• Chinese Optics Letters
  • Vol. 19, Issue 7, 070501 (2021)
Jun Long1, Ping Cai1, Chiyue Liu1, Weijuan Qu2, and Hao Yan1、*
Author Affiliations
  • 1Department of Instrument and Science in Engineering, School of Electronic Information and Electrical Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
  • 2Ngee Ann Polytechnic, Singapore 599489, Singapore
  • show less
    DOI: 10.3788/COL202119.070501 Cite this Article Set citation alerts
    Jun Long, Ping Cai, Chiyue Liu, Weijuan Qu, Hao Yan. Aperture synthesis based solely on phase images in digital holography[J]. Chinese Optics Letters, 2021, 19(7): 070501 Copy Citation Text show less

    Abstract

    Aperture synthesis is an important approach to improve the lateral resolution of digital holography (DH) techniques. The limitation of the accuracy of registration positions between sub-holograms affects the quality of the synthesized image and even causes the failure of aperture synthesis. It is a major issue in aperture synthesis of DH. Currently intensity images are utilized to find the registration positions of sub-holograms in aperture synthesis. To improve the accuracy of registration positions, we proposed a method based on similarity calculations of the phase images between sub-holograms instead of intensity images. Furthermore, a quantitative indicator, degree of image distortion, was applied to evaluate the synthetic results. Experiments are performed and the results verify that the proposed phase-image-based method is better than the state-of-the-art intensity-image-based techniques in the estimation of registration positions and provides a better synthesized final three-dimensional shape image.

    1. Introduction

    Digital holography (DH)[1] provides access to the quantitative measurement of in-depth information of an object. Due to its advantages, DH has been applied widely in various fields, such as microscopic three-dimensional (3D) measurement and imaging, including 3D cell imaging[2,3], 3D micro-structure measurement[4], super-resolution measurement[5], 3D micro-particles tracking[6,7], etc. Although DH has many advantages compared with conventional holography, its lateral resolution and field of view (FOV) are its major deficiencies due to digital recording devices.

    Nowadays, various aperture synthesis methods aiming at increasing DH’s numerical aperture have been proposed to improve DH’s lateral resolution[821]. Most of them sacrifice or merely maintain the FOV instead of enlarging it. Only a few of them in DH achieved simultaneous improvements of both lateral resolution and FOV by adopting the hologram stitching method[13,16]. Furthermore, the hologram stitching method has proved to be able to improve the axial measurement accuracy of DH as well[22].

    In the hologram stitching method, multiple holograms are recorded by shifting a DH system or a measured object[13,16]. Then, multiple sub-holograms are stitched together to form an enlarged synthesized hologram. Two key issues in hologram stitching are (1) inconsistent phase errors among the sub-holograms; (2) errors of registration positions of sub-holograms between the real positions and the registered positions. Since the instability of the mechanical movement system is inevitable, the above errors perpetuate, which affect the effectiveness of stitching and could lead to failure of hologram stitching.

    To overcome phase inconsistency, various techniques have been introduced[2328], including digital phase masks[2325], phase error model[29], and physical phases[26,28]. However, we would like to point out that, in hologram stitching, the compensation of phase errors is relevant to registration errors in the second issue. Instead, the compensation of phase errors is affected by registration errors. If the registration positions of sub-holograms are incorrect, the phase errors could not be correctly compensated. Hence, both registration errors and phase errors need to be correctly compensated.

    To eliminate the registration position errors in the second issue, registration positions among sub-holograms are re-evaluated. Various methods have been proposed to re-evaluate the registration position. Cross-correlation (CC) methods based on similarity tests of overlapping areas in the reconstructed intensity images[3033] have been used to re-evaluate the registration position. An image-based method adopting the sharpness as the metric to search for correct registration positions has been proposed in Ref. [34]. Mutual information (MI) has also been proposed as a metric to automatically re-evaluate the registration positions in Ref. [35]. Dai et al. in Ref. [36] mentioned an image profile technique to obtain the location between adjacent sub-holograms, but no experiments were presented, and the images used have not been indicated as phase images or intensity images.

    Currently, these above solutions are all based solely on intensity images to estimate registration errors. However, the intensity image is sensitive to the intensity distribution of the illumination light, which affects the correctness of the registration error estimation. Furthermore, the first issue is based on phase images, and the final 3D image of DH is derived solely from phase images. Hence, solving the two issues solely by phase images may have advantages and benefit the final 3D topography image result.

    Therefore, in this work, a hologram stitching method solely based on phase images is proposed for aperture synthesis in DH, by which both registration position errors and phase errors among sub-holograms are corrected. Both the simulation and experiment are performed. The results show that, compared with the state-of-the-art intensity-based techniques, the phase-based scheme performs better in compensation of the registration errors among sub-aperture holograms and provides a better synthesized final 3D topography image result.

    2. Experiment and Method

    2.1. Experiment setup

    In this work, a lensless off-axis DH geometry is adopted, shown in Fig. 1. A HeNe laser of 633 nm wavelength and 10 mW power is used as a light source. The light is split into an object illumination beam and a reference beam by a beam splitter. The object beam Es reflected by the object interferes with the reference beam R reflected by a plane mirror to generate a hologram I on the CCD (with a pixel size of 4.65μm×4.65μm and a pixel number of 1280×960) plane. The CCD records the holograms.

    Experimental setup of off-axis DH for aperture synthesis.

    Figure 1.Experimental setup of off-axis DH for aperture synthesis.

    After hologram recording, hologram reconstruction is performed. We apply Fourier filtering to the recorded hologram I and extract the object wavefront Es at the hologram plane. The reconstructed image of the object Eo can be calculated from Es by wave propagation from the hologram plane to the image plane by the angular spectrum method[37] with a reconstruction distance z=119.7mm in this work. Then, the intensity Io and the phase φo of the reconstructed image at the image plane can be derived. Ideally, once the sub-holograms of all sub-apertures are collected by moving the stage, which is a two-axis translation stage (XR50C/M, Thorlabs) in this work, left to right, bottom to top, we can stitch them together by placing them at the designated positions to get a stitched hologram of a full synthetic aperture. After that, we follow the same reconstruction procedures as above.

    2.2. Error sources of synthetic aperture holography

    However, in practice, the designated positions have errors, which make it difficult to generate a correct synthetic aperture hologram by directly synthesizing them. Furthermore, the temporal and spatial variations in the acquisition process also cause errors[38]. Three major errors focused in this work are: Phase error in each sub-aperture hologram: relative changes of the reference field to the object field caused by vibration and temperature fluctuations.Registration error: the difference between real positions of sub-holograms and their designated positions in in-plane directions (x and y directions) caused by the inaccuracy of the mechanical movement of the translation stage.Phase error between sub-holograms: phase error between different sub-holograms caused by the instability of mechanical movement of the translation stage in the out-of-plane direction (z-direction).

    3. Estimation of Registration Error and Evaluation of Estimated Registration Error

    3.1. Estimation of registration error

    The procedures of the proposed aperture synthesis are presented as follows. Firstly, multiple sub-aperture holograms with overlapping regions are recorded, and the optical field Eo of each sub-aperture is acquired by numerical reconstruction. Secondly, the first major error discussed above, the phase errors φap in each sub-hologram, is estimated and compensated by the method described in Refs. [1115,28] in this work. After that, the second major error, the registration error (e^up,e^vp), is estimated to get the correct registration position of adjacent sub-apertures [note that the spatial coordinates of the object plane and detector plane are defined as (x, y) and (u, v), respectively]. Then, the third major error, the phase error φepq between adjacent pth and qth sub-apertures, is estimated and compensated. We want to emphasize that the compensation accuracy of phase error φepq depends on the estimation accuracy of the registration error (e^up,e^vp). After the compensation of all three types of errors above, a large hologram is synthesized. Then, we reconstruct the synthesized optical field Eo, the intensity image Io, and phase image φo.

    The procedures of the proposed method aiming to estimate the registration errors (e^up,e^vp) include four steps: Step 1. Set up a set of all possible registration errors (eupk,evpl) for the pth sub-aperture, where k and l are integer indices, satisfying KkK and LlL. K and L are the maximum errors in the u direction and v direction, respectively, which are determined by the movement accuracy of the translation stage.Step 2. The phase images φop(ui,vj) and φoq(ui,vj) of two neighboring sub-apertures (note that in the overlapping area with I×J pixels) are used to calculate the constant phase error φIJpq between the pth and qth sub-apertures by Eq. (1): φIJpq=1I×J[i=1Ij=1JφIJp(ui,vj)φIJq(ui,vj)].The phase of the qth sub-aperture is compensated as φoq(ui,vj)=φoq(ui,vj)+φIJpq.Step 3. Assuming the qth sub-aperture moves with the designated distances du and dv with registration errors (eupk, evpl) relative to the pth sub-aperture in the u direction and v direction, respectively, we have ui=ui+du+eupk,vj=vj+dv+evpl,where ui and vj are positions of the pth sub-aperture in the u direction and v direction. The similarities between phase images φIJp(ui,vj) and φoq(ui,vj) are calculated by various methods listed in Refs. [35,39,40], such as CC criterion, normalized CC (NCC) criterion, zero-NCC (ZNCC) criterion, MI criterion, sum of squared differences (SSD) criterion, normalized SSD (NSSD) criterion, and zero-NSSD (ZNSSD) criterion. We want to emphasize that the correlation coefficient CNSSD calculated by NSSD criterion can be derived from CNSSD=2(1CNCC), and the correlation coefficient CZNSSD calculated by ZNSSD criterion can be derived from CZNSSD=2(1CZNCC). By traversing all possible registration errors (eupk,evpl), a set of the corresponding similarity indicators Ckl is obtained.Step 4. The optimal estimation of the registration error (e^up,e^vp) is found by minimum or maximum values of the set Ckl obtained in step 3. For CC criterion, NCC criterion, ZNCC criterion, and MI criterion, the indicator of the similarity Ckl is proportional to the similarity extent. For SSD criterion, NSSD criterion, and ZNSSD criterion, the indicator is inversely proportional to the similarity extent.

    The above algorithm with steps 1 to 4 is extended to the next neighboring sub-aperture one by one for aperture synthesis of a final whole aperture.

    The process of state-of-the-art methods is based on the similarity of intensity images in the overlapping regions of nearby sub-apertures. It normally includes all steps except step 3 in the proposed method. Nevertheless, there are limitations in the similarity estimation by intensity images. The most important one is the uneven intensity distribution caused by the uneven distribution of laser light, non-parallel between the coordinate system of the CCD and the measured object, and the off-axis angle between the reference and object lights. These factors affect the estimation accuracy.

    3.2. Evaluation of estimated registration error

    In holographic synthesis, if the registration position is correct, the difference between the synthesized and the pre-synthesized images should be minimal. To quantitatively evaluate the registration error estimated by different methods, an indicator, degree of image distortion (DID), is introduced herein, which is defined as DID=10log10(XX¯)2(XY)2,where X and Y are the pre-synthesized images and the synthesized images, respectively. X¯ is a mean of X. This indicator evaluates the difference between these two images. A larger DID value indicates a smaller difference between X and Y: when X=Y, DID.

    The procedure to verify the effectiveness of DID is shown in Fig. 2. A hologram measured in a real experiment is utilized to perform the simulation. Two sub-holograms with overlapping areas are truncated from it. We assign different registration positions and stitch them together. Then, we reconstruct the object field extracted from the synthesized hologram. Finally, the DID of the phase images before and after stitching is calculated. The DID of different assumed registration errors is shown in Fig. 3. The DID indicator reaches its maximum when the registration error is (0,0). It decreases with the increase of the registration error. Thus, DID proves to be an effective indicator to evaluate registration errors.

    Procedure to verify the effectiveness of the DID indicator.

    Figure 2.Procedure to verify the effectiveness of the DID indicator.

    Relationship between the DID indicator and the assumed registration error.

    Figure 3.Relationship between the DID indicator and the assumed registration error.

    4. Results

    The experiment setup is shown in Fig. 1. The United States Air Force (USAF) target is used as a specimen to examine the proposed method. Nine sub-holograms (3×3) covering different parts of the USAF target are recorded by moving the translation stage and transferred into the computer. The registration errors are estimated by the proposed method discussed in subsection 3.1 as well as the state-of-the-art technique[3032,34,35]. Sub-holograms are synthesized to form a large hologram with error compensated registration positions. Then, they are reconstructed by the angular spectrum method[37] with a reconstruction distance of z=119.7mm. Then, the intensity Io and phase φo of the reconstructed image at the image plane can be derived. The DID values of the synthesized phase images obtained by the state-of-the-art techniques and the proposed method are compared in Table 1 with different correlation criteria[3032,34,35].

    MethodCCNCCZNCCSSDMI
    State-of-the-art method46.1847.6147.6144.2046.98
    Proposed method49.1449.1449.1449.3449.25

    Table 1. The DID Values of Synthetic Phase Imagesa

    The synthesized phase images (the center part with 1280×1280 pixels) by different methods are presented in Fig. 4, with their corresponding DID values shown as well. From Figs. 4(a)4(f), the DID increases from 44.20, 46.18, 47.61, 49.14, and 49.34 to the maximum value of 49.49 (found by traversing all possible positions for all holograms). The quality of the synthesized image increases as well. Note that the phase jumps (as the upper right corner image of each image shows) become smaller, and the high-resolution bars located in the center of the images (as the upper left corner image of each image shows) are better resolved. It is consistent with the results of the simulation performed above. Therefore, the rule that the higher the DID, the better the synthesized quality, is still applicable.

    Comparison of synthesized phase images (the center part with 1280 × 1280 pixels) by different methods. (a)–(c) are the synthesized 3D phase images by the state-of-the-art method with SSD, CC, and ZNSSD criteria, respectively; (d), (e) are the synthesized 3D phase images by the proposed method with ZNSSD and SSD criteria, respectively; (f) is the synthesized 3D phase image by traversing all of the positions for the optimal registration positions with the maximum DID value.

    Figure 4.Comparison of synthesized phase images (the center part with 1280 × 1280 pixels) by different methods. (a)–(c) are the synthesized 3D phase images by the state-of-the-art method with SSD, CC, and ZNSSD criteria, respectively; (d), (e) are the synthesized 3D phase images by the proposed method with ZNSSD and SSD criteria, respectively; (f) is the synthesized 3D phase image by traversing all of the positions for the optimal registration positions with the maximum DID value.

    The synthesized images of the proposed methods shown in Figs. 4(d) and 4(e) are better than the synthesized images of the state-of-the-art methods shown in Figs. 4(a)4(c) in phase jumps and in lateral resolutions in the center of images, where high-resolution bars are located. It is hard to find any difference between Figs. 4(e) and 4(f) through the method of subjective observation. Thus, the registration positions estimated by the proposed method could be considered as the optimal registration positions (the maximum DID value). To highlight the above difference between the proposed and state-of-the-art methods, their best results are compared in Fig. 5. Figure 5(a) is the highlight of Fig. 4(c), while Fig. 5(b) is the highlight of Fig. 4(e). In the right corner of Figs. 5(a) and 5(b), the high-resolution bars are better resolved in Fig. 5(b) than in Fig. 5(a). Let us further focus on their profiles in detail in Figs. 5(a) and 5(b), where only part of the profile was selected (marked with a dotted line). As shown in Fig. 5(c), the red dotted profiles of the state-of-the-art method present more severe phase jumps than the blue solid profiles of the proposed method. Hence, through subjective observations, the synthesized images of the proposed methods are better than the synthesized images of the state-of-the-art methods in both phase jumps and resolution.

    Comparison of square dotted line areas in Figs. 4(c) and 4(e) by different methods. (a) is the best result of the state-of-the-art method (with ZNCC criterion), (b) is the best result of the proposed method (with SSD criterion), (c) comparison of profiles of (a) and (b) in their 450th rows (marked with dotted line).

    Figure 5.Comparison of square dotted line areas in Figs. 4(c) and 4(e) by different methods. (a) is the best result of the state-of-the-art method (with ZNCC criterion), (b) is the best result of the proposed method (with SSD criterion), (c) comparison of profiles of (a) and (b) in their 450th rows (marked with dotted line).

    The state-of-the-art method with SSD yields the worst result with the lowest DID value. This may be due to the inconsistency of the light intensity in the overlapping regions, which is caused by the uneven intensity distribution of the illumination light. By using the state-of-the-art method with NCC and ZNCC, the results improved as shown in Table 1. This is because, in the latter two criteria, where the intensities are normalized, they are less sensitive to the scale changes of intensity in the overlapping area. On the contrary, the phase images do not suffer such a problem, because the phase image is the reflection of the object surface shape unrelated to the distribution of the light illumination. Thus, the results obtained via the proposed method with all criteria are extremely similar and better than that obtained by the state-of-the-art method.

    5. Conclusions and Discussion

    Therefore, based on the above discussions, we conclude that both subjective observations and objective quantitative DID evaluations verify that our proposed methods based solely on phase images are better than the state-of-the-art methods. Also, we note that though the off-axis DH was applied in this work, any geometry of DH, which is capable of obtaining an accurate phase map from a hologram, is suitable for the proposed method.

    References

    [1] E. Cuche, F. Bevilacqua, C. Depeursinge. Digital holography for quantitative phase-contrast imaging. Opt. Lett., 24, 291(1999).

    [2] A. Anand, B. Javidi. Digital holographic microscopy for automated 3D cell identification. Chin. Opt. Lett., 12, 060012(2014).

    [3] B. Rappaz, P. Marquet, E. Cuche, Y. Emery, C. Depeursinge, P. J. Magistretti. Measurement of the integral refractive index and dynamic cell morphometry of living cells with digital holographic microscopy. Opt. Express, 13, 9361(2005).

    [4] F. Charrière, J. Kühn, T. Colomb, F. Montfort, E. Cuche, Y. Emery, K. Weible, P. Marquet, C. Depeursinge. Characterization of microlenses by digital holographic microscopy. Appl. Opt., 45, 829(2006).

    [5] W. Mingjun, F. Shaodong, W. Jigang. Pixel super-resolution lensless in-line holographic microscope with hologram segmentation. Chin. Opt. Lett., 17, 110901(2019).

    [6] H. Liu, J. Yu, T. Wang, Y. Yang, J. Wang, R. Zheng. Digital holography experiment of 3D detection of underwater bubble fields. Chin. Opt. Lett., 11, S20901(2013).

    [7] L. Miccio, P. Memmolo, F. Merola, S. Fusco, V. Embrione, A. Paciello, M. Ventre, P. A. Netti, P. Ferraro. Particle tracking by full-field complex wavefront subtraction in digital holography microscopy. Lab Chip, 14, 1129(2014).

    [8] L. Martínez-León, B. Javidi. Synthetic aperture single-exposure on-axis digital holography. Opt. Express, 16, 161(2008).

    [9] Y. Zhang, Y. Sheng, X. Lu, D. Hsu, C. Yu, Y. Luo, L. Zhong, B. Lee, C. She. Synthetic aperture holography by movement of object. Proc. SPIE, 5636, 581(2005).

    [10] H. Jiang, J. Zhao, J. Di. Numerical correction of splicing dislocation between sub-holograms in synthetic aperture digital holography using convolution approach. Chin. Opt. Lett., 10, 090901(2012).

    [11] J. H. Massig. Digital off-axis holography with a synthetic aperture. Opt. Lett., 27, 2179(2002).

    [12] F. Le Clerc, M. Gross, L. Collot. Synthetic-aperture experiment in the visible with on-axis digital heterodyne holography. Opt. Lett., 26, 1550(2001).

    [13] Y. Hao, A. Asundi. Studies on aperture synthesis in digital Fresnel holography. Opt. Lasers Eng., 50, 556(2012).

    [14] P. Feng, X. Wen, R. Lu. Long-working-distance synthetic aperture Fresnel off-axis digital holography. Opt. Express, 17, 5473(2009).

    [15] C. Yuan, H. Zhai, H. Liu. Angular multiplexing in pulsed digital holography for aperture synthesis. Opt. Lett., 33, 2356(2008).

    [16] J. Di, J. Zhao, H. Jiang, P. Zhang, Q. Fan, W. Sun. High resolution digital holographic microscopy with a wide field of view based on a synthetic aperture technique and use of linear CCD scanning. Appl. Opt., 47, 5654(2008).

    [17] V. Mico, Z. Zalevsky, P. García-Martínez, J. García. Synthetic aperture superresolution with multiple off-axis holograms. J. Opt. Soc. Am. A, 23, 3162(2006).

    [18] T. R. Hillman, T. Gutzler, S. A. Alexandrov, D. D. Sampson. High-resolution, wide-field object reconstruction with synthetic aperture Fourier holographic optical microscopy. Opt. Express, 17, 7873(2009).

    [19] L. Granero, V. Micó, Z. Zalevsky, J. García. Synthetic aperture superresolved microscopy in digital lensless Fourier holography by time and angular multiplexing of the object information. Appl. Opt., 49, 845(2010).

    [20] R. Binet, J. Colineau, J.-C. Lehureau. Short-range synthetic aperture imaging at 633 nm by digital holography. Appl. Opt., 41, 4775(2002).

    [21] D. Claus. High resolution digital holographic synthetic aperture applied to deformation measurement and extended depth of field method. Appl. Opt., 49, 3187(2010).

    [22] Y. Hao, A. Asundi. Impact of charge-coupled device size on axial measurement error in digital holographic system. Opt. Lett., 38, 1194(2013).

    [23] T. Colomb, E. Cuche, F. Charrière, J. Kühn, N. Aspert, F. Montfort, P. Marquet, C. Depeursinge. Automatic procedure for aberration compensation in digital holographic microscopy and applications to specimen shape compensation. Appl. Opt., 45, 851(2006).

    [24] F. Montfort, F. Charrière, T. Colomb, E. Cuche, P. Marquet, C. Depeursinge. Purely numerical compensation for microscope objective phase curvature in digital holographic microscopy: influence of digital phase mask position. J. Opt. Soc. Am. A, 23, 2944(2006).

    [25] S. Liu, Q. Lian, Y. Qing, Z. Xu. Automatic phase aberration compensation for digital holographic microscopy based on phase variation minimization. Opt. Lett., 43, 1870(2018).

    [26] J. K. T. Colomb, F. Charrière, C. Depeursinge, P. Marquet, N. Aspert. Total aberrations compensation in digital holographic microscopy with a reference conjugated hologram. Opt. Express, 14, 4300(2006).

    [27] C. J. Mann, L. Yu, C.-M. Lo, M. K. Kim. High-resolution quantitative phase-contrast microscopy by digital holography. Opt. Express, 13, 8693(2005).

    [28] W. Qu, C. O. Choo, V. R. Singh, Y. Yingjie, A. Asundi. Quasi-physical phase compensation in digital holographic microscopy. J. Opt. Soc. Am. A, 26, 2005(2009).

    [29] Y. Wen, W. Qu, H. Cheng, H. Yan, A. Asundi. Further investigation on the phase stitching and system errors in digital holography. Appl. Opt., 54, 266(2015).

    [30] L. Romero, F. Caldero. A Tutorial on Parametric Image Registration(2007).

    [31] A. E. Tippie, A. Kumar, J. R. Fienup. High-resolution synthetic-aperture digital holography with digital phase and pupil correction. Opt. Express, 19, 12027(2011).

    [32] L. G. Brown. A survey of image registration techniques. ACM Comput. Surv., 24, 325(1992).

    [33] X. Xu, G. Lu, G. Han, F. Gao, Z. Jiao, D. Li. Phase stitching and error correction in aperture synthesis for generalized phase-shifting interferometry. Appl. Opt., 52, 4864(2013).

    [34] S. Lim, K. Choi, J. Hahn, D. L. Marks, D. J. Brady. Image-based registration for synthetic aperture holography. Opt. Express, 19, 11716(2011).

    [35] A. Pelagotti, M. Paturzo, M. Locatelli, A. Geltrude, R. Meucci, A. Finizio, P. Ferraro. An automatic method for assembling a large synthetic aperture digital hologram. Opt. Express, 20, 4830(2012).

    [36] P. Schelkens, C. Dai, T. Ebrahimi, Y. Yu, G. Chen, G. Cristóbal, F. Truchetet, A. Asundi. Study of the holographic phase stitching technique. Proc. SPIE, 7000, 70001T(2008).

    [37] J. W. Goodman. Introduction to Fourier Optics(1996).

    [38] Y. Hao, C. Liu, J. Long, P. Cai, Q. Kemao, A. Asundi. Investigation of the systematic axial measurement error caused by the space variance effect in digital holography. Opt. Lasers Eng., 112, 16(2019).

    [39] B. Pan, K. Qian, H. Xie, A. Asundi. Two-dimensional digital image correlation for in-plane displacement and strain measurement: a review. Meas. Sci. Technol., 20, 062001(2009).

    [40] A. Giachetti. Matching techniques to compute image. Image Vis. Comput., 18, 247(2000).

    Jun Long, Ping Cai, Chiyue Liu, Weijuan Qu, Hao Yan. Aperture synthesis based solely on phase images in digital holography[J]. Chinese Optics Letters, 2021, 19(7): 070501
    Download Citation