# Blind Crosstalk Correction for Color Fringe Projection Profilometry Combining the Sine and Cosine Phase-Shifting Patterns

Authors: Jun Yang; Guoliang Li; Yang Zhang; Jun Wu; Min Zhang; Jin Li; Wenxiang Liu
Journal: Journal of Sensors
 (2022)

---

## Abstract

Color crosstalk is challenging for color fringe projection profilometry, which will cause color imbalance and phase-shift deviation, thus resulting in phase error. To address this challenge, this paper presents a blind crosstalk correction method for color phase-shifting fringe projection profilometry without crosstalk calibration. Three channels of the color fringe are firstly normalized to resolve color imbalance. The phase error induced by phase-shift deviation between three channels doubles the frequency of the color fringe. Therefore, the original cosine color fringe is shifted by quarter period to generate another sine color fringe to compensate the periodic phase error. Both simulations and experiments were carried out, and their results confirmed the effectiveness and accuracy of the proposed method.


---

## Body

## 1. Introduction

Fringe projection profilometry (FPP) based on triangulation principle is one of the most commonly used methods in optical three-dimensional (3D) shape measurement due to its low cost, fast speed, and high-resolution capabilities [1–6]. Normally, a FPP system contains a computer, a projector, and a camera. The projector and the camera are directed along different directions to the measured objects, all of which are controlled by the computer. Firstly, the fringe patterns should be designed and generated by the computer. Then, these patterns are projected onto the measured objects by the projector; at the same time, the deformed patterns are recorded by the camera. Finally, these captured images are transferred to the computer and analyzed by phase extraction algorithms to obtain the geometric information of the measured objects [7]. There are several phase extraction algorithms, including Fourier transform [8], wavelet transform [9], and phase-shifting profilometry (PSP) [10]. Among them, PSP is the most popular method, because it can achieve pixel-wise phase extraction, which can reconstruct the measured objects with abrupt changes more robustly. However, it generally requires more fringe patterns (normally at least three) than transform methods [11].In order to overcome this problem, color phase-shifting technique has been developed [12–19]. Three phase-shifting patterns are assigned to red, green, and blue channels, respectively. Thus, only one color fringe pattern is required to reconstruct 3D shapes. However, the spectra of red, green, and blue channels are usually made to have some overlaps so that there would be no color-blind areas on the spectrum [20–22]. Thence, phase error will be introduced if conventional phase-shifting algorithms are employed directly. To suppress the phase error caused by color crosstalk, Huang et al. [23] proposed a method to compensate the coupling effects between channels. Since the degrees of overlaps are fixed, they projected pure red, green, and blue fringe patterns onto a white-painted flat board to calculate the degrees of coupling effects between three channels and then compensated it. This method has been explored when different numbers of fringes are projected on each channel by Zhang [24]. On the other hand, Wan et al. [25] presented an active phase error suppression method based on hue precorrection. The hue is precorrected to be linear by measuring the function of the combined effect and modifying the computer-generated hue according to the source-distorted relationship. Then, the new modified color fringe pattern can be generated for measuring objects. However, all these techniques need preprocessing before measurement, which makes it inconvenient to use. Hu et al. [26] proposed an algorithm for estimating the color remixing matrix based on the color fringe patterns captured from the reference plane or the surface of objects. Ma et al. [27] proposed a blind phase error suppression approach based on isotropic n-dimensional fringe pattern normalization [28] and carrier squeezing interferometry [29]. Wu et al. [30] proposed a point-to-point compensation method for color coupling and imbalance. The proposed method achieves point-to-point compensation based on the correspondence between camera pixels and projector pixels. Recently, Wang et al. [31] found the phase error introduced by crosstalk is approximately twice the frequency of color fringe. Therefore, Hilbert transform was applied on the color fringe to compensate the phase error. However, due to the fact that Hilbert transform belongs to frequency domain methods, this approach cannot handle the measured objects that contain discontinuities or separations well.In this paper, we present a blind crosstalk correction method for color phase-shifting fringe projection profilometry. Due to the fact that phase error introduced by color crosstalk is twice the frequency of color fringe, the other color fringe with phase-shift is used to compensate the phase error. In addition, the fringe normalization was also designed to resolve the color imbalance problem. We describe the principle of this method in Section2. Simulations and experiments are illustrated in Sections 3 and 4, respectively. Finally, Section 5 summarizes this work.


## 2. Principle

### 2.1. Color Phase-Shifting Fringe

Phase-shifting algorithms are widely used in optical metrology because of their high accuracy, among which three-step phase-shifting algorithm is the most commonly chosen algorithm since it requires the least number of patterns. Three phase-shifting patterns can be separately encoded into three channels of one color fringe. Ideally, the three channels of the color fringe can be written as(1)IRx,y=A0x,y+B0x,ycosϕx,y−2π3,IGx,y=A0x,y+B0x,ycosϕx,y,IBx,y=A0x,y+B0x,ycosϕx,y+2π3,wherex,y denotes the image pixel, A0x,y is the average intensity, B0x,y is the intensity modulation, and ϕx,y is the wrapped phase. Simultaneously solving the above equations will lead to
(2)ϕx,y=tan−13IR−IB2IG−IR−IB.The arctangent function gives the wrapped phase ranging from−π to +π with 2π phase jumps. Therefore, phase unwrapping should be performed to recover the absolute phase map. The relationship between the wrapped phase ϕx,y and the absolute phase Φx,y can be described as
(3)Φx,y=ϕx,y+2πkx,y,where kx,y represents the fringe order, which can be determined by phase unwrapping algorithms.


### 2.2. Color Crosstalk Model

Most color cameras and projectors are designed to have spectrum overlaps between three channels to avoid color-blind regions in the spectrum. As a result, it is difficult to completely separate the three encoded fringes because of color crosstalk. In general, the linear color crosstalk can be described as [32]
(4)IR′IG′IB′=MIRIGIB=m11m12m13m21m22m23m31m32m33IRIGIB,where IR′,IG′,IB′ are the RGB channels of the captured fringe pattern, M represents the crosstalk matrix, and miji,j=1,2,3 are crosstalk factors. In general, the crosstalk between red and blue channels is smaller than that between red and green channels or between blue and green channels. Thus, crosstalk factors generally satisfy the following constraints:
(5)m11>m12>m13,m22>m21,m23,m31<m32<m33.Taking the red channel as an example,IR′ can be further expressed as
(6)IR′=AR+XRcosϕ−2π3+YRsinϕ−2π3,(7)AR=m11+m12+m13×A0,XR=2m11−m12−m13×B02,YR=−m12−m13×3B02.The above equations can be rewritten as(8)IR′=AR+BRcosϕ−2π3+εR,(9)BR=XR2+YR2,cosεR=XRBR,sinεR=−YRBR.SinceXR>−YR>0, εR is positive and usually very small. The phase-shift deviations of green channel and blue channel can be analyzed in the equivalent way; εG and εB are also usually very small and εB is negative. As a result, the green channel IG′ and the blue channel IB′ can be further expressed as
(10)IG′=AG+BGcosϕ+εG,(11)IB′=AB+BBcosϕ+2π3+εB.Through above analysis, it can be known that the values ofAR, AG, and AB are different (the same for BR, BG, and BB), which will cause color imbalance between three channels of color fringe. To avoid introducing phase error caused by color imbalance, fringe normalization should be performed, which will be discussed in Section 2.5. After fringe normalization, the three channels of the color fringe can be rewritten as
(12)IR′¯=cosϕ−2π3+εR,IG′¯=cosϕεG,IB′¯=cosϕ+2π3+εB.


### 2.3. Phase Error Analysis

Although the color imbalance problem is addressed after fringe normalization, there are still phase-shift deviations between three channels of color fringe due to the existence ofεR,εG,εB. If we continue to calculate the phase as
(13)ϕ′x,y=tan−13IR′¯−IB′¯2IG′¯−IR′¯−IB′¯.Then, phase error will be introduced, which can be expressed as(14)Δϕ′x,y=ϕ′x,y−ϕx,y=tan−1cosϕ∑In′¯sinδn+sinϕ∑In′¯cosδnsinϕ∑In′¯sinδn−cosϕ∑In′¯cosδn=tan−1cos2ϕ∑sinθn+sin2ϕ∑cosθn−∑sinεnsin2ϕ∑sinθn−cos2ϕ∑cosθn−∑cosεn,where n=R,G,B and θn=2δn+εn; δn=−2π/3,0,+2π/3. It can be seen that Δϕ′ highly relates to 2ϕ. To simplify the above equation, we assume that εG=0 and εB=−εR=ε⟶0−. Based on this assumption, Δϕ′ will be simplified as
(15)Δϕ′x,y=tan−1Xsin2ϕY−Xcos2ϕ,(16)X=cosε+π3−12≈−3ε2,Y=cosε+12≈32.Therefore, Equation (15) can be further approximated as [29]
(17)Δϕ′x,y≈tan−1−3εsin2ϕ3+3εcos2ϕ≈tan−1−3εsin2ϕ3≈−3εsin2ϕ3.This indicates that the phase error approximately doubles the frequency of the color fringe.


### 2.4. Phase Error Compensation

Based on this observation, if the original phase map is shifted by quarter period orπ/2, the shifted phase map can be used to compensate the quasisine phase error by averaging it with the original phase map. Therefore, the other color fringe with π/2 phase shift to the original one is also obtained as
(18)IR″¯=sinϕ−2π3+εR,IG″¯=sinϕεG,IB″¯=sinϕ+2π3+εB.Then, we can calculate the carried phase as(19)ϕ″x,y=tan−12IG″¯−IR″¯−IB″¯3IB″¯−IR″¯.The phase error can be expressed as(20)Δϕ″x,y=ϕ″x,y−ϕx,y≈3εsin2ϕ3.Clearly,Δϕ′ and Δϕ″ have the same amplitude but opposite distributions. Thence, we can compensate the phase error as
(21)ϕx,y=ϕ′x,y+ϕ″x,y2.


### 2.5. Fringe Normalization

As mentioned above, two color fringes withπ/2 phase shift are used to compensate the phase error which caused by phase-shift deviation. And the color crosstalk matrix of two color fringes can be considered the same. Namely, for the same channel, the values of An and Bn of two color fringes are equal. Thus, a background pattern is utilized to obtain the values of An. By subtracting An of each channel, the fringe intensities of two color fringes can be rewritten as
(22)In′~=Bncosϕ+δn+εn,In″~=Bnsinϕ+δn+εn.Then, they can be normalized by the following equations:(23)In′¯=In′~In′~2+In″~2,In″¯=In″~In′~2+In″~2.As shown in Figure1, the framework of the proposed method contains three steps. Firstly, a background pattern is used to remove color imbalance in different channels of cosine and sine fringes. Then, two phase maps with phase errors caused by phase-shift deviation are extracted, respectively. Finally, the phase error caused by phase-shift deviation is compensated by averaging two phase maps.Figure 1
Framework of the proposed method.


## 2.1. Color Phase-Shifting Fringe

Phase-shifting algorithms are widely used in optical metrology because of their high accuracy, among which three-step phase-shifting algorithm is the most commonly chosen algorithm since it requires the least number of patterns. Three phase-shifting patterns can be separately encoded into three channels of one color fringe. Ideally, the three channels of the color fringe can be written as(1)IRx,y=A0x,y+B0x,ycosϕx,y−2π3,IGx,y=A0x,y+B0x,ycosϕx,y,IBx,y=A0x,y+B0x,ycosϕx,y+2π3,wherex,y denotes the image pixel, A0x,y is the average intensity, B0x,y is the intensity modulation, and ϕx,y is the wrapped phase. Simultaneously solving the above equations will lead to
(2)ϕx,y=tan−13IR−IB2IG−IR−IB.The arctangent function gives the wrapped phase ranging from−π to +π with 2π phase jumps. Therefore, phase unwrapping should be performed to recover the absolute phase map. The relationship between the wrapped phase ϕx,y and the absolute phase Φx,y can be described as
(3)Φx,y=ϕx,y+2πkx,y,where kx,y represents the fringe order, which can be determined by phase unwrapping algorithms.


## 2.2. Color Crosstalk Model

Most color cameras and projectors are designed to have spectrum overlaps between three channels to avoid color-blind regions in the spectrum. As a result, it is difficult to completely separate the three encoded fringes because of color crosstalk. In general, the linear color crosstalk can be described as [32]
(4)IR′IG′IB′=MIRIGIB=m11m12m13m21m22m23m31m32m33IRIGIB,where IR′,IG′,IB′ are the RGB channels of the captured fringe pattern, M represents the crosstalk matrix, and miji,j=1,2,3 are crosstalk factors. In general, the crosstalk between red and blue channels is smaller than that between red and green channels or between blue and green channels. Thus, crosstalk factors generally satisfy the following constraints:
(5)m11>m12>m13,m22>m21,m23,m31<m32<m33.Taking the red channel as an example,IR′ can be further expressed as
(6)IR′=AR+XRcosϕ−2π3+YRsinϕ−2π3,(7)AR=m11+m12+m13×A0,XR=2m11−m12−m13×B02,YR=−m12−m13×3B02.The above equations can be rewritten as(8)IR′=AR+BRcosϕ−2π3+εR,(9)BR=XR2+YR2,cosεR=XRBR,sinεR=−YRBR.SinceXR>−YR>0, εR is positive and usually very small. The phase-shift deviations of green channel and blue channel can be analyzed in the equivalent way; εG and εB are also usually very small and εB is negative. As a result, the green channel IG′ and the blue channel IB′ can be further expressed as
(10)IG′=AG+BGcosϕ+εG,(11)IB′=AB+BBcosϕ+2π3+εB.Through above analysis, it can be known that the values ofAR, AG, and AB are different (the same for BR, BG, and BB), which will cause color imbalance between three channels of color fringe. To avoid introducing phase error caused by color imbalance, fringe normalization should be performed, which will be discussed in Section 2.5. After fringe normalization, the three channels of the color fringe can be rewritten as
(12)IR′¯=cosϕ−2π3+εR,IG′¯=cosϕεG,IB′¯=cosϕ+2π3+εB.


## 2.3. Phase Error Analysis

Although the color imbalance problem is addressed after fringe normalization, there are still phase-shift deviations between three channels of color fringe due to the existence ofεR,εG,εB. If we continue to calculate the phase as
(13)ϕ′x,y=tan−13IR′¯−IB′¯2IG′¯−IR′¯−IB′¯.Then, phase error will be introduced, which can be expressed as(14)Δϕ′x,y=ϕ′x,y−ϕx,y=tan−1cosϕ∑In′¯sinδn+sinϕ∑In′¯cosδnsinϕ∑In′¯sinδn−cosϕ∑In′¯cosδn=tan−1cos2ϕ∑sinθn+sin2ϕ∑cosθn−∑sinεnsin2ϕ∑sinθn−cos2ϕ∑cosθn−∑cosεn,where n=R,G,B and θn=2δn+εn; δn=−2π/3,0,+2π/3. It can be seen that Δϕ′ highly relates to 2ϕ. To simplify the above equation, we assume that εG=0 and εB=−εR=ε⟶0−. Based on this assumption, Δϕ′ will be simplified as
(15)Δϕ′x,y=tan−1Xsin2ϕY−Xcos2ϕ,(16)X=cosε+π3−12≈−3ε2,Y=cosε+12≈32.Therefore, Equation (15) can be further approximated as [29]
(17)Δϕ′x,y≈tan−1−3εsin2ϕ3+3εcos2ϕ≈tan−1−3εsin2ϕ3≈−3εsin2ϕ3.This indicates that the phase error approximately doubles the frequency of the color fringe.


## 2.4. Phase Error Compensation

Based on this observation, if the original phase map is shifted by quarter period orπ/2, the shifted phase map can be used to compensate the quasisine phase error by averaging it with the original phase map. Therefore, the other color fringe with π/2 phase shift to the original one is also obtained as
(18)IR″¯=sinϕ−2π3+εR,IG″¯=sinϕεG,IB″¯=sinϕ+2π3+εB.Then, we can calculate the carried phase as(19)ϕ″x,y=tan−12IG″¯−IR″¯−IB″¯3IB″¯−IR″¯.The phase error can be expressed as(20)Δϕ″x,y=ϕ″x,y−ϕx,y≈3εsin2ϕ3.Clearly,Δϕ′ and Δϕ″ have the same amplitude but opposite distributions. Thence, we can compensate the phase error as
(21)ϕx,y=ϕ′x,y+ϕ″x,y2.


## 2.5. Fringe Normalization

As mentioned above, two color fringes withπ/2 phase shift are used to compensate the phase error which caused by phase-shift deviation. And the color crosstalk matrix of two color fringes can be considered the same. Namely, for the same channel, the values of An and Bn of two color fringes are equal. Thus, a background pattern is utilized to obtain the values of An. By subtracting An of each channel, the fringe intensities of two color fringes can be rewritten as
(22)In′~=Bncosϕ+δn+εn,In″~=Bnsinϕ+δn+εn.Then, they can be normalized by the following equations:(23)In′¯=In′~In′~2+In″~2,In″¯=In″~In′~2+In″~2.As shown in Figure1, the framework of the proposed method contains three steps. Firstly, a background pattern is used to remove color imbalance in different channels of cosine and sine fringes. Then, two phase maps with phase errors caused by phase-shift deviation are extracted, respectively. Finally, the phase error caused by phase-shift deviation is compensated by averaging two phase maps.Figure 1
Framework of the proposed method.


## 3. Simulations

In this section, simulations were carried out to verify the performance of the proposed method. The first simulated object is calledpeaks function in MATLAB, and its resolution is 400×400 pixels. The fringe pitch of the patterns we used is 20 pixels. As mentioned in Section 2.2, the crosstalk matrix should satisfy Equation (5). Therefore, we choose the crosstalk matrix as
(24)M=1.000.360.080.151.000.200.060.251.00.Meanwhile, to simulate the uneven illumination, a 2D Gaussian functionGx,y with a size of 400×400 pixels and a standard deviation of 100 pixels was applied on the simulated patterns. And the uneven illumination can be expressed as
(25)Rx,y=0.5+0.5×Gx,ymaxGx,y.Figures2(a) and 2(b) plot the 300th cross sections of cosine and sine fringes before fringe normalization, respectively. Obviously, the problem of color imbalance was existing in different channels. Then, the color imbalance problem was eliminated by fringe normalization, and the fringes after fringe normalization are shown in Figures 2(c) and 2(d). Figure 3(a) shows the 3D shape obtained before fringe normalization; periodic errors can be seen clearly. As shown in Figure 3(b), these errors were decreased after fringe normalization, but the errors caused by phase-shift deviation were still existing. Figure 3(c) shows the 3D shape retrieved after phase error compensation. Apparently, periodic errors were effectively suppressed by the proposed method. The reconstructed surface is very smooth, with almost no visible errors. In addition, their corresponding root mean square errors (RMSE) are shown in Figures 3(d)–3(f). Before fringe normalization, the RMSE in Figure 3(d) is 0.1298 rad. Figure 3(e) shows that the RMSE was reduced to 0.1099 rad after fringe normalization. And the RMSE, as shown in Figure 3(f), was reduced to 0.0550 rad after phase error compensation. This simulation illustrates the effectiveness of the proposed method.Figure 2
The 300th cross sections of cosine and sine fringes. (a, b) Before fringe normalization; (c, d) after fringe normalization.Figure 3
Simulation results of thepeaks function. (a) Before fringe normalization; (b) after fringe normalization; (c) after phase error compensation; (d–f) RMSE of (a–c), respectively.To compare the proposed method with Wang’s method [31] which uses frequency domain method, a staircase with four steps was simulated. The 3D shapes reconstructed by the proposed method and Wang’s method are shown in Figures 4(a) and 4(b), respectively. Obviously, Wang’s method performed worse than the proposed method because it uses frequency domain method, which is not suitable for processing objects with abrupt surfaces. Figure 5 illustrates a part of the 300th cross section, which contains the abrupt change between the third and fourth steps. The third and fourth stairs are marked in green and yellow in Figure 4, respectively. The distortion of the abrupt change retrieved by Wang’s method is more serious than the proposed method. This simulation demonstrates that the proposed method can handle objects with abrupt surfaces well.Figure 4
Reconstructions: (a) the proposed method; (b) Wang’s method.Figure 5
A part of the 300th cross section, which contains the abrupt change between the third and fourth steps.


## 4. Experiments

To further verify the proposed method, a FPP system was set up comprising a digital-light-processing (DLP) projector (Light-Crafter 4500) and a complementary-metal-oxide-semiconductor (COMS) camera (Point grey Chameleon3). The resolution of the camera is1280×1024 pixels, and the optical lens (Kowa LM16JCM) mounted on the camera has a focal length of 16 mm. The resolution of the projector is 912×1140 pixels. It is well known that high-frequency fringes can effectively resist noise and provide high-precision measurement [10]. However, using high-frequency fringes also introduces phase ambiguities in the reconstruction process that need to be addressed by means of phase unwrapping. In these experiments, three-frequency method [7] was employed to unwrap the wrapped phase. Here, Th, Tm, and Tl represent periods of high-, medium-, and low-frequency fringes, respectively. And they satisfy the following equation:
(26)N2Th=NTm=Tl,where N is a constant. Corresponding wrapped phases of three frequency are ϕh, ϕm, and ϕl. In addition, Tl is usually greater than the projector’s horizontal and vertical resolution to ensure complete elimination of phase ambiguity, which means the wrapped phase of low-frequency fringes has no phase ambiguity. It can be mathematically expressed as
(27)Φl=ϕl,where Φl is the absolute phase of low-frequency fringes. Then, the fringe order km and absolute phase Φm of medium-frequency fringes can be calculated as
(28)km=roundΦl∗N−ϕm2π,Φm=ϕm+2π∗km,where round returns the closest integer. After determining Φm, the fringe order kh and absolute phase Φh of high-frequency fringes can be obtained similarly:
(29)kh=roundΦm∗N−ϕh2π,Φh=ϕh+2π∗kh.Certainly, the errors caused by noise will inevitably be introduced in the actual measurement. Let us assume that the noise errors of low- and medium-frequency fringes areΔΦl and ΔΦm, respectively. Then, the actual fringe orders obtained are
(30)km′=roundΦl+ΔΦl∗N−ϕm2π=km+roundΔΦl∗N−ϕm2π,kh′=roundΦm+ΔΦm∗N−ϕh2π=kh+roundΔΦm∗N−ϕh2π.Obviously, if the value ofN is too large, the fringe orders of medium- and high-frequency fringes will be calculated wrong, which will result in unwrapping errors. In our measurement system, N=8 is a suitable value. As mentioned above, the period of low-frequency fringes should be greater than the projector’s horizontal and vertical resolution. Therefore, the period of high-frequency fringes should satisfy
(31)Th=TlN2≥1140N2≈17.8.In this paper, we chooseTh=18 pixels. Thus, Tm=NTh=144 pixels and Tl=N2Th=1152 pixels.In the first experiment, a ball was measured. Total seven patterns were projected onto the ball and captured by the camera, as shown in Figure6. Figures 6(a)–6(c) show the deformed cosine fringe patterns with periods of 18, 144, and 1152 pixels, respectively. Their corresponding sine fringe patterns are shown in Figures 6(d)–6(f). Figure 6(g) shows the background pattern that is used for fringe normalization. Before phase error compensation, we need to normalize these fringes to address the color imbalance problem. Figure 7 shows the 500th cross sections of six color patterns. Obviously, the imbalance problem existing in Figures 7(a)–7(f) was effectively resolved after fringe normalization, and the normalized fringes are shown in Figures 7(g)–7(l). Afterwards, the reconstructed 3D shapes are shown in Figure 8. Among them, Figure 8(a) shows the 3D shape reconstruction before fringe normalization; periodic phase errors are clearly visible. Then, the fringe normalization method was employed, and the 3D shape retrieved after fringe normalization is shown in Figure 8(b). The phase errors caused by color imbalance were suppressed, but the phase errors caused by phase-shift deviation were still existing. Thus, Wang’s [31] and the proposed methods were both applied to compensate these errors, and the reconstructions obtained are shown in Figures 8(c) and 8(d), respectively. Visually, both Wang’s and the proposed methods effectively suppressed the phase errors in the center of the ball. However, the frequency domain method used by Wang caused large errors in the edge areas, especially the upper edge of the ball. In order to quantify errors, we also used 18-step phase-shifting method [8] to measure this ball and compared the measurement result with the above four processing results. The corresponding deviation maps are shown in Figures 8(e)–8(h). The RMSE of the reconstruction before fringe normalization is 0.5682 rad, and it is reduced to 0.1255 rad after fringe normalization. As shown in Figures 8(g) and 8(h), after phase error compensation, the RMSEs of Wang’s and the proposed methods are 0.0550 and 0.0368 rad, respectively. And the maximum phase error of Wang’s method reaches 0.8 rad in edge areas which is 5.3 times that of the proposed method. In addition, Figure 9 plots the 800th cross section of the center of the ball. This experiment shows that both Wang’s method and the proposed method are reliable in processing smooth object, but the proposed method handles edge areas better than Wang’s method.Figure 6
The ball. (a–c) Deformed cosine fringe patterns with periods of 18, 144, and 1152 pixels, respectively; (d–f) deformed sine fringe patterns with periods of 18, 144, and 1152 pixels, respectively; (g) the background pattern.Figure 7
The 500th cross sections. (a–c) Cosine fringes with periods of 18, 144, and 1152 pixels before fringe normalization, respectively; (d–f) sine fringes with periods of 18, 144, and 1152 pixels before fringe normalization, respectively; (g–i) cosine fringes with periods of 18, 144, and 1152 pixels after fringe normalization, respectively; (j–l) sine fringes with periods of 18, 144, and 1152 pixels after fringe normalization, respectively.Figure 8
Reconstructions. (a) Before fringe normalization; (b) after fringe normalization; (c) Wang’s method; (d) the proposed method; (e–h) deviation maps of (a–d), respectively.Figure 9
The 800th cross section.In order to further compare the proposed method with Wang’s method in dealing with object with abrupt surface, a staircase with four steps was measured. All captured patterns are shown in Figure10. Similarly, the 18-step phase-shifting method was used as a reference. Figures 11(a)–11(c) show the 3D reconstructions recovered by 18-step phase-shifting, Wang’s, and the proposed methods, respectively. Visually, the 18-step phase-shifting method performed best while Wang’s method worked the worst. Due to the application of frequency domain method, the abrupt changes between steps were smoothed, resulting in incorrect measurement. The abrupt changes between steps reconstructed by the proposed method are obviously better than Wang’s method. Furthermore, we performed plane fitting on each step of the reconstruction retrieved by 18-step phase-shifting method. Figures 11(d)–11(f) show the deviation maps between the three reconstructions and the fitting results. Markedly, the accuracy of the proposed method is similar to that of the 18-step phase-shifting method, with only a difference of 0.01 rad. However, the RMSE of Wang’s method is approximately twice that of the proposed method. Figure 12 plots the 800th cross section of the staircase. The proposed method almost coincides with 18-step phase-shifting method, while Wang’s method performs poorly. This experiment confirms that the proposed method is more robust and accurate than Wang’s method in dealing with object with abrupt surfaces.Figure 10
The staircase. (a–c) Deformed cosine fringe patterns with periods of 18, 144, and 1152 pixels, respectively; (d–f) deformed sine fringe patterns with periods of 18, 144, and 1152 pixels, respectively; (g) the background pattern.Figure 11
Reconstructions. (a) 18-step phase-shifting method; (b) Wang’s method; (c) the proposed method; (d–f) deviation maps of (a–c), respectively.Figure 12
The 800th cross section.Two isolated sculptures were also measured to verify the actual measurement effect of the proposed method. All captured images are shown in Figure13. Figure 14 shows the reconstructed 3D shapes, and the reconstruction in Figure 14(c) is significantly better than that in Figures 14(a) and 14(b). Figures 14(d)–14(f) show their corresponding deviation maps. Compared with the result before fringe normalization, the proposed method significantly improved the measurement accuracy, and the RMSE of the proposed method is only one-fifth of that before fringe normalization. This experiment demonstrates that the proposed method has sufficient robustness and accuracy in the actual measurement.Figure 13
Two sculptures. (a–c) Deformed cosine fringe patterns with periods of 18, 144, and 1152 pixels, respectively; (d–f) deformed sine fringe patterns with periods of 18, 144, and 1152 pixels, respectively; (g) the background pattern.Figure 14
Reconstructions. (a) Before fringe normalization; (b) after fringe normalization; (c) after phase error compensation; (d–f) deviation maps of (a–c), respectively.To further validate the dynamic measurement performance of the proposed method, one moving hand was measured. Figures15(a)–15(d) show the 3D shapes reconstructed without color crosstalk compensation, and phase errors can be seen clearly. In contrast, the corresponding frames processed by the proposed method have no obvious errors, as shown in Figures 15(e)–15(h). This experimental result further confirms that the proposed method can effectively compensate the phase errors and improve measurement accuracy in dynamic measurement.Figure 15
The moving hand. (a–d) Without color crosstalk compensation; (e–h) the proposed method.


## 5. Conclusion

This paper presents an effective color crosstalk correction method for color FPP without crosstalk calibration. Due to the fact that phase error induced by color crosstalk is approximately twice the frequency of color fringe, the other sine color fringe was used to compensate the phase error. Both simulations and experiments confirmed the effectiveness of the proposed method. Experimental data shows that the proposed method has the ability to deal with objects with abrupt surfaces, and phase error caused by color crosstalk can be reduced to one-fifth after applying the proposed method.



---

*Source: 7333291-2022-01-28.xml*