Temperature Measurement Method Based on Riesz Transform Method

A method to measure the temperature of a metal plate is presented using a Riesz transform method and the monogenic signal to extract the optical phase distribution from a fringe pattern from which one can get the unknown temperature. The performance of this method is evaluated by the RFSIM metric obtained from the 2nd-order Riesz transform. A phase distribution with a good accuracy is provided.


Introduction
The measurement of temperature field is of great significance to scientific research and national economic development.This is because the measurement of thermal physical quantities is widely used in the implementation and control of industrial processes such as power, aerospace, chemical industry, and oil refining, among others [1].Temperatures are usually measured by inserting thermocouples into tested fields and entire temperature fields are reconstructed using readings from the thermocouples.However, this method has drawbacks.Temperatures at different points are not obtained simultaneously and tested fields are disturbed by the thermocouples themselves.In addition, obtained readings need to be compensated against the effect brought about by radiation.
Digital speckle pattern interferometry is a whole field optical method for noncontact and nondestructive surface analysis.It is now considered as a powerful tool for industrial measurements.It enables full-field measurement of optical phase changes via the acquisition of speckle patterns [7][8][9].
After acquisition, a simple subtraction is usually performed to obtain a correlation fringe pattern.The greatest challenges in speckle interferometry focus on relating fringe patterns to phase mapping, permitting the direct determination of surface displacement.However, as DSPI fringes are characterized by a strong speckle noise background, a denoising method [10][11][12] must be used before the phase evaluation.The removal of speckle noise in DSPI is a complex problem which has generated an active area of research.
The aim of this study is the application of the Riesz transform and the monogenic signal to extract the optical phase distribution from a single fringe pattern without the step of speckle denoising from which one can get the unknown temperature.Felsberg and Sommer [13] proposed two-dimensional generalization of the analytic signal called monogenic signal.This representation based on the Riesz transform [14] instead of the Hilbert transform preserves the core properties of the analytic signal that decomposes an image into information about structure.The monogenic signal is characterized by its three independent sizes as local amplitude, local phase, and local orientation.It has proved its usefulness in several applications [15][16][17].
The paper is organized as follows.Section 2 is devoted to the description of the monogenic phase extraction technique and the RFSIM metrics.In Section 3, a brief description of speckle correlation fringe obtained by digital speckle pattern interferometry (DSPI) technique for temperature measurement is presented.The numerical simulations and results of temperature evaluation are presented in Section 4. Finally a conclusion is given in Section 5.

The Riesz Transform
The Riesz transform for an image  is defined by where * stands for the convolution and (, ) is the image to be analyzed.  and   are the spatial representation of the Riesz kernel, respectively, defined by ( The two output images  1 and  2 are the Riesz transform components on -axis and -axis.The Riesz transform is the natural 2D extension of the Hilbert transform.
2.1.The Monogenic Signal.For an image, the monogenic signal is defined as the combination of the image  and its Riesz transform components: From this formulation, the monogenic signal is 3D vector and we can get the local features as local amplitude, local orientation, local phase, and phase vector.The local amplitude is defined by The local orientation with  ∈ [−/2, /2] is defined by The local phase with  ∈ [−, ] is defined by It is seen that the local orientation, the local phase, and the local amplitude of the image can be conveniently represented via the Riesz transform based on the monogenic signal which decomposes an image into information at each location.
Due to the fact that structural information is related to a given orientation, the phase vector is obtained by the local phase associated with a given local orientation.The phase vector can be interpreted as a rotation vector, in which magnitude corresponds to the phase angle between the image and the monogenic signal.

RFSIM Metric.
The RFSIM metric (Riesz transform based feature similarity) [18] allows judgeing the image quality assessment.This metric is computed by comparing 2nd-order Riesz transform features at key locations between two images (an image with its reference).
The 2nd-order Riesz transform is obtained from the 1storder Riesz transform components  1 and  2 by with Let ( 1 ,  2 ,  3 ,  4 ,  5 ), the first two are 1st-order Riesz transform coefficients and the last three are 2nd-order coefficients of an image .Similarly, let ( 1ref ,  2ref ,  3ref ,  4ref ,  5ref ): the first two are 1st-order Riesz transform coefficients and the last three are 2nd-order coefficients of the reference image  ref .
The coefficient RFSIM is defined as where   is the feature similarity between the two images given by where   , the distance of comparison at the corresponding location (, ) between the two features   and  ref ( = 1, . . ., 5), is formulated by The feature mask (, ) is defined as the edge locations of the each image by where ⊕ is the logical "OR" operation. 1 denotes the result of edge detection performed on the image  and  2 the result of edge detection on the reference image  ref .
The RFSIM metric quantifies how well image details are preserved.

Correlation Fringe Pattern for Temperature Measurement
3.1.Correlation Fringe Pattern Formation.The basic principle DSPI technique is that the speckle pattern intensity distribution is a function of the relative phases of two interfering plane waves inside each resolution cell of an imaging setup.Displacement of the surface affects the intensity received in each speckle cell on the image.The intensity distribution of a reference speckled image (before displacement) is where (, ) is the bias intensity, (, ) is the visibility, and   (, ) is the original phase from the speckle that appears as the high frequency and apparently random pixel-by-pixel intensity variation.After displacement, the intensity distribution becomes  2 (, ) =  (, ) +  (, ) cos (  (, ) + Δ (, )) , where Δ is the phase change in the light resulting from the displacement.The speckle fringe patterns are obtained by subtraction of a reference speckled image from image of displaced surface.In the subtraction mode, the squared   difference between the fringe pattern before and after deformation leads to a secondary fringe pattern which is very sensitive to microstructural changes of the studied surface; if the deformation is too large, the fringe becomes too dense and may vanish due to speckle decorrelation.Speckle decorrelation limits the range of deformation that can be recorded by speckle interferometry.The effects of decorrelation are a distortion and a corresponding reduction of fringe visibility in the observation plane.We assume that the displacements are sufficiently small that speckle decorrelation effects can be ignored.
The intensity distribution in the speckle correlogram is given by  (, ) =  2 (, ) −  1 (, ) , ( 16) × sin (  (, ) + Δ (, ) 2 ) . ( Fringes obtained are characterized by a speckle noise background; a noise reduction method must be used before the phase distribution is evaluated.The desired information containing the sin(Δ/2) fringe term (shape of envelope modulating the random speckle term) may be rectified and filtered by appropriate computer image processing in order to remove the high frequency sin(  + Δ/2) noise.Squaring ( Since   (, ) changes rapidly across the speckle pattern, the ensemble average of the second sine squared term in (18) across the whole measurement area leads to Substituting ( 19) into (18) yields Equation ( 20) is equivalent to the classical interferometry equation.It describes perfect correlated fringes without speckles.In practice, due to speckle decorrelation, fringe visibility is always <1. Fringes visualization on a monitor is sufficient for qualitative investigation.For quantitative evaluation, the temperature information has to be numerically extracted.
where  is the wavelength.The change in the refractive index of air is usually related to a temperature variation through a factor approximately constant; if the temperature changes are small, the so-called Gladstone-Dale constant.The refractive index of air at 632.8 nm can be evaluated with greater precision by [19]  − 1 = 0.292015 ⋅ 10 −3 (1 + 0.368184 ⋅ 10 −2 ) , (1 + 0.368184 ⋅ 10 Owing to the linear relationship between Δ and Δ, measurements of  may be made automatically from the fringe pattern phase distribution using the monogenic signal phase vector.

Numerical Simulation.
The numerical simulation consists in generating digitally a fringe pattern to verify the ability of the method to determine the phase distribution.The test phase function we used is The intensity distribution of fringe pattern shown in Figure 1(a) is  (, ) = 1 + 0.5 cos  (, ) .
We illustrated the Riesz component of the fringe pattern on -axis and on -axis in Figures 1(b) and 1(c), respectively.The wrapped retrieved phase from Riesz Transform method is plotted in Figure 2. The phase-unwrapping process which provides a continuous phase distribution over its definition domain implies suitable mathematical techniques [20], which are then needed.The unwrapped procedure consists of finding the correct field number for each phase measurement.In this work we used a fast phase unwrapping algorithm [21] that is largely immune to the presence of noise and residues and does not require strong user input.The difference between the unwrapped test phase and the unwrapped retrieved phase is depicted in Figure 3.The obtained RFSIM coefficient is 0.8250.

Temperature Evaluation Method.
To validate this phase evaluation method, real speckle correlation fringe pattern is used.The fringes pattern above a central section of a horizontal aluminum plate (0.7 cm × 25 cm × 25 cm) is shown in Figure 4.The correlation fringe is digitized in 512 × 512 × 8 bit data, with a magnification factor such that 1 pixel = 40.10−6 m.
The Riesz components on -axis and on -axis of the correlation fringe relative to a central section of a horizontal aluminum plate are presented in Figure 5.The wrapped and unwrapped retrieved phases are depicted in Figures 6(a) and 6(b).The obtained results are in very close agreement with those obtained by the phase-shifting method.The obtained RFSIM coefficient is 0.5302.
The temperature variation with distance from the plate as calculated from the corresponding unwrapped phase is illustrated in Figures 7 and 8 which show the temperature variation with distance from the plate obtained by [19].

Conclusion
In this work, we present a Riesz transform method for the temperature evaluation of a metal plate.The technique is based on optical phase evaluation method from which one can get the unknown temperature.The Riesz transform method provides phase distribution from a single fringe pattern with a good accuracy.This method is adequate for real time acquisition and easy to implement.We made a comparison with the results of Spagnolo and Al; we obtained results that agree relatively well with the margin of experimental error which was about 0.4 K near the plate surface and 0.57 K at a distance of about 3 cm from the plate.We note that our results render fairly experimental reality knowing that our method is relatively simple and inexpensive.

Figure 1 :
Figure 1: Simulation results: (a) fringe pattern, (b) Riesz component of the fringe pattern on -axis, and (c) Riesz component of the fringe pattern on -axis.

Figure 2 :
Figure 2: The wrapped retrieved phase from Riesz transform method.

Figure 4 :
Figure 4: Correlation fringes pattern relative to the central section of the plate.The temperature of the plate was 325 ± 0.2 K. Room temperature was 295 ± 0.2 K.

Figure 5 :Figure 6 :
Figure 5: (a) Riesz component of the fringe pattern on -axis; (b) Riesz component of the fringe pattern on -axis.

3. 2 .
Temperature Measurement.Now, let us consider a light beam through a test section of length .The unwrapped phase is related to the variation of the refractive index of air on top of the plate surface by

where
is in Celsius.This equation is based on the Gladstone-Dale is relation, with wavelength dependence calculated according to Meggers and Peters and small corrections due to Tilton.Hence

Figure 7 :
Figure 7: Temperature variation with distance from the plate.

Figure 8 :
Figure 8: Temperature variation with distance from the plate obtained by Spagnolo and Al.