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

Generation of patterned-phase-only holograms (PPOHs)

Open Access Open Access

Abstract

A fast and non-iterative method for generating a phase-only hologram, hereafter referred to as the patterned-phase-only hologram (PPOH), is reported in this paper. Briefly, a phase mask with a periodic phase pattern is added to the source image, and converted into a hologram. Subsequently, only the phase component is retained as a phase-only hologram. Experimental evaluation reveals that the visual quality of the reconstructed images of the PPOH generated with our proposed method is favorable, and superior to that obtained with existing methods.

© 2017 Optical Society of America

1. Introduction

In the past 2 decades, numerous efforts have been conducted in the generation and processing of digital holograms. Despite the success of these methods, they will not be applicable in practice unless the holograms can be displayed and observed visually. Although a digital hologram can be printed on a physical medium with a fringe printer, the content is static and cannot be changed afterwards. For digital holograms representing animated video signal, they are usually displayed with high-resolution devices such as a Spatial Light Modulator (SLM). However, a hologram is a complex-valued image, but existing SLMs are only capable of presenting either the magnitude, or the phase component of a hologram. A straightforward to address the problem is through the optical integration of a pair of SLMs to the real and imaginary components, or the magnitude and phase components of a hologram [1–3]. It is also possible to convert a complex hologram into a double phase hologram, and display it with 2 phase-only SLM [4]. Despite the effectiveness of these approaches, the optical setup is rather cumbersome and precise alignment between different elements in the display system is needed. The use of a single SLM for displaying a complex hologram, has been studied and reported in [5–8]. In these methods, the orthogonal components of a hologram are displayed in non-overlapping sections in a single SLM. The reconstructed images are combined with a high-resolution grating. On the downside, the usable area of the display is reduced and precise alignment between different elements in the optical setup is required. A more promising solution is to convert the complex hologram into a phase-only hologram (POH), and display it with a phase-only SLM. Notably Buckley and his co-authors have proposed a method, which is referred to as One-Step-Phase-Retrieval (OSPR) [9], that by adding a random phase mask to the source image prior to the generation of the hologram, the latter can be reconstructed with the phase component alone. The method is faster than traditional approach that requires numerous rounds of iterations to derive a phase-only or an amplitude-only hologram for a given target source image (for example [10–12]). On the downside, the lines and edges are fragmented, and the reconstructed images are rather noisy. Degradation on the lines and edges can be alleviated by confining the addition of the phase noise to the smooth region of the source image [13]. However, for the noise contamination, multiple hologram frames, each representing the same source image and which has been added with different random phase noise, have to be generated and presented in rapid succession to reduce the noisy appearance of the reconstructed image. As such, the computation loading on the hologram generation process could be overwhelming. Recently, Tsang et al. have proposed the Bi-directional Error Diffusion (BERD) method [14,15] which enables a complex Fresnel hologram to be converted swiftly into a POH through the Floyd-Steinberg error diffusion [16]. They have shown that a BERD hologram is capable of preserving high fidelity on the object image. Despite the success of this method, error diffusion is a recursive process that is rather time-consuming, and the diffraction efficiency of the reconstructed image is quite low. In another method, known as the sampled-phase-only hologram (SPOH) [17,18], the intensity profile of the object is down-sampled with a uniform lattice prior to the generation of the complex-valued hologram. Subsequently, only the phase component is retained to be the phase-only hologram, and the magnitude is forced to a constant value. This method is capable of preserving good visual quality on the object scene, but resulted in a sparse reconstructed image (i.e., lots of empty holes) that is masked with the texture of the down-sampling lattice. The quality of the reconstructed image can be enhanced by rapidly displaying 2 SPOHs, each deriving from a different down-sampling lattice [19], a method known as the complementary phase-only hologram (CSPOH). The method, though effective, doubles the computational loading and requires a fast SLM to switch the pair of holograms at higher frame-rate. In this paper, we propose a method to overcome the above mentioned problems. Details of our proposed method and experimental evaluation will be provided in Sections 2 and 3, respectively.

2. Proposed method for generating patterned-phase-only hologram (PPOH)

The concept of our proposed method in generating a PPOH of a source image is outlined as follows. To begin with, a phase mask that is comprising of repetitive blocks of phase pattern is added to the source image. The repetitive phase mask effectuates a diffuser action that is capable of smoothing the magnitude of the diffracted waves to a near uniform distribution. However, as identical phase pattern is added to uniformly spaced local regions of the source image, the reconstructed image of the phase hologram will bear a fine pixelated rather than a noisy appearance. In the rest of this section, we shall describe the realization of our proposed method. Consider an intensity source image I(x,y) that is consisting of X columns and Y rows. Our proposed method for generating a PPOH of the source image can be divided into 4 steps and described as follows.

  • Stage 1: The intensity image I(x,y) is partitioned into non-overlapping image blocks of size M×M as shown in Fig. 1.
     figure: Fig. 1

    Fig. 1 Partitioning the intensity image into non-overlapping square image blocks

    Download Full Size | PDF

  • Stage 2: A M×M phase mask P(x,y) is generated with each pixel being assigned a uniformly distributed random phase value θ(x,y) (generated with a random number generator with uniform probability distribution) within the range [0,2π), i.e.,
    P(x,y)|0x,y<M=exp[iθ(x,y)].
  • Stage 3: For each image block, its pixel value is multiplied with the corresponding phase term in the phase mask to give a modified image IP(x,y) given by
    IP(x,y)=I(x,y)×P(m,n),

    where m=xmodM, n=ymodM, and modis the modulus operator. Note that the addition of the phase mask does not change the magnitude of the source image.

  • Step 4: The intensity image is assumed to be parallel to, and at an axial distance zo from the hologram plane. A complex-valued hologram H(u,v) is generated as given by
    H(u,v)=x=0X1y=0Y1Ip(x,y)exp{i2πλ1(xu)2δ2+(yv)2δ2+zo2}.

    Subsequently, the patterned-phase-only hologram HP(u,v) is obtained by retaining the phase component of the hologram:

    HP(u,v)=exp{iarg[H(u,v)]}.

3. Experimental evaluation

To demonstrate our proposed method, we apply it to generate the PPOHs for the source image “Lenna” in Fig. 2(a), based on λ=633nm, δ=8.1um, and M = 4,8, and 12. The source image and the hologram, which are separated at an axial distance of zo=0.03m, have sizes of 512×320 and 1920×1080. The numerical reconstruction of the three PPOHs at the focused plane at 0.03m is shown in Figs. 2(b)-2(d). The peak-signal-to-noise-ratios (PSNR) of the reconstructed images based on M = 4,8, and 12 are 23.97dB, 25.62dB, and 24.7dB, respectively. A larger value of M effectuates a stronger diffuser that smooth out the magnitude to a more uniform distribution, but also cause more distortion to the reconstructed image. Although the visual quality of the 3 reconstructed images are similar, the PSNR is highest for M = 8. As such, we have employed M = 8 in our subsequent evaluation. Next we compare our proposed method with 3 existing non-iterative methods, namely the magnitude removal, noise addition [13], and sampling [17,18] methods, as well as the Gerchberg-Saxton algorithm (GSA). For each method, we generate the phase-only holograms for axial distances ranging from zo=0.02m to 0.1m in step of 0.01m. In the magnitude removal method, only the phase component of the hologram is retained as a phase-only hologram. For the sampling method, the image is down-sampled prior to generating the phase-only hologram, and the sampling factor is set to 8 to provide good visual quality. In the GSA method, 5 rounds of back and forth Fresnel transform [20] are conducted so that the computation time is about 10 times of, and not too far-off from the non-iterative methods. The plots of the PSNRs (based on the source image “Lenna” as the reference image) versus the axial distance of the reconstructed images of the phase-only holograms obtained from each method are shown in Fig. 3. We have observed that the PSNRs of our proposed method is at least 2dB higher than that of the magnitude removal, sampling, and the GSA methods. The PSNR of the noise addition method is rather similar to that of our proposed method. However, as we shall show later, the reconstructed image of a phase-only hologram derived from our proposed method is less noisy than the one obtained with noise addition. Next in Figs. 4(a)-4(e), we display the reconstructed image of the phase-only holograms obtained with different methods. The axial distance is selected to be 0.05m, at which medium performance in terms of PSNR is noted for each method. It can be seen that the visual quality of the reconstructed image corresponding to our proposed method is obvious superior than those of considered existing methods. In comparison with the noise addition method, the reconstructed image of our proposed method has a pixelated appearance, but the noise contamination is less severe.

 figure: Fig. 2

Fig. 2 (a) Source image “Lenna”, (b),(c), and (d) Numerical reconstructed images of the PPOHs of the source image “Lenna” based on M = 4,8, and 12, respectively.

Download Full Size | PDF

 figure: Fig. 3

Fig. 3 PSNRs of reconstructed images of phase-only holograms obtained by the proposed method and 4 other existing methods.

Download Full Size | PDF

 figure: Fig. 4

Fig. 4 (a)-(e) reconstructed images of phase-only holograms obtained by the proposed method, noise addition method, magnitude removal method, sampling method, and the GSA, respectively.

Download Full Size | PDF

To further demonstrate our proposed method, we apply it to generate the phase-only hologram of the image “USAF” in Fig. 5(a). The reconstructed images at 0.03m and 0.06m are shown in Figs. 5(b) and 5(c), respectively. We observe that the reconstructed images are similar to the original one and all the fine bar patterns are clearly resolved. Finally, we would like to show the optical reconstructed image of the PPOHs prepared by our proposed method, representing the “Lenna” and the “USAF” image at axial distance of 0.03m, on a phase-only SLM. The SLM has a size of 1920×1080, and a pixel size of δ=8.1um. The pair of optical reconstructed images are shown in Figs. 6(a) and 6(b). We can see that the quality of the reconstructed images is favorable, and the images resemble the numerical reconstructed images shown in Figs. 4(a) and 5(b).

 figure: Fig. 5

Fig. 5 (a) Source image “UASF, (b) and (c) reconstructed images of the phase-only holograms obtained by the proposed method at 0.03m and 0.06m, respectively.

Download Full Size | PDF

 figure: Fig. 6

Fig. 6 (a),(b) optical reconstructed images of the phase-only holograms obtained by the proposed method for source images “Lenna” and “USAF”, respectively.

Download Full Size | PDF

4. Conclusion

In this paper, we have proposed a non-iterative method for generating a patterned phase-only hologram (PPOH). By illuminating the PPOH with a coherent beam, the source image can be reconstructed with favorable quality, and without the need of optical filters or other optical accessories. In comparison with existing methods, our proposed method has three major advantages. First, the additional computation load (on top of the generation of the hologram based on Fresnel diffraction) is negligible. This is an obvious advantage over methods that require multiple rounds of iterations for deriving the phase-only hologram. Second, as the reconstructed image is obtained with a single PPOH, it is unnecessary, as in the case of the OSPR and the CSPOH methods, to display 2 or more holograms in rapid succession on the SLM. Third, the visual quality of the reconstructed image of the PPOH is higher than that of the magnitude removal, sampling, and GSA methods, and also less noisy than the noise addition method.

References and links

1. M. M. M. Makowski, A. S. A. Siemion, I. D. I. Ducin, K. K. K. Kakarenko, M. S. M. Sypek, A. M. S. A. M. Siemion, J. S. J. Suszek, D. W. D. Wojnowski, Z. J. Z. Jaroszewicz, and A. K. A. Kolodziejczyk, “Complex light modulation for lensless image projection,” Chin. Opt. Lett. 9, 120008 (2011). [CrossRef]  

2. M.-L. Hsieh, M.-L. Chen, and C.-J. Cheng, “Improvement of the complex modulated characteristic of cascaded liquid crystal spatial light modulators by using a novel amplitude compensated technique,” Opt. Eng. 46(7), 070501 (2007). [CrossRef]  

3. R. Tudela, E. Martín-Badosa, I. Labastida, S. Vallmitjana, I. Juvells, and A. Carnicer, “Full complex fresnel holograms displayed on liquid crystal devices,” J. Opt. A5, S189–S194 (2003).

4. C. K. Hsueh and A. A. Sawchuk, “Computer-generated double-phase holograms,” Appl. Opt. 17(24), 3874–3883 (1978). [CrossRef]   [PubMed]  

5. X. Li, Y. Wang, J. Liu, J. Jia, Y. Pan, and J. Xie, “Color holographic display using a phase-only spatial light modulator,” in Dig. Holo. and 3D Img., OSA Technical Digest, OSA (2013), paper DTh2A.3.

6. S. Reichelt, R. Häussler, G. Fütterer, N. Leister, H. Kato, N. Usukura, and Y. Kanbayashi, “Full-range, complex spatial light modulator for real-time holography,” Opt. Lett. 37(11), 1955–1957 (2012). [CrossRef]   [PubMed]  

7. H. Song, G. Sung, S. Choi, K. Won, H. S. Lee, and H. Kim, “Optimal synthesis of double-phase computer generated holograms using a phase-only spatial light modulator with grating filter,” Opt. Express 20(28), 29844–29853 (2012). [CrossRef]   [PubMed]  

8. J.-P. Liu, W. Y. Hsieh, T.-C. Poon, and P. Tsang, “Complex fresnel hologram display using a single SLM,” Appl. Opt. 50(34), H128–H135 (2011). [CrossRef]   [PubMed]  

9. I. Naydenova, Advanced Holography - Metrology and Imaging (InTech, 2011).

10. P. Zhou, Y. Bi, M. Sun, H. Wang, F. Li, and Y. Qi, “Image quality enhancement and computation acceleration of 3D holographic display using a symmetrical 3D GS algorithm,” Appl. Opt. 53(27), G209–G213 (2014). [CrossRef]   [PubMed]  

11. T. Shimobaba, M. Makowski, T. Kakue, M. Oikawa, N. Okada, Y. Endo, R. Hirayama, and T. Ito, “Lensless zoomable holographic projection using scaled Fresnel diffraction,” Opt. Express 21(21), 25285–25290 (2013). [CrossRef]   [PubMed]  

12. Y. Nagahama, T. Shimobaba, T. Kawashima, T. Kakue, and T. Ito, “Holographic multi-projection using the random phase-free method,” Appl. Opt. 55(5), 1118–1123 (2016). [CrossRef]   [PubMed]  

13. P. W. M. Tsang, Y. T. Chow, and T.-C. Poon, “Generation of edge-preserved noise-added phase-only hologram,” Chin. Opt. Lett. 14, 100901 (2016). [CrossRef]  

14. P. W. M. Tsang and T.-C. Poon, “Novel method for converting digital Fresnel hologram to phase-only hologram based on bidirectional error diffusion,” Opt. Express 21(20), 23680–23686 (2013). [CrossRef]   [PubMed]  

15. J.-P. Liu, S. Y. Wang, P. W. M. Tsang, and T.-C. Poon, “Nonlinearity compensation and complex-to-phase conversion of complex incoherent digital holograms for optical reconstruction,” Opt. Express 24(13), 14582–14588 (2016). [CrossRef]   [PubMed]  

16. R. W. Floyd and L. Steinberg, “An adaptive algorithm for spatial grey scale,” Proc Soc. Info. Disp. 17, 75–77 (1976).

17. P. W. M. Tsang, Y. T. Chow, and T.-C. Poon, “Generation of phase-only Fresnel hologram based on down-sampling,” Opt. Express 22(21), 25208–25214 (2014). [CrossRef]   [PubMed]  

18. P. W. M. Tsang, Y. T. Chow, and T.-C. Poon, “Enhancement on the generation of sampled phase-only holograms,” Chin. Opt. Lett. 13, 060901 (2015). [CrossRef]  

19. P. W. M. Tsang, Y. T. Chow, and T.-C. Poon, “Generation of complementary sampled phase-only holograms,” Opt. Express 24(20), 23390–23395 (2016). [CrossRef]   [PubMed]  

20. P. W. M. Tsang, T. C. Poon, and W. Osten, “Review on the state-of-the-art technologies for acquisition and display of digital holograms,” IEEE Trans. Ind’l. Info. 12(3), 886–901 (2016). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1
Fig. 1 Partitioning the intensity image into non-overlapping square image blocks
Fig. 2
Fig. 2 (a) Source image “Lenna”, (b),(c), and (d) Numerical reconstructed images of the PPOHs of the source image “Lenna” based on M = 4,8, and 12, respectively.
Fig. 3
Fig. 3 PSNRs of reconstructed images of phase-only holograms obtained by the proposed method and 4 other existing methods.
Fig. 4
Fig. 4 (a)-(e) reconstructed images of phase-only holograms obtained by the proposed method, noise addition method, magnitude removal method, sampling method, and the GSA, respectively.
Fig. 5
Fig. 5 (a) Source image “UASF, (b) and (c) reconstructed images of the phase-only holograms obtained by the proposed method at 0.03m and 0.06m, respectively.
Fig. 6
Fig. 6 (a),(b) optical reconstructed images of the phase-only holograms obtained by the proposed method for source images “Lenna” and “USAF”, respectively.

Equations (4)

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

P ( x , y ) | 0 x , y < M = exp [ i θ ( x , y ) ] .
I P ( x , y ) = I ( x , y ) × P ( m , n ) ,
H ( u , v ) = x = 0 X 1 y = 0 Y 1 I p ( x , y ) exp { i 2 π λ 1 ( x u ) 2 δ 2 + ( y v ) 2 δ 2 + z o 2 } .
H P ( u , v ) = exp { i arg [ H ( u , v ) ] } .
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.