Experimental Study of Twostep Phaseshifting Digital Holography based on the Calculated Intensity of a Reference Wave
 Author: Li Jun, Pan Yang yang, Li Jiao sheng, Li Rong, Zheng Tao
 Organization: Li Jun; Pan Yang yang; Li Jiao sheng; Li Rong; Zheng Tao
 Publish: Journal of the Optical Society of Korea Volume 18, Issue3, p230~235, 25 June 2014

ABSTRACT
Twostep quadrature phaseshifting digital holography based on the calculated intensity of a reference wave is proposed. In the MachZehnder interferometer (MZI) architecture, the method only records two quadraturephase holograms, without referencewave intensity or objectwave intensity measurement, to perform object recoding and reconstruction. When the referencewave intensity is calculated from the 2D correlation coefficient (CC) method that we presented, the clear reconstruction image can be obtained by some specific algorithm. Its feasibility and validity were verified by a series of experiments with 2D objects and 3D objects. The presented method will be widely used in realtime or dynamic digital holography applications.

KEYWORD
Twostep phaseshifting digital holography , Phaseshifting interferometry , Computational imaging

I. INTRODUCTION
Digital holography is a technique that permits digital capture of holograms and subsequent reconstruction of the original object on a digital computer. It is also a convenient form for data transmission and object recognition [1, 2]. The major advantages of digital holography, easy and highly sensitive recording of holograms and quick reconstruction, have been realized by rapidly developing CCD and computer technologies. Digital holography has been applied to microscopy [35], interferometry [67], image encryption [8, 9], and recognition of 3D objects [10, 11], and so on.
In recent years, twostep phaseshifting digital holography has been used in optical image reconstruction because its requires fewer hologram recording [12, 13]. Tatsuki Taharae et.al experimentally demonstrated the parallel twostep phaseshifting digital holography which requires only the intensity distribution of the reference wave and spatial two phaseshifted holograms [14]. JungPing Liu and TingChung Poon adopted a correlation method to calculate the reference wave intensity [15]. These methods are more efficient than three or fourstep phaseshifting holography, but they still need to record additionally the reference wave intensity or the object wave intensity. In this paper, we proposed an image reconstruction method with twostep quadrature phaseshifting digital holography, in which there is no need to record the reference wave intensity and the object wave intensity. In the system only two holograms are required to reconstruct the original object so it is more suitable for a realtime situation. We performed experiments to verify that our method can be used in 2D and 3D objects holography and there will not be convergence errors in CC calculation. We have applied this method well to holographic compressive imaging [16]. Section 2 introduced the principles of our method and Section 3 demonstrated the experimental results.
II. FUNDAMENTAL PRINCIPLES
2.1. Hologram Recording
A typical quadrature phaseshifting holographic setup is shown in Fig. 1., where a linearly polarized laser beam is expanded, collimated, and then divided into an object beam and a reference beam. The object beam illuminates plane P, where the original image O is put on the plane P. The phase of the reference wave is controlled by an electrooptic phase modulator. Then the two waves overlap to form interferograms on the CCD plane H.
Let us assume that the plane reference wave is simply given by the real amplitude
A_{r} , and a complex object field in the plane P isO (x _{0},y _{0}), the complex object field in the recording plane H isu_{H} (x_{H} ,y_{H} ). The distance between planes P and H isZ . Based on the traditional twostep quadrature phaseshifting holography [12, 13], when we regulate the electrooptic phase modulator and set the phases of the reference wave in the first and second exposure to 0 andπ / 2 , respectively, then two onaxis quadraturephase hologramsI_{H1} andI_{H2} in the CCD plane H are recorded sequentially and expressed aswhere
I_{0} is the zeroorder light given byOnly with the two recorded interferograms
I_{H1} andI_{H2} , we can reconstruct the original object. In the process, neither the reference wave intensity nor the object wave intensity needs to be recorded for reconstruction of the original image.2.2. Hologram Reconstruction
According to above image recording method, if the reference wave intensity
A_{r}^{2} and the object wave intensity are measured, the original image can be reconstructed correctly with the elimination of the twinimage noise and the zeroorder light. We show that such a goal can be achieved from an algorithm, which determines the reference wave intensityA_{r}^{2} directly from the hologram without the actual need to record the reference wave intensity at all.We find the reference wave intensity by searching all of the most possible values. From holographic imaging principle,
A_{r}^{2} should be within 0 and the maximum ofI_{H1} , we then construct a parameter called the 2D correlation coefficient to evaluate the correlation of the reconstruction results by single hologram and twostep phaseshifting holograms. Single hologram contains zeroorder light and twin image, but based on the twostep phaseshifting holography method, twostep phaseshifting holograms can restrain zeroorder light and twin image [12,13]. So if the reference wave intensity we determined is correct, there will be just one minimum point in the CC curve. Consequently, we can reconstruct the original image without convergence errors.Where abs (
E ) denotes the absolute value of (E );E_{T} is the image reconstructed byI_{H1} ;E denotes the reconstructed image at different intensity values of the reference waveA_{r}^{2} (the value is gotten by searching from 0 to max (I_{H1} ) without actual measurements) with two hologramsI_{H1} ,I_{H2} , which is calculated by Eqs. (9)(10) with the assumed valuesA_{rc}^{2} of the reference waveA_{r}^{2} ，⊗ denotes the 2D correlation operation.We then draw the curve of CC versus
A_{rc} and set the criterion such that the minimum point of the curve locates the actual value ofA_{r}^{2} (see Fig. 2(c)). With this found value ofA_{r}^{2} we can reconstruct the original image based on Eqs. (9)(11).Now in this way only two quadraturephase holograms are needed to reconstruct the original image without zeroorder image and twinimage, and the clear reconstructed image can be obtained at high speed by a certain algorithm.
The algorithm can be expressed as follows:
Replace the
A_{r} in Eqs. (1)(3) with the valueA_{rc} found in the above method, so we can construct a complex hologramI_{C} according toWhere
u_{H} is the complex amplitude of the object light on the CCD camera, which is given byu_{H} = Re(u_{H} )+i Im(u_{H} ). By taking the square of the absolute value of both sides of Eq. (5), we obtain the solution to a quadratic equation inI_{0} .Moreover, from Eqs. (1)(2), we calculate
from which, we can write
Comparing Eq. (6) with Eq. (8), we find that if
A_{rc} + Re(u_{H} ) + Im(u_{H} ) > 0, the minus sign should be chosen for Eq. (6). The quantity ofA_{rc} + Re(u_{H} ) + Im(u_{H} ) will be positive everywhere if the intensity of reference wave is larger than that of the object wave intensity in the CCD plane, which is generally true in practice [12,13]. This will ensure that Eq. (6) takes the minus sign and then we can calculateThen from Eq. (5) we can calculate the complex amplitude in the CCD plane
According to Eqs. (9)(10), we can obtain zeroorderand twinimagefree hologram
u_{H} (x_{H} ,y_{H} ) with the two captured quadraturephase interferograms and the reference wave intensity acquired with above the CC method. Onceu_{H} (x_{H} ,y_{H} ) is known, we can retrieve the complex field in the original object plane P asWhere
IFR_{z} denotes the inverse Fresnel transformation of distanceZ . That is to say, original images have been well reconstructed.Previously, Liu
et al . [15] proposed a seemingly similar correlation method to calculate the reference wave intensity. However, first they defined a different amplitude ratio A from us aswhere
R is the amplitude of the reference light on the CCD plane andϑ denotes the complex amplitude of the object light on the CCD plane, such that when A≥1, only two quadrature phase holograms are needed to obtain a twinimage and zeroorderfree object wave in their method. However, according to the equation, the value of A is related to object wave intensity on the CCD plane and diffraction distance, so in order to reconstruct different original images or dynamic scenes, different object wave intensity on the CCD plane is required to be recorded to judge whether A≥1. In contrast, in our method, we just need to evaluate the amplitude ratio of reference wave intensity and illumination light intensity in the light path. Based on CC calculation, when the amplitude ratio of reference wave intensity and illumination light intensity in the light path is larger than 1.329, the correlation coefficient of the reconstruction object and the original object is 1, thus we can reconstruct original images accurately. So our method can reach the goal to record only two quadrature phase holograms. Secondly, the correlation coefficient defined in their paper, is obviously different from ours.III. EXPERIMENTAL RESULTS
A series of experiments has been made to verify the feasibility of our proposed method. We shall now show a series of results based on the following conditions. The wavelength of the HeNe laser is 632.8 nm. The quantitative bit depth of the MVC1000 type CCD is 8 bit and the largest number of pixels is 1280 (H) by 1024 (V) with pixel size of 5.2
μ m×5.2μ m. The quantitative bit depth of the Pike F421 type CCD is 14 bit and the largest number of pixels is 2048 (H) by 2048 (V) with pixel size of 7.4μ m × 7.4μ m. The NEW FOCUS 4002M type electrooptic phase modulator is controlled by a voltage amplifier and a lag wave plate is used for phase modulation.First we verify the validity of this method by 2D objects reconstruction experiments. The object we use is the transmission type object “SCNU”, as shown in Fig. 2(a) with the size of 1.25 cm × 1.25 cm. The parameter that we use is
d = 46.1 cm. We regulate the electrooptic phase modulator and set the phases of the reference wave in the first and second exposure to 0 andπ / 2 , respectively. We get two quadraturephase interferograms on the CCD plane, one interferogramI_{H1} is shown in Fig. 2(b). Based on the idea of CC discussed previously, Fig. 2(c) shows the CC curve for differentA_{rc} values. From the figure, if the amplitude of the reference light is smaller than 0.34, there will exist convergence errors. However, this situation will not appear in the experiments. With the increase of the reference wave amplitude, there will be just one minimum point in the CC curve. And the true amplitude of the reference wave is 5.55. Fig. 2(e) shows the reconstruction by using the value found in the CC method corresponding to those of Fig. 2(c). Reconstruction results of single hologram and reconstruction results by subtracting the mean value and mean filtering are shown in Fig. 2(d), (f) and (g). Compared with Fig. 2(d), Fig. 2(e) using our method shows a zeroorderand twinimagefree reconstruction image and we can get a clear reconstruction result by subtracting the mean value (see Fig. 2(f)) and mean filtering (see Fig. 2(g)).Then, the 3D objects reconstruction has been tested. The object we use in this experiment is a plastic ring is as shown in Fig. 3(a) with size of 6 mm × 6 mm × 6 mm. The parameter is
d = 34.7 cm. We regulate the electrooptic phase modulator for 3D objects. According to 2D method, the true value of the reference wave is 8.37 (see Fig. 3(c)). Finally, we obtain the reconstruction image, as shown in Fig. 3(e). Reconstruction results of single hologram and reconstruction result by subtracting the mean filtering are shown in Fig. 3(d) and (g). From the pictures we can see that we can also reconstruct the 3D original object with zeroorderand twinimagefree.To further investigate the effectiveness of our method, we reconstruct the object by fourstep phaseshifting holography, as shown in Fig. 3(f). Due to the limitation of our experimental instrument, the experimental result for fourstep phaseshifting digital holography is not ideal. However, with this limitation, the reconstruction result by the twostep phaseshifting holography method is better than the reconstruction result by fourstep phaseshifting holography. Therefore, in the lowlevel optical equipment conditions, twostep phaseshifting holography can work better than fourstep phaseshifting holography. And our method will also get ideal results in the high precision equipment conditions as well as in fourstep phaseshifting holography. So the method proposed in this paper is effective not only in the high precision equipment conditions but also in the lowlevel optical equipment conditions. The results of the experiments show that our method not only can be applied to 2D object reconstruction, but also can be applied to 3D object reconstruction.
IV. CONCLUSIONS
In this paper the twostep phaseshifting digital holography based on the calculated intensity of reference wave algorithm is proposed. With only two captured quadraturephase holograms known, in addition to the acquired referencewave intensity, which is not directly recorded in advance but calculated from CC method with two holograms, the clear retrieved image can be reconstructed at high speed. We can obtain clear reconstruction images without referencewave intensity or objectwave intensity measurement. In the traditional twostep phaseshifting digital holography method, object wave intensity and the reference wave intensity need to be recorded to reconstruct the original image. Due to the high storage requirement of image data and long recording time, these methods have some difficulty particularly in realtime dynamic measurement. Compared with traditional twostep phaseshifting digital holography, our method just need to calculate the reference wave intensity by the CC method instead of recording object wave intensity and reference wave intensity. In this method, our system need to record only two quadrature phase holograms and we don't need to record other data on the CCD plane. It greatly reduces the recording time. For realtime dynamic measurement, we just need to modulate the phase rapidly by phase shifter and get synchronization measurement on the CCD plane to achieve realtime dynamic imaging. In addition, our method can also be used in onestep parallel phaseshifting digital holography to simplify the system. So, this proposed method can improve the speed of holographic imaging and can be considered as a realtime solution. For 3D objects, twostep phaseshifting holography works better than fourstep phaseshifting holography. It reduces by half the storage requirements of the digital hologram. The results of the experiments show that the algorithm we proposed is effective. We believe that the method proposed in this paper will be widely used in parallel holography [17,18], and will provide an effective solution to the realtime and dynamic holographic field.

[FIG. 1.] Setup of image recording with twosteponly quadrature phaseshifting digital holography.

[]

[]

[]

[]

[FIG. 2.] The experimental results of twostep algorithm based on the calculated intensity of reference beam for 2D objects: (a) original image, (b) one of the two holograms, (c) the curve of correlation coefficient versus reference wave amplitude Arc, (d) reconstruction result of single hologram, (e) reconstruction result using our method, (f) minus the mean value of (e), (g) mean filtering of (f).

[]

[]

[]

[]

[]

[]

[]

[FIG. 3.] The experimental results of twostep algorithm based on the calculated intensity of reference beam for 3D objects: (a) small plastic ring using as 3D object, (b) one of the two holograms, (c) the curve of correlation coefficient versus reference wave amplitude Arc, (d) reconstruction result of single hologram, (e) reconstruction result using our method, (f) reconstruction result using fourstep phaseshifting holography, (g) mean filtering of (e).