Adaptive monogenic filtering and normalization of ESPI fringe patterns

Fringe pattern demodulation using Zernike polynomials and a l1-norm regularized extended Kalman filter

A novel algorithm for closed fringe demodulation for an absolute phase estimation, to the best of our knowledge, is proposed. The two-dimensional phase is represented as a weighted linear combination of a certain number of Zernike polynomials (ZPs). Essentially, the problem of phase estimation is converted into the estimation of ZP coefficients. The task of ZP coefficient estimation is performed based on a state space model. Due to the nonlinear dependence of the fringe intensity measurement model on the ZP coefficients, the extended Kalman filter (EKF) is used for the state estimation. A pseudo-measurement model is considered based on the state vector sparsity constraint to improve the convergence performance of the EKF. Simulation and experimental results are provided to demonstrate the noise robustness and the practical applicability of the proposed method.

Fringe analysis: single-shot or two-frames? Quantitative phase imaging answers

Conditions of the digital recording of the fringe pattern determine the phase reconstruction procedure, which in turn directly shapes the final accuracy and throughput of the full-field (non-scanning) optical measurement technique and defines the system capabilities. In this way, the fringe pattern analysis plays a crucial role in the ubiquitous optical measurements and thus is under constant development focused on high temporal/spatial resolution. It is especially valuable in the quantitative phase imaging technology, which emerged in the high-contrast label-free biomedical microscopy. In this paper, I apply recently blossomed two-frame phase-shifting techniques to the QPI and merge them with advanced adaptive interferogram pre-filtering algorithms. Next, I comprehensively test such frameworks against classical and adaptive single-shot methods applied for phase reconstruction in dynamic QPI enabling highest phase time-space-bandwidth product. The presented study systematically tackles important question: what is the gain, if any, in QPI realized by recording two phase-shifted interferograms? Counterintuitively, the results show that single-shot demodulation exhibited higher phase reconstruction accuracy than two-frame phase-shifting methods in low and medium interferogram signal-to-noise ratio regimes. Thus, the single-shot approach is promoted due to not only high temporal resolution but also larger phase-information throughput. Additionally, in the majority of scenarios, the best option is to shift the paradigm and employ two-frame pre-filtering rather than two-frame phase retrieval. Experimental fringe analysis in QPI of LSEC/RWPE cell lines successfully corroborated all novel numerical findings. Hence, the presented numerical-experimental research advances the important field of fringe analysis solutions for optical full-field measurement methods with widespread bio-engineering applications.

Automatic fringe pattern enhancement using truly adaptive period-guided bidimensional empirical mode decomposition

Fringe patterns encode the information about the result of a measurement performed via widely used optical full-field testing methods, e.g., interferometry, digital holographic microscopy, moiré techniques, structured illumination etc. Affected by the optical setup, changing environment and the sample itself fringe patterns are often corrupted with substantial noise, strong and uneven background illumination and exhibit low contrast. Fringe pattern enhancement, i.e., noise minimization and background term removal, at the pre-processing stage prior to the phase map calculation (for the measurement result decoding) is therefore essential to minimize the jeopardizing effect the mentioned error sources have on the optical measurement outcome. In this contribution we propose an automatic, robust and highly effective fringe pattern enhancement method based on the novel period-guided bidimensional empirical mode decomposition algorithm (PG-BEMD). The spatial distribution of the fringe period is estimated using the novel windowed approach and then serves as an indicator for the truly adaptive decomposition with the filter size locally adjusted to the fringe pattern density. In this way the fringe term is successfully extracted in a single (first) decomposition component alleviating the cumbersome mode mixing phenomenon and greatly simplifying the automatic signal reconstruction. Hence, the fringe term is dissected without the need for modes selection nor summation. The noise removal robustness is ensured employing the block matching 3D filtering of the fringe pattern prior to its decomposition. Performance validation against previously reported modified empirical mode decomposition techniques is provided using numerical simulations and experimental data verifying the versatility and effectiveness of the proposed approach.

Impact of background and modulation on parameter estimation using fractional Fourier transform and its solutions

Analysis of closed-fringe patterns with quadratic phase that are often encountered plays an important role in optical interferometry. But because the frequency spectra of the two exponential signals that compose closed-fringe patterns overlap in the Fourier domain while one is clearly distinct from the other in the fractional Fourier domain, fractional Fourier transform (FRFT) is a useful method for analyzing the images to provide parameter estimation. However, when the fringe pattern has varying background and/or modulation due to non-uniform illumination, parameter estimation accuracy based on FRFT is affected, which lacks theoretical justification. Thus, the impact of varying background and/or modulation on the FRFT is studied with theoretic analysis and presented in this paper. Key factors that contribute to the optimal results are discussed when employing three kinds of fringe normalization methods to eliminate the impact. Here, the fringe pattern is first processed by the normalization technique. Then the cosine-only term is used to estimate parameter by use of the FRFT-based method. Physical quantities are then obtained by parameter estimation. In comparison with our previous method based on FRFT, more accurate results are achieved. The feasibility and applicability of the proposed approach are demonstrated using simulation and experimental results.

Phase retrieval in two-shot phase-shifting interferometry based on phase shift estimation in a local mask

Fringe analysis in two-shot phase-shifting interferometry is important but meets challenges due to a limited number of images, corrupting noise, and background modulation. Here we propose an effective algorithm for phase retrieval from two interferograms with unknown phase shifts. The algorithm first evaluates the phase shift in a local mask through phase fitting and global optimization and then computes a full-field phase map using an arctangent function. Since the phase shift evaluation is performed within a local mask, the algorithm is fast compared with conventional optimization-based algorithms and typically needs tens of seconds to complete the processing. Computer simulation and experimental results show that the proposed algorithm has excellent performance compared with state-of-the-art algorithms. A complete software package of the algorithm in MATLAB is available at http://two-shot.sourceforge.io/.

Demodulation of two-shot fringe patterns with random phase shifts by use of orthogonal polynomials and global optimization

We propose a simple and robust phase demodulation algorithm for two-shot fringe patterns with random phase shifts. Based on a smoothness assumption, the phase to be recovered is decomposed into a linear combination of finite terms of orthogonal polynomials, and the expansion coefficients and the phase shift are exhaustively searched through global optimization. The technique is insensitive to noise or defects, and is capable of retrieving phase from low fringe-number (less than one) or low-frequency interferograms. It can also cope with interferograms with very small phase shifts. The retrieved phase is continuous and no further phase unwrapping process is required. The method is expected to be promising to process interferograms with regular fringes, which are common in optical shop testing. Computer simulation and experimental results are presented to demonstrate the performance of the algorithm.

Single shot fringe pattern phase demodulation using Hilbert-Huang transform aided by the principal component analysis

Hybrid single shot algorithm for accurate phase demodulation of complex fringe patterns is proposed. It employs empirical mode decomposition based adaptive fringe pattern enhancement (i.e., denoising, background removal and amplitude normalization) and subsequent boosted phase demodulation using 2D Hilbert spiral transform aided by the Principal Component Analysis method for novel, correct and accurate local fringe direction map calculation. Robustness to fringe pattern significant noise, uneven background and amplitude modulation as well as local fringe period and shape variations is corroborated by numerical simulations and experiments. Proposed automatic, adaptive, fast and comprehensive fringe analysis solution compares favorably with other previously reported techniques.

Theory and algorithms of an efficient fringe analysis technology for automatic measurement applications

Some advances in fringe analysis technology for phase computing are presented. A full scheme for phase evaluation, applicable to automatic applications, is proposed. The proposal consists of: a fringe-pattern normalization method, Fourier fringe-normalized analysis, generalized phase-shifting processing for inhomogeneous nonlinear phase shifts and spatiotemporal visibility, and a phase-unwrapping method by a rounding-least-squares approach. The theoretical principles of each algorithm are given. Numerical examples and an experimental evaluation are presented.

Two-shot fringe pattern phase-amplitude demodulation using Gram-Schmidt orthonormalization with Hilbert-Huang pre-filtering

Gram-Schmidt orthonormalization is a very fast and efficient method for the fringe pattern phase demodulation. It requires only two arbitrarily phase-shifted frames. Images are treated as vectors and upon orthogonal projection of one fringe vector onto another the quadrature fringe pattern pair is obtained. Orthonormalization process is very susceptible, however, to noise, uneven background and amplitude modulation fluctuations. The Hilbert-Huang transform based preprocessing is proposed to enhance fringe pattern phase demodulation by filtering out the spurious noise and background illumination and performing fringe normalization. The Gram-Schmidt orthonormalization process error analysis is provided and its filtering-expanded capabilities are corroborated analyzing DSPI fringes and performing amplitude demodulation of Bessel fringes. Synthetic and experimental fringe pattern analyses presented to validate the proposed technique show that it compares favorably with other pre-filtering schemes, i.e., Gaussian filtering and continuous wavelet transform.

Adaptive enhancement of optical fringe patterns by selective reconstruction using FABEMD algorithm and Hilbert spiral transform

Presented method for fringe pattern enhancement has been designed for processing and analyzing low quality fringe patterns. It uses a modified fast and adaptive bidimensional empirical mode decomposition (FABEMD) for the extraction of bidimensional intrinsic mode functions (BIMFs) from an interferogram. Fringe pattern is then selectively reconstructed (SR) taking the regions of selected BIMFs with high modulation values only. Amplitude demodulation and normalization of the reconstructed image is conducted using the spiral phase Hilbert transform (HS). It has been tested using computer generated interferograms and real data. The performance of the presented SR-FABEMD-HS method is compared with other normalization techniques. Its superiority, potential and robustness to high fringe density variations and the presence of noise, modulation and background illumination defects in analyzed fringe patterns has been corroborated.

A generalized regularized phase tracker for demodulation of a single fringe pattern

The regularized phase tracker (RPT) is one of the most powerful approaches for demodulation of a single fringe pattern. However, two disadvantages limit the applications of the RPT in practice. One is the necessity of a normalized fringe pattern as input and the other is the sensitivity to critical points. To overcome these two disadvantages, a generalized regularized phase tracker (GRPT) is presented. The GRPT is characterized by two novel improvements. First, a general local fringe model that includes a linear background, a linear modulation and a quadratic phase is adopted in the proposed enhanced cost function. Second, the number of iterations in the optimization process is proposed as a comprehensive measure of fringe quality and used to guide the demodulation path. With these two improvements, the GRPT can directly demodulate a single fringe pattern without any pre-processing and post-processing and successfully get rid of the problem of the sensitivity to critical points. Simulation and experimental results are presented to demonstrate the effectiveness and robustness of the GRPT.

Regularized frequency-stabilizing method for single closed-fringe interferogram demodulation

We present a simple and fast regularized frequency-stabilizing method for single open- or closed-fringe interferogram demodulation. The proposed method recovers the phase maps of interferograms by establishing a cost function, according to prior knowledge. Because only the phase field to be estimated is employed in the cost function, the optimization process is fast. Moreover, the recovered phase is continuous, and no further phase unwrapping is necessary. Computer simulation and experimental results have demonstrated both the rapidity and the efficiency of the method.

Normalization of fringe patterns using the bidimensional empirical mode decomposition and the Hilbert transform

We evaluate a data-driven technique to perform bias suppression and modulation normalization of fringe patterns. The proposed technique uses a bidimensional empirical mode decomposition method to decompose a fringe pattern in a set of intrinsic frequency modes and the partial Hilbert transform to characterize the local amplitude of the modes in order to perform the normalization. The performance of the technique is tested using computer simulated fringe patterns of different fringe densities and illumination defects with high local variations of the modulation, and its advantages and limitations are discussed. Finally, the performance of the normalization approach in processing real data is also illustrated.

Multiresolution monogenic signal analysis using the Riesz-Laplace wavelet transform

The monogenic signal is the natural 2-D counterpart of the 1-D analytic signal. We propose to transpose the concept to the wavelet domain by considering a complexified version of the Riesz transform which has the remarkable property of mapping a real-valued (primary) wavelet basis of L(2) (R(2)) into a complex one. The Riesz operator is also steerable in the sense that it give access to the Hilbert transform of the signal along any orientation. Having set those foundations, we specify a primary polyharmonic spline wavelet basis of L(2) (R(2)) that involves a single Mexican-hat-like mother wavelet (Laplacian of a B-spline). The important point is that our primary wavelets are quasi-isotropic: they behave like multiscale versions of the fractional Laplace operator from which they are derived, which ensures steerability. We propose to pair these real-valued basis functions with their complex Riesz counterparts to specify a multiresolution monogenic signal analysis. This yields a representation where each wavelet index is associated with a local orientation, an amplitude and a phase. We give a corresponding wavelet-domain method for estimating the underlying instantaneous frequency. We also provide a mechanism for improving the shift and rotation-invariance of the wavelet decomposition and show how to implement the transform efficiently using perfect-reconstruction filterbanks. We illustrate the specific feature-extraction capabilities of the representation and present novel examples of wavelet-domain processing; in particular, a robust, tensor-based analysis of directional image patterns, the demodulation of interferograms, and the reconstruction of digital holograms.

Estimation of phase derivatives using discrete chirp-Fourier-transform-based method

Estimation of phase derivatives is an important task in many interferometric measurements in optical metrology. This Letter introduces a method based on discrete chirp-Fourier transform for accurate and direct estimation of phase derivatives, even in the presence of noise. The method is introduced in the context of the analysis of reconstructed interference fields in digital holographic interferometry. We present simulation and experimental results demonstrating the utility of the proposed method.

Windowed high-order ambiguity function method for fringe analysis

This paper introduces a windowed high-order ambiguity function (WHAF) method for the demodulation of fringe patterns recorded in holographic interferometry. It first obtains the analytic signal of the fringe pattern and models it as a piecewise polynomial phase signal. A parametric estimation procedure based on HAF is then employed to calculate the polynomial coefficients of the phase over each window of the segmented analytic signal. A salient feature of the proposed method is that it provides an accurate and direct estimation of the unwrapped phase distribution from a single fringe pattern, even when the pattern's phase is rapidly varying. WHAF's application to both digital and classical holographic interferometry is demonstrated by simulation and experimental results.

Fast phase recovery from a single closed-fringe pattern

A new framework for phase recovery from a single fringe pattern with closed fringes is proposed. Our algorithm constructs an unwrapped phase from previously computed phases with a simple open-fringe-analysis algorithm, twice applied for analyzing horizontal and vertical oriented fringes, respectively. That reduces the closed-fringe-analysis problem to that of choosing the better phase between the two oriented computed phases and then of estimating the local sign. By propagating the phase sign [and a tilewise constant (DC) term] by regions [here named tiles] instead of a pixelwise phase propagation, our analysis of closed-fringe patterns becomes more robust and faster. Additionally, we propose a multigrid refinement for improving the final computed phase.

Exploration of event-induced EEG phase synchronization patterns in cognitive tasks using a time–frequency-topography visualization system

In this paper, we present a method for the study of synchronization patterns measured from EEG scalp potentials in psychophysiological experiments. This method is based on various techniques: a time-frequency decomposition using sinusoidal filters which improve phase accuracy for low frequencies, a Bayesian approach for the estimation of significant synchrony changes, and a time-frequency-topography visualization technique which allows for easy exploration and provides detailed insights of a particular experiment. Particularly, we focus on in-phase synchrony using an instantaneous phase-lock measure. We also discuss some of the most common methods in the literature, focusing on their relevance to long-range synchrony analysis; this discussion includes a comparison among various synchrony measures. Finally, we present the analysis of a figure categorization experiment to illustrate our method.