Next Article in Journal
Inner Profile Measurement for Pipes Using Penetration Testing
Next Article in Special Issue
A Novel Method for Early Gear Pitting Fault Diagnosis Using Stacked SAE and GBRBM
Previous Article in Journal
An Improved Recognition Approach for Noisy Multispectral Palmprint by Robust L2 Sparse Representation with a Tensor-Based Extreme Learning Machine
Previous Article in Special Issue
A Switched-Element System Based Direction of Arrival (DOA) Estimation Method for Un-Cooperative Wideband Orthogonal Frequency Division Multi Linear Frequency Modulation (OFDM-LFM) Radar Signals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ultrasonic Flaw Echo Enhancement Based on Empirical Mode Decomposition

1
State Key Lab of Fluid Power and Mechatronic Systems, Zhejiang University, Hangzhou 310027, China
2
CRRC Zhuzhou Institute Co. Ltd., Zhuzhou 412001, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(2), 236; https://doi.org/10.3390/s19020236
Submission received: 3 December 2018 / Revised: 29 December 2018 / Accepted: 6 January 2019 / Published: 9 January 2019
(This article belongs to the Special Issue Sensor Signal and Information Processing II)

Abstract

:
The detection of flaw echoes in backscattered signals in ultrasonic nondestructive testing can be challenging due to the existence of backscattering noise and electronic noise. In this article, an empirical mode decomposition (EMD) methodology is proposed for flaw echo enhancement. The backscattered signal was first decomposed into several intrinsic mode functions (IMFs) using EMD or ensemble EMD (EEMD). The sample entropies (SampEn) of all IMFs were used to select the relevant modes. Otsu’s method was used for interval thresholding of the first relevant mode, and a window was used to separate the flaw echoes in the relevant modes. The flaw echo was reconstructed by adding the residue and the separated flaw echoes. The established methodology was successfully employed for simulated signal and experimental signal processing. For the simulated signals, an improvement of 9.42 dB in the signal-to-noise ratio (SNR) and an improvement of 0.0099 in the modified correlation coefficient (MCC) were achieved. For experimental signals obtained from two cracks at different depths, the flaw echoes were also significantly enhanced.

1. Introduction

The ultrasonic technique has been widely used in nondestructive testing. Usually, the backscattered signal is complex due to the existence of electronic noise and backscattering noise. Consequently, flaw echo detection may be challenging. Numerous methods have been proposed to enhance flaw echoes, such as split spectrum processing [1,2,3,4], wavelet transforms [5,6,7,8,9,10], the Stockwell transform [11,12,13,14], and empirical mode decomposition (EMD) [15,16,17,18,19,20,21,22] (including the so-called ensemble EMD, i.e., EEMD [23]).
Split spectrum processing has significant advantages in processing ultrasonic signals with scattered noise. Split-spectrum analysis separates the spectrum of the signals to obtain several sub-bands, and uses some nonlinear de-noising criteria (such as thresholding method, etc.) to process the signals in each sub-band to achieve the purpose of de-noising. The difficulty of split spectrum processing is how to determine the filter type, central frequency, bandwidth and other parameters. In addition, split spectrum processing lacks the capability of multiresolution analysis.
Wavelet transform is a classical multiresolution analysis method. The difficulty of wavelet transform is how to choose the appropriate wavelet base function and decomposition layer. The disadvantage of conventional wavelet transform is that the phase information of signals is lost.
S transform is the development of short-time Fourier transform and wavelet transform. S transform combines the multiresolution analysis ability of wavelet transform and the phase retention ability of short-time Fourier transform. Meanwhile, the S transform adopts Gaussian window function, which satisfies the normalization characteristic, so the S transform is invertible, that is, the original signal can be obtained from the converted time spectrum. However, due to the fact that the standard deviation of Gaussian window function in S transform is inversely proportional to the frequency and lacks flexibility, S transform may output the result of poor time-frequency resolution. At present, to compensate for the limitations of S transform, researchers introduce additional parameters to control the window function morphology, so that the generalized S transform has the ability to flexibly adjust the time-frequency resolution.
For a single component signal, Hilbert transform can be applied to obtain its analytical signal, and then the envelope spectrum and instantaneous frequency of the signal can be obtained. Detection of a flaw signal from the envelope spectrum is a common method of ultrasonic nondestructive testing. For a multi-component signal, it is necessary to decompose them into single component signals, and then obtain their analytical signals separately. Empirical mode decomposition (EMD) is an adaptive decomposition method. The original signal is decomposed into a series of intrinsic mode functions. The main disadvantages of EMD are the possibility of endpoint effect and mode mixing. Ensemble empirical mode decomposition (EEMD) can solve the mode mixing issue. EEMD decomposes the original signal by adding Gaussian white noise, and takes the result of the lumped average as the mode function. The main difficulty of EEMD is to select the intensity of Gaussian white noise and the times of lumped average.
In this article, an EMD-based methodology for ultrasonic flaw echo enhancement was established. The proposed methodology enhanced the flaw echo through six steps. First, the backscattered signal was adaptively decomposed into several intrinsic mode functions (IMFs) by EMD or EEMD. Second, the sample entropies (SampEn) [24,25] of all IMFs were calculated, and the differences in consecutive SampEn values were studied. Third, those IMFs with a large SampEn were considered to be irrelevant modes and were discarded, significantly suppressing the electronic noise. Fourth, the intervals containing the flaw echo were determined based on IMF interval thresholding and mode cell merging. Otsu’s method [26] was used to search the threshold in IMF interval thresholding. Fifth, a Turkey-Hanning window was used to separate the flaw echo for each relevant mode. Finally, the denoised signal was reconstructed by combining the separated flaw echoes and the residue.
The remainder of this article is organized as follows. In Section 2, reviews of the required tools, including EMD and EEMD, EMD-based denoising methods, SampEn, and Otsu’s method for thresholding, are presented. An analysis of the modes extracted from the backscattered signal, including the mixing of noise and flaw echoes, and the SampEn is given in Section 3. The proposed EMD-based methodology for ultrasonic flaw echo enhancement is given in Section 4 and tested using a simulated signal in Section 5. In Section 6, experimental validations of the proposed EEMD-based methodology are presented. Finally, conclusions are drawn in Section 7.

2. Required Tools

2.1. EMD and EEMD

In EMD, a signal is adaptively decomposed into a collection of IMFs. The EMD results can be presented as
x ( t ) = i = 1 L h ( i ) ( t ) + r ( t )
where x ( t ) is the observed signal, h ( i ) ( t ) ( i L ) are the extracted IMFs, and r ( t ) is the residue.
Unfortunately, EMD is susceptible to mode-mixing. EEMD is an effective technique for alleviating mode-mixing in EMD by repeatedly adding Gaussian white noise and finding the mean of the individual ensemble IMFs as the final IMF.

2.2. Denoising Strategies Based on EMD

Partial reconstruction, direct thresholding, and interval thresholding are three typical strategies adopted in EMD-based denoising.
Partial reconstruction removes noise from the observed signal by discarding irrelevant modes, which can be expressed as
x ^ ( t ) = i = M 1 L h ( i ) ( t ) + r ( t ) = x ( t ) i = 1 M 1 1 h ( i ) ( t )
where h ( i ) ( t ) ( i < M 1 ) are the irrelevant modes.
Direct thresholding is a direct application of wavelet thresholding in the EMD case. For the hard thresholding case, the denoised IMF is given by
h ˜ ( i ) ( t ) = { h ( i ) ( t ) , | h ( i ) ( t ) | > T i 0 , | h ( i ) ( t ) | T i
where T i is the threshold of h ( i ) ( t ) . The denoised signal can be given by
x ˜ ( t ) = i = M 1 M 2 h ˜ ( i ) ( t ) + i = M 2 + 1 L h ( i ) ( t ) + r ( t )
The interval thresholding divides an IMF into several mode cells and treats each mode cell as a whole to perform thresholding. Generally, a mode cell is defined as the signal between two adjacent zero-crossings. For the interval z j ( i ) = [ z j ( i ) z j + 1 ( i ) ] defined by two zero-crossings z j ( i ) and z j + 1 ( i ) , the denoised IMF in the hard thresholding case is given as
h ˜ ( i ) ( z j ( i ) ) = { h ( i ) ( z j ( i ) ) , | h ( i ) ( r j ( i ) ) | > T i 0 , | h ( i ) ( r j ( i ) ) | T i
where h ( i ) ( z j ( i ) ) are all of the samples from z j ( i ) to z j + 1 ( i ) , and h ( i ) ( r j ( i ) ) is the single extremum of h ( i ) ( t ) in the interval z j ( i ) .
Interval thresholding generally outperforms direct thresholding as it avoids catastrophic consequences for the continuity of the reconstructed signal, which are inevitable in direct thresholding.

2.3. Sample Entropy

Approximate entropy [27] and SampEn are two popular metrics for signal complexity measurement. Entropy values increase with increased signal complexity. It has been reported that SampEn outperforms approximate entropy in many aspects, such as reduced bias, independence from the signal, and relative consistency. SampEn was used here for signal complexity assessment.
The SampEn of a specified time series { u ( i ) , i = 1 , 2 , , N } can be obtained through the following steps.
Step 1: Form m dimensional vectors as
U m ( i ) = [ u ( i ) , u ( i + 1 ) , , u ( i + m 1 ) ]
where i = 1 , 2 , , N m + 1 .
Step 2: Define the distance of two such vectors:
d [ U m ( i ) , U m ( j ) ] = max k = 0 ( m 1 ) | u ( i + k ) u ( j + k ) |
Step 3: Consider the first N m vectors of length m so that for i = 1 , , N m , both U m ( i ) and U m + 1 ( i ) can be defined.
Step 4: Given a threshold r > 0 , define
B i m ( r ) = 1 N m 1 Θ ( r d [ U m ( i ) , U m ( j ) ] ) A i m ( r ) = 1 N m 1 Θ ( r d [ U m + 1 ( i ) , U m + 1 ( j ) ] )
where j = 1 , 2 , , N m , j i , and Θ ( · ) is defined as
Θ ( x ) = { 1 , x 0 0 , x < 0
Step 5: Calculate B m ( r ) and A m ( r ) :
B m ( r ) = 1 N m i = 1 N m B i m ( r ) A m ( r ) = 1 N m i = 1 N m A i m ( r )
Step 6: For a limited series, the SampEn (denoted by s ) is estimated as
s = S a m p E n ( m , r ) = ln A m ( r ) B m ( r )
In general, the dimension and threshold are often set to m = 2 and r = ( 0.1 0.25 ) S D u , where S D u is the standard deviation of the time series { u ( i ) } .

2.4. Otsu’s Method for Thresholding

Otsu’s method determines a threshold by maximizing the between-class variance σ B 2 . For a histogram with H levels (i.e., bins), the probability at each level can be first obtained:
p i = n i N , i = 1 , 2 , , H
where n i is the number of elements in the i th level, and N is the number of elements in the histogram. Obviously, p i 0 and p i = 1 are satisfied.
The histogram can be divided into two classes, C 1 and C 2 , with a threshold. σ B 2 is defined as
σ B 2 = ω 1 ( μ 1 μ 0 ) 2 + ω 2 ( μ 2 μ 0 ) 2
where
ω 1 = C 1 p i , ω 2 = C 2 p i , ω 1 + ω 2 = 1
μ 0 , μ 1 , μ 2 are the means of the histogram, class C 1 and class C 2 , respectively. Therefore, Equation (15) can be obtained:
μ 0 = ω 1 μ 1 + ω 2 μ 2
According to Equations (13)–(15), σ B 2 can also be expressed as
σ B 2 = ω 1 ω 2 ( μ 1 μ 2 ) 2

3. Analysis of Modes from Ultrasonic Signals

3.1. The Clutter Model

For metallic materials, when an incident ultrasonic wave propagates into the specimen, the backscattered signal will be primarily composed of three components: (1) the flaw echo signal s ( t ) , (2) backscattering noise v ( t ) due to the grains, and (3) electronic noise n ( t ) due to the instruments and the environment. n ( t ) can be approximated as Gaussian white noise. The frequency spectra of s ( t ) and v ( t ) can be expressed as [2]
V ( ω ) = H t 2 ( ω ) k = 1 K β k ω 2 x k e α s 2 x k ω 4 e i ω 2 x k c 0
S ( ω ) = H t 2 ( ω ) exp ( α s 2 d flaw ω ) exp ( i 2 d flaw c 0 )
where H t ( ω ) is the frequency response of the ultrasonic transducer, and d flaw is the location of the flaw. α s is the material attenuation coefficient, c 0 is the velocity of the longitudinal waves, and K is the total number of scatterers. β k and x k are the scattering coefficient and the position of the k th scatterer, respectively. The amplitudes of s ( t ) and v ( t ) are often normalized for brevity.
The similarity function obtained by deconvolution can be used to distinguish the flaw signals in ultrasonic inspection from other no flaw signals, such as specimen geometric reflection. It was found that the deconvolution patterns of the geometric reflection were impulse-like patterns, whereas those of flaws were bipolar patterns [28]. Therefore, geometric reflection was not considered in the model.
The observed backscattered signal x ( t ) can be expressed as
x ( t ) = s ( t ) + μ v ( t ) + σ n n 0 ( t )
where n 0 ( t ) is standard Gaussian white noise. μ and σ n are scale factors.
Typical simulated results are depicted in Figure 1. The centre frequency of the ultrasonic transducer was 5 MHz, and the sampling frequency was 100 MHz. The scale factors were set to μ = 0.3 and σ n = 0.2 . It can be seen in Figure 1a that the flaw echo had been polluted by intense noise. In addition, frequency aliasing of the flaw echo, backscattering noise, and wide-band Gaussian white noise can be found in Figure 1b.

3.2. Signal Decomposition

The EMD results for the observed signal x ( t ) are shown in Figure 2.
As shown in Figure 2, intense white noise was found in the low-order IMFs; in particular, IMF 1 resembled pure noise. The flaw echo was clearly detected in IMF 3. Consequently, the flaw echo was significantly enhanced by discarding IMF 1 and IMF 2 from the observed signal. However, partial reconstruction is often inadequate; further processing is required. Here, we take IMF 3 as an example.
The flaw echo is an instant signal. Ideally, oscillations can only be detected in the interval in which the flaw echo is located in IMF 3, in contrast to Figure 2. The difference mainly arises from frequency aliasing of the flaw echo, backscattering noise and wide-band Gaussian white noise. In addition, even in the noiseless case, the IMFs still contain false oscillations, as they resemble AM-FM modulated sinusoids. Consequently, further denoising is required to suppress the mixed noise lying outside the interval in which the flaw echo is located.

3.3. SampEn Values of IMFs

The SampEn values of all of the IMFs were calculated and are listed in Table 1, where the threshold r in the SampEn calculations was set to ( 0.1 , 0.15 , 0.2 ) S D u .
It can be noted that the SampEn tends to decrease with increasing IMF order, which indicates that the noise intensity in each IMF decreases as IMF order increases. It is noteworthy that the SampEn values of the first two IMFs, i.e., s 1 and s 2 , were much higher than the others. In addition, a sharp drop between s 2 and s 3 was detected.

4. Proposed Methodology

According to Section 3.2, flaw echo enhancement can be achieved by discarding irrelevant modes and suppressing the mixed noise in the remaining relevant modes. Specifically, the relevant modes are determined by the SampEn of all IMFs or their differences, and the mixed noise is suppressed by separating the flaw echo from those relevant modes by windowing.

4.1. Relevant Mode Selection

The IMFs with intense noise were of much higher complexity than other IMFs. Consequently, those relevant modes were determined according to the SampEn of the IMFs. Specifically, the parameter M 1 , which determined the first relevant mode, was determined by s i or the difference of s i .
Given a predefined SampEn threshold T s , M 1 can be determined as
M 1 = ( max i ) + 1 , s . t .   s i > T s
Using the difference of s i , i.e., d s i , M 1 can also be determined:
M 1 = arg min i { d s i } = arg min i { s i s i 1 }
To determine T s , we should note that the SampEn is dependent on the threshold r . For example, Table 1 shows that s 1 and s 2 decreased rapidly with increasing r . Consequently, the selection of T s was dependent on r . For this article, r = 0.15 S D u and T s = 1 were selected.

4.2. Mixed Noise Suppression

Two steps are required to suppress mixed noise: determine the location of the flaw echo in the first relevant mode and separate the flaw echo from all of the relevant modes.
For the first relevant mode, i.e., h ( M 1 ) ( t ) , a collection of mode cells is built according to the zero-crossings. The set of absolute values of the extrema in h ( M 1 ) ( t ) , which is denoted | r j ( i ) | = [ | r 1 ( i ) | , | r 2 ( i ) | , ] , can be determined. Otsu’s method is used for searching the threshold of | r j ( i ) | , i.e., T i . It is performed on the first relevant mode, and all of the mode cells in which the flaw echoes are located can be detected.
All of the adjacent mode cells in which the flaw echo is located are further merged, yielding an interval in which the flaw echo is located. For example, three adjacent mode cells, z j ( i ) = [ z j ( i ) z j + 1 ( i ) ] , z j + 1 ( i ) = [ z j + 1 ( i ) z j + 2 ( i ) ] and z j + 2 ( i ) = [ z j + 2 ( i ) z j + 3 ( i ) ] , can be merged into one interval [ z j ( i ) z j + 3 ( i ) ] . In other words, the interval in which the flaw echo is located is often composed of several mode cells containing the flaw echo.
Thus far, the location of the flaw echo in the first relevant mode has been detected. Next, the flaw echoes in all relevant modes will be separated.
A window was determined and used for separating the flaw echoes in the relevant modes. Herein, the Turkey-Hanning window was used. The window was defined as
ω ( m ) = { ( 1 cos 2 π m M + 1 ) 2 , 1 m M 2 1 , M 2 < m M 2 + W 1 cos 2 π m W M + 1 2 , M 2 + W < m M + W
where W determined the width of the pass zone and M determined the width of the transition zone. As an example, Figure 3 shows the waveform of a Turkey-Hanning window.
The parameters W and M were determined according to the width of the interval in which the flaw echo was located. To preserve the flaw echo, W was set to the width of this interval. M was selected more flexibly. In this article, M = [ W / 4 ] was adopted, where the operator “ [ ] ” indicated rounding.

4.3. Summarization of the Proposed Methodology

In this methodology, the flaw echo was enhanced in six steps.
Step 1. Signal decomposition. Decompose the observed signal x ( t ) into a collection of the residue r ( t ) and the IMFs’ h ( i ) ( t ) ( i = 1 , 2 , , L ) using EMD or EEMD.
Step 2. SampEn calculation. Obtain the SampEn of all of the IMFs.
Step 3. Relevant mode selection. Determine the first relevant mode using the SampEn or the differences between them.
Step 4. Determine the interval in which the flaw echo is located in the first relevant mode. Otsu’s method is used for threshold selection. Any two adjacent mode cells in which the flaw echo is located are merged into one interval.
Step 5. Separate the flaw echo using the Turkey-Hanning window in each relevant mode, which yields a collection of denoised modes h ˜ ( i ) ( t ) ( M 1 i L ).
Step 6. Reconstruct the denoised signal:
x ˜ ( t ) = i = M 1 L h ˜ ( i ) ( t ) + r ( t )

5. Simulated Signal Processing

5.1. Performance Assessment

To assess the performance of the EMD-based methodology, two metrics are introduced.
The signal-to-noise ratio (SNR) is the first metric. Suppose that s ( t ) is a noiseless signal; the SNR of signal x ( t ) is then defined as
S N R = 10 lg i = 1 N s 2 ( i ) i = 1 N ( x ( i ) s ( i ) ) 2
Well-preserved flaw echoes are expected in practice. Specifically, the amplitudes and shapes of flaw echoes are expected to be unchanged. This can be assessed by the modified correlation coefficient (MCC), which is defined as
M C C = | A x A s A s | ( 1 i = m n ( s ( i ) s ¯ ) ( x ( i ) x ¯ ) i = m n ( s ( i ) s ¯ ) 2 i = m n ( x ( i ) x ¯ ) 2 )
where m and n are instants defining the interval in which the flaw echo is located. A x and A s are the amplitudes of the flaw echoes in x ( t ) and s ( t ) , respectively.
A high SNR and low MCC are expected.

5.2. Signal Processing Results

The relevant modes were first determined. On one hand, with the predefined parameters r = 0.15 S D u and T s = 1 , it can be seen in Table 1 that s 1 , s 2 > T s and s 3 < T s and were satisfied. Consequently, M 1 = 3 was determined according to Equation (20). On the other hand, it could be inferred from Table 1 that d s 3 = s 3 s 2 < d s k ( k 3 ) was satisfied, which also indicates M 1 = 3 according to Equation (21). Thus, these two methods yielded the same relevant mode selection results. Consequently, the electronic noise could be significantly suppressed if we reconstructed the flaw echo via partial reconstruction by discarding h ( 1 ) ( t ) and h ( 2 ) ( t ) . The corresponding reconstructed signal, x ^ ( t ) , is depicted in Figure 4a.
h ( 3 ) ( t ) was used to determine the interval in which the flaw echo was located. The corresponding result is depicted in Figure 4b, for which the threshold given by Otsu’s method was 0.2902. The flaw echoes in all relevant modes h ( i ) ( t ) ( 3 i 8 ) were separated by the Turkey-Hanning window, and the denoised signal was reconstructed. The denoised signal is depicted in Figure 4c, which indicated that the electronic noise and backscattering noise were significantly suppressed, and the flaw echo were preserved well.
The performance of the EMD-based methodology was assessed using the SNR and the MCC. The results are provided in Table 2, which shows that the methodology achieved a high SNR and low MCC as expected. The improvements in SNR and MCC were 12.89 dB and 0.0099, respectively.

6. Experimental Study

The ultrasonic testing system was primarily composed of an ADVANTECH industrial personal computer (IPC), an Olympus ultrasonic probe with a centre frequency of 5 MHz, and a PCIUT3100 ultrasonic acquisition card installed on the IPC. A 6061 aluminium alloy specimen with two artificial cracks was used for ultrasonic testing. The specimen material density was 2.7 × 103 kg/m3 and the sound wave velocity was 6300 m/s. The two cracks, denoted F1 and F2, were machined using wire electrical discharge machining. The buried depths of F1 and F2 were 60 mm and 35 mm, respectively, and the sampling frequency was 100 MHz. The ultrasonic testing system and the specimen are shown in Figure 5.
The ultrasonic signals acquired from cracks F1 and F2 are shown in Figure 6. The backscattered signals from instant 1001 to instant 3000 were used for further analysis.
The processing results of the backscattered signal from crack F1 are shown in Figure 7. The backscattered signal was first decomposed into a collection of IMFs, as shown in Figure 7a. For brevity, only the first five IMFs are shown. Correspondingly, the SampEn values of the first five IMFs were calculated and are listed in Table 3. As the first relevant mode, the third IMF was used to determine the interval in which the flaw echo was located based on Otsu’s method for thresholding and merging of adjacent mode cells, as shown in Figure 7b. The denoised signal is shown in Figure 7c, indicating that both the backscattering noise and the electronic noise were significantly suppressed.
As supplementary information, Figure 8 shows the denoised signal of the backscattered signal from F2, with a significant reduction of the noise.

7. Conclusions

An EMD-based methodology was proposed for ultrasonic flaw echo enhancement. The observed signal was decomposed into IMFs by EMD or EEMD. The relevant modes were determined according to the SampEn of the IMFs. Otsu’s method was used for interval thresholding of the first relevant mode, obtaining the interval in which the flaw echo was located. The flaw echoes in all IMFs were separated by the Turkey-Hanning window. The separated flaw echo and the residue were added together, yielding the denoised signal. Simulation results demonstrated that the EMD-based methodology achieved a high SNR and low MCC. Applications of the EMD-based methodology in flaw echo enhancement in experimental signals was also presented.

Author Contributions

Conceptualization, W.F., X.Z. (Xiaojun Zhou) and C.Y.; formal analysis, W.F.; funding acquisition, C.Y.; investigation, W.F., X.Z. (Xiang Zeng) and C.Y.; methodology, W.F. and X.Z. (Xiang Zeng); resources, X.Z. (Xiaojun Zhou) and C.Y.; supervision, X.Z. (Xiaojun Zhou); writing—original draft, W.F.; writing—review and editing, W.F., X.Z. (Xiaojun Zhou), X.Z. (Xiang Zeng) and C.Y.

Funding

This research was funded by 1: The Fundamental Research Funds for the Central Universities under Grant No.2018QNA4001; 2: Zhejiang Provincial Natural Science Foundation of China under Grant No.LY18E050002.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shankar, P.M.; Karpur, P.; Newhouse, V.L.; Rose, J.L. Split-spectrum processing: Analysis of polarity threshold algorithm for improvement of signal-to-noise ratio and detectablity in ultrasonic signals. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1989, 36, 101–108. [Google Scholar] [CrossRef] [PubMed]
  2. Gustafsson, M.G.; Stepinski, T. Studies of split spectrum processing, optimal detection, and maximum likelihood amplitude estimation using a simple clutter model. Ultrasonics 1997, 35, 31–52. [Google Scholar] [CrossRef]
  3. Rodriguez, A.; Miralles, R.; Bosch, I.; Vergara, L. New analysis and extensions of split-spectrum processing algorithms. NDT E Int. 2012, 45, 141–147. [Google Scholar] [CrossRef]
  4. Jafar, S.; Daniel, T.N.; Kevin, D.D. Analysis of order statistic filters applied to ultrasonic flaw detection using split-spectrum processing. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1991, 38, 133–140. [Google Scholar] [CrossRef]
  5. Lazaro, J.C.; San Emeterio, J.L.; Ramos, A.; Fernandez-Marron, J.L. Influence of thresholding procedures in ultrasonic grain noise reduction using wavelets. Ultrasonics 2002, 40, 263–267. [Google Scholar] [CrossRef]
  6. Matz, V.; Smid, R.; Starman, S.; Kreidl, M. Signal-to-noise ratio enhancement based on wavelet filtering in ultrasonic testing. Ultrasonics 2009, 49, 752–759. [Google Scholar] [CrossRef] [Green Version]
  7. Peng, C.Y.; Gao, X.R.; Wang, A. Novel wavelet self-optimisation threshold denoising method in axle press-fit ultrasonic defect detection. Insight 2016, 58, 145–151. [Google Scholar] [CrossRef]
  8. Liu, Y.; Li, Z.; Zhang, W. Crack detection of fibre reinforced composite beams based on continuous wavelet transform. Nondestruct. Test. Eval. 2010, 25, 25–44. [Google Scholar] [CrossRef]
  9. Chen, Y.C.; Zhou, X.J.; Yang, C.L.; Li, Z. The ultrasonic evaluation method for the porosity of variable-thickness curved CFRP workpiece: Using a numerical wavelet transform. Nondestruct. Test. Eval. 2014, 29, 195–207. [Google Scholar] [CrossRef]
  10. Praveen, A.; Vijayarekha, K.; Abraham, S.T.; Venkatraman, B. Signal quality enhancement using higher order wavelets for ultrasonic TOFD signals from austenitic stainless steel welds. Ultrasonics 2013, 53, 1288–1292. [Google Scholar] [CrossRef]
  11. Stockwell, R.G.; Mansinha, L.; Lowe, R.P. Localization of the complex spectrum: The S transform. IEEE Trans. Signal Process. 1996, 44, 998–1001. [Google Scholar] [CrossRef]
  12. Ari, S.; Das, M.K.; Chacko, A. ECG signal enhancement using S-Transform. Comput. Biol. Med. 2013, 43, 649–660. [Google Scholar] [CrossRef] [PubMed]
  13. Muhammad, A.M.; Jafar, S. S-transform applied to ultrasonic nondestructive testing. In Proceedings of the IEEE Ultrasonics Symposium, Beijing, China, 2–5 November 2008; IEEE: Piscataway, NJ, USA, 2008; pp. 184–187. [Google Scholar] [CrossRef]
  14. Benammar, A.; Drai, R.; Guessoum, A. Ultrasonic flaw detection using threshold modified S-transform. Ultrasonics 2014, 54, 676–683. [Google Scholar] [CrossRef] [PubMed]
  15. Huang, N.E.; Shen, Z.; Long, S.R.; Wu, M.C.; Shih, H.H.; Zheng, Q.N.; Yen, N.C.; Tung, C.C.; Liu, H.H. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proc. R. Soc. A-Math. Phys. Eng. Sci. 1998, 454, 903–995. [Google Scholar] [CrossRef]
  16. Gabriel, R.; Patrick, F. One or Two Frequencies? The empirical mode decomposition answers. IEEE Trans. Signal Process. 2008, 56, 85–95. [Google Scholar] [CrossRef]
  17. Kazys, R.; Tumsys, O.; Pagodinas, D. Ultrasonic detection of defects in strongly attenuating structures using the Hilbert–Huang transform. NDT E Int. 2008, 41, 457–466. [Google Scholar] [CrossRef]
  18. Abdel, O.B.; Jean, C.C. EMD-based signal filtering. IEEE Trans. Instrum. Meas. 2007, 56, 2196–2202. [Google Scholar] [CrossRef]
  19. Zhang, Q.; Que, P.W.; Liang, W. Applying sub-band energy extraction to noise cancellation of ultrasonic NDT signal. J. Zhejiang Univ. Sci. A 2008, 9, 1134–1140. [Google Scholar] [CrossRef]
  20. Yannis, K.; Stephen, M. Development of EMD-based denoising methods inspired by wavelet thresholding. IEEE Trans. Signal Process. 2009, 57, 1351–1362. [Google Scholar] [CrossRef]
  21. Yang, G.L.; Liu, Y.Y.; Wang, Y.Y.; Zhu, Z.L. EMD interval thresholding denoising based on similarity measure to select relevant modes. Signal Process. 2015, 109, 95–109. [Google Scholar] [CrossRef]
  22. Bouden, T.; Djerfi, F.; Nibouche, M. Adaptive split spectrum processing for ultrasonic signal in the pulse echo test. Russ. J. Nondestruct. Test. 2015, 51, 245–257. [Google Scholar] [CrossRef]
  23. Sharma, G.K.; Kumar, A.; Jayakumar, T.; Rao, B.P.; Mariyappa, N. Ensemble empirical mode decomposition based methodology for ultrasonic testing of coarse grain austenitic stainless steels. Ultrasonics 2015, 57, 167–178. [Google Scholar] [CrossRef] [PubMed]
  24. Richman, J.S.; Moorman, J.R. Physiological time-series analysis using approximate entropy and sample entropy. Am. J. Physiol. Heart Circ. Physiol. 2000, 278, H2039–H2049. [Google Scholar] [CrossRef]
  25. Hu, X.S.; Jiang, J.C.; Cao, D.P.; Bo, E. Battery Health Prognosis for Electric Vehicles Using Sample Entropy and Sparse Bayesian Predictive Modeling. IEEE Trans. Ind. Electron. 2016, 63, 2645–2656. [Google Scholar] [CrossRef]
  26. Otsu, N. A threshold selection method from gray-level histograms. IEEE Trans. Syst. Man Cybern. 1979, 9, 62–66. [Google Scholar] [CrossRef]
  27. Steven, M.P. Approximate entropy as a measure of system complexity. Proc. Natl. Acad. Sci. USA 1991, 88, 2297–2301. [Google Scholar] [CrossRef]
  28. Jung, H.J.; Kim, H.J.; Song, S.J.; Kim, Y.H. Model-based enhancement of the TIFD for flaw signal identification in ultrasonic testing of welded joints. In Proceedings of the 29th Annual Review of Progress in Quantitative Nondestructive Evaluation, Western Washington University, Bellingham, WA, USA, 14–19 July 2002; Amer Inst Physics: College Park, MD, USA, 2002; pp. 628–634. [Google Scholar] [CrossRef]
Figure 1. Simulated results. (a) Waveforms and (b) frequency spectra.
Figure 1. Simulated results. (a) Waveforms and (b) frequency spectra.
Sensors 19 00236 g001
Figure 2. Intrinsic mode functions (IMFs) extracted. (a) IMF 1 ~ IMF 4; (b) IMF 5 ~ IMF 8.
Figure 2. Intrinsic mode functions (IMFs) extracted. (a) IMF 1 ~ IMF 4; (b) IMF 5 ~ IMF 8.
Sensors 19 00236 g002
Figure 3. Waveform of a Turkey-Hanning window.
Figure 3. Waveform of a Turkey-Hanning window.
Sensors 19 00236 g003
Figure 4. Processing results of the simulated signal. (a) The reconstructed signal given by partial reconstruction; (b) the detection interval in which the flaw echo was located; (c) the denoised signal.
Figure 4. Processing results of the simulated signal. (a) The reconstructed signal given by partial reconstruction; (b) the detection interval in which the flaw echo was located; (c) the denoised signal.
Sensors 19 00236 g004
Figure 5. The ultrasonic test system and specimen. (a) The ultrasonic testing system. (b) The specimen. (c) Geometry of the specimen. (d) Schematic diagram of flaws location and size, with dimensions in millimeters.
Figure 5. The ultrasonic test system and specimen. (a) The ultrasonic testing system. (b) The specimen. (c) Geometry of the specimen. (d) Schematic diagram of flaws location and size, with dimensions in millimeters.
Sensors 19 00236 g005
Figure 6. The ultrasonic signals. (a) Acquired from F1; (b) acquired from F2.
Figure 6. The ultrasonic signals. (a) Acquired from F1; (b) acquired from F2.
Sensors 19 00236 g006
Figure 7. Processing results of the backscattered signal from F1. (a) The IMFs extracted; (b) the interval in which the flaw echo was located; (c) the denoised signal.
Figure 7. Processing results of the backscattered signal from F1. (a) The IMFs extracted; (b) the interval in which the flaw echo was located; (c) the denoised signal.
Sensors 19 00236 g007
Figure 8. The denoised signal of the ultrasonic signal from F2.
Figure 8. The denoised signal of the ultrasonic signal from F2.
Sensors 19 00236 g008
Table 1. SampEn of the IMFs.
Table 1. SampEn of the IMFs.
r/SDu0.10.150.2
s12.01211.79431.535
s21.88231.48721.213
s30.72240.62420.5682
s40.6460.59840.5616
s50.55010.49210.4472
s60.47170.35890.2708
s70.10950.08210.065
s80.05260.0360.0276
Table 2. Performance assessment.
Table 2. Performance assessment.
SNR/dBMCC/ × 10−2
observed−3.471.05
partial reconstruction denoised1.81
9.42
0.06
0.06
Table 3. SampEn of the first five IMFs.
Table 3. SampEn of the first five IMFs.
s1s2s3s4s5
1.49021.17390.56840.55080.4764

Share and Cite

MDPI and ACS Style

Feng, W.; Zhou, X.; Zeng, X.; Yang, C. Ultrasonic Flaw Echo Enhancement Based on Empirical Mode Decomposition. Sensors 2019, 19, 236. https://doi.org/10.3390/s19020236

AMA Style

Feng W, Zhou X, Zeng X, Yang C. Ultrasonic Flaw Echo Enhancement Based on Empirical Mode Decomposition. Sensors. 2019; 19(2):236. https://doi.org/10.3390/s19020236

Chicago/Turabian Style

Feng, Wei, Xiaojun Zhou, Xiang Zeng, and Chenlong Yang. 2019. "Ultrasonic Flaw Echo Enhancement Based on Empirical Mode Decomposition" Sensors 19, no. 2: 236. https://doi.org/10.3390/s19020236

APA Style

Feng, W., Zhou, X., Zeng, X., & Yang, C. (2019). Ultrasonic Flaw Echo Enhancement Based on Empirical Mode Decomposition. Sensors, 19(2), 236. https://doi.org/10.3390/s19020236

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop