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

Wide-viewing holographic stereogram based on self-interference incoherent digital holography

Open Access Open Access

Abstract

We propose a holographic stereogram synthesis method which uses holograms that are optically captured by self-interference incoherent digital holography (SIDH) based on a geometric phase lens. SIDH is a promising solution for hologram acquisition under low-coherence lighting condition. A mechanical scanning system is constructed to acquire multiple perspective holograms. Numerical simulations and experimental analyses conducted using high-resolution diffractive optical element demonstrate that the proposed method can produce a wide-viewing hologram which can realize realistic 3D scenarios with depth cues such as accommodation and motion parallax. The future objectives include the implementation of a multiple-camera system for holographic videos.

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

1. Introduction

A holographic display reconstructs optical field and provides complete depth cues. Various approaches have been proposed over several decades to efficiently modulate the optical field and produce high-quality three-dimensional (3D) images [14]. However, a high space-bandwidth product of a spatial light modulator (SLM) must be developed along with an efficient holographic data generation approach to present a realistic 3D visual experience using a holographic display.

The existing hologram generation methods are primarily based on computer-generated holograms (CGHs). The CGHs synthesis algorithm aims to determine the solution to propagation equation in a specific 3D formation, which is an ill-posed problem. CGHs can be categorized into direct and iterative methods. The direct method calculates the optical field from the physical condition using a propagation equation [5]. Various sources of 3D information have been employed, including point [68], layer [911], polygon [1214], and light field [1521]. However, the ideally calculated result does not sufficiently match the actual optical configuration, which leads to degraded image quality. Additionally, the computation time increases with the increase complexity of 3D information. An iterative method such as the Gerchberg-Saxton (GS) algorithm propagates a complex wave back and forth between the image and the SLM planes [2224]. This approach is simple and straightforward, but the iteration requires optimization and incurs a higher computational cost.

Data-driven approaches are one of the latest solutions in modern computer science [2427]. The development in deep neural networks (DNNs) and general-purpose graphics processing units (GPGPUs) has made deep learning as a powerful solution for conventional computer vision problems. Various studies have proposed the application of deep learning in the field of general holographic tasks such as holographic generation and reconstruction. However, machine learning requires numerous datasets for high-precision learning and high computing power to train the models.

However, generating the CGH in a real-world scenario requires additional 3D information. The speed of generation of 3D information from virtual scenarios has rapidly increased owing to the development of computer graphics. However, supplementary devices are required to extract 3D information from the real scenarios. Even if precise 3D information is obtained by using an accurate depth sensing device, complicated calculations are required during the hologram synthesis process.

This paper proposes a method for synthesizing wide-angle holographic stereograms using optically captured holograms. Holographic stereograms are ray-based approaches which convert light fields into holograms by adding phase and depth factors [1922]. They comprise holographic elements (hogels) which consider the wavefield as a summation of a spatio-angular distribution of light rays. A ray-based approach enables the approximation of the optical field, which helps perform effective calculations. Additionally, a holographic stereogram is used as the basic principle of holographic printing [28]. Self-interference incoherent digital holography (SIDH) or Fresnel incoherent correlation holography (FINCH) is a promising incoherent hologram recording method [2936], which splits the light emitted from an object into two waves and modulates each light wave differently. The two light waves are emitted from the same point, they have mutual coherence which enables interference on the image sensor even under low-coherence lighting conditions. SIDH or FINCH can be implemented using SLM [30], birefringence lenses [31], or metalens [32]. This study proposed the implementation of SIDH using a Pancharatnam-berry phase lens or a geometric phase (GP) lens, which modulate light based on the polarization state [3335]. Because the proposed holographic camera has a compact configuration and acquires a phase-shifted hologram with a single exposure, we constructed a rotational scanning holographic camera system. Each directional holographic scenario was multiplied by the corresponding plane carrier wave. To enhance the visibility of synthetic holographic stereograms, holographic post-processing which compensates a wavefront by using optimization, was applied to each hologram. Hence, the proposed method can provide precise depth cues with a smooth parallax and occlusion. The digital reconstruction of a 25k resolution hologram and the optical reconstruction of 25k and 100k high-resolution holograms demonstrate the effectiveness of the angular synthesis approach and realistic 3D scenarios.

2. Proposed holographic stereogram synthesis method

2.1 Holographic stereogram synthesis

The stereogram produces 3D depth cues by composing images of different perspectives, and the holographic stereogram consists of the perspective view by hologram element or ‘hogel’. In this study, the hogels were replaced by angular-multiplexed holograms. Figure 1 presents a schematic of the concept. The finally synthesized hologram is called ‘master-hologram,’ and each hologram obtained through rotation is called ‘sub-hologram.’ Each perspective reproduced by the master-hologram comprises a sub-hologram multiplied by the plane carrier wave corresponding to the capturing direction. The plane carrier wave can be expressed as a k-vector.

 figure: Fig. 1.

Fig. 1. Schematic diagram of holographic stereogram.

Download Full Size | PDF

The proposed holographic stereogram synthesis approach consists of three steps. Firstly, the sub-hologram is acquired. Figure 2 illustrates the sub-hologram capturing scheme and the flow chart of the holographic stereogram synthesis procedure. Figure 2(a) presents the optical configuration of the GP lens based SIDH. In the sub-hologram acquisition, the GP lens splits the incident wave into two waves, based on the circular polarization state. According to the principle of self-interference, these two waves are coherent because they originate from the same point. Additionally, the parallel phase shift method is used to generate the complex hologram since the GP lens simultaneously induces phase shift with the wavefront modulation [3335]. The light passes through the polarizer again after passing through the GP lens, to ensure that both the waves have the same polarization; the phase difference is determined by the relative angle difference between the two polarizers. If the angle variance is four steps with 45° per step, four phase-shifted images with 90° per step are obtained and recombined to the complex hologram without bias and twin image noise. The polarized image sensor consists of the micro-polarizer array on the pixel array, which can acquire the intensity value of the polarization component of 0°, 45°, 90°, and 135°, and each extracted intensity hologram is recombined into the complex hologram.

 figure: Fig. 2.

Fig. 2. (a) Schematic of the optical configuration of geometric phase based SIDH system. Solid and dashed lines illustrate the optical path of the marginal rays emitted from the rotation center, (b) the flow chart of the proposed method.

Download Full Size | PDF

Secondly, free space propagation is performed to the position of the master-hologram corresponding to the rotation center. The sub-hologram must match the geometric configuration of the capturing scheme. Therefore, the angular spectrum propagator is used to prevent the sampling problem [37]. Since the sampling period of the angular spectrum propagation remains constant during the propagation, the propagation result exhibit same magnification with the raw hologram. After the backpropagation, to match the resolution of the master-hologram, the sub-hologram is resized to the resolution of the master-hologram using the bicubic interpolation method.

Finally, the plane carrier wave corresponding to the capturing direction is multiplied and then added to form the master-hologram. The resolution and sampling period of the carrier wave are matched with that of the master-hologram. The field of the master-hologram, which is synthesized from the sub-hologram of $M \times N$ is given as follows:

$${\textrm{U}_{Master}}({x,\; y;{z_c}} )= \mathop \sum \limits_M \mathop \sum \limits_N {U_{Sub}}({x,\; y;{z_c}} )\cdot \textrm{exp}({j({{k_{x;m,\; n}}{x_{sub}} + {k_{y;m,\; \; n}}{y_{sub}}} )} ),$$
where $U({x,\; y;{z_c}} )$. is an optical field at the center of rotation, ${z_c}$., with the sampling point $({x,\; y} )$. . The k-vector of the plane carrier wave $({{k_{x;m,\; n}},\; {k_{y;m,\; \; n}}} )$. is given by:
$$({k_{x;m,n}},\; {k_{y;m,n}}) = \frac{{2\pi }}{\lambda }({\sin {\theta_{x;m,n}},\sin {\theta_{y;m,n}}} ), $$
where $\lambda $ is the wavelength, and ${\theta _{x;m,n}}$ and ${\theta _{y;m,n}}$ are the horizontal and vertical angles of the sub-hologram $({m,\; n} )$, respectively. After the clation, the holographic encoding method is employed to convert complex hologram datbased on the optical reconstruction criteria. In the case of phase modulation, the phase-ane part of the complex hologram is extracted and converted into a binary hologram.

The angular sampling distance of each sub-hologram is determined by its maximum spatial frequency, which is obtained based on the specifications of the image sensor and the Nyquist theorem. The largest spatial frequency is ${f_{max\; }} = 1/2\varDelta x$, which can determine the synthetic interval of the sub-hologram as follows:

$$\sin \mathrm{\Delta }\theta = \frac{\lambda }{{2\mathrm{\Delta }x}} = \lambda \cdot {f_{xmax}} = \frac{{{N_x}}}{{N_x^{\prime}}} \cdot \frac{\lambda }{{\mathrm{\Delta }x^{\prime}}} $$
The synthetic angular period, $\mathrm{\Delta }\theta $, changes depending on the space-bandwidth product of the master-hologram. Therefore, the actual synthetic angular period, $\mathrm{\Delta }{\theta ^{\prime}}$, must be larger than $\mathrm{\Delta }\theta $. Since the bandwidth of the sub-hologram is fixed by the specifications of the image sensor, an empty view occurs in the reconstruction if the sampling interval is larger than the optimal period, $\mathrm{\Delta }\theta $.

2.2 Sub-hologram preprocessing

The proposed holographic stereogram carries holographic information on the plane carrier waves. Therefore, the image quality of the sub-hologram determines the image quality of the master-hologram. The proposed SIDH system is focused on the field-of-view (FOV), depth resolution, and sacrifice beam condition matching of the conventional FINCH, which presents high resolution, to capture a macroscopic object. The FOV of the SIDH system is determined by the ratio of the distances between the GP lens and object and the image sensor and GP lens. However, during the process of adjusting the distance condition to satisfy the FOV, the image quality degrades due to the illumination conditions, spectral bandwidth of the light source, and defects in the GP lens. The illumination condition is directly related to the clarity of the holographic fringe pattern; the narrower the spectral bandwidth of the incoherent light source, the higher the resolution of the self-interference hologram that is generated. Defects in the GP lens considerably affect the phase modulation accuracy. The GP lens is fabricated through the photo-alignment of liquid crystal cells. If the liquid crystal alignment is distorted by an external factor during the process of fabricating the GP lens, the phase modulation efficiency decreases, and it remains as an additional error term of a specific polarization component. This is not removed as a bias or twin image noise in the phase shift hologram, but remains as a residual bias and degrades the quality of the reconstructed image.

Among the various approaches used to calibrate holograms [3840], the black box approach is applied in this study to numerically suppress the noise of the holographic signal. The aforementioned optical aberration factors ultimately construct one aberration term on the sub-hologram. The consequential aberration field is defined and the aberration term modeled as a quadratic surface is determined using a numerical optimization technique.

Figure 3 illustrates the wavefront compensation procedure. While obtaining a phase-shift hologram using polarization, the aberration terms generated by the various factors introduce errors in the polarization component, as shown in Figs. 3(a) and 3(d). Therefore, the numerical aberration compensation process aims to approximate the aberration field as the real and imaginary parts of the hologram. The compensated complex amplitude, ${U_{com}}$., can be expressed as

$${U_{com}}({x,\; y} )= {U_{ori}} - {U_{abe}},$$
where ${U_{ori}}$ is the original complex amplitude, and ${U_{abe}}$. is the aberration complex amplitude. Figures 3(b) and 3(e) present the modeled aberration surface. According to the sparse optimization technique described in Ref. [40], ${U_{abe}}$ is considered for the surface model in the second-order as
$${U_{abe}}({x,\; y} )= {\beta _1}{x^2} + {\beta _2}{y^2} + {\beta _3}xy + {\beta _4}x + {\beta _5}y + {\beta _6},$$
where $\beta $ denotes the fitting coefficient to be estimated. Note that for mathematical simplification, second-order model is used. Second-ord approximation can compensate the most common aberrations such as tilt, astigmatism, and parabolic surface. The extension of the fitting model to a higher order is straightforward. Conventionally, this problem can be solved by the least square method; the regularization technique is used to prevent overfitting. Optimization is performed by a nonlinear programming solver integrated into MATLAB. The optimization process involves identifying the waveont which minimizes the objective function for the input wavefront. The objective loss function, $\mathrm{{\cal L}}$, is
$$\mathrm{{\cal L}} = \mathop {\textrm{argmin}}\limits_{{U_{abe}}} {|{{U_{abe}} - U} |_1} + \mathrm{\Lambda }{|{{U_{abe}}} |_1},$$
where ${|\cdot |_1}$. denotes the ${\ell _1}$. -norm, and $\mathrm{\Lambda }$ denotes the regularization parameter, which gives the sparsity of optimization. This optimization problem is equivalent to minimizing energy in the system. Furthermore, the enforced sparsity prevents the overfitting problem and improves the image quality enhancement. Figures 3(c) and 3(f) present the compensated results obtained by subtracting the aberration field acquired using the proposed method. The overall average value approached zero, and the irregulasurface was corrected accordingly.

 figure: Fig. 3.

Fig. 3. Wavefront optimization procedure applied at real and imaginary parts of raw hologram. (a–c) Real part compensation procedure, (d–f) imaginary part compensation procedure.

Download Full Size | PDF

Figure 4 presents the results of applying wavefront compensation to SIDH. The removal of background noise in the digital reconstruction of the ensated hologram is demonstrated through the comparison with the raw hologram. Particularly, it is eliminated that the window-like noise at the peripheral region which is similar to the ‘ringing’ noise of the coherent image. The focus of the reconstruction is maintained, indicating that the proposed wavefront correction preserves the phase information. This is because the average of the complex data was adjusted to a value close to zero through wavefront compensation. The preserved holographic information leads to a precise accommodation cue during reconstruction. However, there is fluctuation in the phase-angle part after the compensation, as shown in Fig. 4(h). There are two resultant aberration surfaces because the compensation is utilized to the real and imaginary parts of the hologram. Therefore, these two surfaces act independently to compensate the wavefront; the continuity for an arbitrary local phase appears to be lost. Nevertheless, the reconstruction result is not significantly affected by this phase fluctuation.

 figure: Fig. 4.

Fig. 4. Wavefront compensation result comparison and digital reconstruction. (a–d) Raw hologram, (e–h) Compensated hologram (i) digital reconstruction of raw hologram, (j) digital reconstruction of compensated hologram.

Download Full Size | PDF

The brightness of the peripheral region in the digital reconstruction of the compensated hologram was reduced, as shown in Fig. 4(j). The brightness of the image was determined by the aperture limitation. In SIDH and FINCH, the overall hologram is composed as a point spread hologram from a point source. This approach has two main limitations: the optical path length difference and the Nyquist frequency. The optical path length difference is based on the focal length of the GP lens and optical configuration of the objective illumination and must be smaller than the coherence length. The Nyquist frequency is related to the specifications of the image sensor; an aliasing problem occurs if the highest spatial frequency of the point spread hologram is smaller than the Nyquist frequency. Compared with the Nyquist frequency, the optical path difference is the dominant factor of the aperture limitation, which defines the minimum size of the point spread hologram. The aperture limitation owing to the size of the image sensor results in the cropping of the periphery of the hologram, which makes the periphery of the reconstructed image appear dark.

3. Experiment

3.1 Sub-hologram generation

The proposed method was verified by performing a numerical simulation and a proof-of-concept (PoC) demonstration using a diffractive optical element (DOE). A multi-angle rotation system was set up and integrated with a holographic camera based on our previous research [3335]. The GP lens based SIDH can be installed in a motorized scanning system because it is robust against vibrations. The implemented motorized stage has a range of motion of 90° and 45° in horizontal and vertical directions, respectively. Figure 5 illustrates the capture conditions. The holograms of the emissive and the reflective objects placed 50 mm apart from the rotation center, were captured. The emissive target consisted of a LED surface light source and patterned aperture of letters ‘ETRI’ and ‘KHU’, and the reflective target was a white statue. The distance between the rotation center and GP lens was 200 mm for the emissive target and 250 mm for the reflective target. A polarization image sensor was used to implement the single-shot SIDH system, with a resolution of 1024 × 1024 pixels and a pixel pitch of 6.9 $\mathrm{\mu m}$. In our experiment, a 550 nm bandpass filter with a bandwidth of ±50 nm and a GP lens with focal length 265 mm were used. Furthermore, the distance from the GP lens to the image sensor was 8 mm, which determines the field-of-view (FOV) of the holographic camera, which was set as 23.6°.

 figure: Fig. 5.

Fig. 5. (a) Concept of the sub-hologram capture system. The experimental setup of (b) emissive target, (c) reflective target.

Download Full Size | PDF

3.2 Numerical simulation

The digital reconstruction was verified by synthesizing the master-hologram with a resolution of 25,000 × 25,000 and the sampling period of 1 $\mathrm{\mu m}$, with the sub-hologram of the emissive target. It is easier to confirm the holographic signal than the reflective target as the brightness of the incident light determines the quality of the holographic fringe pattern. According to Eq. (3), the angular sampling distance was 1.29°, and the number of sub-holograms was 25 × 25.

Figure 6 presents the spatial frequency domain and digital reconstruction of the master-hologram. The signal distribution was arranged periodically in the frequency domain, as shown in Fig. 6(a). The periodic arrangement of holographic signals in the Fourier domain was obtained through the multiplication of the plane carrier wave of each k-vector. The view corresponding to the angle was reconstructed by selecting a specific area, as shown in Figs. 6(b)–6(g). Each reconstructed view image can provide accommodation according to the propagation equation and a smooth parallax owing to the viewing directions.

 figure: Fig. 6.

Fig. 6. Numerical simulation result of synthesized hologram. (a) Spatial frequency domain, digital reconstruction results when focusing on the (b–d) backward letter ‘KHU’ and (e-g) forward letter ‘ETRI’, viewing from (b, e) -10°, (c, f) center, (d, g) +10°.

Download Full Size | PDF

3.3 Optical experiment

The holographic stereogram made of silica DOE was fabricated using lithography to demonstrate optical reconstruction. Off-axis processing, phase encoding, and binarization techniques were applied to the complex holographic stereogram. Off-axis processing on the angular domain was employed at 20° to avoid the illumination bias. The conjugate signal induced by binary modulation was deleted in the off-axis processing. Subsequently, phase angle encoding and binarization based on zero were performed by applying the inverse Fourier transform.

Holographic stereograms with two resolutions, i.e., 25k and 100k were produced. The 25k hologram synthesis conditions were identical to those in the simulation in terms of the resolution and pixel pitch. The synthesis condition of the 100k high-resolution hologram included a sampling period of 0.5 $\mathrm{\mu m}$, synthesis interval of 0.64°, and perspective of 99 × 75. Figure 7 depicts the fabricated holographic stereogram of a 25k hologram without illumination. The dimensions of the holographic stereogram were 25 mm × 25 mm. The hologram was reconstructed by using tilted LED illumination at an incident angle of 12°. The wavelength was 550 nm and the field of view was 31.9°.

 figure: Fig. 7.

Fig. 7. Fabricated 25k holographic stereogram without illumination.

Download Full Size | PDF

Figure 8 presents each result of the overall synthesis procedure and optical reconstruction comparison. Figure 8(a) depicts the phase-angle part of the compensated sub-hologram which has the resolution of $({{N_x},\; {N_y}} )= ({1024,\; 1024} ).$ Fig. 8(b) depicts the result of applying backpropagation to the rotation center and resizing to the resolution of $({N_x^{\prime},\; N_y^{\prime}} )= ({25,000,\; 25,000} ).$ The bicubic algorithm was used as the interpolation method. Figure 8(c) depicts the plane carrier wave, and Fig. 8(d) is a single perspective sub-hologram which multiplied with the resized hologram data and carrier wave. The master-hologram shown in Fig. 8(e) was obtained by adding all the sub-holograms obtained by repeating this process at all the measured angles. The phase-angle part of the hologram was extracted and binarized to reconstruct the synthesized holographic stereogram, as shown in Fig. 8(f).

 figure: Fig. 8.

Fig. 8. Step-by-step phase-angle part hologram visualization of the proposed method. (a) Raw hologram with resolution of $({{N_x},\; {N_y}} )$, and sampling period of $({\varDelta x,\; \varDelta y} )$, (b) back propagation and resizing result, (c) arbitrary plane carrier wave, (d) plane carrier multiplied result, (e) the result of summation of sub-hologram, (f) binarization of phase-angle part of master-hologram, (g) optical reconstruction without wavefront compensation, (h) backward focused optical reconstruction with wavefront compensation, (i) frontward focused optical reconstruction with wavefront compensation.

Download Full Size | PDF

The effectiveness of the wavefront compensation was demonstrated through a comparison between the compensation adapted proof of concept and the uncompensated one, as shown in Figs. 8(g), 8(h), and 8(i). While the optical reconstruction of the uncompensated data can record the holographic signal, it is disrupted by noise and cannot be observed properly, as shown in Fig. 8(g). Conversely, in the case of the corrected data, clear images and focus cues can be observed, as shown in Figs. 8(h) and 8(i).

Figure 9 presents the optical reconstruction result of the 25k holographic stereogram of the emissive target. The same angle as that of the digital reconstruction result shown in Fig. 6, was observed. The synthesized conditions were a viewing zone of 25 × 25 and an angular distance of 1.29°. The theoretical diffraction angle was 15.9°. However, the observation was conducted from -10° to +10° because the diffraction efficiency decreases as the angle of observation increases, which is the characteristic of DOE. Variation of focus and motion parallax are observed in the results of Fig. 9. The dark regions at the periphery of the letters are induced by overlapping images of different angles.

 figure: Fig. 9.

Fig. 9. Optical reconstruction of the 25k holographic stereogram. (a-c) Backward focus, (d-g) Frontward focus, viewing from (a, d) -10°, (b, f) center, (c, g) 10°.

Download Full Size | PDF

Figure 10 and visualizations 1 and 2 depict the fabricated 100k high-resolution holographic stereogram and the resultant optical reconstruction. The hologram of the reflective target was noisier than the emissive target due to the scattering effect on the surface. The same wavefront compensation was applied to the sub-holograms, and the master-hologram of the 100k resolution was synthesized. The dimensions of the produced holographic stereogram were 50 mm × 50 mm, and the viewing angle of the holographic stereogram was approximately 66.7° from -33° to 33°. However, the diffractive efficiency was relatively low compared with the close view of the central axis. Figures 10(b)–10(g) present the holographic stereogram observed at a viewing angle of approximately 25° where the diffraction efficiency was relatively preserved. The optical reconstruction results presented a natural 3D image even when viewed with the naked eye, as shown in Figs. 10(h) and 10(i). In the case of the periphery, a dim image was observed owing to the diffractive efficiency of the DOE, but the depth cue and occlusion were also confirmed.

 figure: Fig. 10.

Fig. 10. Experimental result of the 100k holographic stereogram result. (a) Fabricated 100k hologram without illumination. (b–d) Backward focus, (e–g) frontward focus, viewing from (b, e) -25°, (c, f) center, (d, g) 25°, (h–i) close view of (c) and (f) (see Visualization 1 and Visualization 2).

Download Full Size | PDF

The proposed method can be used to develop a high-resolution hologram synthesis algorithm which uses holographic cameras. A higher resolution is required to provide a more natural 3D visual experience using holographic displays. However, the current CGHs generation techniques which are based on the FFT algorithm require more computing power as the resolution of the hologram increases to more than 100k. The proposed holographic stereogram approach is a potential solution for high-resolution holograms, and a hologram can be generated without requiring additional 3D information such as a depth sensor. However, the numerical approximation-based solution is used in this case because the noise suppression is required. The characteristics of the optical element must be improved along with the optical structure of the holographic camera to acquire high-quality holographic stereograms.

4. Conclusion

In summary, this study presented a self-interference incoherent digital hologram-based holographic stereogram synthesis method. The method was used to synthesize a high space-bandwidth hologram which excludes additional 3D information for separate CGH generation. Wavefront correction through an optimization technique was applied to SIDH to address the arbitrary noise in the hologram, thereby improving the image quality of the holographic stereogram. The effectiveness of the proposed method was demonstrated through numerical simulation and optical reconstruction, and accommodation and motion parallax were produced. The future objectives include the implementation of a multiple-camera system for holographic videos. Therefore, careful assessment of the holographic camera is required, and the physical limitations of the image sensor must be overcome. We believe that this attempt inspires various fields of the incoherent holographic camera.

Funding

Institute for Information and Communications Technology Promotion (IITP) grant funded by the Korea government (MSIT)(2019-0-00001).

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. S. A. Benton and V. M. Bove Jr, Holographic imaging (John Wiley & Sons, 2008).

2. C. Chang, K. Bang, G. Wetzstein, B. Lee, and L. Gao, “Toward the next-generation VR/AR optics: a review of holographic near-eye displays from a human-centric perspective,” Optica 7(11), 1563–1578 (2020). [CrossRef]  

3. Y. Lim, K. Hong, H. Kim, H.-E. Kim, E.-Y. Chang, S. Lee, T. Kim, J. Nam, H.-G. Choo, J. Kim, and J. Hahn, “360-degree tabletop electronic holographic display,” Opt. Express 24(22), 24999–25009 (2016). [CrossRef]  

4. J. An, K. Won, Y. Kim, J.-Y. Hong, H. Kim, Y. Kim, H. Song, C. Choi, Y. Kim, J. Seo, A. Morozov, H. Park, S. Hong, S. Hwang, K. Kim, and H.-S. Lee, “Slim-panel holographic video display,” Nat. Commun. 11(1), 5568 (2020). [CrossRef]  

5. J.-H. Park, “Recent progress in computer-generated holography for three-dimensional scenes,” J. Inf. Disp. 18(1), 1–12 (2017). [CrossRef]  

6. M. E. Lucente, “Interactive computation of holograms using a look-up table,” J. Electron. Imaging 2(1), 28–34 (1993). [CrossRef]  

7. C. Gao, J. Liu, X. Li, G. Xue, J. Jia, and Y. Wang, “Accurate compressed look up table method for CGH in 3D holographic display,” Opt. Express 23(26), 33194–33204 (2015). [CrossRef]  

8. D. Yasuki, T. Shimobaba, M. Makowski, J. Suszek, T. Kakue, and T. Ito, “Hologram computation using the radial point spread function,” Appl. Opt. 60(28), 8829–8837 (2021). [CrossRef]  

9. J.-S. Chen and D. P. Chu, “Improved layer-based method for rapid hologram generation and real-time interactive holographic display applications,” Opt. Express 23(14), 18143–18155 (2015). [CrossRef]  

10. Y. Zhao, L. Cao, H. Zhang, D. Kong, and G. Jin, “Accurate calculation of computer-generated holograms using angular-spectrum layer-oriented method,” Opt. Express 23(20), 25440–25449 (2015). [CrossRef]  

11. H. Zhang, L. Cao, and G. Jin, “Computer-generated hologram with occlusion effect using layer-based processing,” Appl. Opt. 56(13), F138–F143 (2017). [CrossRef]  

12. H. Kim, J. Hahn, and B. Lee, “Mathematical modeling of triangle-mesh-modeled three-dimensional surface objects for digital holography,” Appl. Opt. 47(19), D117–D127 (2008). [CrossRef]  

13. D. Im, J. Cho, J. Hahn, B. Lee, and H. Kim, “Accelerated synthesis algorithm of polygon computer-generated holograms,” Opt. Express 23(3), 2863–2871 (2015). [CrossRef]  

14. Y.-P. Zhang, F. Wang, T.-C. Poon, S. Fan, and W. Xu, “Fast generation of full analytical polygon-based computer-generated holograms,” Opt. Express 26(15), 19206–19224 (2018). [CrossRef]  

15. M. Yamaguchi, H. Hoshino, T. Honda, and N. Ohyama, “Phase-added stereogram: calculation of hologram using computer graphics technique,” Proc. SPIE 1914(9), 25–31 (1993). [CrossRef]  

16. H. Kang, T. Yamaguchi, and H. Yoshikawa, “Accurate phase-added stereogram to improve the coherent stereogram,” Appl. Opt. 47(19), D44–D54 (2008). [CrossRef]  

17. Q. Y. Smithwick, J. Barabas, D. E. Smalley, and V. M. Bove Jr, “Interactive holographic stereograms with accommodation cues,” in Practical Holography XXIV: Materials and Applications, vol. 7619 (International Society or Optics and Photonics, 2010), p. 761903.

18. K. Wakunami and M. Yamaguchi, “Calculation for computer generated hologram using ray-sampling plane,” Opt. Express 19(10), 9086–9101 (2011). [CrossRef]  

19. H. Zhang, Y. Zhao, L. Cao, and G. Jin, “Fully computed holographic stereogram based algorithm for computer-generated holograms with accurate depth cues,” Opt. Express 23(4), 3901–3913 (2015). [CrossRef]  

20. J.-H. Park and M. Askari, “Non-hogel-based computer generated hologram from light field using complex field recovery technique from wigner distribution function,” Opt. Express 27(3), 2562–2574 (2019). [CrossRef]  

21. H. Zhang, L. Cao, and G. Jin, “Three-dimensional computer-generated hologram with Fourier domain segmentation,” Opt. Express 27(8), 11689–11697 (2019). [CrossRef]  

22. N. Padmanaban, Y. Peng, and G. Wetzstein, “Holographic near-eye displays based on overlap-add stereograms,” ACM Trans. Graph. 38(6), 1–13 (2019). [CrossRef]  

23. R. W. Gerchberg, “A practical algorithm for the determination of phase from image and diffraction plane pictures,” Optik 35, 237–246 (1972).

24. R. G. Dorsch, A. W. Lohmann, and S. Sinzinger, “Fresnel ping-pong algorithm for two-plane computer-generated hologram display,” Appl. Opt. 33(5), 869–875 (1994). [CrossRef]  

25. Y. Peng, S. Choi, N. Padmanaban, and G. Wetzstein, “Neural holography with camera-in-the-loop training,” ACM Trans. Graph. 39(6), 1–14 (2020). [CrossRef]  

26. Y. Rivenson, Y. Wu, and A. Ozcan, “Deep learning in holography and coherent imaging,” Light: Sci. Appl. 8(1), 85 (2019). [CrossRef]  

27. L. Shi, B. Li, C. Kim, P. Kellnhofer, and W. Matusik, “Towards real-time photorealistic 3d holography with deep neural networks,” Nature 591(7849), 234–239 (2021). [CrossRef]  

28. J. Su, X. Yan, Y. Huang, X. Jiang, Y. Chen, and T. Zhang, “Progress in the synthetic holographic stereogram printing technique,” Appl. Sci. 8(6), 851 (2018). [CrossRef]  

29. J. Rosen, S. Alford, V. Anand, J. Art, P. Bouchal, Z. Bouchal, M. U. Erdenebat, L. Huang, A. Ishii, S. Juodkazis, N. Kim, P. Kner, T. Koujin, Y. Kozawa, D. Liang, J. Liu, C. Mann, A. Marar, A. Matsuda, T. Nobukawa, T. Nomura, R. Oi, M. Potcoava, T. Tahara, B. L. Thanh, and H. Zhou, “Roadmap on recent progress in FINCH technology,” J. Imaging 7(10), 197 (2021). [CrossRef]  

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

31. N. Siegel, V. Lupashin, B. Storrie, and G. Brooker, “High-magnification super-resolution finch microscopy using birefringent crystal lens interferometers,” Nat. Photonics 10(12), 802–808 (2016). [CrossRef]  

32. H. Zhou, L. Huang, X. Li, X. Li, G. Geng, K. An, Z. Li, and Y. Wang, “All-dielectric bifocal isotropic metalens for a single-shot hologram generation device,” Opt. Express 28(15), 21549–21559 (2020). [CrossRef]  

33. K. Choi, J. Yim, S. Yoo, and S.-W. Min, “Self-interference digital holography with a geometric-phase hologram lens,” Opt. Lett. 42(19), 3940–3943 (2017). [CrossRef]  

34. K. Choi, K.-I. Joo, T.-H. Lee, H.-R. Kim, J. Yim, H. Do, and S.-W. Min, “Compact self-interference incoherent digital holographic camera system with real-time operation,” Opt. Express 27(4), 4818–4833 (2019). [CrossRef]  

35. J. Lee, Y. Kim, K. Choi, J. Hahn, S.-W. Min, and H. Kim, “Digital incoherent compressive holography using a geometric phase metalens,” Sensors 21(16), 5624 (2021). [CrossRef]  

36. A. Vijayakumar, T. Katkus, S. Lundgaard, D. Linklater, E. P. Ivanova, S. H. Ng, and S. Juodkazis, “Fresnel incoherent correlation holography with single camera shot,” Opto-Electron. Adv. 3(8), 200004 (2020). [CrossRef]  

37. J. W. Goodman, Introduction to Fourier Optics, 2nd ed. (McGraw-Hill, 1996).

38. J. Zhang, N. Pégard, J. Zhong, H. Adesnik, and L. Waller, “3D computer-generated holography by non-convex optimization,” Optica 4(10), 1306–1313 (2017). [CrossRef]  

39. S. Choi, J. Kim, Y. Peng, and G. Wetzstein, “Optimizing image quality for holographic near-eye displays with michelson holography,” Optica 8(2), 143–146 (2021). [CrossRef]  

40. Z. Ren, J. Zhao, and E. Y. Lam, “Automatic compensation of phase aberrations in digital holographic microscopy based on sparse optimization,” APL Photonics 4(11), 110808 (2019). [CrossRef]  

Supplementary Material (2)

NameDescription
Visualization 1       Optical reconstruction of the 100k high-resolution hologram of frontward focused.
Visualization 2       Optical reconstruction of the 100k high-resolution hologram of backward focused.

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

Fig. 1.
Fig. 1. Schematic diagram of holographic stereogram.
Fig. 2.
Fig. 2. (a) Schematic of the optical configuration of geometric phase based SIDH system. Solid and dashed lines illustrate the optical path of the marginal rays emitted from the rotation center, (b) the flow chart of the proposed method.
Fig. 3.
Fig. 3. Wavefront optimization procedure applied at real and imaginary parts of raw hologram. (a–c) Real part compensation procedure, (d–f) imaginary part compensation procedure.
Fig. 4.
Fig. 4. Wavefront compensation result comparison and digital reconstruction. (a–d) Raw hologram, (e–h) Compensated hologram (i) digital reconstruction of raw hologram, (j) digital reconstruction of compensated hologram.
Fig. 5.
Fig. 5. (a) Concept of the sub-hologram capture system. The experimental setup of (b) emissive target, (c) reflective target.
Fig. 6.
Fig. 6. Numerical simulation result of synthesized hologram. (a) Spatial frequency domain, digital reconstruction results when focusing on the (b–d) backward letter ‘KHU’ and (e-g) forward letter ‘ETRI’, viewing from (b, e) -10°, (c, f) center, (d, g) +10°.
Fig. 7.
Fig. 7. Fabricated 25k holographic stereogram without illumination.
Fig. 8.
Fig. 8. Step-by-step phase-angle part hologram visualization of the proposed method. (a) Raw hologram with resolution of $({{N_x},\; {N_y}} )$, and sampling period of $({\varDelta x,\; \varDelta y} )$, (b) back propagation and resizing result, (c) arbitrary plane carrier wave, (d) plane carrier multiplied result, (e) the result of summation of sub-hologram, (f) binarization of phase-angle part of master-hologram, (g) optical reconstruction without wavefront compensation, (h) backward focused optical reconstruction with wavefront compensation, (i) frontward focused optical reconstruction with wavefront compensation.
Fig. 9.
Fig. 9. Optical reconstruction of the 25k holographic stereogram. (a-c) Backward focus, (d-g) Frontward focus, viewing from (a, d) -10°, (b, f) center, (c, g) 10°.
Fig. 10.
Fig. 10. Experimental result of the 100k holographic stereogram result. (a) Fabricated 100k hologram without illumination. (b–d) Backward focus, (e–g) frontward focus, viewing from (b, e) -25°, (c, f) center, (d, g) 25°, (h–i) close view of (c) and (f) (see Visualization 1 and Visualization 2).

Equations (6)

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

U M a s t e r ( x , y ; z c ) = M N U S u b ( x , y ; z c ) exp ( j ( k x ; m , n x s u b + k y ; m , n y s u b ) ) ,
( k x ; m , n , k y ; m , n ) = 2 π λ ( sin θ x ; m , n , sin θ y ; m , n ) ,
sin Δ θ = λ 2 Δ x = λ f x m a x = N x N x λ Δ x
U c o m ( x , y ) = U o r i U a b e ,
U a b e ( x , y ) = β 1 x 2 + β 2 y 2 + β 3 x y + β 4 x + β 5 y + β 6 ,
L = argmin U a b e | U a b e U | 1 + Λ | U a b e | 1 ,
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.