Open access peer-reviewed chapter

Fringe Pattern Analysis in Wavelet Domain

Written By

Yassine Tounsi, Abdulatef Ghlaifan, Manoj Kumar, Fernando Mendoza-Santoyo, Osamu Matoba and Abdelkrim Nassim

Submitted: 14 March 2019 Reviewed: 10 June 2019 Published: 01 July 2019

DOI: 10.5772/intechopen.87943

From the Edited Volume

Holographic Materials and Applications

Edited by Manoj Kumar

Chapter metrics overview

1,191 Chapter Downloads

View Full Metrics


We present a full-field technique for single fringe pattern analysis based on wavelet transform. Wavelets technique is a powerful method that quantifies at different scales how spatial energy is distributed. In the wavelets domain, fringe pattern analysis requires spatial modulation by a high-frequency carrier. We realize the modulation process numerically by combining the fringe pattern and its quadrature generated analytically by spiral phase transform. The first application concerns the speckle denoising by thresholding the two-dimensional stationary wavelet transform (2D-swt) coefficients of the detail sub-bands. In the second application, the phase derivatives are estimated from the 1D-continuous wavelet transform (1D-cwt) and 2D-cwt analysis of the modulated fringe pattern by extracting the extremum scales from the localized spatial frequencies. In the third application, the phase derivatives distribution is evaluated from the modulated fringe pattern by the maximum ridge of the 2D-cwt coefficients. The final application concerns the evaluation of the optical phase map using two-dimensional discrete wavelet transform (2D-dwt) decomposition of the modulated fringe pattern. The optical phase is computed as the arctangent function of the ratio between the detail components (high-frequency sub-bands) and the approximation components (low-frequency sub-bands). The performance of these methods is tested on numerical simulations and experimental fringes.


  • continuous wavelets transform (cwt)
  • stationary wavelets transform (swt)
  • discrete wavelets transform (dwt)
  • fringe pattern analysis
  • residual speckle denoising

1. Introduction

The single fringe pattern is an emerging technique for the analysis of full-field measurements in optical metrology. Fringe pattern analysis becomes a key technique in interferometric metrology which is more suitable for dynamic processes. The main purpose of this chapter is to exploit wavelet transform for single fringe pattern analysis to extract useful information. Wavelets are a powerful method allowing to know the spatial energy distribution at several scales.

Fringe pattern analysis using wavelets transform requires a spatial modulation by a high-frequency carrier as the case of Fourier domain [1]. We realize the modulation process digitally [2] by combining the fringe pattern and its quadrature generated analytically by spiral phase transform SPT [3, 4]. The advantage of this method lies in its implementation of numerical algorithms to reduce experimentation burden present in the phase shifting methods [5], namely, multiple fringe pattern generation and the experimental carrier introduction.

In the first application, we describe a residual speckle noise reduction technique by 2D-stationary wavelet decomposition of the speckle correlation fringe pattern [6]. The speckle noise is removed by thresholding the 2D-swt detail sub-bands coefficients [7].

The second application presents a method for phase derivative estimation based on 1D-continuous wavelet analysis of the modulated fringe pattern using Paul’s function as the mother wavelet [8]. The phase derivative is calculated by extracting the extremum scales from the localized spatial frequencies. The third application is phase derivatives evaluation by 2D-continuous wavelet analysis of the modulated fringe pattern using complex Morlet’s function as the mother wavelet [9]. The phase derivative is estimated by a maximum ridge of wavelets coefficients. Finally, we introduce an algorithm for optical phase extraction by 2D-Discrete wavelet decomposition of the modulated fringe pattern using Gabor’s function as the mother wavelet [10]. The optical phase is extracted as the ratio between the detail components (high-frequency sub-bands) and the approximation components (low-frequency sub-bands).

The remaining part of the chapter proceeds as follow: divided into three sections: the second section will examine a brief description of continuous, discrete and stationary wavelet transforms. In the third section, we briefly introduce the speckle correlation fringes obtained by digital speckle pattern interferometry (DSPI). Finally, the fourth section is devoted to applications of wavelets to fringe pattern analysis.


2. Wavelet transform analysis

The concept of the wavelets and its numerous fields of the application make them a useful tool in several studies concerning localized variations analysis for non-stationary or transient signal analysis. This concept is based on the multiresolution analysis that represents signal variations at different scales. A detailed review of wavelets theory has been published by Ingrid Daubechies in [11]. However, we present here a brief description of three wavelets families for completeness.

2.1. Continuous wavelet transform (cwt)

The analysis of a given signal by continuous wavelets transform concerns to decompose it into several basic functions named wavelets. They are oscillatory functions with a finite duration and having zero average value, also, they are characterized by irregularity and the good localization These properties of the continuous wavelets make them a superior basis for signals analysis with discontinuities. The wavelets are constructed by translating and dilating them other wavelet functions resulting in a self-similar wavelet’s families as follows.

ψ s , ξ x = ψ s 1 x ξ E1

where s and ξ are respectively the scale and translation parameters. The wavelet coefficients obtained by 1D-cwt decomposition of the signal f(x) are given by:

w s ξ = + f x . ψ s , ξ x . dx = s 1 / 2 + f x . ψ s 1 x ξ . dx E2

where ψ*s,ξ(x) represents the conjugate of the wavelet ψs,ξ(x) defined for each shift ξ and each scale s. Wavelet coefficients are the output of the correlation product between a signal and the mother wavelet for different values of dilatation, and this is interpreted as a measure of the local similarity between them.

The 2D-continuous wavelet coefficients are obtained by the computation of the correlation product between the input image (signal 2D) and the mother wavelet. In the 2D case, a parameter of orientation angle is added, and the 2D-cwt wavelet coefficients of an input image f(x,y) are defined as:

w ξ η s θ = s 2 f x y ψ s 1 R θ x ξ y η dxdy E3

where ξ and η are respectively the translation parameters along x- and y- directions, s and θ are respectively the scale vector and rotation angle; ψ denotes the 2D mother wavelet and Rθ represents the conventional 2×2 rotation matrix corresponding to θ given by

R θ = x . cos θ + y . sin θ y . cos θ x . sin θ E4

2.2. Discrete wavelets transform (dwt)

The discrete wavelet transform denoted by dwt is a fast computing algorithm; it was introduced by Mallat [12] in 1989 for signal or image decomposition into two important components called approximation and details.

The 1D-dwt decomposition is based on two functions called scaling and wavelets. These functions, i.e., the scaling and the mother wavelet functions are related to the impulse response of h[n] and g[n], respectively, as illustrated in Figure 1. The approximation (low-frequency) and details (high-frequency) obtained respectively by scaling and wavelets functions are the down-sampled outputs of the first filters h[n] and g[n].

Figure 1.

Block diagram of filter analysis.

The relation between the two functions and the two filters is expressed as:

ϕ x = 2 n h n φ 2 x n E5
ψ x = 2 n g n ψ 2 x n E6

g is the conjugate mirror of h:

g n = 1 n . h 1 n E7

In Eq. (7), conjugation and mirror effects are represented respectively by 1 n and (−n). The internal orthogonality relation is satisfied by the low-pass filter h as follows:

n h n h n + 2 j = 0 E8

and to have

n h n 2 = 1 E9

The filter g satisfies h, the same internal orthogonality and both obey the mutual orthogonality relation expressed as following:

n h n g n + 2 j = 0 E10

In 2D-dwt, the wavelet transform of an image involves recursive filtering and sub-sampling process. Figure 2 describes the procedure of dwt analysis, three detail images, and one approximation image are obtained at each level. Concerning detail images, they contain the high-frequency information we denote horizontal image sub-band by HD, vertical image sub-band by VD and the diagonal image sub-band by DD and the approximation image sub-band is denoted by AI which accommodates the low-frequency information.

Figure 2.

Diagrams of dwt image decomposition.

Figure 3.

Block diagram of filter analysis.

The 2D scaling function denoted by ϕ(x,y), and the 3D wavelet ψH(x,y), ψV(x,y), and ψD(x,y) are computed by the algebraic product between the one dimensional scaling and wavelet function as expressed in the following Eq. [13],

ϕ x y = ϕ x . ϕ y ψ H x y = ψ x . ϕ y ψ V x y = ϕ y . ψ x ψ D x y = ψ x . ψ y E11

The three functions ψV, ψH, and ψD provide respectively the vertical, horizontal and diagonal variations, from there, the three details images are obtained.

2.3. Stationary wavelet transform (swt)

Stationary wavelets transform (swt) was presented in [14]. Swt has the same principle of decomposition as dwt transform, but the process of down-sampling is eliminated which means the swt is translation-invariant. The 2D-swt is founded on the idea of no down-sampling. Specifically, this transform is applied at each pixel of the image and saves the detail coefficients and exploits the low-frequency data at each level. A clear description of this technique is detailed in [15].

The principle of swt analysis is schematized in Figure 3. In brief, the swt technique consists to decompose an input sequence at each level and the given output is a new sequence that has the same length as the original sequences. In order to implement this transform, the process of original image decimation is removed, nonetheless, the size of the filter is modified at each level by zero paddings.

We introduce an operator denoted Z ^ that alternates an input sequence with zeros so that for all the integers j, we have ( Z ^ x )2j = xj and ( Z ^ x )2j + 1 = 0. It is also assumed the filters Hr and Gr to have weights Z ^ r h and Z ^ r g , respectively. Thus, the filter Hr is characterized by a weight hr2rj = hj and hrkj = 0 in the case where k is not a multiple of 2r. The filter Hr is attained by introducing a zero between each adjacent pair of the filter H(r−1) elements, and an identical trend should be followed for Gr. The following relations satisfy the above statements:

D 0 r H r = HD 0 r and D 0 r G r = GD 0 r E12

We start the swt definition by setting uJ to be the original input sequence, the stationary wavelet transform can be described for j = J, J-1, …,1,

u j 1 = H J j u j and v j 1 = G J j u j E13

The vectors uj and vj acquire the same length for the vector uJ with a length of 2J.


3. Speckle correlation fringe analysis

3.1. Speckle effect

The speckle appears as a granular structure and it results from the self-interference of a large number of coherent waves randomly scattered from and/or transmitted through a rough object surface [16]. When we illuminate a porous surface with coherent light, the scattered light intensity has a random spatial variation called “speckle effect.” The “speckle pattern” appears chaotic and disordered; it is described by using statistics and probability. The speckle pattern structure is dependent on the coherence properties of the used illumination light and also on the characteristics of the object’s surface diffusion. Locally in space, the electric field of the speckle pattern is given by computing the contribution sum of all illuminated scattering elements of the rough object’s surface and is given by [16]

E x y z = k = 1 n a k exp j ϕ k E14

In Eq. (14), n represents the scattering number of elements; ak and ϕk are respectively the amplitude and phase of the kth scatterer contribution. The theorem of the central limit state that the random variable resulting from the sum of several independent random variables results in a Gaussian distribution when their number (number of an independent random variable) approaches infinity.

In the case where the number of scatterers is large, by applying the central Limit we deduce that:

  1. The two components (real and imaginary) of the field are identically distributed Gaussian variables, independent and with zero means.

  2. The intensity I(x,y,z) is characterized by a negative exponential probability distribution expressed as [16].

p I = exp I / < I > / < I > E15

The associated phase ϕ is uniformly distributed between −π and π [16]

p ϕ = 1 / 2 π E16

3.2. Digital speckle pattern interferometry (DSPI)

Digital/electronic speckle pattern interferometry (DSPI/ESPI) is an optical interferometric technique demonstrated simultaneously by Macovski et al. [17], and Butters and Leendertz [18, 19] at the beginning of the 1970s. Later on, the technique was enriched by Beidermann and Ek [20] and Lokberg and Hogmoen [21] with several new investigations. Figure 4 shows the typical schematic of the DSPI system with a sensitivity vector responding to out-of-plane deformation/displacement measurements (Figure 4a), and the digital electronics and image processing unit (Figure 4b). The DSPI technique measures the phase changes introduced by the speckle intensity changes. In this technique, a spatially filtered reference beam is added to the speckle pattern which is scattered from or transmitted through the test object, to code its phase. Therefore, a speckled interferogram, resulting from the superposition of the speckle pattern and the reference beams, contains essential information about the random phase. In DSPI, as shown in Figure 4, two speckle interferograms, corresponding to two different states of the object are recorded and stored in the frame grabber card. For example, one speckle interferogram, A, is recorded when the object is in its initial state, say the reference state, while another speckle interferogram, B, is recorded in the deformed state of the object, say by a small distance “d,” as shown in Figure 4a. The DSPI fringe pattern is observed after subtraction of these two speckle interferograms A and B by using the appropriate software. The useful information about the object is encoded in the DSPI fringe pattern like the one shown in Figure 4b.

Figure 4.

Schematic of the DSPI system: (a) Speckle interferometer and (b) Digital electronics and image processing unit.

The DSPI system can be made sensitive to out-of-plane or in-plane deformations, or both, depending on the design of the optical setup. Numerous measurement applications of DSPI have been investigated with different proposed experimental configurations [22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48] which have immense importance in the scientific, engineering, and industrial fields. The popularity of the technique, in optical metrology, is by virtue of the several advantages it offers such as: it is a non-contact/invasive type technique; it provides full-field of view information; it is faster in operation as almost real-time observations can be obtained; investigations on large-size objects can be observed with the use of proper optics, and the technique is less sensitive to environmental perturbations in comparison with its counterparts. The DSPI technique has proven its capability in many investigations, e.g., for the measurement of displacement/deformation with variable sensitivity to in-plane and out-of-plane direction [22, 23, 24, 25, 26, 27], measurement of three-dimensional shape [28, 29, 30], surface roughness measurement [31], vibration measurement/monitoring [32, 33, 34], measurement of material properties [35, 36, 37], flow visualization [38, 39], measurement of refractive index and temperature distribution [40, 41, 42], and the investigation of the magnetic fields on the temperature profile of gaseous flames [43, 44]. A Doctoral Thesis describing the comprehensive study of DSPI and some of its novel investigations is recommended for interested readers working in this field [48].

Mathematically, the fringe formation in DSPI is demonstrated here. The intensity distribution recorded before the object deformation is expressed as follows:

f B x y = b x y + a x y . cos ϕ s x y E17

In this equation, the term b(x,y) represents the background or bias, a(x,y) the visibility and term ϕs(x,y) is the original speckle phase appearing as a high frequency, and the intensity of pixel by pixel is randomly distributed. When the object is deformed, the intensity distribution of speckle-pattern becomes:

f A x y = b x y + a x y . cos ϕ s + φ . E18

where φ(x,y) represents the phase variation stemming from the object deformation. Assuming that the introduced deformations are sufficiently small, we can ignore the speckle decorrelation effects.

The speckle fringe correlation intensity distribution in the subtraction mode is expressed as:

f = f A f B = 2 . a . sin φ / 2 . sin ϕ s + φ / 2 E19

The above equation can be further expressed as:

f 2 = 4 a 2 . sin 2 φ / 2 . sin 2 ϕ s + φ / 2 E20

The speckle phase ϕs(x,y) changes faster across the speckle pattern, the second sine squared term has an ensemble average expressed as:

sin 2 ϕ s + φ / 2 1 / 2 πn 0 2 π sin 2 ϕ s + φ / 2 d ϕ s + φ / 2 = 1 / 2 E21

Hence, the speckle fringe correlation intensity distribution becomes:

f C = f 2 = 2 a 2 . sin 2 φ / 2 = a 2 . 1 cos φ E22

The obtained intensity distribution fc is rendered as the classical form of a cosine fringe pattern.

3.3. Fringe pattern analysis

The analysis of the fringe patterns concerns to evaluate the coded phase distribution related to the physical magnitude. We classified the phase extraction techniques to methods that explore a shifted fringe pattern as the phase shifting techniques and methods that require the introduction of a spatial carrier such as the Fourier transform method and the wavelet method.

The phase shifting (see Ref. [5]) and the Fourier transform methods (see Ref. [1]) are the most common methods exploited for fringe pattern analysis. The goal of the next section is to present algorithms based on wavelets transform to analyze a single frame speckle fringe pattern.


4. Applications to speckle correlation fringe analysis

As we stated previously, the fringe pattern analysis in wavelet domain requires a modulation process by introducing a spatial carrier. In our case, the spatial carrier in generated numerically and the entire process is carried out directly using a computer. Before presenting the four applications, we present in the next subsection the modulation process required by wavelets transform.

4.1. Fringe pattern modulation

In order to modulate a given fringes pattern in a chosen direction, we combine the fringe pattern and its quadrature with a cos(m.x) and sin(m.x) matrix respectively where m means the modulation rate. We consider the fringe pattern intensity distribution expressed as follow:

f x y = b x y + a x y cos φ x y E23

A low-pass filter applied to the intensity distribution removes the background illumination and Eq. (23) becomes:

f ˜ x y = a x y cos φ x y E24

Larkin proposed recently a transform called spiral phase quadrature denoted -SPT- for the 2D fringe pattern (see Refs. [3, 4]). The SPT of f is computed as:

SPT f ˜ = j exp j D x y a x y sin φ E25

where, a(x,y).sinφ is the quadrature term, j2 = −1, and D(x,y) represents the direction map. Then, we obtain the sine fringe pattern as:

a x y sin φ = j exp j D x y SPT f ˜ E26

The direction map is computed by the ratio between the horizontal and vertical gradient of the f as and analytically expressed as:

tan D = y f ˜ / x f ˜ E27

So, we define the quadrature map as:

q x y = b x y sin φ = j exp j D SPT f ˜ E28

Then, the intensity distribution of the modulated fringe pattern is defined as:

f m = f ˜ cos m x q sin m x = b x y cos φ + m x E29

4.2. Speckle noise reduction using various wavelet-based techniques: a simulation study

4.2.1. Speckle noise reduction using 2D-swt

As a reminder, the speckle fringe pattern intensity distribution in the subtraction mode is expressed in (Eq. (19)) as f = 2 . a . sin φ / 2 . sin ϕ s + φ / 2 . The speckle fringes, presented in Eq. (19), are characterized by multiplicative residual speckle noise, therefore, a denoising step is necessary before the evaluation of the phase distribution and their derivatives. The high-frequency sin(ϕs + φ/2) noise should be removed with an appropriate filtering technique. After decomposition by swt of the speckle fringe pattern, the noise is reduced by thresholding the detail sub-bands coefficients using a soft threshold function. The choice of this type of thresholding stems from its characteristic to obtain near-optimal minimax rate. Generally, the threshold estimation requires the knowledge of the noise variance and the optimal threshold is estimated by taking into consideration the variance (σ2) of the coefficients in the highest wavelets decomposition. Figure 5 illustrates the capability of stationary wavelet transform thresholding technique to reduce the speckle noise. The left and the right halves present the intensity distribution of speckle fringe correlation before and after denoising. It is clearly shown in the line profile plot in Figure 5b along the line AB, the influence of swt to remove the frequency from the noised image.

Figure 5.

(a) Speckle fringe correlation before and after denoising; and (b) the 1D intensity profile along line AB.

4.2.2. Phase derivatives extraction using 1D-cwt

The idea is the phase extraction by spatial frequencies using the1D-continuous wavelet transform and Paul’s function as the mother wavelet. The phase derivatives are evaluated from the modulus of the cwt coefficients by extracting the extremum scales from the localized spatial frequencies and the phase is obtained by numerical integration.

Generally, the modulated fringe pattern intensity distribution, as represented by Eq. (18), can also be expressed as:

f m x y = b x y + a x y cos m . x + φ x y E30

Using Eq. (2), the wavelet transform coefficients of the modulated fringe pattern, given in Eq. (30), can be represented as

w x s ξ = s 1 / 2 + b x y + a ( x y ) cos m . x + φ x y ψ y ξ / s dy E31

The localization property of the wavelet is taken into account in order to write the Taylor series of the phase φ(x,y) near the central value ξ as:

φ x y = φ x ξ + y ξ φ x y / y + E32

Assuming that the bias a and the visibility b have a slow variation, so the higher order of (y-ξ) is neglected with respect to the phase-modulated carrier because of the localization of the wavelet. With these considerations, the 1D-cwt coefficients are now computed as:

w x s ξ = s 1 / 2 b x ξ a x ξ + cos mx + φ x ξ + y ξ φ ( x ξ ) / y ψ y ξ / s dy E33

and the Parseval’s identity gives

w x s ξ = s 2 π + f a ̂ x k ψ ̂ sk e iξk dk E34


f a ̂ x k = b x ξ a x ξ πh x k E35


h x k = δ k m 1 exp i φ x ξ ξ φ ( x ξ ) / y + δ k + m 1 exp i φ x ξ ξ φ ( x ξ ) / y E36
m 1 = m + φ x ξ / x E37

The wavelet transform reduces to.

w x s ξ = 0.5 . a x ξ b x ξ s ψ ̂ sm 1 exp i ξm + φ x ξ + ψ ̂ sm 1 exp ( i ( ξm + φ ( x ξ ) ) ) E38

Introducing Paul’s mother wavelet defined as:

ψ x = 2 n n ! 1 ix n + 1 / 2 π 2 n ! / 2 E39

we get

w x s ξ = 2 n ! 1 a x ξ b x ξ s n + 1 / 2 m 1 n exp sm 1 × exp i ξm + φ x ξ E40

whose modulus is computed as

w x s ξ = 2 n ! 1 a x ξ b x ξ m 1 n s n + 1 / 2 exp sm 1 E41

The extremum scale denoted by S is represented by

S x ξ = 2 n + 1 / 2 m 1 E42

Equations (37) and (42) provide the phase derivative as

φ x ξ / y = 2 n + 1 / 2 S x ξ m E43

where m is the modulation ratio.

The performance of this algorithm is illustrated in Figure 6, where the phase derivatives of the speckle fringe correlation along x- and y-directions are presented.

Figure 6.

The evaluated phase derivatives along (a) x-direction and (b) y-direction using 1D-cwt technique.

4.2.3. Phase derivatives extraction using 2D-cwt

The idea is to evaluate the phase derivative by ridge point using2D-continuous wavelet transform and Morlet’s function as the mother wavelet. The phase derivative is computed from maximums scales relating to the ridge point of the wavelet coefficient modules. From the modulated fringe pattern fm(x, y), we compute the 2D-cwt wavelet coefficients as defined in Eq. (3) and represented again as

w t d s θ = s 2 f m x y ψ M s 1 R θ x t y d dxdy E44

The 2D complex Morlet is used as the mother wavelet: it is essentially a plane wave modulated by a Gaussian window, and is expressed as:

ψ M = exp x 2 + y 2 / 2 exp ( jk 0 x cos θ + y sin θ E45

where k0 is the specific spatial frequency that should be in the range 5–6 in order to satisfy the admissibility condition [49], and j2 = −1.

A new matrix containing the maximum value of each column of the wavelet coefficient modulus array defines the wavelet ridge, and then, the corresponding scale value is determined from the ridge wavelet [50, 51]. The maximum scales smax correspond to the maximum ridge of the obtained wavelet coefficients modulus:

s max θ = arg max s R + , θ 0 2 π w t d s θ E46

where smax represents the scale value for maxima. The phase derivative is obtained by

φ = k 0 + k 0 2 + 2 1 / 2 / 2 s max m E47

The horizontal and vertical phase derivatives evaluated from the speckle fringe correlation are presented in Figure 7.

Figure 7.

The evaluated phase derivatives along (a) x-direction and (b) y-direction using 2D-cwt technique.

4.2.4. Optical phase extraction using 2D-dwt

This application concerns to extract the optical phase distribution by dwt components using the 2D discrete wavelets transform and Gabor’s function as the mother wavelet. The 2D-dwt method decomposes the fringe pattern into two components: approximation (low-frequency) and details (high-frequency) using two principal quadrature mirror decomposition filters: low-pass and high-pass filters.

We consider the modulated fringes expressed in Eq. (30) as f m x y = b x y + a x y . cos m . x + φ x y . It depends on the desired phase φ(x,y). Figure 8 presents the four sub-bands of the modulated speckle fringe correlation after 2D-dwt decomposition.

Figure 8.

The four sub-bands images after modulated fringe 2D-dwt decomposition.

The optical phase distribution coded in the modulated fringe pattern is evaluated by computing the arctangent function of the ratio between the detail and approximation images as following:

φ x y = atan 2 f m h / f m g E48

Inside the atan2 function, the numerator and denominator present the convolution product between fm and the two filters h and g, and give the detail and approximation images respectively. An illustration of the two images is shown in Figure 9.

Figure 9.

Approximation and detail image of speckle fringe correlation.

The detailed image in the numerator is chosen according to the direction of modulation. The low-pass filter h and high-pass filters g are respectively the real and imaginary part of Gabor function expressed as:

h = exp x 2 / σ x 2 + y 2 / σ y 2 / 2 . cos 2 π f 0 x cos θ + y sin θ E49
g = exp x 2 / σ x 2 + y 2 / σ y 2 / 2 . sin 2 π f 0 x cos θ + y sin θ E50

The 2D Gabor’s function [52] is the extension of the 1D Gabor function introduced by Gabor [53]. It is defined by a Gaussian window that modulates a complex exponential centered at a fixed frequency It is formulated by:

G = exp x 2 / σ x 2 + y 2 / σ y 2 / 2 . exp j .2 π f 0 x cos θ + y sin θ E51

Central frequency, orientation, spatial extent, and aspect are four essential parameters characterize Gabor function, and they can be adjusted by varying respectively f0, θ, σx and σy. The use of the atan2 function in Eq. (48) provides the phase distribution modulo 2π as shown in Figure 10a. To remove this discontinuity, a phase unwrapping step is obligatory [54]; for this reason, we have implemented the PUMA algorithm [55]. The unwrapped phase illustration in three-dimensional representations is presented in Figure 10b.

Figure 10.

(a) Wrapped phase distribution retrieved by 2-dwt; and (b) continuous phase distribution volume visualization.

4.3. Discussion of wavelet techniques on the numerically simulated fringe correlations

In order to study the performance and validity of the algorithms, an image quality assessment (Q) is used [56]. This metric model any image distortion by combining three factors. The first factor is the loss of correlation between the ideal and the obtained images denoted O and E respectively, which measures the linear correlation degree between the two input images. This factor is defined as:

Q 1 = σ OE / σ O σ E E52

The second factor is luminance distortion, it measures the mean luminance between O and E, which its value is computed as:

Q 2 = 2 O ¯ . E ¯ / O ¯ 2 + E ¯ 2 E53

The third factor is contrast distortion that measures the contrast similarity between the two images and it is defined as:

Q 3 = 2 σ O σ E / σ O 2 + σ E 2 E54

And then, the Q coefficient is computed by the product of three factors:

Q = Q 1 × Q 2 × Q 3 E55

where O ¯ and E ¯ present average of O and E, respectively, σ O and σ E are respectively the standard deviations of the two images. The Q index takes values between −1 and 1 where 1 means that the retrieval characteristic is exact. The Table 1 below summarizes the computed Q index values.

Retrieved characteristic map Q index value Time in seconds
Horizontal phase gradient 0.95 (1-cwt) & 0.96 (2-cwt) 50.23 (1-cwt) & 20 (2-cwt)
Vertical phase gradient 0.95 (1-cwt) & 0.96 (2-cwt) 50.26 (1-cwt) & 20.11 (2-cwt)
Recovered phase distribution by 2-dwt 0.97 10.55

Table 1.

The measured Q index values.

As a note, the 1D-cwt analyses images line by line, whereas the 2D-cwt scans images with the help of the parameter of orientation. On the analysis time side, we used Matlab on a 2.93 GHz Intel Pentium processor machine with 4 GB RAM. According to the Q index values, the presented algorithms give the desired information with good accuracy (Q between 0.95 and 0.97).

4.4. Application on an experimental speckle fringe correlation

After the presentation of the four applications of the stationary, continuous and discrete wavelet transform for fringe pattern analysis, and the favorable effectiveness validated by computer simulation, we exploit this algorithmic arsenal to analyze an experimental speckle fringe pattern recorded by using the DSPI setup showed in Figure 11a.

Figure 11.

(a) Schematic of DSPI setup [57, 58], (b) Recorded DSPI fringe pattern, and (c) Denoised DSPI fringe pattern by using swt thresholding technique.

In DSPI, two specklegrams corresponding to the two different states (before and after deformation) of the object are recorded with an image sensor and then subtracted in order to get the speckle fringe correlation as shown in Figure 11b. Figure 11c shows the denoised fringe resulting from the swt thresholding technique. The wrapped phase distribution and corresponding three-dimensional continuous phase distribution are shown in Figures 12a and b, respectively. The horizontal and vertical phase derivatives are shown in Figure 13.

Figure 12.

(a) Wrapped phase distribution extracted using dwt algorithm. (b) Three-dimensional representation of continuous phase distribution.

Figure 13.

The estimated phase gradient distribution using 2-cwt technique along: (a) x-direction and (b) y-direction.


5. Conclusion

This chapter has introduced an important thematic in interferometric metrology namely, fringe pattern analysis. We have discussed the capability of wavelets transform families (stationary, continuous and discrete wavelet) to analyze the fringe patterns.

This analysis consists in extracting the optical phase distribution coded in the fringes and its horizontal and vertical derivatives. We have presented with the help of computer simulations four applications of the wavelet transform, in order to analyze the speckle fringes correlation obtained in digital speckle pattern interferometry (DSPI). We summarize these applications as:

  • Speckle noise reduction using stationary wavelet transform thresholding technique.

  • Phase derivative extraction using one and two-dimensional continuous wavelet transform.

  • Optical phase distribution extraction using the two-dimensional discrete wavelet transform



Manoj Kumar and Osamu Matoba acknowledge the support by Japan Society for the Promotion of Science KAKENHI (17F17370, 18H03888) and JST CREST Grant Number JPMJCR1755. All the authors wish to thank the anonymous reviewers for their great insight, comments, and suggestions.


  1. 1. Takeda M, Ina H, Kobayashi S. Fourier-transform method of fringe-pattern analysis for computer-based topography and interferometry. Journal of the Optical Society of America. 1982;72(1, 1):156-160
  2. 2. Barj EM, Afifi M, Idrissi AA, Rachafi S, Nassim K. A digital spatial carrier for wavelet phase extraction. Optik. 2005;116(11):507-510
  3. 3. Larkin KG, Bone DJ, Oldfield MA. Natural demodulation of two-dimensional fringe patterns. I. General background of the spiral phase quadrature transform. Journal of the Optical Society of America A. 2001;18(8):1862-1870
  4. 4. Larkin KG. Natural demodulation of two-dimensional fringe patterns. II. Stationary phase analysis of the spiral phase quadrature transform. Journal of the Optical Society of America A. 2001;18(8):1871-1881
  5. 5. Creath K. Phase-shifting speckle interferometry. Applied Optics. 1985;24(18):3053-3058
  6. 6. Barj EM, Afifi M, Idrissi AA, Nassim K, Rachafi S. Speckle correlation fringes denoising using stationary wavelet transform. Application in the wavelet phase evaluation technique. Optics and Laser Technology. 2006;38(7):506-511
  7. 7. Donoho DL, Johnstone JM. Ideal spatial adaptation by wavelet shrinkage. Biometrika. 1994;81(3):425-455
  8. 8. Afifi M, Fassi-Fihri A, Marjane M, Nassim K, Sidki M, Rachafi S. Paul wavelet-based algorithm for optical phase distribution evaluation. Optics Communication. 2002;211(1-6):47-51
  9. 9. Ghlaifan A, Tounsi Y, Muhire D, Nassim A. Phase gradient retrieval from fringes pattern by using of two-dimensional continuous wavelet transforms. International Journal of Optics and Applications. 2017;7(4):69-74
  10. 10. Ghlaifan A, Tounsi Y, Zada S, Muhire D, Nassim A. Two-dimensional discrete wavelets transform for optical phase extraction: Application on speckle correlation fringes. Optical Engineering. 2016;55(12):121708
  11. 11. Daubechies I. Ten Lectures on Wavelets. Philadelphia, Pennsylvania: Society for Industrial and Applied Mathematics (SIAM); 1992. ISBN: 0-89871-274-2.
  12. 12. Mallat SG. A theory for multiresolution signal decomposition: The wavelet representation. IEEE Transactions on Pattern Analysis and Machine Intelligence. 1989;(7):674-693
  13. 13. Gonzalez RC, Woods RE. Digital Image Processing. 3rd ed. Upper Saddle River, N.J: Pearson; 2007
  14. 14. Pesquet JC, Krim H, Carfantan H. Time-invariant orthonormal wavelet representations. IEEE Transactions on Signal Processing. 1996;44(8):1964-1970
  15. 15. Nason GP, Silverman BW. Lecture Notes in Statistics. Wavelets and Statistics. Vol. 103. New York: Springer; 1995. p. 281
  16. 16. Goodman JW. Laser speckle and related phenomenon. In: Topics in Applied Physics. Vol. 9. Berlin: Springer; 1975
  17. 17. Macovski A, Ramsey SD, Schaefer LF. Time-lapse interferometry and contouring using television systems. Applied Optics. 1971;10(12):2722-2727
  18. 18. Butters JN, Leendertz JA. Speckle pattern and holographic techniques in engineering metrology. Optics and Laser Technology. 1971;3:26-30
  19. 19. Butters JN, Leendertz JA. Holographic and video techniques applied to engineering measurement. Measurement and Control. 1971;4(12):349-354
  20. 20. Biedermann K, Ek L. A recording and display system for hologram interferometry with low resolution imaging devices. Journal of Physics E: Scientific Instruments. 1975;8(7):571
  21. 21. Løkberg OJ, Høgmoen K. Vibration phase mapping using electronic speckle pattern interferometry. Applied Optics. 1976;15(11):2701-2704
  22. 22. Wykes C. Use of electronic speckle pattern interferometry (ESPI) in the measurement of static and dynamic surface displacements. Optical Engineering. 1982;21(3):213400
  23. 23. Joenathan C, Franze B, Haible P, Tiziani HJ. Large in-plane displacement measurement in dual-beam speckle interferometry using temporal phase measurement. Journal of Modern Optics. 1998;45(9):1975-1984
  24. 24. Fricke-Begemann T, Burke J. Speckle interferometry: Three-dimensional deformation field measurement with a single interferogram. Applied Optics. 2001;40(28):5011-5022
  25. 25. Kumar M, Agarwal R, Bhutani R, Shakher C. Measurement of strain distribution in cortical bone around miniscrew implants used for orthodontic anchorage using digital speckle pattern interferometry. Optical Engineering. 2016;55(5):054101. DOI: 10.1117/1.OE.55.5.054101
  26. 26. Kumar M, Agarwal R, Bhutani R, Shakher C. Deformation measurements in cortical bone-miniscrew interface in human maxilla by using digital speckle pattern interferometry. In: Speckle 2018, VII International Conference on Speckle Metrology; 7 September 2018; Vol. 10,834. International Society for Optics and Photonics; 2018. p. 1083412. DOI: 10.1117/12.2317714
  27. 27. Pedrini G, Zou YL, Tiziani HJ. Simultaneous quantitative evaluation of in-plane and out-of-plane deformations by use of a multidirectional spatial carrier. Applied Optics. 1997;36(4):786-792
  28. 28. Creath K, Cheng YY, Wyant JC. Contouring aspheric surfaces using two-wavelength phase-shifting interferometry. Optica Acta: International Journal of Optics. 1985;32(12):1455-1464
  29. 29. Moore AJ, Tyrer JR. Two-dimensional strain measurement with ESPI. Optics and Lasers in Engineering. 1996;24(5-6):381-402
  30. 30. Yang L, Zhang P, Liu S, Samala PR, Su M, Yokota H. Measurement of strain distributions in mouse femora with 3D-digital speckle pattern interferometry. Optics and Lasers in Engineering. 2007;45(8):843-851
  31. 31. Lehman MM, Pomarico JA, Torroba RD. Digital speckle pattern interferometry applied to a surface roughness study. Optical Engineering. 1995;34(4):1148-1153
  32. 32. Creath K, Slettemoen GÅ. Vibration-observation techniques for digital speckle-pattern interferometry. Journal of the Optical Society of America A. 1985;2(10):1629-1636
  33. 33. Wang WC, Hwang CH, Lin SY. Vibration measurement by the time-averaged electronic speckle pattern interferometry methods. Applied Optics. 1996;35(22):4502-4509
  34. 34. Chen F, Griffen CT, Allen TE. Digital speckle interferometry: Some developments and applications for vibration measurement in the automotive industry. Optical Engineering. 1998;37(5):1390-1398
  35. 35. Ganesan AR. Measurement of poisson’s ratio using real-time digital speckle pattern interferometry. Optics and Lasers in Engineering. 1989;11(4):265-269
  36. 36. Kumar M, Gaur KK, Shakher C. Measurement of material constants (Young’s modulus and Poisson’s ratio) of polypropylene using digital speckle pattern interferometry (DSPI). Journal of the Japanese Society for Experimental Mechanics. 2015;15(Special_Issue):s87-s91. DOI: 10.11395/jjsem.15.s87
  37. 37. Kumar M, Khan GS, Shakher C. Measurement of elastic and thermal properties of composite materials using digital speckle pattern interferometry. In: SPECKLE 2015, VI International Conference on Speckle Metrology; 24 August 2015; Vol. 9660. International Society for Optics and Photonics; 2015. p. 966011. DOI: 10.1117/12.2196390
  38. 38. Song Y, Wu Y, Kulenovic R, Guo Z, Groll M. Digital speckle pattern interferometry (DSPI) using optical polarization phase shift for measurement of temperature field in thermal flow. In: Optical Technology in Fluid, Thermal, and Combustion Flow III; 21 November 1997; Vol. 3172. International Society for Optics and Photonics; 1997, pp. 442-453
  39. 39. Song Y, Zhang W, Wu Y, Yao X. Digital speckle technique applied to flow visualization. Tsinghua Science and Technology. 2000;5(1):89-95
  40. 40. Barbosa EA, dos Santos SC. Refractive and geometric lens characterization through multi-wavelength digital speckle pattern interferometry. Optics Communication. 2008;281(5):1022-1029
  41. 41. Kumar M, Kumar V, Shakher C. Measurement of temperature and temperature distribution in diffusion flames using digital speckle pattern interferometry. In: Eleventh International Conference on Correlation Optics; 17 December 2013; Vol. 9066. International Society for Optics and Photonics; 2013. p. 90660Y. DOI: 10.1117/12.2047634
  42. 42. Kumar M, Shakher C. Measurement of temperature and temperature distribution in gaseous flames by digital speckle pattern shearing interferometry using holographic optical element. Optics and Lasers in Engineering. 2015;73:33-39. DOI: 10.1016/j.optlaseng.2015.04.002
  43. 43. Kumar M, Agarwal S, Kumar V, Khan GS, Shakher C. Experimental investigation on butane diffusion flames under the influence of magnetic field by using digital speckle pattern interferometry. Applied Optics. 2015;54(9):2450-2460. DOI: 10.1364/AO.54.002450
  44. 44. Kumar M, Shakher C. Holographic optical element based digital speckle pattern shearing interferometer. In: Proceedings of SPIE—The International Society for Optical Engineering (SPECKLE); 2018. Vol. 10834. p. 1083414. DOI: 10.1117/12.2318342
  45. 45. Kumar M, Birhman AS, Kannan S, Shakher C. Measurement of initial displacement of canine and molar in human maxilla under different canine retraction methods using digital holographic interferometry. Optical Engineering. 2018;57(9):094106. DOI: 10.1117/1.OE.57.9.094106
  46. 46. Kumar M, Shakher C. Experimental characterization of the hygroscopic properties of wood during convective drying using digital holographic interferometry. Applied Optics. 2016;55(5):960-968. DOI: 10.1364/AO.55.000960
  47. 47. Kumar M, Shakher C. Measurement of hygroscopic strain in deodar wood during convective drying using lensless Fourier transform digital holography. In: Optical Micro-and Nanometrology VI2016. Vol. 9890. International Society for Optics and Photonics. 2016. p. 989007. DOI: 10.1117/12.2227464
  48. 48. Kumar M. Some novel applications of digital speckle pattern interferometry and digital holographic interferometry (doctoral dissertation). Delhi: Indian Institute of Technology; 2016. Available from:
  49. 49. Torrence C, Compo GP. A practical guide to wavelet analysis. Bulletin of the American Meteorological Society. 1998;79(1):61-78
  50. 50. Liu H, Cartwright AN, Basaran C. Moire interferogram phase extraction: A ridge detection algorithm for continuous wavelet transforms. Applied Optics. 2004;43(4):850-857
  51. 51. Li S, Su X, Chen W. Wavelet ridge techniques in optical fringe pattern analysis. Journal of the Optical Society of America A. 2010;27(6):1245-1254
  52. 52. Daugman JG. Uncertainty relation for resolution in space, spatial frequency, and orientation optimized by two-dimensional visual cortical filters. Journal of the Optical Society of America A. 1985;2(7):1160-1169
  53. 53. Gabor D. Theory of communication. Electrical Engineers-Part I: General. Journal of the Institution of. 1947;94(73):58
  54. 54. Pritt MD, Ghiglia DC. Two-Dimensional Phase Unwrapping: Theory, Algorithms, and Software. New York: Wiley; 1998. ISBN-10: 0471249351
  55. 55. Bioucas-Dias JM, Valadao G. Phase unwrapping via graph cuts. IEEE Transactions on Image Processing. 2007;16(3):698-709
  56. 56. Wang Z, Bovik AC. A universal image quality index. IEEE Signal Processing Letters. 2002;9(3):81-84
  57. 57. Tounsi Y, Kumar M, Nassim A, Santoyo FM. Speckle noise reduction in digital speckle pattern interferometric fringes by nonlocal means and its related adaptive kernel-based methods. Applied Optics. 2018;57(27):7681-7690. DOI: 10.1364/AO.57.007681
  58. 58. Zada S, Tounsi Y, Kumar M, Mendoza-Santoyo F, Nassim A. Contribution study of monogenic wavelets transform to reduce speckle noise in digital speckle pattern interferometry. Optical Engineering. 2019;58(3):034109. DOI: 10.1117/1.OE.58.3.034109

Written By

Yassine Tounsi, Abdulatef Ghlaifan, Manoj Kumar, Fernando Mendoza-Santoyo, Osamu Matoba and Abdelkrim Nassim

Submitted: 14 March 2019 Reviewed: 10 June 2019 Published: 01 July 2019