A new two step onaxis digital holography (DH) is proposed without any assumptions, phase shifting,or complicated optical components. A dualchannel MachZehnder interferometer was employed. Using that setup, the object field can be reconstructed requiring only two step measurements. To eliminate position difference between two chargecoupled device (CCD) cameras, a matched filter algorithm was used.Experimental results are compared to those of the traditional phase shifting technique. The proposed approach can also be applied to singleexposure onaxis DH for real time measurement.
I. INTRODUCTION
DH is a novel technique for recording and reconstructing the fully complex field of a wavefront
[1

2]
. Two processes are involved in DH: the recording process, and the hologram reconstruction process. Due to the nature of holography,three diffraction components: the dc term, object field and its conjugate, are produced in the reconstruction process.In offaxis DH, these three diffraction components are spatially separated. In contrast, in onaxis DH, the diffraction components are overlapped. However, onaxis DH has a larger field of view and space bandwidth product (SBP)than offaxis DH, features that could be very useful in some applications
[3]
. To solve the overlap problem of onaxis DH, various efforts have been made.
The use of multiple exposures with phase shifting is a simple solution and could eliminate both the dc term and conjugate field
[4

5]
. However, it requires a phase shifting mechanism during measurements. Maleki and Devaney suggested the recording of two holograms at two different distances
[6]
. By processing the holograms in the spatial frequency domain, they extracted the object image without the conjugate image. However, because singularities occur during the computation, they had to use an inverse filter and regularization, which requires more extensive computation.Kim and Javidi
[7]
also used two holograms. To reduce the dc term, the average intensity of the hologram was subtracted from the hologram. However, this method is limited to an object which has a uniform intensity distribution.Situ
[8]
reported a method which records four holograms from two recording planes. This method still requires a phase delay of reference beams during measurements,and it has a singularity problem.
Zhang and Tiziani
[9

10]
assumed a weak object wave so that the hologram signal could be approximated as an exponential function. Although this method involves simple computation and does not suffer from singularities, the assumption of the weak object wave is not always valid and favorable in terms of the signal to noise ratio. Liu and Poon
[11]
also assumed a weak object wave, although they retrieved the object field by two step measurements.
There are other purely iterative methods that use the GerchbergSaxton algorithm
[12

13]
. These methods, however,may strongly depend on initial guesswork and requires computations to obtain an optimized result.
Another approach using a tilted reference wave was proposed
[14

16]
. It has a very simple optical setup and requires no phase shifting devices. However, the isotropic phase shifting is not possible in that method. As a phaseshifting free method like the tilted reference approach, a singleexposure technique using fractional Talbot effect was reported
[17]
. However, the accurate fabrication and alignment of the gratings are difficult, and have a great effect on the resultant image.
Recently, pixelated optical components have been used
[18

22]
. These approaches can directly record phaseshifted holograms with specially designed components. However,they are difficult to fabricate and align. As a related work,Nomura
[23]
suggested a method using a randomphase reference wave. This method requires only ordinary imaging devices. However, they assumedthe same object and reference intensity on any neighboring four pixels.
In this paper, a new two step onaxis DH is proposed.Unlike the above mentioned methods, it does not require any phase shifting during the recording process, weak object wave assumption, and any other complicated optical components. The proposed method also provides a simple reconstruction algorithm and system setup.
A dualchannel MachZehnder interferometer was employed to verify our method. The matched filter algorithm
[24]
was also employed to align two CCDs. Simulation and experiment results are presented as a demonstration.
II. BASIC PRINCIPLE
FIG. 1.
is the schematic diagram of the experimental setup. The light source is a singlemode laser with 635 nm wavelength. Laser light is collimated by a spatial filter and lens, and then passes through a linear polarizer (P3) whose polarization direction is 45
^{0}
with respect to the xaxis.This beam is then divided by a beam splitter (BS1) into two paths, one of which illuminates the transparent object,and the other, which passes through a quarterwave plate,acts as the reference beam. The quarterwave plate has the fastaxis along the yaxis. Therefore,
π / 2
phase shift occurs only on the xaxis component (0
^{0}
) of reference beam.Using two polarizers (P1, P2) whose polarization directions are 90
^{0}
, 0
^{0}
respectively, we can make a
π / 2
phaseshifted reference wave on CCD2 with respect to that of the CCD1.That makes it possible to obtain two holograms simultaneously with
π / 2
phase shift between them.
Generally, in onaxis DH, the hologram can be represented as follows,
Dualchannel MachZehnder interferometer with aquarterwave plate. SF spatial filter; P polarizer; BS beamsplitter; M mirror; WP wave plate.
where, x, y are the rectangular coordinates of the recording plane,
θ
is the phase shift, O (x, y), R (x, y) are the object and reference field, and * denotes the complex conjugate. To retrieve the object field, a two hologram method
[7]
was chosen among many reconstruction approaches.The method consists of two procedures. First, a complex field is obtained at the recording plane using Eq. (2),
where, H (x, y, 0), H (x, y,
π /
2) are the hologram and its
π
/ 2 phaseshifted one. Second, object field is reconstructed by the Fresnel transformation,
where,
x ' , y '
are the rectangular coordinates of the reconstruction plane, k is the wavenumber, d is the distance between recording and reconstruction planes,
λ
is the wavelength of the light, and
F
denotes the two dimensional fast Fourier transformation.
Note that the object field can be obtained by two step measurements, if we use the MachZehnder interferometer as shown in the
FIG. 1.
First, two holograms can be obtained simultaneously, having a
π /
2 phase shift between them. Second, object intensity can be measured by blocking the reference arm. Substituting these two holograms and the object intensity into Eq. (2), and then using Eq.(3), the object field can be reconstructed. Here, note that reference intensity is constant in time. Therefore, there is no necessity for measuring reference intensity each time.
III. ANALYSIS AND RESULT
However, if the two CCDs have different positions relative to the object, the position effect should also be considered in the two hologram approach. For example, a CCD can have position information relative to the object,which are x, y, z, roll, pitch, and yaw, as shown in the
FIG. 2.
In general, the position information of two CCDs are different from each other. Therefore, a resolution target was employed as an object as shown in the
FIG. 3
(a) to examine the effect of position difference between two CCDs on the reconstructed image.
If the two CCDs have the same position relative to the object, the two holograms obtained using
FIG. 1
have a difference of
π / 2
phase shift between them. Therefore, the object field is reconstructed without any noise components as shown in the
FIG. 3
(b). However, if there is any difference between position information of the two CCDs,we cannot retrieve a clean object image as shown in the
FIG. 3
(c) and (d).
FIG. 3
(c) was obtained when there was a lateral shift of two pixels between two CCDs in both x and y directions. And,
FIG. 3
(d) was obtained when an axial shift of 500 ㎛ existed between two CCDs in the z axis. As you can see in the
FIG. 3
, a small position difference can badly influence the reconstructed image.Therefore, it is in practice important that two CCDs be aligned each other to have the same positions relative to the object, especially before any measurements.
Six geometric variables: x y z roll pitch and yaw.
In the previous paper
[24]
, we proposed an alignment method between two CCDs using a matched filter algorithm.To put it simply, the algorithm is a kind of correlation method, which compares intensities of two CCDs. When two CCDs have the same position for an object, their intensity patterns may be very similar. Therefore,the correlation between the intensities is almost one.On the other hand, when the two CCDs have different positions, the correlation between the intensities may be a value lower than one. Therefore, in the process of making the correlation approach one, the two CCDs should have the same position for an object. For detailed descriptions about matched filter algorithm, refer to the previous paper.
Simulation result: (a) original object (b) reconstructedobject without position error (c) reconstructed object withlateral position error (d) reconstructed object with axialposition error.
Experimental result: (a) proposed method (b) conventionalphaseshifting method (c) onaxis DH without phase shifting.
Threechannel MachZehnder interferometer forsingleexposure DH. NF neutral density filter.
FIG. 4
(a) is the reconstructed object intensity obtained by the proposed dualchannel onaxis DH. A part of USAF resolution target (38257, Edmund) was employed as an object as shown in the
FIG. 4
.
FIG. 4
(b) is reconstructed by the conventional phase shifting method which has only one CCD. Therefore, there are not any noises caused by the position difference between detectors.
FIG. 4
(c) is obtained by onechannel onaxis DH without any phase shifting process. As you can see, the conjugate component of the object degrades the reconstructed image.
IV. DISCUSSION
As shown in the
FIG. 4
, the proposed dualchannel onaxis DH method provides almost the same quality result compared to the conventional phase shifting scheme.Therefore, we can retrieve the object field much faster than by conventional approaches, if the position difference between two CCDs is eliminated. The proposed method consists of two processes: measurement of two phaseshifted holograms and the object intensity. In addition,there are not any limitations in our method.
Furthermore, singleexposure onaxis DH can be realized in practice based on our scheme. To realize singleexposure onaxis DH, we should use one more channel as shown in the
FIG. 5
.
FIG. 5
is almost identical to
FIG. 1
, except the added detector, CCD3. CCD3 measures object intensity,and CCD1 and 2 measure two phaseshifted holograms.Then, using Eq. (2), (3), singleexposure onaxis DH can be realized in practice without any limitations. The neutral density filter decreases the reference intensity approximately in half to match the intensity level of reference and object.
V. CONCLUSION
In this paper, a new two step onaxis DH was proposed.A dualchannel MachZehnder interferometer was employed to implement the scheme. Object field was retrieved in two step measurements. We eliminated position differences between two CCDs using matched filter algorithm, and then provided a reconstruction image which was of almost the same quality as the result from the conventional phaseshifting method. In addition, there are not any limitations in our method. Furthermore, if we use one more channel, singleexposure onaxis DH can be realized in a practical way. We believe that our scheme can be widely applied to realtime applications like monitoring dynamics of biological samples or particles in flows.
Acknowledgements
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MEST) (20100016064).
Schnars U
,
Juptner W
1994
Direct recording of holograms by a CCD target and numerical reconstruction
Appl. Opt.
33
179 
181
Schnars U
,
Juptner W
2002
Digital recording and numerical reconstruction of holograms
Meas. Sci. Technol.
13
R85 
Stern A
,
Javidi B
2007
Theoretical analysis of threedimensional imaging and recognition of microorganisms with a singleexposure online holographic microscope
J. Opt. Soc. Am.A
24
163 
168
Takaki Y
,
Kawai H
,
Ohzu H
1999
Hybrid holographic microscopy free of conjugate and zeroorder images
Appl.Opt.
38
4990 
4996
Zhang I
1997
Phaseshifting digital holography
Opt. Lett.
22
1268 
1270
Maleki M
,
Devaney A
1994
Noniterative reconstruction of complexvalued objects from two intensity measurements
Opt. Eng.
33
3243 
3253
Kim D
,
Javidi B
2006
3D object recognition using single exposure online digital holography
Proc. SPIE
6027
602703 
602709
Situ G
,
Ryle J
,
Gopinathan U
,
Sheridan J
2008
Generalized inline digital holographic technique based on intensity measurements at two different planes
Appl. Opt.
47
711 
717
Zhang Y
,
Zhang X
2003
Reconstruction of a complex object from two inline holograms
Opt. Express
11
572 
578
Zhang Y
,
Pedrini G
,
Osten W
,
Tiziani H
2004
Reconstruction of inline digital holograms from two intensity measurements
Opt. Lett.
29
1787 
1789
Liu J
,
Poon T
2009
Twosteponly quadrature phaseshifting digital holography
Opt. Lett.
34
250 
252
Zhang Y
,
Pedrini G
,
Osten W
,
Tiziani H
2003
Whole optical wave field reconstruction from double or multi inline holograms by phase retrieval algorithm
Opt. Express
11
3234 
3241
Latychevskaia T
,
Fink H
2007
Solution to the twin image problem in holography
Phys. Rev. Lett.
98
233901 
Nakadate S
,
Isshiki M
1997
Realtime vibration measurement by a spatial phaseshifting technique with a tilted holographic interferogram
Appl. Opt.
36
281 
284
Yasuno Y
,
Makita S
,
Endo T
,
Aoki G
,
Sumimura H
,
Itoh M
,
Yatagai T
2004
Oneshotphaseshifting Fourier domain optical coherence tomography by reference wavefront tilting
Opt. Express
12
6184 
6191
Toge H
,
Fujiwara H
,
Sato K
2008
Oneshot digital holography for recording color 3D images
Proc. SPIE
6912
69120U 
MartinezLeon L
,
AraizaE M
,
Javidi B
,
Andres P
,
Climent V
,
Lancis J
,
Tajahuerce E
2009
Singleshot digital holography by use of the fractional Talbot effect
Opt.Express
17
12900 
12909
Awatsuji Y
,
Sasada M
,
Kubota T
2004
Parallel quasiphaseshifting digital holography
Appl. Phys. Lett.
85
1069 
Awatsuji Y
,
Fujii A
,
Kubota T
,
Matoba O
2006
Parallel threestep phaseshifting digital holography
Appl. Opt.
45
2995 
3002
Awatsuji Y
,
Koyama T
,
Tahara T
,
Ito K
,
Shimozato Y
,
Kaneko A
,
Nishio K
,
Ura S
,
Kubota T
,
Matoba O
2009
Parallel opticalpathlengthshifting digital holography
Appl. Opt.
48
160 
167
Nomura T
,
Murata S
,
Nitanai E
,
Numata T
2006
Phaseshifting digital holography with a phase difference between orthogonal polarizations
Appl. Opt.
45
4873 
4877
Suzuki H
,
Nomura T
,
Nitanai E
,
Numata T
2010
Dynamic recording of a digital hologram with single exposure by a wavesplitting phaseshifting method
Opt. Rev.
17
176 
180
Nomura T
,
Imbe M
2010
Singleexposure phaseshifting digital holography using a randomphase reference wave
Opt. Lett.
35
2281 
2283
Lee H
,
Kim D
,
Kim S
2009
A simple and quantitative alignment procedure between solid state cameras
Opt.Express
17
23947 
23952