the whole SLM for one color component. More effective method is to
use multiplexing by encoding three color components into one SLM
with whole area [18]. However, most of these methods only consider
the projection of two dimensional color images. For dynamic holo-
graphic 3D projection, Gaolei Xue proposed a method to reconstruct
full color 3D object by modeling its complex amplitude [19]. In this
method, the CGH of each color component is calculated using polygon-
based methods and three CGHs are synthesized into one CGH by
multiplexing encoding method. However, the calculation time of the
CGH will be proportional to the numbers of the modeled polygons.
In this paper we propose a method to simply calculate CGH for
color holographic 3D projection, based on the combination of the
nonuniform sampled WRP method and the multiplexing encoding
technique. The CGH calculation is simple and fast, the calculation time
is independent on the number of object points. The CGH is capable to
reconstruct zoomable color 3D object in holographic projection system.
2. Method of CGH generation in color holographic 3D
projection
2.1. CGH calculation from 3D object using nonuniform sampled WRP
The nonuniform sampled WRP (NS-WRP) based method [10]
presents the effectiveness to simply and fast calculate a CGH that can
project 3D object with perceptive depth and zoomable function. The
schematic diagram of this method of calculating CGH from a mono-
chrome 3D object is shown in Fig. 1(a). First, a virtual wavefront
recording plane (WRP) which is sufficiently close to the 3D object is
established. The 3D object is represented by a 2D intensity map I(u, v)
and a depth map d(u, v). The intensity map I(u, v) is located at distance
d
min
from the WRP where d
min
is the minimum value of d(u, v). Both of
the intensity map and the WRP is sampled in the same way and their
sampling pitch is given according to the depth map as:
Wuv d duv du(, )= / (, )⋅
min
(1)
Where du is the transversal pixel pitch of the 3D object. In this way,
the different depth value d(u, v) of each object point give rise to the
nonuniform sampled (NS) WRP and the intensity map. Next, the
calculation of the NS-WRP from the NS intensity map is performed by
using the nonuniform fast Fourier transform (NUFFT) based angular
spectrum method expressed as:
xw yw NUFFT NUFFT I u v T f f( , ) = 2{ 1[ ( , )]⋅ ( , )}
xy
(2)
Where W(xw, yw) is the light field of the WRP and T(f
x
, f
y
) is the
transfer function of the angular spectrum method. NUFFT1 and
NUFFT2 denote the first and second type of the NUFFT respectively.
The detailed explanation of NUFFT can be found in Ref. [20–22]. After
that the calculated WRP W(xw, yw) is forcely re-sampled to an
uniform sampled WRP W(xw’, yw’) with the samping pitch of du. In
the last step of the method, the hologram is obtained by calculating the
diffraction from the W(xw’ , yw’) employing the ARSS –Fresnel
diffraction algorithm [23]. The calculation of the hologram H(x, y)is
given by:
x y FFT FFT W xw yw iφ FFT iφ ct( , ) = { [ ( , )⋅exp( )]⋅ [exp( )⋅Re ]}
−1
12
(3)
Where f
1
and f
2
are the quadratic phase term and Rect is a rectangular
function that reduces aliasing noise. Their detailed expression can be
found in Ref. [23].
Fig. 1(b) shows the reconstruction of the hologram calculated by the
above described NS-WRP based method. Due to the resampling
operation from nonuniform sampled WRP to uniform sampled WRP
in the calculation step, each object point of the intensity map will be
reconstructed to a new position at its original depth distance according
to the relations between the reconstruction distance and the sampling
pitch change [10]. Consequently the whole monochrome 3D object
with inherent depth cue can be reconstructed.
2.2. CGH synthesize from color 3D object by multiplexing technique
The calculation process of the CGH from a color 3D object is shown
in Fig. 2. The color 3D object can be divided into RGB component.
Each divided component is represented by a 2D intensity map and
distributed transversely in space as shown in Fig. 2. Then the intensity
map of each color component is propagated toward to the three NS-
WRPs directly and separately, using the NUFFT based angular
spectrum method. After resampling three NS-WRPs to nonuniform
sampled WRPs according to the depth map of the object, we calculate
the sub-hologram of RGB component from each re-sampled WRP.
Here the propagation from each WRP to the hologram is off-axis
diffraction, hence we employ the ARSS- Fresnel diffraction algorithm
[23] to calculate the off-axis diffraction and generate three complex
sub-holograms of RGB component. Finally we superimpose the three
complex sub-holograms of RGB component in order to synthesize one
final CGH.
The reconstruction process of the CGH is shown in Fig. 3. In the
reconstruction, the CGH is illuminated simultaneously by three
reference lights of RGB colors with different angles. By properly
adjusting the angle θ of the three illumination lights, the RGB
component of the target object can be reconstructed at the desired
Fig. 1. (a) Hologram calculation by the NS-WRP based method. (b) Reconstruction of
the hologram..
C. Chang et al.
Optics Communications 387 (2017) 267–274
268