Annular multi-focal-phase mask multiplexing based large depth of field imaging by interferenceless coded aperture correlation holography

Extending depth-of-field (DOF) of the imaging system without modifying the structure and sacrificing imaging performances of the optical system is of great significance to broaden the capability and application of the imaging system. In this paper, the interferenceless coded aperture correlation holography(I-COACH) is developed to be a large-depth incoherent imaging system by employing an annular multi-focal coded phase mask (AM-CPM). Based on the analyses of axial defocus characteristics in I-COACH, the defocus compensation function is defined, the AM-CPM is designed and multiplexed on the system optical pupil, which plays the role of a gradual lens. In AM-CPM, multi-annular zones with different focal lengths are used to compensate different axial defocus aberrations and adjacent annular zones have symmetric axial defocus aberration correction capability according to the imaging characteristics of the system. The simulations and experimental results fully demonstrate that the axial point spread function distribution of the system obtained by AM-CPM is continuous and the development method enables the extension of the DOF of the I-COACH system by only single exposure point spread hologram. This solution is expected to provide great potential in the field of microscopic imaging and other fields of that based on I-COACH system.

Optical imaging technologies provide the most intuitive and effective way to observe the morphology, structure and even functional information of the samples. Large depth imaging can offer the whole information of expected volume once, thus both large enough imaging depth of field and high spatial resolution are pursued for imaging the structure of a whole samples at the same time or 3D single particle tracking in thick samples. Various methods of depth-of-field (DOF) extension have been developed in imaging systems, such as wave-front coding technology 1-3 , image fusion methods 4 , multi-plane imaging based on digital holography 5 and scattering based methods 6 , and DOF engineering 7 . Compared to other DOF extension imaging methods with complex matching algorithms and expensive diffractive optics, wave-front encoding technique achieves large DOF imaging by placing a phase mask on the optical pupil of the conventional imaging system to encode and modulate light waves, making the point spread function and optical transfer function of the system are not sensitive to the defocus aberration over the DOF range 8 . Therefore, wave-front encoding enables DOF extension without changing any structure of the imaging system, thus enhancing the system imaging performance.
As one of the incoherent 3D imaging techniques, interferenceless coded aperture correlation holography (I-COACH) technique attracted the focus of increasing research interest in recent years for its advantages such as wide spectrum, non-scanning 3D imaging and higher axial resolution compare to Fresnel incoherent correlation holography [9][10][11] , and various imaging enhancement methods [12][13][14][15][16][17][18][19] and potential application sceneries have been development [20][21][22][23][24][25] . This technique modulates the incident light waves by a wave-front pseudo-random coded phase mask (CPM) to record a series of point spread holograms (PSHs) along z-axial depth and one object hologram (OH) and enables thick 3D objects to be imaged within a calibrated axial range. Based on the imaging scheme of the I-COACH, it is also advantageous to combine wave-front encoding technique for modifying I-COACH as one of incoherent large depth imaging technique. The axial imaging properties of I-COACH system were

Results
The mutual constraint analysis of parameters m and α. According to Eqs. (5a) and (5b), the AM-CPM adjacent different ring zones have axially symmetric defocused aberration quadratic phase functions, which play the role of a gradual lens to extend the imaging DOF. Therefore, by modulating the parameters m and α, it is possible to obtain multifocal encoded phase masks with different DOF extensions.
In order to enable the system response of the object point in the CCD plane at any axial position to satisfy the conjugate relationship, it is necessary to calculate the correction of the system axial position d 0 with respect to the parameters m and α. As shown in Fig. 1, the front focal plane of the input lens L 0 is defined as the coordinate origin, with the positive direction to the left and the negative direction to the opposite. According to Eq. (4), the d 0 and α curve are calculated for the axial range − 0.1 m ≤ d 0 ≤ 0.1 m when m is taken as 2,4,7,8,10,13 and 15, respectively, and the results are shown in Fig. 1. According to the Fig. 1, d0 is a non-linear relationship with parameter α when m is a certain determined value. Therefore, the imaging of the two sides of the input lens L 0 in the I-COACH system is a non-symmetrical relationship, and in the same axial range, the axial defocus aberration of the lens L 0 in the negative direction is greater than that in the positive direction, which leads to different axial resolution of the two sides of the lens L 0 in the system. Meanwhile, when m is fixed at a certain value, different axial positions correspond to different values of α. However, the spatial light modulator (SLM) is loaded with a phase mask that can only satisfy the α value at a certain axial position. As a result, the object points in other axial positions are not imaged due to the non-linear relationship, which makes the object points quickly defocus when imaging in the CCD plane. Furthermore, as shown in the blue rectangular wireframe in Fig. 1 that α values approximately satisfy a certain linear relationship for a larger choice of m in a smaller axial range. The larger the m value, the smaller the fluctuation range of α values. Therefore, in the I-COACH system, the larger the axial modulation range corresponding to a certain α value is chosen, the larger the DOF ranges of the corresponding system. However, with the increasing of m, the effective area of each defocus compensation ring on the AM-CPM decreases. As the results, the equivalent numerical aperture of the system decreases, which causes the system to lose a certain degree of transverse resolution. At the same time, another important parameter (one of the very key points in the I-COACH-like systems), the signal to noise (SNR) ratio of the images on each foci within the extended DOF decreases, as well. With the increasing of m, the SNR of the in-focused images on each axial plane within the extended DOF gradually decreases. In this paper, the SNRs of the reconstructed images with m > 10 were defined to be unacceptable in the following experiments. Therefore, m = 10 is chosen as the upper limit of the extension of the DOF in the proposed method. Taking Fig. 3a, the full width at half maxima (FWHM) of the axial direction PSF at 0 mm shows that the PSF is distributed asymmetrically, with lower axial resolution in the positive direction than in the negative direction, which is consistent with the above theoretical analysis that I-COACH has similar imaging properties to regular lens imaging 10 . The axial imaging resolution of the proposed method, for m = 10 for example, can be quantified using the FWHM of the corresponding curve in the Fig. 3a. The PSF normalized intensities of point objects at different axial positions are shown in Fig. 3c-g, and it can be clearly seen that a larger m corresponding to an AM-CPM with a better DOF extension capability. Further, we verify the DOF extension capability of AM-CPM in the I-COACH system by simulations. The simulation parameters are the same as above, with the SLM loaded synthetic coded phase mask (SCPM) synthesized from CPM, AM-CPM and QPM with a focal length f slm = 150 mm, where the annular sparse CPM is obtained by the modified GS algorithms and the scattering degree σ = 0.167 [16][17][18] . The PSHs of single exposures at different axial positions were obtained, and the non-linear reconstruction (NLR) 15 was used to obtain the I-COACH axial PSF distribution corresponding to different m in the range of − 10 mm ≤ d 0 ≤ 10 mm for point objects. The distribution of the axial PSF of the point object corresponding to different m is shown in Fig. 4. The PSF obtained by SCPM has a continuous distribution over a certain range.
The experiment of DOF extension in I-COACH system. As shown in Fig. 5, a two-channel experimental system was setup to verify the DOF extension capability of the proposed AM-CPM for the I-COACH system. In experiments, a pinhole or object was illuminated by the incoherent light emitting diode (LED) (Thorlabs LED 625L4, 700mW, center wavelength of λ = 625 nm, Δλ = 17 nm), and the light from beam splitter BS 1 was collimated by the lens L 0 and passed through a polarizer P. The polarizer P polarizes the light along the orientation of the active axis of the spatial light modulator (SLM, Holoeye PLUTO, 1080 × 1920 pixels, 8 μm pixel pitch, phase-only modulation) located at a distance of 55 mm from the L 0 . On the SLM, a SCPM phase mask is displayed whereas its phase is the combination of CPM with the scattering degree σ = 0.167, AM-CPM and QPM. The light modulated by the SLM was collected by a Charge Coupled Device (CCD, Thorlabs CS235MU, 1200 × 1920pixels, 5.86 μm pixel pitch, and monochrome) located at a distance of d s = 171 mm from the SLM.
Corresponding to the simulation results in Fig. 4, in the first experiments, only channel 1 worked and the AM-CPM corresponding to the different m was used to record the PSHs. Taking m = 2 as an example, the SCPM corresponding to m = 2 was loaded into the SLM and the PSHs of 20 μm pinhole at different axial positions were recorded in steps of 1 mm within − 6 mm ≤ d 0 ≤ 0 mm. The PSH recorded at the axial position of 0 mm was separately reconstructed with the PSH at other axial positions using an NLR method with a modulation parameter of p = − 0.3 and o = 0.8 to obtain pinhole reconstructed images to demonstrate the DOF extension capability of the AM-CPM for the I-COACH system. The corresponding PSHs were obtained in the similar way for different m and the normalized reconstructed image of the pinhole was obtained from the NLR, the results of which are shown in Fig. 6. Comparison with the lens imaging in Fig. 6a shows that different m has different DOF extension range and can obtain a better pinhole reconstruction image within the effective DOF modulation www.nature.com/scientificreports/    www.nature.com/scientificreports/ loaded into the SLM, fixing the channel 2 object at the 0 mm position and moving the channel 1 object within − 6 mm ≤ d 0 ≤ 0 mm in steps of 1 mm to record the object hologram OHs of the two-channel object at different axial positions, and the holograms of the PSHs and OHs are shown in Fig. 7a and b. The single exposure PSH of the pinhole at the 0 mm position was obtained separately from the OHs at different axial positions using NLR method with a modulation parameter of p = − 0.3 and o = 0.8 to obtain a reconstructed image. The corresponding OHs and object reconstruction images were obtained in the same way for different m. The results are shown in Fig. 7c. It can be seen that different m can achieve simultaneous imaging of objects at different depth planes in the I-COACH system. Thus, by adjusting the AM-CPM, the I-COACH system can quickly achieve DOF imaging of samples with different thicknesses.

Discussion
In this paper, we have proposed a design method for implementing a depth-of-field phase mask AM-CPM for the I-COACH system to further extend the imaging performance of the system. The source of axial defocus characteristics is demonstrated in I-COACH system and the defocus compensation function is determined, thus the AM-CPM with axial defocus aberration correction capability is proposed. The simulations and experimental results fully demonstrate the ability of the proposed method to extend the DOF by adjusting the parameters m and α to obtain different AM-CPM, thus enabling the I-COACH system to image samples at different depths rapidly and simultaneously. In conclusion, the proposed method improves the imaging performance of the I-COACH technology, and provides an important alternative method for further promoting the practical application in synthetic aperture imaging, scattered media imaging, microscopic imaging, and other fields. www.nature.com/scientificreports/

Methods
The optical setup of the I-COACH system with DOF extension and the coordinates of each plane are shown in Fig. 8. The object or pinhole is illuminated critically with incoherent light to guarantee complete incoherence between any two points on the object 10 . The diffracted light of the object is collected and collimated by the lens L 0 with focal length f L0 located at a distance of d 0 from the object. Then, the light is polarized by a polarizer P along the active axis of the SLM and totally modulated by a mask which is synthesized by a CPM, an annular multi-focal coded phase mask and a quadratic phase mask (QPM) with focal length f slm loaded onto the SLM. The light modulated by the SLM is recorded by a charge coupled device (CCD) located at a distance of d s from the SLM. In the general I-COACH setup, the light scattered from the mask is projected on the sensor plane by the QPM to satisfy the Fourier-transform relations between the CPM plane and the CCD plane, which means to d s = − f slm in the optical setup. Hence, the system is a regular lens imaging system when the SLM is only loading QPM, and an I-COACH system when a mixed mask is loaded. When a point object is illuminated in the above optical setup, the light field of the point object on the SLM and CCD planes given as respectively, where ℑ 2D is the 2D Fourier transform. C 1 = −1 2 d i d 0 , C 2 = 1 2 d s , r i = x i , y i and r s = x s , y s are transverse location vector on the SLM and CCD, respectively. r 0 = x 0 , y 0 and r = x ′ , y ′ are transverse location vector on the object and rear plane of lens L 0 , respectively. Q is a quadratic phase function and given by Q(ar) = exp iaπ −1 x 2 + y 2 . A(r) is the aperture function of the rear plane of the lens L 0 , where the value less than the lens L 0 diameter R is 1 otherwise 0. Φ is a phase synthesized by a CPM, an AM-CPM and a QPM. According to the above equations and optical setup, it can be seen the existence of quadratic function Q r which is defined as axial defocus aberrations make the relationship between the object and image plane not satisfy the conjugate relationship of a well-focused of the lens. Therefore, the distance d 0 of point-object varies along z-axis will result in imaging defocus which should be compensated for extending imaging DOF. Unconventionally, we proposed a method of annular multi-focal coded phase mask that is multiplexed on the imaging system pupil plane to compensate simultaneously the various of defocus aberration. In our proposal, the axial defocus compensation function ξ is defined as, where m and α are the designed parameters which can be setting respectively to compensate the different defocus. Let the conjugate relationship between axial defocus aberrations and axial defocus compensation function be fulfilled, and they satisfy the relationship of Eq. (4) at the same time, In Eq. (4), the parameters d i and f L0 are determined according to the system, and d 0 is related to the axial distribution depth of the 3D object, so that different axial locations d 0 correspond to different m and α. In order to achieve continuous correction of axial defocus aberration within the depth of field d 0 , a space division multiplexing method is designed to divide the phase mask into annular equal-area zones. The ring zones with different focal lengths achieve the correction of the corresponding the different axial defocus aberrations. Thus, according to Eq. (4), the AM-CPM is designed as, (5a) AM−CPM= arg exp j2π −1 αf (r i ) , In other words, f (r i ) realizes the different ring zones with different focal lengths, and αf (r i ) works together to realize further correction for different focal lengths with different ring zones. The adjacent annular zones of the AM-CPM have axisymmetric defocus aberration compensation functions to extend the DOF of the imaging system. Therefore, by modulating the parameters m and α, The AM-CPM with different DOF extension capabilities can be obtained.

Data availability
The datasets used and/or analyzed during the current study available from the corresponding author on reasonable request. (5b) f (r i ) =      n 1 γ 2 , n 1 m ≤ γ 2 < 2n 1 + 1 2m