Photonics Research, 2014, 2 (2): 02000064, Published Online: Nov. 5, 2014   

Enhanced depth resolution in optical scanning holography using a configurable pupil Download: 1045次

Author Affiliations
1 Department of Electrical and Electronic Engineering, University of Hong Kong, Pokfulam, Hong Kong, China
2 Institute of Applied Physics, University of Electronic Science and Technology of China, 610054, Chengdu, China
3 Bradley Department of Electrical and Computer Engineering, Virginia Tech, Blacksburg, Virginia 24061, USA
Abstract
The optical scanning holography (OSH) technique can capture all the three-dimensional volume information of an object in a hologram via a single raster scan. The digital hologram can then be processed to reconstruct individual sectional images of the object. In this paper, we present a scheme to reconstruct sectional images in OSH with enhanced depth resolution, where a spatial light modulator (SLM) is adopted as a configurable point pupil. By switching the SLM between two states, different Fresnel zone plates (FZPs) are generated based on the same optical system. With extra information provided by different FZPs, a depth resolution at 0.7 μm can be achieved.

1. INTRODUCTION

Digital holography is a technique that seeks to capture the three-dimensional (3D) volume information of an object and to record it on a two-dimensional (2D) hologram electronically [1,2]. As one of the digital holographic techniques, optical scanning holography (OSH) can record the 2D holographic information of the 3D object by a single raster scan [3]. The idea of OSH was generated by Poon et al. when he was dealing with bipolar point spread functions in incoherent image processing [4,5]. The theory of OSH was then developed, and the first experiment was carried out in 1990 [6]. Since then, OSH has undergone great development and has found various applications ranging from biological microscopy [7,8] to remote sensing [9,10] to 3D cryptography [11], etc.

Unlike laser scanning [12,13] or confocal microscopy techniques [14,15], where multiple 2D scans are needed to capture all the 3D information, OSH can greatly save on data acquisition time. While it is so powerful in capturing the 3D information, there is also a need in OSH to view the individual 2D planes of the object from the acquired digital hologram, which is known as sectioning or sectional image reconstruction. In the sectioning process, the challenge lies in the suppressing of the defocus noise, which is the undesired residue signal from neighboring sections. The conventional method of achieving this involves computing the convolution of the conjugate impulse response at the desired section with the hologram, but this approach suffers from large defocus noise [3]. To reduce the defocus noise, many methods including the Wiener filter [16], Wigner distribution [17], and inverse imaging [18,19] have been demonstrated. Further development related to inverse imaging includes the development of a blind edge detection technique to locate the sections [20], compressed sensing [21], and edge-preserving regularization for sharper sections [22,23]. The use of a random-phase pupil is also shown to achieve better sectioning resolution [24], which requires the averaging of several independent holograms of the same object. Another method using a dual-wavelength laser source has been developed, with the depth resolution improved up to 2.5 μm [25,26]. More recently, Ou et al. have proposed a double-detection method, which has led to a depth resolution of 1 μm by capturing holographic information of the same object at two different depth locations [27].

In this paper, we propose a sectioning method based on a configurable point pupil for achieving enhanced depth resolution. A spatial light modulator (SLM) is used as the switchable pupil, which generates different Fresnel zone plates (FZPs) for the same OSH system. We can then make use of the extra information provided by the different FZPs to achieve a better depth resolution. The proposed method is easier to realize in practice compared to some of the methods mentioned above.

2. PRINCIPLE

The OSH system setup is shown in Fig. 1. In the optical system, the laser source centered at ω is divided into two parts by a beam splitter (BS1), of which one passes through a pupil function p1(x,y) to have a spherical wavefront on the object. Meanwhile, the other path first has a frequency shift Ω via an acousto-optic frequency shifter (AOFS) and then arrives at a point pupil p2(x,y) to result in a planar wavefront on the object. The two coherent beams of different frequencies ω and ω+Ω are then combined by a second beamsplitter (BS2) and are used to scan the object located at a distance z away from the scanning mirror. Lens 3 is used to collect the light from the object. A photodetector collects the transmitted and scattered light from the object and converts it into an electronic signal. To recover the complex scattered wavefront from the object, an electronic demodulator and a low-pass filter are used to process the measured holograms, the details of which can be found in [28].

Fig. 1. Schematic of OSH.

下载图片 查看所有图片

The optical transfer function (OTF) of the OSH system can be expressed as [28] where x,y,z are the spatial coordinates, kx and ky are the corresponding lateral spatial frequency coordinates, f is the focal length of lens 1 and lens 2, and k0=2π/λ is the wavenumber of the optical source, with λ representing the wavelength of the optical source.

The description and the imaging equation so far are similar to the conventional OSH setup. However, our scheme actually differs in two important aspects. First, we make use of an SLM to realize the pupil function p2(x,y), making it much more flexible and programmable. For our purpose here, we switch the SLM to two configurations, denoted as p2a(x,y) and p2b(x,y). In principle, it is also possible to design the SLM for more pupil functions, which may be combined with further computation for better sectioning capability. Second, the imaging process is also different, in that we take one measurement for each pupil configuration. Specifically, we first set p2a=δ(x,y) and record the first hologram g1(x,y); then, we tune the SLM to a second configuration and obtain the other hologram g2(x,y). The requirement for the second configuration is that p2b(x,y) is sufficiently different from p2a(x,y) to maximize the information content from the two measurements. In what follows, we will show that even with a simple shift of the impulse such that p2b(x,y)=δ(xx0,yy0), where x0 and y0 are certain nonzero constants, we can increase the depth resolution of the OSH system significantly.

For the first hologram, the impulse response of the OSH system [with p1(x,y)=1 and p2a(x,y)=δ(x,y)] can be derived from Eq. (1) as where the impulse response h(x,y;z) is the Fourier transform of the OTF H(kx,ky;z). Given an input object with an intensity distribution φ(x,y,z), the resulting hologram is therefore where * denotes the 2D convolution operation and zi represents the depth location of the ith section, where we have assumed that we can discretize the object into n sections.

For the second measurement, with p1(x,y)=1 and p2b(x,y)=δ(xx0,yy0), the shifted OTF can be simplified as the impulse response thus can be expressed as and the corresponding hologram is

We can combine the two hologram recordings into a single imaging equation as follows. Let g1 and g2 denote the lexicographical ordering of the holograms g1(x,y) and g2(x,y), respectively. In a similar way, the object at depth i becomes φi. Its convolution with the impulse response at zi is denoted by the matrix–vector product H1(zi)φi for the first hologram and H2(zi)φi for the second one. The overall imaging is then represented by where ϵ1 to ϵn represent the respective observation noise values. H1 and H2 denote the matrices of the first and second measurement, respectively, and H encapsulates all the impulse responses, while the object at various sections is represented by φ.

Our objective in sectioning is to recover the object vector φ from the observation g. This process can be viewed as a minimization problem [18]: where ·2 denotes the 2 norm, λ>0 is a penalty parameter, and C stands for the Laplacian operator. Other regularization methods, such as 1 norm or total variation norm, have shown to have better results [23,29,30]; however, we choose 2-norm regularization here to have a better comparison with both the inverse imaging method [19] and the double-detection technique [27].

The solution to this problem can be obtained as [26] where H+ indicates the conjugate transpose of H. We use the conjugate gradient method to solve this problem, the details of which are given in [31].

3. SIMULATION AND ANALYSIS

In this section, we analyze the sectioning results of the proposed method with an object having two and three dominant sections. We also make a comparison of the proposed method with double detection and study the effect of measurement error.

3.2 A. Sectioning with Two Sections

The proposed method is first demonstrated with an object that contains two sections at z1=34mm and z2=34.0007mm, with the section distance equal to 0.7 μm. Each section has a size of 1mm×1mm in the lateral directions and is sampled to 512×512 pixels, as shown in Figs. 2(a) and 2(b). The object is illuminated by a He–Ne laser with the wavelength centered at 632 nm. The configurable pupil p2(x,y) is realized by a liquid-crystal-on-silicon (LCoS) transmissive SLM with an 8 μm pixel size, which can switch between two configurations: (1) p2a(x,y)=δ(x,y) and (2) p2b(x,y)=δ(xx0,yy0), with x0=0.48mm and y0=0mm.

Fig. 2. Two sections together with the real part of the FZPs for each scan.

下载图片 查看所有图片

In the simulation, we use an OSH simulator to evaluate the performance of the proposed method, in which the finite diameter of the apertures and different kinds of noise are taken into consideration. The details of the OSH simulator can be found in [26]. Here, in the proposed method, the diameter of the collimated beam is D=40mm and the focal length is f=50mm, which lead to a system numerical aperture [NA, where NA=sin(D/2f)] equal to 0.3894. Both the thermal noise and the Poisson noise are evaluated in the simulator: the former is a white Gaussian noise with an energy per bandwidth at σ02 (σ0=12.75), while the other has a Gaussian distribution N(p¯,cp¯2), where p¯ denotes the mean value of the detected hologram signal amplitude and c=0.0025.

The measured FZPs with each configuration for the two sections are shown in Fig. 2. One can observe that as the pupil switches to p2b(x,y)=δ(xx0,yy0), the center of the FZP shifts 0.54 mm along the x axis as well. As a result, the opaque and transparent zones are different in the observation area, which indicates that the spatial sampling rate of the FZPs will change as they shift away from the center. Accordingly, extra information can be captured during the second measurement, with which we can expect to allow for better sectioning. Convolving these FZPs with the object results in the holograms shown in Fig. 3.

Fig. 3. Holograms containing two-sectional images of the object with different pupil functions (a), (b) for x0=y0=0 and (c), (d) for x0=0.48mm, y0=0.

下载图片 查看所有图片

These holograms become the input to the sectional image reconstruction described in the previous section. For comparison, we have implemented three different methods for sectioning: (1) the conventional method, (2) the conjugate gradient method with a single measurement, and (3) the proposed method with two measurements. The results are shown in Fig. 4. It can be observed that neither the conventional method nor the conjugate gradient method with a single measurement provides acceptable sectioning results to distinguish the different sections. On the other hand, with the proposed method, thanks to the extra information captured during the second measurement, each section is clearly separated with a small amount of defocus noise. This indicates that the depth resolution can reach 0.7 μm in this case. The measured signal-to-noise rate (SNR) is 38.5 dB in the simulation.

Fig. 4. Sectioning results using (a), (b) the conventional method; (c), (d) the conjugate-gradient-based method with a single measurement; (e), (f) the conjugate-gradient-based method with the proposed method.

下载图片 查看所有图片

The extra information captured during the second measurement is highly related to the correlation between the matrices H1 and H2 according to Eq. (7). The correlation of H1 and H2 can be evaluated using the following equation [32]: N=x=1Ny=1Mi=1Nj=1MH1(x+i,y+j)·H2(i,j)i=1Nj=1MH1(x+i,y+j)2·i=1Nj=1MH2(i,j)2.

As the difference between H1 and H2 is mainly caused by the lateral shift x0, it is therefore useful to study the relationship between x0 and the correlation of H1 and H2, as well as the relationship between x0 and the system resolution. The results are shown in Fig. 5. There are two important features one can observe from this plot: (1) the correlation between the two matrices H1 and H2 decreases as x0 increases, which is in accordance with our expectation, as a larger x0 would lead to less similarity between H1 and H2; (2) as x0 increases, the depth resolution would improve from 7.5 to 0.7 μm. This is because as the center of the FZP shifts away, the transition between the opaque and transparent zones becomes faster, which indicates that the spatial sampling rate of the FZPs becomes higher. This would provide us with more extra information for the hologram, which results in an enhanced depth resolution.

Fig. 5. Relationship between x0 and the resolution and the relationship between x0 and the correlation of H1 and H2.

下载图片 查看所有图片

3.3 B. Comparison with the Double-Detection Method

From the analysis in the previous section, we can see that there is a trade-off between the lateral shift x0 and the depth resolution. This is also the case for the double-detection method, where the trade-off exists between the shift in depth (z axis) and the sectioning resolution [27]. In order to evaluate the performance of the proposed method and the double-detection method, it is necessary to make a comparison between the two.

We first compare the relationship between the spatial shift and the correlation of H1 and H2, as is shown in Fig. 6. It can be seen that with the same amount of shift, the lateral displacement with the proposed method would introduce more difference (i.e., less similarity and smaller correlation) between the two matrices representing the two OTFs compared to the double-detection method with displacement along the z axis.

Fig. 6. Relationship between spatial shift and the correlation of H1 and H2.

下载图片 查看所有图片

We have also compared the relationship between the spatial shift and the resolution. The result is shown in Fig. 7. One can observe that the two curves hold the same trend, except that the one with the proposed method changes much faster, which is in accordance with the analysis of the correlation above. To achieve a depth resolution up to 1 μm, one needs to move the object along the z axis for about 15 to 20 mm by using the double-detection technique, while with the proposed method, the lateral shift of the pupil p2 is only around 0.1 mm. It is worth mentioning that due to the SLM used in the system, there is no need to move the pupil p2 itself in the actual optical system; one only needs to adjust the power distribution applied on the SLM along the x axis for several pixels, which is more time efficient.

Fig. 7. Relationship between spatial shift and resolution.

下载图片 查看所有图片

3.4 C. Sectioning with Three Sections

We then evaluate the performance of the proposed method with an object with three sections, as is shown in Fig. 8. The location of the three sections are z1=34mm, z2=34.0007mm, and z3=34.0014mm, respectively. The OSH simulator is also used to perform the simulation, in which the system NA and the noise are kept the same as before. The pupil function also switches between two states: (1) p2a(x,y)=δ(x,y) and (2) p2b(x,y)=δ(xx0,yy0), with x0=0.48mm and y0=0mm. The holograms from each scan of the object are shown in Fig. 9.

Fig. 8. Object with three sections.

下载图片 查看所有图片

Fig. 9. Holograms containing three-section images of the object with different pupil functions (a), (b) for x0=y0=0 and (c), (d) for x0=0.48mm, y0=0.

下载图片 查看所有图片

The sectioning results for this three-section object are shown in Fig. 10, in which the sectioning effect of the traditional method, the conjugate gradient method with a single measurement, and our proposed method with two measurements are compared. The proposed method also outperforms the other two in suppressing the defocus noise, which indicates that it can handle objects of multiple sections as well.

Fig. 10. Sectioning results using (a)–(c) the conventional method; (d)–(f) the conjugate-gradient-based method with a single measurement; (g)–(i) the conjugate-gradient-based method with the proposed method.

下载图片 查看所有图片

3.5 D. Measurement Error Analysis

Because we need to switch the SLM between two states, there would be inevitable measurement error in practice. However, unlike the double-detection method in [27], where the object itself is displaced along the z axis, with configurable pupils there is no need to move any part of the optical setup, which should give rise to a smaller amount of the measurement error. The error mostly comes from the manufacturing defect of the pixel size of the SLM.

We use the same two-section object to test the effect of measurement error, with x0=0.48mm for the second measurement. Here we define Δx as the measurement error along the x axis, so the pupil function p2b(x,y) actually becomes p2b(x,y)=δ(xx0,yy0), with x0=x0+Δx and y0=0mm. Three different situations are considered in the simulation: (1) Δx=0mm, z1=34mm, z2=34.0007mm; (2) Δx=100nm, z1=34mm, z2=34.0007mm; (3) Δx=100nm, z1=34mm, z2=34.01mm. The corresponding sectioning results using the proposed method are shown in Fig. 11. One can see that it is hard to distinguish the two sections when the measurement error reaches 100 nm. However, as the section distance increases up to 10 μm, this measurement error would be tolerated, leaving only a small amount of defocus noise. This indicates that the depth resolution has been degraded from 0.7 to around 10 μm when measurement error Δx=100nm.

Fig. 11. Sectioning results using the proposed method with (a)–(b) Δx=0nm, z1=34mm, z2=34.0007mm; (c)–(d) Δx=100nm, z1=34mm, z2=34.0007mm; (e)–(f) Δx=100nm, z1=34mm, z2=34.01mm.

下载图片 查看所有图片

We have also measured the relationship between the measurement error of the proposed method and the system depth resolution, which is shown in Fig. 12. As one would expect, the depth resolution would degrade as the measurement error increases. The depth resolution would degrade to 10 μm when the measurement error reaches 100 nm. For comparison, we also show the relationship between the measurement error and the system depth resolution of the double-detection method here. One can see that the proposed method with a configurable pupil is more sensitive to measurement error compared to the double-detection method. This is due to the fact that the matrix H changes more quickly with the lateral spatial shift than the z axis spatial shift, as discussed before.

Fig. 12. Relationship between measurement error and the depth resolution.

下载图片 查看所有图片

4. CONCLUSIONS

A method has been proposed to reconstruct a sectional image using a configurable point pupil, which can be realized with a commercially available LCoS SLM. By controlling the voltage applied on each pixel of the SLM, the point pupil can be switched between two configurations. This can create two different sets of FZPs, with which extra information about the object can be extracted. By combining the holograms recorded during each scan process, we have shown that better sectioning results are possible.

References

[1] J. W. Goodman, R. W. Lawrence. Digital image formation from electronically detected holograms. Appl. Phys. Lett., 1967, 11: 77-79.

[2] GoodmanJ. W., “Digital image formation from holograms: early motivations and modern capabilities,” in OSA Topical Meeting in Digital Holography and Three-Dimensional Imaging (2007), paper JMA1.

[3] PoonT.-C., Optical Scanning Holography with MATLAB (Springer, 2007).

[4] T.-C. Poon, A. Korpel. Optical transfer function of an acousto-optic heterodyning image processor. Opt. Lett., 1979, 4: 317-319.

[5] T.-C. Poon. Scanning holography and two-dimensional image processing by acousto-optic two-pupil synthesis. J. Opt. Soc. Am. A, 1985, 2: 521-527.

[6] T.-C. Poon, B. D. Duncan, M. H. Wu, K. Shinoda, Y. Suzuki. Real-time optical holography using a spatial light modulator. Jpn. J. Appl. Phys., 1990, 29: L1840-L1842.

[7] J. Swoger, M. Martínez-Corral, J. Huisken, E. H. K. Stelzer. Optical scanning holography as a technique for high-resolution three-dimensional biological microscopy. J. Opt. Soc. Am. A, 2002, 19: 1910-1918.

[8] G. Indebetouw, W. Zhong. Scanning holographic microscopy of three-dimensional fluorescent specimens. J. Opt. Soc. Am. A, 2006, 23: 1699-1707.

[9] B. W. Schilling, G. C. Templeton. Three-dimensional remote sensing by optical scanning holography. Appl. Opt., 2001, 40: 5474-5481.

[10] T. Kim, T.-C. Poon, G. Indebetouw. Depth detection and image recovery in remote sensing by optical scanning holography. Opt. Eng., 2002, 41: 1331-1338.

[11] T.-C. Poon, T. Kim, K. Doh. Optical scanning cryptography for secure wireless transmission. Appl. Opt., 2003, 42: 6496-6503.

[12] W. Denk, J. Strickler, W. Webb. Two-photon laser scanning fluorescence microscopy. Nature, 1990, 248: 73-76.

[13] J. S. Ploem. Laser scanning fluorescence microscopy. Appl. Opt., 1987, 26: 3226-3231.

[14] S. Nie, D. Chiu, R. Zare. Probing individual molecules with confocal fluorescence microscopy. Science, 1994, 266: 1018-1021.

[15] J. Korlach, P. Schwille, W. Webb, G. Feigenson. Characterization of lipid bilayer phases by confocal microscopy and fluorescence correlation spectroscopy. Proc. Natl. Acad. Sci. USA, 1999, 96: 8461-8466.

[16] T. Kim. Optical sectioning by optical scanning holography and a Wiener filter. Appl. Opt., 2006, 45: 872-879.

[17] H. Kim, S.-W. Min, B. Lee, T.-C. Poon. Optical sectioning for optical scanning holography using phase-space filtering with Wigner distribution functions. Appl. Opt., 2008, 47: D164-D175.

[18] E. Y. Lam, X. Zhang, H. Vo, T.-C. Poon, G. Indebetouw. Three-dimensional microscopy and sectional image reconstruction using optical scanning holography. Appl. Opt., 2009, 48: H113-H119.

[19] X. Zhang, E. Y. Lam, T.-C. Poon. Reconstruction of sectional images in holography using inverse imaging. Opt. Express, 2008, 16: 17215-17226.

[20] X. Zhang, E. Y. Lam, T. Kim, Y. S. Kim, T.-C. Poon. Blind sectional image reconstruction for optical scanning holography. Opt. Lett., 2009, 34: 3098-3100.

[21] ZhangX.LamE. Y., “Sectional image reconstruction in optical scanning holography using compressed sensing,” in IEEE International Conference on Image Processing (2010), pp. 33493352.

[22] ZhangX.LamE. Y., “Edge detection of three-dimensional object by manipulating pupil functions in optical scanning holography system,” in IEEE International Conference on Image Processing (2010), pp. 36613664.

[23] X. Zhang, E. Y. Lam. Edge-preserving sectional image reconstruction in optical scanning holography. J. Opt. Soc. Am. A, 2010, 27: 1630-1637.

[24] Z. Xin, K. Dobson, Y. Shinoda, T.-C. Poon. Sectional image reconstruction in optical scanning holography using a random-phase pupil. Opt. Lett., 2010, 35: 2934-2936.

[25] KeJ.LamE. Y., “Using a dual-wavelength source for depth resolution enhancement in optical scanning holography,” in OSA Topical Meeting in Digital Holography and Three-Dimensional Imaging (2011), paper DTuC30.

[26] J. Ke, T.-C. Poon, E. Y. Lam. Depth resolution enhancement in optical scanning holography with a dual-wavelength laser source. Appl. Opt., 2011, 50: H285-H296.

[27] H. Ou, T.-C. Poon, K. K. Wong, E. Y. Lam. Depth resolution enhancement in double-detection optical scanning holography. Appl. Opt., 2013, 52: 3079-3087.

[28] T.-C. Poon. On the fundamentals of optical scanning holography. Am. J. Phys., 2008, 76: 738-745.

[29] H. Di, K. Zheng, X. Zhang, E. Y. Lam, T. Kim, Y. S. Kim, T.-C. Poon, C. Zhou. Multiple-image encryption by compressive holography. Appl. Opt., 2012, 51: 1000-1009.

[30] D. J. Brady, K. Choi, D. L. Marks, R. Horisaki, S. Lim. Compressive holography. Opt. Express, 2009, 17: 13040-13049.

[31] Z.-Z. Bai, Z.-Q. Wang. Restrictive preconditioners for conjugate gradient methods for symmetric positive definite linear systems. J. Comput. Appl. Math., 2006, 187: 202-226.

[32] S. Wei, S. Lai. Fast template matching based on normalized cross correlation with adaptive multilevel winner update. IEEE Trans. Image Process., 2008, 17: 2227-2235.

Haiyan Ou, Ting-Chung Poon, Kenneth K. Y. Wong, Edmund Y. Lam. Enhanced depth resolution in optical scanning holography using a configurable pupil[J]. Photonics Research, 2014, 2(2): 02000064.

本文已被 1 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!