Computer-generated holograms for three-dimensionalsurface objects with shade and texture
Kyoji Matsushima

Digitally synthetic holograms of surface model objects are investigated for reconstructing threedimensional objects with shade and texture. The objects in the proposed techniques are composed ofplanar surfaces, and a property function defined for each surface provides shape and texture. The fieldemitted from each surface is independently calculated by a method based on rotational transformationof the property function by use of a fast Fourier transform (FFT) and totaled on the hologram. This
technique has led to a reduction in computational cost: FFT operation is required only once for calculating
a surface. In addition, another technique based on a theoretical model of the brightness of the reconstructed surfaces enables us to shade the surface of a reconstructed object as designed. Optical reconstructions of holograms synthesized by the proposed techniques are demonstrated. © 2005 OpticalSociety of America
OCIS codes: 090.1760, 090.2870, 090.1970.

1. Introduction

Computer-generated
hologramsforthreedimensional (3-D) displays, sometimes called digitally synthetic holograms, are desired media forcreating 3-D autostereoscopic images of virtual objects. However, the technology suffers from two problems: the necessity for extremely high spatial
resolution to fabricate or display the holograms, andlong computation times for the creation, especially infull parallax holograms.
For the past decade, techniques using point sourcesof light have been widely used to calculate object
waves.1,2 This point source method is simple in principle and potentially the most flexible for synthesizing holograms of 3-D objects. However, because it is
too time consuming to create full parallax holograms,3 many methods to reduce the computation
time, including geometric symmetry,4 look-up tables,3difference formulas,5 recurrence formulas,6 employing
computer-graphics hardware,7 and constructing special CPUs,8 have been attempted.
Point source methods for calculating spherical

The author (matsu@kansai-u.ac.jp) is with the Department ofElectrical Engineering and Computer Science, Kansai University,
3-3-35 Yamate-cho, Suita, Osaka 564-8680, Japan.
Received 17 June 2004; revised manuscript received 8 March
2005; accepted 21 March 2005.
0003-6935/05/224607-08$15.00/0© 2005 Optical Society of America

waves emitted from point sources are commonly ray
oriented. As they trace the ray from a point source toa sampling point on the hologram, the procedure issometimes referred to as ray tracing.2 However, there
are also wave-oriented methods to calculate object
fields in which fields emitted from objects defined asplanar segments9,10 or 3-D distributions of fieldstrength11 are calculated by methods based on wave
optics. The major advantage of wave-oriented methods is that they can use a fast Fourier transform(FFT) for numerical calculations. Therefore the computation time is shorter than for point source methods, especially in full parallax holograms. However,
the optical reconstruction of accurately rendered 3-D
objects such as a shaded cube, as reported for waveoriented methods, was not discussed in the papers
cited above. This is so because of a lack of well-definedprocedures to generate object fields for arbitrarily
shaped surfaces that are diffusive and sometimes
have texture. The technique for shading the reconstructed object according to such design parameters
as the position of the illumination light and the ratioof the surrounding light is also important in creating
real 3-D images by wave-oriented methods.
In wave-oriented methods, calculating fields arecommonly based on coordinate transformation in
Fourier space.10,11 A similar method based on theRayleigh–Sommerfeld integral has been reported
within the context of free-space beam propagation.12
Recently, the author reported a more precise formulation and numerical consideration13 as an extension
1 August 2005 兾 Vol. 44, No. 22 兾 APPLIED OPTICS

4607

Fig. 2. Fields emitted from surfaces with (a) a constant phase, (b)a diffusive phase, and (c) a diffusive phase multiplied by the phaseof a plane wave propagating to a hologram.

Fig. 1. Geometry and definitions of global coordinates and tilted
local coordinates defined for a planar surface.

of an angular spectrum of plane waves14 in which
remapping the angular spectrum plays an importantrole. The remapping also eases the difficulty of creating object fields in wave-oriented methods.
In this paper two techniques for synthesizing object
fields in surface models are presented for creating
3-D images by use of computer-generated holograms.
The first technique, based on the rotational transformation of wave fields presented in Ref. 13 and on
remapping of the angular spectrum, provides amethod for synthesis of the object fields. This technique makes it possible to create diffusive fields ofarbitrarily tilted planar surfaces that have an arbitrary shape and texture. Furthermore, another technique is also presented for avoiding unexpectedchanges in brightness of the surfaces of objects. The
technique enables us to render surface objects as thedesigners intended.
2. Object Model and the Property Function of Surfaces

The coordinate systems and geometry used in thisstudy are shown in Fig. 1. Objects consist of planar
surfaces that are diffusive and luminous by reflecting
virtual illumination. Each surface has its own two
local coordinates, called tilt and parallel. The tilted
local coordinates defined for the nth planar surface
are denoted rn ⫽ 共xn, yn, zn兲, defined such that theplanar surface is laid on the 共xn, yn, 0兲 plane. A complex function hn共xn, yn兲 is defined on the plane to givethe nth surface such properties as shape, brightness,diffusiveness, and texture. Thus these complex functions are referred to as the property functions of thesurface.
Parallel local coordinates r̂n ⫽ 共x̂n, ŷn, ẑn兲 are also
defined for each surface. They share their origin withthe tilted coordinates, but the axes are parallel tothose of the global coordinates. In the global coordinates, denoted r ⫽ 共x̂, ŷ, ẑ兲, the hologram is placed on
the 共x̂, ŷ, 0兲 plane. All property functions of surfaces
are defined in the following form:hn(xn, yn) ⫽ an(xn, yn)⌿(xn, yn)pn(xn, yn),
4608

APPLIED OPTICS 兾 Vol. 44, No. 22 兾 1 August 2005

(1)

where an共xn, yn兲 is a real function that provides amplitudes of the property function to keep the shapeand the texture of the nth surface.
If the property function is defined only as the amplitude distribution, the surface yields little diffusiveness, as shown in Fig. 2(a). For example, an共xn, yn兲
in Fig. 1 is a simple rectangular function; i.e., theamplitude is constant within the rectangular surface.
This situation is similar to the optical diffraction of aplane wave by a rectangular aperture; therefore, if
the surface is visible to the naked eye, the light has
not been sufficiently diffracted by the aperture. To
give surfaces large diffusiveness, the amplitude functions must be multiplied by a given diffusive phase:⌿(xn, yn) ⫽ exp[ik␾d(xn, yn)],(2)where ␾d共xn, yn兲 is a phase that behaves as a numerical diffuser. Random functions are candidates for the
diffusive phase, but full random functions are notappropriate to the diffusive phase because the random phases are discontinuous and have a large Fourier frequency. Thus the random phases cause
speckles in the reconstruction and problems in numerical calculation. In the research reported in thispaper, a digital diffuser proposed for Fourier holograms15 is used for phase function ␾d共xn, yn兲.
If a property function is given by the product of theamplitude function and the diffusive phase, the carrier frequency of the field on the tilted 共xn, yn, 0兲 planeis zero. This forces the surface to emit light perpendicularly to the surface, as shown in Fig. 2(b). If thesurface is sufficiently diffusive, a portion of the emitted field may reach the hologram, but high diffusiveness results in high computational costs such as theneed for a great number of sampling points. Therefore the phase of a plane wave propagating perpendicularly to the hologram should be multiplied by thetwo factors given above. This plane-wave factor
causes the field to propagate into the hologram, expressed bypn(xn, yn) ⫽ exp[i(kx, nxn ⫹ ky, nyn)],(3)

where kx, n and ky, n are the xn and yn components,
respectively, of the wave vector of the plane wave.
The property function given by hn共xn, yn兲 is transformed into the complex amplitude ĥn共x̂n, ŷn兲 in theparallel coordinates by the method described in Sec-tion 5 below. When this transformation is written asĥn(x̂n, ŷn) ⫽ ᏾␪x␪y␪z{hn(xn, yn)},Ĥ(û, v̂) ⫽ H[␣(û, v̂), ␤(û, v̂)].

(4)

fields from all surfaces are superimposed upon thehologram plane as follows:ĥ(x̂, ŷ) ⫽ 兺 ᏼdn{ĥn(x̂n, ŷn)},given byComplex amplitudes of the field are obtained inparallel coordinates by inverse Fourier transformation of the spectrum in the paraxial condition13:(5)ĥ(x̂, ŷ) ⫽ Ᏺ⫺1{Ĥ(û, v̂)}nwhere ␪x, ␪y, and ␪z are rotation angles on each axisand ᏼdn兵 其 represents translational propagation
through distance dn.
3. Rotational Transformation of Property Functions

The details of rotational transformation were alreadyreported in Ref. 13. However, I summarize it for convenience in this section, because rotational transformation of wave fields is the core of techniques
proposed in this paper. Note that rotation of just asingle planar surface is presented in this section;
therefore the subscript n is omitted in this section.
A.

Rotation of Coordinates in Fourier SpaceThe Fourier spectrum of a property function is givenin the tilted local coordinates asH(u, v) ⫽ Ᏺ{h(x, y)}
⫽

冕冕

⬁

⫽

(6)
where u and v denote the Fourier frequencies in thex and y axes, respectively. The frequency in the z axis
is not independent of u and v and is given by
w共u, v兲 ⫽ 共␭⫺2 ⫺ u2 ⫺ v2兲1兾2, where ␭ is the wavelength. The relation among Fourier frequencies in theparallel local coordinates 关û, v̂, ŵ共û, v̂兲兴 is analogous
to this and given by ŵ共û, v̂兲 ⫽ 共␭⫺2 ⫺ û2 ⫺ v̂2兲1兾2. In
Fourier space, the frequencies 关u, v, w共u, v兲兴 can betransformed into 关û, v̂, ŵ共û, v̂兲兴 by ordinary
coordinate-transformation procedures and vice versa.
Suppose that r is transformed into r̂ by a rotation
matrix T, i.e., r̂ ⫽ Tr and r ⫽ T⫺1r̂. The frequencies
in the parallel coordinates are given asu ⫽ ␣(û, v̂) ⫽ a1û ⫹ a2v̂ ⫹ a3ŵ(û, v̂),v ⫽ ␤(û, v̂) ⫽ a4û ⫹ a5v̂⫹a6ŵ(û, v̂),(7)

where the rotation matrix is defined as冋册

a1 a2 a3
⫺1
T ⫽ a4 a5 a6 .
a7 a8 a9冕冕⬁Ĥ(û, v̂)exp[i2␲(x̂û ⫹ ŷv̂)]dûdv̂.

⫺⬁

(10)
B.

Remapping the Fourier Spectrum

In the actual numerical calculation, to use a FFT one
must sample the spectrum as well as the field at an
equidistant sampling grid within a given samplingarea. The coordinate transformation of Eq. (9), however, causes distortion and a shift of the samplingarea and points. This distortion and shift are dependent on the way the tilted local coordinates are defined.16
Figure 3(a) shows an example of the definition inwhich the parallel coordinates r̂ are transformed into
the tilted coordinates r by rotation of the coordinates
on the ẑ axis before the y axis. In this case, rotation
matrix T⫺1 is given byh(x, y)exp[⫺i2␲(ux ⫹ vy)]dxdy,

⫺⬁

T⫺1 ⫽

冋

Therefore the spectrum in the parallel coordinates is册cos ␪y cos ␪z cos ␪y cos ␪z ⫺sin ␪y⫺sin ␪zcos ␪z0.
sin ␪y cos ␪z sin ␪y sin ␪z cos ␪y(11)

Figure 3(b) shows the sampling area of Ĥ共û, v̂兲 forseveral rotation angles when h共x, y兲 is sampled atintervals of ␦x ⫽ ␦y ⫽ 2 ␮m in ␭ ⫽ 633 nm. The
sampling area of H共u, v兲, a ␦x⫺1 ⫻ ␦y⫺1 square, isdistorted and shifted. Therefore, resampling accompanied with interpolation is necessary for using an
inverse FFT in calculating the complex amplitudes ofthe field. However, simple interpolation is not sufficient because the FFT algorithm does not work effectively for spectra sampled far from zero frequency.
This shift can be interpreted to be a carrier frequency
observed in the parallel local coordinates.
Let us reverse the procedure of rotation of the coordinates and shift the origin of Fourier space 共u, v兲
in the tilted coordinates. The origin of Fourier space
共û, v̂兲 in the parallel coordinates is inversely projectedto the frequencies 共u0, v0兲 in the tilted coordinates bymatrix (8) as follows:u0 ⫽ ␣(0, 0) ⫽ a3/␭,

(8)

(9)

v0 ⫽ ␤(0, 0) ⫽ a6/␭.

(12)

To ensure that the center of the spectrum in theparallel coordinates is located at the origin of theFourier space after rotation of the coordinates, thefollowing new shifted Fourier space should be intro1 August 2005 兾 Vol. 44, No. 22 兾 APPLIED OPTICS

4609

The exponential factor in brackets in Eq. (16) is attributed to the carrier frequency observed in the parallel coordinates, whereas factor p共x, y兲 of theproperty function was introduced to force the emitted
field toward the hologram, canceling the carrier frequency in the parallel coordinates. In fact, the exponential factors of Eq. (16) and p共x, y兲 cancel eachother out. Equation (16) is rewritten by substitutionof Eq. (3) as follows:H⬘(u⬘, v⬘) ⫽ Ᏺ{a(x, y)⌿(x, y)exp{i[(kx ⫺ 2␲u0)x⫹ (ky ⫺ 2␲v0)y]}},(17)
where the subscript n is omitted again. The wave
vector of a plane wave propagating along the ẑ axis isexpressed by 共0, 0, 2␲兾␭兲 in parallel coordinates.
Thus the plane wave in the tilted coordinates is obtained by coordinates rotation by use of matrix (8) as
follows:
kx ⫽ 2␲a3兾␭,ky ⫽ 2␲a6兾␭.

(18)

The spectrum of relation (15) is rewritten by substitution of Eqs. (18) and (12):H⬘(u⬘, v⬘) ⫽ Ᏺ{a(x, y)⌿(x, y)}.

(19)

As a result, the factor p共x, y兲 is no longer required in
the property function if the spectrum is calculated in
shifted Fourier space 共u⬘, v⬘兲. Therefore let us redefine the property function as
Fig. 3. Schematic of rotation upon two axes: (a) a plane rotatedupon the ẑ axis before the x axis and (b) resampling areas of theFourier spectrum at several rotation angles in the rotation scheme.

h(x, y) ⬅ a(x, y)⌿(x, y),(20)

and its spectrum asH(u, v) ⬅ Ᏺ{h(x, y)}.

duced in the tilted coordinates:u⬘ ⫽ u ⫺ u0,v⬘ ⫽ v ⫺ v0.

(13)

The spectrum expressed in shifted Fourier space
共u⬘, v⬘兲 is written as
H⬘(u⬘, v⬘) ⫽ H(u⬘ ⫹ u0, v⬘ ⫹ v0).

(14)

The spectrum in the parallel coordinates is obtainedby remapping spectrum H⬘共u⬘, v⬘兲 onto Fourier space
共û, v̂兲 as follows:Ĥ(û, v̂) ⬵ H⬘(u ⫺ u0, v ⫺ v0)⫽ H⬘(␣(û, v̂) ⫺ u0, ␤(û, v̂) ⫺ v0),(15)

where the sign for nearly equal means that an interpolation is required.
The Fourier spectrum in the shifted Fourier space
is obtained by application of the shift theorem of theFourier-transform theory to Eq. (14):H⬘(u⬘, v⬘) ⫽ Ᏺ{h(x, y)exp[⫺i2␲(u0x ⫹ v0y)]}. (16)4610

APPLIED OPTICS 兾 Vol. 44, No. 22 兾 1 August 2005

(21)

Consequently, the rotational transformation issummarized as follows: First, one obtains spectrum
H共u, v兲 of the property function of Eq. (20) by fast
Fourier transformation. The center of the spectrum isplaced at the origin in the Fourier space. Next, thespectrum in the parallel coordinates is obtained byremapping spectrum H共u, v兲, expressed by substituting Eq. (12) into relation (15) as follows:Ĥ(û, v̂) ⬵ H(␣(û, v̂) ⫺ ␣(0, 0), ␤(û, v̂) ⫺ ␤(0, 0)).
(22)
Finally, the complex amplitudes of the field are obtained in the parallel coordinates by an inverse Fourier transformation of Eq. (10).
4. Holograms of a Single Surface with Texture
A.

Single Axis Rotation

The hologram of a single planar surface with texture
was fabricated for verifying the technique describedin Section 5. The planar surface and the hologram areFig. 4. (a) Planar object used for fabricating the hologram of aplane rotated upon a single axis. (b) Geometry for capturing thereconstruction. The dimensions of texture of a checker embeddedin the property function are 16.4 mm ⫻ 8.2 mm.

sampled at intervals of 2 ␮m in the x axis and 4 ␮min the y axis. The planar surface has sampling pointsof 16,384 ⫻ 4096 and a binary texture. Amplitude
distribution a共x, y兲 is shown in Fig. 4(a). The surface
object was placed at ẑ ⫽ ⫺10 cm and rotated only on
the y axis at an angle of 80°, as shown in Fig. 4(b). The
hologram with 8192 ⫻ 4096 pixels was encoded by apoint-oriented method and fabricated by a special
printer constructed for printing synthetic fringes.17
The hologram is reconstructed by a 633 nm He–Ne
laser, and the reconstructed image is captured by adigital camera placed at ẑ ⫽ 10 cm. Figure 5 shows
photographs of a reconstructed image. The camera,
whose focus is fixed, was moved from left to right for
Figs. 5(a)–5(c) and back and forth for Figs. 5(d)–5(f).
The apparent size of the texture pattern changes
when the viewpoint is moved in Figs. 5(a)–5(c), while
the defocused position of the texture pattern changes
when the camera is moved because the focal plane of
the camera moves in Figs. 5(d)–5(f).
B.

Two-Axis Rotation

Fig. 6. Planar object used for fabricating a hologram in two-axis
rotation.

The number of sampling points and sampling pitchesof the hologram is the same as in the single-axis
rotation. Optical reconstructions of the holograms areshown in Fig. 7. Four holograms with different rotation angles were fabricated and reconstructed. The
appearance of texture on the planar surfaces varies
according to the rotation angle.
5. Holograms of a Three-Dimensional Object and
Its Shading

Three-dimensional objects such as cubes and pyramids can be built from planar surfaces. Thus, one cansynthesize the fields of these objects by superimposing the fields from the planar surfaces onto the hologram. However, a problem that does not exist when
one is synthesizing a single-surface object is that thebrightness of the reconstructed surface varies, depending on the angle of the surfaces. As a result,objects are shaded as if an unexpected illumination
were throwing light. In a single-surface object, theHolograms of a single planar surface rotated on twoaxes were also fabricated and optically reconstructed.
The method of rotation of two axes is the same asshown schematically in Fig. 3. The amplitude image
with 8192 ⫻ 4096 sampling points is shown in Fig. 6.

Fig. 5. Optically reconstructed images of a hologram captured bymoving a camera (a)–(c) from left to right and (d)–(f) back andforth.

Fig. 7. Optical reconstructions of holograms of planar surfaces
rotated at several angles.
1 August 2005 兾 Vol. 44, No. 22 兾 APPLIED OPTICS

4611

Fig. 8. Model of brightness of a planar surface expressed by aproperty function sampled at an equidistant grid.

change of brightness of a surface is not perceived
because there is nothing to compare with the single
surface in a piece of hologram. This unexpected andunwanted change of brightness must be resolved if
one is to shade the object as intended.
A.

Theory of Brightness of Reconstructed Surfaces

To compensate for unexpected shading it is necessary
to investigate which parameters govern the brightness of the surface in reconstruction. Figure 8 is atheoretical model that predicts the brightness of asurface represented by sampled property function
h共x, y兲. Suppose that the amplitude of a property
function of a surface is a constant, i.e., that a共x, y兲⬅ a, and suppose that a2 provides optical intensity on
the surface. In such cases, the radiant flux ⌽ of asmall area ␦A on the surface is given by⌽⫽冕冕|h(x, y)|2dxdy

␦A⬵ ␦A␴a2,(23)

where ␴ is the surface density of the sampling points.
Assuming that the small area emits light within adiffusion angle in a direction at ␪v to the normal
vector, the solid angle corresponding to the diffusion
cone is given as ⍀ ⫽ A兾r2, where A ⫽ ␲共r tan ␺d兲2 isthe section of the diffusion cone at a distance r and ␺dis the diffusion angle of light that depends on diffuser
function ⌿共x, y兲 of Eq. (1).
According to photometry, brightness of the surface,observed in a direction at an angle ␪v, is given by
L⫽

d⌽兾d⍀
.
cos ␪v␦A(24)

Assuming that light is diffused almost uniformly, i.e.,that d⌽兾dA ⯝ ⌽兾A, the brightness is rewritten by
substitution of d⌽ ⯝ 共⌽兾A兲dA, d⍀ ⫽ dA兾r2, and relation (23) into Eq. (24) as follows:
L⯝
4612

␴a2
␲ tan2 ␺d cos ␪v

.

(25)

APPLIED OPTICS 兾 Vol. 44, No. 22 兾 1 August 2005

Fig. 9. Curves of the angle factor for several values of ␥.

As a result, the brightness of the surface depends
on the surface density of sampling, the diffusivenessof the diffuser function, and the amplitude of thesurface property function. In addition, the brightnessof the surface is governed by observation angle ␪v. In
other words, if several surfaces with the same property function are reconstructed from a hologram, thebrightness varies according to the direction of thenormal vector of the surface. This phenomenon
causes unexpected shading.
Inasmuch as only a simple theoretical model has
been discussed so far, relation (25) is only partiallyappropriate for measuring the brightness of opticallyreconstructed surfaces of real holograms. The brightness given in relation (25) diverges in the limit
␪v → ␲兾2, but an actual hologram cannot produce infinite brightness for its reconstructed surface. Thus
relation (25) is not sufficient to compensate for thebrightness. To avoid the divergence of brightness
in relation (25), one should introduce angle factor
共1 ⫹ ␥兲兾共cos ␪v ⫹ ␥兲 shown in Fig. 9, instead of1兾cos ␪v, a priori. This angle factor is unity in ␪v⫽ 0 and 1 ⫹ 1兾␥ in ␪v ⫽ ␲兾2. Consequently, thebrightness is given as an expression ofL⫽␴a2(1 ⫹ ␥),(cos
␪v ⫹ ␥)␲ tan ␸d2

(26)

where ␥ is a parameter that plays a role in preventing
the divergence of brightness and in preventing overcompensation. Because ␥ is dependent on actual
methods for fabricating holograms, such as encoding
the field or the property of recording materials, it
should be determined experimentally.
B.

Compensation for Brightness and Shading Objects

The amplitude of a property function that reconstructs a surface in a given brightness L is obtainedby solution of Eq. (26) for a as follows:冋册L␲ tan2 ␸d (cos ␪v ⫹ ␥) 1兾2a⫽.
␴
(1 ⫹ ␥)(27)

Fig. 10. Optical reconstructions of unshaded hexagonal prisms (a) without brightness compensation and (b), (c) with compensation in␥ ⫽ 0, 0.5, respectively.

However, angle ␪v is unknown in synthesizing theobject field, and therefore it seems impossible to compensate for the change of brightness. But holograms
are observed in a direction along the ẑ axis, i.e., perpendicular to the hologram, because the hologram isusually observed at a distance of more than several
tens of centimeters. Hence it is possible to approximate ␪v by an angle ␪n formed between the nth surface and the hologram. Objects are shaded by amethod based on Lambert’s law and the diffused reflection model. The brightness of the nth surface, of
which the normal vector forms angle ^␪n with the vector of illumination, is given by
Ln ⫽ L0(cos ^␪n ⫹ le),(28)

where le is the ratio of the surrounding light to theillumination and L0 is brightness in ^
␪n ⫽ 0 and le⫽ 0. By substitution of Ln of Eq. (28) into L of Eq. (27)amplitude an of the nth surface is given as follows:冋(cos ^
␪n ⫹ le)(cos ␪n ⫹ ␥)an ⫽ a0
1⫹␥

a0 ⬅

冋

册

L0␲ tan2 ␸d 1兾2
.
␴

册

1兾2

,

(29)

(30)

Here, observation angle ␪v is replaced by the angle ofthe normal vector, ␪n.
6. Optical Reconstruction of Three-Dimensional
Objects

First, I fabricated several holograms of the same hexagonal prism with which to determine the value ofparameter ␥. Figure 10 shows the optical reconstruction of three holograms. The reconstructed image ofthe hologram without compensation for brightness isshown in Fig. 10(a). The left-hand surface of the hexagonal prism, which has the largest angle ␪n, is thebrightest of the object surfaces. As shown in Fig.
10(b), the hologram with compensation in ␥ ⫽ 0 is
contrasted to that in Fig. 10(a). Here, remember thatcompensation in ␥ ⫽ 0 leads to unlimited compensation. Therefore the surface that forms a large angle
with the hologram is dark as a result of overcompensation. Figure 10(c) is also applicable to a hexagonal
prism whose brightness is compensated for by␥ ⫽ 0.5. Differences of brightness disappear by propercompensation for brightness, which dissolves borders
between surfaces.
Figure 11 shows optical reconstruction of 3-D objects whose brightness is completely compensated forat ␥ ⫽ 0.5. In addition, the surfaces are shaded; theamplitudes of the surfaces are determined by use of
Eq. (29) in given virtual illumination and surrounding light. Arrows and numbers in Fig. 11 indicate

Fig. 11. Optical reconstructions of 3-D objects shaded with illumination light. Cubes are illuminated from the upper right in(a) le ⫽ 0 and from the upper left in (b) le ⫽ 0.7; a hexagonal prism 共le ⫽ 0.5兲 is shown in (c). Brightnesses of objects are all compensatedfor at ␥ ⫽ 0.5. Arrows and numbers in parentheses define the illumination vector in global coordinates.
1 August 2005 兾 Vol. 44, No. 22 兾 APPLIED OPTICS

4613

illumination vectors in global coordinates. As expected from the vectors, object surfaces are shaded in
the reconstructions.
7. Discussion

The computation time in the proposed techniques isgiven by rotational transformation of the surfaces ofthe object. According to Ref. 13, the computation timeof rotational transformation is dominated by FFTs,
and FFT operation must be executed twice to rotate aplanar surface. However, most of the inverse FFTs of
Eq. (10) can be omitted from calculating the total
field; just an inverse FFT operation is necessary tocreate a hologram because the translational propagation of the field ᏼd兵 其 can be carried out in Fourier
space. In the synthesis of holograms described in previous sections, the method of the angular spectrum ofplane waves14 is used for the operation of the propagation. Therefore the total field of Eq. (5) on the hologram is expressed by兵其ĥ(x̂, ŷ) ⫽ Ᏺ⫺1 兺 Ĥn(ûn, v̂n)exp[i2␲ŵ(ûn, v̂n)dn] ,n(31)
where dn is the distance between the 共x̂n, ŷn, 0兲 planeof the parallel coordinates and the hologram. Thus
the number of times a FFT is executed is N ⫹ 1, to
calculate the total field of an object composed of N
pieces of planar surface. As a result, one FFT兾surface
is approximately estimated as the computational costin the proposed techniques.
8. Conclusion

Full parallax computer-generated holograms of
three-dimensional surface objects were synthesized
by use of a wave-optical method. In this method, anobject is composed of some planar surfaces, and acomplex function defined for each surface retains
such properties as shape, texture, and brightness.
The fields emitted from the tilted surfaces are calculated by use of the rotational transformation of theproperty function and totaled on the hologram.
When surfaces build an object, the change ofbrightness that depends on the angle of view causes
unexpected shading of the surface. A theoreticalmodel with which to predict the brightness of thereconstructed surface and prevent unexpected shading has been proposed. This technique allows theobject to be shaded as one intends. Finally, opticalreconstructions of holograms synthesized by use of
the proposed techniques have been demonstrated toverify the validity of the methods.

4614

APPLIED OPTICS 兾 Vol. 44, No. 22 兾 1 August 2005

This study is partly supported by the Kansai University High Technology Research Center and inpart by Kansai University research grants, including a grant-in-aid for encouragement of scientists,
in 2003.
References
1. J. P. Waters, “Holographic image synthesis utilizing theoretical methods,” Appl. Phys. Lett. 9, 405– 407 (1966).
2. A. D. Stein, Z. Wang, and J. J. S. Leigh, “Computer-generated
holograms: a simplified ray-tracing approach,” Comput. Phys.
6, 389 –392 (1992).
3. M. Lucente, “Interactive computation of holograms using alook-up table,” J. Electron. Imag. 2, 28 –34 (1993).
4. J. L. Juárez-Pérez, A. Olivares-Pérez, and R. Berriel-Valdos,“Nonredundant calculation for creating digital Fresnel holograms,” Appl. Opt. 36, 7437–7443 (1997).
5. H. Yoshikawa, S. Iwase, and T. Oneda, “Fast computation of
Fresnel holograms employing difference,” in Practical Holography XIV and Holographic Materials VI, S. A. Benton, S. H.
Stevenson, and J. T. Trout, eds., Proc. SPIE 3956, 48 –55(2000).
6. K. Matsushima and M. Takai, “Recurrence formulas for fastcreation of synthetic three-dimensional holograms,” Appl. Opt.
39, 6587– 6594 (2000).
7. A. Ritter, J. Böttger, O. Deussen, M. König, and T. Strothotte,“Hardware-based rendering of full-parallax synthetic holograms,” Appl. Opt. 38, 1364 –1369 (1999).
8. T. Ito, H. Eldeib, K. Yoshida, S. Takahashi, T. Yabe, and T.
Kunugi, “Special purpose computer for holography HORN-2,”
Comput. Phys. Commun. 93, 13–20 (1996).
9. D. Leseberg and C. Frère, “Computer-generated holograms of
3-D objects composed of tilted planar segments,” Appl. Opt. 27,3020 –3024 (1988).
10. T. Tommasi and B. Bianco, “Computer-generated holograms oftilted planes by a spatial frequency approach,” J. Opt. Soc. Am.
A 10, 299 –305 (1993).
11. D. Leseberg, “Computer-generated three-dimensional image
holograms,” Appl. Opt. 31, 223–229 (1992).
12. N. Delen and B. Hooker, “Free-space beam propagation between arbitrarily oriented planes based on full diffraction theory: a fast Fourier transform approach,” J. Opt. Soc. Am. A 15,
857– 867 (1998).
13. K. Matsushima, H. Schimmel, and F. Wyrowski, “Fast calculation method for optical diffraction on tilted planes by use of
the angular spectrum of plane waves,” J. Opt. Soc. Am. A 20,
1755–1762 (2003).
14. J. W. Goodman, Introduction to Fourier Optics, 2nd ed.
(McGraw-Hill, 1996), Chap. 3.10.
15. R. Bräuer, F. Wyrowski, and O. Bryngdahl, “Diffusers in digital holography,” J. Opt. Soc. Am. A 8, 572–578 (1991).
16. K. Matsushima and A. Kondoh, “Wave optical algorithm forcreating digitally synthetic holograms of three-dimensionalsurface objects,” in Practical Holography XVII and Holographic Materials IX, T. H. Jeong and S. H. Stevenson, eds.,
Proc. SPIE 5005, 190 –197 (2003).
17. K. Matsushima and A. Joko, “A high-resolution printer for
fabricating computer-generated display holograms (in Japanese),” J. Inst. Image Inf. Television Eng. 56, 1989 –1994(2002).