**Single-Shot Phase-Shifting Digital Holography Based on the Spatial Carrier Interferometry and Its Tolerance Analysis**

Yasuhiro Harada, Aizuddin Wan and Hiroyasu Sone *Kitami Institute of Technology Japan*

### **1. Introduction**

90 Advanced Holography – Metrology and Imaging

Shemirani, M.B., Wei Mao Panicker, R.A. & Kahn, J.M. (2009). Principal Modes in Graded-

Skarman, B.,Becker, J. &Wozniak, K. (1996). Simultaneous 3D-PIV and temperature

Sokkar, T. Z. N., El-Morsy, M. A. & Wahba, H. H. (2008). Automatic fringe analysis of the induced anisotropy of bent optical bres, *Opt. Commun.*, Vol. 281, pp. 1915–1923 Sweeney, D. W. & Vest, C. M. (1973). Reconstruction of three dimensional refractive index elds from multidirectional interferometric data, *Appl. Opt.*, Vol. 12, pp. 2649–2663 Vollertsen, F. & Wochnowski, C. (2004). UV-laser assisted Fabrication of integrated-optical Waveguides, *CIRP Annals - Manufacturing Technology*, Vol. 53, pp. 199-202 Wahba, H. H. & Kreis, Th. (2009c). Digital holographic interferometric characterization of

Wahba, H. H. & Kreis, Th. (2009a). Characterization of graded index optical fibers by digital

Wahba, H. H. & Kreis, Th. (2009b). Characterization of optical fibers by digital holographic interferometry'', *Proc. SPIE 7389*, 73890K, Munich, Germany, June, 2009 Wochnowski, C. , Metev, S. and Sepold, G.(2000). UV-laser-assisted Modification of the

Yamaguchi, I. & Zhang, T. (1997). Phase-shifting digital holography, *Opt. Lett.*, Vol. 22, pp.

Yassien, K., Agour, M., V. Koyplow, C. & Dessouky, H. (2010). On the Digital Holographic

optical Properties of Polymethylmethacrylate,"*Appl. Surf. Sci.*, Vol.154–155 , pp.

Interferometry of Fibrous Materials, I: Optical Properties of Polymer and Optical

bent optical fibers, *J. Opt. A: Pure Appl. Opt.*, Vol. 11, pp.105407

holographic interferometry, *Appl. Opt.*, Vol. 48, pp. 1573-1582

Fibers, *Opt. and Lasers Eng.*, Vol. 48, pp.555

*Lightwave Tech.*, Vol. 27, pp. 1248

*Instrum.*, Vol. 7, pp. 1-6

706-711

1268-1270

Index Multimode Fiber in Presence of Spatial- and Polarization-Mode Coupling, *J.* 

measurements using a new CCD-based holographic interferometer, *Flow Meas.* 

In digital holography (Kreis, 2005; Schnars & Jueptner, 2005), holograms are recorded using an image sensor such as CCD cameras, and reconstruction of images of the object is performed by a computer. From this, digital holography has following attractive features: A wet process for developing recording media is not required; qualitative evaluation is easy for three-dimensional images of objects; and focused images of three-dimensional objects at a desired depth can be instantaneously acquired without a mechanical focusing process. However, this technology suffers from lower resolution of image sensors than photographic materials, and therefore, object size to be recorded is restricted to be large and far from the image sensor in off-axis recording geometry. These limitations can be moderated by introducing the phase-shifting method (Schreiber et al., 2007) that enables us to derive complex amplitude of the object wave on the image sensor directly from three or more phase-shifted holograms recorded in the in-line geometry (Yamaguchi & Zhang, 1997). However, since this method requires to record at least three or four phase-shifted hologram, it is useless for instantaneous measurement of the dynamic object.

In order to resolve this problem, various single-shot technologies with phase-shifting features in subsequent holograms are proposed(Awatsuji et al., 2004; Nomura et al., 2006; Toge et al., 2008; Wyant, 2003). All of these methods are based on space-division multiplexing of the phase-shifted holograms in the single recording but device and technique used are different each other. In the former three methods, special micro-element arrays of phase retarders (Awatsuji et al., 2004) or polarizers (Nomura et al., 2006; Wyant, 2003) suitable to a CCD array for recording are strictly required, and therefore, any satisfactory result for practical instrumentations has not been reported until now. On the other hand, the last method only require an introduction of off-axis plane wave illumination to an image sensor for recording single holograms (Toge et al., 2008). From its simplicity, this method seems to be one of most practical candidate for single-shot digital holography. However, it should be point out that we must align the incident angle of the reference wave plane wave so that the relative phase difference between adjacent pixels of the image sensor becomes 2*π*/3 or *π*/2 radian. This means that very rigorous alignment in reference wave illumination is required. Furthermore, this requirement may be more serious problem in the case of color digital holography where multicolor laser light is utilized. In this situation, optical alignment should be independently

**2.1 Hologram intensity distribution and its spectral distributions**

on the Spatial Carrier Interferometry and Its Tolerance Analysis

*I*(*x*, *y*) = |*UR*(*x*, *y*) + *UO*(*x*, *y*)|

where

regard only to *x* axis gives

shows its schematic representation.

given in ideal case (Gaskill, 1978) by

*<sup>S</sup>*˜(*ξ*, *<sup>y</sup>*) = <sup>F</sup>*<sup>x</sup>* [*s*(*x*, *<sup>y</sup>*)] <sup>=</sup> ∑

Now let us consider mathematical expression of interference fringe on the image sensor and its sampled version –digital hologram intensity– obtained by sampling detection by the image

<sup>93</sup> Single-Shot Phase-Shifting Digital Holography Based

2 <sup>=</sup> *<sup>A</sup>*<sup>2</sup> <sup>+</sup> *<sup>a</sup>*2(*x*, *<sup>y</sup>*) + <sup>2</sup>*Aa*(*x*, *<sup>y</sup>*) cos {2*πξ*0*<sup>x</sup>* <sup>+</sup> *<sup>φ</sup>*(*x*, *<sup>y</sup>*)}

respectively. *ξ*<sup>0</sup> means a carrier frequency in the interferogram introduced by the oblique incidence of the plane reference wave. One-dimensional Fourier transform of *I*(*x*, *y*) with

where ˜*I*(*ξ*, *<sup>y</sup>*) = <sup>F</sup>*<sup>x</sup>* [*I*(*x*, *<sup>y</sup>*)], *<sup>B</sup>*(*ξ*, *<sup>y</sup>*) = <sup>F</sup>*<sup>x</sup>* [*B*(*x*, *<sup>y</sup>*)], and *<sup>C</sup>*(*ξ*, *<sup>y</sup>*) = <sup>F</sup>*<sup>x</sup>* [*c*(*x*, *<sup>y</sup>*)], respectively. From this expression, it is clear that information of the object wave is mainly stored in the term *C*(*ξ*, *y*) shifted by the carrier frequency *ξ*<sup>0</sup> = sin *θ*/*λ* in the spectral domain. Figure 2(a)

In digital holography, hologram intensity distribution is detected by an image sensor where the detection of light intensity is performed as spatial sampling by the photodiode in the image sensor. Therefore image sensor can be considered as the sampling array with interval same as the pixel interval. Therefore a mathematical expression for this sampling array can

exp (2*πndξ*) *<sup>δ</sup>*(*<sup>y</sup>* − *md*) = ∑

*<sup>δ</sup>*(*<sup>y</sup>* − *md*) = ∑

*UO*(*x*, *y*) = *a*(*x*, *y*) exp {*iφ*(*x*, *y*)} , (2) *UR*(*x*, *y*) = *A* exp {*i*2*πξ*0*x*} , (3)

= *b*(*x*, *y*) + *c*(*x*, *y*) exp (*i*2*πξ*0*x*) + *c*∗(*x*, *y*) exp (−*i*2*πξ*0*x*) (4)

˜*I*(*ξ*, *<sup>y</sup>*) ≡ F*<sup>x</sup>* [*I*(*x*, *<sup>y</sup>*)] <sup>=</sup> *<sup>B</sup>* (*ξ*, *<sup>y</sup>*) <sup>+</sup> *<sup>C</sup>*(*<sup>ξ</sup>* <sup>−</sup> *<sup>ξ</sup>*0, *<sup>y</sup>*) + *<sup>C</sup>*∗(*<sup>ξ</sup>* <sup>+</sup> *<sup>ξ</sup>*0, *<sup>y</sup>*) (8)

*b*(*x*, *y*) = *A*<sup>2</sup> + *a*2(*x*, *y*), (5)

<sup>2</sup> *Aa*(*x*, *<sup>y</sup>*) exp {*iφ*(*x*, *<sup>y</sup>*)} , (6)

*<sup>λ</sup>* , (7)

*δ*(*x* − *nd*)*δ*(*y* − *md*), (9)

exp *π n ξN ξ* 

*δ*(*y* − *md*),

*δ* (*ξ* − *n*2*ξN*) *δ*(*y* − *md*), (10)

*<sup>m</sup>* <sup>∑</sup>*<sup>n</sup>*

*<sup>m</sup>* <sup>∑</sup>*<sup>n</sup>*

sensor. If we represent object and reference waves at the image sensor plane by

respectively, the hologram intensity, *I*(*x*, *y*), at the image sensor is given by

*<sup>c</sup>*(*x*, *<sup>y</sup>*) = <sup>1</sup>

**2.2 Image sensor as a sampling array and its spectral distributions**

*s*(*x*, *y*) = ∑

and its one-dimensional Fourier transform is given by

= ∑ *<sup>m</sup>* <sup>∑</sup>*<sup>n</sup> δ <sup>ξ</sup>* <sup>−</sup> *<sup>n</sup> d* 

*<sup>m</sup>* <sup>∑</sup>*<sup>n</sup>*

*<sup>m</sup>* <sup>∑</sup>*<sup>n</sup>*

*<sup>ξ</sup>*<sup>0</sup> <sup>=</sup> sin *<sup>θ</sup>*

performed for each wavelength but there are no description about care of the angle of incidence of each reference wave (Toge et al., 2008).

In the present chapter, we concentrate our attention on the problem of the tolerance of the incidence angle of the plane reference wave in the method proposed by (Toge et al., 2008). We theoretically analyze the problem in the Fourier domain using the principle of single-shot phase-shifting digital holography based on the spatial carrier interferometry (Harada et al., 2008; 2009). From our theoretical analysis, importance in the alignment of the plane reference wave is brought to light and special procedure to be applied under the situation, where the required condition for the incident angle of the plane reference wave is not satisfied, is given.

#### **2. Theoretical analysis of the single-shot scheme**

Figure 1 (a) shows a schematic diagram of optical geometry for recoding holograms in the method proposed by (Toge et al., 2008). It is just a off-axis geometry of digital holography where *θ* is an incident angle of the plane reference wave, *λ* the wavelength of the light in the holographic recording. We assume the interval of photodiodes in an image sensor is *d*. If the incident angle, *θ*, of the reference plane wave satisfies the equation given by

$$\frac{\lambda}{\sin \theta} = 3d \quad \text{or} \quad \sin \theta = \lambda / 3d \tag{1}$$

phase of the reference wave along a surface of the image sensor changes periodically with spatial interval of 3*d*. Therefore, the single-shot fringe pattern –hologram– recorded by the image sensor can be considered as a sum of three phase-shifted sub-holograms with constant phase difference of 2*π*/3 in the reference waves as shown in Fig. 1 (b). Therefore, after interpolating the missing data point in each sub-hologram, complex amplitude of the object wave at the image sensor can be reconstructed by using three-step algorithm of phase-shifting interferometry (Schreiber et al., 2007).

Fig. 1. Schematic diagrams of recording geometry (a) and resultant phase distributions at the image sensor under the condition of *λ*/ sin *θ* = 3*d* (b).

### **2.1 Hologram intensity distribution and its spectral distributions**

Now let us consider mathematical expression of interference fringe on the image sensor and its sampled version –digital hologram intensity– obtained by sampling detection by the image sensor. If we represent object and reference waves at the image sensor plane by

$$\mathcal{U}\_{\mathcal{O}}(\mathbf{x}, \mathcal{y}) = a(\mathbf{x}, \mathcal{y}) \exp\left\{i\phi(\mathbf{x}, \mathcal{y})\right\},\tag{2}$$

$$\mathcal{U}\_{\mathbb{R}}(\mathbf{x}, \mathbf{y}) = A \exp\left\{i 2\pi \mathfrak{f}\_0 \mathbf{x}\right\},\tag{3}$$

respectively, the hologram intensity, *I*(*x*, *y*), at the image sensor is given by

$$\begin{split} I(\mathbf{x},\mathbf{y}) &= |\mathcal{U}\_{R}(\mathbf{x},\mathbf{y}) + \mathcal{U}\_{O}(\mathbf{x},\mathbf{y})|^{2} \\ &= A^{2} + a^{2}(\mathbf{x},\mathbf{y}) + 2Aa(\mathbf{x},\mathbf{y})\cos\left\{2\pi\mathfrak{T}\_{0}\mathbf{x} + \phi(\mathbf{x},\mathbf{y})\right\} \\ &= b(\mathbf{x},\mathbf{y}) + c(\mathbf{x},\mathbf{y})\exp\left(i2\pi\mathfrak{T}\_{0}\mathbf{x}\right) + c^{\*}(\mathbf{x},\mathbf{y})\exp\left(-i2\pi\mathfrak{T}\_{0}\mathbf{x}\right) \end{split} \tag{4}$$

where

2 Will-be-set-by-IN-TECH

performed for each wavelength but there are no description about care of the angle of

In the present chapter, we concentrate our attention on the problem of the tolerance of the incidence angle of the plane reference wave in the method proposed by (Toge et al., 2008). We theoretically analyze the problem in the Fourier domain using the principle of single-shot phase-shifting digital holography based on the spatial carrier interferometry (Harada et al., 2008; 2009). From our theoretical analysis, importance in the alignment of the plane reference wave is brought to light and special procedure to be applied under the situation, where the required condition for the incident angle of the plane reference wave is not satisfied, is given.

Figure 1 (a) shows a schematic diagram of optical geometry for recoding holograms in the method proposed by (Toge et al., 2008). It is just a off-axis geometry of digital holography where *θ* is an incident angle of the plane reference wave, *λ* the wavelength of the light in the holographic recording. We assume the interval of photodiodes in an image sensor is *d*. If the

phase of the reference wave along a surface of the image sensor changes periodically with spatial interval of 3*d*. Therefore, the single-shot fringe pattern –hologram– recorded by the image sensor can be considered as a sum of three phase-shifted sub-holograms with constant phase difference of 2*π*/3 in the reference waves as shown in Fig. 1 (b). Therefore, after interpolating the missing data point in each sub-hologram, complex amplitude of the object wave at the image sensor can be reconstructed by using three-step algorithm of phase-shifting

sin *<sup>θ</sup>* <sup>=</sup> <sup>3</sup>*<sup>d</sup>* or sin *<sup>θ</sup>* <sup>=</sup> *<sup>λ</sup>*/3*<sup>d</sup>* (1)

incident angle, *θ*, of the reference plane wave satisfies the equation given by

Image sensor

pixel

*d*

Fig. 1. Schematic diagrams of recording geometry (a) and resultant phase distributions at the

*z*

*d*

*λ*/ sin *θ* = 3*d*

0

2π/3 4π/3 3*d*

3*d*

4π/3

4π/3

4π/3

2π/3

2π/3

2π/3

2π/3

0

0

= ++

3*d*

(b) resultant phase distributions at the image sensor under the condition of

0

0

*x*

*λ*

incidence of each reference wave (Toge et al., 2008).

**2. Theoretical analysis of the single-shot scheme**

interferometry (Schreiber et al., 2007).

Object wave

θ

λ

(a) Recording geometry

image sensor under the condition of *λ*/ sin *θ* = 3*d* (b).

Reference plane wave

Object

$$b(\mathbf{x}, y) = A^2 + a^2(\mathbf{x}, y),\tag{5}$$

$$\mathcal{L}(\mathbf{x}, y) = \frac{1}{2} A a(\mathbf{x}, y) \exp\left\{i\phi(\mathbf{x}, y)\right\},\tag{6}$$

$$
\xi\_0 = \frac{\sin \theta}{\lambda},
\tag{7}
$$

respectively. *ξ*<sup>0</sup> means a carrier frequency in the interferogram introduced by the oblique incidence of the plane reference wave. One-dimensional Fourier transform of *I*(*x*, *y*) with regard only to *x* axis gives

$$\tilde{I}(\tilde{\boldsymbol{\xi}},\boldsymbol{y}) \equiv \mathcal{F}\_{\boldsymbol{x}}\left[I(\boldsymbol{x},\boldsymbol{y})\right] = B\left(\tilde{\boldsymbol{\xi}},\boldsymbol{y}\right) + \mathcal{C}(\tilde{\boldsymbol{\xi}} - \tilde{\boldsymbol{\xi}}\_{0},\boldsymbol{y}) + \mathcal{C}^\*(\tilde{\boldsymbol{\xi}} + \tilde{\boldsymbol{\xi}}\_{0},\boldsymbol{y})\tag{8}$$

where ˜*I*(*ξ*, *<sup>y</sup>*) = <sup>F</sup>*<sup>x</sup>* [*I*(*x*, *<sup>y</sup>*)], *<sup>B</sup>*(*ξ*, *<sup>y</sup>*) = <sup>F</sup>*<sup>x</sup>* [*B*(*x*, *<sup>y</sup>*)], and *<sup>C</sup>*(*ξ*, *<sup>y</sup>*) = <sup>F</sup>*<sup>x</sup>* [*c*(*x*, *<sup>y</sup>*)], respectively. From this expression, it is clear that information of the object wave is mainly stored in the term *C*(*ξ*, *y*) shifted by the carrier frequency *ξ*<sup>0</sup> = sin *θ*/*λ* in the spectral domain. Figure 2(a) shows its schematic representation.

#### **2.2 Image sensor as a sampling array and its spectral distributions**

In digital holography, hologram intensity distribution is detected by an image sensor where the detection of light intensity is performed as spatial sampling by the photodiode in the image sensor. Therefore image sensor can be considered as the sampling array with interval same as the pixel interval. Therefore a mathematical expression for this sampling array can given in ideal case (Gaskill, 1978) by

$$\delta s(\mathbf{x}, y) = \sum\_{m} \sum\_{n} \delta(\mathbf{x} - nd) \delta(y - md)\_{\prime} \tag{9}$$

and its one-dimensional Fourier transform is given by

$$\begin{aligned} \tilde{S}(\boldsymbol{\xi}, \boldsymbol{y}) &= \mathcal{F}\_{\boldsymbol{X}} \left[ s(\boldsymbol{x}, \boldsymbol{y}) \right] = \sum\_{m} \sum\_{n} \exp \left( 2 \pi n d \boldsymbol{\xi} \right) \delta(\boldsymbol{y} - md) = \sum\_{m} \sum\_{n} \exp \left( \pi \frac{n}{\tilde{\xi}\_{N}} \boldsymbol{\xi} \right) \delta(\boldsymbol{y} - md), \\ &= \sum\_{m} \sum\_{n} \delta \left( \tilde{\xi} - \frac{n}{d} \right) \delta(\boldsymbol{y} - md) = \sum\_{m} \sum\_{n} \delta \left( \tilde{\xi} - n \mathfrak{L}\_{N}^{\mathsf{T}} \right) \delta(\boldsymbol{y} - md), \end{aligned} \tag{10}$$

respectively. Schematic representation of absolute value of the spectral distribution *S*˜(*j*) is shown in Fig. 2(c). Please note that, in addition to the widen interval of 3*d*, the shift of the detecting point *jd* is introduced in the spatial domain representation of Eq. (11). These two differences result in reduced interval of spectral array with 1/3*d* = 2*ξN*/3 and new phase term whose phase value depends on the number *j* and spectral order *n*� in the spectral domain representation of Eqs. (12) and (13). For the sub-detector array *s*(0)(*x*, *y*), the additional phase for spectral array is 0 for all impulsive spectrum of the order *n*� since *j* = 0. In contrast to this, for the sub-detector array *s*(1)(*x*, *y*) and *s*(2)(*x*, *y*) result in the impulsive spectral arrays

<sup>95</sup> Single-Shot Phase-Shifting Digital Holography Based

be shown in the next subsection this phase term plays very important role in the method for

Let us consider the digital hologram detected by the image sensor of Eq. (9) and its sub-holograms which correspond to the hologram considered to be detected by individual sparse sampling array, *s*(*j*)(*x*, *y*), given by Eq. (11). Mathematical expression of the digital hologram, *Is*(*x*, *y*), detected by the image sensor is simply given by a product of the hologram intensity, *I*(*x*, *y*), given by Eq. (4), and the sampling array, *s*(*x*, *y*), given by Eq. (9). Therefore, spectral expression of the digital hologram, ˜*I*(*ξ*, *y*) can be given by the convolution integral

In the same manner and with help of linearity shown in Eq. (11), mathematical expressions for

*<sup>s</sup>* (*x*, *<sup>y</sup>*) = *<sup>I</sup>*(*x*, *<sup>y</sup>*) · *<sup>s</sup>*(*j*)

*<sup>s</sup>* (*ξ*, *<sup>y</sup>*) = ˜*I*(*ξ*, *<sup>y</sup>*) <sup>∗</sup>*<sup>ξ</sup> <sup>S</sup>*˜(*j*)

By substituting Eqs. (8) and (13) into Eq. (16), we obtain important expression for spectral

<sup>∑</sup> *m*

This equation implies that spectral distribution of the hologram intensity is replicated at

method Toge et al. (2008), we simply observe spectra falling into the Nyquist base-band of the sub-sampling array, i.e., |*ξ*| < *ξN*/3. Now we are a position to investigate the properties of

2*ξN*/3 and modulated in phase by −(2*πn*�

*<sup>ξ</sup>* <sup>−</sup> *<sup>n</sup>*� <sup>2</sup>*ξ<sup>N</sup>*

<sup>3</sup> , *<sup>y</sup>* + *C* 

. According to the interpolating procedure for missing data points in the

˜*Is*(*ξ*, *<sup>y</sup>*) = ˜*I*(*ξ*, *<sup>y</sup>*) <sup>∗</sup>*<sup>ξ</sup> <sup>S</sup>*˜(*ξ*, *<sup>y</sup>*), (14)

(*x*, *y*), (15)

*<sup>ξ</sup>* <sup>−</sup> *<sup>n</sup>*� <sup>2</sup>*ξ<sup>N</sup>*

*δ*(*y* − *md*) (*i* = 0, 1, 2) (17)

<sup>3</sup> <sup>−</sup> *<sup>ξ</sup>*0, *<sup>y</sup>*

/3)*j* depending on the

(*ξ*, *y*). (16)

*<sup>s</sup>* (*x*, *y*) (*i* = 0, 1, 2), and their one-dimensional Fourier transform,

, respectively. As will

*<sup>S</sup>*˜(1)(*ξ*, *<sup>y</sup>*) and *<sup>S</sup>*˜(2)(*ξ*, *<sup>y</sup>*) with the phase of <sup>−</sup>(2*π*/3)*n*� and <sup>−</sup>(4*π*/3)*n*�

on the Spatial Carrier Interferometry and Its Tolerance Analysis

single-shot phase-shifting digital holography (Harada et al., 2008; 2009).

where ∗*<sup>ξ</sup>* stands the convolution integral with regard to the variable *ξ*.

(*j*)

exp

*<sup>ξ</sup>* <sup>−</sup> *<sup>n</sup>*� <sup>2</sup>*ξ<sup>N</sup>*

the spectrum for various condition of the reference plane wave.

*I* (*j*)

˜*I* (*j*)

<sup>−</sup>*i*2*πn*� *<sup>j</sup>* 3 *<sup>B</sup>* 

<sup>3</sup> <sup>+</sup> *<sup>ξ</sup>*0, *<sup>y</sup>*

between their Fourier spectra as follows:

distribution of the sub-holograms as follows:

<sup>3</sup> ∑ *n*�

+*C*∗ 

the digital sub-holograms, *I*

*<sup>s</sup>* (*ξ*, *<sup>y</sup>*) = <sup>2</sup>*ξ<sup>N</sup>*

the spatial frequency of *n*�

replication order *n*�

*<sup>s</sup>* (*ξ*, *y*) are given by

˜*I* (*j*)

˜*I* (*j*)

**2.3 Digital holograms detected by the image sensor and its spectral distributions**

where *ξ<sup>N</sup>* = 1/2*d* is the Nyquist frequency of the image sensor. Figure 2(b) shows its schematic representation.

Fig. 2. Spectral distribution of interference fringes and pixel arrays of the image sensor. (a) spectral distribution of the interference fringe at the image sensor, (b) spectral distribution of a pixel array of the image sensor with spatial interval of *d*, and (c) that of a pixel array with the spatial interval of 3*d*, respectively. Here *ξ<sup>N</sup>* = 1/2*d* is Nyquist frequency of the original pixel array of the image sensor.

Since the original sampling array of the image sensor with the interval of *d* in *x* direction can be recognized as a sum of sparse sampling arrays with an extended interval of 3*d* and spatial shift of detecting point by *d*, Eqs.(9) and (10) can be rewritten as

$$s(\mathbf{x}, y) = \sum\_{j=0}^{2} s^{(j)}(\mathbf{x}, y) = \sum\_{j=0}^{2} \sum\_{n'} \sum\_{m} \delta(\mathbf{x} - n' \mathbf{3}d - jd) \delta(y - md), \tag{11}$$

and

$$\tilde{S}(\tilde{\xi}, y) = \sum\_{j=0}^{2} \tilde{S}^{(j)}(\tilde{\xi}, y) = \frac{1}{3d} \sum\_{j=0}^{2} \sum\_{n'} \sum\_{m} \exp\left[-i2\pi n' \frac{j}{3}\right] \delta\left(\tilde{\xi} - \frac{n'}{3d}\right) \delta\left(y - md\right), \tag{12}$$

$$=\frac{2\xi\_N}{3}\sum\_{j=0}^{2}\sum\_{n'}\sum\_{m}\exp\left[-i2\pi n'\frac{j}{3}\right]\delta\left(\xi-n'\frac{2\xi\_N}{3}\right)\delta\left(y-md\right),\tag{13}$$

4 Will-be-set-by-IN-TECH

where *ξ<sup>N</sup>* = 1/2*d* is the Nyquist frequency of the image sensor. Figure 2(b) shows its

A (ξ, y)



Fig. 2. Spectral distribution of interference fringes and pixel arrays of the image sensor. (a) spectral distribution of the interference fringe at the image sensor, (b) spectral distribution of a pixel array of the image sensor with spatial interval of *d*, and (c) that of a pixel array with the spatial interval of 3*d*, respectively. Here *ξ<sup>N</sup>* = 1/2*d* is Nyquist frequency of the original

Since the original sampling array of the image sensor with the interval of *d* in *x* direction can be recognized as a sum of sparse sampling arrays with an extended interval of 3*d* and spatial

*δ*(*x* − *n*�

<sup>−</sup>*i*2*πn*� *<sup>j</sup>* 3 *δ <sup>ξ</sup>* <sup>−</sup> *<sup>n</sup>*� 3*d* 

<sup>−</sup>*i*2*πn*� *<sup>j</sup>* 3 *δ* 

2 ∑ *j*=0 ∑ *n*� ∑ *m*

shift of detecting point by *d*, Eqs.(9) and (10) can be rewritten as

3*d*

<sup>=</sup> <sup>2</sup>*ξ<sup>N</sup>* 3

(*x*, *y*) =

2 ∑ *j*=0 ∑ *n*� ∑ *m* exp 

> 2 ∑ *j*=0 ∑ *n*� ∑ *m* exp

2 ∑ *j*=0 *s* (*j*)

(*ξ*, *<sup>y</sup>*) = <sup>1</sup>

0 ξ<sup>0</sup>

C (ξ-ξ0, y)

ξ

ξ

ξ

3*d* − *jd*)*δ*(*y* − *md*), (11)

*<sup>ξ</sup>* <sup>−</sup> *<sup>n</sup>*� <sup>2</sup>*ξ<sup>N</sup>* 3

*δ* (*y* − *md*), (12)

*δ* (*y* − *md*), (13)

C\* (ξ+ξ0, y)


schematic representation.

(a)

(b)

(c)

pixel array of the image sensor.

and

*S*˜(*ξ*, *y*) =

*s*(*x*, *y*) =

2 ∑ *j*=0 *S*˜(*j*) respectively. Schematic representation of absolute value of the spectral distribution *S*˜(*j*) is shown in Fig. 2(c). Please note that, in addition to the widen interval of 3*d*, the shift of the detecting point *jd* is introduced in the spatial domain representation of Eq. (11). These two differences result in reduced interval of spectral array with 1/3*d* = 2*ξN*/3 and new phase term whose phase value depends on the number *j* and spectral order *n*� in the spectral domain representation of Eqs. (12) and (13). For the sub-detector array *s*(0)(*x*, *y*), the additional phase for spectral array is 0 for all impulsive spectrum of the order *n*� since *j* = 0. In contrast to this, for the sub-detector array *s*(1)(*x*, *y*) and *s*(2)(*x*, *y*) result in the impulsive spectral arrays *<sup>S</sup>*˜(1)(*ξ*, *<sup>y</sup>*) and *<sup>S</sup>*˜(2)(*ξ*, *<sup>y</sup>*) with the phase of <sup>−</sup>(2*π*/3)*n*� and <sup>−</sup>(4*π*/3)*n*� , respectively. As will be shown in the next subsection this phase term plays very important role in the method for single-shot phase-shifting digital holography (Harada et al., 2008; 2009).

### **2.3 Digital holograms detected by the image sensor and its spectral distributions**

Let us consider the digital hologram detected by the image sensor of Eq. (9) and its sub-holograms which correspond to the hologram considered to be detected by individual sparse sampling array, *s*(*j*)(*x*, *y*), given by Eq. (11). Mathematical expression of the digital hologram, *Is*(*x*, *y*), detected by the image sensor is simply given by a product of the hologram intensity, *I*(*x*, *y*), given by Eq. (4), and the sampling array, *s*(*x*, *y*), given by Eq. (9). Therefore, spectral expression of the digital hologram, ˜*I*(*ξ*, *y*) can be given by the convolution integral between their Fourier spectra as follows:

$$
\tilde{I}\_s(\xi, y) = \tilde{I}(\xi, y) \*\_{\tilde{\xi}} \tilde{S}(\xi, y),
\tag{14}
$$

where ∗*<sup>ξ</sup>* stands the convolution integral with regard to the variable *ξ*. In the same manner and with help of linearity shown in Eq. (11), mathematical expressions for the digital sub-holograms, *I* (*j*) *<sup>s</sup>* (*x*, *y*) (*i* = 0, 1, 2), and their one-dimensional Fourier transform, ˜*I* (*j*) *<sup>s</sup>* (*ξ*, *y*) are given by

$$I\_s^{(j)}(\mathbf{x}, y) = I(\mathbf{x}, y) \cdot \mathbf{s}^{(j)}(\mathbf{x}, y), \tag{15}$$

$$
\tilde{I}\_s^{(j)}(\xi, y) = \tilde{I}(\xi, y) \ast\_{\xi} \tilde{S}^{(j)}(\xi, y). \tag{16}
$$

By substituting Eqs. (8) and (13) into Eq. (16), we obtain important expression for spectral distribution of the sub-holograms as follows:

$$\begin{split} I\_{\mathsf{s}}^{(j)}(\xi,\mathsf{y}) &= \frac{2\xi\_{\mathsf{N}}}{3} \sum\_{n'} \exp\left[-i2\pi n' \frac{j}{3}\right] \left[B\left(\xi - n' \frac{2\xi\_{\mathsf{N}}}{3}, \mathsf{y}\right) + \mathcal{C}\left(\xi - n' \frac{2\xi\_{\mathsf{N}}}{3} - \xi\_{\mathsf{0}}, \mathsf{y}\right)\right.\\ &\left. + \mathcal{C}^{\*}\left(\xi - n' \frac{2\xi\_{\mathsf{N}}}{3} + \xi\_{\mathsf{0}}, \mathsf{y}\right)\right] \sum\_{m} \delta(\mathsf{y} - md) \quad (i = 0, 1, 2) \end{split} \tag{17}$$

This equation implies that spectral distribution of the hologram intensity is replicated at the spatial frequency of *n*� 2*ξN*/3 and modulated in phase by −(2*πn*� /3)*j* depending on the replication order *n*� . According to the interpolating procedure for missing data points in the method Toge et al. (2008), we simply observe spectra falling into the Nyquist base-band of the sub-sampling array, i.e., |*ξ*| < *ξN*/3. Now we are a position to investigate the properties of the spectrum for various condition of the reference plane wave.

*I*

*I*

following phase:

(1) *<sup>s</sup>* (*ξ*, *<sup>y</sup>*) = <sup>2</sup>*ξ<sup>N</sup>*

(2) *<sup>s</sup>* (*ξ*, *<sup>y</sup>*) = <sup>2</sup>*ξ<sup>N</sup>*

+*C*∗ *ξ* + 

+*C*∗ *ξ* + 

3

on the Spatial Carrier Interferometry and Its Tolerance Analysis

3

2*π* 

after the measurement of object wave.

technique Takeda et al. (1982).

considered here is given by

**3. Verification by computer simulations**

*UO*(*x*, *<sup>y</sup>*) = exp

*<sup>ξ</sup>*<sup>0</sup> <sup>−</sup> <sup>2</sup>*ξ<sup>N</sup>* 3 

*B*(*ξ*, *y*) + *C*

*<sup>ξ</sup>*<sup>0</sup> <sup>−</sup> <sup>2</sup>*ξ<sup>N</sup>* 3

*B*(*ξ*, *y*) + *C*

*<sup>ξ</sup>*<sup>0</sup> <sup>−</sup> <sup>2</sup>*ξ<sup>N</sup>* 3 , *y* exp −*i* 4*π* 3

 *ξ* − 

<sup>97</sup> Single-Shot Phase-Shifting Digital Holography Based

 , *y* exp −*i* 2*π* 3

 *ξ* − 

As shown in Eqs. (23)–(25), phase-shift relationship between these three sub holograms still remains but the components *C*(*ξ*, *y*) and *C*∗(*ξ*, *y*) are modulated by the spatial linear carrier with spatial frequency of *ξ*<sup>0</sup> − 2*ξN*/3. It comes from the misalignment of the reference plane wave. Therefore, simple application of the three-step algorithm given by Eqs. (21) and (22) will lead the wrong result. For the case of the color digital holography based on the single-shot method proposed by Sato *et al.* Toge et al. (2008), the angle of the incident satisfying Eq. (1) is inevitably different with different color of the light. Therefore, in such case, one must align the optical components for each color of the reference waves independently. This suggests that the color digital holography by Sato *et al.* Toge et al. (2008) seems to be done with very careful and technical procedures. The simple way to obtain the correct complex amplitude of the object wave at the image sensor without any careful alignment is to compensate the phase due to linear carrier with the frequency of *ξ*<sup>0</sup> − 2*ξN*/3 after unwrapping the phase obtained by using Eq. (22). More concretely, correct phase distribution is obtained by subtracting the

from the unwrapped phase at each pixel. For linear-phase compensation, we have to measure additional hologram without object and derive the frequency of the linear carrier before or

Another extreme case where Eq. (1) is not satisfied is shown in Fig. 3(c). In this case the incident angle, *θ*, is too small for components *C*s of the higher order of *n*� to fall into the Nyquist base-band, |*ξ*| ≤ *ξN*/3. Therefore, there is no phase-shifting relationship between three sub holograms obtained from single off-axis hologram, and therefore, reconstruction of object waves using Eqs. (21) and (22) will be in vain both amplitude and phase. One solution to derive the correct complex amplitude in the present case is to use of the Fourier transform

In order to demonstrate the present theoretical analysis, computer simulations are performed for each case discussed in the previous section. The object wave at the image sensor

*i.e.*, the phase object is assumed. Total number of pixels of the image sensor in the horizontal and vertical directions is 256. Therefore pixel size, *d*, is 7.8*μ*m. Figures 4 (a) and (b) show a

i2*π*(*x*<sup>2</sup> + *y*2)<sup>2</sup>

*<sup>ξ</sup>*<sup>0</sup> <sup>−</sup> <sup>2</sup>*ξ<sup>N</sup>* 3 , *y* exp *i* 2*π* 3 

*<sup>ξ</sup>*<sup>0</sup> <sup>−</sup> <sup>2</sup>*ξ<sup>N</sup>* 3 , *y* exp *i* 4*π* 3 

<sup>∑</sup> *m*

<sup>∑</sup> *m*

*nd* (*n* = 0, 1, 2, ···). (26)

(|*x*|, |*y*| ≤ 1[mm]), (27)

(28)

*δ*(*y* − *md*), (24)

*δ*(*y* − *md*). (25)

#### **2.4 Considerations**

In the present subsection, we investigate effects of incident angle of the reference plane wave in the single-shot phase-shifting digital holography from spectral point of view based on Eq. (17).

First case is the case where Eq. (1) is satisfied. Corresponding spectral distribution is shown in Fig. 3(a). In this case, spectral components of the hologram falling into the Nyquist base-band, |*ξ*| ≤ *ξN*/3, are found to be *B* of *n*� = 0, *C*<sup>∗</sup> of *n*� = +1, and *C* of *n*� = −1, and their locations are all at the origin of spectral coordinate. Therefore, the three sub hologram obtained in this case are give by

$$\overline{I\_{s}^{(0)}}(\xi,y) = \frac{2\tilde{\xi}\_{N}}{3} \left[ B(\xi,y) + \mathcal{C}(\xi,y) + \mathcal{C}^\*(\xi,y) \right] \sum\_{m} \delta(y-md),\tag{18}$$
  $\overline{I\_{s}^{(1)}}(\xi,y) = \frac{2\tilde{\xi}\_{N}}{3} \left[ B(\xi,y) + \mathcal{C}(\xi,y) \exp\left(i\frac{2\pi}{3}\right) + \mathcal{C}^\*(\xi,y) \exp\left(-i\frac{2\pi}{3}\right) \right] \sum\_{m} \delta(y-md),$ 

$$\overline{I\_s^{(2)}}(\xi, y) = \frac{2\xi\_N}{3} \left[ B(\xi, y) + \mathbb{C}(\xi, y) \exp\left(i\frac{4\pi}{3}\right) + \mathbb{C}^\*(\xi, y) \exp\left(-i\frac{4\pi}{3}\right) \right] \sum\_m \delta(y - md), \tag{20}$$

where symbol ··· denotes filter out only low-frequency component less than *ξ* = *ξN*/3. This is, in fact, performed by interpolating missing data due to sparse sampling in the space domain. As shown in Eqs.(18)—(20), important component *C* in the each sub hologram has the relative phase difference of 2*π*/3, which corresponds relative phase shift in the space domain representation of the sub holograms. Therefore we are able to simply apply the three-step algorithm in phase-shifting interferometry Schreiber et al. (2007) for calculating complex amplitude of the object wave at the image sensors, *i.e.*, amplitude *a*(*x*, *y*) and phase *φW*(*x*, *y*) (subscript *<sup>W</sup>* denotes *wrapped* within −*π* to +*π*),

$$a(\mathbf{x}, \mathbf{y}) = \frac{1}{A} \left[ \left\{ \overline{I\_s^{(1)}(\mathbf{x}, \mathbf{y})} - \overline{I\_s^{(2)}(\mathbf{x}, \mathbf{y})} \right\}^2 + \left\{ \overline{2I\_s^{(0)}(\mathbf{x}, \mathbf{y})} - \overline{I\_s^{(1)}(\mathbf{x}, \mathbf{y})} - \overline{I\_s^{(2)}(\mathbf{x}, \mathbf{y})} \right\}^2 \right]^{1/2} \tag{21}$$
 
$$\left[ \begin{array}{c} \sqrt{3} \left\{ \overline{I\_s^{(1)}(\mathbf{x}, \mathbf{y})} - \overline{I\_s^{(2)}(\mathbf{x}, \mathbf{y})} \right\} \end{array} \right]$$

$$\phi\_W(\mathbf{x}, y) = \tan^{-1} \left[ \frac{\mathbf{v} \cdot \mathbf{v} \left\{ I\_s^{(1)}(\mathbf{x}, y) - I\_s^{(1)}(\mathbf{x}, y) \right\}}{2I\_s^{(0)}(\mathbf{x}, y) - \overline{I\_s^{(1)}(\mathbf{x}, y)} - \overline{I\_s^{(2)}(\mathbf{x}, y)}} \right] \cdot \text{ (22)}$$

Next, we consider the case where Eq.(1) is not satisfied. Figure 3(b) shows the case where sin *θ* > *λ*/3*d*. Spectral components of the hologram falling into the Nyquist base-band, |*ξ*| ≤ *ξN*/3, are the same as the previous case but their locations are different from the origin. Taking account this fact, the three sub holograms obtained in this case are given by

$$\overline{I\_s^{(0)}}(\xi, y) = \frac{2\tilde{\xi}\_N}{3} \left[ B(\xi, y) + C \left\{ \xi - \left( \tilde{\xi}\_0 - \frac{2\tilde{\xi}\_N}{3} \right), y \right\} + C^\* \left\{ \xi + \left( \tilde{\xi}\_0 - \frac{2\tilde{\xi}\_N}{3} \right), y \right\} \right] \tag{23}$$

$$\times \sum\_m \delta(y - md), \quad \text{(23)}$$

6 Will-be-set-by-IN-TECH

In the present subsection, we investigate effects of incident angle of the reference plane wave in the single-shot phase-shifting digital holography from spectral point of view based on

First case is the case where Eq. (1) is satisfied. Corresponding spectral distribution is shown in Fig. 3(a). In this case, spectral components of the hologram falling into the Nyquist base-band, |*ξ*| ≤ *ξN*/3, are found to be *B* of *n*� = 0, *C*<sup>∗</sup> of *n*� = +1, and *C* of *n*� = −1, and their locations are all at the origin of spectral coordinate. Therefore, the three sub hologram obtained in this

<sup>3</sup> [*B*(*ξ*, *<sup>y</sup>*) + *<sup>C</sup>*(*ξ*, *<sup>y</sup>*) + *<sup>C</sup>*∗(*ξ*, *<sup>y</sup>*)]∑

� −*i* 2*π* 3

� −*i* 4*π* 3

(1) *<sup>s</sup>* (*x*, *<sup>y</sup>*) <sup>−</sup> *<sup>I</sup>*

(1) *<sup>s</sup>* (*x*, *<sup>y</sup>*) <sup>−</sup> *<sup>I</sup>*

(1) *<sup>s</sup>* (*x*, *<sup>y</sup>*) <sup>−</sup> *<sup>I</sup>*

*<sup>ξ</sup>*<sup>0</sup> <sup>−</sup> <sup>2</sup>*ξ<sup>N</sup>* 3 � , *y* ��

× ∑ *m*

(2) *<sup>s</sup>* (*x*, *<sup>y</sup>*)

(2) *<sup>s</sup>* (*x*, *<sup>y</sup>*)

�<sup>2</sup> �1/2

�

(2) *<sup>s</sup>* (*x*, *<sup>y</sup>*)

+ *C*∗(*ξ*, *y*) exp

+ *C*∗(*ξ*, *y*) exp

(0) *<sup>s</sup>* (*x*, *<sup>y</sup>*) <sup>−</sup> *<sup>I</sup>*

√3 � *I*

(0) *<sup>s</sup>* (*x*, *<sup>y</sup>*) <sup>−</sup> *<sup>I</sup>*

⎡ ⎢ ⎢ ⎣

2*I*

*m*

�� ∑ *m*

�� ∑ *m*

*δ*(*y* − *md*), (18)

(19)

(20)

, (21)

⎤ ⎥ ⎥ ⎦ . (22)

*δ*(*y* − *md*), (23)

*δ*(*y* − *md*),

*δ*(*y* − *md*),

*I*

*B*(*ξ*, *y*) + *C*(*ξ*, *y*) exp

*B*(*ξ*, *y*) + *C*(*ξ*, *y*) exp

*φW*(*x*, *y*) (subscript *<sup>W</sup>* denotes *wrapped* within −*π* to +*π*),

(2) *<sup>s</sup>* (*x*, *<sup>y</sup>*)

�<sup>2</sup> + � 2*I*

*φW*(*x*, *y*) = tan−<sup>1</sup>

account this fact, the three sub holograms obtained in this case are given by

� *ξ* − �

Next, we consider the case where Eq.(1) is not satisfied. Figure 3(b) shows the case where sin *θ* > *λ*/3*d*. Spectral components of the hologram falling into the Nyquist base-band, |*ξ*| ≤ *ξN*/3, are the same as the previous case but their locations are different from the origin. Taking

> *<sup>ξ</sup>*<sup>0</sup> <sup>−</sup> <sup>2</sup>*ξ<sup>N</sup>* 3 � , *y* � + *C*∗ � *ξ* + �

(1) *<sup>s</sup>* (*x*, *<sup>y</sup>*) <sup>−</sup> *<sup>I</sup>*

(0) *<sup>s</sup>* (*ξ*, *<sup>y</sup>*) = <sup>2</sup>*ξ<sup>N</sup>*

� *i* 2*π* 3 �

� *i* 4*π* 3 �

where symbol ··· denotes filter out only low-frequency component less than *ξ* = *ξN*/3. This is, in fact, performed by interpolating missing data due to sparse sampling in the space domain. As shown in Eqs.(18)—(20), important component *C* in the each sub hologram has the relative phase difference of 2*π*/3, which corresponds relative phase shift in the space domain representation of the sub holograms. Therefore we are able to simply apply the three-step algorithm in phase-shifting interferometry Schreiber et al. (2007) for calculating complex amplitude of the object wave at the image sensors, *i.e.*, amplitude *a*(*x*, *y*) and phase

**2.4 Considerations**

case are give by

(1) *<sup>s</sup>* (*ξ*, *<sup>y</sup>*) = <sup>2</sup>*ξ<sup>N</sup>*

(2) *<sup>s</sup>* (*ξ*, *<sup>y</sup>*) = <sup>2</sup>*ξ<sup>N</sup>*

*<sup>a</sup>*(*x*, *<sup>y</sup>*) = <sup>1</sup>

*I*

*A*

(0) *<sup>s</sup>* (*ξ*, *<sup>y</sup>*) = <sup>2</sup>*ξ<sup>N</sup>*

3

�

*B*(*ξ*, *y*) + *C*

�� *I*

3

3

�

�

Eq. (17).

*I*

*I*

$$\begin{split} I\_{s}^{(1)}(\xi,\mathsf{y}) &= \frac{2\mathsf{f}\_{N}}{3} \left[ B(\xi,\mathsf{y}) + \mathsf{C} \left\{ \xi - \left( \xi\_{0} - \frac{2\mathsf{f}\_{N}}{3} \right), \mathsf{y} \right\} \exp\left( i\frac{2\pi}{3} \right) \right. \\ &\left. + \mathsf{C}^{\*} \left\{ \xi + \left( \xi\_{0} - \frac{2\mathsf{f}\_{N}}{3} \right), \mathsf{y} \right\} \exp\left( -i\frac{2\pi}{3} \right) \right] \sum\_{m} \delta(\mathsf{y} - md), \end{split} \tag{24}$$

$$\begin{split} \overline{I\_{s}^{(2)}}(\xi,y) &= \frac{2\tilde{\xi}\_{N}}{3} \left[ B(\tilde{\xi},y) + \mathcal{C} \left\{ \tilde{\xi} - \left( \tilde{\xi}\_{0} - \frac{2\tilde{\xi}\_{N}}{3} \right), y \right\} \exp\left( i\frac{4\pi}{3} \right) \right. \\ &+ \mathcal{C}^{\*} \left\{ \tilde{\xi} + \left( \tilde{\xi}\_{0} - \frac{2\tilde{\xi}\_{N}}{3} \right), y \right\} \exp\left( -i\frac{4\pi}{3} \right) \right] \sum\_{m} \delta(y-md). \end{split} \tag{25}$$

As shown in Eqs. (23)–(25), phase-shift relationship between these three sub holograms still remains but the components *C*(*ξ*, *y*) and *C*∗(*ξ*, *y*) are modulated by the spatial linear carrier with spatial frequency of *ξ*<sup>0</sup> − 2*ξN*/3. It comes from the misalignment of the reference plane wave. Therefore, simple application of the three-step algorithm given by Eqs. (21) and (22) will lead the wrong result. For the case of the color digital holography based on the single-shot method proposed by Sato *et al.* Toge et al. (2008), the angle of the incident satisfying Eq. (1) is inevitably different with different color of the light. Therefore, in such case, one must align the optical components for each color of the reference waves independently. This suggests that the color digital holography by Sato *et al.* Toge et al. (2008) seems to be done with very careful and technical procedures. The simple way to obtain the correct complex amplitude of the object wave at the image sensor without any careful alignment is to compensate the phase due to linear carrier with the frequency of *ξ*<sup>0</sup> − 2*ξN*/3 after unwrapping the phase obtained by using Eq. (22). More concretely, correct phase distribution is obtained by subtracting the following phase:

$$2\pi \left[ \mathfrak{f}\_0 - \frac{2\mathfrak{f}\_N}{3} \right] \text{ } nd \qquad \left( n = 0, 1, 2, \cdots \right). \tag{26}$$

from the unwrapped phase at each pixel. For linear-phase compensation, we have to measure additional hologram without object and derive the frequency of the linear carrier before or after the measurement of object wave.

Another extreme case where Eq. (1) is not satisfied is shown in Fig. 3(c). In this case the incident angle, *θ*, is too small for components *C*s of the higher order of *n*� to fall into the Nyquist base-band, |*ξ*| ≤ *ξN*/3. Therefore, there is no phase-shifting relationship between three sub holograms obtained from single off-axis hologram, and therefore, reconstruction of object waves using Eqs. (21) and (22) will be in vain both amplitude and phase. One solution to derive the correct complex amplitude in the present case is to use of the Fourier transform technique Takeda et al. (1982).

### **3. Verification by computer simulations**

In order to demonstrate the present theoretical analysis, computer simulations are performed for each case discussed in the previous section. The object wave at the image sensor considered here is given by

$$\mathcal{U}\_{\mathcal{O}}(\mathbf{x}, \mathbf{y}) = \exp\left[\mathbf{i}2\pi(\mathbf{x}^2 + \mathbf{y}^2)^2\right] \quad (|\mathbf{x}|, |\mathbf{y}| \le \mathbf{1}[\mathbf{m}\mathbf{m}]),\tag{27}$$

$$(28)$$

*i.e.*, the phase object is assumed. Total number of pixels of the image sensor in the horizontal and vertical directions is 256. Therefore pixel size, *d*, is 7.8*μ*m. Figures 4 (a) and (b) show a

phase distribution of the object wave under test in gray level and its intensity distribution of

<sup>99</sup> Single-Shot Phase-Shifting Digital Holography Based

(a) (b)

Figures 5(a)–(c) show intensity distributions of holograms detected by the image sensor in the off-axis recording geometry shown in Fig. 1 (a) with an angle of (a) sin *θ* = *λ*/3*d*, (b) sin *θ* = 1.4*λ*/3*d* > *λ*/3*d*, and (c) sin *θ* = 0.8*λ*/6*d* < *λ*/6*d*, respectively. Corresponding spatial carrier frequency introduced by tilted reference plane waves are (a) *ξ*<sup>0</sup> = 1/3*d* = 2*ξN*/3, (b) *ξ*<sup>0</sup> = 1.4/3*d* = 2.8*ξN*/3, and (c) *ξ*<sup>0</sup> = 0.8/6*d* = 0.8*ξN*/3, respectively. Corresponding spectral relations are schematically given in Figs. 3(a), (b), and (c), respectively. In the case of the Figs.5 (a) and (b), where side-band spectra of the higher order (*n*� = ±1) falls in to the Nyquist base-band of |*ξ*| ≤ *ξN*/3 = 1/6*d* [see Figs.3(a) and (b)], intensity distributions of the detected holograms contain aliased or moiré fringes. In contrast to them, Fig. 5 (c) shows clear interference fringes with linear carrier because there are no higher order side-band spectra in

Figures 5(a1)–(a3) denote sub-holograms obtained from the Fig. 5(a) by sparsely sampling with the interval of 3*d*. Each of them shows similar coarse fringe pattern those obtained in the in-line geometry as shown in Fig. 4(b). Furthermore, their coarse fringes are shifted in their phase 2*π*/3 each other. This fact confirms that phase distribution of the object wave at the

Figures 5(b1)–(b3) show sub-holograms obtained from the Fig. 5(b) by using the same manner. The coarse fringes in these figures are different shape from those of the in-line holograms

2*ξN*/3| = |2.8*ξN*/3 − 2*ξN*/3| = 0.8*ξN*/3 due to large amount of tilt in the reference plane wave. However, it must be noted that there is the phase-shift relation between these coarse fringes with the phase value of 2*π*/3 each other. Therefore, after interpolating process we obtain three phase-shift fringes which are available to calculate a phase distribution where

Final phase distribution of the object wave at the image sensor can be obtained by subtracting

<sup>0</sup> = |*ξ*<sup>0</sup> −

<sup>0</sup> *nd* = (0.6*ξN*/3)*nd* (*n* = 1, 2, 3, . . .) are added.

<sup>0</sup> *nd* = (0.6*ξN*/3)*nd* (*n* = 1, 2, 3, . . .). This result supports our derived

image sensor can be evaluated directly by using Eq.(22) after interpolating process.

shown in Figs. 5(a1)–(a3) since there are carrier fringes with a spatial frequency of *ξ*�

Fig. 4. Object wave at the image sensor investigated in computer simulation: (a) phase distribution in gray level and (b) intensity distribution of a hologram obtained by the in-line

the Nyquist base-band [see Fig. 3(c)], and therefore, no moiré fringes.

the phase of object wave and linear phase of *ξ*�

the spatial linear phase *ξ*�

a hologram recorded by in-line geometry, respectively.

on the Spatial Carrier Interferometry and Its Tolerance Analysis

geometry.

(c) sin *θ* < *λ*/6*d*, *ξ*<sup>0</sup> < 1/6*d* = *ξN*/3

Fig. 3. Spectral distributions of resultant sub-holograms for various conditions of the incident angle of the reference plane wave. (a) sin *θ* = *λ*/3*d* and *ξ*<sup>0</sup> = 1/3*d* = 2*ξN*/3, (b) sin *θ* > *λ*/3*d* and *ξ*<sup>0</sup> > 1/3*d* = 2*ξN*/3, and (c) sin *θ* < *λ*/6*d* and *ξ*<sup>0</sup> < 1/6*d* = *ξN*/3, respectively.

8 Will-be-set-by-IN-TECH



(a) sin *θ* = *λ*/3*d*, *ξ*<sup>0</sup> = 1/3*d* = 2*ξN*/3

<sup>2</sup>ξ<sup>N</sup> <sup>ξ</sup><sup>0</sup> <sup>4</sup>ξN/3 -2ξ<sup>N</sup> -4ξN/3 -2ξN/3 <sup>0</sup> 2ξN/3 -ξ<sup>0</sup> -ξN/3 <sup>ξ</sup>N/3

(b) *λ*/2*d* > sin *θ* > *λ*/3*d*, 1/2*d* = *ξ<sup>N</sup>* > *ξ*<sup>0</sup> > 1/3*d* = 2*ξN*/3


Fig. 3. Spectral distributions of resultant sub-holograms for various conditions of the incident angle of the reference plane wave. (a) sin *θ* = *λ*/3*d* and *ξ*<sup>0</sup> = 1/3*d* = 2*ξN*/3, (b) sin *θ* > *λ*/3*d* and *ξ*<sup>0</sup> > 1/3*d* = 2*ξN*/3, and (c) sin *θ* < *λ*/6*d* and *ξ*<sup>0</sup> < 1/6*d* = *ξN*/3,



respectively.

ξ<sup>0</sup> ξ<sup>0</sup>

ξ

ξ

ξ

2ξ<sup>N</sup>

4ξN/3 ξ<sup>0</sup>

ξN/3

(c) sin *θ* < *λ*/6*d*, *ξ*<sup>0</sup> < 1/6*d* = *ξN*/3

phase distribution of the object wave under test in gray level and its intensity distribution of a hologram recorded by in-line geometry, respectively.

Fig. 4. Object wave at the image sensor investigated in computer simulation: (a) phase distribution in gray level and (b) intensity distribution of a hologram obtained by the in-line geometry.

Figures 5(a)–(c) show intensity distributions of holograms detected by the image sensor in the off-axis recording geometry shown in Fig. 1 (a) with an angle of (a) sin *θ* = *λ*/3*d*, (b) sin *θ* = 1.4*λ*/3*d* > *λ*/3*d*, and (c) sin *θ* = 0.8*λ*/6*d* < *λ*/6*d*, respectively. Corresponding spatial carrier frequency introduced by tilted reference plane waves are (a) *ξ*<sup>0</sup> = 1/3*d* = 2*ξN*/3, (b) *ξ*<sup>0</sup> = 1.4/3*d* = 2.8*ξN*/3, and (c) *ξ*<sup>0</sup> = 0.8/6*d* = 0.8*ξN*/3, respectively. Corresponding spectral relations are schematically given in Figs. 3(a), (b), and (c), respectively. In the case of the Figs.5 (a) and (b), where side-band spectra of the higher order (*n*� = ±1) falls in to the Nyquist base-band of |*ξ*| ≤ *ξN*/3 = 1/6*d* [see Figs.3(a) and (b)], intensity distributions of the detected holograms contain aliased or moiré fringes. In contrast to them, Fig. 5 (c) shows clear interference fringes with linear carrier because there are no higher order side-band spectra in the Nyquist base-band [see Fig. 3(c)], and therefore, no moiré fringes.

Figures 5(a1)–(a3) denote sub-holograms obtained from the Fig. 5(a) by sparsely sampling with the interval of 3*d*. Each of them shows similar coarse fringe pattern those obtained in the in-line geometry as shown in Fig. 4(b). Furthermore, their coarse fringes are shifted in their phase 2*π*/3 each other. This fact confirms that phase distribution of the object wave at the image sensor can be evaluated directly by using Eq.(22) after interpolating process.

Figures 5(b1)–(b3) show sub-holograms obtained from the Fig. 5(b) by using the same manner. The coarse fringes in these figures are different shape from those of the in-line holograms shown in Figs. 5(a1)–(a3) since there are carrier fringes with a spatial frequency of *ξ*� <sup>0</sup> = |*ξ*<sup>0</sup> − 2*ξN*/3| = |2.8*ξN*/3 − 2*ξN*/3| = 0.8*ξN*/3 due to large amount of tilt in the reference plane wave. However, it must be noted that there is the phase-shift relation between these coarse fringes with the phase value of 2*π*/3 each other. Therefore, after interpolating process we obtain three phase-shift fringes which are available to calculate a phase distribution where the phase of object wave and linear phase of *ξ*� <sup>0</sup> *nd* = (0.6*ξN*/3)*nd* (*n* = 1, 2, 3, . . .) are added. Final phase distribution of the object wave at the image sensor can be obtained by subtracting the spatial linear phase *ξ*� <sup>0</sup> *nd* = (0.6*ξN*/3)*nd* (*n* = 1, 2, 3, . . .). This result supports our derived

equations and procedures, *i.e.*, Eqs.(23)–(25) and compensation of the linear phase given by

<sup>101</sup> Single-Shot Phase-Shifting Digital Holography Based

Figures 5(c1)–(c3) show sub-holograms from the Fig. 5(c). In contrast to previous two cases, there are no phase-shifting relation between each sub-hologram. This result comes from the fact that there are no aliasing, in other word, no moiré fringe in the obtained hologram shown in Fig. 5(c). As predicted by our theoretical analysis we are not able not derived the phase map of the object wave at the image sensor by using the phase-shifting algorithm in this case. It will be done by using Fourier transform method Takeda et al. (1982) with the single-shot

We have investigated theoretically tolerance of the incident angle of plane reference wave in space-division multiplexed single-shot phase-shifting digital holography presented by Sato *et al* Toge et al. (2008). It is found from our analysis that the rigorous alignment of the incident angle of the plane reference wave is very important requirement to reconstruct the correct complex amplitude of the object wave by direct application of the phase-shifting algorithm. It is also found that even if the incident angle of the plane reference wave departs from the condition given by Eq. (1), we have phase-shifted sub holograms from the single digital hologram by sparsely sampling. However, the direct application of the phase-shifting algorithm for these phase-shifted holograms results in wrong complex amplitude especially in phase distribution. We described the way to obtain correct phase distribution in this case from our theoretical analysis. That is to subtract the linear phase distribution introduced by misalignment of the reference plane wave. This solution was first derived and brought to light by our theoretical analysis within the authors knowledge, and widen field of the single-shot

Awatsuji, Y.; Sasada, M.; Kubota, T. (2004). Parallel quasi-phase-shifting digital holography,

Gaskill, J. (1978). *Linear Systems, Fourier Transforms, and Optics*, John Wiley & Sons,

Harada, Y.; Goto, H.; Iwaki, T.; Wakaguri, S. (2008). Single-shot phase-shifting digital

Harada, Y.; Goto, H.; Wan, A. (2009). Spatially-multiplexed single-shot phase-shifting digital

Kreis, T. (2005). *Handbook of Holographic Interferometry*, Wiley VCH, ISBN:3-527-40546-1,

Nomura, T.; Murata, S.; Nitanai, E; Numata, T. (2006). Phase-shifting digital holography

Schnars, U. & Jueptner, W. (1978). *Digital Holography*, Springer, ISBN:3-540-21934-x, Berlin.

holography based on the spatial carrier interferometry, *Extended Abstract on Optics & Photonics Japan 2008*, Optical Society of Japan (Japan Society of Applied Physics),

holography based on the spatial carrier interferometry, *Holographic Display Artists and Engineers Club (HODIC) Circular*, Optical Society of Japan (Japan Society of Applied

with a phase difference between orthogonal polarization, *Applied Optics*, Vol.45, No.

phase-shifting digital holography using off-axis plane wave illuminations.

*Applied Physics Letters*, Vol.85, No.6, 1069-1071.

Tsukuba, Japan, pp.284-285. (in Japanese).

Physics), Vol.20, No. 29, 35-44. (in Japanese).

ISBN:0-471-29288-5, New York.

Eq. (26), for the case where the condition Eq. (1) is not satisfied.

on the Spatial Carrier Interferometry and Its Tolerance Analysis

hologram of Fig. 5(c).

**4. Conclusions**

**5. References**

Weinheim.

4873–4877.

Fig. 5. Typical results of the computer simulations: (a) interference fringes, *Is*(*x*, *y*), for sin *θ* = *λ*/3*d* and *ξ*<sup>0</sup> = 1/3*d* = 2*ξN*/3, (b) interference fringe, *Is*(*x*, *y*), for sin *θ* = 1.4*λ*/3*d* > *λ*/3*d* and *ξ*<sup>0</sup> = 1.4/3*d* = 2.8*ξN*/3, (c) interference fringe, *Is*(*x*, *y*), for sin *θ* = 0.8*λ*/6*d* < *λ*/6*d* and *ξ*<sup>0</sup> = 0.8/6*d* = 0.8*ξN*/3, (a1)–(a3) corresponding sub-holograms *I* (0) *<sup>s</sup>* , *<sup>I</sup>* (1) *<sup>s</sup>* , *<sup>I</sup>* (2) *<sup>s</sup>* for sin *<sup>θ</sup>* <sup>=</sup> *<sup>λ</sup>*/3*d*, (b1)–(b3) corresponding sub-holograms *<sup>I</sup>* (0) *<sup>s</sup>* , *<sup>I</sup>* (1) *<sup>s</sup>* , *<sup>I</sup>* (2) *<sup>s</sup>* for sin *θ* = 1.4*λ*/3*d* > *λ*/3*d*, and (c1)–(c3) corresponding sub-holograms *I* (0) *<sup>s</sup>* , *<sup>I</sup>* (1) *<sup>s</sup>* , *<sup>I</sup>* (2) *<sup>s</sup>* for sin *θ* = 0.8*λ*/6*d* < *λ*/6*d*.

equations and procedures, *i.e.*, Eqs.(23)–(25) and compensation of the linear phase given by Eq. (26), for the case where the condition Eq. (1) is not satisfied.

Figures 5(c1)–(c3) show sub-holograms from the Fig. 5(c). In contrast to previous two cases, there are no phase-shifting relation between each sub-hologram. This result comes from the fact that there are no aliasing, in other word, no moiré fringe in the obtained hologram shown in Fig. 5(c). As predicted by our theoretical analysis we are not able not derived the phase map of the object wave at the image sensor by using the phase-shifting algorithm in this case. It will be done by using Fourier transform method Takeda et al. (1982) with the single-shot hologram of Fig. 5(c).
