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

Multiwavelength digital holography with wavelength-multiplexed holograms and arbitrary symmetric phase shifts

Open Access Open Access

Abstract

We propose multiwavelength in-line digital holography with wavelength-multiplexed phase-shifted holograms and arbitrary symmetric phase shifts. We use phase-shifting interferometry selectively extracting wavelength information to reconstruct multiwavelength object waves separately from wavelength-multiplexed monochromatic images. The proposed technique obtains systems of equations for real and imaginary parts of multiwavelength object waves from the holograms by introducing arbitrary symmetric phase shifts. Then, the technique derives each complex amplitude distribution of each object wave selectively and analytically by solving the two systems of equations. We formulate the algorithm in the case of an arbitrary number of wavelengths and confirm its validity numerically and experimentally in the cases where the number of wavelengths is two and three.

© 2017 Optical Society of America

1. Introduction

Holography [1] is a technique used to record and reconstruct a wavefront diffracted from an object. Both amplitude and phase information of an object wave are recorded by utilizing interference of light, and a three-dimensional (3D) image is reconstructed by the obtained amplitude and phase distributions. The recorded information is called a “hologram”. Digital holography [2,3] is a technique used to record a digital hologram that contains an object wave, and to reconstruct both the 3D and quantitative phase images of an object using a computer. This technique has potential application to the fields of microscopy [4], particles and flow measurements [5], quantitative phase imaging [6], multimodal imaging [7], and encryption [8].

Phase-shifting digital holography [9,10] is a representative technique to capture an object wave with a full space-bandwidth product of an image sensor [11]. This digital holography has performed 3D shape measurement with multiwavelength phase unwrapping [12] and holographic 3D color imaging [13,14]. In phase-shifting color/multiwavelength digital holography, there are two types of representative implementations: time division [12] and space-division multiplexing [13,14] of multiple wavelengths. In the former, wavelength information is recorded sequentially by changing the wavelengths of light for measurement. Therefore mechanical shutters to select the measurement wavelength or operations for turning the lasers on and off are required. Furthermore, at least three phase-shifted holograms are required at each wavelength. As a result, this implementation is time consuming. In the latter, simultaneous recording of three wavelengths is possible by using a color image sensor with a Bayer color-filter array, and three phase-shifted holograms are sufficient. However, the crosstalk between object waves at multiple wavelengths occurs if the wavelength selectivity of the array is low [15]. Furthermore, recordable spatial bandwidth is restricted by the array. Space-division multiplexing decreases the field of view (FOV) and resolution, and the FOV becomes a quarter in square compared to single-wavelength phase-shifting technique. There is a tradeoff between the number of recordable wavelengths and spatial information. Therefore, it is difficult to conduct multispectral 3D imaging with a wide FOV.

Since 2013, we have presented phase-shifting digital holography selectively extracting wavelength information using wavelength-multiplexed phase-shifted holograms to obtain multiwavelength information without a color-filter array [16–19]. In this technique, multiple wavelengths are superimposed on each other in the space and spatial frequency domains and separated in the polar coordinate plane by using phase shifts, which is called phase-division multiplexing (PDM) of wavelengths. With the technique, a full space-bandwidth product can be available at each wavelength and there is no need for changing the wavelengths of the light sources. Where the number of wavelengths is N, 2N + 1 wavelength-multiplexed holograms are sufficient for 3D multiwavelength imaging, while 3N holograms are necessary in the time-division technique. Furthermore, we have proposed a technique that requires only 2N wavelength-multiplexed holograms by making the best use of 2π ambiguity of the phase [19]. Recordable wavelength bandwidth is determined by the spectral sensitivity of a monochromatic image sensor and is wider than an implementation with a color image sensor. However, up to now, phase shifts with integral multiples of 2π were needed to separately extract object waves with multiple wavelengths, and experimental demonstrations for only two wavelengths were reported in the digital holography. It is important to present a concrete image-reconstruction algorithm applicable to multicolor and multispectral 3D measurements.

In this paper, we propose multiwavelength digital holography with wavelength-multiplexed holograms and arbitrary symmetric phase shifts. The proposed technique is based on PDM of wavelengths, and arbitrary phase shifts can be adopted. We present the image reconstruction algorithm applicable when the measurement wavelengths are more than two. The proposed algorithm makes the best use of wavelength separation in the polar coordinate plane. By the technique introducing arbitrary symmetric phase shifts, multiple object waves are analytically derived from the systems of equations for real and imaginary parts of the waves. Formulation in the case of an arbitrary number of wavelengths helps to make possible multispectral holographic imaging without approximation. Our motivation of the proposal is to achieve multispectral holographic imaging analytically with small number of holograms, small phase shifts which are given by a piezo-driven mirror, no approximation, and full space-bandwidth product of a monochromatic image sensor at each wavelength, which are impossible for the conventional PDM techniques. We demonstrate two- and three-wavelength digital holographic imaging ability of the proposed technique.

2. Phase-shifting interferometry selectively extracting wavelength information

The concept of the proposed digital holography is based on phase-shifting interferometry selectively extracting wavelength information [16–19]. Figure 1 illustrates the schematic of the basic concept. Multiple wavelength information is superimposed in the space and spatial frequency domains, and wavelength-multiplexed phase-shifted holograms are continuously recorded by introducing different phase shifts per wavelength as shown in Fig. 1(a). Wavelength information is separated in the polar coordinate plane as illustrated in Fig. 1(b) and this separation is derived by introducing different phase shifts for respective wavelengths. Object waves with multiple wavelengths are obtained separately by phase-shifting interferometry selectively extracting wavelength information and a multiwavelength 3D image is reconstructed with diffraction integrals. It is noted that, where the number of wavelengths is N, 2N + 1 holograms are required for the initially presented PDM [16–18].

 figure: Fig. 1

Fig. 1 Basic concept of phase-shifting interferometry selectively extracting wavelength information. (a) Flow from recording to reconstruction. (b) Wavelength separation in the polar coordinate plane when applying arbitrary symmetric phase shifts.

Download Full Size | PDF

However, up to now, utilization of 2π ambiguity was needed to separately extract object waves with multiple wavelengths. As a result, large phase shifts were needed and therefore multicolor 3D imaging was not demonstrated experimentally because it was difficult to prepare a phase shifter that satisfies the need. In the case of recording holograms with three wavelengths of633, 532, and 473 nm by the conventional phase-division technique, a piezo actuator must move a mirror at 336.756 μm distance in total with one-nanometer positioning accuracy. It is difficult to prepare such a mirror with a piezo actuator to satisfy the required specification due to the effect of hysteresis, 105 required dynamic range (from 10−4 to 10−9 meter) with 1 nm-order resolution, and high cost of the piezo-driven mirror. Moreover, there is a poor applicability to incoherent and low-coherence digital holography conventionally because the required moving distance becomes longer than the coherence length. Although the total movement distance is decreased by choosing the wavelengths carefully, in the case of three wavelengths of 640, 532, and 488 nm, whole movement distance is 85.120 μm with one-nanometer accuracy. Furthermore, regular phase shifts were not possible by the previously reported algorithms.

3. Multiwavelength digital holography with wavelength-multiplexed holograms and arbitrary symmetric phase shifts

Figure 2 illustrates an optical implementation of the proposed multiwavelength digital holography in the case where the number of wavelengths is N. A phase modulator such as a mirror with a piezo actuator, wave plates, or a spatial light modulator sequentially generates different phase shifts per wavelength, and a monochromatic image sensor records 2N + 1 wavelength-multiplexed phase-shifted holograms I(x,y:α11,α21, ..., αN1), ..., I(x,y:α1 2N+1,α2 2N+1, ..., αN 2N+1), where α1, α2, ..., and αN are phase shifts at the wavelengths of λ1, λ2, ..., and λN. Where I0th(x,y) is the summation of 0th-order diffraction waves, Ao(x,y) and Ar(x,y) are amplitude distributions of object and reference waves, ϕo(x,y) is a phase distribution of an object wave, then a recorded wavelength-multiplexed hologram is expressed as follows:

 figure: Fig. 2

Fig. 2 An optical implementation and obtained holograms in the proposed digital holography.

Download Full Size | PDF

I(x,y:α11,α21,,αN1)= I0th(x,y)+2i=1NAoi(x,y)Ari(x,y)cos[ϕoi(x,y)αi1].

By using a matrix equation, 2N + 1 holograms I(x,y:α11,α21, ..., αN1), ..., I(x,y:α1 2N+1,α2 2N+1, ..., αN 2N+1) are expressed as

(I(x,y:0,,0)I(x,y:α11,,αN1)I(x,y:α12,,αN2)I(x,y:α12N,,αN2N))=(1101cosα11sinα111cosα12sinα12...10...cosαN1sinαN1...cosαN2sinαN21cosα12Nsinα12N......cosαN2NsinαN2N)×(I0th(x,y)2Ao1(x,y)Ar1(x,y)cosϕo1(x,y)2Ao1(x,y)Ar1(x,y)sinϕo1(x,y)2AoN(x,y)ArN(x,y)sinϕoN(x,y)).
The proposed technique adopts arbitrary symmetric phase shifts and does not need integral multiples of 2π. By introducing symmetric phase shifts to the holograms I(x,y:0, 0, ..., 0), I(x,y:α11,α21, ..., αN1), I(x,y:-α11,-α21, ..., -αN1), ..., I(x,y:α1N,α2N, ..., αNN), and I(x,y:-α1N,-α2N, ..., -αNN), object waves are selectively extracted with simple mathematical expressions. Equation (2) is rewritten as
(I1(x,y)I2(x,y)I3(x,y)I2N(x,y)I2N+1(x,y))=(I(x,y:0,,0)I(x,y:α11,,αN1)I(x,y:α11,,αN1)I(x,y:α1N,,αNN)I(x,y:α1N,,αNN))=(1101cosα11sinα111cosα11sinα11...10...cosαN1sinαN1...cosαN1sinαN11cosα1Nsinα1N......cosαNNsinαNN1cosα1Nsinα1N...cosαNNsinαNN)×(I0th(x,y)2Ao1(x,y)Ar1(x,y)cosϕo1(x,y)2Ao1(x,y)Ar1(x,y)sinϕo1(x,y)2AoN(x,y)ArN(x,y)cosϕoN(x,y)2AoN(x,y)ArN(x,y)sinϕoN(x,y)).
Using the wavelength-multiplexed holograms obtained by the proposed technique, object waves at multiple wavelengths are selectively extracted with the following procedures. Firstly, we remove I0th(x,y) from the holograms with the equations below,
Rek(x,y)=i=1NAoi(x,y)Ari(x,y)cosϕoi(x,y)(1cosαik)(k=1,...,N)=2I1(x,y)[I2k(x,y)+I2k+1(x,y)]4,
Imk(x,y)=i=1NAoi(x,y)Ari(x,y)sinϕoi(x,y)sinαik=I2k(x,y)I2k+1(x,y)4.
As described in Eqs. (4) and (5), real and imaginary parts of N object waves are clearly separated. This is derived by symmetric phase shifts and arbitrary shifts are applicable. When Ar(x,y) is assumed as constant and phase shifts are known, N variables are contained in Eqs. (4) and (5), respectively. N systems of equations are obtained for real part Re1(x,y), Re2(x,y), ..., and ReN(x,y), and imaginary one Im1(x,y), Im2(x,y), ..., and ImN(x,y), respectively, and therefore N complex amplitude distributions of object waves U1(x,y) = Ao1(x,y)cosϕo1(x,y) + jAo1(x,y)sinϕo1(x,y), ..., and UN(x,y) are selectively, analytically, and rigorously derived, where j = (−1)1/2. LU decomposition or other mathematical techniques can solve 2N systems of equations analytically with no approximations. Thus, by the proposal of PDM technique with regular, small, and arbitrary symmetric phase shifts, it is easy to conduct multispectral digital holographic imaging using wavelength-multiplexed phase-shifted holograms. As examples, we describe applications to two- and three-wavelength digital holography.

3-1. Two-wavelength digital holography

Five wavelength-multiplexed phase-shifted holograms I1(x,y), I2(x,y), I3(x,y), I4(x,y), and I5(x,y) are recorded with phase shifts (0,0), (α11,α21), (-α11,-α21), (α12,α22), and (-α12,-α22) at the wavelengths of (λ1,λ2). From Eqs. (3)-(5) we obtain the formulas below.

Re1(x,y)=Ao1(x,y)Ar1(x,y)cosϕo1(x,y)(1cosα11)+Ao2(x,y)Ar2(x,y)cosϕo2(x,y)(1cosα21),
Re2(x,y)=Ao1(x,y)Ar1(x,y)cosϕo1(x,y)(1cosα12)+Ao2(x,y)Ar2(x,y)cosϕo2(x,y)(1cosα22),
Im1(x,y)=Ao1(x,y)Ar1(x,y)sinϕo1(x,y)sinα11+Ao2(x,y)Ar2(x,y)sinϕo2(x,y)sinα21,
Im2(x,y)=Ao1(x,y)Ar1(x,y)sinϕo1(x,y)sinα12+Ao2(x,y)Ar2(x,y)sinϕo2(x,y)sinα22.
Real and imaginary parts of two object waves are derived by solving two systems of equations. As a result, object waves at two wavelengths U1(x,y) and U2(x,y) are selectively extracted with the expressions:

U1(x,y)=Re1(x,y)(1cosα22)Re2(x,y)(1cosα21)Ar1(x,y)[(1cosα11)(1cosα22)(1cosα12)(1cosα21)]+jIm1(x,y)sinα22Im2(x,y)sinα21Ar1(x,y)[sinα11sinα22sinα12sinα21],
U2(x,y)=Re1(x,y)(1cosα12)Re2(x,y)(1cosα11)Ar2(x,y)[(1cosα12)(1cosα21)(1cosα11)(1cosα22)]+jIm1(x,y)sinα12Im2(x,y)sinα11Ar2(x,y)[sinα12sinα21sinα11sinα22].

3-2. Three-wavelength digital holography

From seven wavelength-multiplexed phase-shifted holograms I1(x,y), I2(x,y), I3(x,y), I4(x,y), I5(x,y), I6(x,y), and I7(x,y) and Eqs. (3)-(5), we obtain object waves at the wavelengths of λ1, λ2, and λ3.

Re1(x,y)=Ao1(x,y)Ar1(x,y)cosϕo1(x,y)(1cosα11)+Ao2(x,y)Ar2(x,y)cosϕo2(x,y)(1cosα21)+Ao3(x,y)Ar3(x,y)cosϕo3(x,y)(1cosα31),
Re2(x,y)=Ao1(x,y)Ar1(x,y)cosϕo1(x,y)(1cosα12)+Ao2(x,y)Ar2(x,y)cosϕo2(x,y)(1cosα22)+Ao3(x,y)Ar3(x,y)cosϕo3(x,y)(1cosα32),
Re3(x,y)=Ao1(x,y)Ar1(x,y)cosϕo1(x,y)(1cosα13)+Ao2(x,y)Ar2(x,y)cosϕo2(x,y)(1cosα23)+Ao3(x,y)Ar3(x,y)cosϕo3(x,y)(1cosα33),
Im1(x,y)=Ao1(x,y)Ar1(x,y)sinϕo1(x,y)sinα11+Ao2(x,y)Ar2(x,y)sinϕo2(x,y)sinα21+Ao3(x,y)Ar3(x,y)sinϕo3(x,y)sinα31,
Im2(x,y)=Ao1(x,y)Ar1(x,y)sinϕo1(x,y)sinα12+Ao2(x,y)Ar2(x,y)sinϕo2(x,y)sinα22+Ao3(x,y)Ar3(x,y)sinϕo3(x,y)sinα32,
Im3(x,y)=Ao1(x,y)Ar1(x,y)sinϕo1(x,y)sinα13+Ao2(x,y)Ar2(x,y)sinϕo2(x,y)sinα23+Ao3(x,y)Ar3(x,y)sinϕo3(x,y)sinα33.
Three unknown variables are contained in each equation described above and therefore six unknown variables in two systems of equations are solved analytically and rigorously from seven wavelength-multiplexed holograms as follows:
Ao1(x,y)cosϕo1(x,y)=C4Re'(x,y)C3Re''(x,y)Ar1(x,y)(C1C4C2C3),
Ao2(x,y)cosϕo2(x,y)=(1cosα32)Re1(x,y)(1cosα31)Re2(x,y)Ar2(x,y)C3C1Ao1(x,y)Ar1(x,y)cosϕo1(x,y)Ar2(x,y)C3,
Ao3(x,y)cosϕo3(x,y)=Re3(x,y)(1cosα13)Ao1(x,y)Ar1(x,y)cosϕo1(x,y)Ar3(x,y)(1cosα33)(1cosα23)Ao2(x,y)Ar2(x,y)cosϕo2(x,y)Ar3(x,y)(1cosα33),
Ao1(x,y)sinϕo1(x,y)=S4Im'(x,y)S3Im''(x,y)Ar1(x,y)(S1S4S2S3),
Ao2(x,y)sinϕo2(x,y)=sinα32Im1(x,y)sinα31Im2(x,y)Ar2(x,y)S3S1Ao1(x,y)Ar1(x,y)sinϕo1(x,y)Ar2(x,y)S3,
Ao3(x,y)sinϕo3(x,y)=Im3(x,y)sinα13Ao1(x,y)Ar1(x,y)sinϕo1(x,y)Ar3(x,y)sinα33sinα23Ao2(x,y)Ar2(x,y)sinϕo2(x,y)Ar3(x,y)sinα33,
where,
Re'(x,y)=(1cosα32)Re1(x,y)(1cosα31)Re2(x,y) ,
Re''(x,y)=(1cosα33)Re1(x,y)(1cosα31)Re3(x,y),
C1=(1cosα11)(1cosα32)(1cosα12)(1cosα31),
C2=(1cosα11)(1cosα33)(1cosα13)(1cosα31),
C3=(1cosα21)(1cosα32)(1cosα22)(1cosα31),
C4=(1cosα21)(1cosα33)(1cosα23)(1cosα31),
Im'(x,y)=sinα32Im1(x,y)sinα31Im2(x,y),
Im''(x,y)=sinα33Im1(x,y)sinα31Im3(x,y),
S1=sinα11sinα32sinα12sinα31,
S2=sinα11sinα33sinα13sinα31,
S3=sinα21sinα32sinα22sinα31,
S4=sinα21sinα33sinα23sinα31.
In the same manner, N object waves are derived from 2N + 1 wavelength-multiplexed phase-shifted holograms by using Eqs. (3)-(5) and solving systems of equations. In comparison with time-division technique, we can obtain a color 3D image more quickly with the PDM technique. The acceleration of the measurement speed is due to the following reasons. Phase-shifting color digital holography with time-division technique requires 3N holograms. Furthermore, N mechanical shutters to select a wavelength from N wavelengths should be set between lasers and beam combiners. A computer to synchronize the shutters and exposure of an image sensor is required to obtain holograms without crosstalk between wavelengths. In contrast, the proposed technique can obtain multiwavelength object waves separately with 2N + 1 holograms. Therefore, measurement speed is 3N/(2N + 1) times improved by the proposed technique in principle. The proposed technique can obtain four-wavelength object waves during obtaining three wavelengths by the time-division technique theoretically. As increasing the number of wavelengths recorded, the degree of acceleration approaches to 1.5. Moreover, no mechanical shutter is needed to separate wavelength information and measurement speed is not affected by the synchronization and shutters. Therefore, the proposed technique is nearly 1.5 times faster than the time-division technique theoretically, and can further accelerate the measurement speed practically.

4. Numerical simulations

We conducted numerical simulations to verify the effectiveness of the proposed technique in the case of N = 3. We assumed an optical system containing a mirror with a piezo actuator used as a phase modulator. Figure 3 shows the amplitude and phase distributions and color-synthesized image of the object wave. λ1 = 633 nm, λ2 = 532 nm, and λ3 = 473 nm were assumed as the red-, green-, and blue-wavelengths of the light sources. Respective color components of a standard image “peppers” were used as amplitude images at the wavelengths, respectively. In these simulations, distance between the object and image sensor was assumed as 250 mm, pixel pitch was 2.2 μm, number of pixels was 512 × 512, and dynamic range of the image sensor was 12 bits. We set the condition that piezo-driven mirror moved at regular distances sequentially and obtained wavelength-multiplexed phase-shifted holograms sequentially. We obtained wavelength-multiplexed holograms with various optical-path shifts of the piezo-driven mirror numerically and reconstructed the object images with the proposed technique. Then we calculated root-mean-square errors (RMSEs) and cross-correlation coefficients (CCs) of the reconstructed amplitude images to evaluate the image quality quantitatively.

 figure: Fig. 3

Fig. 3 The optical system and object wave set for numerical simulations, and numerically obtained hologram. (a) Optical setup with a piezo-driven mirror. (b) Amplitude and (c) phase distributions of the object wave. (d) Red-, (e) green-, and (f) blue-color components of (b). (g) One of the numerically obtained wavelength-multiplexed phase-shifted holograms.

Download Full Size | PDF

Figure 4 shows the numerical results. Numerical results clarified that moving a distance of 50 nm of the mirror with a piezo actuator per phase shift was sufficient for faithful holographic multiwavelength 3D imaging. From the color-synthesized image seen in Fig. 4(d), multicolor 3D imaging ability was verified numerically. However, random noise can be seen slightly on the phase image in Fig. 4(f) because minor noise was not removed due to the small phase shifts. Such noise was successfully removed by introducing sufficient optical path shifts to separate wavelengths, as shown in Figs. 4(h)-4(j). Figure 5 indicates the graphs of the calculated RMSEs and CCs. These graphs showed its color 3D imaging ability quantitatively with less than a 100 nm shift of the mirror with a piezo actuator for each phase shift. Numerical results indicate that whole movement distance to achieve three-wavelength digital holographic imaging with high accuracy is decreased to 300 nm when recording three wavelengths of 633, 532, and 473 nm. The whole distance is further shortened to 240 nm although a little noise appears in Fig. 5. From the results, it is concluded that the whole distance of the mirror can be decreased to 1/1123 by the proposed technique in comparison to the conventional PDM. Furthermore, the minimum regular phase shifts can be estimated from the numerical results. Figures 5(a) and 5(b) indicate small phase shifts can be applied because phase shifts are calculated from the wavelengths used and moving distance of the mirror with a piezo actuator. The minimum shifts to achieve high-quality three-wavelength digital holographic imaging were estimated as 56.9, 67.7, and 76.1 degrees at the wavelengths of 633 nm, 532 nm, and 473 nm, respectively, in the case where Fresnel domain was assumed and dynamic range of a monochromatic image sensor was 12 bits. This estimation is derived from the numerical results with 50 nm regular phase shifts, as shown in Figs. 4(a)-4(g), 5(a), and 5(b). As a result, the minimum difference between the phase shifts at the wavelengths of 633 nm and 532 nm can be estimated as 10.8 degrees and that between the shifts at 532 nm and 473 nm can be as 8.44 ones. Thus, the minimum difference of phase shifts is estimated as 8.44 degrees in the case described above, which is less than π/20. The proposed PDM technique requires the wavelength separations in the polar coordinate plane and therefore the phase shifts should be different. However, the phase-shift difference of π/20 enabled selective extractions of multiple object waves. Additionally, the numerical simulations clarify that the proposed technique has the potential of the applicability to incoherent and low-coherence color digital holography because a total moving distance required for the proposed technique can be less than the coherence length of incoherent and low-coherence light. Therefore, the problems described for the previous publications are solved by the proposal. It is worthy of note that RMSEs increases and CCs decreases when the optical-path shifts of the reference arm are close to integral multiples of the halves of the respective wavelengths, which correspond to regular π or 2π phase shifts. When adopting regular phase shifts as performed in the numerical simulations, regular π phase shifts at a wavelength prevent to extract the imaginary part of the object wave at the corresponding wavelength, in the same manner of single-wavelength phase-shifting digital holography. Equation (5) indicates that sinα at the corresponding wavelength becomes zero when α is an integral multiple of π. Therefore, there is no way to retrieve imaginary part of the object wave at the wavelength and the conjugate image superimposes on the object image. Moreover, when α is an integral multiple of 2π, the hologram intensity at the wavelength is not changed and the object wave at the corresponding wavelength is removed by the subtraction procedures based on Eqs. (4) and (5). As a result, object information cannot be obtained. Therefore, π or 2π phase shifts should be avoided for their respective wavelengths and integral multiple of π is not suitable for regular phase shifts in the proposed technique. Thus, the validity of the proposed technique was numerically confirmed, and the guideline for phase shifting was indicated.

 figure: Fig. 4

Fig. 4 Numerical results. Reconstructed amplitude images at the wavelengths of (a) 633 nm, (b) 532 nm, and (c) 473 nm and (d) the color-synthesized image. Phase images at (e) 633 nm, (f) 532 nm, and (g) 473 nm obtained from the holograms with a 50 nm shift of the mirror for each phase shift. Phase images at (h) 633 nm, (i) 532 nm, and (j) 473 nm obtained from the holograms with a 190 nm shift.

Download Full Size | PDF

 figure: Fig. 5

Fig. 5 (a) RMSEs and (b) CCs of the reconstructed amplitude images, which are calculated for quantitative evaluations of the image quality.

Download Full Size | PDF

5. Experiments

We have experimented to demonstrate and validate the proposed technique. We constructed a reflection-type two-wavelength digital holography system and a transmission-type three-wavelength digital holographic microscopy system to set the cases of N = 2 and 3. In these experiments, we adopted a mirror with a piezo actuator as a phase modulator.

5-1. Two-wavelength digital holography

Five wavelength-multiplexed phase-shifted holograms were recorded sequentially by changing the optical path of the reference arm with an in-line configuration. The oscillation wavelengths of the lasers used as light sources λ1 and λ2 were 671 nm (MSL-FN-671-100-1, CNI) and 532 nm (GLK-32300, LASOS), respectively. Coherence lengths of the lasers are longer than 10 meters. A monochromatic CMOS image sensor was used to record the grayscale holograms. The sensor has 12 bits, 2592 × 1944 pixels, and a pixel pitch of 2.2 μm. The mirror with a piezo actuator moved 70 nm for each phase shift sequentially. A miniature model of baby birds with green-colored grass was set as a color 3D object. The front side of the object was set at the distance of 115 mm from the image sensor plane. For comparison, object images were also reconstructed from a wavelength-multiplexed hologram obtained with the same in-line configuration to highlight the ability of the proposal for the wavelength separation.

Figure 6 shows the experimental results. In the same manner presented by experts in color digital holography with coherent light [20–23], we showed a photograph of the object taken with white light, monochromatic images at respective wavelengths, and its color-synthesized images in Fig. 6. As seen in Figs. 6(c) and 6(d), not only the 0th-order diffraction wave and the conjugate image, but also image components given by the crosstalk between the intensity distributions at λ1 and λ2 were reconstructed from a wavelength-multiplexed hologram. In contrast, object waves at the respective wavelengths were clearly reconstructed by the proposed technique as shown in Figs. 6(e) and 6(f). Moreover, depth information was obtained and focused images of the objects placed at different depths were successfully reconstructed as seen in Figs. 6(e)-6(h). As a result, color images at arbitrary depths were obtained, as shown in Figs. 6(i) and 6(j). Figure 6 visualizes obviously the color difference between the miniature models of the baby birds and grass, which is seen under the birds and has different color from the birds. Green-color grass was reconstructed as green color by the proposed technique with two wavelengths. Thus, the proposed technique has performed two-color 3D imaging of scattering objects. Thus, the effectiveness and color 3D imaging ability were experimentally demonstrated. Furthermore, these results clarified that the proposed technique was applicable to a scattering colored 3D object with rough surface.

 figure: Fig. 6

Fig. 6 Experimental results. (a) A photograph of the object that is illuminated by natural light. (b) One of wavelength-multiplexed phase-shifted holograms and its magnified image. Images reconstructed with Fig. 6(b) at the wavelengths of (c) 640 mm and (d) 532 mm. Images reconstructed by the proposed technique at the wavelengths of (e) 640 mm and (f) 532 mm, which are focused on the depth of 115 mm from the image sensor plane. Those at the wavelengths of (g) 640 mm and (h) 532 mm, and 145 mm distance from the image sensor plane. Color-synthesized images obtained from (i) Figs. 6(e) and 6(f), and (j) Figs. 6(g) and 6(h).

Download Full Size | PDF

5-2. Three-wavelength digital holographic microscopy

We constructed a three-wavelength in-line digital holographic microscope based on the proposed technique. Figure 7 illustrates a schematic of the system. A Telecentric magnification system was designed and its magnification and theoretical resolution at the wavelength of 488 nm were 25 and 774 nm, which were calculated from the ratio of focal lengths of the lenses in the telecentric system and the numerical aperture of the microscope objective N.A. = 0.8. CW lasers with the oscillation wavelengths λ1 = 640 nm (RLK-4075, LASOS), λ2 = 532 nm (GLK-32300, LASOS), and λ3 = 488 nm (Sapphire 488SF-100-HS, Coherent) were used as light sources. A monochromatic CMOS image sensor used was the same as the previously described experiment. A preparation of stained mouse kidney cells was set as a specimen. A piezo-driven mirror was moved a distance of 183 nm for each phase shift. We recorded one hologram per 1.5 seconds. Seven wavelength-multiplexed phase-shifted holograms with three wavelengths were obtained sequentially. Photographs of the specimen at respective wavelengths were taken with the lasers and CMOS image sensor sequentially to investigate the quality of the images reconstructed by the proposed technique. For comparison, we reconstructed three-wavelength images using a wavelength-multiplexed hologram obtained with the same in-line configuration and using calculation of diffraction integral alone. Figure 8 shows the experimental results. In the same manner of refs [20–23], we showed monochromatic images at the respective wavelengths and the color-synthesized images of the specimen illuminated by laser beams. As seen in Figs. 8(a)-8(l), wavelength information was successfully separated and similar images at the wavelengths and its color-synthesized one of the specimen were obtained by the proposed technique, while wavelength dependency of transparency cannot be visualized from a wavelength-multiplexed in-line hologram. Furthermore, 3D imaging of stained cell nuclei with a 1 μm-order radius were achieved as shown in Figs. 8(m) and 8(n). If the object images at the respective wavelengths were successfully retrieved from wavelength-multiplexed phase-shifted holograms, intensity difference between the photographs and reconstructed intensity images at these wavelengths on the image sensor plane approaches to zero. As described in ref [24], it is important to conduct quantitative evaluations of experimental results. Therefore, we calculated root-mean-square errors (RMSEs) to evaluate the reconstructed images quantitatively. Intensity images seen in Figs. 8(a)-8(c) were used for the evaluation as the original images in this experiment. Table 1 showed the results and errors were notably decreased by using the proposed technique in comparison to the method using a wavelength-multiplexed in-line hologram and diffraction-integral calculation alone. The results indicated the ability of the proposed technique for wavelength separation and removals of unwanted image components. Thus, the effectiveness of the proposed microscopy was experimentally demonstrated.

 figure: Fig. 7

Fig. 7 Schematic of the constructed three-wavelength digital holographic microscopy system.

Download Full Size | PDF

 figure: Fig. 8

Fig. 8 Experimental results. Photographs of a specimen at (a) λ1, (b) λ2, and (c) λ3, and (d) its color-synthesized image. (e) – (g) Intensity images reconstructed from a wavelength-multiplexed in-line hologram and (h) its color-synthesized image. (i) – (k) Images obtained by the proposed microscopy and (l) its color-synthesized image. Magnified color-synthesized images in which focused planes are (m) 0 and (n) 10 mm from the image sensor one, respectively. Arrows shown in (n) indicate the focused images of the stained cell nuclei on a preparation of a mouse kidney cells. Rectangles shown in (i) – (l) correspond with the areas of (m) and (n).

Download Full Size | PDF

Tables Icon

Table 1. RMSEs of the reconstructed images.

6. Conclusion

We have proposed multiwavelength in-line phase-shifting digital holography that uses wavelength-multiplexed phase-shifted holograms and arbitrary symmetric phase shifts, based on phase-shifting interferometry selectively extracting wavelength information. Its effectiveness was numerically confirmed and experimentally demonstrated. Numerical results clarified that it is possible to achieve noiseless holographic multiwavelength 3D imaging with less than 100 nm movement of the mirror per phase shift while it has been difficult with such small movements by the previously reported PDM techniques [25]. Experimental results clarified its ability of color 3D imaging of scattering objects with rough surfaces and its applicability to multicolor 3D microscopy. In comparison to the time-division technique, the proposed technique is valid for implementing a simple setup and expediting the measurement because the number of images recorded decreases.

The next step of the proposed PDM technique is the extension to multispectral holographic 3D image sensing, simultaneous imaging of color and 3D shape with multiwavelength phase unwrapping, dispersion imaging of a 3D specimen by using phase images at multiple wavelengths, multicolor incoherent digital holographic imaging, and multidimensional holographic imaging. This technique has prospective applications to multispectral microscopy to observe 3D specimens with a wide field of view, color 3D image sensing, multidimensional holographic image sensors, and other multiwavelength 3D imaging applications.

Funding

PRESTO, Japan Science and Technology Agency (JST) (JPMJPR16P8), The Nakajima Foundation (2017-2018), Konica Minolta Science and Technology Foundation (2013-2014), Research Foundation for Opt-Science and Technology (2013-2015), the Japan Society for the Promotion of Science (JSPS) (Grant-in-Aid for Young Scientists (B) 15K17474), and MEXT-Supported Program for the Strategic Research Foundation at Private Universities (2013-2018).

Acknowledgments

We appreciate Ms. Kris Cutsail for checking the grammar of this article. We thank Dr. Sumio Nakahara and Dr. Shigeyoshi Hisada to give us optical elements and anti-vibration tables.

References and links

1. D. Gabor, “A new microscopic principle,” Nature 161(4098), 777–778 (1948). [CrossRef]   [PubMed]  

2. J. W. Goodman and R. W. Lawrence, “Digital image formation from electronically detected holograms,” Appl. Phys. Lett. 11(3), 77–79 (1967). [CrossRef]  

3. T.-C. Poon and J.-P. Liu, eds., Introduction to Modern Digital Holography with MATLAB (Cambridge University, 2014).

4. Y. Takaki, H. Kawai, and H. Ohzu, “Hybrid holographic microscopy free of conjugate and zero-order images,” Appl. Opt. 38(23), 4990–4996 (1999). [CrossRef]   [PubMed]  

5. S. Murata and N. Yasuda, “Potential of digital holography in particle measurement,” Opt. Laser Technol. 32(7-8), 567–574 (2000). [CrossRef]  

6. E. Watanabe, T. Hoshiba, and B. Javidi, “High-precision microscopic phase imaging without phase unwrapping for cancer cell identification,” Opt. Lett. 38(8), 1319–1321 (2013). [CrossRef]   [PubMed]  

7. X. Quan, K. Nitta, O. Matoba, P. Xia, and Y. Awatsuji, “Phase and fluorescence imaging by combination of digital holographic microscopy and fluorescence microscopy,” Opt. Rev. 22(2), 349–353 (2015). [CrossRef]  

8. W. N. Li, C.-X. Shi, M.-L. Piao, and N. Kim, “Multiple-3D-object secure information system based on phase shifting method and single interference,” Appl. Opt. 55(15), 4052–4059 (2016). [CrossRef]   [PubMed]  

9. J. H. Bruning, D. R. Herriott, J. E. Gallagher, D. P. Rosenfeld, A. D. White, and D. J. Brangaccio, “Digital wavefront measuring interferometer for testing optical surfaces and lenses,” Appl. Opt. 13(11), 2693–2703 (1974). [CrossRef]   [PubMed]  

10. I. Yamaguchi and T. Zhang, “Phase-shifting digital holography,” Opt. Lett. 22(16), 1268–1270 (1997). [CrossRef]   [PubMed]  

11. A. Stern and B. Javidi, “Improved-resolution digital holography using the generalized sampling theorem for locally band-limited fields,” J. Opt. Soc. Am. A 23(5), 1227–1235 (2006). [CrossRef]   [PubMed]  

12. Y.-Y. Cheng and J. C. Wyant, “Two-wavelength phase shifting interferometry,” Appl. Opt. 23(24), 4539–4543 (1984). [CrossRef]   [PubMed]  

13. I. Yamaguchi, T. Matsumura, and J. Kato, “Phase-shifting color digital holography,” Opt. Lett. 27(13), 1108–1110 (2002). [CrossRef]   [PubMed]  

14. J. Kato, I. Yamaguchi, and T. Matsumura, “Multicolor digital holography with an achromatic phase shifter,” Opt. Lett. 27(16), 1403–1405 (2002). [CrossRef]   [PubMed]  

15. J. M. Desse, P. Picart, and P. Tankam, “Sensor influence in digital 3λ holographic interferometry,” Meas. Sci. Technol. 22(6), 064005 (2011). [CrossRef]  

16. T. Tahara, S. Kikunaga, Y. Arai, and Y. Takaki, “Phase-shifting interferometry capable of selectively extracting multiple wavelength information and color three-dimensional imaging using a monochromatic image sensor,” in Proceedings of Optics and Photonics Japan (OPJ, 2013), paper 13aE9.

17. T. Tahara, S. Kikunaga, Y. Arai, and Y. Takaki, “Phase-shifting interferometry capable of selectively extracting multiple wavelength information and its applications to sequential and parallel phase-shifting digital holography,” in Digital Holography and Three-Dimensional Imaging 2014 (DH), OSA Technical Digest (online) (Optical Society of America, 2014), paper DM3B.4.

18. T. Tahara, R. Mori, S. Kikunaga, Y. Arai, and Y. Takaki, “Dual-wavelength phase-shifting digital holography selectively extracting wavelength information from wavelength-multiplexed holograms,” Opt. Lett. 40(12), 2810–2813 (2015). [CrossRef]   [PubMed]  

19. T. Tahara, R. Mori, Y. Arai, and Y. Takaki, “Four-step phase-shifting digital holography simultaneously sensing dual-wavelength information using a monochromatic image sensor,” J. Opt. 17, 125707 (2015).

20. P. Ferraro, S. Grilli, L. Miccio, D. Alfieri, S. De Nicola, A. Finizio, and B. Javidi, “Full color 3-D imaging by digital holography and removal of chromatic aberrations,” J. Disp. Technol. 4(1), 97–100 (2008). [CrossRef]  

21. M. A. Araiza-Esquivel, L. Martínez-León, B. Javidi, P. Andrés, J. Lancis, and E. Tajahuerce, “Single-shot color digital holography based on the fractional Talbot effect,” Appl. Opt. 50(7), B96–B101 (2011). [CrossRef]   [PubMed]  

22. T. Kiire, D. Barada, J. Sugisaka, Y. Hayasaki, and T. Yatagai, “Color digital holography using a single monochromatic imaging sensor,” Opt. Lett. 37(15), 3153–3155 (2012). [CrossRef]   [PubMed]  

23. M. Leclercq and P. Picart, “Method for chromatic error compensation in digital color holographic imaging,” Opt. Express 21(22), 26456–26467 (2013), https://www.osapublishing.org/oe/abstract.cfm?uri=oe-21-22-26456. [CrossRef]   [PubMed]  

24. A. Alfalou and C. Brosseau, “Optical image compression and encryption methods,” Adv. Opt. Photonics 1(3), 589–636 (2009). [CrossRef]  

25. T. Tahara, R. Otani, Y. Arai, and Y. Takaki, “Multiwavelength digital holography and phase-shifting interferometry selectively extracting wavelength information: phase-division multiplexing (PDM) of wavelengths,” in Holographic Materials and Optical Systems, I. Naydenova Ed. (InTech, 2017).

Cited By

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

Alert me when this article is cited.


Figures (8)

Fig. 1
Fig. 1 Basic concept of phase-shifting interferometry selectively extracting wavelength information. (a) Flow from recording to reconstruction. (b) Wavelength separation in the polar coordinate plane when applying arbitrary symmetric phase shifts.
Fig. 2
Fig. 2 An optical implementation and obtained holograms in the proposed digital holography.
Fig. 3
Fig. 3 The optical system and object wave set for numerical simulations, and numerically obtained hologram. (a) Optical setup with a piezo-driven mirror. (b) Amplitude and (c) phase distributions of the object wave. (d) Red-, (e) green-, and (f) blue-color components of (b). (g) One of the numerically obtained wavelength-multiplexed phase-shifted holograms.
Fig. 4
Fig. 4 Numerical results. Reconstructed amplitude images at the wavelengths of (a) 633 nm, (b) 532 nm, and (c) 473 nm and (d) the color-synthesized image. Phase images at (e) 633 nm, (f) 532 nm, and (g) 473 nm obtained from the holograms with a 50 nm shift of the mirror for each phase shift. Phase images at (h) 633 nm, (i) 532 nm, and (j) 473 nm obtained from the holograms with a 190 nm shift.
Fig. 5
Fig. 5 (a) RMSEs and (b) CCs of the reconstructed amplitude images, which are calculated for quantitative evaluations of the image quality.
Fig. 6
Fig. 6 Experimental results. (a) A photograph of the object that is illuminated by natural light. (b) One of wavelength-multiplexed phase-shifted holograms and its magnified image. Images reconstructed with Fig. 6(b) at the wavelengths of (c) 640 mm and (d) 532 mm. Images reconstructed by the proposed technique at the wavelengths of (e) 640 mm and (f) 532 mm, which are focused on the depth of 115 mm from the image sensor plane. Those at the wavelengths of (g) 640 mm and (h) 532 mm, and 145 mm distance from the image sensor plane. Color-synthesized images obtained from (i) Figs. 6(e) and 6(f), and (j) Figs. 6(g) and 6(h).
Fig. 7
Fig. 7 Schematic of the constructed three-wavelength digital holographic microscopy system.
Fig. 8
Fig. 8 Experimental results. Photographs of a specimen at (a) λ1, (b) λ2, and (c) λ3, and (d) its color-synthesized image. (e) – (g) Intensity images reconstructed from a wavelength-multiplexed in-line hologram and (h) its color-synthesized image. (i) – (k) Images obtained by the proposed microscopy and (l) its color-synthesized image. Magnified color-synthesized images in which focused planes are (m) 0 and (n) 10 mm from the image sensor one, respectively. Arrows shown in (n) indicate the focused images of the stained cell nuclei on a preparation of a mouse kidney cells. Rectangles shown in (i) – (l) correspond with the areas of (m) and (n).

Tables (1)

Tables Icon

Table 1 RMSEs of the reconstructed images.

Equations (35)

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

I ( x , y : α 11 , α 21 , , α N 1 ) =   I 0 t h ( x , y ) + 2 i = 1 N A o i ( x , y ) A r i ( x , y ) cos [ ϕ o i ( x , y ) α i 1 ] .
( I ( x , y : 0 , , 0 ) I ( x , y : α 11 , , α N 1 ) I ( x , y : α 12 , , α N 2 ) I ( x , y : α 12 N , , α N 2 N ) ) = ( 1 1 0 1 cos α 11 sin α 11 1 cos α 12 sin α 12 ... 1 0 ... cos α N 1 sin α N 1 ... cos α N 2 sin α N 2 1 cos α 12 N sin α 12 N ... ... cos α N 2 N sin α N 2 N ) × ( I 0 t h ( x , y ) 2 A o 1 ( x , y ) A r 1 ( x , y ) cos ϕ o 1 ( x , y ) 2 A o 1 ( x , y ) A r 1 ( x , y ) sin ϕ o 1 ( x , y ) 2 A o N ( x , y ) A r N ( x , y ) sin ϕ o N ( x , y ) ) .
( I 1 ( x , y ) I 2 ( x , y ) I 3 ( x , y ) I 2 N ( x , y ) I 2 N + 1 ( x , y ) ) = ( I ( x , y : 0 , , 0 ) I ( x , y : α 11 , , α N 1 ) I ( x , y : α 11 , , α N 1 ) I ( x , y : α 1 N , , α N N ) I ( x , y : α 1 N , , α N N ) ) = ( 1 1 0 1 cos α 11 sin α 11 1 cos α 11 sin α 11 ... 1 0 ... cos α N 1 sin α N 1 ... cos α N 1 sin α N 1 1 cos α 1 N sin α 1 N ... ... cos α N N sin α N N 1 cos α 1 N s i n α 1 N ... cos α N N sin α N N ) × ( I 0 t h ( x , y ) 2 A o 1 ( x , y ) A r 1 ( x , y ) cos ϕ o 1 ( x , y ) 2 A o 1 ( x , y ) A r 1 ( x , y ) sin ϕ o 1 ( x , y ) 2 A o N ( x , y ) A r N ( x , y ) cos ϕ o N ( x , y ) 2 A o N ( x , y ) A r N ( x , y ) sin ϕ o N ( x , y ) ) .
R e k ( x , y ) = i = 1 N A o i ( x , y ) A r i ( x , y ) cos ϕ o i ( x , y ) ( 1 cos α i k ) ( k = 1 , ... , N ) = 2 I 1 ( x , y ) [ I 2 k ( x , y ) + I 2 k + 1 ( x , y ) ] 4 ,
I m k ( x , y ) = i = 1 N A o i ( x , y ) A r i ( x , y ) sin ϕ o i ( x , y ) sin α i k = I 2 k ( x , y ) I 2 k + 1 ( x , y ) 4 .
R e 1 ( x , y ) = A o 1 ( x , y ) A r 1 ( x , y ) cos ϕ o 1 ( x , y ) ( 1 cos α 11 ) + A o 2 ( x , y ) A r 2 ( x , y ) cos ϕ o 2 ( x , y ) ( 1 cos α 21 ) ,
R e 2 ( x , y ) = A o 1 ( x , y ) A r 1 ( x , y ) cos ϕ o 1 ( x , y ) ( 1 cos α 12 ) + A o 2 ( x , y ) A r 2 ( x , y ) cos ϕ o 2 ( x , y ) ( 1 cos α 22 ) ,
I m 1 ( x , y ) = A o 1 ( x , y ) A r 1 ( x , y ) sin ϕ o 1 ( x , y ) sin α 11 + A o 2 ( x , y ) A r 2 ( x , y ) sin ϕ o 2 ( x , y ) sin α 21 ,
I m 2 ( x , y ) = A o 1 ( x , y ) A r 1 ( x , y ) sin ϕ o 1 ( x , y ) sin α 12 + A o 2 ( x , y ) A r 2 ( x , y ) sin ϕ o 2 ( x , y ) sin α 22 .
U 1 ( x , y ) = R e 1 ( x , y ) ( 1 cos α 22 ) R e 2 ( x , y ) ( 1 cos α 21 ) A r 1 ( x , y ) [ ( 1 cos α 11 ) ( 1 cos α 22 ) ( 1 cos α 12 ) ( 1 cos α 21 ) ] + j I m 1 ( x , y ) sin α 22 I m 2 ( x , y ) sin α 21 A r 1 ( x , y ) [ sin α 11 sin α 22 sin α 12 sin α 21 ] ,
U 2 ( x , y ) = R e 1 ( x , y ) ( 1 cos α 12 ) R e 2 ( x , y ) ( 1 cos α 11 ) A r 2 ( x , y ) [ ( 1 cos α 12 ) ( 1 cos α 21 ) ( 1 cos α 11 ) ( 1 cos α 22 ) ] + j I m 1 ( x , y ) sin α 12 I m 2 ( x , y ) sin α 11 A r 2 ( x , y ) [ sin α 12 sin α 21 sin α 11 sin α 22 ] .
R e 1 ( x , y ) = A o 1 ( x , y ) A r 1 ( x , y ) cos ϕ o 1 ( x , y ) ( 1 cos α 11 ) + A o 2 ( x , y ) A r 2 ( x , y ) cos ϕ o 2 ( x , y ) ( 1 cos α 21 ) + A o 3 ( x , y ) A r 3 ( x , y ) cos ϕ o 3 ( x , y ) ( 1 cos α 31 ) ,
R e 2 ( x , y ) = A o 1 ( x , y ) A r 1 ( x , y ) cos ϕ o 1 ( x , y ) ( 1 cos α 12 ) + A o 2 ( x , y ) A r 2 ( x , y ) cos ϕ o 2 ( x , y ) ( 1 cos α 22 ) + A o 3 ( x , y ) A r 3 ( x , y ) cos ϕ o 3 ( x , y ) ( 1 cos α 32 ) ,
R e 3 ( x , y ) = A o 1 ( x , y ) A r 1 ( x , y ) cos ϕ o 1 ( x , y ) ( 1 cos α 13 ) + A o 2 ( x , y ) A r 2 ( x , y ) cos ϕ o 2 ( x , y ) ( 1 cos α 23 ) + A o 3 ( x , y ) A r 3 ( x , y ) cos ϕ o 3 ( x , y ) ( 1 cos α 33 ) ,
I m 1 ( x , y ) = A o 1 ( x , y ) A r 1 ( x , y ) sin ϕ o 1 ( x , y ) sin α 11 + A o 2 ( x , y ) A r 2 ( x , y ) sin ϕ o 2 ( x , y ) sin α 21 + A o 3 ( x , y ) A r 3 ( x , y ) sin ϕ o 3 ( x , y ) sin α 31 ,
I m 2 ( x , y ) = A o 1 ( x , y ) A r 1 ( x , y ) sin ϕ o 1 ( x , y ) sin α 12 + A o 2 ( x , y ) A r 2 ( x , y ) sin ϕ o 2 ( x , y ) sin α 22 + A o 3 ( x , y ) A r 3 ( x , y ) sin ϕ o 3 ( x , y ) sin α 32 ,
I m 3 ( x , y ) = A o 1 ( x , y ) A r 1 ( x , y ) sin ϕ o 1 ( x , y ) sin α 13 + A o 2 ( x , y ) A r 2 ( x , y ) sin ϕ o 2 ( x , y ) sin α 23 + A o 3 ( x , y ) A r 3 ( x , y ) sin ϕ o 3 ( x , y ) sin α 33 .
A o 1 ( x , y ) cos ϕ o 1 ( x , y ) = C 4 R e ' ( x , y ) C 3 R e ' ' ( x , y ) A r 1 ( x , y ) ( C 1 C 4 C 2 C 3 ) ,
A o 2 ( x , y ) cos ϕ o 2 ( x , y ) = ( 1 cos α 32 ) R e 1 ( x , y ) ( 1 cos α 31 ) R e 2 ( x , y ) A r 2 ( x , y ) C 3 C 1 A o 1 ( x , y ) A r 1 ( x , y ) cos ϕ o 1 ( x , y ) A r 2 ( x , y ) C 3 ,
A o 3 ( x , y ) cos ϕ o 3 ( x , y ) = R e 3 ( x , y ) ( 1 cos α 13 ) A o 1 ( x , y ) A r 1 ( x , y ) cos ϕ o 1 ( x , y ) A r 3 ( x , y ) ( 1 cos α 33 ) ( 1 cos α 23 ) A o 2 ( x , y ) A r 2 ( x , y ) cos ϕ o 2 ( x , y ) A r 3 ( x , y ) ( 1 cos α 33 ) ,
A o 1 ( x , y ) sin ϕ o 1 ( x , y ) = S 4 I m ' ( x , y ) S 3 I m ' ' ( x , y ) A r 1 ( x , y ) ( S 1 S 4 S 2 S 3 ) ,
A o 2 ( x , y ) sin ϕ o 2 ( x , y ) = sin α 32 I m 1 ( x , y ) sin α 31 I m 2 ( x , y ) A r 2 ( x , y ) S 3 S 1 A o 1 ( x , y ) A r 1 ( x , y ) sin ϕ o 1 ( x , y ) A r 2 ( x , y ) S 3 ,
A o 3 ( x , y ) sin ϕ o 3 ( x , y ) = I m 3 ( x , y ) sin α 13 A o 1 ( x , y ) A r 1 ( x , y ) sin ϕ o 1 ( x , y ) A r 3 ( x , y ) sin α 33 sin α 23 A o 2 ( x , y ) A r 2 ( x , y ) sin ϕ o 2 ( x , y ) A r 3 ( x , y ) sin α 33 ,
R e ' ( x , y ) = ( 1 cos α 32 ) R e 1 ( x , y ) ( 1 cos α 31 ) R e 2 ( x , y )   ,
R e ' ' ( x , y ) = ( 1 cos α 33 ) R e 1 ( x , y ) ( 1 cos α 31 ) R e 3 ( x , y ) ,
C 1 = ( 1 cos α 11 ) ( 1 cos α 32 ) ( 1 cos α 12 ) ( 1 cos α 31 ) ,
C 2 = ( 1 cos α 11 ) ( 1 cos α 33 ) ( 1 cos α 13 ) ( 1 cos α 31 ) ,
C 3 = ( 1 cos α 21 ) ( 1 cos α 32 ) ( 1 cos α 22 ) ( 1 cos α 31 ) ,
C 4 = ( 1 cos α 21 ) ( 1 cos α 33 ) ( 1 cos α 23 ) ( 1 cos α 31 ) ,
I m ' ( x , y ) = sin α 32 I m 1 ( x , y ) sin α 31 I m 2 ( x , y ) ,
I m ' ' ( x , y ) = sin α 33 I m 1 ( x , y ) sin α 31 I m 3 ( x , y ) ,
S 1 = sin α 11 sin α 32 sin α 12 sin α 31 ,
S 2 = sin α 11 sin α 33 sin α 13 sin α 31 ,
S 3 = sin α 21 sin α 32 sin α 22 sin α 31 ,
S 4 = sin α 21 sin α 33 sin α 23 sin α 31 .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.