A Wavelet Threshold Function for Treatment of Partial Discharge Measurements

*Caio F.F.C. Cunha, Mariane R. Petraglia, André T. Carvalho and Antonio C.S. Lima*

#### **Abstract**

Based on the wavelet transform filtering theory, the chapter will describe the elaboration of a wavelet threshold function intended for the denoising of the partial discharge phenomenon measurements. This new function, conveniently named Fleming threshold, is based on the logistic function, which is well known for its utility in several important areas. In the development is shown some variations in the application of the Fleming function, in an attempt to identify the decomposition levels where the thresholding process must be more stringent and those where it can be more lenient, which increases its effectiveness in the removal of noisy coefficients. The proposed function and its variants demonstrate excellent results compared to other wavelet thresholding methods already described in the literature, including the famous Hard and Soft functions.

**Keywords:** wavelet transform, threshold function, partial discharge, signal denoising

#### **1. Introduction**

The analysis of the Partial Discharges (PD) phenomenon, which manifests itself in the existing imperfections into the insulation of high voltage equipment have received global acceptance as an important tool for the predictive diagnosis of the operational conditions of these, allowing taking measures that can safeguard both the material and the power supply quality of the electrical system.

PD are short duration impulsive signals and, consequently, these can be detected in a wide frequency range, from a few kHz to GHz. Normally, there is a direct relationship between the frequency range where there is a higher incidence of PD pulses and the type of high-voltage equipment evaluated, e.g. transformers and generators usually emit pulses from a few tens of kHz up to about 30 MHz [1], whereas Gas Insulated Substations (GIS) are affected by very fast pulses ranging from 300 MHz to 3 GHz and for cables the spectrum covers frequencies from 300 MHz to 1 GHz.

**Figure 1** shows two examples of measured PD pulses in two different HV equipment, a GIS and a hydro generator. Note the marked white noise presence. The pulses normally have an exponentially damped oscillatory shape or only an exponentially damped shape [2].

and Amplitude Modulated (AM) noise. The results showed that our proposal is able to

Noise degrades the accuracy and precision of analysis, in addition to reducing the detection limit of the instrument applied in the PD measurements. Often the WT is a tool designed to attenuate continuous random noise (white noise), because after the decomposition of a signal in the wavelet domain can be noted that the average density of the coefficients is inversely proportional to the dyadic scale 1*=*2 *<sup>j</sup>* (*j* indicates the level of decomposition), i.e., half of the number of extreme local coefficients do not spread from a 1*=*2 *<sup>j</sup>* scale to the next 1*=*2 *<sup>j</sup>*þ<sup>1</sup> scale, distributing it uniformly across the scales. As the wavelet coefficients distribution pattern of the PD signal (which tends to have its energy concentrated in few decomposition levels) differs from the noise pattern, it becomes easy to identify and separate the PD signals from the noise [30–32]. However, in wavelet denoising, the noise attenuation occurs not only to the white noise but also to the noise with frequency components that do not match the frequency components of the PD pulses.

Basically, the wavelet shrinkage denoising process involves three steps [13, 31]:

1.Determine the WT decomposition tree (discrete WT, wavelet packet transform, stationary WT or dual-tree complex WT) to be applied, the number of decomposition levels *J* and the wavelet function that will be employed on each of the *j* levels (where j ¼ 1, 2, ⋯, *J*), and then perform the

decomposition of the analyzed signal into its wavelet coefficients;

3.Reconstruct the signal by applying the Inverse Wavelet Transform

2.Calculate the threshold values using one of the threshold selection rules, which depend on statistical estimation of the noise level present in the signal. Apply the calculated value in a threshold function to thus reduce the coefficients of the noise figure and preserve the signal coefficients of interest, in our case the PD pulse;

(corresponding to the decomposition tree selected in the first step) in the threshold coefficients, to obtain the filtered signal in the time domain.

overcome, qualitatively, and quantitatively, all the confronted functions.

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements*

**2. Wavelet domain detection**

*DOI: http://dx.doi.org/10.5772/intechopen.94115*

**Figure 2.**

**313**

*Signal denoising steps by wavelet transform.*

**Figure 1.** *PD pulses measured in: (a) GIS; (b) Hydro generator.*

The proper diagnosis of equipment is closely related to the peak amplitude and shape of the pulses detected. Therefore, it is important to preserve the amplitude characteristics of the signal (especially the peaks), providing higher Signal to Noise Ratio (SNR) and lower Amplitude Error (EA).

The application of FFT and STFT filtering is not as effective in the treatment of non-stationary, transient, and stochastic signals as the PD [3], since these transforms do not allow a location in the time and frequency domain in the same way as the wavelet transform does [4, 5] (with better resolution in frequency and worse resolution in time for the low frequency components of the signal; and worse resolution in frequency and better resolution in time for the high frequency components of the signal). Therefore, the performance of these methods becomes limited in comparison with the wavelet denoising, which presents a capacity of self-adaptation to the signal.

Partial discharges, almost entirely, are electrically detected and quantified, exposing them to the extensive noise interferences that may compromise the PD signals measurement, limiting the diagnosis accuracy. Different signal processing tools have been used to extract the PD signals from these noise sources; among them, it is possible to highlight the Wavelet Transform (WT). The filtering by wavelet processing is recommended in the extraction of PD signals immersed in Gaussian noise [1, 6].

An efficient application of wavelet processing depends on the careful selection of the parameters that will concentrate the coefficients on the most suitable decomposition levels to minimize the PD signal information loss. Among these, we have the applied WT, the number of decomposition levels, the wavelet functions used in each of these levels, the method of estimating the threshold value of the obtained coefficients and the threshold function.

The choice of most of these parameters has already been widely explored in several works meant to PD processing [2, 6–13]. However, with respect to threshold functions, most studies do not focus on PD signals denoising but on audio signals [14–18], Electrocardiogram (ECG) [19–21] or images [22–28]. Therefore, there is a lack of a dedicated threshold function to improve the PD pulses denoising, in order to increase the precision in the diagnosis of High Voltage (HV) equipment.

Based on the WT filtering theory, this chapter will be described the development of a wavelet threshold function aiming to improve the noise reduction in PD measurements. The logistic function serves as an inspiration to this new function [29], which is well known for its usefulness in numerous areas. Since it is customary to associate functions of this type with something that refers them to the name of their developers (e.g., [25]), it was designated as Fleming threshold function.

The denoising performance of the proposed threshold function was compared with the traditional Hard and Soft functions and with twelve other thresholding functions. For a fair analysis of the filtering results, were used 2064 simulated and measured PD pulses contaminated with uniform white noise, Gaussian white noise, and Amplitude Modulated (AM) noise. The results showed that our proposal is able to overcome, qualitatively, and quantitatively, all the confronted functions.

## **2. Wavelet domain detection**

Noise degrades the accuracy and precision of analysis, in addition to reducing the detection limit of the instrument applied in the PD measurements. Often the WT is a tool designed to attenuate continuous random noise (white noise), because after the decomposition of a signal in the wavelet domain can be noted that the average density of the coefficients is inversely proportional to the dyadic scale 1*=*2 *<sup>j</sup>* (*j* indicates the level of decomposition), i.e., half of the number of extreme local coefficients do not spread from a 1*=*2 *<sup>j</sup>* scale to the next 1*=*2 *<sup>j</sup>*þ<sup>1</sup> scale, distributing it uniformly across the scales. As the wavelet coefficients distribution pattern of the PD signal (which tends to have its energy concentrated in few decomposition levels) differs from the noise pattern, it becomes easy to identify and separate the PD signals from the noise [30–32]. However, in wavelet denoising, the noise attenuation occurs not only to the white noise but also to the noise with frequency components that do not match the frequency components of the PD pulses.

Basically, the wavelet shrinkage denoising process involves three steps [13, 31]:


**Figure 2.** *Signal denoising steps by wavelet transform.*

**Figure 2** illustrates each of these steps involved in the wavelet denoising processing of a digitalized signal.

As several parameters are involved, they should be carefully selected according to the signal characteristics, in order to maximize their wavelet coefficients above the noise level. Thus, filtering performance is closely related to each of these parameters and some of these will have a greater influence on the quality of the result. The determination of these parameters shows to be an optimization challenge [33]. In this chapter, we will focus our attention on the improvement of the threshold function applied in the second step.

probability, statistics, chemistry, etc. Eq. (1) represents the logistic function, where *H* is the maximum value of the curve, *α* controls the slope of this curve and *x* corresponds to the value of *x* at the midpoint of the sigmoid curve dictated by the numerator value (**Figure 3**). When the *x* value tends to þ∞ the curve approaches *H*

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements*

*f x*ð Þ¼ *<sup>H</sup>*

1 þ *e*�*α*ð Þ *<sup>x</sup>*�*x*<sup>0</sup>

Eq. (1) enabled us to develop the threshold function for filtering signals in such a way that it circumvents the problems previously described. As the objective in the thresholding process is to preserve the coefficients above the threshold value, it is easy to see that the maximum value *H* will be the decomposed wavelet coefficient *wj*,*<sup>k</sup>* (which corresponds to the variable *x*). Thus, the function will maintain symmetry when we vary the inclination constant *c* of the curve. Finally, it is necessary to move the function along the abscissa axis so that the graph shown in **Figure 2** leaves the ordinate axis and stays over the threshold value, this is done by

subtracting the variable *x* from the value where we want to move the function (*x*0), i.e., the value of the coefficients *wj*,*<sup>k</sup>* must be subtracted from the threshold value *λ*.

For a more efficient implementation, in which it is not necessary to worry about

�*<sup>c</sup>*ð Þ *<sup>ω</sup>j*,*k*�*<sup>λ</sup>* <sup>¼</sup> *<sup>ω</sup><sup>j</sup>*,*<sup>k</sup>*

1 þ *e*

*<sup>c</sup>*ð Þ �*ωj*,*k*þ*<sup>λ</sup>* (2)

*<sup>c</sup>*ð Þ ð Þ �*sign*ð Þ *<sup>ω</sup>j*,*<sup>k</sup>* �*ωj*,*<sup>k</sup>* <sup>þ</sup>*<sup>λ</sup>* (3)

By making these adaptations, we obtain the following threshold function:

1 þ *e*

*<sup>η</sup> <sup>f</sup> <sup>ω</sup><sup>j</sup>*,*<sup>k</sup>*, *<sup>λ</sup>*,*<sup>c</sup>* <sup>¼</sup> *<sup>ω</sup><sup>j</sup>*,*<sup>k</sup>*

wavelet denoising in relation to the Hard and the Soft functions.

the fact that the coefficient *wj*,*<sup>k</sup>* is positive or negative, Eq. (2) can be rewritten using the signum (*sign*) function, which returns þ1 if the value is positive and �1 if

1 þ *e*

With high *c* values, the curve inclination on the threshold point is such that it approaches the Hard function, but with a smoother (differentiable) transition. For low *c* values, the inclination of the function will act with less intensity on the coefficients below the threshold value and with greater intensity on the coefficients above this value, i.e., a large part of noisy coefficients may pass and there will be information losses on those coefficients that represent the signal of interest, in our case the PD pulse. With the appropriate choice of the *c* value for each processed signal, it is possible to obtain a significant improvement in the result of the PD

Considering a threshold value *λ* ¼ 1, **Figure 4** shows the behavior of the proposed threshold function (called the Fleming function) for different *c* values (3, 10, 20, 30, 50, 80, 100 e 200). With very low values of *c* there is the possibility of

*<sup>η</sup> <sup>f</sup> <sup>ω</sup><sup>j</sup>*,*<sup>k</sup>*, *<sup>λ</sup>*,*<sup>c</sup>* <sup>¼</sup> *<sup>ω</sup><sup>j</sup>*,*<sup>k</sup>*

the value is negative. Thus, we have:

**Figure 3.** *Logistic function.*

**315**

(1)

and when it tends to �∞ it approaches to zero.

*DOI: http://dx.doi.org/10.5772/intechopen.94115*

#### **3. Fleming threshold function**

In most of the wavelet denoising literature, especially those focused in the treatment of PD signals, the choice of the threshold function normally falls between the Hard and the Soft functions. Moreover, it is well known that for PD pulse filtering the Hard function tends to preserve more of the signal information, providing a higher SNR and a lower Amplitude Error (AE). However, the Hard estimate has discontinuities, being not differentiable, which ends up causing instability problems and sensitivity to small changes in the data pseudo-Gibbs effect. The Soft function is weakly differentiable and produces a high attenuation of the coefficients and, therefore, the reduction of the amplitude in the resulting signal.

In an attempt to get around these problems, many alternatives are being proposed. The main idea is to generate a high-derivative order thresholding function, which contributes to its use in optimization algorithms that look for the optimal parameters to be applied in the thresholding of each signal [34]. Therefore, the function becomes adaptable to the signal to be processed, improving the quality of the denoised signal.

When analyzing the threshold functions applied, whether in the area of PD, audio, ECG, or image processing, it is remarkable that those seek improvements by combining both, the preservation properties of the coefficients and magnitudes provided by the Hard function, as well as the differentiation and smoothness provided by the Soft function. In image processing, the smoothness property is interesting so that the resulting image shows more pleasant contours. In signal processing, such as audio, ECG, and PD pulse processing, it is important to achieve better preservation of signal magnitude (peak) and signal noise ratio.

For this reason, many authors have explored functions that correspond to an interpolation of the Soft and Hard alternatives. As an example, it is possible to mention functions such as: the Garrote described by Nasiri et al. in [19]; the Non Negative Garrote described in [12]; the Adaptive Shirinkage showed by Partha Ray in [35]; the Liu developed by Shan Liu in [16]; the Hui presented in [36]; Stein and Semi-Soft shown in [12]; and the functions described by Zhang et al. in [37, 38]. However, the majority of the functions cannot adapt to the different signals due to the fixed transition curve on the threshold value. In these functions, there is still a greater tendency to smooth the coefficients than to preserve them, not realizing that for PD signals it is appropriate for the function to be closer to the Hard them to the Soft threshold function, but still preserving some of the smoothness (differentiability) in the transition of the threshold value, which will allow an improvement in the EQM and CC.

Following this line of reasoning, we propose a new threshold function similar to the Hard but being differentiable for higher orders and being able to adjust to each signal. This proposal is based on the well-known logistic function, shown in **Figure 2**, widely used in artificial neural networks, demography, economics,

#### *A Wavelet Threshold Function for Treatment of Partial Discharge Measurements DOI: http://dx.doi.org/10.5772/intechopen.94115*

probability, statistics, chemistry, etc. Eq. (1) represents the logistic function, where *H* is the maximum value of the curve, *α* controls the slope of this curve and *x* corresponds to the value of *x* at the midpoint of the sigmoid curve dictated by the numerator value (**Figure 3**). When the *x* value tends to þ∞ the curve approaches *H* and when it tends to �∞ it approaches to zero.

$$f(\mathbf{x}) = \frac{H}{\mathbf{1} + e^{-a(\mathbf{x} - \mathbf{x}\_0)}} \tag{1}$$

Eq. (1) enabled us to develop the threshold function for filtering signals in such a way that it circumvents the problems previously described. As the objective in the thresholding process is to preserve the coefficients above the threshold value, it is easy to see that the maximum value *H* will be the decomposed wavelet coefficient *wj*,*<sup>k</sup>* (which corresponds to the variable *x*). Thus, the function will maintain symmetry when we vary the inclination constant *c* of the curve. Finally, it is necessary to move the function along the abscissa axis so that the graph shown in **Figure 2** leaves the ordinate axis and stays over the threshold value, this is done by subtracting the variable *x* from the value where we want to move the function (*x*0), i.e., the value of the coefficients *wj*,*<sup>k</sup>* must be subtracted from the threshold value *λ*. By making these adaptations, we obtain the following threshold function:

$$\eta\_f(o\_{j,k}, \lambda, c) = \frac{o\_{j,k}}{\mathbf{1} + e^{-c\left(o\_{j,k} - \lambda\right)}} = \frac{o\_{j,k}}{\mathbf{1} + e^{\xi\left(-o\_{j,k} + \lambda\right)}}\tag{2}$$

For a more efficient implementation, in which it is not necessary to worry about the fact that the coefficient *wj*,*<sup>k</sup>* is positive or negative, Eq. (2) can be rewritten using the signum (*sign*) function, which returns þ1 if the value is positive and �1 if the value is negative. Thus, we have:

$$\eta\_f(o\_{j,k}, \lambda, \mathfrak{c}) = \frac{o\_{j,k}}{\mathbf{1} + \mathbf{e}^{\mathfrak{c}\left(\left(-\operatorname{sign}\left(o\_{j,k}\right) \times a\_{j,k}\right) + \lambda\right)}}\tag{3}$$

With high *c* values, the curve inclination on the threshold point is such that it approaches the Hard function, but with a smoother (differentiable) transition. For low *c* values, the inclination of the function will act with less intensity on the coefficients below the threshold value and with greater intensity on the coefficients above this value, i.e., a large part of noisy coefficients may pass and there will be information losses on those coefficients that represent the signal of interest, in our case the PD pulse. With the appropriate choice of the *c* value for each processed signal, it is possible to obtain a significant improvement in the result of the PD wavelet denoising in relation to the Hard and the Soft functions.

Considering a threshold value *λ* ¼ 1, **Figure 4** shows the behavior of the proposed threshold function (called the Fleming function) for different *c* values (3, 10, 20, 30, 50, 80, 100 e 200). With very low values of *c* there is the possibility of

**Figure 3.** *Logistic function.*

**Figure 4.** *Behavior of the Fleming threshold function to λ* ¼ 1*.*

*<sup>η</sup> <sup>f</sup> <sup>ω</sup><sup>j</sup>*,*<sup>k</sup>*, *<sup>λ</sup>*,*<sup>c</sup>* � � <sup>¼</sup>

*DOI: http://dx.doi.org/10.5772/intechopen.94115*

*Behavior of the modified Fleming threshold function to λ* ¼ 0, 05*.*

**Figure 6.**

**317**

**3.1 Relevant wavelet coefficient identification**

or PD components (high kurtosis ≥3).

*ω<sup>j</sup>*,*<sup>k</sup>*

*ω<sup>j</sup>*,*<sup>k</sup>*

Therefore, in Eq. (4) we have a function capable of adapting to different types of wavelet coefficients, varying between the Soft and the Hard threshold functions according to the *c* inclination value defined. Thus, there is a need to define how

From the idea of identifying the most important coefficients to form the PD signal, used for the SNRBWS method, we were able to perform a variant on the threshold function. In this case, we chose to use kurtosis (*Ku*) as a statistical measure of the probability distribution's flatness [39] of the *ω<sup>j</sup>*,*<sup>k</sup>* wavelet coefficient, because the tapered this curve, the farther from the Normal probability distribution (Gaussian), which is characteristic of the white noise presence. Therefore, kurtosis will serve as an indicator to know if we have noisy coefficients (kurtosis close to 3)

**Figure 7** shows the detail coefficients at level *j* ¼ 1 and the detail coefficients at level *j* ¼ 6 with their respective histograms. Notice that in **Figure 7(a)**, formed almost exclusively by noise components, the histogram is very close to the Normal probability distribution, a fact that is confirmed by the kurtosis value equal to 2.9687; in the **Figure 7(b)** the coefficients have significant information about the PD pulse and the histogram is more tapered (leptokurtic), moving away from the

Then, to fulfill the task of identifying the most relevant coefficients to form the PD signal, it is enough to assume the following the condition regarding the kurtosis value: if the kurtosis of the coefficient is greater than 4, it must be considered important and the threshold function will make use of a lower *c* inclination

thus, when *λ*< 0, 5 the lower the *λ* threshold value, the greater the rigor in discarding the coefficients (closer to the Hard function), with a significant

*<sup>c</sup>*ð Þ ð Þ �*sign*ð Þ *<sup>ω</sup>j*,*<sup>k</sup>* �*ωj*,*<sup>k</sup>* <sup>þ</sup>*<sup>λ</sup> if <sup>λ</sup>*<sup>≥</sup> 0, 5

(4)

h i� *<sup>c</sup>=<sup>λ</sup>* ð Þ ð Þ *sign*ð Þ *<sup>ω</sup>j*,*<sup>k</sup>* �*ωj*,*<sup>k</sup>* <sup>þ</sup>*<sup>λ</sup> if <sup>λ</sup>*<sup>&</sup>lt; 0, 5

1 þ *e*

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements*

(and which) the inclination value should be applied to the coefficients.

Normal distribution, as indicated by the kurtosis value of 9.9612.

1 þ *e*

8 >><

>>:

improvement in the function's behavior, as shown in **Figure 6**.

**Figure 5.** *Behavior of the Fleming threshold function to λ* ¼ 0, 05*.*

passing a large number of noisy coefficients, so it is indicated that the value of the constant be greater than or equal to 5.

In a PD evaluation, most measurements provide signals with amplitude around mV. Thus, if the WT technique is applied to filter the signals, its decomposed coefficients will also be in the mV range and by using a threshold rule (in our case *scaledep*) the threshold value *λ* will be small and usually smaller than 1, mainly for coefficients that contain more noise than the PD components. When we evaluate the threshold function for a small threshold value (e.g., *λ* ¼ 0, 05), the accuracy with which the coefficients are attenuated becomes lower, as illustrated in **Figure 5**. Note that even for *c* ¼ 200, most of the noisy coefficients can pass, different than what was seen for the threshold value *λ* ¼ 1.

One solution to overcome this problem was to adapt Eq. (3) according to the threshold value when it is considered small (understand as small as *λ*< 0, 5), by simply changing the *c* constant that controls the inclination proportionally to the *λ* threshold values. With this, we can rewrite Eq. (3) as follows:

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements DOI: http://dx.doi.org/10.5772/intechopen.94115*

**Figure 6.** *Behavior of the modified Fleming threshold function to λ* ¼ 0, 05*.*

$$\eta\_f(o\_{j,k}, \lambda, c) = \begin{cases} \frac{o\_{j,k}}{\mathbf{1} + \mathbf{c}^\epsilon \big( \left( -\operatorname{sign} \left( {o\_{j,k}} \right) \times {}\_{j,k} \big) + \lambda \big)}} & \text{if } \lambda \ge \mathbf{0}, \mathbf{5} \\\\ \frac{o\_{j,k}}{\mathbf{1} + \mathbf{c}^{\left\langle \iota/\lambda \big\rangle} \big( \left( {-\operatorname{sign} \left( {o\_{j,k}} \right) \times {}\_{j,k} \big) + \lambda \big)} + \lambda \big)} & \text{if } \lambda < \mathbf{0}, \mathbf{5} \end{cases} \end{cases} \tag{4}$$

thus, when *λ*< 0, 5 the lower the *λ* threshold value, the greater the rigor in discarding the coefficients (closer to the Hard function), with a significant improvement in the function's behavior, as shown in **Figure 6**.

Therefore, in Eq. (4) we have a function capable of adapting to different types of wavelet coefficients, varying between the Soft and the Hard threshold functions according to the *c* inclination value defined. Thus, there is a need to define how (and which) the inclination value should be applied to the coefficients.

#### **3.1 Relevant wavelet coefficient identification**

From the idea of identifying the most important coefficients to form the PD signal, used for the SNRBWS method, we were able to perform a variant on the threshold function. In this case, we chose to use kurtosis (*Ku*) as a statistical measure of the probability distribution's flatness [39] of the *ω<sup>j</sup>*,*<sup>k</sup>* wavelet coefficient, because the tapered this curve, the farther from the Normal probability distribution (Gaussian), which is characteristic of the white noise presence. Therefore, kurtosis will serve as an indicator to know if we have noisy coefficients (kurtosis close to 3) or PD components (high kurtosis ≥3).

**Figure 7** shows the detail coefficients at level *j* ¼ 1 and the detail coefficients at level *j* ¼ 6 with their respective histograms. Notice that in **Figure 7(a)**, formed almost exclusively by noise components, the histogram is very close to the Normal probability distribution, a fact that is confirmed by the kurtosis value equal to 2.9687; in the **Figure 7(b)** the coefficients have significant information about the PD pulse and the histogram is more tapered (leptokurtic), moving away from the Normal distribution, as indicated by the kurtosis value of 9.9612.

Then, to fulfill the task of identifying the most relevant coefficients to form the PD signal, it is enough to assume the following the condition regarding the kurtosis value: if the kurtosis of the coefficient is greater than 4, it must be considered important and the threshold function will make use of a lower *c* inclination

**4. Fleming threshold function**

*DOI: http://dx.doi.org/10.5772/intechopen.94115*

thresholding functions mentioned in the Section 4.

compare the results for different values of this constant.

criteria were described in [42].

**Table 1.**

**319**

**4.1 Investigating the better inclination value** *c*

In order to perform the evaluation of the Fleming thresholding functions, we took 2064 signals and submitted to the wavelet denoising processes. Among these signals, we included real PD measurements from HV equipment and PD simulated with different levels of uniform white, Gaussian white and AM noise (created of the same way described in [40]). For each data, we compare the performance of our proposal against the classical Hard and Soft thresholding, along with 12 other

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements*

In addition to thresholding, the wavelet shrinkage process also requires the choice of the decomposition tree, the mother wavelet, the decomposition levels number and the threshold value ð Þ*λ* estimation method. As our goal is to evaluate only the thresholding functions performance, we change only these and keep fixed the other wavelet parameters necessary to the signal filtering. We chose to use the FWT structure, due to the ease of its implementation and because it is widely applied in the treatment of PD signals. We use the SNRBWS method to select the mother wavelet and the NWDLS method to find the decomposition levels number. In the threshold value estimative, we chose the *scaledep* method [3, 40, 41].

Since the Fleming function depends on an *c* inclination constant, which controls

The comparisons were done using statistical parameters as Absolute Mean Error (AME), Mean Square Error (*MSE*), Root Mean Square Error (*RMSE*), Correlation Coefficient (*CC*), Normalized Correlation Coefficient (*NCC*), Energy Difference (*EnD*), Signal to Noise Ratio (*SNR*), Signal to Noise Ratio Difference (*DSNR*), Noise Level Reduction (NLR), kurtosis difference (*Δk*); and local similarity criteria that involve maximum Magnitude Error (*MEmax*), minimum Magnitude Error (*MEmin*), maximum Peak Time Variation (*PTVmax*), minimum Peak Time Variation (*PTVmin*) and Rise Time Variation (*RTV*). Some of these parameters are used to form a fitness function (*JApt*), composed by global similarity criteria (*csg*) and local similarity criteria (*csl*), that can determine the best filtering result. All these

In a first analysis, it was investigated, through the *JApt* fitness criterion, what is

thresholding. **Table 1** evinced that *c* ¼ 5 produces the highest amount of best results per threshold function. In **Table 2** both methods produce best results with a lower constant, in case *c* ¼ 10 to Fleming and *c* ¼ 5 to Fleming 2. In this way, it is possible

Fleming **47,09** 25,00 8,96 4,12 3,63 1,50 1,07 1,16 7,46 Fleming 2 **65,16** 11,82 5,52 2,37 1,79 1,55 0,78 1,16 9,84

*Best results percentage (by JApt) comparison between Fleming functions to various inclination constants.*

**Inclination constant** *c* **5 10 30 50 100 200 300 500 1000**

the best *c* inclination value to be used in each alternative of the Fleming

to recommend not to use inclination values higher than 10.

**Function Best** *JApt* **results percentage**

how the decomposed wavelet coefficients are eliminated or attenuated, we also

**Figure 7.** *Wavelet coefficients and histograms of a real PD pulse: (a) first detail coefficient (kurtosis = 2.9687); (b) sixth detail coefficient (kurtosis = 9.9612).*

constant, allowing the passage of more coefficients, otherwise it will be considered as noisy coefficients and a much higher inclination constant must be assigned (in case *<sup>c</sup>* <sup>¼</sup> <sup>10</sup>20), eliminating a greater amount of noise, which approximates our function of the Hard. In equational terms, we have the Eq. (5):

$$\eta\_{f\_{1}}(o\_{\boldsymbol{j},k},\boldsymbol{\lambda},\boldsymbol{c})=\begin{cases} \frac{o\_{\boldsymbol{j},k}}{\mathbf{1}+\boldsymbol{e}^{10^{\lambda}\big(\left(-\operatorname{sign}\left(a\_{\boldsymbol{j},k}\right)\times a\_{\boldsymbol{j},k}\right)+\boldsymbol{\lambda}\big)}} & \text{if } K\_{u}\big(o\_{\boldsymbol{j},k}\big)\leq 4\\\begin{cases} \frac{o\_{\boldsymbol{j},k}}{\mathbf{1}+\boldsymbol{e}^{\ell\big(\left(-\operatorname{sign}\left(a\_{\boldsymbol{j},k}\right)\times a\_{\boldsymbol{j},k}\right)+\boldsymbol{\lambda}\big)}} & \text{if } \boldsymbol{\lambda}\geq \boldsymbol{0},\boldsymbol{5} \\\\ \frac{o\_{\boldsymbol{j},k}}{\mathbf{1}+\boldsymbol{e}^{\left(c/\boldsymbol{\lambda}\right)\big(\left(-\operatorname{sign}\left(a\_{\boldsymbol{j},k}\right)\times a\_{\boldsymbol{j},k}\big)+\boldsymbol{\lambda}\big)}} & \text{if } \boldsymbol{\lambda}<\boldsymbol{0},\boldsymbol{5} \end{cases} & \text{if } K\_{u}\succeq 4 \end{cases} \tag{5}$$

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements DOI: http://dx.doi.org/10.5772/intechopen.94115*

#### **4. Fleming threshold function**

In order to perform the evaluation of the Fleming thresholding functions, we took 2064 signals and submitted to the wavelet denoising processes. Among these signals, we included real PD measurements from HV equipment and PD simulated with different levels of uniform white, Gaussian white and AM noise (created of the same way described in [40]). For each data, we compare the performance of our proposal against the classical Hard and Soft thresholding, along with 12 other thresholding functions mentioned in the Section 4.

In addition to thresholding, the wavelet shrinkage process also requires the choice of the decomposition tree, the mother wavelet, the decomposition levels number and the threshold value ð Þ*λ* estimation method. As our goal is to evaluate only the thresholding functions performance, we change only these and keep fixed the other wavelet parameters necessary to the signal filtering. We chose to use the FWT structure, due to the ease of its implementation and because it is widely applied in the treatment of PD signals. We use the SNRBWS method to select the mother wavelet and the NWDLS method to find the decomposition levels number. In the threshold value estimative, we chose the *scaledep* method [3, 40, 41].

Since the Fleming function depends on an *c* inclination constant, which controls how the decomposed wavelet coefficients are eliminated or attenuated, we also compare the results for different values of this constant.

The comparisons were done using statistical parameters as Absolute Mean Error (AME), Mean Square Error (*MSE*), Root Mean Square Error (*RMSE*), Correlation Coefficient (*CC*), Normalized Correlation Coefficient (*NCC*), Energy Difference (*EnD*), Signal to Noise Ratio (*SNR*), Signal to Noise Ratio Difference (*DSNR*), Noise Level Reduction (NLR), kurtosis difference (*Δk*); and local similarity criteria that involve maximum Magnitude Error (*MEmax*), minimum Magnitude Error (*MEmin*), maximum Peak Time Variation (*PTVmax*), minimum Peak Time Variation (*PTVmin*) and Rise Time Variation (*RTV*). Some of these parameters are used to form a fitness function (*JApt*), composed by global similarity criteria (*csg*) and local similarity criteria (*csl*), that can determine the best filtering result. All these criteria were described in [42].

#### **4.1 Investigating the better inclination value** *c*

In a first analysis, it was investigated, through the *JApt* fitness criterion, what is the best *c* inclination value to be used in each alternative of the Fleming thresholding. **Table 1** evinced that *c* ¼ 5 produces the highest amount of best results per threshold function. In **Table 2** both methods produce best results with a lower constant, in case *c* ¼ 10 to Fleming and *c* ¼ 5 to Fleming 2. In this way, it is possible to recommend not to use inclination values higher than 10.


**Table 1.**

*Best results percentage (by JApt) comparison between Fleming functions to various inclination constants.*


**Table 2.**

*Mean value results (by JApt) comparison between Fleming functions to various inclination constants.*

#### **4.2 Comparison between Fleming, Hard and Soft threshold functions**

The main objective of building a dedicated threshold function is to make it able to produce results superior to those of conventional functions. As seen, the most applied functions in wavelet coefficient filtering are Hard and Soft, not only for PD signals, but also for image processing, audio signals, etc.

First, we show in **Figure 8** the results of the comparison between the first proposed alternative using *c* ¼ 5 against Hard and Soft functions. According to the *JApt*, we find that the proposed function achieves a higher percentage of better results than the Hard and Soft. As expected, due to its simplicity, the Soft thresholding is the fastest in runtime.

In **Figure 10** is shown a signal consisting of 3 simulated PD pulses wrapped in white noise and in AM noise, which was created as performed in [3]. In addition, note the filtering results for the Hard, Soft, Fleming and Fleming 2 thresholding. The Soft function tends to considerably attenuate the pulses peak amplitudes; the Hard function shows greater preservation of these amplitudes; the Fleming function allows the passage of a little more noise with negligible amplitudes, but it achieves better preservation of the amplitudes than the Hard and Soft, while the Fleming 2 function is able to solve Fleming's problem by identifying the coefficients of greater importance. In this way, the Fleming 2 method presents better amplitudes preservation than the other functions and still manages to eliminate the low amplitude noise seen with the use of the Fleming. The filtering improvement is also indicated by the fitness function, with higher value (JApt = 16.4607) for the filtering result

*Comparison of the better denoising results obtained between the Hard, Soft and Fleming 2 threshold functions.*

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements*

*DOI: http://dx.doi.org/10.5772/intechopen.94115*

*Comparison of the better denoising results obtained between the Hard, Soft, Fleming and Fleming 2 threshold*

using the Fleming 2 thresholding.

**Figure 10.**

**321**

**Figure 9.**

*functions to a simulated PD signal.*

We then compare in **Figure 9** the second alternative proposed with the Hard and Soft functions. Note that there is a significant improvement in the number of better results, achieving superior performance in the *EMA* and *Δk* criteria, which did not occur with the first alternative of our function.

Therefore, is evidenced by the superiority of the proposed alternatives in relation to the amount of better results obtained compared to the usual Hard and Soft methodologies. The only drawback is that our second proposal needs a little more time to be processed, but it is a relatively low price to be paid to achieve better results in reducing noisy components of PD signals. Also, note that, compared with the Soft function, the Hard thresholding tends to provide a better preservation of the PD pulses amplitudes and of the SNR, which confirms the statements made in the literature [3, 22].

**Figure 8.** *Comparison of the better denoising results obtained between the Hard, Soft and Fleming threshold functions.*

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements DOI: http://dx.doi.org/10.5772/intechopen.94115*

**Figure 9.** *Comparison of the better denoising results obtained between the Hard, Soft and Fleming 2 threshold functions.*

In **Figure 10** is shown a signal consisting of 3 simulated PD pulses wrapped in white noise and in AM noise, which was created as performed in [3]. In addition, note the filtering results for the Hard, Soft, Fleming and Fleming 2 thresholding. The Soft function tends to considerably attenuate the pulses peak amplitudes; the Hard function shows greater preservation of these amplitudes; the Fleming function allows the passage of a little more noise with negligible amplitudes, but it achieves better preservation of the amplitudes than the Hard and Soft, while the Fleming 2 function is able to solve Fleming's problem by identifying the coefficients of greater importance. In this way, the Fleming 2 method presents better amplitudes preservation than the other functions and still manages to eliminate the low amplitude noise seen with the use of the Fleming. The filtering improvement is also indicated by the fitness function, with higher value (JApt = 16.4607) for the filtering result using the Fleming 2 thresholding.

#### **Figure 10.**

*Comparison of the better denoising results obtained between the Hard, Soft, Fleming and Fleming 2 threshold functions to a simulated PD signal.*

#### **4.3 Comparison between all threshold functions**

With the results described in the previous subsection, we have a quantitative idea of the used method's capacity, but only with the average results is possible to have a real sense of the quality of each one. Taking advantage of the opportunity, we implement the various wavelet thresholding methods (mentioned in Section 4), including the: Adapt Shrink, Garrote, Hui, Liu, Non Negative Garrote (NNG), Semi Soft (SS), Stein, Zhang 1 (Z1), Zhang 2 (Z2), Zhang 3 (Z3), Zhang 4 (Z4), and Zhang 5 (Z5). The required variables for each of these alternatives were designated according to the specifications provided by the respective authors in the works that describe them.

**Parameter Threshold function**

*DOI: http://dx.doi.org/10.5772/intechopen.94115*

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements*

*Percentage of best results by evaluation parameters for all threshold functions.*

**Parameter Threshold function Hard Soft Fleming Fleming**

**2**

*AME* 0,021 0,020 0,021 0,021 0,021 0,023 0,023 0,019 *MSE* 0,020 0,017 0,020 0,020 0,017 0,021 0,025 0,016 *RMSE* 0,029 0,031 0,029 0,029 0,030 0,030 0,038 0,028 *CC* 0,770 0,736 0,790 0,787 0,761 0,773 0,713 0,775 *NCC* 0,769 0,733 0,789 0,786 0,759 0,773 0,706 0,773 *EnD* 0,302 0,525 0,267 0,283 0,472 0,364 0,571 0,294 *SNR* 6,47 1,22 6,82 6,83 2,57 6,51 0,00 5,40 *DSNR* 3,05 2,20 3,40 3,41 0,85 3,09 3,42 1,98 *NLR* 37,33 36,57 37,55 37,52 36,86 38,26 36,40 37,16 *MEmax* 14,13 39,75 12,97 14,77 35,08 13,78 43,78 20,53 *MEmin* 16,11 35,57 13,20 16,24 30,09 14,03 40,20 17,37 *PTVmax* 15,71 11,44 11,66 16,28 13,97 13,29 12,45 9,08 *PTVmin* 23,06 20,14 19,93 23,93 23,45 22,07 23,58 16,49

**Adapt shrink** **Garrote Hui Liu**

**Table 4.**

**323**

**NNG SS Stein Z1 Z2 Z3 Z4 Z5**

*AME* 0,97 1,02 6,35 1,31 4,36 4,12 7,03 3,05 *MSE* 2,28 1,11 6,83 1,11 5,14 3,05 2,37 3,29 *RMSE* 1,70 1,11 7,41 1,11 5,14 2,96 2,37 3,39 *CC* 1,70 0,39 5,14 1,60 1,55 2,03 0,00 4,46 *NCC* 1,31 0,39 5,28 1,55 1,07 1,89 0,00 4,75 *EnD* 2,03 0,48 5,43 2,86 3,59 4,12 1,55 5,57 *SNR* 1,31 0,78 5,33 0,87 0,05 0,97 0,00 17,34 *DSNR* 1,21 0,78 5,43 0,87 0,05 0,97 0,00 17,34 *NLR* 0,00 0,00 0,00 26,26 0,00 0,00 0,00 **73,74** *MEmax* 1,70 0,15 5,81 9,06 1,36 4,36 0,24 13,28 *MEmin* 1,11 0,05 11,05 5,18 2,71 5,28 0,15 6,40 *PTVmax* 0,00 1,16 1,79 6,15 5,33 20,16 2,52 **41,38** *PTVmin* 0,00 1,16 1,26 6,06 4,99 18,90 1,31 **42,05** DM 0,00 1,79 0,68 7,36 2,42 12,74 37,26 **20,74** *RTV* 0,00 2,66 3,34 4,65 4,65 14,24 18,94 **21,08** *csg* 4,02 2,57 7,95 0,24 1,07 3,59 0,00 4,26 *csl* 1,55 1,02 8,62 4,55 1,41 4,89 0,78 6,73 *Δk* 3,49 3,44 6,59 2,81 2,96 6,06 0,29 5,18 *JApt* 4,41 2,33 8,19 0,44 0,87 4,36 0,00 4,02 *texec* 0,24 2,37 9,93 18,02 0,00 3,05 0,58 **33,38**

Similarly to what was done in **Table 1**, we made a percentage evaluation of the amount of best results considering all threshold functions and the proposed Fleming functions (compared for a constant c = 5). From **Tables 3** and **4**, the bold values evidence that the Fleming threshold had a superior performance when compared to the other alternatives. In terms of fitness, the one with the highest amount of better filtering results was Fleming. The Stein function outperforms the others in execution time. The Soft function ends up losing space in practically all the evaluated criteria, confirming that it is not suitable to treat PD signals, due to the high attenuation generated in the wavelet coefficients processing.


**Table 3.** *Percentage of best results by evaluation parameters for all threshold functions.*


*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements DOI: http://dx.doi.org/10.5772/intechopen.94115*

#### **Table 4.**

*Percentage of best results by evaluation parameters for all threshold functions.*



#### **Table 5.**

*Average results by evaluation parameters for all threshold functions.*


#### **Table 6.**

*Average results by evaluation parameters for all threshold functions.*

Also was evaluated the average results of the evaluation parameters, according to the **Tables 5** and **6**. Note that the fitness JApt = 3,20 of the Fleming function (using c = 10) is the highest among all the others, being followed by the Garrote and the Hard functions. Thus, the proposed alternatives achieve the objective of

**Figure 11.**

**325**

*PD signal from a hydro generator.*

*Comparison of the better denoising results obtained between the all evaluated threshold functions to a measured*

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements*

*DOI: http://dx.doi.org/10.5772/intechopen.94115*

#### *A Wavelet Threshold Function for Treatment of Partial Discharge Measurements DOI: http://dx.doi.org/10.5772/intechopen.94115*

#### **Figure 11.**

*Comparison of the better denoising results obtained between the all evaluated threshold functions to a measured PD signal from a hydro generator.*

overcoming other methods, also providing a better qualitative result in the treatment of PD pulses.

**Figure 11** exemplifies the wavelet shrinkage process using each of the thresholding functions discussed above for a PD signal measured from a hydro generator. In this case, note that the functions we have created are superior in preserving the amplitudes of the signals and eliminating the present noise, especially the filtering using the Fleming 2 function, which obtained the highest level of JApt compared to the other functions, followed by Fleming, Garrote and Hard functions. The Soft and the other thresholding alternatives end up causing deformations of the pulses waveforms and greater attenuation of these peak amplitudes. The Zhang 1 and Zhang 5 functions allow most of the noise to pass through the denoising process and the Zhang 4 function ends up eliminating the PD signal that we are interested in obtaining.

#### **5. Conclusions**

Was presented a new threshold function called Fleming, which combines the quality of a strongly differentiable function and a more flexible alternative, enabling its optimization to provide better results in the PD signals treatment, in order to preserve its important characteristics for the diagnosis of the HV equipment subjected to the partial discharge analysis. The proposal inspired by the wellknown logistic function [29], which depends on a parameter that controls the inclination of the curve in the threshold value (calculated a priori). Also was created a variant of this same function, using a simple idea, but little investigated in the literature: identifying the decomposed coefficients with the greatest contribution in the desired signal recovering [3].

With the results described in Section 5, in which hundreds of signals (measured and simulated) were evaluated, the ability of the Fleming function and its Fleming 2 variant to overcome the most common functions such as Hard and Soft, as well as twelve other alternatives presented in some publications [15, 19, 22, 35–38]. The Fleming function can be applied with different inclination values, but for PD signals, the ideal is that these values are limited between 5 and 10 to provide the best results.

The Fleming 2 alternative showed the highest percentage of the best results and the Fleming alternative showed the highest average value in terms of amplitude. Thus, if the goal is to achieve a higher number of better results, the indicated is to threshold the wavelet coefficients using the Fleming 2 function, but if the idea is to achieve better average results, consider using the Fleming function. As for the average processing time, these functions are relatively fast when compared to the other evaluated functions, not falling far behind the classic ones Hard and Soft.

**Author details**

**327**

and Antonio C.S. Lima<sup>1</sup>

Caio F.F.C. Cunha1,2\*, Mariane R. Petraglia1

provided the original work is properly cited.

\*Address all correspondence to: caioflemin@gmail.com

1 COPPE/UFRJ, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements*

*DOI: http://dx.doi.org/10.5772/intechopen.94115*

2 CEPEL/ELETROBRAS, Electrical Energy Research Center, Rio de Janeiro, Brazil

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

, André T. Carvalho<sup>2</sup>

The application of the developed thresholding functions is extensible to other types of signals, such as acoustic emissions, electrocardiogram signals, image processing, among others. However, in each case it would be necessary to investigate the appropriate values of the inclination parameter *c*.

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements DOI: http://dx.doi.org/10.5772/intechopen.94115*

## **Author details**

Caio F.F.C. Cunha1,2\*, Mariane R. Petraglia1 , André T. Carvalho<sup>2</sup> and Antonio C.S. Lima<sup>1</sup>

1 COPPE/UFRJ, Federal University of Rio de Janeiro, Rio de Janeiro, Brazil

2 CEPEL/ELETROBRAS, Electrical Energy Research Center, Rio de Janeiro, Brazil

\*Address all correspondence to: caioflemin@gmail.com

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## **References**

[1] Carvalho A T, Lima A C, Cunha C F, Mariane R P. Identification of partial discharges immersed in noise in large hydro-generators based on improved wavelet selection methods. Measurement 2015; 75: 122–133. DOI: 10.1016/j.measurement.2015.07.050

[2] Cunha C F, Carvalho A T, Petraglia M R, Lima A C. A new wavelet selection method for partial discharge denoising. Electric Power Systems Research 2015; 125: 184–195. DOI: 10.1016/j. epsr.2015.04.005.

[3] Gaouda A, El-Hag A, Abdel-Galil T, Salama MM, Bartinikas R. On-line detection and measurement of partial discharge signals in a noisy environment. IEEE Transactions on Dielectrics and Electrical Insulation. 2008;**15**:1162-1173. DOI: 10.1109/ TDEI.2008. 4591239

[4] Cunha CF, Carvalho AT, Petraglia MR, Lima AC. An improved scale dependent wavelet selection for data denoising of partial discharge measurement. IEEE International Conference on Solid Dielectrics (ICSD). 30 June-4 July 2013. Bologna. IEEE. 2013:100-104

[5] Hongxia L. Xuefeng Z. A method of second wavelet transform automated threshold for partial discharge signal extraction. In: Second International Conference on Digital Manufacturing and Automation (ICDMA). 5–7 August 2011. Zhangjiajie. IEEE. 2011. pp. 42-45. DOI: 10.1109/ICDMA.2011.19.

[6] Shin I, Soraghan J, Siew W. Detection of PD utilizing digital signal processing methods. Part 3: open-loop noise reduction. Electrical Insulation Magazine. IEEE. 2001;**17**:6-13. DOI: 10.1109/57.901611

[7] Li J, Sun C, Yang J. Adaptive denoising for PD online monitoring based on wavelet transform. In: Proceedings of the SoutheastCon; 31 March 2006; Memphis: IEEE. 2006. pp. 71-74

Electrical Insulation. 2005;**12**:586-594. DOI: 10.1109/TDEI.2005.1453464

*DOI: http://dx.doi.org/10.5772/intechopen.94115*

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements*

[20] He H, Wang Z, Tan Y. Noise reduction of ECG signals through genetic optimized wavelet threshold filtering. In: IEEE International Conference on Computational

Intelligence and Virtual Environments

June 2015. Shenzhen. IEEE. 2015. pp. 1-6

[21] Tulsani H. Gupta R. Wavelet based optimized polynomial threshold function for ECG signal denoising. In: International Conference on Computing for Sustainable Global Development (INDIACom). 11-13 March 2015; New

for Measurement Systems and Applications (CIVEMSA). 12–14

Delhi. IEEE;**2015**:1563-1566

1012-1025. DOI: 10.1016/j. neucom.2008.04.016

Sciences. 2012;**1**:60-65

International Conference on Computational Intelligence and Computing Research (ICCIC). 28–29 December 2010. Coimbatore. IEEE.

2011. pp. 1-4

IEEE. 2011. pp. 721-724

DOI: 10.1109/83.862633

[22] Nasri M, Nezamabadi-pour H. Image denoising in the wavelet domain using a new adaptive thresholding function. Neurocomputing. 2009;**72**:

[23] Jangra S. Kumar S. A new threshold function for image denoising based on wavelet transform. International Journal of Engineering and Mathematical

[24] Mohideen S, Perumal S, Krishnan N. Selvakumar R. A novel approach for image denoising using dynamic tracking with new threshold technique. In: IEEE

[25] Norouzzadeh Y, Rashidi M. Image denoising in wavelet domain using a new thresholding function. In: International Conference on Information Science and Technology. 26–28 March 2011. Nanjing.

[26] Chang SG, Yu B, Vetterli M. Adaptive wavelet thresholding for image denoising and compression. IEEE Transactions on Image Processing. 2000;**9**:1532-1546.

[14] Vishwakarma D K, Kapoor R, Dhiman A, Goyal A, Jamil D. De-noising of audio signal using heavy tailed distribution and comparison of wavelets

and thresholding techniques. In:

Delhi. IEEE; 2015. p. 755–760.

[15] Tulsani H, Gupta R. 1-D signal denoising using wavelets based optimization of polynomial threshold function. In: 3rd International Conference on Reliability, Infocom Technologies and Optimization (ICRITO) (Trends and Future

Directions). 8–10 October 2014; Noida.

[16] Zhang WQ. Song G X. A translationinvariant wavelet de-noising method based on a new thresholding function. In: International Conference on Machine Learning and Cybernetics. 5 November 2003; Xi'an. IEEE. 2003. pp. 2341-2345

[17] Dong X, Yue Y, Qin X, Wang X, Tao Z. Signal denoising based on improved wavelet packet thresholding function. In: 2010 International Conference on Computer,

Mechatronics, Control and Electronic Engineering. 24-26 August 2014; Changchun. IEEE; 2014. p. 382–385.

[18] Guo Q, Yu S, Leng M, Liu W. Signal denoising using simplified SURE based linear expansion of thresholds in the wavelet domain. In: Congress on Image and Signal Processing. 27–30 May 2008;

Sanya. IEEE. 2008. pp. 592-596

wavelet transform. Journal of

jics20102845

**329**

[19] Liu S, Li Y, Hu X, Liu L, Hao D. A novel thresholding method in removing noises of electrocardiogram based on

Information & Computational Science. 2013;**10**:5031-5041. DOI: 10.12733/

IEEE. 2014. pp. 1-5

International Conference on Computing for Sustainable Global Development (INDIACom); 11–13 March 2015; New

[8] Li J, Jiang T, Grzybowski S, et al. Scale dependent wavelet selection for de-noising of partial discharge detection. Dielectrics and Electrical Insulation. IEEE Transactions; 2010; 17, 1705–1714. DOU: 10.1109/ TDEI.2010.5658220.

[9] Evagorou D, Kyprianou A, Lewin P, et al. Feature extraction of partial discharge signals using the wavelet packet transform and classification with a probabilistic neural network. IET Science Measurement & Technology. 2010;**4**:177-192. DOI: 10.1049/ietsmt.2009.0023

[10] Zhaoheng D, Shanghe L, Lei W. Selection of the optimal wavelet bases for wavelet de-noising of partial discharge signal. In: 2nd International Conference on Signal Processing Systems (ICSPS); 5–7 July 2010; Dalian. IEEE. 2010. pp. 400-404

[11] Hao L, Lewin P, Hunter J, et al. Discrimination of multiple PD sources using wavelet decomposition and principal component analysis. IEEE Transactions on Dielectrics and Electrical Insulation. 2011; 18: 1702– 1711. DOI: 10.1109/TDEI.2011. 6032842.

[12] Xiaorong W, Zongjun G, Yong S, et al. Extraction of partial discharge pulse via wavelet shrinkage. In: Proceedings of the 6th International Conference on Properties and Applications of Dielectric Materials; 21– 26 June 2000; Xi'an. IEEE. 2002. pp. 685-688. DOI: 10.1109/ ICPADM.2000.876113.

[13] Zhou X, Zhou C, Kemp IJ. An improved methodology for application of wavelet transform to partial discharge measurement denoising. IEEE Transactions on Dielectrics and

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements DOI: http://dx.doi.org/10.5772/intechopen.94115*

Electrical Insulation. 2005;**12**:586-594. DOI: 10.1109/TDEI.2005.1453464

[14] Vishwakarma D K, Kapoor R, Dhiman A, Goyal A, Jamil D. De-noising of audio signal using heavy tailed distribution and comparison of wavelets and thresholding techniques. In: International Conference on Computing for Sustainable Global Development (INDIACom); 11–13 March 2015; New Delhi. IEEE; 2015. p. 755–760.

[15] Tulsani H, Gupta R. 1-D signal denoising using wavelets based optimization of polynomial threshold function. In: 3rd International Conference on Reliability, Infocom Technologies and Optimization (ICRITO) (Trends and Future Directions). 8–10 October 2014; Noida. IEEE. 2014. pp. 1-5

[16] Zhang WQ. Song G X. A translationinvariant wavelet de-noising method based on a new thresholding function. In: International Conference on Machine Learning and Cybernetics. 5 November 2003; Xi'an. IEEE. 2003. pp. 2341-2345

[17] Dong X, Yue Y, Qin X, Wang X, Tao Z. Signal denoising based on improved wavelet packet thresholding function. In: 2010 International Conference on Computer, Mechatronics, Control and Electronic Engineering. 24-26 August 2014; Changchun. IEEE; 2014. p. 382–385.

[18] Guo Q, Yu S, Leng M, Liu W. Signal denoising using simplified SURE based linear expansion of thresholds in the wavelet domain. In: Congress on Image and Signal Processing. 27–30 May 2008; Sanya. IEEE. 2008. pp. 592-596

[19] Liu S, Li Y, Hu X, Liu L, Hao D. A novel thresholding method in removing noises of electrocardiogram based on wavelet transform. Journal of Information & Computational Science. 2013;**10**:5031-5041. DOI: 10.12733/ jics20102845

[20] He H, Wang Z, Tan Y. Noise reduction of ECG signals through genetic optimized wavelet threshold filtering. In: IEEE International Conference on Computational Intelligence and Virtual Environments for Measurement Systems and Applications (CIVEMSA). 12–14 June 2015. Shenzhen. IEEE. 2015. pp. 1-6

[21] Tulsani H. Gupta R. Wavelet based optimized polynomial threshold function for ECG signal denoising. In: International Conference on Computing for Sustainable Global Development (INDIACom). 11-13 March 2015; New Delhi. IEEE;**2015**:1563-1566

[22] Nasri M, Nezamabadi-pour H. Image denoising in the wavelet domain using a new adaptive thresholding function. Neurocomputing. 2009;**72**: 1012-1025. DOI: 10.1016/j. neucom.2008.04.016

[23] Jangra S. Kumar S. A new threshold function for image denoising based on wavelet transform. International Journal of Engineering and Mathematical Sciences. 2012;**1**:60-65

[24] Mohideen S, Perumal S, Krishnan N. Selvakumar R. A novel approach for image denoising using dynamic tracking with new threshold technique. In: IEEE International Conference on Computational Intelligence and Computing Research (ICCIC). 28–29 December 2010. Coimbatore. IEEE. 2011. pp. 1-4

[25] Norouzzadeh Y, Rashidi M. Image denoising in wavelet domain using a new thresholding function. In: International Conference on Information Science and Technology. 26–28 March 2011. Nanjing. IEEE. 2011. pp. 721-724

[26] Chang SG, Yu B, Vetterli M. Adaptive wavelet thresholding for image denoising and compression. IEEE Transactions on Image Processing. 2000;**9**:1532-1546. DOI: 10.1109/83.862633

[27] Biswas M, Om H. An adaptive wavelet thresholding image denoising method. In: National Conference on Communications (NCC). 15-17 February 2013. New Delhi. IEEE; 2013. p. 1–5.

[28] Begum AS, Poornachandra S. Curvelet based image de-noising using beta-trim shrinkage for magnetic resonance images. 2014 International Conference on Science Engineering and Management Research (ICSEMR). 27-29 November 2014. Chennai. IEEE. 2014: 1-8

[29] Lumen. Build a logistic model from data [Internet]. Available from: https:// courses.lumenlearning.com/ivytechcollegealgebra/chapter/build-a-logisticmodel-from-data/.

[30] Ma X, Zhou C, Kemp I. Interpretation of wavelet analysis and its application in partial discharge detections. IEEE Transactions on Dielectrics and Electrical Insulation. 2002;**9**:446-457. DOI: 10.1109/ TDEI.2002.1007709.

[31] Donoho DL, Johnstone IM. Adapting to unknown smoothness via wavelet shrinkage. Journal of the American Statistical Association. 1995; **90**:1200-1224. DOI: 10.1080/ 01621459.1995.10476626

[32] Li J, Sun C, Yang J. Adaptive denoising for PD online monitoring based on wavelet transform. IEEE Southeast Conference. Memphis. 31 March – 2 April 2006. IEEE; 2006. p. 71–74.

[33] Satish L, Zaengl WS. Artificial neural networks for recognition of 3-d partial discharge patterns. IEEE Transactions on Dielectrics and Electrical Insulation. 1994;**1**:265-275. DOI: 10.1109/94.300259

[34] Li J, Cheng C, Jiang T, Grzybowski S. Wavelet de-noising of partial discharge signals based on

genetic adaptive threshold estimation. IEEE Transactions on Dielectrics and Electrical Insulation. 2012;**19**:543-549. DOI: 10.1109/TDEI.2012.6180248

1550016. DOI: 10.1142/ S0219691315500162.

[42] Cunha CF, Carvalho AT, Mariane RP, Hélio PA, Lima AC. Proposal of a novel fitness function for evaluation of wavelet shrinkage parameters on partial discharge

10.1049/iet-smt.2017.0201

**331**

denoising. IET Science, Measurement & Technology. 2018;**12**:283-289. DOI:

*DOI: http://dx.doi.org/10.5772/intechopen.94115*

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements*

[35] Ray P, Maitra AK, Basuray A. Extract Partial Discharge signal using wavelet for on-line measurement. In: International Conference on Communications and Signal Processing (ICCSP). 3–5 April 2013. Melmaruvathur. IEEE. 2013. pp. 888-892. DOI: 10.1109/ iccsp.2013.6577184.

[36] Huimin C, Ruimei Z, Yanli H. Improved threshold denoising method based on wavelet transform. Physics Procedia. 2012;**33**:1354-1359. DOI: 10.1016/j.phpro.2012. 05.222

[37] Zang H, Wang Z, Zheng Y. Analysis of signal de-noising method based on an improved wavelet thresholding. In: 9th International Conference on Electronic Measurement Instruments (ICEMI). 16– 19 August 2009. Beijing. IEEE: 2009; p. 1–987–1–990.

[38] Zhang X. Thresholding neural network for adaptive noise reduction. IEEE Transactions on Neural Networks. 2001;**12**:567-584. DOI: 10.1109/ 72.925559

[39] Celso G, Steven WM, Edward O, James AY. Final state sensitivity: An obstruction to predictability. Physics Letters A. 1983;**99**:415-418. DOI: 10.1016/0375-9601(83)90945-3

[40] Cunha CF. Aplicação da Transformada Wavelet na redução de ruídos em medições de descargas parciais [dissertation]. Rio de Janeiro: Federal University of Rio de Janeiro (UFRJ; 2013

[41] Omari T, Bereksi-Reguig F. An automatic wavelet denoising scheme for heart sounds. International Journal of Wavelets, Multiresolution and Information Processing. 2015; 13:

*A Wavelet Threshold Function for Treatment of Partial Discharge Measurements DOI: http://dx.doi.org/10.5772/intechopen.94115*

1550016. DOI: 10.1142/ S0219691315500162.

[42] Cunha CF, Carvalho AT, Mariane RP, Hélio PA, Lima AC. Proposal of a novel fitness function for evaluation of wavelet shrinkage parameters on partial discharge denoising. IET Science, Measurement & Technology. 2018;**12**:283-289. DOI: 10.1049/iet-smt.2017.0201

**Chapter 15**

**Abstract**

**1. Introduction**

**333**

Functions

*Paulo César Linhares da Silva*

Use of Daubechies Wavelets in the

This chapter aims to use Daubechies' wavelets as basis functions to generate analytical functions, thus being able to rewrite the Taylor series using these wavelets. This makes it possible to analyze functions with a high degree of complexity, in problems that require a high degree of precision in their solution. Wavelet analysis can be applied to practical problems that require a high degree of precision, for example, in the study and analysis of electromagnetic propagation in optical fibers, solutions of differential equations involving engineering problems, in the transmission of WiFi signals, in the treatment and analysis of biomedical

**Keywords:** wavelets, Daubechies, analytical functions, basis functions, Taylor series

Wavelets [1] were born from the need to generate functions, especially those that present singularities, high gradients, discontinuities both in the time domain and in the frequency domain. Wavelets enable the high-resolution analysis of functions with these characteristics. An example of a problem that occurs when generating functions with a Fourier base is the Gibbs phenomenon. Such a phenomenon occurs because there is no way to represent functions that present discontinuities, even adding more elements in the base that will generate the function. A characteristic of wavelets is that they do not produce such an effect. Wavelets are widely used in the solution of numerical problems in several areas

of knowledge such as image compression, Numerical Harmonic Analysis [2], financial analysis, oil detection, differential Equations [3, 4], biomedical signals, analysis of electromagnetic integral Equations [5], optical fibers [6], among others.

Many of these applications use the specific properties of wavelets, such as coefficients that are determined numerically, multi-resolution analysis to decompose a signal, integrals, and derivatives obtained numerically, energy concentrated

For the development of topics presented in this chapter, the reader must have as a prerequisite knowledge of functional analysis, linear algebra, measure theory

in its compact and base with orthogonal elements.

**2. Short introduction to wavelet theory**

Representation of Analytical

images, detection of oil sources through the study of seismic signals.

#### **Chapter 15**
