On window size selection in the windowed Fourier ridges algorithm: Addendum

The windowed Fourier transform is a useful technique for fringe pattern analysis. It has been shown that the proper selection of the window size is a balance between the linear phase approximation error and the influence of noise [Q. Kemao, On widow size selection in windowed Fourier ridges algorithm. Opt Lasers Eng, 2007, accepted for publication]. Since the fringe intensity and noise level usually vary spatially, the window size should also be spatially adaptive in order to reach a good balance for each pixel of the fringe pattern. This addendum first shows that the windowed Fourier transform with a spatially fixed window size (SFWS) is still practically useful and then discusses the window size competition strategies for the windowed Fourier transform with a spatially adaptive window size (SAWS). The windowed Fourier transform with a SAWS is theoretically better than that with a SFWS but it is also more challenging in use. The windowed Fourier ridges algorithm is used for analysis throughout this paper. This analysis is also applicable to the windowed Fourier filtering algorithm.     

Fringe pattern analysis by S-transform

S-transform proposed in 1996 by Stockwell R.G is a simple and popular technique for the time–frequency analysis. It has been introduced in optical three-dimensional shape measurement, recently. In this paper, a study about applications of S-transform in the demodulation of deformed fringe patterns is performed. We focus on discussing not only the S-transform spectrum filtering technique, the S-transform ridge technique and the phase gradient calculation method based on S-transform used in fringe pattern demodulation, but also the phase unwrapping technique. In addition, a generalized S-transform was introduced to analyze fringe patterns, which is helpful to improve the measurement accuracy and flexibility of the method based on S-transform. The reconstruction results based on S-transform were compared with that on wavelet transform and windowed Fourier transform in fringe analysis.     

An improved windowed Fourier transform for fringe demodulation

A modified algorithm of windowed Fourier transform (WFT) for phase retrieval from electronic speckle-shearing fringe patterns with carriers is proposed. The algorithm is based on the introduction of a fast Fourier transform (FFT) in WFT to reduce computation time for fringe demodulation. Since boundary effects in FFT will influence the accuracy of phase retrieval, the Gerchberg method is employed to extrapolate the fringe pattern at the boundaries to reduce boundary effects. A theoretical analysis of the algorithm is presented. Both simulated and experimental results show that the proposed method has reduced the computation time significantly compared with the convolution method of WFT without sacrificing measurement accuracy.     

Two-dimensional multiscale windowed Fourier transform based on two-dimensional wavelet transform for fringe pattern demodulation

A two-dimensional multiscale windowed Fourier transform (2D-MWFT), based on two-dimensional Gabor wavelet transform (2D-GWT), for the phase extraction from a spatial fringe pattern in fringe projection profilometry is presented. First, the instantaneous frequencies on x and y direction of the modulated fringe pattern are determined by 2D-GWT, and then the local stationary lengths are obtained. The 2D-MWFT with different two-dimensional Gaussian windows whose width is set according to the local stationary length is preformed for each section of the modulated fringe pattern to achieve multiresolution analysis and phase demodulation. Comparing the result of the phase demodulated by 2D-GWT and two-dimensional windowed Fourier transform (2D-WFT) with that by 2D-MWFT in a numerical simulation, we show that the 2D-MWFT method is superior to these methods, especially for the local non-stationary signal with low frequency. The theory and the results of a simulation and experiment are shown.     

Method for eliminating zero spectrum in Fourier transform profilometry

Windowed Fourier transform (FT) is introduced into Fourier transform profilometry (FTP) to extract zero component included in deformed fringe pattern without capturing two fringe patterns with π phase difference. First, we extract zero component of deformed fringe image using windowed FT and normalize it, then calculate the normalized Fourier spectra of deformed fringe image, finally subtract the zero component from Fourier spectra of deformed fringe image. We give the theoretical analysis of the proposed method. The experiment is carried out to demonstrate that the method can eliminate zero component and extend the measurable slope of height variation to nearly three times of the conventional FTP.     

三维测量中一种新的自适应窗口傅里叶相位提取法

针对多尺度窗口傅里叶变换中,窗口尺寸的自适应选取及提取基频时的频谱混叠等问题,提出基于希尔伯特-黄变换(HHT)的自适应窗口傅里叶相位提取法。对变形条纹信号进行HHT后,通过谱分析,自适应确定能够准确描述条纹信号变化情况的瞬时频率及条纹图的背景分量。根据所得的瞬时频率,给出自适应定位条纹信号局部平稳区域的步骤,进而确定窗口尺寸。不需额外计算,可有效去除背景分量以减少基频提取过程中零频频谱的干扰。与现有的用最大脊法确定窗口尺寸的方法相比,本方法不受被测相位必须线性逼近且变化缓慢的前提约束。实验证明本方法有效、可行,且对测量携带陡峭边缘或面形复杂的物体也能进行较为精确有效的测量。     

Two-dimensional windowed Fourier transform for fringe pattern analysis: Principles, applications and implementations

Fringe patterns from optical metrology systems need to be demodulated to get the desired parameters. Two-dimensional windowed Fourier transform is chosen for the determination of phase and phase derivatives. Two algorithms, one based on filtering and the other based on similarity measure, are developed. Some applications based on these two algorithms are explored, including strain determination, phase unwrapping, phase-shifter calibration, fault detection, edge detection and fringe segmentation. Various examples are given to demonstrate the ideas. Finally implementations of these algorithms are addressed. Most of the work has appeared in various papers and its originality is not claimed. Instead, this paper gives an overview and more insights of our work on windowed Fourier transform.     

Generalized three-dimensional windowed Fourier transform for fringe analysis

A 3D windowed Fourier transform is proposed for fringe sequence analysis, which processes the joint spatial and temporal information of the fringe sequence simultaneously. The 2D windowed Fourier transform in the spatial domain and the 1D windowed Fourier transform in the temporal domain are two special cases of the proposed method. The principles of windowed Fourier filtering and windowed Fourier ridges are developed. Experimental verification shows encouraging results despite a longer processing time.     

基于Gabor滤波的散斑条纹图平滑方法

提出了一种基于伽博(Gabor)滤波的散斑条纹图平滑方法.通过加窗傅里叶运算提取散斑条纹图的条纹频率和条纹梯度方向,并利用它们确定具有频率和方向选择性的Gabor滤波器对散斑条纹图进行Gabor滤波.数值模拟和实验结果表明,该方法在滤除散斑噪声的同时能够有效的保留散斑条纹图的条纹结构,为从单幅散斑条纹图中进一步提取条纹相位场奠定了良好的基础.     

窗口傅里叶变换轮廓术中窗口尺度选取的改进

多尺度窗口傅里叶变换法根据条纹信号的瞬时频率梯度来确定信号的局域平稳长度,再由局域平稳长度来控制窗口的尺度,即窗口的尺度和局域平稳长度成正比。使用多尺度窗口傅里叶变换法来使条纹信号的频谱局域化,可以在条纹信号的频率分辨率和空间分辨率之间达到一种较佳的调和。针对多尺度窗口傅里叶变换三维形貌测量技术中局域平稳长度提取算法的不足进行了改进,使窗口尺度的选取更具合理性,对变形光栅基频提取更精确,进一步减小了位相测量的误差。给出了理论分析、计算机模拟以及实验结果。     

On window size selection in the windowed Fourier ridges algorithm

The windowed Fourier ridges (WFR) algorithm can be used to extract the frequency and phase information from an exponential or carrier fringe pattern. The selection of the window size for this algorithm is investigated in this paper. For exponential phase fringe patterns, the window size does not need to be adjusted according to the signal frequency. Hence the scaling strategy of the wavelet transform is not necessary. Instead, the window size should be selected according to the balance between the linear phase approximation error and the noise level. For carrier fringe patterns, the influence of the conjugate component should also be considered. However it can be ignored when the lowest frequency of the fringe pattern is higher than 1.5/σ (σ is used to represent the window size).     

Parallel computing for fringe pattern processing: A multicore CPU approach in MATLAB environment

In the process of measurements such as optical interferometry and fringe projection, an important stage is fringe pattern analysis. Many advanced fringe analysis algorithms have been proposed including regularized phase tracking (RPT), partial differential equation based methods, wavelet transform, Wigner–Ville distribution, and windowed Fourier transform. However, most of those algorithms are computationally expensive. MATLAB^® is a general algorithm development environment with powerful image processing and other supporting toolboxes. It is also commonly used in photomechanical data analysis. With rapid development of multicore CPU technique, using multicore computer and MATLAB^® is an intuitive and simple way to speed up the algorithms for fringe pattern analysis. The paper introduces two acceleration approaches for fringe pattern processing. The first approach is task parallelism using multicore computer and MATLAB^® parallel computing toolbox. Since some algorithms are embarrassing problems, our first approach makes use of this characteristic to parallelize these algorithms. For this approach, parallelized windowed Fourier filtering (WFF) algorithm serves as an example to show how parallel computing toolbox accelerates the algorithm. Second, data parallelism using multicore computer and MATLAB^® parallel computing toolbox is proposed. A high level parallel wrapping structure is designed, which can be used for speeding up any local processing algorithms. WFF, windowed Fourier ridges (WFR), and median filter are used as examples to illustrate the speedup. At last, the results show that the parallel versions of former sequential algorithm with simple modifications achieve the speedup up to 6.6 times.     

基于S变换的改进窗口傅里叶三维测量法

讨论了一维S变换在处理条纹图中的应用,推导了S变换表达式,提出一种基于S变换脊的自适应窗口傅里叶三维测量法。利用S变换对于频率的敏感性以及S变换脊思想,求得任意点的瞬时频率。推导了在相位跳变剧烈区域相位二阶导φ″(b)对S变换脊的影响,并求得去除此影响后更为准确的瞬时频率。然后将瞬时频率倒数作为经典窗口傅里叶变换(WFT)的窗口大小,解出条纹图的精确相位。最后对去除相位二阶导影响前后基于S脊的WFT进行对比仿真;实验验证其与传统基于小波脊的WFT在抗噪声、解相位精度上的优势。     

Phase retrieval of optical fringe patterns from the ridge of a wavelet transform

A new method for phase retrieval of optical fringe patterns is presented. This method is based on a wavelet transform and is capable of extracting the full 2D phase distribution from a single fringe pattern. An important conclusion that the phase of the optical fringe pattern is equal to the phase of its wavelet transform on the ridge of the wavelet transform is theoretically clarified. The method is compared with the Fourier transform and the integration methods. A numerical simulation and an experimental example of phase retrieval are shown.     

Windowed Fourier transform profilometry based on improved S-transform

A novel method for windowed Fourier transform (WFT) profilometry is presented. This method is based on improved S-transform. The impact of the second order derivative of the phase (φ^'(b)) to the ridge of S-transform is derived, and how to estimate this deviation is discussed. An important conclusion that more accurate instantaneous frequency can be obtained after removing this deviation is shown. Thus, an accurate phase map of the fringe pattern is obtained by using the WFT based on the window size map, and this map is related to the instantaneous frequency. The method is compared with the WFT based on the wavelet transform. A numerical simulation and experimental example have shown its validity in practical applications.     

Single fringe projection profilometry based on sinusoidal intensity normalization and subpixel fitting

A novel fringe projection profilometry using a single sinusoidal fringe pattern projected is proposed. Computer-generated sinusoidal fringe and uniform intensity patterns are firstly projected on a testing object by a liquid crystal display projector. The variable reflection intensity of a fringe pattern is then roughly normalized by division operation applied to the grabbed fringe and uniform intensity patterns projected. Fringe intensity is further normalized by employing an interpolation algorithm. The deformed sinusoidal pattern encoding object shape is converted to a wrapped phase map without using phase-shifting or Fourier transform. Computer simulation and experimental performance are evaluated to demonstrate the validity of the proposed method. The experimental results compared with those of the four-step phase-shifting and fast Fourier transform methods are also presented.     

基于窗口傅里叶变换剪切干涉法波前检测

提出了一种利用二维窗口傅里叶变换从径向剪切干涉条纹中准确得到波前的重建技术。首先对剪切干涉条纹做二维窗口傅里叶变换,设置阈值和频率积分范围后,进行二维窗口傅里叶逆变换,然后对包裹相位做去载频和相位展开处理得到相位差分布,最后使用波前迭代算法从相位差中复原实际波前。模拟计算表明,使用该方法最大相位复原误差为0.82%,均方根值为0.020 9 rad,实验结果验证了该方法的有效性。同时也对窗口傅里叶变换的关键参数,如窗函数的选择、窗口大小的确定以及阈值的选取等进行了简要讨论。与传统傅里叶变换法(FFT)相比,基于窗口傅里叶变换的剪切干涉波前检测法有更高的精度和稳定性,为波前检测提供一种新的处理方法。     

Computer aided evaluation of fringe patterns

The crucial step in computer aided evaluation of interferometric fringe patterns is the determination of the interference phase distribution from the recorded and stored intensity pattern. Methods for determination of the interference phase distributions from fringe patterns are presented and the most important methods, which are skeletonising, temporary heterodyning, phase shifting and Fourier transform evaluation are compared with regard to experimental requirements, achievable resolution and precision, as well as inherent noise suppression and image enhancement. The comparison shows that whenever phase shifting is possible, it is the best choice. If only one interference pattern is offered, Fourier transform evaluation is recommendable. For demodulation of the wrapped phase a path-independent algorithm is presented. Practical examples are given from holographic interferometric measurements.     

Phase-shifting windowed Fourier ridges for determination of phase derivatives

Determination of the phase or phase derivative from interferometric fringe patterns is an important task in optical interferometry. The use of wavelet ridges was recently shown to be an effective method for phase retrieval from a single fringe pattern. One necessary requirement in this method is the need for carrier frequency. In cases when carrier frequency is not available, the novel phase-shifting windowed Fourier ridges method can be used. Phase derivatives with the proper sign can be directly retrieved even in the presence of noise. An application for curvature determination from speckle shearographic fringes demonstrates the effectiveness of the method.