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

Homodyne scanning holography

Open Access Open Access

Abstract

We have developed a modified version of a scanning holography microscope in which the Fresnel Zone Plates (FZP) are created by a homodyne rather than a heterodyne interferometer. Therefore, during the scanning the projected pattern on the specimen is frozen rather than varied as previously. In each scanning period the system produces an on-axis Fresnel hologram. The twin image problem is solved by a linear combination of at least three holograms taken with three FZPs with different phase values.

©2006 Optical Society of America

1. Introduction

Fluorescence microscopy is widely used in many areas of biology. Digital holography [1] is a well established technique for three-dimensional (3D) imaging of micro structures. However, digital holograms are recorded by interfering two mutual coherent beams and therefore cannot be applied to the incoherent light which is the characteristic of the emitted light in a fluorescent microscope. Scanning holography [24], on the other hand, is the only method so far that has demonstrated the ability to produce a hologram of the fluorescent emission distributed in a 3D structure. However, scanning holography, in its previous form, has been recorded by a heterodyne interferometer in which the holographic information has been coded on a high carrier frequency. Such method suffers from several drawbacks. On one hand, trying to keep the scanning time as short as possible requires using carrier frequencies which are higher than the bandwidth limit of some, or all of the electronic devices in the system. On the other hand, working with a carrier frequency, that is lower than the system limitation, extends the scanning time far beyond the minimal time needed to capture the holographic information according to the sampling theorem. Long scanning times limit the system from recording dynamical scenes.

In this paper we combine the two holographic techniques scanning and digital holography to a single method of recording on-line Fresnel holograms of microscopic fluorescence samples after laser excitation. In this modified system, the hologram is recorded without temporal carrier frequency using homodyne interferometer. By doing this we offer an improved method of 3D imaging which can be applied to fluorescence microscopy. To the best of our knowledge, it is first time that scanning holography is demonstrated by a homodyne rather than heterodyne interferometer.

2. Principles of homodyne scanning holography

In classic heterodyne scanning holography a pattern of Fresnel Zone Plates (FZP) scans the object and at each and every scanning position the light intensity is integrated by a detector. The overall process is a 2D convolution operation between the object and the FZP patterns. Our modified scanning holography is also a method of convolution by scanning the object with a set of frozen-in-time FZP patterns. The FZP is created by interference of two mutually coherent spherical waves. As shown in Fig. 1, the interference pattern is projected on the specimen, scans it in 2D, and the reflected light from the specimen is integrated on an area detector. Due to the line-by-line scanning by the FZP along the specimen, the one dimensional signal detected is composed of the entire lines of the convolution matrix between the object function and the FZP. In the computer, the detected signal is reorganized in the shape of a 2D matrix the values of which actually represent the Fresnel hologram of the specimen. The specimen we consider is 3D, and its 3D structure is stored in the hologram by the effect that during the convolution, the number of cycles of the FZP (its Fresnel number) contributed from a distant object point is slightly smaller than the number of cycles of the FZP contributed from closer object points.

 figure: Fig. 1.

Fig. 1. Optical setup: EOPM, electro-optic phase modulator introducing a phase difference between the two beams; BS, beam splitter; PI, piezo XY stage, OBJ, objective; PM, photomultiplier tube detector; LPF, lowpass filter; PC, personal computer.

Download Full Size | PDF

As mentioned above the FZP is the intensity pattern of the interference between two spherical waves given by,

F(x,y,z)=Ap(x,y){2+exp[iπλ(γ+z)(x2+y2)+iθ]+exp[iπλ(γ+z)(x2+y2)iθ]},

where p(x,y) is a disk function with the diameter D that indicates the limiting aperture on the projected FZP, A is a constant, θ is the phase difference between the two spherical waves and λ is the wavelength of the light source. The constant γ indicates that at a plane z=0 there is effectively interference between two spherical waves, one emerging from a point at z=-γ and the other converging to a point at z=γ. This does not necessarily imply that these particular spherical waves are exclusively needed to create the FZP. For a 3D specimen S(x,y,z) the convolution with the FZP of Eq. (1) is,

O(x,y;γ)=S(x,y,z)*F(x,y,z)=2AS(x,y,z)*p(x,y)dz
+AS(x,y,z)p(xx,yy)exp{iπ[(xx)2+(yy)2]λ(γ+z)+iθ}dxdydz
+AS(x,y,z)p(xx,yy)exp{iπ[(xx)2+(yy)2]λ(γ+z)iθ}dxdydz,

where the asterisk denotes a 2D convolution. Note that O(x,y;γ) is a 2D function which is different for different values of the parameter γ. This convolution result is similar to a conventional Fresnel on-axis digital hologram, and therefore, it suffers from the same problems. Specifically, O(x,y;γ) of Eq. (2) contains three terms which represent the information on three images namely the 0th diffraction order, the virtual and the real images. Trying to reconstruct the image of the specimen directly from a hologram of the form of Eq. (2) would fail because of the disruption originated from two images out of the three. This difficulty is solved here with the same solution applied in an on-axis digital holography. Explicitly, at least three holograms of the same specimen are recoded, where for each one of them a FZP with a different phase value is introduced. A linear combination of the three holograms cancels the two undesired terms and the remaining is a complex valued on-axis Fresnel hologram which contents only the information of the single desired image, either the virtual, or the real one, according to our choice. A possible linear combination of the three holograms to extract a single convolution between the object and one of the quadratic phase function of Eq. (2), is

 figure: Fig. 2.

Fig. 2. Three recorded holograms with phase difference between the two interferometers arms of (a) 0 (b) π/2 and (c) π.

Download Full Size | PDF

OF(x,y;γ)=O1(x,y;γ)[exp(±iθ3)exp(±iθ2)]+O2(x,y;γ)[exp(±iθ1)exp(±iθ3)]
+O3(x,y;γ)[exp(±iθ2)exp(±iθ1)],

where Oi(x,y;γ) is the ith recorded hologram of the form of Eq. (2) and θi is the phase value of the ith FZP used during the recording. The choice between the signs in the exponents of Eq.v(3) determines which image, virtual or real, is kept in the final hologram. If for instance the virtual image is kept, OF(x,y;γ) is the final complex valued hologram of the form,

OF(x,y;γ)=S(x,y,z)*p(x,y)exp[iπλ(γ+z)(x2+y2)]dz.

The function OF(x,y;γ) is the final hologram which contains the information of only one image - the 3D virtual image of the specimen in this case. Such image S’(x,y,z) can be reconstructed from OF(x,y;γ) by calculating in the computer the inverse operation to Eq. (4), as follows,

S(x,y,z)=OF(x,y;γ)*exp[iπ2λz(x2+y2)].

The resolution properties of this imaging technique are determined by the properties of the FZP. More specifically the diameter D and the constant γ characterize the system resolution in a similar way to the effect of an imaging lens [5]. Suppose the image is a single infinitesimal point at z=0, then OF(x,y;γ) gets the shape of a quadratic phase function limited by a finite aperture. The reconstructed point image has a transverse diameter of 1.22λγ/D, which defines the transverse resolution, and an axial length of 8λγ2/D2 which defines the axial resolution. Note also that the width of the FZP’s last ring along its perimeter is about λγ/D, and therefore the size of the specimen’s smallest distinguishable detail is approximately equal to the width of this ring.

 figure: Fig. 3.

Fig. 3. (a) The magnitude and (b) The phase of the final hologram. [Media 1]

Download Full Size | PDF

3. Experimental results

To demonstrate the proposed technique, the setup shown in fig.1 was built on a standard wide-field fluorescence microscope (ZEISS, Axiovert S100). The specimen was a slide with several pollen grains [Carolina Biological Slide No. 30-4264(B690)] positioned at different distances from the microscope objective. The microscope objective was a Zeiss infinity corrected Fluar 20×, NA=0.75. The slide was illuminated by the FZP created by the interferometer. A DPSS-532 laser beam (λ=532nm) was split in two beams with beam expanders consisting each of a microscope objective and a 12-cm focal-length achromat as a collimating lens. One of the beams passed through an electro-optic phase modulator (New Focus 4002) driven by three (or more) constant voltage values which induce three (or more) phase difference values between the interfering beams. Note that unlike previous studies [24] there is no frequency difference between the two interfering waves since this time we record a hologram with a homodyne interferometer. The two waves were combined by the beam splitter to create an interference pattern in the space of the specimen. The pattern was then reduced in size and projected through the objective onto the specimen. The sample was scanned in a 2D raster with an X-Y piezo stage (Physik Instrument P-527). The data were collected by a GageScope CS1602 acquisition system, and data manipulation was performed by programs written in MATLAB.

 figure: Fig. 4.

Fig. 4. Movie of the image reconstruction along the light propagation axis.

Download Full Size | PDF

The three recorded holograms of the specimen taken with phase difference values of θ1,2,3=0, π/2, and π are shown in Fig. 2, respectively. In this figures it clearly appears that the dominant term is the low frequency term [the first in Eq. (2)], and therefore without mixing the three holograms in the linear combination that eliminates the low frequency along with the twin image term, there is no possibility to recover the desired image with a reasonable quality. These three holograms are substituted into Eq. (3) and yield a complex valued hologram shown in Fig. 3. This time the grating lines are clearly revealed in the phase pattern.

The computer reconstruction of two pollen grains along the z axis is shown in the movie of Fig. 4. As can be seen in this movie different parts of the pollen grains are in focus at different transverse planes.

4. Conclusions

We have presented and verified experimentally a modified version of a scanning holography system. This setup has some advantages over the previous designs but there are problems that should be solved in the future. The main technical difficulty is that during the serial scanning of each hologram there is a slight shift in the phase value of each FZP. This phase shift introduces some error and uncertainty in the measured phase value of each hologram. In future work, we intend more closely temporally correlate each point in each of the three holograms so as to minimize this phase error. Nevertheless, the proposed system is more immune from noise and it operates faster than previous scanning holography systems.

Acknowledgments

The authors thank Maria DeBernardi for valuable comments, Wenwei Zhong for technical help with the programming, and Brian Storrie for his role in envisioning the construction of a holographic microscope. This research is supported by the National Science Foundation grant DBI-0420382. J. Rosen’s research is supported in part by the Israel Science Foundation grant 119/03. G. Indebetouw’s research is supported in part by the National Institute of Health grant 5R21 RR18440.

References and links

1. I. Yamaguchi and T. Zhang, “Phase-shifting digital holography,” Opt. Lett. 22, 1268–1269 (1997). [CrossRef]   [PubMed]  

2. B. W. Schilling, T.-C. Poon, G. Indebetouw, B. Storrie, K. Shinoda, Y. Suzuki, and M. H. Wu, “Threedimensional holographic fluorescence microscopy,” Opt. Lett. 22, 1506–1508 (1997). [CrossRef]  

3. G. Indebetouw, A. El Maghnouji, and R. Foster, “Scanning holographic microscopy with transverse resolution exceeding the Rayleigh limit and extended depth of focus,” J. Opt. Soc. Am. A 22, 892–898 (2005). [CrossRef]  

4. G. Indebetouw, P. Klysubun, T. Kim, and T.-C. Poon, “Imaging properties of scanning holographic microscopy,” J. Opt. Soc. Am. A 17, 380–390 (2000). [CrossRef]  

5. J. W. Goodman, Introduction to Fourier Optics, 2nd ed. (McGraw-Hill, New York, 1996), 126–165.

Supplementary Material (1)

Media 1: AVI (4480 KB)     

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 (4)

Fig. 1.
Fig. 1. Optical setup: EOPM, electro-optic phase modulator introducing a phase difference between the two beams; BS, beam splitter; PI, piezo XY stage, OBJ, objective; PM, photomultiplier tube detector; LPF, lowpass filter; PC, personal computer.
Fig. 2.
Fig. 2. Three recorded holograms with phase difference between the two interferometers arms of (a) 0 (b) π/2 and (c) π.
Fig. 3.
Fig. 3. (a) The magnitude and (b) The phase of the final hologram. [Media 1]
Fig. 4.
Fig. 4. Movie of the image reconstruction along the light propagation axis.

Equations (8)

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

F ( x , y , z ) = Ap ( x , y ) { 2 + exp [ i π λ ( γ + z ) ( x 2 + y 2 ) + i θ ] + exp [ i π λ ( γ + z ) ( x 2 + y 2 ) i θ ] } ,
O ( x , y ; γ ) = S ( x , y , z ) * F ( x , y , z ) = 2 A S ( x , y , z ) * p ( x , y ) dz
+ A S ( x , y , z ) p ( x x , y y ) exp { i π [ ( x x ) 2 + ( y y ) 2 ] λ ( γ + z ) + i θ } dx dy dz
+ A S ( x , y , z ) p ( x x , y y ) exp { i π [ ( x x ) 2 + ( y y ) 2 ] λ ( γ + z ) i θ } dx dy dz ,
O F ( x , y ; γ ) = O 1 ( x , y ; γ ) [ exp ( ± i θ 3 ) exp ( ± i θ 2 ) ] + O 2 ( x , y ; γ ) [ exp ( ± i θ 1 ) exp ( ± i θ 3 ) ]
+ O 3 ( x , y ; γ ) [ exp ( ± i θ 2 ) exp ( ± i θ 1 ) ] ,
O F ( x , y ; γ ) = S ( x , y , z ) * p ( x , y ) exp [ i π λ ( γ + z ) ( x 2 + y 2 ) ] dz .
S ( x , y , z ) = O F ( x , y ; γ ) * exp [ i π 2 λ z ( x 2 + y 2 ) ] .
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.