A fast-Fourier-transform method of topography and interferometry is proposed. By computer processing of a noncontour type of fringe pattern, automatic discrimination is achieved between elevation and depression of the object or wave-front form, which has not been possible by the fringe-contour-generation techniques. The method has advantages over moire topography and conventional fringe-contour interferometry in both accuracy and sensitivity. Unlike fringe-scanning techniques, the method is easy to apply because it uses no moving components.

Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry

Synthetic aperture radar (SAR) is a coherent active microwave imaging method. In remote sensing it is used for mapping the scattering properties of the Earth's surface in the respective wavelength domain. Many physical and geometric parameters of the imaged scene contribute to the grey value of a SAR image pixel. Scene inversion suffers from this high ambiguity and requires SAR data taken at different wavelength, polarization, time, incidence angle, etc. Interferometric SAR (InSAR) exploits the phase differences of at least two complex-valued SAR images acquired from different orbit positions and/or at different times. The information derived from these interferometric data sets can be used to measure several geophysical quantities, such as topography, deformations (volcanoes, earthquakes, ice fields), glacier flows, ocean currents, vegetation properties, etc. This paper reviews the technology and the signal theoretical aspects of InSAR. Emphasis is given to mathematical imaging models and the statistical properties of the involved quantities. Coherence is shown to be a useful concept for system description and for interferogram quality assessment. As a key step in InSAR signal processing two-dimensional phase unwrapping is discussed in detail. Several interferometric configurations are described and illustrated by real-world examples. A compilation of past, current and future InSAR systems concludes the paper.

https://iopscience.iop.org/article/10.1088/0266-5611/14/4/001/pdf

Synthetic aperture radar interferometry

A self-scanned 1024 element photodiode array and minicomputer are used to measure the phase (wavefront) in the interference pattern of an interferometer to lambda/100. The photodiode array samples intensities over a 32 x 32 matrix in the interference pattern as the length of the reference arm is varied piezoelectrically. Using these data the minicomputer synchronously detects the phase at each of the 1024 points by a Fourier series method and displays the wavefront in contour and perspective plot on a storage oscilloscope in less than 1 min (Bruning et al. Paper WE16, OSA Annual Meeting, Oct. 1972). The array of intensities is sampled and averaged many times in a random fashion so that the effects of air turbulence, vibrations, and thermal drifts are minimized. Very significant is the fact that wavefront errors in the interferometer are easily determined and may be automatically subtracted from current or subsequent wavefrots. Various programs supporting the measurement system include software for determining the aperture boundary, sum and difference of wavefronts, removal or insertion of tilt and focus errors, and routines for spatial manipulation of wavefronts. FFT programs transform wavefront data into point spread function and modulus and phase of the optical transfer function of lenses. Display programs plot these functions in contour and perspective. The system has been designed to optimize the collection of data to give higher than usual accuracy in measuring the individual elements and final performance of assembled diffraction limited optical systems, and furthermore, the short loop time of a few minutes makes the system an attractive alternative to constraints imposed by test glasses in the optical shop.

Digital wavefront measuring interferometer for testing optical surfaces and lenses

We describe what is to our knowledge a novel technique for phase unwrapping. Several algorithms based on unwrapping the most-reliable pixels first have been proposed. These were restricted to continuous paths and were subject to difficulties in defining a starting pixel. The technique described here uses a different type of reliability function and does not follow a continuous path to perform the unwrapping operation. The technique is explained in detail and illustrated with a number of examples.

Fast two-dimensional phase-unwrapping algorithm based on sorting by reliability following a noncontinuous path.

Phase shifting algorithms for fringe projection profilometry: A review

We have been reporting several new techniques of analysis and synthesis applied to Phase Shifting Interferometry (PSI). These works are based upon the Frequency Transfer Function (FTF) and how this new tool of analysis and synthesis in PSI may be applied to obtain very general results, among them; rotational invariant spectrum; complex PSI algorithms synthesis based on simpler first and second order quadrature filters; more accurate formulae for estimating the detuning error; output-power phase noise estimation. We have made our cases exposing these aspects of PSI separately. Now in the light of a better understanding provided by our past works we present and expand in a more coherent and holistic way the general theory of PSI algorithms. We are also providing herein new material not reported before. These new results are on; a well defined way to combine PSI algorithms and recursive linear PSI algorithms to obtain resonant quadrature filters.

The general theory of phase shifting algorithms

Isotropic n-dimensional fringe pattern normalization

Algorithm for fringe pattern normalization

Quadrature operators are useful for obtaining the modulating phase ϕ in interferometry and temporal signals in electrical communications. In carrier-frequency interferometry and electrical communications, one uses the Hilbert transform to obtain the quadrature of the signal. In these cases the Hilbert transform gives the desired quadrature because the modulating phase is monotonically increasing. We propose an n-dimensional quadrature operator that transforms cos(ϕ) into -sin(ϕ) regardless of the frequency spectrum of the signal. With the quadrature of the phase-modulated signal, one can easily calculate the value of ϕ over all the domain of interest. Our quadrature operator is composed of two n-dimensional vector fields: One is related to the gradient of the image normalized with respect to local frequency magnitude, and the other is related to the sign of the local frequency of the signal. The inner product of these two vector fields gives us the desired quadrature signal. This quadrature operator is derived in the image space by use of differential vector calculus and in the frequency domain by use of a n-dimensional generalization of the Hilbert transform. A robust numerical algorithm is given to find the modulating phase of two-dimensional single-image closed-fringe interferograms by use of the ideas put forward.

/pdf/general-n-dimensional-quadrature-transform-and-its-49key0gtq9.pdf

General n-dimensional quadrature transform and its application to interferogram demodulation.

In recent years phase-measuring techniques have been applied to the problem of extracting information of photoelastic data. We present a new phase-measuring algorithm for extraction of the isochromatics of photoelastic fringe patterns. The algorithm permits the extraction of the isochromatic phase with almost no influence from the isoclinics, thus avoiding the usual problems of low-modulation areas associated with isoclinics. The isochromatic phase map obtained with this algorithm is well suited for a full separation of the stress components in a sample. The algorithm can be used with any commercial diffuse-light circular polariscope.

/pdf/phase-measuring-algorithm-for-extraction-of-isochromatics-of-4c5ymq125i.pdf

Juan Antonio Quiroga

Papers

The general theory of phase shifting algorithms

Isotropic n-dimensional fringe pattern normalization

Algorithm for fringe pattern normalization

General n-dimensional quadrature transform and its application to interferogram demodulation.

Phase measuring algorithm for extraction of isochromatics of photoelastic fringe patterns