A New Classification Method of Infrasound Events Using Hilbert-Huang Transform and Support Vector Machine

Infrasound is a type of low frequency signal that occurs in nature and results fromman-made events, typically ranging in frequency from 0.01Hz to 20Hz. In this paper, a classificationmethod based onHilbert-Huang transform (HHT) and support vectormachine (SVM) is proposed to discriminate between three different natural events. The frequency spectrum characteristics of infrasound signals produced by different events, such as volcanoes, are unique, which lays the foundation for infrasound signal classification. First, the HHT method was used to extract the feature vectors of several kinds of infrasound events from the Hilbert marginal spectrum. Then, the feature vectors were classified by the SVM method. Finally, the present of classification and identification accuracy are given.The simulation results show that the recognition rate is above 97.7%, and that approach is effective for classifying event types for small samples.


Introduction
Infrasound is a type of low frequency signal that is undetectable to the human ear, ranging in frequency from 0.01 to 20 Hz.Many events produce infrasound signals, including natural phenomena such as tsunami, aurora, volcanic eruptions, and earthquakes and human activities such as nuclear blasts, rockets, and artillery shooting [1][2][3][4][5].Because various events produce infrasound by different mechanisms, the energy of the signals is also distributed in different frequency [6][7][8][9].We can analyze the characteristics of detected infrasound signals to complete classification, which is the technical basis of the global infrasound detection system, one of the four types of monitoring used by the Comprehensive Nuclear-Test-Ban Treaty (CTBT) International Monitoring System (IMS).
Infrasound classification consists of two parts: feature extraction and classification recognition.The key to classification of infrasonic events is the extraction of effective feature vectors from a signal.Feature vectors must be distinguished from other targets.Effective feature extraction techniques are the foundation of the classification of infrasound events.After feature extraction of the target signal, feature vectors must be processed by a classifier.Although the feature vectors largely determine classification, the performance of the classifier also directly affects the classification results.
Both natural events and human activities produce nonlinear and nonstationary infrasound signals.Much research has focused on extracting feature vectors accurately from such signals and using them effectively for classification.Tsybul'skaya et al. [10] describe possible methods for the classification of infrasonic signals from different sources.By using wavelet packet decomposition, Wang and Li classified the acoustic signals produced by fish and ships [11,12]; Ham and others applied the cepstrum coefficient method, the process used in speaker recognition, to classify infrasound signals produced by volcanoes and waves associated with mountains [13][14][15][16]; and Chilo et al. compared three feature extraction techniques [17,18]discrete wavelet transform (DWT), time scale spectrum using continuous wavelet transforms (TSS), and cepstral coefficients and their derivatives-and used two neural networks for classification.Their results show that the TSS method reflects the characteristics of infrasound signals effectively.
HHT, a new method for time-frequency analysis, performs better than other signal processing methods for analyzing nonlinear and nonstationary signals [19][20][21].In this paper, we introduce HHT into the field of infrasound classification.The feature vectors extracted by HHT are high-dimensional.Because SVM is an efficient classifier for analyzing highdimensional data, we use SVM as a classifier.We propose a classification method for infrasound signals based on HHT and SVM.A comparison experiment was performed to test the efficacy of the proposed method compared with the advanced hybrid methodology by José Chilo et al.For this comparison experiment, we performed different classification methods on the same data.The simulation results show that the feature vectors extracted by the method in this paper have a relative good degree of differentiation and the classification accuracy of infrasound signals is improved over other methods.

Feature Extraction Based on HHT.
In 1998, Huang and his team proposed a new signal analysis method, Hilbert-Huang transforms (HHT) [20].HHT is based on concept of the intrinsic mode function (IMF) and the empirical mode decomposition (EMD) method.After the EMD process, the time series signal () is decomposed into several IMF imf  () components and a residue, as shown in the following formula: Here,   () is the residue component.Using the Hilbert transform, the instantaneous frequency of a signal can be calculated from the IMF.After performing the Hilbert transform on every IMF component, the Hilbert spectrum, (, ), may be obtained.The marginal spectrum of the signal, ℎ(), can be obtained by integration: HHT is a set of complete signal processing methods.The EMD technique allows adaptive decomposition according to the characteristics of the signal without any basis function setting in advance.This is an essential difference compared to Fourier transform and wavelet transforms.Hilbert spectrum can clearly reflect the energy distribution of a signal as a function of time and frequency.The Hilbert marginal spectrum, which is the integral of the instantaneous frequency, shows the total energy contribution of each frequency.It represents the cumulative probabilistic amplitude over the entire dataset and the actual frequency distribution of the signal.
In this study, the original signal is decomposed into several IMF components by applying the EMD technique to several kinds of typical infrasound signals.The Hilbert transform of each IMF component was performed, the marginal spectrum was calculated, and the real energy distribution  characteristics of the signal in the frequency domain were obtained.The feature vectors from the energy distribution of several kinds of events were extracted to serve the basis of postclassification.
2.2.Classification Using SVM.SVM, which was first proposed by Corinna Cortes and Vapnik in 1995 [22], has many advantages for classification and regression prediction in solving small sample, nonlinear, and high dimension problems [23,24].The first issue to be considered is the regression function for data fitting when the SVM is applied for classification.Let (  ,   ) be the sample data, where  = 1, 2, . .., .Here  represents the number of samples, and  is the desired output.Through training and learning, a regression function can be obtained from the sample data: Here  is the normal vector of the classification hyperplane and  is the offset of the hyperplane.The optimal regression function is the solution of quadratic optimization.The objective function that requires optimization is as follows: The constraints are Here only a small portion of  and  * will be nonzero, and they are the support vectors to the corresponding samples.To solve the equations above, replace (  ,   ) with the kernel function (  ,   ) in the above formula to obtain  and the regression function: The key to the SVM is the kernel function.Different kernel function leads to different SVM algorithm [25].Kernel functions that are used commonly are the linear kernel function, the polynomial kernel function, the radial basis function (RBF), and the sigmoid kernel function [26].Here we use the RBF kernel function.RBF kernel function is widely used because SVM with RBF kernel function has less RBF parameters to be determined and can get a very smooth estimate and relatively good generalization ability [26][27][28][29].The RBF kernel function is shown as follows, where  2 represents the variance of the test sample: 3. Simulations   around the world.All codes used in this study were written in MATLAB.

Simulations Setup.
The volcano, tsunami, and earthquake signals exhibit strong spectral characteristics in the 0.01 Hz-1.0 Hz range [14].The signals presented to the classifier are restricted to this frequency range by a fourth order Butterworth bandpass filter with cutoff frequencies of 0.01 Hz and 1.0 Hz.After filtering and denoising the signals, the HHT method was used for signal processing and acquisition of the marginal spectrum.The marginal spectrum obtained by HHT reflects the real distribution of infrasound signals in the frequency domain.The feature vector comprised of the amplitude of the marginal spectrum demonstrates a clear distinction between infrasound events.Because the energy of the signal after 120 is nearly zero, as shown in Figure 2, the feature vectors extracted from the marginal spectrum have 120 dimensions.By feature extraction, three classes of a total of 139 feature vectors are obtained.Each class is randomly divided into two groups: the training group and the testing group.The proportion of training group data to test group data is about three to one.First, the training group is used to train the SVM classification model.Then, the testing group is used to test the SVM classification model.Finally, the classification results and accuracy are given.This process is shown in Figure 1.

Simulation Results and Discussion
. Figure 2 shows infrasound feature vectors of three different classes: earthquake, tsunami, and volcano.They were obtained using the method described previously.These results show that there is a high degree of similarity among the feature vectors in each class, while there are clear differences between classes.These features may be helpful for classification and may improve the accuracy of classification.
The infrasound classes used for training and testing are shown in Table 2.The SVM model was trained with the training set data, which provided the optimal parameters of the SVM classifier.The residual vector was used as a test set to test the performance of the SVM classifier.The final classification result is shown in Figure 3. On the graph, the abscissa represents feature vectors for each test and the ordinate represents the category of test feature vector: "∘" instead of the actual category of test feature vector; " * " instead of the prediction category of test feature vector.The results show that the rate of correct classification is 97.7% for the three infrasound events.
To compare SVM with other methods, we extracted feature vectors using the DWT method for the same data [14] and performed SVM for classification, as shown in Figure 4. Compared with DWT algorithm, feature vectors extracted by HHT can better reflect the characteristics and frequency distribution of the infrasound signal.The method proposed has relatively good results of identification, as shown in Table 3.

Conclusion and Future Work
By performing a large number of tests, the method described above has proved to be effective, and the classification accuracy is above 97%.The results show that the method proposed in this study performs better than other methods, such as feature extraction using discrete wavelet transforms and classification using radial basis function neural networks when applied to the classification of infrasound signals.We conclude that this new method has a certain reference value for feature extraction and classification recognition of infrasound signals.
Due to the limitation of current conditions, only small sample sizes and a small number of infrasound types were used for these tests.These factors will influence the reliability of the results.In order to obtain more accurate results, more infrasound data and infrasound event types should be analyzed.

5 Figure 2 :
Figure 2: Infrasound feature vectors of the three different classes extracted by HHT.

Table 1 :
Infrasound data summary a .
a Data from Comprehensive Nuclear-Test-Ban Treaty Beijing National Data Center.b The station codes of IMS.

Table 2 :
Infrasound classes used for training and testing.

Table 3 :
The comparison of infrasound classification system.