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

Generalized phase unwrapping method that avoids jump errors for fringe projection profilometry

Open Access Open Access

Abstract

Jump errors easily occur on the discontinuity of the wrapped phase because of the misalignment between wrapped phase and fringe order in fringe projection profilometry (FPP). In this paper, a phase unwrapping method that avoids jump errors is proposed for FPP. By building two other staggered wrapped phases from the original wrapped phase and dividing each period of fringe order into three parts, the proposed generalized tripartite phase unwrapping (Tri-PU) method can be used to avoid rather than compensatorily correct jump errors. It is suitable for the phase unwrapping method assisted by fringe order with a basic wrapped phase and fringe order, no matter which method is used to recover them. The experimental results demonstrate the effectiveness and generality of the proposed method, which is simple to implement and superior to measure complex objects with sharp edges.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Phase unwrapping is a key issue for three-dimensional (3D) shape measurement based on fringe projection profilometry (FPP) [13]. In FPP, the retrieved phase is wrapped between -π and π due to the inverse trigonometric operation, which is ambiguous for further phase-to-height mapping. To eliminate phase ambiguity, spatial and temporal phase unwrapping methods have been developed. The spatial phase unwrapping method detects and removes the 2π phase jumps by judging the phase difference between adjacent pixels, while the temporal phase unwrapping (TPU) method eliminates phase ambiguity by projecting additional patterns to uniquely label each period of the wrapped phase. Because of the pixel-by-pixel phase unwrapping, the TPU method is more suitable to be applied in complex and isolated scenes.

TPU methods can be further classified into the phase-assisted method and the intensity-assisted method depending on the different ways to modulate and retrieve the information of the fringe order. The former mainly includes multiple-frequency or multiple-wavelength method [2] and phase-coding method [4], and the latter mainly includes De Bruijn coding method [5], Grey-level coding method [6] and Gray code coding method [7]. All the TPU methods eliminate phase ambiguity by matching the wrapped phase and the corresponding fringe order. But in actual measurement, TPU methods suffer two types of phase unwrapping errors, which seriously affect the final reconstructed 3D result. The first one is the phase unwrapping error caused by noise, which is uniformly distributed and can be easily eliminated by simple filtering. The second one is the error caused by arctangent calculation and misalignment between the wrapped phase and fringe order. This kind of errors is periodically distributed and gather on the discontinuities of the wrapped phase, so it is named the jump errors. Jump errors are in zonal distribution and the width of the error region further increases in dynamic measurement due to the motion of the object and the active defocusing of the projector [8]. Therefore, it is more challenging to detect and remove the jump errors compared with the noise errors. It should be noted that jump errors do not occur in multiple-frequency methods or multiple-wavelength methods because solving phase order in these methods do not only use the assisted low-frequency wrapped phase but utilize the high-frequency wrapped phase which will be unwrapped. So, the unstable jump of the phase order in the discontinuity is well-matched with the jump of wrapped phase. Therefore, multiple-frequency methods or multiple-wavelength methods are not included in the scope of this work.

The existing methods to overcome this challenge can be classified into post-corrected methods and pre-avoided methods. Among these methods, spatial phase unwrapping [9], median filtering [10], monotonicity detection [1113], and robust principal component analysis [14] are the post-corrected methods to detect and correct jump errors. These kinds of methods perform well in static scenes. But it is hard to apply them in dynamic measurement because the width of the error region is larger in dynamic scenes, causing difficulty to distinguish the jump errors from the abrupt change of the object. Pre-avoided methods directly avoid the jump errors by staggering the error regions of wrapped phase and the fringe order. Complementary Gray-coded (CGC) method [15] obtain two staggered fringe orders to overcome the jump errors by projecting an additional pattern, which decreases the measuring efficiency. Cyclic complementary Gray-coded (CCGC) and shifting Gray-coded (SGC) methods [16,17] achieved robust phase unwrapping for complex dynamic scenes by utilizing Gray codes in adjacent projected sequences, which is not suitable for static scenes. Tripartite phase unwrapping (Tri-PU) method [18] solves the problem of jump errors by changing the built-up sequence of sinusoidal patterns to obtain staggered wrapped phases and calculating the wrapped phase to divide the fringe order into three parts, which must satisfy the precondition of three-step phase shifting and obtaining reference phase.

To sum up, post-corrected methods have better adaptability and can be applied in all phase unwrapping methods in FPP, but not suitable for dynamic scenes. Pre-avoided methods fundamentally avoid jump errors and can be applied in different measurement scenes but they are not generalized approaches for FPP regardless of phase-shifting steps and phase unwrapping methods. Therefore, an adaptive and generalized pre-avoided phase unwrapping method is required to solve the problem of jump errors in FPP.

In this paper, a phase unwrapping method that avoids jump errors is proposed for FPP. By building two other staggered wrapped phases from the original wrapped phase and dividing each period of fringe order into three parts, the proposed generalized Tri-PU method can be used to avoid rather than compensatorily correct the jump errors. No additional information is required to assist to divide region and no restrictions on the way of obtaining the wrapped phase and fringe order. So, this method is generalized for phase unwrapping, no matter which method is used to obtain wrapped phase and its corresponding fringe order.

2. Principle

2.1 Typical phase unwrapping errors in TPU methods

In TPU methods, sinusoidal patterns are projected to modulate the height information of a measured object, and the modulated wrapped phase can be extracted using different methods. But all the wrapped phase are wrapped in (-π, π] due to arctangent function. To eliminate the phase ambiguity, fringe order is required and can be solved from coding patterns.

In actual measurement, two types of phase unwrapping errors easily occur due to different reasons as shown in Fig. 1. The first one is the error caused by noise, which is uniformly distributed and discrete. So, it can be easily detected and eliminated by the reviewed post-eliminated methods. The second one is the jump errors caused by two main reasons. Firstly, arctangent calculation in phase extraction leads to the unstable jump near the discontinuity of the wrapped phase, therefore the jump edge of the wrapped phase is unreliable for phase unwrapping. In addition, the motion of the object and the defocusing of the projector cause the misalignment between the wrapped phase and fringe order, which leads to the obvious errors in phase unwrapping. The jump errors are periodically distributed and in zonal distribution, so it is hard for the post-eliminated methods to distinguish the jump errors from the abrupt edges of the object when the width of error region is large. Therefore, the jump errors are more difficult to be eliminated compared with the noise-caused errors in FPP. So, in this work, we focus on solving the problem of the latter error (jump errors) and finding a generalized method which can be adapted to different measuring system (focusing and defocusing measuring system) and widely applied in different methods (coding methods and phase extraction methods) in FPP.

 figure: Fig. 1.

Fig. 1. Illustration of typical phase unwrapping errors in FPP.

Download Full Size | PDF

2.2 Proposed generalized phase unwrapping method

In our recent work, the pre-avoided Tri-PU method [18] is proposed to avoid jump errors in Gray-code-based phase unwrapping, in which three parts of divided fringe orders are used to unwrap three corresponding staggered wrapped phases. But the precondition of three-step phase shifting and reference wrapped phase in Tri-PU method dramatically limits its generality in FPP.

In this work, a generalized Tri-PU phase unwrapping method that avoided jump errors is proposed for FPP. The framework of the proposed method is shown in Fig. 2. In TPU methods, sinusoidal patterns are projected to modulate the height information of a measured object, and the modulated wrapped phase can be extracted using different methods such as phase-shifting profilometry (PSP) [19], Fourier transform profilometry (FTP) [20], window Fourier transform profilometry (WFTP) [21] and so on. To eliminate the phase ambiguity, fringe order is required and can be calculated by phase coding [4], De Bruijn coding [5], Gray-level coding [6], Gray code coding [7] and so on. In traditional TPU methods, the discontinuities of the wrapped phase and the fringe order should be well aligned to avoid the jump errors, which is hard and even impossible in the actual measurement. In the proposed method, the discontinuity of wrapped phase and fringe order are pre-staggered by building other two staggered wrapped phases from the original wrapped phase and dividing each period of fringe order into three parts. Three parts of the fringe order are used to unwrap the middle part of the corresponding three staggered wrapped phases which has no discontinuity, so the proposed method can be used to avoid rather than compensatorily correct the jump errors.

 figure: Fig. 2.

Fig. 2. Framework of the proposed generalized phase unwrapping method.

Download Full Size | PDF

Taking the phase shifting method as an example to illustrate the principle of the proposed method, the ith deformed fringes can be described as in N-step PSP:

$${I_i}(x,y) = a(x,y) + b(x,y)\cos [\Phi (x,y) - {\delta _i}],i = 1,2,3\ldots N$$
where, a(x,y) and b(x,y) are the average and modulation of patterns; Φ(x,y) is the modulated continuous phase; δi=2π(i-1)/N is the phase-shift amount. And the phase shifting algorithm is used to calculate the wrapped phase ϕ(x,y):
$$\phi (x,y) = \textrm{artan} [\frac{{\sum {_{i = 1}^N(I_i^{}(x,y)\sin {\delta _i})} }}{{\sum {_{i = 1}^N(I_i^{}(x,y)\cos {\delta _i})} }}].$$
The periodic natures of the sinusoidal patterns and arctangent function introduce phase ambiguity, so the fringe (phase) order k(x,y) is required to eliminate phase ambiguity.

In actual measurement, jump errors easily occur on the discontinuity as shown in Fig. 3(a). The essential cause of edges errors is the misaligment of jump edge between wrapped phase and fringe order, so phase unwrapping in the edge region of wrapped phase is unreliable. But in the middle region of wrapped phase, both fringe order and wrapped phase are continuous so phase unwrapping on this region is realiable. Based on this point, a generalized phase unwrapping method is proposed to avoid jump errors.

 figure: Fig. 3.

Fig. 3. Principle of the proposed generalized Tri-PU phase unwrapping method. (a) Schematic diagram of region division. (b) Process of tripartite phase unwrapping.

Download Full Size | PDF

Firstly, each period of fringe order A(i) is divided into three parts. The middle part of A(i) is judged by realiable phase value of wrapped phase ϕ(x,y):

$${k_m}(x,y) = k(x,y),\textrm{ }where\textrm{ }|{\phi (x,y)} |< \pi /3.$$
And the left and right parts of A(i) are distinguished by the minmum point of wrapped phase:
$${k_l}(x,y) = k(x,y),\textrm{ }where\textrm{ }y < {y_{\left\langle {\min |{\phi (x)} |} \right\rangle }},\textrm{ }(x,y) \in A(i)\textrm{ }and\textrm{ }(x,y) \notin {k_m}$$
$${k_r}(x,y) = k(x,y),\textrm{ }where\textrm{ }y > {y_{\left\langle {\min |{\phi (x)} |} \right\rangle }},\textrm{ }(x,y) \in A(i)\textrm{ }and\textrm{ }(x,y) \notin {k_m}$$
in which, y denotes the direction perpendicular to fringe and <min|ϕ(x)|> denotes the position of the minmum wrpped phase in each row.

After region division of k(x,y), three staggered wrapped phases are calculated using following equation:

$$\left\{ {\begin{array}{l} {{\phi_1}(x,y)\textrm{ = }\boldsymbol{wrap}[\phi (x,y) + 2\pi /3]}\\ {{\phi_2}(x,y) = \phi (x,y)\textrm{ }}\\ {{\phi_3}(x,y)\textrm{ = }\boldsymbol{wrap}[\phi (x,y) - 2\pi /3]} \end{array}} \right.,$$
in which,
$$\boldsymbol{wrap}[\phi ] = \left\{ {\begin{array}{l} {\phi - 2\pi ,\textrm{ }\phi > \pi }\\ {\phi + 2\pi ,\textrm{ }\phi \le - \pi } \end{array}} \right..$$
Finally, three parts of fringe order kl(x,y), km(x,y) and kr(x,y) can be used to unwrap the middle part of the corresponding wrapped phases ϕ1(x,y), ϕ2(x,y) and ϕ3(x,y) as follows:
$$\Phi (x,y) = \left\{ {\begin{array}{l} {{\phi_1}(x,y) + 2\pi {k_l}(x,y) - 2\pi /3}\\ {{\phi_2}(x,y) + 2\pi {k_m}(x,y)\textrm{ }}\\ {{\phi_3}(x,y) + 2\pi {k_r}(x,y) + 2\pi /3} \end{array}} \right..$$
As shown in Fig. 3(b), the discontinuity of wrapped phase and fringe order are staggered in each subregion of k(x,y), so the jump errors are pre-avoided in the proposed phase unwrapping method rather than post-corrected in the existing methods.

Compared with our recently proposed Tri-PU method [18], this generalized Tri-PU method doesn’t need additional reference plane for region division of fringe order and has no limit for the way of phase extraction. The proposed method is effective only if the wrapped phase and fringe order are properly obtained using any approach. Therefore, it is a generalized phase unwrapping method in FPP.

To further illustrate the proposed method, the algorithm flow chart and the MATLAB codes are shown in Fig. 4. The symbols and acquirement of key arguments are described as follows:

  • [m,n] is the size of the image and n is the direction perpendicular to fringe direction;
  • k is the number of fringe period;
  • KK is the fringe order, which can be calculated by different coding methods;
  • Mo is the binary modulation template, which can be solved from the sinusoidal fringe pattern [22,23];
  • Mo_s and Mo_ss are binary templates of each class of fringe order before and after background removing;
  • phase_H is the original wrapped phase, which can be calculated by different phase analysis methods and phase_H_part is the prat of wrapped phase in each class of fringe order;
  • phase_H1 and phase_H2 are the staggered wrapped phase using Eq. (6);
  • thre denotes the position of the minmum wrpped phase in each row;
  • MMo_l, MMo_m, and MMo_r are the binary templates for the left, middle and right parts of the fringe order.

 figure: Fig. 4.

Fig. 4. Algorithm flow chart and the MATLAB codes of generalized Tri-PU method.

Download Full Size | PDF

3. Experiments and results

3.1 Experiments on a 1-bit defocusing projecting measuring system

To verify the effectiveness of the proposed method. We firstly developed a 1-bit binary defocusing projecting measuring system which easily suffers jump errors due to the defocusing of the projector. This system includes a DMD projector (LightCrafter4500) and a Photron FASTCAM Mini UX100 camera with a 16-mm imaging lens. The resolutions of the projector and the camera are 912×1140 pixels (896×1140 pixels are used to generate fringes) and 1280×800 pixels, respectively. This system is developed for high-speed measurement by projecting 1-bit patterns at 2000 fps.

In the first experiment, a complex honeycomb metal specimen is tested using the traditonal Gray coded method and the proposed generalized Tri-PU method. Three phase-shifting sinsoidal patterns with 32 fringe periods are projected for phase modulation and one of these is shown in Fig. 5(a). In addition, 5 Gray code pattens are projected to label the fringe order. Three-step phase-shifitng algorithm is used to obtain the traditional wrapped phase ϕ2, as shown in Fig. 5(b). And other staggered wrapped phases ϕ1 and ϕ3 shown in Fig. 5(c) and 5(d) are calculated using Eq. (6) and Eq. (7). Then, the decoding fringe order is divided into three parts kl, km, kr using Eqs. (3) - (5) and shown in Figs. 5(e) and 5(h). Finally, kl, km and kr are used to unwrap the middle part of the corresponding wrapped phase ϕ1, ϕ2 and ϕ3 to avoid the jump errors. The reconstructed results using the traditional and this proposed method are shown in Figs. 5(f) and 5(g) and the enlarged detailed information of results are shown in Figs. 5(i) and 5(j). Experimental results indicate the proposed method can well avoid the jump errors during the complex object 3D measuring using binary defocusing projecting system.

 figure: Fig. 5.

Fig. 5. Experiments on a 1-bit defocusing projecting measuring system. (a) One of the captured deformed fringe pattern. (b) Original wrapped phase. (c)-(d) Staggered wrapped phases. (e) Fringe order. (f) Reconstructed absolute phase using the traditional phase unwrapping method. (g) Reconstructed absolute phase using the proposed generalized Tri-PU method. (h) Divided tripartite regions of fringe order. (i) Enlarged subimage in (f). (j) Enlarged subimage in (g).

Download Full Size | PDF

To compare the performance the proposed pre-avoided method with the traditional post-corrected method in the complex measuring scenes, the measuring result of the honeycomb metal specimen with jump errors as shown in Fig. 6(a) are corrected using the post-corrected method. Different from the proposed pre-avoided method, the post-corrected method directly detects and corrects the jump errors using different algorithms, in which median filtering is the typical one. So, in this experiment, median filtering is applied to smooth jump errors and jump errors is detected and corrected by substracting two unwrapping phases before and after filtering. The size of the filtering window increases from 3×3 to 9×9 and the corresponding filtering results and corrected results are shown in Figs. 6(b)–6(d) and Figs. 6(f)–6(h), respectively.

 figure: Fig. 6.

Fig. 6. Comparative experiments using pre-avoided and post-corrected methods. (a) Reconstructed result using the traditional phase unwrapping method. (b)-(d) Filtering results using a filtering window with increasing size. (e) Reconstructed result using the generalized Tri-PU method. (f)-(h) Corresponding corrected results using post-corrected method.

Download Full Size | PDF

Results indicate that jump errors are gradually eliminated with the increasing size of the filtering window, but the jump errors region with large width still exist even using relatively large filtering window (9×9) as shown in Fig. 6(h). In addition, sharp edges of the measured object will be also eliminated after median filtering with a large window as shown in Fig. 6(d), and it causes the emerging errors as shown in Fig. 6(h). Therefore, it is hard to distinguish sharp edges from the phase jump errors and might suffer undesired emerging errors due to the strategy of dectecting and correcting in the post-corrected methods. But for the proposed method, it can directly avoid the jump errors while well perserving the sharp edges in the measured scene as shown in Fig. 6(e). The comparative experiments demonstrate that the proposed method is more adaptive and robust compared with the traditional post-corrected methods, so the proposed method can be applied in the complex measuring scene.

3.2 Experiments on an 8-bit focusing projecting measuring system

To verify the generality of the proposed method. We also developed a 8-bit focusing projecting measuring system including a DMD projector (LightCrafter4500) and a Baumer HXC40 camera with a 16-mm imaging lens, which is more common to be applied in FPP. The resolutions of the projector and camera are the same as that of the first experiment. This system is developed for high-accuracy measurement by projecting high-quality 8-bit patterns at 60 fps.

In the second experimrnt, a complex scene including three isolated objects are measured using different phase calculation and unwrapping methods with different fringe frequency. PSP and FTP are two representative phase analysis methods. In this experiment, PSP and π-shift FTP [24] are used to retrive wrapped phase and several typical fringe coding methods are used to calculate the fringe order to verify the generality of the proposed method. Firstly, phase coding method is used to unwrap the PSP-calculated wrapped phase (f=8, 4 patterns for phase retrieval and 3 for fringe order calculation); Grey-level coding method is used to unwrap the PSP-calculated wrapped phase (f=16, 4 patterns for phase retrieval and 2 for fringe order calculation); Gray coding method is used to unwrap the PSP-calculated wrapped phase (f=32, 4 patterns for phase retrieval and 5 for fringe order calculation). In addition, Gray coding method is also used to unwrap the π-shift FTP-calculated wrapped phase (f=64, 2 patterns for phase retrieval and 6 for fringe order calculation). The captured coding patterns of four combining methods for solving fringe order are shown in Fig. 7(a) and one of the sinusoidal patterns in four combining methods is given in Fig. 7(b). Figure 7(c) shows the extracted wrapped phase with different frequency. And the reconstructed results using the traditional phase unwrapping method and the proposed method are shown in Figs. 7(d) and 7(e), respectively. And Fig. 7(f) shows the details of the enlargred subimages Figs. 7(d) and 7(e). The experimental results indicate the proposed generalized Tri-PU method can well avoid the jump errors for different phase-calculated methods and different coding methods with different fringe frequency. And all of them verify the generality of this proposed method.

 figure: Fig. 7.

Fig. 7. Experimental results by different phase calculation and unwrapping methods using a 8-bit focusing projecting measuring system. (a) Captured coding patterns. (b) Captured fringe patterns. (c) Wrapped phase. (d) Reconstructed height using the traditional phase unwrapping method. (e) Reconstructed height using the proposed generalized Tri-PU method. (f) Enlarged subimages of the details in (d) and (e).

Download Full Size | PDF

4. Discussion

Our proposed method has the following advantages compared with the existing correction method for jump errors.

  • Pre-avoided method. The core idea of traditional post-corrected methods such as spatial phase unwrapping [9], median filtering [10], monotonicity detection [1113], and robust principal component analysis [14] is detecting the jump errors and compensating the unwrapping phase. This kind of method is well suitable for the static and slowly varying scene with narrow error region, but it is hard to correct the wide error region caused by the defocusing of the projector and motion of objects. When the measuring scene has abrupt change or discontinuous surface, it is diffcult to distinguish sharp edges from the phase jump errors and unwrapping phase might suffer undesired emerging errors using the post-corrected methods. In addition, the parameters of the post-corrected methods should be carefully chosen to stimultaneously ensure accurate detection of the jump errors and preserve abrupt change of measuring scene. But for the proposed method in this work, the jump errors are pre-avoided by staggering the discontinuity of the wrapped phase and fringe order, hence it is simple and superior to measure complex scenes with abrupt change.
  • Generalized phase unwrapping method. The previous pre-avoided methods including CGC [15], CCGC [16], SGC [17] and Tri-PU [18] have respective limitation in measuring methods. CGC, CCGC and SGC methods only work for Gray-coded method. CGC method is suitable for statics scenes by projecting one additional pattern, and CCGC and SGC methods are only suitable for dynamic scenes using Gray codes in the adjacent projecting sequences. The Tri-PU method is limited to phase extration using three-step phase shifting algorithm and additional reference wrapped phase is required to divide region of the fringe order. Compared with the mentioned pre-avoided methods, the proposed method is suitable for the different phase unwrapping methods assisted by fringe order in FPP, once the wrapped phase vlaue and the fringe order information is obtained. Therefore, this method can be widely applied to avoid the jump error in FPP as a generalized method.
Of course the proposed method also has a critical condition to ensure the correct phase unwrapping, because the proposed generalized Tri-PU method has a threshold 1/3 fringe period (2π∕3 phase difference). When the mismatch between the wrapped phase and the phase order is no more than 1/3 fringe period, the proposed method works well.

To better illustrate this question, a schematic diagram of the critical value is shown in Fig. 8. Both the shift of the wrapped phase and phase order will widen the mismatch between them. Motion of the object makes phase order have a direct shift of Δorder and also causes an indirect phase error Δϕ in wrapped phase as shown in Fig. 8. The motion-induced phase error also leads to a shift of wrapped phase’s discontinuity and the function D(·) is used to denote this mapping relationship from phase error to shifting distance.

 figure: Fig. 8.

Fig. 8. Critical condition of generalized Tri-PU method to ensure successful phase unwrapping.

Download Full Size | PDF

It can be seen from Fig. 8 that when the mismatch between the wrapped phase and phase order exceeds 1/3 fringe period T, the middle part of phase order km will be used to unwrap the discontinuity of wrapped phase, which will lead to the failure of the proposed method. Therefore, the critical condition of the proposed method to ensure the correct phase unwrapping is concluded in Eq. (9).

$$\max [\frac{{{\Delta _{order}} + D(\Delta \phi )}}{T}] < \frac{1}{3}$$

5. Conclusion

In conclusion, a generalized phase unwrapping method that avoids jump errors is proposed for FPP in this work. By building other two staggered wrapped phases from the original wrapped phase and dividing each period of fringe order into three parts, the proposed generalized Tri-PU method can be used to avoid the jump errors in FPP. Because no additional information is required except for basic wrapped phase and fringe order, this method is suitable for the phase unwrapping method assisted by fringe order in FPP, no matter which method is used to recover the wrapped phase vlaue and the fringe order information.

The experimental results demonstrated the effectiveness and generality of this proposed method. In addition, it is worth noting that the proposed method is the pre-avoid rather than post-corrected method, therefore it is superior to measure complex objects with sharp edges.

Funding

National Natural Science Foundation of China (62075143); National Postdoctoral Program for Innovative Talents (BX2021199).

Disclosures

The authors declare no conflicts of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

1. D. C. Ghiglia and M.D. Pritt, Two-dimensional phase unwrapping: theory, algorithms, and software (John Wiley and Sons, 1998).

2. C. Zuo, L. Huang, M. Zhang, Q. Chen, and A. Asundi, “Temporal phase unwrapping algorithms for fringe projection profilometry: A comparative review,” Opt. Laser Eng. 85, 84–103 (2016). [CrossRef]  

3. X. He and Q. Kemao, “A comparative study on temporal phase unwrapping methods in high-speed fringe projection profilometry,” Opt. Laser Eng. 142, 106613 (2021). [CrossRef]  

4. Y. Wang and S. Zhang, “Novel phase-coding method for absolute phase retrieval,” Opt. Lett. 37(11), 2067–2069 (2012). [CrossRef]  

5. J. Pages, J. Salvi, C. Collewet, and J. Forest, “Optimised De Bruijn patterns for one-shot shape acquisition,” Image Vision Comput. 23(8), 707–720 (2005). [CrossRef]  

6. R. Porras-Aguilar, K. Falaggis, and R. Ramos-Garcia, “Optimum projection pattern generation for grey-level coded structured light illumination systems,” Opt. Laser Eng. 91, 242–256 (2017). [CrossRef]  

7. G. Sansoni, M. Carocci, and R. Rodella, “Three-dimensional vision based on a combination of gray-code and phase-shift light projection: analysis and compensation of the systematic errors,” Appl. Opt. 38(31), 6565–6573 (1999). [CrossRef]  

8. S. Lei and S. Zhang, “Flexible 3-D shape measurement using projector defocusing,” Opt. Lett. 34(20), 3080–3082 (2009). [CrossRef]  

9. Y. Wang, S. Zhang, and J. H. Oliver, “3D shape measurement technique for multiple rapidly moving objects,” Opt. Express 19(9), 8539–8545 (2011). [CrossRef]  

10. D. Zheng, F. Da, Q. Kemao, and H. S. Seah, “Phase-shifting profilometry combined with Gray-code patterns projection: unwrapping error removal by an adaptive median filter,” Opt. Express 25(5), 4700–4713 (2017). [CrossRef]  

11. S. Zhang, “Flexible 3D shape measurement using projector defocusing: extended measurement range,” Opt. Lett. 35(7), 934–936 (2010). [CrossRef]  

12. L. Huang and A. Asundi, “Phase invalidity identification framework with the temporal phase unwrapping method,” Meas. Sci. Technol. 22(3), 035304 (2011). [CrossRef]  

13. S. Feng, Q. Chen, C. Zuo, R. Li, G. Shen, and F. Feng, “Automatic identification and removal of outliers for high-speed fringe projection profilometry,” Opt. Eng. 52(1), 013605 (2013). [CrossRef]  

14. Y. Zhang, J. Tong, L. Lu, J. Xi, Y. Yu, and Q. Guo, “Fringe Order Correction for Fringe Projection Profilometry Based on Robust Principal Component Analysis,” IEEE Access 9, 23110–23119 (2021). [CrossRef]  

15. Q. Zhang, X. Su, L. Xiang, and X. Sun, “3-D shape measurement based on complementary Gray-code light,” Opt. Laser Eng. 50(4), 574–579 (2012). [CrossRef]  

16. Z. Wu, C. Zuo, W. Guo, T. Tao, and Q. Zhang, “High-speed three-dimensional shape measurement based on cyclic complementary Gray-code light,” Opt. Express 27(2), 1283–1297 (2019). [CrossRef]  

17. Z. Wu, W. Guo, and Q. Zhang, “High-speed three-dimensional shape measurement based on shifting Gray-code light,” Opt. Express 27(16), 22631–22644 (2019). [CrossRef]  

18. Z. Wu, W. Guo, Y. Li, Y. Liu, and Q. Zhang, “High-speed and high-efficiency three-dimensional shape measurement based on Gray-coded light,” Photonics Res. 8(6), 819–829 (2020). [CrossRef]  

19. V. Srinivasan, H. C. Liu, and M. Halioua, “Automated phase-measuring profilometry of 3-D diffuse objects,” Appl. Opt. 23(18), 3105–3108 (1984). [CrossRef]  

20. M. Takeda and K. Mutoh, “Fourier transform profilometry for the automatic measurement of 3-D object shapes,” Appl. Opt. 22(24), 3977–3982 (1983). [CrossRef]  

21. Q. Kemao, “Two-dimensional windowed Fourier transform for fringe pattern analysis: principles, applications and implementations,” Opt. Laser Eng. 45(2), 304–317 (2007). [CrossRef]  

22. X. Su, L. Su, W. Li, and L. Xiang, “New 3D profilometry based on modulation measurement,” Proc. SPIE 3558, 1–7 (1998). [CrossRef]  

23. X. Su, L. Su, and W. Li, “New Fourier transform profilometry based on modulation measurement,” Proc. SPIE 3749, 438–439 (1999). [CrossRef]  

24. J. Li, X. Su, and L. Guo, “Improved Fourier transform profilometry for the automatic measurement of three-dimensional object shapes,” Opt. Eng. 29(12), 1439–1444 (1990). [CrossRef]  

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

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. Illustration of typical phase unwrapping errors in FPP.
Fig. 2.
Fig. 2. Framework of the proposed generalized phase unwrapping method.
Fig. 3.
Fig. 3. Principle of the proposed generalized Tri-PU phase unwrapping method. (a) Schematic diagram of region division. (b) Process of tripartite phase unwrapping.
Fig. 4.
Fig. 4. Algorithm flow chart and the MATLAB codes of generalized Tri-PU method.
Fig. 5.
Fig. 5. Experiments on a 1-bit defocusing projecting measuring system. (a) One of the captured deformed fringe pattern. (b) Original wrapped phase. (c)-(d) Staggered wrapped phases. (e) Fringe order. (f) Reconstructed absolute phase using the traditional phase unwrapping method. (g) Reconstructed absolute phase using the proposed generalized Tri-PU method. (h) Divided tripartite regions of fringe order. (i) Enlarged subimage in (f). (j) Enlarged subimage in (g).
Fig. 6.
Fig. 6. Comparative experiments using pre-avoided and post-corrected methods. (a) Reconstructed result using the traditional phase unwrapping method. (b)-(d) Filtering results using a filtering window with increasing size. (e) Reconstructed result using the generalized Tri-PU method. (f)-(h) Corresponding corrected results using post-corrected method.
Fig. 7.
Fig. 7. Experimental results by different phase calculation and unwrapping methods using a 8-bit focusing projecting measuring system. (a) Captured coding patterns. (b) Captured fringe patterns. (c) Wrapped phase. (d) Reconstructed height using the traditional phase unwrapping method. (e) Reconstructed height using the proposed generalized Tri-PU method. (f) Enlarged subimages of the details in (d) and (e).
Fig. 8.
Fig. 8. Critical condition of generalized Tri-PU method to ensure successful phase unwrapping.

Equations (9)

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

I i ( x , y ) = a ( x , y ) + b ( x , y ) cos [ Φ ( x , y ) δ i ] , i = 1 , 2 , 3 N
ϕ ( x , y ) = artan [ i = 1 N ( I i ( x , y ) sin δ i ) i = 1 N ( I i ( x , y ) cos δ i ) ] .
k m ( x , y ) = k ( x , y ) ,   w h e r e   | ϕ ( x , y ) | < π / 3.
k l ( x , y ) = k ( x , y ) ,   w h e r e   y < y min | ϕ ( x ) | ,   ( x , y ) A ( i )   a n d   ( x , y ) k m
k r ( x , y ) = k ( x , y ) ,   w h e r e   y > y min | ϕ ( x ) | ,   ( x , y ) A ( i )   a n d   ( x , y ) k m
{ ϕ 1 ( x , y )  =  w r a p [ ϕ ( x , y ) + 2 π / 3 ] ϕ 2 ( x , y ) = ϕ ( x , y )   ϕ 3 ( x , y )  =  w r a p [ ϕ ( x , y ) 2 π / 3 ] ,
w r a p [ ϕ ] = { ϕ 2 π ,   ϕ > π ϕ + 2 π ,   ϕ π .
Φ ( x , y ) = { ϕ 1 ( x , y ) + 2 π k l ( x , y ) 2 π / 3 ϕ 2 ( x , y ) + 2 π k m ( x , y )   ϕ 3 ( x , y ) + 2 π k r ( x , y ) + 2 π / 3 .
max [ Δ o r d e r + D ( Δ ϕ ) T ] < 1 3
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.