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

Random phase-shifting digital holography based on a self-calibrated system

Open Access Open Access

Abstract

Random phase-shifting digital holography based on a self-calibrated system is proposed. In the proposed method, the hologram and the calibration interference fringes can be recorded simultaneously in a single image based on the space-division-multiplexing technique. Three randomly phase-shifted holograms and corresponding interference fringes are recorded, and the phase-shifting amount between each two adjacent holograms is calculated by the sampling Moiré method from the calibration interference fringes. A reflective object is used to demonstrate the effectiveness of the proposed method in the numerical and experiment.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Digital holography [13] is a powerful technique for three-dimensional (3-D) measurement because of its ability to acquire both amplitude and phase information of the objects. There are essentially two types of ways to overcome an inherent issue in digital holography that the conjugate image and the zeroth-order diffracted wave are superimposed on the object image, which deteriorates the image quality. One way is called off-axis digital holography [4,5]. In this method, a large angle between the reference wave and object wave is introduced when recording a hologram. The spectra of the spatial-frequency for the zeroth-order diffracted wave, the conjugate image, and the object image will be separated in the direction of the off-axis angle. Then the object image can be easily reconstructed by extracting the spectrum in spatial-frequency of the object image. However, only a small part of bandwidth can be used for image reconstruction, so that the off-axis method’s measurable area is extremely narrow. Another way is to utilize the phase-shifting method to remove the conjugate image and the zeroth-order diffracted wave efficiently.

In the general phase-shifting digital holography [69], the phase-shifting amount must be known precisely. A piezoelectric transducer (PZT) is typically used as the phase shifter to record multiple phase-shifted holograms. However, the phase-shifting errors could inevitably be introduced in phase-shifted holograms due to the nonlinearity of the PZT and environmental disturbances. To overcome the problem, numerous phase-retrieval algorithms based on statistics that can estimate the phase-shifting amount from arbitrary phase-shifted interference patterns have been proposed [1012]. Unfortunately, these methods are formidable to use for an object wave, which is not fully random or zero mean in the diffraction field. A method that sets a reference object beside the measured object to detect the phase-shifting amount was proposed [13]. Unfortunately, this method is arduous for digital holography because of the enormous difference in the reflectivity between the reference object and the measured object. Furthermore, although a closed-loop phase control system [14] and dual camera calibrated system [15,16] were developed to calibrate phase-shifting errors, an additional electronic device such as a single photodiode or camera is required in the systems.

To eliminate this inconvenience, we proposed a calibrated phase-shifting digital holography based on the space-division multiplexing technique, which can detect phase-shifting errors without using any additional electronic devices [17]. In the previous study, four phase-shifted holograms were required for the phase-shifting calculation, and a transparent object was used to demonstrate the effectiveness of the method. Note that it is easy to control the range of the object wave of a transparent object on the image sensor to avoid disturbing the calibration wave.

In this paper, we propose a random phase-shifting digital holography based on a self-calibrated system and records a reflective object for the first challenge. It is well known that the investigation of phase-shifting digital holography has entered to two-step phase-shifting method [8]. However, the power fluctuation of the laser could strongly influence the quality of the reconstructed image in the two-step phase-shifting method. Thus, we use the three-step random phase-shifting method [18] to confirm the effectiveness of the proposed method in the numerical simulation and experiment.

2. Principle

2.1 Random phase-shifting method

Assuming that the complex distribution of an object wave at the image sensor plane is $O({x,y} )= A({x,y} )\textrm{exp}[{i\phi ({x,y} )} ]$, the complex distribution of reference wave is ${R_j}({x,y} )= {A_R}\exp [{i{\delta_j}({x,y} )} ]$, $({j = 1,2,3} )$, the intensity distribution of the recorded holograms can be described as

$${I_1}(x,y;{\delta _1}) = {A^2} + A_R^2 + A{A_R}\exp [i(\phi - {\delta _1})] + A{A_R}\exp [i({\delta _1} - \phi )].$$
$${I_2}(x,y;{\delta _2}) = {A^2} + A_R^2 + A{A_R}\exp [i(\phi - {\delta _2})] + A{A_R}\exp [i({\delta _2} - \phi )].$$
$${I_3}(x,y;{\delta _3}) = {A^2} + A_R^2 + A{A_R}\exp [i(\phi - {\delta _3})] + A{A_R}\exp [i({\delta _3} - \phi )].$$
Solving these equations for $\textrm{cos}(\phi )$ and $\textrm{sin}(\phi )$ yields, the complex distribution of the object wave at the image sensor plane can be obtained by the following equation [18]:
$$\begin{aligned} O(x,y) &= \frac{1}{{4{A_R}\sin [({\delta _3} - {\delta _2})/2]}}\left\{ {\frac{{\exp [i({\delta_1} + {\delta_2})/2]}}{{\sin [({\delta_3} - {\delta_1})/2]}}} \right.\\ &\times ({I_1} - {I_3})\left. { - \frac{{\exp [i({\delta_1} + {\delta_3})/2]}}{{\sin [({\delta_2} - {\delta_1})/2]}}({I_1} - {I_2})} \right\}. \end{aligned}$$

Then, the amplitude and phase distributions of the object wave at the object plane are obtainable through the diffraction integral calculation. To calculate the complex distribution of the object wave correctly, all the ${\delta _j}$ should be known precisely. Therefore, we proposed a calibrated phase-shifting digital holography based on the space-division multiplexing technique, which can detect the phase-shifting amount precisely.

2.2 Calibrated phase-shifting digital holography

The principle of the calibrated phase-shifting digital holography based on the space-division multiplexing technique is shown in Fig. 1. Make the reference wave and object wave cover the entire area of the image sensor, and the calibration wave impinges a small area of the image sensor. Furthermore, the intensities of the object and calibration waves are adjusted to low and high, respectively. Such an adjustment creates a recorded image including two types of information: one is the hologram that is used to reconstruct the object image, and the other is the interference fringes that are used to detect the phase-shifting amount, as indicated in the left-hand side of Fig. 1.

 figure: Fig. 1.

Fig. 1. Principle of the calibrated phase-shifting digital holography based on the space-division multiplexing technique.

Download Full Size | PDF

The phase-shifting amount only depends on the change in the optical length of the reference wave. Consequently, the phase-shifting amounts of the image-reconstruction and calibration areas will be changed synchronously. Then the phase-shifting amount between two consecutive phase-shifted holograms can be calculated from the calibration areas by certain image-processing methods, such as the sampling Moiré method [1921] and windowed Fourier transform method [22]. The accuracy of the sampling Moiré method can theoretically reach 1/1000 of the interference fringe pitch [21], and recently, we confirmed that the sampling Moiré method has the powerful ability to resist random noise [23]. Therefore, we applied the sampling Moiré method to detect the phase-shifting amount from the interference fringes. The maximum measurable displacement by the sampling Moiré method is half the pitch of the fringe patterns. Therefore, the phase-shifting amount between any two consecutive images should be small than π. First, the interference fringes in the calibration areas are extracted. Next, multiple phase-shifted Moiré fringes can be obtained by down-sampling and intensity interpolation processing. Then three phase distributions that are required for Eq. (4) can be obtained through the discrete Fourier transform algorithm [19].

3. Numerical simulation

A numerical simulation was conducted to verify the effectiveness of the proposed method. In the general phase-shifting digital holography, the phase-shifting amount is set as a constant. To achieve a convenient calculation, we set the phase-shifting amount to π/2, and the phase value of the reference wave in the first image was set as zero. It means that the theoretical phase-shifting amounts in the second and third images were π/2 and π, respectively. However, the phase-shifting amount is formidable to achieve precisely due to the precision of the phase shifter. Then the phase-shifting errors at a maximum of 20% were randomly introduced in the holograms and the interferograms. Here, the introduced true phase-shifting amount of the reference wave in the second and third images were 1.8585 and 2.9927 radians, respectively.

Besides, some random noise at a maximum of 20% was also added in the calibration area, because it is inevitable for a part of the object wave to enter the calibration area. The images with 1024×1024 pixels of 3.45µm pixel pitch shown in Figs. 2(a) and 2(b) were treated as the object. The pixel value of the amplitude image was from 10 (black area) to 255 (white area), and the phase distribution was a random pattern. The phase value was ranged from -π to π, and the histogram of the phase image is indicated in Fig. 2(c). The wavelength of the laser and the distance between the object and image sensor were assumed to be 473 nm and 100 mm, respectively. Then three phase-shifted images at the image sensor plane were generated, as described in Figs. 2(d)–2(f).

 figure: Fig. 2.

Fig. 2. Simulated object. (a) and (b) are amplitude and phase distributions, (c) is the histogram of the phase distribution, (d)-(f) are the phase-shifted images, respectively.

Download Full Size | PDF

The detected phase values of the reference wave by the sampling Moiré technique were 1.8581 and 2.9919. Here, we call the method that calculates the object image by the Eq. (4) with the assumed the phase-shifting amounts π/2 and π, as the conventional method. On the other hand, the method that reconstructs the object image by the Eq. (4) with the detected the phase-shifting amounts is called as the proposed method. Then the object images were reconstructed by the conventional method and the proposed method. The numerical simulation results are represented in Fig. 3.

 figure: Fig. 3.

Fig. 3. Simulation results. (a) and (b) are the amplitude and phase distributions obtained by the conventional method; (c) and (d) are the amplitude and phase distributions reconstructed by the proposed method; (e)-(h) are the magnified images of the areas indicated in (a)-(d), respectively.

Download Full Size | PDF

Obviously, the quality of the reconstructed images by the proposed method is much better than that obtained by the conventional method. In addition, the root-mean-square errors (RMSEs) of the reconstructed images were calculated and the results are presented in Table 1. The effectiveness of the proposed method is numerically confirmed by the simulation.

Tables Icon

Table 1. RMSE results of the amplitude and phase obtained by the conventional and the proposed methods.

4. Experiment

4.1 Optical setup

The experimental setup of the random phase-shifting digital holography measurement system is presented in Fig. 4. The collimated beam was divided into two beams by a beam splitter (BS) 1, one of which passed through BS 2 for illuminating the object. The other beam worked as the reference wave whose optical length was shifted by a PZT (including a mirror). The beam reflected from BS 2 worked as a calibration wave that passed through a slit to make it impinge a small area of the image sensor. By adjusting the intensity of the object wave, the reference wave, and the calibration wave, an image data which included the hologram and interference fringes were recorded, as described in Fig. 5. In the experiment, a Nd:YAG laser with a wavelength of 473 nm and output power of 50 mW was utilized as the light source. A complementary metal-oxide semiconductor camera (VCXU-50, Baumer, Inc.) with a resolution of 2448 × 2048 pixels and a 3.45-µm pixel pitch was used for recording three phase-shifted holograms and interference fringes. To compare with the conventional method, the phase-shifting amount was set to π/2.

 figure: Fig. 4.

Fig. 4. Optical setup of the proposed system. BS1–5: beam splitter, BE: beam expander, M1-4: Mirror, L: lens.

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. A typical example of the recorded image in the experiment.

Download Full Size | PDF

4.2 Experiment results and discussion

A Japanese 100-yen coin was set as the object at the position 280 mm away from the image sensor, and fixed on a rotation stage, as presented in Fig. 6. A PZT (PAZ005, Thorlabs, Inc.) with a resolution of 5 nm in the closed-loop mode was used as the phase shifter. The shift amount of the PZT was set to 60 nm, equivalent to 1.594 radian in theory. Then the phase-shifting of ${\delta _2},{\delta _3}$ were 1.594 and 3.188 radians. We recorded three phase-shifted images and reconstructed the object images by using the conventional method and the proposed method. The detected ${\delta _2},{\delta _3}$ were 1.982 and 3.917 radians. The reason for causing the difference between the set values and detected values mainly has two respects. On the one hand caused by the precision of the PZT, on the other hand was that the surface of the mirror mounted on the PZT was not strictly perpendicular to the incoming beam. The amplitude images of the entire object are represented in Figs. 7(a) and 7(b). Figures 7(c) and 7(d) show the magnified images of the areas indicated in Figs. 7(a) and 7(b). The right side of the reconstructed image is the conjugate image area. In theory, the conjugate image should be eliminated completely by phase-shifting calculation. However, the signal of the conjugate image was residual on the right side of the image that reconstructed by the conventional method. The profiles of the lines located in the conjugate image area in Figs. 7(c) and 7(d) were compared and the results are indicated in Fig. 7(e). Furthermore, we reconstructed the object image in the opposite direction (−280 mm), as described in Fig. 8. As can be seen, the conjugate image which symmetrically appears with the object image in the opposite direction could be observed clearly in the conventional reconstruction. On the other hand, the conjugate image was completely removed because of the correct phase-shifting calculation. The experiment provides a reasonable capability for high-quality image reconstruction by the proposed method.

 figure: Fig. 6.

Fig. 6. Photograph of the measurement object.

Download Full Size | PDF

 figure: Fig. 7.

Fig. 7. Experiment results: (a) and (b) are the amplitude images reconstructed by the conventional method and the proposed method, respectively; (c) and (d) are the magnified images in the areas indicated in (a) and (b), respectively; (e) is the profiles of amplitude values along lines indicated in (c) and (d), respectively.

Download Full Size | PDF

 figure: Fig. 8.

Fig. 8. Inverse reconstruction of the amplitude image: (a) the conventional method, (b) the proposed method.

Download Full Size | PDF

It is well known that the phase measurement is an important application of digital holography in experimental mechanics. Here, we simply manually rotated the stage to give a bitty rigid-body displacement to the object, and recorded phase-shifted holograms before and after rotations. The amplitude images reconstructed by the conventional method and the proposed method, and the contrast of the reconstructed amplitude images was adjusted under the same conditions to observe the conjugate image clearly, as shown in Figs. 9(a) and 9(b). The conjugate image can be observed clearly in the reconstructed image by the conventional method due to the phase-shifting errors. An example of the reconstructed phase image is shown in Fig. 9(c). The phase of the conjugate image is arduous to observe directly from the phase image because of the rough surface of the object and speckle noise. To reduce the speckle noise, we utilized a windowed phase-shifting digital holography [24] to calculate the phase difference between before and after rotations by conventional method and the proposed method. The extracted phase difference images obtained by each method are displayed in Figs. 9(d) and 9(f). Besides, the magnified images in the areas indicated in Figs. 9(d) and 9(f) are indicated in Figs. 9(e) and 9(g), respectively. The phase difference image obtained by the proposed method is more precise and has little phase variation than the conventional method. There results demonstrated the effectiveness of the proposed method in both the image reconstruction and the phase distribution calculation.

 figure: Fig. 9.

Fig. 9. Experiment results of phase difference measurement: (a) and (b) are the amplitude images reconstructed by the conventional method and the proposed method before the object rotated; (c) is the phase image reconstructed proposed method before the object rotated; (d) and (f) are the phase difference images obtained by the conventional method and the proposed method, respectively; (e) and (g) are the magnified images in the areas indicated in (d) and (f).

Download Full Size | PDF

5. Conclusions

A random phase-shifting digital holography based on a self-calibrated system was firstly presented to record a reflective object. In the proposed method, only three randomly phase-shifted holograms and interference fringes were required, and the phase-shifting amount could be detected precisely without utilizing any additional electronic devices such as a single photodiode or camera. By the proposed system, high-quality image reconstruction and phase measurement can be achieved conveniently even if using an inexpensive PZT that cannot perform accurate phase-shifting. Therefore, it is expected to contribute to various fields such as micro-displacement measurement, 3D object recognition, cell measurement, and 3D surface measurement, etc.

Funding

Japan Society for the Promotion of Science (19K15003, 20H02038).

Disclosures

The authors declare no conflicts of interest.

References

1. J. W. Goodman and R. W. Lawrence, “Digital image formation from electronically detected holograms,” Appl. Phys. Lett. 11(3), 77–79 (1967). [CrossRef]  

2. B. Javidi and E. Tajahuerce, “Three-dimensional object recognition by use of digital holography,” Opt. Lett. 25(9), 610–612 (2000). [CrossRef]  

3. D. Lebrun, A. M. Benkouider, S. Coëtmellec, and M. Malek, “Particle field digital holographic reconstruction in arbitrary tilted planes,” Opt. Express 11(3), 224–229 (2003). [CrossRef]  

4. S. Chen, C. Li, C. Ma, T.-C. Poon, and Y. Zhu, “Phase sensitivity of off-axis digital holography,” Opt. Lett. 43(20), 4993–4996 (2018). [CrossRef]  

5. X. He, C. V. Nguyen, M. Pratap, Y. Zheng, Y. Wang, D. R. Nisbet, R. J. Williams, M. Rug, A. G. Maier, and W. M. Lee, “Automated Fourier space region-recognition filtering for off-axis digital holographic microscopy,” Biomed. Opt. Express 7(8), 3111–3123 (2016). [CrossRef]  

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

7. I. Yamaguchi, S. Ohta, and J. Kato, “Image formation in phase-shifting digital holography and applications to microscopy,” Appl. Opt. 40(34), 6177–6186 (2001). [CrossRef]  

8. X. F. Meng, L. Z. Cai, X. F. Xu, X. L. Yang, X. X. Shen, G. Y. Dong, and Y. R. Wang, “Two-step phase-shifting interferometry and its application in image encryption,” Opt. Lett. 31(10), 1414–1416 (2006). [CrossRef]  

9. T. Nomura and B. Javidi, “Object recognition by use of polarimetric phase-shifting digital holography,” Opt. Lett. 32(15), 2146–2148 (2007). [CrossRef]  

10. L. Z. Cai, Q. Liu, Y. R. Wang, X. F. Meng, and M. Z. He, “Experimental demonstrations of the digital correction of complex wave errors caused by arbitrary phase-shift errors in phase-shifting interferometry,” Appl. Opt. 45(6), 1193–1202 (2006). [CrossRef]  

11. N. Yoshikawa, T. Shiratori, and K. Kajihara, “Robust phase-shift estimation method for statistical generalized phase-shifting digital holography,” Opt. Express 22(12), 14155–14165 (2014). [CrossRef]  

12. N. Yoshikawa, S. Namiki, and A. Uoya, “Object wave retrieval using normalized holograms in three-step generalized phase-shifting digital holography,” Appl. Opt. 58(5), A161–A168 (2019). [CrossRef]  

13. G. Mauvoisin, F. Brémand, and A. Lagarde, “Three-dimensional shape reconstruction by phase-shifting shadow moiré,” Appl. Opt. 33(11), 2163–2169 (1994). [CrossRef]  

14. C. R. Mercer and G. Beheim, “Fiber optic phase stepping system for interferometry,” Appl. Opt. 30(7), 729–734 (1991). [CrossRef]  

15. P. Xia, Q. Wang, S. Ri, and H. Tsuda, “Calibrated phase-shifting digital holography based on dual-camera system,” Opt. Lett. 42(23), 4954–4957 (2017). [CrossRef]  

16. P. Xia, Q. Wang, S. Ri, and H. Tsuda, “Nanometer-order thermal deformation measurement by a calibrated phase-shifting digital holography system,” Opt. Express 26(10), 12594–12604 (2018). [CrossRef]  

17. P. Xia, Q. Wang, S. Ri, and H. Tsuda, “Calibrated phase-shifting digital holography based on space-division multiplexing,” Opt. Laser Eng. 123, 8–13 (2019). [CrossRef]  

18. L. Z. Cai, Q. Liu, and X. L. Yang, “Phase-shift extraction and wave-front reconstruction in phase-shifting interferometry with arbitrary phase steps,” Opt. Lett. 28(19), 1808–1810 (2003). [CrossRef]  

19. S. Ri, T. Muramatsu, M. Saka, K. Nanbara, and D. Kobayashi, “Accuracy of the sampling moiré method and its application to deflection measurements of large-scale structures,” Exp. Mech. 52(4), 331–340 (2012). [CrossRef]  

20. S. Ri, M. Fujigaki, and Y. Morimoto, “Sampling moiré method for accurate small deformation distribution measurement,” Exp. Mech. 50(4), 501–508 (2010). [CrossRef]  

21. S. Ri and T. Muramatsu, “Theoretical error analysis of the sampling moiré method and phase compensation methodology for single-shot phase analysis,” Appl. Opt. 51(16), 3214–3223 (2012). [CrossRef]  

22. K. Qian, “Two-dimensional windowed Fourier transform for fringe pattern analysis: Principles, applications and implementations,” Opt. Laser Eng. 45(2), 304–317 (2007). [CrossRef]  

23. S. Ri, Q. Wang, P. Xia, and H. Tsuda, “Spatiotemporal phase-shifting method for accurate phase analysis of fringe pattern,” J. Opt. 21(9), 095702 (2019). [CrossRef]  

24. Y. Morimoto, M. Toru, M. Fujigaki, and N. Kawagishi, “Subnanometer displacement measurement by averaging of phase difference in windowed digital holographic interferometry,” Opt. Eng. 46(2), 025603 (2007). [CrossRef]  

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

Fig. 1.
Fig. 1. Principle of the calibrated phase-shifting digital holography based on the space-division multiplexing technique.
Fig. 2.
Fig. 2. Simulated object. (a) and (b) are amplitude and phase distributions, (c) is the histogram of the phase distribution, (d)-(f) are the phase-shifted images, respectively.
Fig. 3.
Fig. 3. Simulation results. (a) and (b) are the amplitude and phase distributions obtained by the conventional method; (c) and (d) are the amplitude and phase distributions reconstructed by the proposed method; (e)-(h) are the magnified images of the areas indicated in (a)-(d), respectively.
Fig. 4.
Fig. 4. Optical setup of the proposed system. BS1–5: beam splitter, BE: beam expander, M1-4: Mirror, L: lens.
Fig. 5.
Fig. 5. A typical example of the recorded image in the experiment.
Fig. 6.
Fig. 6. Photograph of the measurement object.
Fig. 7.
Fig. 7. Experiment results: (a) and (b) are the amplitude images reconstructed by the conventional method and the proposed method, respectively; (c) and (d) are the magnified images in the areas indicated in (a) and (b), respectively; (e) is the profiles of amplitude values along lines indicated in (c) and (d), respectively.
Fig. 8.
Fig. 8. Inverse reconstruction of the amplitude image: (a) the conventional method, (b) the proposed method.
Fig. 9.
Fig. 9. Experiment results of phase difference measurement: (a) and (b) are the amplitude images reconstructed by the conventional method and the proposed method before the object rotated; (c) is the phase image reconstructed proposed method before the object rotated; (d) and (f) are the phase difference images obtained by the conventional method and the proposed method, respectively; (e) and (g) are the magnified images in the areas indicated in (d) and (f).

Tables (1)

Tables Icon

Table 1. RMSE results of the amplitude and phase obtained by the conventional and the proposed methods.

Equations (4)

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

I 1 ( x , y ; δ 1 ) = A 2 + A R 2 + A A R exp [ i ( ϕ δ 1 ) ] + A A R exp [ i ( δ 1 ϕ ) ] .
I 2 ( x , y ; δ 2 ) = A 2 + A R 2 + A A R exp [ i ( ϕ δ 2 ) ] + A A R exp [ i ( δ 2 ϕ ) ] .
I 3 ( x , y ; δ 3 ) = A 2 + A R 2 + A A R exp [ i ( ϕ δ 3 ) ] + A A R exp [ i ( δ 3 ϕ ) ] .
O ( x , y ) = 1 4 A R sin [ ( δ 3 δ 2 ) / 2 ] { exp [ i ( δ 1 + δ 2 ) / 2 ] sin [ ( δ 3 δ 1 ) / 2 ] × ( I 1 I 3 ) exp [ i ( δ 1 + δ 3 ) / 2 ] sin [ ( δ 2 δ 1 ) / 2 ] ( I 1 I 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.