1 February 2009 Three-dimensional light modulation using phase-only spatial light modulator
Author Affiliations +
Optical Engineering, 48(2), 020502 (2009). doi:10.1117/1.3076211
The simultaneous modulation of coherent light distributions on several closely spaced planes by one phase-only spatial light modulator is described to create 3-D light distribution. The pure phase distribution on the phase-only spatial light modulator is retrieved by a new multiplane iteration method based on the Gerchberg-Saxton iteration method, which consecutively propagates the coherent light on several parallel planes with amplitude constraints. A comparison between our proposed method and the complex distribution addition method is presented by numerical simulation. The result indicates that a sharp grayscale image is realized on each closely spaced plane by our proposed method, implying potential applications both in 3-D holographic displays and optical tweezers.
Xia and Yin: Three-dimensional light modulation using phase-only spatial light modulator

The phase distribution of coherent light is difficult to measure directly, because phase information is lost in the measured intensity. Phase retrieval is used to determine the phase distribution on one plane by the intensity distribution measured from several other planes, and thereby is widely used both for wavefront sensing and wavefront modulation. For wavefront sensing, as in x-ray crystallography, electron microscopy,1, 2 and astronomy,3 the intensity distributions on several output planes are directly measured, while the coherent light distribution on the input plane is “sensed” by a retrieval algorithm. For wavefront modulation, as in optical tweezers4 and 2-D video frame projection,5, 6 the coherent light on the input plane is directly “modulated” to generate a special image pattern on the output plane.

A phase-only spatial light modulator (POSLM) is widely used for wavefront modulation. Since only the phase distribution is modulated, leaving the amplitude intact, no energy is lost and thus very high diffraction efficiency is achieved. Recently, liquid crystal on silicon (LCOS) with 2π range of phase modulation was developed to be used as POSLM.7 By modulating the “pure” phase on LCOS, the amplitude on the output plane with arbitrary distribution was realized, even at high video frame rates.6

To generate a special image pattern on the output plane, the required phase distribution on the LCOS modulator is deduced by either iteration algorithms8, 9 or direct algorithms.10 These algorithms are based on scalar diffraction theory,11 in which the coherent light propagates from one plane to another, and therefore they are efficient for deducing the light distribution in a 2-D plane. Modulation of coherent light in 3-D space will extend the application field of phase retrieval, for example, 3-D holographic displays12 and 3-D manipulation of biomedical particles.4 3-D light distribution can be realized by simultaneously modulating coherent light in several parallel 2-D planes. Two kinds of multiplane iteration methods have been introduced to modulate 3-D light distribution. One is based on consecutive propagation between multiplanes;13, 14 the other is based on the addition of complex distribution directly propagated from multiplanes.15, 16 However, to converge a sharp image, the distance between planes cannot be very small in the first category method, and the reconstructed image in the second category method sometimes looks dirty due to the addition of complex distribution propagated from several planes. The purpose of this study is to develop a multiplane iteration method that can modulate the intensity distribution on several closely spaced planes with sharp images.

As mentioned before, a 3-D light distribution can be discretized into several 2-D planes. The light distribution on each plane is efficiently calculated by fast Fourier transform (FFT) based on the scalar diffraction theory. Our iteration method is based on Gerchberg-Saxton’s iteration method.8, 9 To retrieve the pure phase distribution on a POSLM, forward and backward propagation is calculated from plane to plane, and the amplitude distribution on each output plane is modified according to the amplitude constraint, which is the image pattern that will be modulated on each plane.

Figure 1 shows the proposed iteration method to modulate images on two output planes by a POSLM. The choice of two output planes is arbitrary and does not affect the generality of the method and results. The method consists of four steps: 1. the light starts at the POSLM and propagates to the first plane (P1) ; 2. the amplitude distribution on plane P1 is modified according to the amplitude constraint and then propagates to plane P2 ; 3. the amplitude distribution on the final plane P2 is inversely propagated to plane P1 , after applying the amplitude constraint; 4. the light distribution on plane P1 is inversely propagated to the POSLM, and finally the amplitude distribution on POSLM is modified to a uniform light distribution. During each step, the phase distribution on each plane is calculated by scalar diffraction from the previous plane, combined with the modified amplitude distribution, and is then propagated to the next plane.

Fig. 1

Scheme of proposed multiplane iteration algorithm.


The angular spectrum method,11 based on the scalar diffraction theory, is used to propagate light between output planes, because it is suitable for light propagation over a short distance, and thus high resolution 3-D light distributions can be realized by discretizing in small steps. Moreover, the sampling space is invariant between two propagation planes. The distance between the POSLM and the first plane (P1) is not limited to very short distances. For propagation over a long distance, single FFT Fresnel diffraction is used. It should be noted that the sampling space during Fresnel diffraction is changed, and thus forward and inverse propagation should take place between the same two planes. Therefore, during inverse propagation, the light is first propagated back to plane P1 in step 3 and then finally propagated back to the POSLM plane in step 4.

The numerical simulation is based on the scheme depicted in Fig. 1. The wavelength used in the simulation is λ=5×107m . The distance between the POSLM and the first plane (P1) is L1=0.5m . The distances between the two output planes is L2=1×102m . The sampling space on the output plane is dx=dy=2×105m , and the sampling numbers in horizontal and vertical directions are both N=512 . The sampling space on the POSLM plane is determined by Fresnel diffraction.

Figures 2 and 2 show the amplitude constraints that are used in the simulation on two output planes, i.e., planes P1 and P2 , respectively. They are generated by segmenting the original image in Fig. 2. The resolution of these three images is 256×256 , and is padded to 512×512 by zeros during the calculations to avoid circular convolution. Figure 2 shows the constraint area on plane P2 . When light propagates to this plane, the amplitude distribution inside the constraint area (the white area) is substituted by the image content, and the amplitude distribution in the dark area is intact during propagation. The constraint area on plane P1 is an inverse one that is the dark area in Fig. 2.

Fig. 2

(a) The original image used to generate constraint. (b) and (c) are amplitude constraints on planes P1 and P2 , respectively. (d) The constraint area on plane P2 .


Figures 3 and 3 show the simulated amplitude distribution on planes P1 and P2 after 512 iterations based on our proposed method. Inside the constraint area on each plane, the amplitude distribution converges to a sharp image. At the same area on the adjacent planes, this sharp image becomes a smeared image due to diffraction. For comparison, Figs. 3 and 3 show the simulated amplitude distribution based on the method in Ref. 16. The edges look blurred due to the light distribution being close to the edge on one plane propagated into the other side on other plane. The convergence of two algorithms presented in Fig. 3 is also different. Our proposed algorithm shows a fast convergence due to the do-not-care area outside the constraint.6 The root mean square error (rms) of the calculated amplitude distribution on plane P2 is lower than 10 after 4 times of iteration. In contrast, the rms for the method in Ref. 16 is higher than 35 after 512 times iteration on the same plane due to the blurred edges.

Fig. 3

The simulated amplitude distributions on planes P1 and P2 , respectively. (a) and (b) are the results by our method. (c) and (d) are the results by the complex distribution addition method.


As mentioned before, the output planes are placed close to each other in our simulations; this enables the generation of high resolution 3-D light distribution. However, they also display smeared images on the other planes. Fortunately, this is not an issue in some applications. For instance, when the 3-D light distribution is discretized into several 2-D planes and placed close to human eyes,12 people will clearly perceive a sharp image on one plane, and the images on other planes look blurred due to accommodation. Another example is a multiplane projection: when we use a POSLM to project an image on an arbitrarily shaped surface, e.g., a sphere, this 3-D surface can be discretized into several 2-D planes. On each plane a sharp image is projected on one area and stops further propagation.

In conclusion, we develop a new iteration method to project a 3-D image on several output planes by one POSLM. The numerical simulation implies that this method is very suitable for modulation of 3-D images with very fine structure. The method described can be used to generate 3-D holographic displays, and also has potential application in optical tweezers.


This work is supported by the national high technology research and development program of China under the grant number 2007AA01Z303. The authors wish to thank Daniel den Engelsen for fruitful comments and help to prepare this work.


1.  J. S. Wu, U. Weierstall, J. C. H. Spence, and C. T. Koch, “Iterative phase retrieval without support,” Opt. Lett.0146-9592 10.1364/OL.29.002737 29, 2737–2739 (2004). Google Scholar

2.  V. Elser, “Phase retrieval by iterated projections,” J. Opt. Soc. Am. A Opt. Image Sci. Vis1084-7529 20, 40–55 (2003). Google Scholar

3.  C. Roddier and F. Roddier, “Combined approach to Hubble space telescope wave-front distortion analysis,” Appl. Opt.0003-6935 32, 2992–3008 (1993). Google Scholar

4.  D. G. Grier, “A revolution in optical manipulation,” Nature (London)0028-0836 10.1038/nature01935 424, 810–816 (2003). Google Scholar

5.  M. Makowski, M. Sypek, and A. Kolodziejczyk, “Colorful reconstructions from a thin multi-plane phase hologram,” Opt. Express1094-4087 16, 11618–11623 (2008). Google Scholar

6.  A. Georgiou, J. Christmas, N. Collings, J. Moore, and W. A. Crossland, “Aspects of hologram calculation for video frames,” J. Opt. A, Pure Appl. Opt.1464-4258 10, 035302 (2008). Google Scholar

7.  H. T. Dai, K. S. Xu, Y. J. Liu, X. Wang, and J. Liu, “Characteristics of LCoS Phase-only spatial light modulator and its applications,” Opt. Commun.0030-4018 238, 269–276 (2004). Google Scholar

8.  R. W. Gerchberg and W. O. Saxton, “A practical algorithm for the determination of phase from image and diffraction plane pictures,” Optik (Jena)0030-4026 35, 237–246 (1972). Google Scholar

9.  J. R. Fienup, “Phase retrieval algorithms: a comparison,” Appl. Opt.0003-6935 21, 2758–2769 (1982). Google Scholar

10.  T. E. Gureyev, A. Pogany, D. M. Paganin, and S. W. Wilkins, “Linear algorithms for phase retrieval in the Fresnel region,” Opt. Commun.0030-4018 10.1016/j.optcom.2003.12.020 231, 53–70 (2004). Google Scholar

11.  J. W. Goodman, Introduction to Fourier Optics, 3rd ed., Chap. 3, McGraw-Hill, New York (2005). Google Scholar

12.  A. Schwerdtner, N. Leister, and R. Häussler, “A new approach to electro-holography for TV and projection display,” SID Intl. Symp. Digest Tech. Papers, Vol. 38, pp. 1224–1227, Academic, New York (2007). Google Scholar

13.  M. Makowski, M. Sypek, A. Kolodziejczyk, and G. Mikula, “Three-plane phase-only computer hologram generated with iterative Fresnel algorithm,” Opt. Eng.0091-3286 10.1117/1.2148980 44, 125805 (2005). Google Scholar

14.  M. Makowski, M. Sypek, A. Kolodziejczyk, G. Mikula, and J. Suszek, “Iterative design of multiplane holograms: experiments and applications,” Opt. Eng.0091-3286 10.1117/1.2727379 46, 045802 (2007). Google Scholar

15.  T. Haist, M. Schonleber, and H. J. Tiziani, “Computer-generated holograms from 3D-objects written on twisted-nematic liquid crystal display,” Opt. Commun.0030-4018 10.1016/S0030-4018(97)00192-2 140, 299–308 (1997). Google Scholar

16.  G. Sinclair, J. Leach, P. Jordan, G. Gibson, E. Yao, Z. J. Laczik, M. J. Padgett, and J. Courtial, “Interactive application in holographic optical tweezers of a multi-plane Gerchberg-Saxton algorithm for three-dimensional light shaping,” Opt. Express1094-4087 10.1364/OPEX.12.001665 12, 1665–1670 (2004). Google Scholar

© (2009) Society of Photo-Optical Instrumentation Engineers (SPIE)
Jun Xia, Hanchun Yin, "Three-dimensional light modulation using phase-only spatial light modulator," Optical Engineering 48(2), 020502 (1 February 2009). http://dx.doi.org/10.1117/1.3076211


Phase shift keying

3D displays


3D image processing


Algorithm development

Back to Top