Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Coherent noise reduction in digital holographic microscopy by averaging multiple holograms recorded with a multimode laser

Open Access Open Access

Abstract

In digital holographic microscopy (DHM), it is undesirable to observe coherent noise in the reconstructed images. The sources of the noise are mainly the parasitic interference fringes caused by multiple reflections and the speckle pattern caused by the optical scattering on the object surface. Here we propose a noise reduction approach in DHM by averaging multiple holograms recorded with a multimode laser. Based on the periodicity of the temporal coherence of a multimode semiconductor laser, we acquire a series of holograms by changing the optical path length difference between the reference beam and object beam. Because of the use of low coherence light, we can remove the parasitic interference fringes caused by multiple reflections in the holograms. In addition, the coherent noise patterns change in this process due to the different optical paths. Therefore, the coherent noise can be reduced by averaging the multiple reconstructions with uncorrelated noise patterns. Several experiments have been carried out to validate the effectiveness of the proposed approach for coherent noise reduction in DHM. It is shown a remarkable improvement both in amplitude imaging quality and phase measurement accuracy.

© 2017 Optical Society of America

1. Introduction

Digital holographic microscopy (DHM) [1] is a well-known technique that can optically record and numerically reconstruct 3D information of an object. In the reflection mode, the recovered phase provides surface topography information and permits direct measurement of the object shape and the associated parameters. In the transmission mode, the recovered phase of the optical path length describes the phase function of the object. Therefore, DHM has led to the development of a huge set of applications in many fields, including quantitative microscopy in biology [2], nondestructive measurement in the industry [3], and etc. However, due to the coherent nature of the light source, the digital holograms are corrupted by a mixture of coherent noises when a laser beam illuminates an object with a rough surface or inhomogeneous interior. These noises include the parasitic interference fringes originating from the multiple reflections from different interfaces in the optical path and the unwanted scattering on the micro structures with sizes larger than the wavelength [4]. As such, the image quality and the measurement accuracy are seriously degraded in the holographic reconstructions. In order to reduce coherent noises, a number of techniques have been proposed and they can be classified into four groups: low coherence methods [5–12], multiple holograms methods [13–24], image processing methods [25–34] and hybrid methods [35].

In the first group, some approaches are based on temporally low coherent light emitted from light-emitting diodes [5, 6] or super luminescent diode [7], spatially low coherent light by use of rotating diffusers [8–10], and a random laser illumination [11]. By combining the Gabor’s holography and Zernike’s phase contrast imaging, a speckle-free phase microscopy is presented by using Halogen Lamp illumination with extremely short coherence length [12]. These light sources can guarantee that no multi-reflected-transmitted lights will interfere, and thus, remove deleterious parasitic interference fringes in the recorded holograms. Furthermore, the relationship between the degrees of spatiotemporal coherence of an illumination and the imaging quality of interferometric microscopy has been quantitatively investigated by experiments [13]. In the second group of multiple holograms approach, the reduction of coherent noise is achieved by averaging a number of reconstructed images from a set of holograms, whose noise patterns are decorrelated with each other. The decorrelation of noise can be achieved by using different wavelengths [14], using different polarization states [15, 16], varying the positions of camera [17], shifting the object slightly [18], changing illumination angle [19, 20], and applying a random phase generator, such as an optical diffuser [21, 22], a spatial light modulator [23]. As a related approach, the synthetic aperture imaging technique was also used to reduce coherent noises in DHM [24, 25]. For this strategy, the denoising performance depends on the noise decorrelation extent between the different holograms. The methods based on imaging processing were proposed with the use of digital filters [26], such as the median filter [27], the Wiener filter [28], the discrete Fourier filter [29], the wavelet filtering [30], the non-local means [31], and a multilevel two-dimensional empirical decomposition mode [32]. As an alternative approach, the coherent noise can also be suppressed by averaging a certain number of images, which are reconstructed from a single hologram with resampling masks [33, 34] or numerical jittering [35]. In these way, the coherent noise reduction can be performed at the cost of resolution and contrast of the reconstructed image. More recently, a novel approach is proposed by combining the Multi-Look Digital Hologram (MLDH) techniques and Three-Dimension Block Matching (BM3D) filtering to reduce noise in DH [36]. The method is made of two steps. Firstly the noise is reduced to the theoretical bound by averaging multiple reconstructed images with uncorrelated speckle noises. And then, based on the averaging result, a quasi noise-free reconstruction is obtained by use of the block grouping and collaborative filtering. However, the method’s capability was only proved in holographic intensity imaging, and the results are not provided for holographic phase imaging.

In this paper, we present a new approach to reduce the coherent noise by combining the low coherence recordings and multiple holograms averaging. First, by employing the periodicity of temporal coherence of multimode semiconductor laser, we record a series of holograms by changing the optical path length of the reference beam. Due to the use of low coherent light, the parasitic interference fringes caused by multiple reflections or transmissions are avoided in holograms recording. Second, we numerically reconstruct a set of intensity and phase images with uncorrelated coherent noise patterns. Third, we average the reconstructed images to suppress the coherent noise caused by the optical scattering. By performing these three steps, we can significantly improve the image quality and the measurement accuracy. The experimental results demonstrate the effectiveness of the proposed method for the reconstructed intensity and phase images.

2. Principle and method

We consider the process of digital holography recording: a reference beam R interferences with an object beam O, both beams being emitted from the same light source. The interference pattern takes the expression as

IH=|R|2+|O|2+RO+RO.
where the superscript * denotes the complex conjugation of the wave amplitude.

Supposed that the light source is a multimode semiconductor laser and the oscillation of N longitudinal modes, the complex optical field can be written as

E=n=(N1)/2n=(N1)/2Anexp{i[(k+nΔk)L+φn]}.
where An and φn are the optical power and phase of the nth lasing mode respectively, k = 2π/λ, λ is the center mode wavelength of laser, Δk = 2πΔλ/λ2, Δλ is the wavelength interval between two adjacent modes, L is the optical transmission path.

Thus, the R and O beams can be written as, respectively

R=R0n=(N1)/2n=(N1)/2Anexp{i[(k+nΔk)LR+φn]}.O=O0n=(N1)/2n=(N1)/2Anexp{i[(k+nΔk)LO+φn]}.
where R0 is the complex amplitude of R beam, which usually is a plane or spherical wave, O0 is the complex amplitude of O beam, which represents the surface topography for reflection imaging and the optical path difference for transmission imaging, LR and LO are the optical paths of R and O arms, respectively.

So, the interference term of R*O in Eq. (1), which contains the information of object, can be described as,

R*O=R0O0n=(N1)/2n=(N1)/2m=(N1)/2m=(N1)/2AnAmexp[ik(LOLR)]×exp[iΔk(mLOnLR)]exp[i(φmφn)].

Because the interference of different modes does not mutually coherent, the phase different of φm-φn (n≠m) is likely to vary randomly in time. So, the interference of different modes does not contribute to the interference term, and it follows that

R*O=R0O0exp[ik(LOLR)]n=(N1)/2n=(N1)/2m=(N1)/2m=(N1)/2δn,mAnAm×exp[iΔk(nLOmLR)]exp[i(φmφn)].
where δn,m is Kronecker delta function and described as,

δn,m={1(n=m)0(nm)

According to Eq. (5), the inference of different modes does not pay role to the holographic recording, and then each mode contributes for the recording of its respective hologram. The resulting interference pattern will be a superposition of all holograms, and the interference term can be written as

R*O=R0O0exp[ik(LOLR)]n=(N1)/2n=(N1)/2An2exp[iΔkn(LOLR)].

In order to simplify the following analysis, the value of An for each mode is considered equal to 1. Now, Eq. (7) can be written as

R*O=R0*O0exp[ik(LOLR)]{sin[NΔk(LOLR)/2]sin[Δk(LOLR)/2]}.
where the sin denotes the sine function, the term in the {·} corresponds to the visibility of interference fringes.

Based on Eq. (8), the graph in Fig. 1(b) shows the behavior of fringe visibility as a function of (LO-LR) with Δk = 1.046rad/mm and N = 12, which are close to the experimental condition. It can be seen that the fringe visibility repeats periodically with a spacing of 6.75 mm, which is called the interference period (IP). In order to verify this effect, a series of holograms are recorded by changing the value of LR and keeping the value of LO constant in the experiment setup, which is shown in the Fig. 2. Figure 1(a) represents a hologram with high fringe visibility, and Figs. 1(c), 1(d) and 1(e) show the cutouts of same region (indicated with a white rectangle) in these holograms, which are recorded with different (LO-LR). As can be seen, when the changing of (LO-LR) is non-integer multiple of the IP, the fringe visibility is low, as shown in the Figs. 1(d) and 1(e). In this case, the object information cannot be recorded in the hologram. On the other hand, when the value of (LO-LR) are an integer multiple of the IP, the fringe visibility are maximum, as shown in the Fig. 1(c). At these positions, the hologram has a high signal to noise ratio. Meanwhile, due to the changing of (LO-LR), the coherent noise in the holograms will randomly fluctuate. Thus, by digital reconstruction, a series of intensity and phase images are obtained with uncorrelated coherent noise patterns from these holograms. Moreover, thanks to the shortly coherent light, it is possible to guarantee no parasitic interference fringes caused by multiple reflections in holograms. And then, the averaging process can effectively suppress the artifact speckle pattern caused by the optical scattering on object surface or interior. Consequently, the imaging quality and the measurement accuracy can be significantly improved.

 figure: Fig. 1

Fig. 1 The fringe visibility versus the optical path different between the R and O beams. (a): the digital hologram, (b): the calculation results, (c) (d) and (e): the experimental results, the cutouts (indicated with a white rectangle) of the holograms recorded with different (LO-LR).

Download Full Size | PDF

 figure: Fig. 2

Fig. 2 Schematic of the experimental system. Laser, a multimode semiconductor laser; λ/2, half-wave plate; PBS, polarizing beam splitter; BE, beam expander with spatial filter; M, mirror; L, lens; O, object wave; R, reference wave; MO, microscope objective; BS, beam splitter; CMOS, camera; OPR, optical path retarder to adjust the optical path length of the reference beam.

Download Full Size | PDF

3. Experiments

The experimental setup is a reflective off-axis holographic interferometer with a sample pre-magnifying configuration, as depicted in Fig. 2. The source is a multimode semiconductor laser with the oscillation of 12 longitudinal modes (Laser, Oxxius LBX, with 638nm wavelength and 50mW power). The laser beam is divided into an object beam (O) and a reference beam (R) by a polarized beam splitter (PBS). A half-wave plate (λ/2) in front of the PBS is used to adjust the intensity ratio between the O and R beams. Behind the PBS, another half-wave plate (λ/2) is introduced to obtain the parallel direction of polarization in two beams. Both the beams are collimated as plane waves by the beam expanders (BE1 and BE2, composed of spatial filters and collimating lenses). In the object arm, a microscope objective (MO, Olympus UplanFLN, the magnification of 20 × , the numerical aperture of 0.4) is used to enhance the lateral resolution, which is diffraction limited and depends on the numerical aperture of the MO. The object beam is focused near the back focal plan of the MO by the condenser lens (L2) to deliver a collimated illumination wave. Then, the optical wave reflected by the object is collected by the MO and a magnified image of the object is produced behind the camera with the distance of 8 cm. In the reference arm, the optical path length of the R beam is adjusted by moving an optical path retarder (OPR, composed by the mirror of M1 and M2) to match the optical path length of the O beam. In order to improve the sampling capacity of the camera, the matching lens (L1) is optionally used to create a spherical wave with a curvature very similar to the curvature created by the MO. Finally, by calibrating the position of the OPR, the interference between the object wave emerging from the image and the spherical reference beam is accomplished with a maximal fringe contrast, and recorded by the camera (CMOS, Lumenera Lu125M, the resolution of 1024 × 1024, the pixel pitch of 6.7μm). As mentioned above, due to the periodicity of laser temporal coherence, a sequence of holograms can be recorded by continuously adjusting the OPR and keeping the position of the object with the IP.

The goal of the first experiment is to investigate the fluctuation of coherent noise with changing the value of LR-LO. A 1951 USAF resolution target is used as the observed object, and a sequence of holograms is recorded by keeping the LO and changing the LR by an integer multiple of IP. The reconstructed results are shown in Fig. 3. Figures 3(a) and 3(b) represent the amplitude images with a reference LR and a shifted LR, where the distance between the LR and LO is equal to one IP. Figures 3(c) and 3(d) represent the phase images with the reference LR and the shifted LR with one IP. In these figures, the enlarged cutouts of a uniform region indicated by a white rectangle show the noise patterns with more detail. Obviously, the coherent noise distribution has randomly fluctuated with the changing of LR in both amplitude and phase images. Furthermore, in order to qualitatively discuss the statistical relation of coherent noise, the normalized cross correlations (γ) between two noise patterns is calculated by

γ=n,m[ϕ0(n,m)ϕ¯0][ϕL(n,m)ϕ¯L]n,m[ϕ0(n,m)ϕ¯0]2n,m[ϕL(n,m)ϕ¯L]2.
where ϕ0 and ϕL denotes either the amplitude or phase distribution in the reconstructed image with the reference LR0 and the shifted LRS, ϕ¯0 and ϕ¯Ldenote the mean values, n and m denote the pixel number of the row and column in region of interest. Based on the uniform region as mentioned above, the graph in Fig. 3(e) displays the cross correlation against the value of LR-LO, where the changing of LR-LO is expressed as the multiple of the IP. It can be seen that γ of different positions are less than 0.2. As predicted, the results show a very efficient noise decorrelation between the reconstructed images with different optical path difference between O and R beams.

 figure: Fig. 3

Fig. 3 The experimental results of the noise decorrelation by changing the LR and keeping LO. (a) and (b): the amplitude images with a reference LR and a shifted LR, (c) and (d): the phase images with a reference LR and a shifted LR, (e): the normalized cross correlations (γ) of the noise patterns over the LR-LO, where the horizontal axis is expressed as the multiple of the interference period.

Download Full Size | PDF

In order to verify the proposed method, the different numbers of reconstructed images are used in the averaging process. Figures 4(a) and 4(b) represents the amplitude images from a single hologram and after 15 averaging processes, respectively. Figures 4(c) and 4(d) represent the corresponding phase images before and after the averaging process. The reduction of coherent noise is clearly visible both in the amplitude and phase images. For quantitatively comparing the effect of coherent noise reduction, the spatial standard deviation (σ) of the noise pattern can be evaluated by

σ=n,m[ϕ(n,m)ϕ¯]2/(n×m1).
where ϕ denotes the value distribution of the average image, ϕ¯ denotes the mean values, n and m represent the pixel number of row and column of calculation region, which is indicated with a white rectangle in Fig. 3. Figure 4(e) shows the value of σ as a function of the number of averaged images. As expected, the σ of the noise pattern in the amplitude and phase images rapidly decrease with the number of averaged images. For intensity imaging, the amplitude noise is reduced from 3.3 × 10−3 to 0.6 × 10−3 and the image quality is improved significantly. For phase measurement, the phase noise is suppressed from 11° to 2°, which corresponds to an equivalent height from 16 nm to 3 nm. Figures 5(a) and 5(b) represent the 3D surface profile of the resolution targets by using a single hologram and 15 averaging processes, respectively. Figure 5(c) shows the profiles along the straight line indicated by a white line in Fig. 5(a). These results demonstrate that the proposed method can significantly reduce the coherent noise both in amplitude and phase image.

 figure: Fig. 4

Fig. 4 (a) and (b): the amplitude images from a single hologram and after 15 averaging process, (c) and (d): the phase images from a single hologram and after 15 averaging process, (e): the standard deviation in the amplitude and phase images over the number of reconstructed images for averaging process.

Download Full Size | PDF

 figure: Fig. 5

Fig. 5 (a): 3D phase image from a single hologram, (b): phase image after 15 averaging processes, (c) phase distributions along a line before and after the averaging process.

Download Full Size | PDF

As an example of application, a microlens (the diameter of 90 μm and the height of 21 μm) is measured by using the proposed approach. Figures 6(a) and 6(b) present the amplitude images reconstructed from a single hologram and after averaging 15 holograms together, respectively. The amplitude noise can be seen in Fig. 6(a) and is clearly reduced in Fig. 6(b). Figures 6(c) and 6(d) show the corresponding wrapped phase images. In comparison, the phase noise in Fig. 6(c) appear significantly reduced by the averaging process. To demonstrate topography measurement, the central region of the reconstructed phase image, which is indicated by a rectangular region in Fig. 6(c), are presented in Figs. 7(a) and 7(b) with 3D perspective. And the corresponding unwrapped phase images are obtained by using a least-squares algorithm [37] and are shown in Figs. 7(c) and 7(d). Furthermore, the Figs. 7(e) and 7(f) present the compensated phase images obtained by application of a digital phase mask, which is defined by fitting the unwrapped phase data with the Zernike parameters [38]. It can be noticed that some noise points that appear in Fig. 7(e) are suppressed in Fig. 7(f). Figure 7(g) presents the phase profiles along a straight line indicated by a white line in Fig. (e) before and after the averaging process. The reduction of phase noise can clearly be observed from these phase pictures, and the accuracy of shape measurement can be significantly improved.

 figure: Fig. 6

Fig. 6 The measurement results of a microlens. (a) and (b): the amplitude image from a single hologram and after 15 averaging process, (c) and (d): the phase image before and after the averaging process.

Download Full Size | PDF

 figure: Fig. 7

Fig. 7 (a) and (b): 3D wrapped phase images before and after averaging process, (c) and (d): unwrapped phase images corresponding (a) and (b), (e) and (f): compensated phase images corresponding (c) and (d), (g): phase profiles along a straight line indicated by a white line in (e).

Download Full Size | PDF

4. Conclusion

This paper presents a new approach to reduce coherent noise in DHM by averaging multiple holograms recorded with a multimode laser. First, based on the periodicity of temporal coherence of multimode semiconductor laser, a series of holograms are generated and recorded by changing the optical path difference between reference and object beams. In this step, due to the shortly coherent light, the parasitic interference fringes caused by the multiple reflections or transmissions are prevented from corrupting holograms. Meanwhile, a set of object reconstructed image is obtained, which occupied different noise patterns. Afterward, the average processing can effectively suppress the speckle pattern, which comes from the optical scattering on object surface or interior. Finally, the image quality and the measurement accuracy are improved significantly. Experimental results show that the proposed method is superior to reduce the coherent noise in intensity image and suppress the phase noise in phase image. As we all know, if a set of totally uncorrelated noise patterns are added for averaging process, the standard deviation of noise pattern falls in proportion to 1/N, where N is the number of holograms used to averaging process [39]. According to the experimental results of the proposed method, the values of standard deviation in amplitude and phase images decrease considerably for the first few holograms but at a lesser extent for higher number of holograms. It is easy to see that the noise reduction trend is bounded by the ideal 1/N, as well as other multiple holograms approaches. The limit of denoising performance can be interpreted in terms of the partial correlation between the noise patterns. Although the multiple holograms methods are different from each other in noise decorrelation extents, the value of cross correlations is typically in the range of 01 to 0.2 [15–18]. Compared with other multiple hologram approaches, the proposed method can effectively suppress the parasitic fringes and speckle pattern in both intensity and phase images for DHM. Meanwhile the acquisition of multiple holograms is merely dependent on the linear motion of an optical path retarder, and the extra image registration is not required before averaging process. Although the coherent noise reduction is achieved at cost of time for capturing multiple holograms, the result is fully dependent on optical process without any loss of resolution of imaging. The proposed method is effective and feasible for inspecting the shape and surface quality of the microlens with high precision.

Funding

National Major Scientific Instruments Development Project of China (No. 2013YQ030595).

References and links

1. M. K. Kim, Digital holographic Microscopy: Principles, Techniques, and Applications (Springer, 2011).

2. P. Picart, New Techniques in Digital Holography (ISTE-Wiley, 2015).

3. T. Kreis, Handbook of Holographic Interferometry: Optical and Digital Methods (Wiley, 2005).

4. J. W. Goodman, Optical Methods for Suppressing Speckle, in Speckle Phenomena in Optics (Roberts & Company, 2007).

5. J. Garcia-Sucerquia, “Noise reduction in digital lensless holographic microscopy by engineering the light from a light-emitting diode,” Appl. Opt. 52(1), A232–A239 (2013). [CrossRef]   [PubMed]  

6. R. L. Guo, B. L. Yao, J. W. Min, M. L. Zhou, X. H. Yu, M. Lei, S. H. Yan, Y. L. Yang, and D. Dan, “LED-based digital holographic microscopy with slightly off-axis interferometry,” J. Opt. 125408, 1–8 (2014).

7. P. Langehanenberg, G. V. Bally, and B. Kemper, “Application of partially coherent light in live cell imaging with digital holographic microscopy,” J. Mod. Opt. 57(9), 709–717 (2010). [CrossRef]  

8. F. Dubois, L. Joannes, and J. C. Legros, “Improved three-dimensional imaging with a digital holography microscope with a source of partial spatial coherence,” Appl. Opt. 38(34), 7085–7094 (1999). [CrossRef]   [PubMed]  

9. Y. Choi, T. D. Yang, K. J. Lee, and W. Choi, “Full-field and single-shot quantitative phase microscopy using dynamic speckle illumination,” Opt. Lett. 36(13), 2465–2467 (2011). [CrossRef]   [PubMed]  

10. J. Dohet-Eraly, C. Yourassowsky, A. E. Mallahi, and F. Dubois, “Quantitative assessment of noise reduction with partial spatial coherence illumination in digital holographic microscopy,” Opt. Lett. 41(1), 111–114 (2016). [CrossRef]   [PubMed]  

11. B. R. A. Michael and C. H. Cao, “Speckle-free laser imaging using random laser illumination,” Nat. Photonics 6, 355–359 (2012).

12. Z. Wang, L. Millet, M. Mir, H. Ding, S. Unarunotai, J. Rogers, M. U. Gillette, and G. Popescu, “Spatial light interference microscopy (SLIM),” Opt. Express 19(2), 1016–1026 (2011). [CrossRef]   [PubMed]  

13. S. Shin, K. Kim, K. Lee, S. Lee, and Y. Park, “Effects of spatiotemporal coherence on interferometric microscopy,” Opt. Express 25(7), 8085–8097 (2017). [CrossRef]   [PubMed]  

14. T. Nomura, M. Okamura, E. Nitanai, and T. Numata, “Image quality improvement of digital holography by superposition of reconstructed images obtained by multiple wavelengths,” Appl. Opt. 47(19), D38–D43 (2008). [CrossRef]   [PubMed]  

15. L. Rong, W. Xiao, F. Pan, S. Liu, and R. Li, “Speckle noise reduction in digital holography by use of multiple polarization holograms,” Chin. Opt. Lett. 8(7), 653–655 (2010). [CrossRef]  

16. W. Xiao, J. Zhang, L. Rong, F. Pan, S. Liu, F. Wang, and A. He, “Improvement of speckle noise suppression in digital holography by rotating linear polarization state,” Chin. Opt. Lett. 9(6), 060901 (2011). [CrossRef]  

17. F. Pan, W. Xiao, S. Liu, and L. Rong, “Coherent noise reduction in digital holographic microscopy by laterally shifting camera,” Opt. Commun. 292, 68–72 (2013). [CrossRef]  

18. F. Pan, W. Xiao, S. Liu, F. Wang, L. Rong, and R. Li, “Coherent noise reduction in digital holographic phase contrast microscopy by slightly shifting object,” Opt. Express 19(5), 3862–3869 (2011). [CrossRef]   [PubMed]  

19. C. G. Quan, X. Kang, and C. J. Tay, “Speckle noise reduction in digital holography by multiple holograms,” Opt. Eng. 46(11), 115801 (2007). [CrossRef]  

20. Y. Wang, P. Meng, D. Wang, L. Rong, and S. Panezai, “Speckle noise suppression in digital holography by angular diversity with phase-only spatial light modulator,” Opt. Express 21(17), 19568–19578 (2013). [CrossRef]   [PubMed]  

21. Y. Park, W. Choi, Z. Yaqoob, R. Dasari, K. Badizadegan, and M. S. Feld, “Speckle-field digital holographic microscopy,” Opt. Express 17(15), 12285–12292 (2009). [CrossRef]   [PubMed]  

22. S. Kubota and J. W. Goodman, “Very efficient speckle contrast reduction realized by moving diffuser device,” Appl. Opt. 49(23), 4385–4391 (2010). [CrossRef]   [PubMed]  

23. J. Leng, X. Z. Sang, and B. B. Yan, “Speckle noise reduction in digital holography with spatial light modulator and nonlocal means algorithm,” Chin. Opt. Lett. 12(4), 040301 (2014). [CrossRef]  

24. V. Bianco, P. Memmolo, M. Paturzo, and P. Ferraro, “On-speckle suppression in IR digital holography,” Opt. Lett. 41(22), 5226–5229 (2016). [CrossRef]   [PubMed]  

25. P. Feng, X. Wen, and R. Lu, “Long-working-distance synthetic aperture Fresnel off-axis digital holography,” Opt. Express 17(7), 5473–5480 (2009). [CrossRef]   [PubMed]  

26. S. Montresor and P. Picart, “Quantitative appraisal for noise reduction in digital holographic phase imaging,” Opt. Express 24(13), 14322–14343 (2016). [CrossRef]   [PubMed]  

27. J. G. Garcia-Sucerquia, J. A. H. Ramirez, and D. V. Prieto, “Reduction of speckle noise in digital holography by using digital image processing,” Optik (Stuttg.) 116(1), 44–48 (2005). [CrossRef]  

28. X. O. Cai, “Reduction of speckle noise in the reconstructed image of digital holography,” Optik (Stuttg.) 121(4), 394–399 (2010). [CrossRef]  

29. J. Maycock, B. M. Hennelly, J. B. McDonald, Y. Frauel, A. Castro, B. Javidi, and T. J. Naughton, “Reduction of speckle in digital holography by discrete Fourier filtering,” J. Opt. Soc. Am. A 24(6), 1617–1622 (2007). [CrossRef]   [PubMed]  

30. A. Sharma, G. Sheoran, Z. A. Jaffery, and Moinuddin, “Improvement of signal-to-noise ratio in digital holography using wavelet transform,” Opt. Lasers Eng. 46(1), 42–47 (2008). [CrossRef]  

31. A. Uzan, Y. Rivenson, and A. Stern, “Speckle denoising in digital holography by nonlocal means filtering,” Appl. Opt. 52(1), A195–A200 (2013). [CrossRef]   [PubMed]  

32. M. Leo, R. Piccolo, C. Distante, P. Memmolo, M. Paturzo, and P. Ferraro, “Multilevel bidimensional empirical mode decomposition: a new speckle reduction method in digital holography,” Opt. Eng. 53(11), 112314 (2014). [CrossRef]  

33. V. Bianco, M. Paturzo, P. Memmolo, A. Finizio, P. Ferraro, and B. Javidi, “Random resampling masks: a non-Bayesian one-shot strategy for noise reduction in digital holography,” Opt. Lett. 38(5), 619–621 (2013). [CrossRef]   [PubMed]  

34. D. Hincapie, J. Herrera-Ramírez, and J. Garcia-Sucerquia, “Single-shot speckle reduction in numerical reconstruction of digitally recorded holograms,” Opt. Lett. 40(8), 1623–1626 (2015). [CrossRef]   [PubMed]  

35. M. Haouat, J. Garcia-Sucerquia, A. Kellou, and P. Picart, “Reduction of speckle noise in holographic images using spatial jittering in numerical reconstructions,” Opt. Lett. 42(6), 1047–1050 (2017). [CrossRef]   [PubMed]  

36. V. Bianco, P. Memmolo, M. Paturzo, A. Finizio, B. Javidi, and P. Ferraro, “Quasi noise-free digital holography,” Light Sci. Appl. 5(9), 1–12 (2016). [CrossRef]  

37. D. C. Ghiglia and M. D. Pritt, Two-Dimensional Phase Unwrapping: Theory, Algorithms, and Software (Wiley, 1998).

38. T. Colomb, F. Montfort, J. Kühn, N. Aspert, E. Cuche, A. Marian, F. Charrière, S. Bourquin, P. Marquet, and C. Depeursinge, “Numerical parametric lens for shifting, magnification, and complete aberration compensation in digital holographic microscopy,” J. Opt. Soc. Am. A 23(12), 3177–3190 (2006). [CrossRef]   [PubMed]  

39. J. C. Dainty, Laser Speckle and Related Phenomena (2nd ed.) (Springer-Verlag, 1984).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (7)

Fig. 1
Fig. 1 The fringe visibility versus the optical path different between the R and O beams. (a): the digital hologram, (b): the calculation results, (c) (d) and (e): the experimental results, the cutouts (indicated with a white rectangle) of the holograms recorded with different (LO-LR).
Fig. 2
Fig. 2 Schematic of the experimental system. Laser, a multimode semiconductor laser; λ/2, half-wave plate; PBS, polarizing beam splitter; BE, beam expander with spatial filter; M, mirror; L, lens; O, object wave; R, reference wave; MO, microscope objective; BS, beam splitter; CMOS, camera; OPR, optical path retarder to adjust the optical path length of the reference beam.
Fig. 3
Fig. 3 The experimental results of the noise decorrelation by changing the LR and keeping LO. (a) and (b): the amplitude images with a reference LR and a shifted LR, (c) and (d): the phase images with a reference LR and a shifted LR, (e): the normalized cross correlations (γ) of the noise patterns over the LR-LO, where the horizontal axis is expressed as the multiple of the interference period.
Fig. 4
Fig. 4 (a) and (b): the amplitude images from a single hologram and after 15 averaging process, (c) and (d): the phase images from a single hologram and after 15 averaging process, (e): the standard deviation in the amplitude and phase images over the number of reconstructed images for averaging process.
Fig. 5
Fig. 5 (a): 3D phase image from a single hologram, (b): phase image after 15 averaging processes, (c) phase distributions along a line before and after the averaging process.
Fig. 6
Fig. 6 The measurement results of a microlens. (a) and (b): the amplitude image from a single hologram and after 15 averaging process, (c) and (d): the phase image before and after the averaging process.
Fig. 7
Fig. 7 (a) and (b): 3D wrapped phase images before and after averaging process, (c) and (d): unwrapped phase images corresponding (a) and (b), (e) and (f): compensated phase images corresponding (c) and (d), (g): phase profiles along a straight line indicated by a white line in (e).

Equations (10)

Equations on this page are rendered with MathJax. Learn more.

I H = | R | 2 + | O | 2 + R O + R O .
E = n = ( N 1 ) / 2 n = ( N 1 ) / 2 A n exp { i [ ( k + n Δ k ) L + φ n ] } .
R = R 0 n = ( N 1 ) / 2 n = ( N 1 ) / 2 A n exp { i [ ( k + n Δ k ) L R + φ n ] } . O = O 0 n = ( N 1 ) / 2 n = ( N 1 ) / 2 A n exp { i [ ( k + n Δ k ) L O + φ n ] } .
R * O = R 0 O 0 n = ( N 1 ) / 2 n = ( N 1 ) / 2 m = ( N 1 ) / 2 m = ( N 1 ) / 2 A n A m exp [ i k ( L O L R ) ] × exp [ i Δ k ( m L O n L R ) ] exp [ i ( φ m φ n ) ] .
R * O = R 0 O 0 exp [ i k ( L O L R ) ] n = ( N 1 ) / 2 n = ( N 1 ) / 2 m = ( N 1 ) / 2 m = ( N 1 ) / 2 δ n , m A n A m × exp [ i Δ k ( n L O m L R ) ] exp [ i ( φ m φ n ) ] .
δ n , m = { 1 ( n = m ) 0 ( n m )
R * O = R 0 O 0 exp [ i k ( L O L R ) ] n = ( N 1 ) / 2 n = ( N 1 ) / 2 A n 2 exp [ i Δ k n ( L O L R ) ] .
R * O = R 0 * O 0 exp [ i k ( L O L R ) ] { sin [ N Δ k ( L O L R ) / 2 ] sin [ Δ k ( L O L R ) / 2 ] } .
γ = n , m [ ϕ 0 ( n , m ) ϕ ¯ 0 ] [ ϕ L ( n , m ) ϕ ¯ L ] n , m [ ϕ 0 ( n , m ) ϕ ¯ 0 ] 2 n , m [ ϕ L ( n , m ) ϕ ¯ L ] 2 .
σ = n , m [ ϕ ( n , m ) ϕ ¯ ] 2 / ( n × m 1 ) .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.