Open Access
13 July 2024 Single-beam digital holographic reconstruction: a phase-support enhanced complex wavefront on phase-only function for twin-image elimination
Charlotte Kyeremah, Matthew Weiss, Dila Kandel, Daniel Haehn, Chandra Yelleswarapu
Author Affiliations +
Abstract

Significance

In in-line digital holographic microscopy (DHM), twin-image artifacts pose a significant challenge, and reduction or complete elimination is essential for object reconstruction.

Aim

To facilitate object reconstruction using a single hologram, significantly reduce inaccuracies, and avoid iterative processing, a digital holographic reconstruction algorithm called phase-support constraint on phase-only function (PCOF) is presented.

Approach

In-line DHM simulations and tabletop experiments employing the sliding-window approach are used to compute the arithmetic mean and variance of the phase values in the reconstructed image. A support constraint mask, through variance thresholding, effectively enabled twin-image artifacts.

Results

Quantitative evaluations using metrics such as mean squared error, peak signal-to-noise ratio, and mean structural similarity index show PCOF’s superior capability in eliminating twin-image artifacts and achieving high-fidelity reconstructions compared with conventional methods such as angular spectrum and iterative phase retrieval methods.

Conclusions

PCOF stands as a promising approach to in-line digital holographic reconstruction, offering a robust solution to mitigate twin-image artifacts and enhance the fidelity of reconstructed objects.

1.

Introduction

Digital holographic microscopy (DHM) is a powerful technique that enables capturing and reconstructing three-dimensional (3D) information of phase objects. By leveraging interference patterns created between objects and reference waves, information-rich reconstructed images with high axial and diffraction-limited transverse resolutions are acquired.19 Various configurations such as Gabor’s in-line1,10,11 and Leith-Upatnieks’ off-axis12,13 have been explored to reconstruct the object information from the recorded hologram, each offering its own set of advantages and disadvantages.

In off-axis DHM, the object and reference beams travel separate paths and interfere at the image plane with some angle between them. Even though this scheme eases the reconstruction process, it comes with trade-offs, specifically requiring vibration isolation and underutilization of the recording medium’s space-bandwidth product (SBP). By contrast, in the in-line DHM setup, both object and reference wavefields travel along the same path, allowing for full utilization of the SBP without the need for vibration isolation. However, this setup introduces undesirable twin-image artifacts due to the coincidence of virtual and real images, resulting in an apparent superimposition of out-of-focus images when focusing on either of them, which poses a significant challenge during object reconstruction.

Several techniques have been developed to reduce or eliminate twin-image artifacts. These include conventional diffraction methods such as the Fresnel or the angular spectrum method (ASM),1419 linear filtering,20,21 iterative phase retrieval (ItPR) methods,2228 deep learning,2932 sparsity-based phase retrieval,3335 and complex wavefront filtering.25,3638 ASM involves computing wave propagation through Fourier transforms, enabling analysis of diffraction patterns and reconstructing wavefronts at various distances. It is a widely used numerical approach; however, it cannot effectively eliminate/reduce the unwanted twin-image artifacts, leading to a degradation in the reconstructed object’s quality. Linear filtering techniques alter the spatial or frequency components of the holographic data using filters to selectively remove undesirable components associated with twin images. It is quick and computationally efficient, especially for simple objects. However, it is effective only for simple objects,38 and it may introduce artifacts when applied indiscriminately.

Deep learning employs neural network architectures trained on holographic data to directly learn the mapping between holograms and their corresponding object waves.2932,39,40 This approach effectively mitigates twin-image artifacts but requires a large dataset for training. In addition, complex neural networks may be computationally intensive and require substantial computational resources. By contrast, the sparsity-based method exploits the inherent sparsity within the object domain to improve holographic reconstruction.3335 It can provide high-quality reconstructions with reduced artifacts. However, its performance may vary depending on the sparsity level of the object.

ItPR methods iteratively refine an estimate of the phase until a satisfactory reconstruction is achieved, thus involving an intensive iterative computational process and time. In addition, convergence may not always be guaranteed, and the choice of initialization parameters can significantly impact performance. Complex wavefront filtering techniques manipulate the complex-valued holographic data directly to suppress twin-image artifacts. They can provide high-quality reconstructions without extensive computational overhead. However, designing effective filters may require a deep understanding of the underlying physics and may not always be straightforward for complex scenes.

It is well known that phase is an essential aspect of an image as it contains key features of the object.41 Phase manipulation has been exploited to process images and extract the features of interest.4245 As the phase-only information of the object contains all the required information, here, we report a novel approach called phase-support constraint on phase-only function (PCOF), wherein complex wave filtering is applied to the phase-only function. This technique facilitates object reconstruction using a single hologram, significantly reducing inaccuracies and avoiding iterative processing. Through in-line DHM simulations and tabletop experimental setups involving microbeads, non-centrosymmetric simulated objects, and biological specimens, PCOF’s effectiveness in eliminating the twin-image phenomenon is demonstrated. Quantitative assessments employing metrics such as mean squared error (MSE), peak signal-to-noise ratio (PSNR), and mean structural similarity index (MSSIM) underscore PCOF’s capacity for high-fidelity reconstruction, minimizing deviation from the original image. Furthermore, our method established superior performance on these metrics compared with conventional ASM and ItPR methods. This contribution holds promise for advancing holographic reconstruction methodologies, mitigating artifact-related challenges, and ensuring precise object representation.

2.

Methodology

Figure 1 presents a computer-generated hologram depicting a symbol Ψ as the simulated phase object with a phase step π, along with the algorithmic steps for its reconstruction. Utilizing a sliding-window technique, different parts of the reconstructed image were systematically analyzed by considering subsets or patches of pixels at a time. At each position, the arithmetic means and variance of pixel values within the window were computed, enabling the distinction of areas of significance solely based on phase values. Then, a support constraint mask was formulated through variance thresholding to effectively mitigate twin-image artifacts.

Fig. 1

Workflow of the proposed reconstruction algorithm for twin-image elimination. A computer-generated hologram of a simulated phase object, Ψ with a phase step, π is used as the test object. Step 1 depicts the reconstruction using a conventional diffraction method, also known as the ASM. Steps 2 to 4 illustrate the PCOF, clearly demonstrating the suppression of the twin image.

JBO_29_7_076502_f001.png

Step 1: Initialization and back-propagation

The first step in this proposed method is preprocessing the recorded hologram for contrast and image enhancements. The preprocessed holographic data are then initialized by creating an initial phase of zeros, and the ASM is employed to simulate wave propagation from the hologram plane to the object plane and vice versa. The wavefront Us(x,y,z) at a distance z from the hologram plane is computed using Eq. (1):

Eq. (1)

Us(x,y,z)=F1{F{IH(x,y,0)}·HBP(fx,fy)}Aoeiϕo(x,y),
where IH(x,y,0) is the preprocessed holographic field and F and F1 are the Fourier and inverse Fourier transforms, respectively. The initial phase, ϕo(x,y) obtained from the reconstructed hologram which has the twin image (Aoeiϕo(x,y)) present, is given by

Eq. (2)

ϕo(x,y)=tan1[Im(Us(x,y))Re(Us(x,y))].

Step 2: Segmentation of the twin image

A segmentation mask is created through thresholding of a variance map derived from the reconstructed phase, ϕo(x,y) given by Eq. (2).

2.1.

Calculation of the Variance Map

A variance map, V(x,y), is calculated from the in-focus reconstructed twin image based on the difference between the squared phase values and the arithmetic mean of local pixel values. This variance is a measure of local phase changes and indicates regions with significant features or abrupt phase changes. A specified window size of 5×5 is chosen, and the phase of the reconstructed twin image is processed by calculating the variance in the overlapping blocks using Eq. (3):

Eq. (3)

V=1m2i=1mj=1m[ϕo(x,y)2μ]2,
where μ=(1m2i=1mj=1m[ϕo(x,y)]) is the arithmetic mean of the current block m×m pixels.

2.2.

Thresholding and Masking

A mask that identifies regions of interest within the holographic reconstruction is created. This mask is derived from a threshold applied to the variance map. A threshold, τ, is set based on a percentage of the maximum phase value [using Eq. (2), 50% of the maximum phase value is used here], and pixels with a variance above this threshold are considered part of the segmented region. The threshold value can be adjusted based on the propagation distance, z. A binary mask,

Eq. (4)

ϕmask={1,if  V(i,j)>τ0,otherwise,
is then created based on the thresholding, where pixels above the threshold are set to 1, while the segmented region and pixels below the threshold are set to 0. This enabled the removal/elimination of the unwanted region (twin-image artifacts). The phase mask, ϕmask, is then used to create a new phase (propagated phase mask) using Eq. (5), effectively focusing the reconstruction on regions with high variance. It acts as a constraint, allowing areas with significant phase variations to be emphasized,

Eq. (5)

ϕ(x,y)=ϕo(x,y)·ϕmask.

A Gaussian blur filter is then applied to the updated phase, ϕ(x,y) to reduce noise and smoothen the phase information.

Step 3: Object constraint and forward-propagation

The phase-only function is calculated to preserve key features of the object in the reconstructed image. This approach involves selecting only the phase of, Us(x,y)=Aeiϕ(x,y), while making the amplitude unity. This process eliminates amplitude with negative values as the intensity cannot be negative. Then, the transmission function is updated, and forward propagation is performed to obtain a new holographic field where the amplitude is replaced by the recorded hologram

Eq. (6)

Uo(x,y)=F1{F{(Us(x,y))}·HFP(fx,fy)}Aeiϕo(x,y).

The recorded hologram, IH(x,y,0), is filtered by dividing the hologram by the square of the reference beam, Ur(x,y) using Eq. (7) to suppress the direct current component and create a new hologram

Eq. (7)

IH(x,y)=IH(x,y)Ur2(x,y).

Step 4: Object reconstruction and back-propagation

The new hologram free from the twin-image artifacts is obtained by multiplying the filtered hologram, IH(x,y) with the complex exponential of the phase angle of Uo(x,y)

Eq. (8)

IHnew(x,y)=IH(x,y)×eiϕUo(x,y).

This operation combines the amplitude of the filtered hologram with the phase-only information of the output, Uo(x,y) [Eq. (6)]. The final backpropagation is performed to obtain the reconstructed holographic field for further analysis to extract the true phase (unwrapped phase) map of the object

Eq. (9)

Uobj(x,y)=F1{F{IHnew(x,y)}·HBP(fx,fy)}Aeiϕobj(x,y).

2.3.

Evaluation of the Algorithm and Image Quality

To facilitate a comprehensive assessment of various algorithms and their corresponding reconstruction outcomes, MSE is used as a metric to quantify the error or discrepancy between the ground truth object and its reconstructed image. The MSE is normally defined as follows:

Eq. (10)

MSE(i,j)=1M×Ni=1MMj=1NN[ρρo]2,
where ρ(i,j) is the reconstructed distribution, ρo(i,j) is the initial distribution, and [N, M] is the number of pixels. The MSE, calculated as the average of squared differences among pixel values, serves as an indicator of how closely the pixel values in the reconstructed image resemble those in the original image. A lower MSE signifies a closer match between two images.

PSNR is another metric used to measure the quality of reconstructed images compared with the original images. It is a logarithmic ratio of the maximum possible pixel value to the square root of the MSE.

Eq. (11)

PSNR=10·log10(max(ρ)MSE)2,
where max(ρ) is the maximum possible pixel value of the object’s phase map. Since PSNR is inversely related to the MSE, as the MSE decreases, the PSNR increases, indicating a higher-quality reconstruction. A higher PSNR denotes increased similarity between the original and reconstructed images.

MSSIM is a metric used to quantify the similarity between two images (ρ,ρo). It assesses structural information preservation during image processing. The MSSIM value ranges between 1 and 1, where 1 indicates perfect similarity, 0 denotes no similarity, and negative values suggest dissimilarity. A higher MSSIM value implies better preservation of structural details between the original and reconstructed images.

Eq. (12)

SSIM(ρ,ρo)=[2μρμρO+C1]×[2σρρO+C2][μρ2+μρO2+C1]×[σρ2+σρO2+C2],
where μρ, μρO, σρ, σρO, and σρρO are the average luminance of images, standard deviations of the pixel intensities, and the covariance of the pixel intensities in images ρ and ρo. C1 and C2 are constants added for numerical stability. MSSIM is computed by comparing local patterns of pixel intensities in both images, offering insights into structural similarity. A higher MSSIM signifies a greater resemblance between the structures of the original and reconstructed images.

3.

Results and Discussion

The efficacy of our proposed methodology is demonstrated by recording and reconstructing microscopic objects via both in-line DHM simulations and tabletop experimental setups. Ground truth objects were generated with precise optical properties for the simulation phase. For instance, a configuration of 150 microsphere beads, each measuring a diameter of 30  μm and a phase step of 2π, was simulated [depicted in Fig. 2(a)]. Holograms were generated using the angular spectrum propagation function and then reconstructed using the proposed approach. The reconstructions were compared against established ASM14,17,19 and ItPR26,46 methods. Figures 2(b)2(d) show the reconstructed microspheres using ASM, ItPR, and the PCOF methods, respectively. The phase profiles of the reconstructed images show “ringing” around reconstructed phase object information in the case of ASM and ItPR, whereas no ringing is observed in the PCOF reconstruction. The absence of ringing highlights its superior performance in effectively suppressing twin-image artifacts and ensures accurate representation and true preservation of the object information. The maximum phase value obtained post-reconstruction closely matched the predefined maximum phase shift of 2π radians, differing by only 0.312 radians.

Fig. 2

Reconstruction and comparison of simulated 30  μm beads. (a) Full FOV hologram of 150 beads; panels (b)–(d) are the reconstructed phase of the region of interest (ROI) (red box) by ASM, ItPR after 30 iterations, and PCOF methods. The inset is the phase profile of the selected single bead; panels (e) and (f) are the 2D and 3D profiles of the ROI by PCOF, respectively. The 3D profile shows that the PCOF can be used to reconstruct objects at different depths. The phase profiles of panels (b)–(d) are for observation.

JBO_29_7_076502_f002.png

To provide a deeper understanding of fidelity, quality, and similarity between the original and reconstructed images, quantitative analysis using multiple metrics, such as the MSE, PSNR, and MSSIM, was performed, and the values are presented in Table 1. Both ASM and ItPR have an MSE of 0.02, indicating a greater deviation from the original image. By contrast, PCOF exhibited an exceptionally low MSE of 0.01, signifying a closer match between the pixel values in the reconstructed and original images. The PSNR values of 65.56 dB for the ASM, 66.78 dB for the ItPR method, and 68.40 dB for PCOF reflect varying reconstruction qualities, with the proposed method outperforming both other methods. The MSSIM values further support these findings. The ItPR method achieves substantial structural similarity with an MSSIM of 0.9987 but falls short of near-perfect resemblance. By contrast, the ASM’s MSSIM of 0.9980 indicates relatively lower structural similarity compared with the other methods. Remarkably, PCOF achieves an almost perfect structural similarity, with an MSSIM of 0.9992, highlighting its exceptional ability to reconstruct images closely resembling the original data. Overall, the PCOF method outperforms ASM and ItPR methods, with lower MSE, higher PSNR, and markedly higher MSSIM values. These findings highlight PCOF’s exceptional capability for high-fidelity reconstructions, affirming its effectiveness in digital holographic image reconstruction.

Table 1

Quantitative metrics were obtained for the reconstruction of the simulated microbeads using the named reconstruction methods.

MethodMSEPSNR (dB)MSSIM
ASM0.0265.560.9980
ItPR0.0266.780.9987
PCOF0.0168.400.9992

To demonstrate the application of the PCOF for non-centrosymmetric objects, an extended object resembling a spiral was simulated.46 The reconstruction employed a threshold set at 50% of the maximum phase value. The resulting output, as shown in Fig. 3, indicates high fidelity and substantial structural similarity between the original and reconstructed images. The quantitative evaluation produced an MSE of 0.02, a PSNR of 64.99 dB, and an MSSIM of 0.9970.

Fig. 3

Simulated and reconstructed hologram of an extended object. (a) The ground truth image, (b) a recorded hologram, and (c) a reconstructed image of the extended object using the PCOF method.

JBO_29_7_076502_f003.png

3.1.

Experimental Verification

A tabletop in-line lensless DHM system was built to study the effectiveness of PCOF in analyzing real-phase objects. The setup comprises a partially coherent light emitting diode (LED) light source (Thorlabs, LED525E, Newton, New Jersey, United States, 2.6 mW @ 20 mA, wavelength λ=525  nm), which passes through a 50  μm pinhole for spatial filtering. The sample is positioned 10  cm from the light source by employing detector plane holography, while the detector sits 1  mm away. This arrangement yields a field of view (FOV) equivalent to the active area of the detector, 30  mm2 for a Raspberry Pi camera boasting a pixel size of 1.55  μm and a resolution of 4056×3040  pixels.

The imaging of red blood cells over a large FOV and the subsequent generation of 3D phase maps offer a significant value for blood and disease analyses. Despite the relatively low range of phase values, the sample’s thickness can be accurately calculated following phase characterization using a known reference sample. To illustrate the utility of PCOF for such analysis, phase characterization through the imaging of 2.5  μm polystyrene silica beads (Polysciences, Inc., Warrington, Pennsylvania, United States) with a refractive index of 1.43 was conducted. The beads were suspended in water (nmed=1.33) and sandwiched between two coverslips. The hologram of the beads was captured using the in-line lensless (DHM) setup. The phase reconstruction obtained using PCOF is depicted in Fig. 4(b). Analysis of the selected glass bead within the region of interest revealed an average size of 2.55±0.20  μm, aligning closely with the manufacturer’s specifications.

Fig. 4

Phase characterization and thickness profile of silica beads. (a) Hologram of silica beads of 2.5  μm size, (b) reconstructed phase, (c) reconstructed intensity/amplitude of the white-lined single bead, and (d) the thickness profile of the reconstructed single bead.

JBO_29_7_076502_f004.png

Further imaging and reconstruction of the cheek cell sample are depicted in Fig. 5. Initially, the sample was verified using a standard bright-field microscope, as shown in Fig. 5(a). Subsequently, the sample was imaged using an in-line DHM setup, and the unwrapped phase was obtained, as shown in Fig. 5(b). Intensity (box I) and phase reconstruction (box II) were performed using the ASM [Fig. 5(c)], ItPR [Fig. 5(d)], and PCOF [Fig. 5(e)] methods. The bottom row of box II displays phase profiles. The results demonstrated that the PCOF method significantly reduced twin-image artifacts present in the ASM and ItPR reconstructions. Furthermore, the unwrapped image reveals all the internal structures of the cell, including the nucleus, cytoplasm, and cell membrane.

Fig. 5

Hologram recording and reconstruction of cheek cells. (a) Brightfield image obtained using a 5×0.15 numerical aperture microscope objective, and the intensity profile; (b) in-line hologram and 3D phase maps obtained using PCOF; (c)–(e) reconstructions using ASM, ItPR, and PCOF methods, respectively; I: intensity reconstruction; II: phase reconstruction (top row) and phase profile (bottom row). Intensity and phase profiles of a cheek cell are obtained along the indicated white line.

JBO_29_7_076502_f005.png

It can be observed that minimal residual speckle noise is present in our images. However, it is important to note that such noise is typical in coherent or partially coherent imaging. Although LEDs produce less speckle than lasers, the use of a pinhole in the in-line DHM setup can still introduce some speckle noise due to diffraction and scattering. The phase-support constraint in our method helps to limit the phase values to a defined support region, reducing the influence of noise outside this region. By focusing the reconstruction process within the defined support, the method effectively filtered out noise and artifacts that are not consistent with the expected phase distribution.

4.

Conclusion

The significance of phase information in imaging techniques motivated the development of complex wavefront using PCOF as a novel approach to tackle twin-image artifacts in digital holography microscopy. PCOF’s broad utility was showcased through extensive simulations, tabletop experiments, and quantitative evaluations utilizing various metrics obtained with cheek cells, glass beads, and non-centrosymmetric objects. Comparative analysis with conventional methods such as ASM and ItPR highlighted PCOF’s superiority in terms of MSE, PSNR, and MSSIM values, emphasizing its potential for advancing holographic reconstruction methodologies. PCOF offers a promising solution by swiftly suppressing twin-image artifacts and accurately reconstructing the object phase without iterative processes, thereby ensuring high-quality holographic image reconstructions. Future endeavors could involve exploring the application of the developed system in real-time imaging and diagnostics, which holds promise for further investigation.

Disclosures

The authors declare that there are no conflicts of interest related to this article.

Code and Data Availability

Data underlying the results may be obtained from the authors upon reasonable request.

Acknowledgments

The UMass Boston Internal Proposal Development Grant Program funded this work. C.K. thanks the Proposal Development Grant for the financial support. MW was supported by the John Templeton Foundation (Grant No. 62424). The opinions expressed in this publication are those of the authors and do not necessarily reflect the views of the John Templeton Foundation.

References

1. 

D. Gabor, “A new microscopic principle,” Nature, 161 (4098), 777 –778 https://doi.org/10.1038/161777a0 (1948). Google Scholar

2. 

S. S. Kou and C. J. Sheppard, “Imaging in digital holographic microscopy,” Opt. Express, 15 (21), 13640 –13648 https://doi.org/10.1364/OE.15.013640 OPEXFF 1094-4087 (2007). Google Scholar

3. 

M. K. Kim, “Principles and techniques of digital holographic microscopy,” Proc. SPIE, 1 (1), 018005 https://doi.org/10.1117/6.0000006 PSISDG 0277-786X (2010). Google Scholar

4. 

J. Sheng, E. Malkiel and J. Katz, “Digital holographic microscope for measuring three-dimensional particle distributions and motions,” Appl. Opt., 45 (16), 3893 –3901 https://doi.org/10.1364/AO.45.003893 APOPAI 0003-6935 (2006). Google Scholar

5. 

B. Kemper and G. Von Bally, “Digital holographic microscopy for live cell applications and technical inspection,” Appl. Opt., 47 (4), A52 –A61 https://doi.org/10.1364/AO.47.000A52 APOPAI 0003-6935 (2008). Google Scholar

6. 

M. K. Kim, Digital Holographic Microscopy, 149 –190 Springer, New York (2011). Google Scholar

7. 

U. Schnars and W. Jüptner, “Direct recording of holograms by a CCD target and numerical reconstruction,” Appl. Opt., 33 (2), 179 –181 https://doi.org/10.1364/AO.33.000179 APOPAI 0003-6935 (1994). Google Scholar

8. 

X. Yu et al., “Review of digital holographic microscopy for three-dimensional profiling and tracking,” Opt. Eng., 53 (11), 112306 https://doi.org/10.1117/1.OE.53.11.112306 (2014). Google Scholar

9. 

P. Marquet et al., “Digital holographic microscopy: a noninvasive contrast imaging technique allowing quantitative visualization of living cells with subwavelength axial accuracy,” Opt. Lett., 30 (5), 468 –470 https://doi.org/10.1364/OL.30.000468 OPLEDP 0146-9592 (2005). Google Scholar

10. 

Dennis Gabor, “Microscopy by reconstructed wave-fronts,” Proc. R. Soc. Lond. A Math. Phys. Sci., 197 (1051), 454 –487 https://doi.org/10.1098/rspa.1949.0075 (1949). Google Scholar

11. 

D. Gabor, W. E. Kock and G. W. Stroke, “Holography,” Science, 173 (3991), 11 –23 https://doi.org/10.1126/science.173.3991.11 SCIEAS 0036-8075 (1971). Google Scholar

12. 

E. N. Leith and J. Upatnieks, “Reconstructed wavefronts and communication theory,” J. Opt. Soc. Am., 52 (10), 1123 https://doi.org/10.1364/JOSA.52.001123 JOSAAH 0030-3941 (1962). Google Scholar

13. 

E. N. Leith and J. Upatnieks, “Wavefront reconstruction with continuous-tone objects,” J. Opt. Soc. Am., 53 (12), 1377 https://doi.org/10.1364/JOSA.53.001377 JOSAAH 0030-3941 (1963). Google Scholar

14. 

T. Kreis, Handbook of Holographic Interferometry: Optical and Digital Methods, John Wiley & Sons( (2006). Google Scholar

15. 

J. W. Goodman, “An introduction to the principles and applications of holography,” Proc. IEEE, 59 (9), 1292 –1304 (1971). Google Scholar

16. 

B. Javidi et al., “Roadmap on digital holography [Invited],” Opt. Express, 29 (22), 35078 –35118 https://doi.org/10.1364/OE.435915 OPEXFF 1094-4087 (2021). Google Scholar

17. 

G. Dwivedi et al., “Comparison of numerical reconstruction of digital holograms using angular spectrum method and Fresnel diffraction method,” J. Opt., 49 1 –10 https://doi.org/10.1007/s12596-017-0424-z (2017). Google Scholar

18. 

N. Akhter et al., “A comparative study of reconstruction algorithms in digital holography,” Optik – Int. J. Light Electron Opt., 124 (17), 2955 –2958 https://doi.org/10.1016/j.ijleo.2012.09.002 (2013). Google Scholar

19. 

J. W. Goodman, Introduction to Fourier Optics, Roberts and Company Publishers( (2005). Google Scholar

20. 

K. H. S. Marie, J. C. Bennett and A. P. Anderson, “Digital processing technique for suppressing the interfering outputs in the image from an inline hologram,” Electron Lett., 15 (8), 241 –243 https://doi.org/10.1049/el:19790171 (1979). Google Scholar

21. 

L. Onural, “Digital decoding of in-line holograms,” Opt. Eng., 26 (11), 261124 https://doi.org/10.1117/12.7974205 (1987). Google Scholar

22. 

M. Guizar-Sicairos and J. R. Fienup, “Understanding the twin-image problem in phase retrieval,” J. Opt. Soc. Amer. A, 29 (11), 2367 –2375 https://doi.org/10.1364/JOSAA.29.002367 JOAOD6 0740-3232 (2012). Google Scholar

23. 

J. R. Fienup, “Phase retrieval algorithms: a comparison,” Appl. Opt., 21 (15), 2758 –2769 https://doi.org/10.1364/AO.21.002758 APOPAI 0003-6935 (1982). Google Scholar

24. 

G. Koren, F. Polack and D. Joyeux, “Iterative algorithms for twin-image elimination in in-line holography using finite-support constraints,” J. Opt. Soc. Amer. A, 10 (3), 423 –433 https://doi.org/10.1364/JOSAA.10.000423 JOAOD6 0740-3232 (1993). Google Scholar

25. 

L. Denis et al., “Twin-image noise reduction by phase retrieval in in-line digital holography,” Proc. SPIE, 5914 59140J https://doi.org/10.1117/12.617405 (2005). Google Scholar

26. 

T. Latychevskaia, “Iterative phase retrieval for digital holography: tutorial,” J. Opt. Soc. Amer. A, 36 (12), D31 –D40 https://doi.org/10.1364/JOSAA.36.000D31 JOAOD6 0740-3232 (2019). Google Scholar

27. 

G. Liu and P. D. Scott, “Phase retrieval and twin-image elimination for in-line Fresnel holograms,” J. Opt. Soc. Amer. A, 4 (1), 159 –165 https://doi.org/10.1364/JOSAA.4.000159 JOAOD6 0740-3232 (1987). Google Scholar

28. 

S. Kumar Gautam et al., “Single-shot and twin-image free unique phase retrieval using an aspect of autocorrelation that considers the object asymmetry,” J. Opt. Soc. Amer. A, 37 (11), 1826 –1831 https://doi.org/10.1364/JOSAA.402731 JOAOD6 0740-3232 (2020). Google Scholar

29. 

Y. Rivenson et al., “Phase recovery and holographic image reconstruction using deep learning in neural networks,” Light Sci. Appl., 7 (2), 17141 https://doi.org/10.1038/lsa.2017.141 (2017). Google Scholar

30. 

A. Sinha et al., “Lensless computational imaging through deep learning,” Optica, 4 (9), 1117 –1125 https://doi.org/10.1364/OPTICA.4.001117 (2017). Google Scholar

31. 

H. Li et al., “Deep DIH: single-shot digital in-line holography reconstruction by deep learning,” IEEE Access, 8 202648 –202659 https://doi.org/10.1109/ACCESS.2020.3036380 (2020). Google Scholar

32. 

L. Huang et al., “Holographic image reconstruction with phase recovery and autofocusing using recurrent neural networks,” ACS Photonics, 8 (6), 1763 –1774 https://doi.org/10.1021/acsphotonics.1c00337 (2021). Google Scholar

33. 

Y. Gao and L. Cao, “Iterative projection meets sparsity regularization: towards practical single-shot quantitative phase imaging with in-line holography,” Light: Adv. Manuf., 4 (1), 37 –53 https://doi.org/10.37188/lam.2023.006 (2023). Google Scholar

34. 

C. Gaur, P. Lochab and K. Khare, “Experiments on sparsity assisted phase retrieval of phase objects,” J. Opt., 19 (5), 055703 https://doi.org/10.1088/2040-8986/aa6357 (2017). Google Scholar

35. 

A. S. Galande et al., “Quantitative phase imaging of biological cells using lensless inline holographic microscopy through sparsity-assisted iterative phase retrieval algorithm,” J. Appl. Phys., 132 (24), 243102 https://doi.org/10.1063/5.0123677 JAPIAU 0021-8979 (2022). Google Scholar

36. 

G. Pedrini et al., “In-line digital holographic interferometry,” Appl. Opt., 37 (26), 6262 –6269 https://doi.org/10.1364/AO.37.006262 APOPAI 0003-6935 (1998). Google Scholar

37. 

J. L. de Almeida et al., “Twin-image suppression in digital in-line holography based on wave-front filtering,” Pattern Anal. Appl., 24 (3), 907 –914 https://doi.org/10.1007/s10044-020-00949-7 (2021). Google Scholar

38. 

C. P. McElhinney, B. M. Hennelly and T. J. Naughton, “Twin-image reduction in inline digital holography using an object segmentation heuristic,” J. Phys. Conf. Ser., 139 012014 https://doi.org/10.1088/1742-6596/139/1/012014 JPCSDZ 1742-6588 (2008). Google Scholar

39. 

X. Chen et al., “DH-GAN: a physics-driven untrained generative adversarial network for holographic imaging,” Opt. Express, 31 (6), 10114 https://doi.org/10.1364/OE.480894 OPEXFF 1094-4087 (2023). Google Scholar

40. 

M. Rogalski et al., “Physics-driven universal twin-image removal network for digital in-line holographic microscopy,” Opt. Express, 32 (1), 742 https://doi.org/10.1364/OE.505440 OPEXFF 1094-4087 (2024). Google Scholar

41. 

A. V. Oppenheim and J. S. Lim, “The importance of phase in signals,” Proc. IEEE, 69 (5), 529 –541 https://doi.org/10.1109/PROC.1981.12022 (1981). Google Scholar

42. 

S.-R. Kothapalli et al., “Spectral phase based medical image processing1,” Acad. Radiol., 12 (6), 708 –721 https://doi.org/10.1016/j.acra.2004.09.017 (2005). Google Scholar

43. 

T. L. Nguyen et al., “Quantitative phase imaging: recent advances and expanding potential in biomedicine,” ACS Nano, 16 (8), 11516 –11544 https://doi.org/10.1021/acsnano.1c11507 ANCAC3 1936-0851 (2022). Google Scholar

44. 

J. Birdi et al., “True 3D reconstruction in digital holography,” J. Phys.: Photonics, 2 (4), 044004 https://doi.org/10.1088/2515-7647/abb586 (2020). Google Scholar

45. 

C. Kyeremah et al., “Exploring different textures of a nematic liquid crystal for quantitative Fourier phase contrast microscopy,” Opt. Laser Technol., 147 107631 https://doi.org/10.1016/j.optlastec.2021.107631 OLTCAS 0030-3992 (2022). Google Scholar

46. 

T. Latychevskaia and H.-W. Fink, “Solution to the twin image problem in holography,” Phys. Rev. Lett., 98 (23), 233901 https://doi.org/10.1103/PhysRevLett.98.233901 PRLTAO 0031-9007 (2007). Google Scholar

Biography

Charlotte Kyeremah earned her bachelor’s (2013) and master’s (2017) degrees from the University of Cape Coast, Ghana, inspired by Prof. Buah-Bassuah to pursue applied optics. She then obtained a second master’s degree in applied physics from the University of Massachusetts Boston (2020), under Prof. Chandra Yelleswarapu. Currently, Charlotte is a PhD candidate in biophysics at the same university, focusing on holographic microscopy for biomedical imaging, aiming to enhance infectious disease diagnostics in developing countries such as Ghana.

Matthew Weiss: Biography is not available.

Dila Kandel obtained his MSc degree in physics from Tribhuvan University, Nepal, in 2013 and his second MS degree in applied physics in 2023 from the University of Massachusetts Boston. He has taught high school physics in Nepal. Currently, he is conducting research in particle tracking using inline digital holographic microscopy under the supervision of Prof. Chandra Yelleswarapu. He is proficient in Java, Python, MATLAB, and R, with a strong background in data analysis and image processing.

Daniel Haehn is an assistant professor of computer science at the University of Massachusetts Boston (UMB) and an associate of the Harvard Paulson School of Engineering and Applied Sciences. Prior to joining UMB, he was a researcher at Brigham and Women’s Hospital, the University of Pennsylvania, and Boston Children’s Hospital. Daniel obtained his PhD in computer science from Harvard University and his diploma (BS/MS) in medical computer science from the University of Heidelberg in Germany. He strongly supports open science and open education.

Chandra Yelleswarapu is an associate professor of physics at the University of Massachusetts Boston. He obtained his MPhil degree from the University of Hyderabad and his PhD in applied optics from Alabama A&M University. His research focuses on nonlinear optics, photoacoustics, and quantitative phase imaging. He is committed to diversity in STEM, engaging underrepresented students in research programs and university-wide initiatives. He aims to translate discoveries into real-world solutions.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 International License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Charlotte Kyeremah, Matthew Weiss, Dila Kandel, Daniel Haehn, and Chandra Yelleswarapu "Single-beam digital holographic reconstruction: a phase-support enhanced complex wavefront on phase-only function for twin-image elimination," Journal of Biomedical Optics 29(7), 076502 (13 July 2024). https://doi.org/10.1117/1.JBO.29.7.076502
Received: 14 March 2024; Accepted: 24 June 2024; Published: 13 July 2024
Advertisement
Advertisement
KEYWORDS
3D image reconstruction

Digital holography

Holography

Phase reconstruction

Holograms

Image restoration

Reconstruction algorithms

Back to Top