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

Expanded field of view frequency-selective incoherent holography by using a triple-beam setup

Open Access Open Access

Abstract

We propose a new, to the best of our knowledge, method of incoherent optical frequency selection called three-pack frequency-selective incoherent holography. Compressed holography is reconstructed using phase shift intercepts and spatial transfer function convolution in the form of separation without loss of magnification or resolution. The frequency-selective reconstruction process removes the conjugate and DC terms along with the interception of the object wave. This work attempts three-dimensional reconstruction and selected-frequency phase extraction of axial slices in submicron steps, and the experimental results show the potential of the proposed method in areas such as compressed holography, extended field of view, and slice tomography.

© 2022 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement

1. Introduction

The simultaneous capture of several complex wavefronts by a single CCD is known as multiplexed holography. Each wavefront is holographically encoded with a separate message. These intricate wavefronts are compressed, and all of them are captured simultaneously. Multiplexing with microscopic amplification allows digital microscopy to have significantly improved performance in terms of wavefront acquisition and field of view (FOV). Multiplexed micro-holography can therefore be used for a variety of purposes. For instance, non-scanning holographic lamination [14], depth-of-field multiplexing [59], and field-of-view multiplexing [1016]. Holography is a method of recording and reconstructing the amplitude and phase of wave fields, invented by Dennis Gabor in 1948 [17]. During this period, researchers used wet chemical and thin film materials for recording holograms. It was not until 1967 that the concept of digital holography was first introduced by J. W. Goodman and R. Lawrence [18]. It is through digital processing that objects are digitally reconstructed. From then on, the era of digital holography began.

The laser has strong coherence for the development of holography, and the interferometer makes it simple to create holograms. Coherent optical is the first area where multiplexed holography initially debuted. Frenklach et al. [19] introduced an interferometric technique in 2014, referred to as interferometry with a triple imaging area, in which the camera simultaneously captures four holograms, tripling the amount of quantitative data it gathers during a single exposure. Six off-axis holograms were combined into one multiplexed off-axis hologram by Rubin et al. in 2016 [20] without sacrificing resolution or magnification. Additionally, it is possible to recreate it using spatial filtering. The use of spatial bandwidth is significantly improved. In 2017, Dardikman et al. [21] positioned two pairs of complicated wave fronts with conjugates in the area inhabited by the DC term. In one pass, eight-pack holograms (8 PH) were reused. It is necessary to acquire two phase-shifted holograms in order to eliminate the DC term. W. Zhang et al. [22] produced the first two-pack multiplexed frequency-selective (TPFS) holography in the area of incoherent light in 2022. Longitudinal separation of on-axis overlapping spectra was achieved by combining the phase shift method of fresnel correlation holography. Without sacrificing resolution, TPFS doubled FOV. The phase of 3D objects may be reconstructed more effectively by using multiplexed holographic phase extraction methods. Currently, a maximum of eight holograms can be extended from the coherent optical multiplexing that is experimentally feasible. The hologram's frequency domain reveals that the spectrum of the eight off-axis wave fronts already fills the two-dimensional frequency domain when the cutoff angular frequency necessary for the multiplexing method is four times the sample's maximum angular frequency. If there are more off-axis holographic waves in this situation, spatial filtering in the two-dimensional frequency domain will not be able to pick the spectrum of extra conjugate terms. In the coherent optical domain, this restricts the number of multiplexed holographic envelopes. Fresnel Incoherent Correlation Holography (FINCH) [23], a self-interference reference holography that defies the Lagrangian invariance requirement of conventional optical systems, was proposed by Rosen et al. in 2007. Compared to typical coherent imaging systems with similar numerical apertures (NA), FINCH's optical transfer function (OTF) is twice as broad and reacts uniformly to all OTF frequencies [24,25]. Because of this, FINCH has a resolution that is around two and one and a half times greater than analogous coherent and incoherent imaging systems, respectively. Hence, in incoherent multiplexed frequency-selective holography, the self-interference holography principle of FINCH is necessary. In addition, the theory of multiplexed frequency-selective holography provides a good foundation for the advancement of frequency-selective slice laminography and holographic FOV expansions in the incoherent field.

In this paper, we extend the concept of frequency selection experimentally to the field of incoherent three-pack holography (3 PH). This work not only realized the incoherent three-pack object wave frequency-selective reconstruction. It is also the first demonstration of incoherent optical frequency-selective slice tomography. We then present the design, experimental studies, and results of a triple-beam frequency-selective Fresnel incoherent correlated holography (TBFS-FINCH) structure.

2. Methodology

In triple-beam frequency-selective Fresnel incoherent correlated holography (TBFS-FINCH), three object waves leaving the traditional microscope configuration are introduced into the SLM loaded with a phase mask (PM). Due to the interferometric properties of the single FINCH, the three beams containing different specimen information interfere with each other to form a mixed hologram. In the proposed system, there are two cases of mixed-stack holograms. The first is a triple-wave interference hologram, which is formed when the optical path difference (OPD) between any two of the three branches of the object wave is less than the temporal coherence length (TCL). If the OPD between the object wave is larger than the TCL, it is the second case, where the intensities of the three holograms are superimposed and do not contain inter-correlation terms. In the following, we discuss the frequency selection options for the first scenario. This is because the three-branch interference already includes the case of holographic superposition. The multiplexed holography received by the camera evolves into an intensity superposition as soon as the mixed interferometric correlation term is zero. So the reconstruction for the second scenario can be included in the resolution for the first scenario. According to Fresnel diffraction theory, the complex amplitude distribution of the three object waves reaching the sensor surface through the imaging system can be represented separately. We start with the example of specimen O.

Based on the first branch in Fig. 1, in the case of a point object with the amplitude $\sqrt {{I_s}}$ located at $({\overline r _s},{z_s}) = ({x_s},{y_s},{z_s})$, the light diffracted from the object reaches the left surface of the microscope objective (MO) with complex amplitude is expressed as:

$$\tau (x,y;{r_s},{z_s}) = \sqrt {{I_s}} {C_1}L(\overline r /{z_s})Q(1/{z_s}),$$
where $L(\overline s /z) = \exp [i2\pi {(\lambda z)^{ - 1}}({s_x}x + {s_y}y)]$ and $Q(a) = \exp [i\pi a{\lambda ^{ - 1}}({x^2} + {y^2})]$ are the linear and the quadratic phase functions, respectively, in which $\lambda$ is defined as the central wavelength. ${C_1}$ is a complex valued constant dependent on the position of the point source. The complex amplitude immediately after MO is given as $\tau (x,y,{r_s},{z_s}) \cdot Q( - 1/{f_{mo1}})$. The diffractive optical element $r(x,y) = Q( - {\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 a}} \right.}\!\lower0.7ex\hbox{$a$}}){e^{i\theta }}$ is shown on a spatial light modulator (SLM) oriented at a small angle with respect to the optical axis. The SLM is polarization sensitive, and the beams incident on the SLM need to be modulated by the linear polarizer $L{P_1}$. For the first branch, LP1 leaves half of the incident light unmodulated by the SLM. At this point, the SLM is equivalent to a reflector, and the complex amplitude arriving at the sensor surface is
$${D_R}(x,y;\overline {{r_s}} ,{z_s}) = \tau (x,y,{r_s},{z_s}) \cdot Q( - 1/{f_{mo1}}) \ast Q(1/({d_2} + {z_h})).$$

 figure: Fig. 1.

Fig. 1. TBFS-FINCH optical multiplexing system. BF, bandpass filter; $L{P_1}$, linear polarizer 1; $L{P_2}$, linear polarizer 2; BS, beam splitter; MO, microscope objective, ${\times}$ 4, 0.1-NA, ${f_{mo1}}$=7 mm; HWP, half-wave plate; M, mirror; SLM, spatial light modulator; CCD, charge coupled device.

Download Full Size | PDF

$L{P_1}$ similarly causes half of the incident beam to be modulated by the PM of the SLM, and the modulated wave at the same position is

$${D_O}(x,y;{\overline r _{s,\theta }},{z_s}) = \tau (x,y,{r_s},{z_s}) \cdot Q( - 1/{f_{mo1}}) \ast Q(1/{d_2}) \cdot Q( - 1/a){e^{i\theta }} \ast Q(1/{z_h}),$$
where ${\ast} $ denotes a two-dimensional convolution. $\theta$ denotes the phase shift for complex amplitudes formed with different PMs. The initial phase of PM can be generated randomly, but the phase shift between two PMs must be guaranteed to be $\pi$. Following mathematical justification presented in [23], the recorded intensity over the CCD plane is given by,
$${I_k}(\overline {{r_0}} ;{\overline r _{s,\theta }},{z_s}) = {|{{D_R}(x,y;\overline {{r_s}} ,{z_s}) + {D_O}(x,y;{{\overline r }_{s,\theta }},{z_s})} |^2} = {I_k}(\overline {{r_0}} - \frac{{{z_h}}}{{{z_s}}}{\overline r _{s,\theta }};0,{z_s}).$$
We assume that the MO collimates the beam to the SLM surface. The lateral magnification of the system does not take into account ${d_2}$. The right-hand side of the Eq. (4) represents the intensity received by the sensor as a shifted version of the intensity response of a point object on the optical axis, where the shifted distance is ${z_h}/{z_s}$. k denotes the number of holograms formed by each of the three single channels and the number of possible interferences of the two-channel waves with each other. For a single channel, each point source is only spatially coherent with itself. For multiplexed channels, the OPD produced by the BS beam splitting of quasi-monochromatic light satisfies the interference conditions to produce a hologram. The hologram recorded by many point sources is therefore only the sum of the contributions from all point sources. That is, the recorded hologram of single branch is
$${H_k}({r_0};{z_s}) = \int\!\!\!\int\!\!\!\int {{I_k}(\overline {{r_0}} - \frac{{{z_h}}}{{{z_s}}}{{\overline r }_{s,\theta }};0,{z_s})d{x_s}} d{y_s}d{z_s}.$$

Obviously, the triple multiplexed hologram is ${H_{mn}} = \sum\limits_{k = 1}^k {{H_k}({r_0};{z_s})}$. (m,n) are the markers of the holograms recorded in the different cases, nothing more, and not parameters.

Using demodulation (i.e., filtering in the Fourier domain), the cross-referenced holograms cannot be effectively separated, so that single phase recovery is in principle impossible. The valuable Fourier terms ${D_R}(x,y;\overline {{r_s}} ,{z_s}) \cdot D_{_O}^{\dagger} (x,y;{\overline r _{s,\theta }},{z_s})$ and $D_{_R}^{\dagger} (x,y;\overline {{r_s}} ,{z_s}) \cdot {D_O}(x,y;{\overline r _{s,\theta }},{z_s})$ are modulated at the spatial carrier frequency of Fig. 2(d) and are concentrated in the centre of the two-dimensional frequency domain. The symbol ${\dagger}$ here stands for complex conjugate. The analysis is performed with respect to the first branch. ${D_O}(x,y;{\overline r _{s,\theta }},{z_s})$ the complex amplitude equation generated by the imaging system for the point object, is available for processing. For the remaining arbitrary branches also the valuable self-referenced signal waves are processed.

 figure: Fig. 2.

Fig. 2. (a) A hologram of TPFS-FINCH. (c) A hologram of TBFS-FINCH. (b) and (d) are the spectrum of (a) and (c). (e) is a reconstruction of the angular spectrum of (a). (f) is a reconstruction of the angular spectrum of (c). (f) is a reconstruction of the power spectrum of (c). (g) and (h) are the intensity distributions of the corresponding colour paths in (e) and (f), respectively.

Download Full Size | PDF

In order to implement triple-beam frequency-selective FINCH (TBFS-FINCH), the corresponding experimental steps are as follows:

  • (1) The specimens are placed at symbols O,S,A in Fig. 1(a) and the fast axis of the HWP is adjusted to be parallel to the direction of polarization of the object wave light.
  • (2) Figure 1(b) is loaded onto the SLM. With the HWP placed on the first branch only, the intensity distribution acquired on the CCD is ${H_{11}}$. By placing the HWP on the second branch only, the intensity distribution acquired on the CCD is ${H_{12}}$. Applying the HWP on the third branch only, the intensity distribution acquired on the CCD is ${H_{13}}$. When placing two HWPs on the first and second branches respectively, the intensity distribution acquired on the CCD is ${H_{14}}$. The two HWPs are placed in the second and third branches and the intensity distribution acquired on the CCD is ${H_{15}}$. The two HWPs are placed in the first and third branches and the intensity distribution acquired on the CCD is ${H_{16}}$.
  • (3) Load Fig. 1(c) onto the SLM. Repeat the recording sequence from step (2) (the PM of the SLM has been replaced). The intensity distributions acquired by the CCD are: ${H_{21}}$, ${H_{22}}$, ${H_{23}}$, ${H_{24}}$, ${H_{25}}$, ${H_{26}}$.

In post-processing, a Fresnel modulation and a straight reference wave are used to remove defects in the optical path. The redundant terms in Eq. (4) are eliminated by a triple frequency selection scheme. The frequency selection equation is shown below.

$$FT\{ O\} = \frac{{FT\{\frac{1}{4}({{H_{23}} - {H_{24}} + {H_{22}} - {H_{26}}} )\} {H_{ - {z_1}}} - FT\{\frac{1}{4}({{H_{13}} - {H_{14}} + {H_{12}} - {H_{16}}} )\} {H_{ - {z_2}}}}}{{2({H_{{z_2} - {z_1}}} - {H_{{z_1} - {z_2}}})}}$$
$$FT\{ S\} = \frac{{FT\{\frac{1}{4}({{H_{23}} - {H_{24}} - {H_{25}} + {H_{21}}} )\} {H_{ - {z_1}}} - FT\{ \frac{1}{4}({{H_{13}} - {H_{14}} - {H_{15}} + {H_{11}}} )\} {H_{ - {z_2}}}}}{{2({H_{{z_2} - {z_1}}} - {H_{{z_1} - {z_2}}})}}$$
$$FT\{ A\} = \frac{{FT\{ \frac{1}{4}({{H_{22}} - {H_{26}} - {H_{25}} + {H_{21}}} )\} {H_{ - {z_1}}} - FT\{ \frac{1}{4}({{H_{12}} - {H_{16}} - {H_{15}} + {H_{11}}} )\} {H_{ - {z_2}}}}}{{2({H_{{z_2} - {z_1}}} - {H_{{z_1} - {z_2}}})}}$$
where ${H_{{z_i}}}(u,\nu ,{z_i})$ is the spatial transfer function. This can be expressed here as: ${H_{{z_i}}}(u,\nu ,{z_i})=$$\textrm{exp}\left[ {ik{z_i}\sqrt {1 - {{(\lambda {f_x})}^2} - {{(\lambda {f_y})}^2}} } \right],i = 1,2.\,{f_x}$ and ${f_y}$ are the spatial coordinates in the frequency domain. FT denotes two-dimensional Fourier transform. In the experiments, the term ${D_R}(x,y;\overline {{r_s}} ,{z_s}) \cdot D_R^{\dagger} (x,y;{\overline r _{s,\theta }},{z_s})$ used for aberration correction was obtained from a reference hologram that was used in the absence of a sample using the so-called double exposure method [26,27]. The angular spectrum reproduction requires the ${H_{{z_i}}}(u,\nu ,{z_i})$ and the object wave spectrum to be correlated to yield: $S(x,y,z) = IFT\{ FT\{ O(x,y)\} {H_z}(u,v,z)\}$.

3. Experiments

The principle of the TBFS-FINCH was demonstrated experimentally using the digital microholographic setup shown in Fig. 1(a). The xenon lamp (CET-TCX250, 250W) is mounted on the left side of the entire illumination channel. Narrow bandwidth filter $\varDelta \lambda$ = 20 nm with a central wavelength of 633 nm. The plano-convex lens (CL) allows the conversion of point light sources into parallel beams. The effective focal length of the biconvex lens L is 50.8 mm and the light from the incoherent light source is focused by the lens L onto the object in the critical illumination configuration. Three pairs of identical Mo and USAF 1951 charts were placed in three branches in order to strictly magnify the information of the specimen under test by a factor of four. The light aggregated by the BS is polarized by the linear polarizer $L{P_1}$ oriented to an angle of ${45^\circ }$ with respect to the active axis of SLM (FSLM-2K70-VIS, CAS Microstar, 1920 ${\times}$ 1080 pixels, 8 $\mu m$ pixel pitch, phase-only modulation). The mask, calculated using the Fresnel diffraction principle, is displayed on the SLM in 1920 ${\times}$ 1080 pixels. About half of the light intensity is modulated by the mask, the remaining beam is unaffected and reflected. As for the phase shift, the frequency selector combined with the filtering algorithm uses a two-step phase shift, corresponding to a phase shift difference of $\pi$ between the two PMs displayed in the SLM to eliminate twin images and zero-order terms. In order to enable only similarly oriented components of the modulated and unmodulated waves to propagate and cause interference at the sensor plane (MER-502-79U3M/C, 2448 ${\times}$ 2048 pixels, 3.45$\mu m$ pixel size, 79 fps capturing speed), a second polarizer $L{P_2}$, positioned after the SLM, is orientated at 45 degrees with regard to the active axis of the SLM. The distance ${z_s}$ is variable, and the distance ${z_h}$ = 15 cm between the SLM and the camera is aligned with the target so that elements 3,4 and 5 of group 2 of the USAF charts can be imaged onto the camera.

To further experimentally analyse the axial resolution of the TBFS in real scenarios, three planar objects were constructed at same relative distances using USAF resolution charts mounted on channels 1,2 and 3 respectively.

The two-pack incoherent frequency-selective holography technique has been demonstrated in a previous study [22]. This technique was shown to reconstruct the object waves of both branches without loss of resolution. To demonstrate that the 3PH and 2PH in the incoherent domain are equally capable of frequency selection. The effect of mixing and superimposing holographic experiments of elements 3 (5.04 lp/mm), 4 (5.66 lp/mm) and 5 (6.35 lp/mm) of group 2 of the USAF chart is shown in Fig. 2. Figure 2(e) is a reconstruction of the angular spectrum of Fig. 3(a). The holograms of the overlapping numbers 3 and 5 were separated without eliminating the coherent background. Figure 2(g) shows the intensity distribution of the corresponding path in (e). The background intensity is around 80 arbitrary unit (a.u.) and the intensity of the figures is around 90 a.u. This is a capability that the two-pack frequency-selective-FINCH (TPFS-FINCH) setup possesses and has been demonstrated previously in the literature. The hologram of the experiment we performed with the setup TBFS-FINCH according to Fig. 1 is shown in Fig. 2(c). The background Fresnel diffraction ring in this hologram is still visible, indicating that the hologram is a superposition of three self-interfering holograms. Similarly, 3PH is available for spectral selection without loss of resolution. Figure 2(f) shows the power spectrum reconstructed by stripping elements 3, 4, and 5 of group 2 of the USAF chart from the corresponding branches. Figure 2(h) presents the intensity distribution of the three colour paths in (f). Since the graph (f) is reconstructed as a power spectrum, i.e., as the square of the inverse Fourier transform, the background intensity is essentially zero. It is easy to discern the contour information of the object. Figure 2(f) also demonstrates that the TBFS-FINCH has a triple field of view, which is advantageous in situations where detector resolution is insufficient. In the frequency selection scheme for axial holographic superposition, the ratio of the cut-off corner frequency to the maximum corner frequency is 1. The calculation of the wavefront acquisition efficiency for intensity superposition shows that the efficiency of arbitrary n-branch multiplexed holographic frequency selection is n/8. For a on-axis FINCH, a single acquisition of an object wavefront typically requires a three-step phase shift technique to eliminate the conjugate and DC terms. The acquisition of n wavefronts requires 3n holograms. For the setup described in this paper, the improvement in wavefront acquisition efficiency is a factor of 9/8, breaking the limit of less than 1 for the on-axis holographic frequency selection technique.

 figure: Fig. 3.

Fig. 3. (a) is a hologram of TBFS-FINCH, containing paramecium and water cotton cell. (c), (d), (e) are three-branch inverse Fourier reconstructions of (a). (f)-(h) are the intensity figures of (c)-(e). (b) is the intensity distribution of (f)-(h) corresponding to the colour path.

Download Full Size | PDF

The TBFS-FINCH setup also has excellent capabilities in biomicroscopic imaging. Many new optical systems [2830] use microlens arrays to split a complete laser wavefront into many tiny parts in space, effectively superimposing a number of images on the camera target surface. Due to the limitations of the imaging equipment’s field of view, there are multiple biological cell images that are multiplexed on a single camera. The ability to actively select the desired cells for individual analysis is extremely beneficial in this case to exclude extraneous interference factors. In Fig. 3(a), we have taken multiplexed holograms of a water cotton cell and two paramecium of different sizes using a frequency-selective setup. Figure 3(c),(d),(e) shows the inverse Fourier patterns reconstructed by frequency selection from Eqs. (6)-(8) in Section 2, while also feeding back the object wave information of the three branches. The results show that the three spectral stripping reconstructions of the triplexed holograms have been successfully performed. (f)-(h) are intensity figures corresponding to (c)-(e). The intensity distributions of the selected lines S1, S2, and S3 in (f)-(h) are shown in Fig. 3(b). It is demonstrated that the reconstructed cell images after frequency selection can be highlighted in the complex surrounding environment. The incoherent optical holographic frequency selection technique reconstructs the cellular hologram without destroying the cellular structure. It is valuable for studying the parameters of cellular physiological activity.

Figure 4 provides the average cross section of a 300 lp/mm grating based on the reconstruction of the hologram obtained with the TBFS setup. The grating visibility is 0.84, and the frequency selective setup still provides excellent lateral resolution compared to the grating visibility of conventional imaging, FINCH, and COACH in the literature [23,24,31].

 figure: Fig. 4.

Fig. 4. (a) a hologram of TBFS-FINCH. (b) The reconstruction of grating. (c) Average cross-sections of the 300 lp/mm grating.

Download Full Size | PDF

It is noticeable that the TFBS-FINCH setup demonstrates notable slicing capabilities, which we demonstrate with three paramecium specimens. As shown in Fig. 5, with the first two branches being the front and back of the 3D object and the third being the internal section of the 3D object, the specimens are placed on the stepper motor and moved in 0.3 mm steps. The stripped spectrum inverse Fourier transform (IFT) reconstruction of multiplexed holograms using a frequency selection algorithm is used in the experiment. As can be seen from the reconstructed figure, as the focus position of the MO is moved axially over the specimen, the reconstructed section also corresponds to a slice of the corresponding step inside the 3D object. An intensity distribution was made for the red line cross section in Fig. 5, indicating that the internal refractive index and profile of the object are different at different locations in the axial direction, as shown in Fig. 6.

 figure: Fig. 5.

Fig. 5. Experimental comparison of the results of the third path reconstruction when the position of the third branch relative to the focal point of the MO ($\Delta d$) is varied from 0 to 1.8 mm in steps of 0.3 mm.

Download Full Size | PDF

 figure: Fig. 6.

Fig. 6. Intensity distribution of the red line segments corresponding to the step IFT reconstruction in Fig. 5.

Download Full Size | PDF

To test the numerical focusing capability of the setup, we made the aerosol generator spray alumina particles so that they spread around the sample location in the three branches. The multiplexed aerosol hologram captured on the CCD is shown in 7(a). The ideal plane for focused reconstruction can be determined using an autofocus technique based on weighted spectral analysis [32] for the microscopes discussed here. The phase variation brought on by aberrations during phase extraction is not insignificant. Here, we employ the double exposure method to solely display the object’s phase distribution. The reconstructed phase has been unwrapped, and the background phase noise and Fresnel mask geometry tilt-related phases have been effectively eliminated, as can be seen in Fig. 7(e). The phase distribution was converted to thickness distribution using $d = {\raise0.7ex\hbox{${\Delta \phi }$} \!\mathord{\left/ {\vphantom {{\Delta \phi } {[({n_a} - {n_r})k]}}} \right.}\!\lower0.7ex\hbox{${[({n_a} - {n_r})k]}$}}$, where ${n_a}$=1.7659 and ${n_r}$=1.0003 [33] are the constant mean refractive indices of the monochromatic light at 633 nm and the aluminium oxide aerosol's surrounding medium (in this example, air), respectively.

 figure: Fig. 7.

Fig. 7. (a) a aerosol hologram of TBFS-FINCH. (b), (c), (d) are three-branch reconstructions of (a). (e) is the unwrapped phase of (c). (f) Three-dimensional thickness reconstruction of partially phased. (g) Corresponds to the intensity distribution of the dashed line. (f) Local magnification curve in (g).

Download Full Size | PDF

Figures 7(b)-(d) are based on frequency-selected reconstructions under TBFS-FINCH, representing aerosol lateral cross sections at different axial distances from MO. This has enlightening implications for the determination of aerosol concentrations within a certain space. Figure 7(f) shows the aerosol thickness of a partial phase conversion. The holographic technique can record not only the amplitude of the aerosol particles but also the phase diagram of the particles. The average intensity of the red curve in Fig. 7(g) is approximately 160 a.u. It indicates that the unmodulated background in the reconstructed map of (c) is still present. The average intensity after phase unwrapping is approximately 35 a.u., eliminating the background noise, as shown in Fig. 7(h). We expanded the true phase of the aerosol particles using a double-exposure unwrapping algorithm. The thickness converted by the particle phase reveals the exact size of the particles, which range from approximately 1 $\mu m$ to 10 $\mu m$. Figure 7 also confirms that the TBFS setup has the capability to quantitatively visualize phase. Clearly, technology that can accurately record the phase of an object has important implications for research in the field of holographic tomography.

4. Discussion and summary

In this study, we report on the implementation of incoherent multiplexed frequency-selective holography using a triple-beam setup. Our primary contributions in this area are the first demonstrations of frequency-selective reconstruction of incoherent three-packet object waves and optical frequency-selective slice tomography. It is not a conventional multi-camera independent exposure phase shift reconstruction [3436] or inverse filtered projection laminar reconstruction [37,38]. The reflective SLM can be reloaded twice as needed with the phase mask without moving the camera. The proposed TBFS technique is based on only a few components in a simple and static configuration, allowing it to be implemented as a stand-alone module for fast and accurate quantitative phase imaging (QPI). Compared to other multiplexed imaging techniques, TBFS not only improves efficiency by reducing resolution, but also extends FOV.

Coherent optical frequency-selective holography [19,20,21] is being studied and cited by an increasing number of researchers due to its vastly improved performance in terms of wavefront acquisition efficiency and FOV. So far, the number of multiplexed waves of coherent light has been experimentally achieved up to 8 PH. Coherent optical multiplexing differs from TBFS in theory in that it employs a dual image area interferometry technique. Using this technique, two interferograms are projected onto the camera, each of which captures a distinct area of the sample or a complicated wave front for a separate sample. Coherent optical multiplexing allows for the reconstruction and splicing of two or more FOVs using spatial filters since the FOVs do not overlap in the spatial frequency domain. Equations (6)-(8) demonstrate how the TBFS technique used in this study vertically shears the matching spectrum using a spatial transfer function with phase shift. Compared to the TPFS scheme [22], the TBFS scheme deterministically achieves an extension from two-branch multiplexed frequency selection to three-branch. It illustrates the theoretical possibility of extending the pack number from two to infinite dimensions for incoherent optical frequency selection.

Depending on the form of spectral filtering, the wavefront acquisition efficiency of TBFS is lower than that of the coherent light multiplexing holographic frequency selection scheme. Since the temporal coherence length of incoherent light is very small, in the sub-micron order, it is difficult to achieve off-axis holography. Therefore, incoherent holography is mainly used in on-axis holography, which requires the use of phase shift techniques. For a single FINCH, wavefront acquisition usually requires a three-step phase shift to remove the conjugate and zero-level terms. The TBFS scheme improves the wavefront acquisition efficiency by 3n/8 compared to conventional on-axis holography. n is the number of object waves to be acquired. Clearly, the larger the value of n, the more efficient the incoherent frequency selection, i.e., the higher the reconstruction time resolution. Multiplexed frequency selection offers more options for the development of holographic encrypted transmissions, quantum encryption [39,40] and so on.

The TBFS technique in this paper shows the capability of QPI on the micrometer scale in aerosol separation reconstruction. The setup has a single-channel multiplexing feature, a feature that reduces phase fluctuations that may originate from ambient noise [41,42]. This opens up the option of using intensity superposition states to extract phases in the field of remote sensing detection. Multi-pack holography enables the compressed transmission of many layers of information, which is useful when detector resolution is insufficient. This implies that it is unnecessary to factor each subwavefront’s capture interval into the camera’s frame rate. Our research supports the premise that several camera exposures are no longer necessary for QPI based on TBFS. Comparing the precision and stability of the TBFS imaging method to that of traditional multiplexed holography [2225], we discovered that there was no degradation (see Fig. 2). Our study paves the way for the development of wide-field super-resolution TBFS-FINCH for biological or metrological imaging, as proven by the USAF1951 resolution charts, paramecium sections, and aerosol dynamics frequency-selective reconstructions.

Funding

National Natural Science Foundation of China (41371336, 52177137).

Acknowledgment

We thank the anonymous reviewers for their comments on this paper.

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. B. W. Schilling, T. C. Poon, G. Indebetouw, B. Storrie, K. Shinoda, Y. Suzuki, and M. H., “Wu,Three-dimensional holographic fluorescence microscopy,” Opt. Lett. 22(19), 1506–1508 (1997). [CrossRef]  

2. C. M. Fang-Yen, W. Choi, Y. Sung, C. J. Holbrow, R. R. Dasari, and M. S. Feld, “Video-rate tomographic phase microscopy,” J. Biomed. Opt. 16(1), 1 (2011). [CrossRef]  

3. D. Jin, R. Zhou, Z. Yaqoob, and P. T. C. So, “Dynamic spatial filtering using a digital micromirror device for high-speed optical diffraction tomography,” Opt. Express 26(1), 428–437 (2018). [CrossRef]  

4. K. Kim, J. Yoon, and Y. Park, “Simultaneous 3D visualization and position tracking of optically trapped particles using optical diffraction tomography,” Optica 2(4), 343–346 (2015). [CrossRef]  

5. P. A. Dalgarno, H. I. C. Dalgarno, A. Putoud, R. Lambert, L. Paterson, D. C. Logan, D. P. Towers, R. J. Warburton, and A. H. Greenaway, “Multiplane imaging and three dimensional nanoscale particle tracking in biological microscopy,” Opt. Express 18(2), 877–884 (2010). [CrossRef]  

6. M. McGuire, W. Matusik, H. Pfister, B. Chen, J. F. Hughes, and S. K. Nayar, “Optical splitting trees for high precision monocular imaging,” IEEE Comput. Grap. Appl. 27(2), 32–42 (2007). [CrossRef]  

7. C. Maurer, S. Khan, S. Fassl, S. Bernet, and M. Ritsch-Marte, “Depth of field multiplexing in microscopy,” Opt. Express 18(3), 3023–3034 (2010). [CrossRef]  

8. P. Ferraro, M. Paturzo, P. Memmolo, and A. Finizio, “Controlling depth of focus in 3D image reconstructions by flexible and adaptive deformation of digital holograms,” Opt. Lett. 34(18), 2787–2789 (2009). [CrossRef]  

9. J. Kostencka, T. Kozacki, and K. Lizewski, “Autofocusing method for tilted image plane detection in digital holographic microscopy,” Opt. Commun. 297, 20–26 (2013). [CrossRef]  

10. S. K. Mirsky and N. T. Shaked, “First experimental realization of six-pack holography and its application to dynamic synthetic aperture superresolution,” Opt. Express 27(19), 26708–26720 (2019). [CrossRef]  

11. P. Girshovitz, I. Frenklach, and N. T. Shaked, “Broadband quantitative phase microscopy with extended field of view using off-axis interferometric multiplexing,” J. Biomed. Opt. 20(11), 1 (2015). [CrossRef]  

12. N. Rotman-Nativ, N. A. Turko, and N. T. Shaked, “Flipping interferometry with doubled imaging area,” Opt. Lett. 43(22), 5543–5546 (2018). [CrossRef]  

13. P. Girshovitz and N. T. Shaked, “Doubling the field of view in off-axis low-coherence interferometric imaging,” Light: Sci. Appl. 3(3), e151 (2014). [CrossRef]  

14. D. Roitshtain, N. A. Turko, B. Javidi, and N. T. Shaked, “Flipping interferometry and its application for quantitative phase microscopy in a micro-channel,” Opt. Lett. 41(10), 2354–2357 (2016). [CrossRef]  

15. T. Tahara, Y. Awatsuji, K. Nishio, S. Ura, O. Matoba, and T. Kubota, “Space-bandwidth capacity-enhanced digital holography,” Appl. Phys. Express 6(2), 022502 (2013). [CrossRef]  

16. Shan Shan Kou, Laura Waller, George Barbastathis, and Colin J. R. Sheppard, “Transport-of-intensity approach to differential interference contrast (TI-DIC) microscopy for quantitative phase imaging,” Opt. Lett. 35(3), 447–449 (2010). [CrossRef]  

17. D. Gabor, “A new microscopic principle,” Nature 161(4098), 777–778 (1948). [CrossRef]  

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

19. I. Frenklach, P. Girshovitz, and N. T. Shaked, “Off-axis interferometric phase microscopy with tripled imaging area,” Opt. Lett. 39(6), 1525–1528 (2014). [CrossRef]  

20. M. Rubin, G. Dardikman, S. K. Mirsky, N. A. Turko, and N. T. Shaked, “Six-pack off-axis holography,” Opt. Lett. 42(22), 4611–4614 (2017). [CrossRef]  

21. G. Dardikman and N. T. Shaked, “Is multiplexed off-axis holography for quantitative phase imaging more spatial bandwidth efficient than on-axis holography?” J. Opt. Soc. Am. A 36(2), A1–A11 (2019). [CrossRef]  

22. Wanbin Zhang, Baosheng Li, Chunlin Shi, and Jianquan Li, “Two-pack frequency-selective incoherent holography by using a dual-beam setup,” Opt. Lasers Eng. 156, 107086 (2022). [CrossRef]  

23. J. Rosen and G. Brooker, “Digital spatially incoherent fresnel holography,” Opt. Lett. 32(8), 912–914 (2007). [CrossRef]  

24. J. Rosen, N. Siegel, and G. Brooker, “Theoretical and experimental demonstration of resolution beyond the Rayleigh limit by FINCH fluorescence microscopic imaging,” Opt. Express 19(27), 26249–26268 (2011). [CrossRef]  

25. J. Rosen and R. Kelner, “Modified Lagrange invariants and their role in determining transverse and axial imaging resolutions of selfinterference incoherent holographic systems,” Opt. Express 22(23), 29048–29066 (2014). [CrossRef]  

26. Pietro Ferraro, Sergio De Nicola, Andrea Finizio, Giuseppe Coppola, Simonetta Grilli, Carlo Magro, and Giovanni Pierattini, “Compensation of the inherent wave front curvature in digital holographic coherent microscopy for quantitative phase-contrast imaging,” Appl. Opt. 42(11), 1938–1946 (2003). [CrossRef]  

27. Tristan Colomb, Jonas Kühn, Florian Charrière, Christian Depeursinge, Pierre Marquet, and Nicolas Aspert, “Total aberrations compensation in digital holographic microscopy with a reference conjugated hologram,” Opt. Express 14(10), 4300–4306 (2006). [CrossRef]  

28. B. Lokesh Reddy, Prakash Ramachandran, and Anith Nelleri, “Compressive complex wave retrieval from a single off-axis digital Fresnel hologram for quantitative phase imaging and microlens characterization,” Opt. Commun. 478, 126371 (2021). [CrossRef]  

29. Xiao-Wen Cao, Yi-Ming Lu, Hua Fan, Hong Xia, Lei Zhang, and Yong-Lai Zhang, “Wet-etching-assisted femtosecond laser holographic processing of a sapphire concave microlens array,” Appl. Opt. 57(32), 9604–9608 (2018). [CrossRef]  

30. Melania Paturzo, Vito Pagliarulo, Vittorio Bianco, Pasquale Memmolo, Lisa Miccio, Francesco Merola, and Pietro Ferraro, “Digital Holography, a metrological tool for quantitative analysis: Trends and future applications,” Opt. Lasers Eng. 104, 32–47 (2018). [CrossRef]  

31. A. Vijayakumar, Yuval Kashter, Roy Kelner, and Joseph Rosen, “Coded aperture correlation holography–a new type of incoherent digital holograms,” Opt. Express 24(11), 12430–12441 (2016). [CrossRef]  

32. Y Guo, X Chen, and T. Zhang, “Robust phase unwrapping algorithm based on least squares,” Opt. Lasers Eng. 63, 25–29 (2014). [CrossRef]  

33. Rolf Apetz and Michel P. B. van Bruggen, “Transparent alumina: a light-scattering model,” J. Am. Ceram. Soc. 86(3), 480–486 (2003). [CrossRef]  

34. Carlos Trujillo and Jorge Garcia-Sucerquia, “Phase-shifting digital holographic microscopy by using a multi-camera setup,” Opt. Lett. 42(23), 4841–4844 (2017). [CrossRef]  

35. Firas Abedi, You Yang, and Qiong Liu, “Group geometric calibration and rectification for circular multi-camera imaging system,” Opt. Express 26(23), 30596–30613 (2018). [CrossRef]  

36. Luigi Bruno, Gianfranco Bianco, and Massimo A. Fazio, “A multi-camera speckle interferometer for dynamic full-field 3D displacement measurement: Validation and inflation testing of a human eye sclera,” Opt. Lasers Eng. 107, 91–101 (2018). [CrossRef]  

37. Arkadiusz Kuś, “Real-time, multiplexed holographic tomography,” Opt. Lasers Eng. 149, 106783 (2022). [CrossRef]  

38. Julianna Kostencka, Tomasz Kozacki, and Michał Józwik, “Holographic tomography with object rotation and two-directional off-axis illumination,” Opt. Express 25(20), 23920–23934 (2017). [CrossRef]  

39. Z. J. Zhang, H. Yuan, C. M. Xie, and B. L. Ye, “Quantum state sharing with mixing state from six-qubit entangled pure one,” Mod. Phys. Lett. A 35(32), 2050264 (2020). [CrossRef]  

40. Z. J. Zhang, L. Deng, L. Zhang, B. Zhuge, and B.L Ye, “Efficient Tripartite Quantum Operation Sharing with Five-Qubit Absolutely Maximally Entangled State,” Int. J. Theor. Phys. 60(7), 2583–2591 (2021). [CrossRef]  

41. AK Gupta, R Mahendra, and NK Nishchal, “Single-shot phase imaging based on transport of intensity equation,” Opt. Commun. 477, 126347 (2020). [CrossRef]  

42. Y Li, J Di, C Ma, J Zhang, J Zhong, K Wang, X Teli, and J Zhao, “Quantitative phase microscopy for cellular dynamics based on transport of intensity equation,” Opt. Express 26(1), 586–593 (2018). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

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

Fig. 1.
Fig. 1. TBFS-FINCH optical multiplexing system. BF, bandpass filter; $L{P_1}$, linear polarizer 1; $L{P_2}$, linear polarizer 2; BS, beam splitter; MO, microscope objective, ${\times}$ 4, 0.1-NA, ${f_{mo1}}$=7 mm; HWP, half-wave plate; M, mirror; SLM, spatial light modulator; CCD, charge coupled device.
Fig. 2.
Fig. 2. (a) A hologram of TPFS-FINCH. (c) A hologram of TBFS-FINCH. (b) and (d) are the spectrum of (a) and (c). (e) is a reconstruction of the angular spectrum of (a). (f) is a reconstruction of the angular spectrum of (c). (f) is a reconstruction of the power spectrum of (c). (g) and (h) are the intensity distributions of the corresponding colour paths in (e) and (f), respectively.
Fig. 3.
Fig. 3. (a) is a hologram of TBFS-FINCH, containing paramecium and water cotton cell. (c), (d), (e) are three-branch inverse Fourier reconstructions of (a). (f)-(h) are the intensity figures of (c)-(e). (b) is the intensity distribution of (f)-(h) corresponding to the colour path.
Fig. 4.
Fig. 4. (a) a hologram of TBFS-FINCH. (b) The reconstruction of grating. (c) Average cross-sections of the 300 lp/mm grating.
Fig. 5.
Fig. 5. Experimental comparison of the results of the third path reconstruction when the position of the third branch relative to the focal point of the MO ($\Delta d$) is varied from 0 to 1.8 mm in steps of 0.3 mm.
Fig. 6.
Fig. 6. Intensity distribution of the red line segments corresponding to the step IFT reconstruction in Fig. 5.
Fig. 7.
Fig. 7. (a) a aerosol hologram of TBFS-FINCH. (b), (c), (d) are three-branch reconstructions of (a). (e) is the unwrapped phase of (c). (f) Three-dimensional thickness reconstruction of partially phased. (g) Corresponds to the intensity distribution of the dashed line. (f) Local magnification curve in (g).

Equations (8)

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

τ ( x , y ; r s , z s ) = I s C 1 L ( r ¯ / z s ) Q ( 1 / z s ) ,
D R ( x , y ; r s ¯ , z s ) = τ ( x , y , r s , z s ) Q ( 1 / f m o 1 ) Q ( 1 / ( d 2 + z h ) ) .
D O ( x , y ; r ¯ s , θ , z s ) = τ ( x , y , r s , z s ) Q ( 1 / f m o 1 ) Q ( 1 / d 2 ) Q ( 1 / a ) e i θ Q ( 1 / z h ) ,
I k ( r 0 ¯ ; r ¯ s , θ , z s ) = | D R ( x , y ; r s ¯ , z s ) + D O ( x , y ; r ¯ s , θ , z s ) | 2 = I k ( r 0 ¯ z h z s r ¯ s , θ ; 0 , z s ) .
H k ( r 0 ; z s ) = I k ( r 0 ¯ z h z s r ¯ s , θ ; 0 , z s ) d x s d y s d z s .
F T { O } = F T { 1 4 ( H 23 H 24 + H 22 H 26 ) } H z 1 F T { 1 4 ( H 13 H 14 + H 12 H 16 ) } H z 2 2 ( H z 2 z 1 H z 1 z 2 )
F T { S } = F T { 1 4 ( H 23 H 24 H 25 + H 21 ) } H z 1 F T { 1 4 ( H 13 H 14 H 15 + H 11 ) } H z 2 2 ( H z 2 z 1 H z 1 z 2 )
F T { A } = F T { 1 4 ( H 22 H 26 H 25 + H 21 ) } H z 1 F T { 1 4 ( H 12 H 16 H 15 + H 11 ) } H z 2 2 ( H z 2 z 1 H z 1 z 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.