Multi-STFT for adaptive SCPS phase determination

This paper presents a fast and reliable approach for phase modulo 2π-calculation from a single fringe pattern. It calculates correct phase values even for very complex and variable shape gradients based on a locally variable fringe period determined for the entire image. In the paper, a new two-step method for wrapped phase calculation is proposed. It is performed through the use of a method based on a multiple local fast Fourier transform for estimation of a local fringes period map and a 5-point spatial carrier phase shifting (SCPS) formula for phase modulo 2π-calculation. The described approach is verified by a correct demodulation of a real fringe pattern taken by a 3D-shape measurement system.     

A derivative based simplified phase tracker for a single fringe pattern demodulation

In this paper, a novel fringe demodulation method for the estimation of phase and its first-order derivative from a closed-fringe interferogram is proposed. The proposed method determines the phase derivatives in both x&y directions from fringe orientation and density. The phase derivatives are subsequently used to determine phase values using a novel simplified phase tracker. In the phase tracking model, the complexity of the cost function is reduced using predetermined derivatives so computation time required for phase tracking is reduced considerably. The proposed model is more robust while dealing with saddle points in fringes than the conventional phase tracker model. Hence it does not require any specialized scanning strategy. The proposed method is validated with simulated and experimental fringe patterns (obtained using electronic speckle pattern interferometry and optical holographic interferometry) and a comparison study is carried out with conventional regularized phase tracker. The simulation results show that the proposed method has good accuracy and requires less computation time than existing phase-tracking algorithms. The experimental results demonstrate the robustness of the proposed method against speckle noise and its practical applicability for static and dynamic applications.     

An improved stair phase encoding method for absolute phase retrieval

An improved phase unwrapping method is proposed to reduce the projection fringes in three-dimensional (3D) surface measurement. Color fringe patterns are generated by encoding with sinusoidal fringe and stair phase fringe patterns in red and blue channels. These color fringe patterns are projected onto the tested objects and then captured by a color CCD camera. The recorded fringe patterns are separated into their RGB components. Two groups of four-step phase-shifting fringe patterns are obtained. One group of the stripes are four sinusoidal patterns, which are used to determine the wrapped phase. The other group of stripes are four sinusoidal patterns with the codeword embedded into stair phase, whose stair changes are perfectly aligned with the 2π discontinuities of sinusoidal fringe phase, which are used to determine the fringe order for the phase unwrapping. The experimental results are analyzed and compared with those of the method in Zheng and Da (2012. Opt Express 20(22):24139–24150). The results show that the proposed method needs only four fringe patterns while having less error. It can effectively reduce the number of projection fringes and improve the measuring speed.     

基于格雷码的分区间相位展开方法

格雷码因具有良好的鲁棒性和抗噪性,被广泛应用到结构光投影三维成像方法中。在三维测量过程中,由于设备以及其他环境噪声的影响,格雷码方法解码条纹级次边沿和截断相位边沿通常无法处于理想的对齐状况,使得展开的相位出现跳变现象。为了更好地避免级次跳变误差,使得边沿跳变区域的容错宽度更大,提出基于格雷码的分区间相位展开方法。在互补格雷码基础上增加一幅格雷码图像,利用所有格雷码解得附加码字,通过对附加码字进行不同位移量的条纹级次映射,得到2个辅助条纹级次。利用所有条纹级次,对截断相位进行分区间相位展开,在边沿跳变区域错误大于半个周期时,仍能获取到无跳变现象的展开相位。实验结果表明,当边沿错误区域宽度小于3/4个条纹周期宽度时,可以有效避免级次跳变产生的误差。     

Fractional derivative applied to the retrieval of phase distribution from an interferogram

This paper proposes a new method for determining the phase distribution of a 2D fringe pattern using 2D fractional derivatives. The fractional derivative of order?r (real 0??r??3), in the ?? direction (????[0,2??]) is defined by the Fourier transform and its inverse and gives as result a phase distribution. The fractional derivative of order r produces a phase displacement of r??/2 without significantly changes in normalized amplitude distribution. There are optimum values of r and ?? for which the phase distribution given by the fractional derivative represents the phase distribution for the fringe pattern. Also, we demonstrate the usefulness of a method based on four fractional derivatives to determine the phase and intensity distribution from a single interferogram and highlights small objects located in the fringe pattern. We demonstrate the proposed method using both numerical simulated and experimental interferograms.     

一种三维数字成像系统的多视点姿态估计方法

为校准多视场深度数据,提出基于条纹投影的三维数字成像系统的多视点姿态估计方法。该方法至少在两个视点分别向被测物体投射出一组正交条纹图,利用条纹投影和相位重建技术,将相位图映射为物体的三维空间坐标。进而,利用投影仪的投射过程是摄像机成像过程的逆过程,建立投影仪的投射平面和摄像机的成像平面的对应关系,将“极线几何约束”应用到基于条纹投影的主动三维视觉的姿态估计问题,并在考虑测量数据受噪声影响的条件下,建立了求解视点姿态参量的数学模型。通过优化求解非线性方程可以获得多视点的姿态估计参量。所设计的实验及结果证明了所提出方法的有效性。     

Blind phase error suppression for color-encoded digital fringe projection profilometry

Color-encoded digital fringe projection profilometry (CDFPP) has the advantage of fast speed, non-contact and full-field testing. It is one of the most important dynamic three-dimensional (3D) profile measurement techniques. However, due to factors such as color cross-talk and gamma distortion of electro-optical devices, phase errors arise when conventional phase-shifting algorithms with fixed phase shift values are utilized to retrieve phases. In this paper, a simple and effective blind phase error suppression approach based on isotropic n-dimensional fringe pattern normalization (INFPN) and carrier squeezing interferometry (CSI) is proposed. It does not require pre-calibration for the gamma and color-coupling coefficients or the phase shift values. Simulation and experimental works show that our proposed approach is able to effectively suppress phase errors and achieve accurate measurement results in CDFPP.     

Analysis on fringe pattern demodulation by use of 2-D CWT

We discuss the wavelet transform profilometry based on the continuous wavelet transform technique as viewed from frequency analysis. We deduce the expression of one-dimensional (1-D) and two-dimensional (2-D) wavelet transform in frequency domain and analyze their characteristics in the application of demodulating the fringe patterns. We also compare 1-D CWT and 2-D CWT in demodulating the oblique fringe patterns with dual carrier frequency components. When oblique fringe patterns are processed, the direction normal to the grating line and x axis or y axis is not identical. By 1-D CWT, in which wavelet transform is carried out row by row, we cannot obtain the most similarity between local signal and the wavelet functions with different dilation values. While a fan 2-D continuous wavelet transformation can deal with the fringe pattern as a 2-D unit as well as has multi-directions, its advantage is that it can be used to exact the information in the spatial direction. However, its spatial localization ability is not very good, which leads that it is not suitable for demodulating the fringe patterns with high phase variation. Computer simulations and experiments have verified our analysis.     

High-resolution frequency estimation technique for recovering phase distribution in interferometers

An integral approach to phase measurement is presented. First, the use of a high-resolution technique for the pixelwise detection of phase steps is proposed. Next, the robustness of the algorithm that is developed is improved by incorporation of a denoising procedure during spectral estimation. The pixelwise knowledge of phase steps is then applied to the Vandermonde system of equations for retrieval of phase values at each pixel point. Conceptually, our proposal involves the design of an annihilating filter that has zeros at the frequencies associated with the polynomial that describes the fringe intensity. The parametric estimation of this annihilating filter yields the desired spectral information embedded in the signal, which in our case represents the phase steps. The proposed method offers the advantage of extracting the interference phase of nonsinusoidal waveforms in the presence of miscalibration error of the piezoelectric transducer. In addition, in contrast to previously reported methods, this method does not require the application of selective phase steps between data frames for nonsinusoidal waveforms.     

Patch-wise denoising of phase fringe patterns based on matrix enhancement

We propose a new approach for the denoising of a phase fringe pattern recorded in an optical interferometric setup. The phase fringe pattern which is generally corrupted by high amount of speckle noise is first converted into an exponential phase field. This phase field is divided into a number of overlapping patches. Owing to the small size of each patch, the presence of a simple structure of the interference phase is assumed in it. Accordingly, the singular value decomposition (SVD) of the patch allows us to separate the signal and noise components effectively. The patch is reconstructed only with the signal component. In order to further improve the robustness of the proposed method, an enhanced data matrix is generated using the patch and the SVD of this enhanced matrix is computed. The matrix enhancement results in an increased dimension of the noise subspace which thus accommodates more amount of noise component. Reassignment of the filtered pixels of the preceding patch in the current patch improves the noise filtering accuracy. The fringe denoising capability in function of the noise level and the patch size is studied. Simulation and experimental results are provided to demonstrate the practical applicability of the proposed method.     

一种基于傅里叶变换的分析载波条纹的新方法

针对传统傅里叶变换法处理光载波干涉条纹图时会有边缘效应产生的问题,提出了一种基于傅里叶变换法的外推延拓方法,并从理论上进行了数学推导。为了验证这种方法的正确性,分别对一维数字信号和二维空间载波条纹图进行了数值模拟,进一步分析了误差产生的原因,并与传统的傅里叶变换法对比。结果表明该法可以有效抑制传统傅里叶变换法处理光载波干涉条纹图时边缘效应所造成的较大误差,在基于空间域相位调制技术的波面干涉测量中,对空间载波条纹图进行处理,可以使相位的计算精度达到3.3 mrad。     

基于希尔伯特变换的干涉条纹相位解调新算法

提出了一种新的应用希尔伯特变换解调干涉条纹相位的算法,可以从单幅干涉条纹图中解调出全场相位分布.在实际应用中,常借助傅里叶变换实现希尔伯特变换算法,但是会忽略负频率成分,造成相位信息的丢失.对于相位分布非单调变化的干涉条纹,提出了一种判断函数,用来计算相位信息零频率点的分布.利用相位的零频率点分布构造了一个二元模板,使用该模板对本文提出的两次希尔伯特变换法产生的包裹相位图进行修正.对修正后的包裹相位图进行解包裹处理,可以得到连续的全场相位分布.对该方法用计算机模拟进行了验证.     

Multi-frequency inverse-phase fringe projection profilometry for nonlinear phase error compensation

A new multi-frequency inverse-phase method was proposed to compensate for nonlinear phase errors in fringe projection profilometry and to measure the three-dimensional shape of discontinuous objects. After introducing a phase offset of π/4 into the multi-frequency four-step phase-shifting method the corresponding nonlinear phase error reversed its sign, which allowed the addition of unwrapped phases before and after the phase-offset operation to compensate for the error. For the four-step phase-shifting method, simulation analysis showed that the nonlinear phase error had quadrupled the fringe frequency. Moreover, experimental results verified the feasibility and applicability of the proposed method.     

散斑干涉相位条纹图的频域滤波处理

为了解决在数字散斑干涉技术测量时,散斑干涉相位条纹图像中大量噪声对相位解包裹结果和精度产生严重影响的问题,介绍了一种条纹正余弦分解和频域低通滤波结合的方法,实现了散斑干涉相位条纹图的高精度滤波。该方法的基本思路是在对相位图像进行滤波处理前,先将相位图通过正余弦函数进行映射转换成两幅图,分别经过频域滤波,然后再合成为相位图。这种分解频域滤波方法可以在滤波的同时,有效保留相位跳变信息。实验结果表明：与传统的图像降噪方法相比,该方法能够在保留图像“尖峰”信息的基础上,较好地滤除图像中的散斑噪声,方法简单有效,有效解决了传统滤波方法应用在相位条纹图中,相图灰度信息丢失10%~40%的问题。     

Predictive flow-field estimation

We are considering the problem of real-time prediction of 3D turbulent velocity fields based on a small number of scalar measurements. The method of proper orthogonal decomposition (POD) allows for the decomposition of an ensemble of velocity fields into a set of spatial basis functions and a set of temporal coefficients. The computation of the temporal coefficients is by no means a trivial matter, especially when one is faced with a large number of modes. In this paper we discuss the use of radial basis function (RBF) models to capture the discrete time evolution and nonlinear dynamics of the POD coefficients. Further, we propose the use of regularized regression techniques to generate models that provide mappings between the POD coefficients and scalar measurements. As a final step towards real-time prediction, the state-space RBF models and regression measurement models are combined using unscented Kalman filters to produce optimal solutions such that a balance between the state models and measurement models is achieved.The proposed methods are tested for two specific cases. The classical Lorenz model is chosen to demonstrate the use and effectiveness of RBF models as a potential candidate for state models. Flow around a wall-mounted cube in a channel at Re=20,000 is considered as the second case. The aim for the second case is to be able to accurately predict the POD coefficients outside the ensemble. It is shown that a large number of POD coefficients is required to approximate the velocity fields with sufficient accuracy. The RBF models are created based on only the temporal information available from the initial ensemble, and it is shown that the RBF model is able to correctly approximate the high-dimensional phase space. Combined with the unscented Kalman filter it is indeed possible to track the evolution of the POD coefficients for a long time. The robustness of the filter is demonstrated by considering the presence of noise in measurements and using measurement information at time steps greater than the evolution time step.     

Absolute phase recovery of three fringe patterns with selected spatial frequencies

A new temporal approach is presented for the recovery of the absolute phase maps from their wrapped versions based on the use of fringe patterns of three different spatial frequencies. In contrast to the two-frequency method recently published, the method proposed is characterized by better anti-error capability as measured by phase error tolerance bound. A general rule for the selection of the three frequencies is presented, and its relationship to the phase error tolerance bound is derived. Theoretical analysis and experimental results are also presented to validate the effectiveness of the proposed three frequency technique.     

Shifted-phase calibration for a 3-D shape measurement system based on orthogonal composite grating projection

This paper proposes a novel method for reducing measurement error caused by spectrum overlapping in orthogonal-composite-grating-based 3-D measurement method. For 3-D measurement systems based on orthogonal composite grating projection, spectrum overlapping causes phase of each deformed phase-shifting fringe changed differently, which violates the principle that the shifted phases between adjacent deformed fringes must be equivalent to 2π/3, and therefore results in phase measurement error. The proposed shifted-phase calibration method is based on that phase variation of each deformed fringe is independent of height and reflectivity of the measured object. Three composite gratings are projected on the reference plane, and each carrier channel includes three phase-shifting gratings needed in phase measuring profilometry (PMP). Because the adjacent phase-shifting fringes demodulated from the same carrier channel have the phase difference of 2π/3, we can respectively calculate the reference plane's phases of three carrier channels by the phase algorithm of PMP method, and the shifted phases between them are obtained. When an object is measured, the shifted phases between deformed phase-shifting fringes can be calibrated. A new 3-D measurement mathematical model is set to reconstruct object. Our experiments prove that the proposed method can effectively restrain the effect of spectrum overlapping and improve measurement accuracy almost one times.     

抑制高阶谐波误差的随机相移条纹分析算法

高阶谐波和随机相移误差是影响条纹分析精度的主要因素。为了同时解决这两个问题,提出了基于频域滤波的迭代相移算法。该算法采用巴特沃斯低通滤波器,从频域上滤除条纹的高阶谐波分量,再运用最小二乘迭代方法从三帧随机相移条纹图像中提取相位信息。数值模拟和实验结果表明,该算法可有效地抑制由高阶谐波和随机相移引入的波纹误差,误差PV值和RMS值分别为0.368 8 rad和0.025 3 rad,其精度高于传统的三步相移算法和Wang算法。该方法适合于高精度干涉测量和三维物体表面轮廓测量。     

Robust wrapping-free phase retrieval method based on weighted least squares method

Phase unwrapping is one of the most challenging processes in many profilometry techniques. To sidestep the phase unwrapping process, Perciante et al. (2015) proposed a wrapping-free method to retrieve the phase based on the direct integration of the spatial derivatives of the fringe patterns. However, this method is only applicable to objects with phase continuity, so it may fail to handle fringe patterns containing complicated singularities such as noise, shadows, shears and surface discontinuities. In the light of this problems, a robust wrapping-free phase retrieval method is proposed that is based on the combined use of Perciante's method and the weighted least squares method. Two partial derivatives of the desired phase are obtained from the fringe patterns, while the carrier is eliminated using the direct phase difference method. The phase singularities are determined using a derivative variance correlation map (DVCM), and the weighting coefficient is obtained from the binary mask of the reverse DVCM. Simulations and experiments are conducted to prove the validity of the proposed method. The results are analyzed and compared with those of Perciante's method demonstrating that in addition to maintaining the advantage of sidestepping the phase unwrapping process, the proposed method is available for measuring objects with some types of singularities sources.     

Signal processing method of phase correction for laser heterodyne interferometry

A novel signal processing method of movement direction identification and phase correction is presented for laser heterodyne interferometry. Based on the reference signal, four intervals with phase difference of 90° each other are set up. The real-time movement direction identification and the integer fringe counting are realized by detecting the times that the rising-edge of the measurement signal crosses the intervals. The phase correction approach is proposed in detail to solve the fraction phase compensation when the initial phase difference is not equal to the zero phase difference. Three experiments of the stability test, the nanometer and micrometer displacement tests on bi-directional movement were performed to demonstrate the usefulness and feasibility of the presented signal processing method.