Singularity detection of 2 D signals using fractal dimension analysis of scale information

Fractal dimension (FD) analysis has been widely used in signal processing. The key issue in signal processing is the singularity detection. One of the main problems for FD analysis of signals is its susceptibility to measurement noise, likely obscuring the identification of singularities. To address this deficiency, a new physical quantity, named ‘the scale-window fractal dimension (SWFD)’, is proposed and a SWFD analysis method is formed to identify the singularities in the noisy 2D signal. With this method, the noisy 2D signal first is decomposed into sets of scale signals with the aid of 2D Gabor wavelet transforms; then SWFD estimates are calculated along every scale signals to form the FD surface. The singularities can be localized by the sudden changes in the spatial variation of the FD surface. As an application of the method, the identification of damage singularity for an experimental composite plate is performed with the mode shapes measured by a scanning laser vibrometer as the analyzed 2D signals. The results show that the SWFD analysis method has the prominent features of high accuracy of singularity localization and strong robustness to noise.


Introduction
Singularity, usually caring the most important information of signal, can be used to determine the abnormal phenomenon of the signal.Hitherto, many kinds of singularity detection techniques have been developed [1], among which the fractal dimension (FD) analysis based methods have been the most widely used ones [2][3][4][5][6][7][8].The FD is a ratio that provides a statistical index of complexity, which compares how the detail in a pattern changes with the scale at which it is measured.Many algorithms for estimating the FD have been proposed, such as algorithms by Pickover and Khorasani [2], Katz [3,4], Higuchi [5], Maragos and Sun [6], Petrosian [7].In time series analysis, Katz's fractal dimension (KFD) calculation is derived directly from a 1D waveform signal, eliminating the preprocessing step of creating a binary sequence and is widely used to quantify the singularity or complexity of a waveform.The KFD for a waveform is defined by [3]: where is the number of spatial sampling intervals of the waveform under investigation, is the maximum distance between the beginning point and the th point of the sampling sequence under investigation and is the sum of distances between successive points.One of the successful applications of KFD is to identify damage of beams for health monitoring by the singularity analysis of the beam's vibrational mode shapes.The typical studies can be seen in references [9][10][11][12][13][14].These studies show that the KFD is high sensitivity to the singularities of the mode shapes and can be used to determine the location of the damage.However, an intrinsic deficiency is its susceptibility to measurement noise, likely obscuring the features of damage; besides, few studies applied the KFD to 2D signals for singularity extraction.To address these problems, this study presents an alternative method, which called 'scale-window fractal dimension (SWFD) analysis', to study the subtle singularities from the noisy 2D signals.This method involves two main processes: 1) Multi-resolution analysis relying on a 2D wavelet is performed to decompose a noisy 2D signal into multiple frequency-band components.Owing to the difference in frequency content, the noise component and the effective signal component that carries singularity information commonly occupy distinct frequency bands; 2) Parallel SWFD analysis is performed on every frequency-band component to produce a group of FD surface spectrum, of which the FD surface of the effective signal component accommodates the feature of singularity by eliminating the noise component.The method is experimentally validated using a scanning laser vibrometer (SLV) to acquire mode shape signals of an actual composite plate.
One of the suitable wavelet for the above processing is the 2D Gabor wavelet which can be defined as [15]: where is the maximum frequency of the modulating sinusoidal plane wave and is the orientation of the major axis of the elliptical Gaussian.To obtain the scale 2D signals, a simplified form of Eq. ( 3) can be used here by making = 1, = = 1, = 0. : ( , , , , , ) can be transformed into a series of scale 2D signals : ( , , , , , ) through Gabor wavelet.The lower-convey higher-frequency components of the 2D signal, and the higher-embody lower-frequency components of the 2D signal.Usually, the noise, trend, and singularity components have higher, lower, and in-between frequencies, respectively.Hence, they can be distributed at different-scale spaces.

SWFD analysis
To facilitate the use of KFD for 1D signal, : ( , , , , , ), can be considered as a set of 1D lines approximately, such as: : , , , , , ≈ : , , , , , or : , , , , , ≈ : , , , , , To locate the singularity of : ( , , , , , ), a sliding window technique is adopted to produce FD trajectories for every signal lines of : ( , , , , ) and : ( , , , , ).Let : ( , ) be a general instance of signal lines of : ( , , , , , ).In particular, in the case of = and = , : ( , ) is a signal line of : ( , , , , ) in the direction; in the case of = and = , : ( , ) is a signal line of : ( , , , , ) in the direction.Let : , be a segment of : ( , ), covering sampling points with the midpoint .A retrofitted KFD defined on is given as: where and are the scale parameters corresponding to Eq. ( 1).The retrofitted KFD, which we call "SWFD", defines on the segment of the scale signals, not the whole original signal, quantifying the complexity of : , centered at .The new SWFD provides a way to observe the singularity characteristics from the scale signals.Furthermore, SWFD can calculate FD for each sampling points of the signal, so it has the potential to locate the singularities in space.
By regarding : , as a sliding window, we can move it to cover the signal line : ( , ) step by step, evaluating the SWFD in terms of Eq. ( 6) at each step.After the window slides across the signal line, a SWFD trajectory can be formed that depicts the line's profile of complexity.
Corresponding to Eq. ( 5), fusion of Π and Π by averaging is adopted to form an --SWFD surface: Π provides a full prospect of complexity of : ( , , , , , ) by calculating the estimates of SWFD for every sampling points.Singularity can cause a sharply increased SWFD peak of an --SWFD surface, overwhelming the neighboring elements.This SWFD peak creates energy concentration, resulting in low Shannon entropy for the surface.In contrast, noise or global trend components cannot induce energy concentration, usually giving rise to relatively greater Shannon entropy.This effect provides a Shannon entropy-based method [16] for selecting the proper --SWFD surface at a particular scale from the --SWFD surface spectrum of the scale 2D signals for singularity identification.

Experimental example
Damage detection in laminated composite plates is a research focus in the filed of structural health monitoring.Vibrational mode shape, a 2D signal, conveys spatial characteristics of the plate being inspected and has the potential to be used to locate damage.The damage detection of an experimental glass fiber reinforced polymer (GFRP) square plate is used here to display the effectiveness of the above SWFD analysis method by calculating the singularities of the mode shape.

Experimental setup
The GFRP plate has dimensions of length 400 mm, width 400 mm, and depth 3 mm.This specimen consists of four-ply unidirectional glass fiber mats oriented in the direction, with one out of four layers are removed from the specimen locally to simulate damage, as enclosed by the little square in Fig. 1.Vibration of the specimen is created using a circular 10 mm-diameter piezoelectric lead-zirconate-titanate (PZT) actuator, located at the geometrical center of the plate, to excite the plate harmonically and perpendicularly.When the specimen vibrates stably, a scanning laser vibrometer (SLV) (Polytec PSV-400) is employed to scan its non-heated surface, to register the isochronous transverse velocities at 374×374 sampling points evenly distributed across the surface.The array of isochronous transverse velocities constitutes mode shapes for the plate, as typically illustrated in Fig. 2 for the mode shapes at natural frequencies of 388.13 Hz and 1126.56Hz.The mode shapes were disturbed by experimental uncertainty and measurement noise due to the limitation of measurement accuracy.

Results
For the mode shape at natural frequency of 388.13 Hz, the proper --SWFD surface identified by the Shannon entropy-based method is Π shown in Fig. 3.A cluster of much higher-intensity patches can be clearly observed towards the upper right of the planform, clearly delineating the distribution of the damage.Similar results can be observed under the natural frequency of 1126.56Hz, as shown in Fig. 4 with Π as the proper --SWFD surface.These observations validate that the proposed method can localize damage singularity in actual composite plates with great accuracy.

Conclusions
New physical quantity, the SWFD, is proposed with the aim of detecting the singularity of the signal, while overcoming the weakness of susceptibility to noise of KFD.A SWFD analysis method is then formed based on combining the merits of wavelet transform for decomposing the signal into distinct scale signals to separate the noise and singularity information in different frequency bands and KFD in singularity characterization.The identification of damage for an experimental GFRP plate is used as an application of the proposed method, with the mode shapes measured by SLV as the analyzed 2D signals.The results show that the SWFD analysis method has an excellent ability to identify singularities in noisy 2D signals, and it's a practical technology for structure damage diagnosis in plates.

Fig. 5 ( 2 .Fig. 3 .Fig. 4 . 5 . 6 .
Fig.5(a) and 5(b) present the results of the planforms of --surfaces Π , obtained by directly applying the KFD to mode shape lines at natural frequencies of 388.13 Hz and 1126.56Hz, respectively, while exclusive of the scale concept.Moreover, planforms of the