Research Article  Open Access
Mechanical Fault Diagnosis Using Color Image Recognition of Vibration Spectrogram Based on Quaternion Invariable Moment
Abstract
Automatic extraction of timefrequency spectral image of mechanical faults can be achieved and faults can be identified consequently when rotating machinery spectral image processing technology is applied to fault diagnosis, which is an advantage. Acquired mechanical vibration signals can be converted into color timefrequency spectrum images by the processing of pseudo WignerVille distribution. Then a feature extraction method based on quaternion invariant moment was proposed, combining image processing technology and multiweight neural network technology. The paper adopted quaternion invariant moment feature extraction method and gray levelgradient cooccurrence matrix feature extraction method and combined them with geometric learning algorithm and probabilistic neural network algorithm, respectively, and compared the recognition rates of rolling bearing faults. The experimental results show that the recognition rates of quaternion invariant moment are higher than gray levelgradient cooccurrence matrix in the same recognition method. The recognition rates of geometric learning algorithm are higher than probabilistic neural network algorithm in the same feature extraction method. So the method based on quaternion invariant moment geometric learning and multiweight neural network is superior. What is more, this algorithm has preferable generalization performance under the condition of fewer samples, and it has practical value and acceptation on the field of fault diagnosis for rotating machinery as well.
1. Introduction
Quality of products and safety of entire production system are affected by working condition of rotary machineries directly. The fault diagnosis is performed and the working condition is monitored so that faults can be detected as early as possible, which is a key issue in the field of engineering. Moreover, signal processing, extractions of the fault features, and pattern recognition play a key role in the process of realization. At present, working state analysis based on vibration signal of machine parts is a common research method of the rotary machinery diagnosis and monitoring. There are many deficiencies of traditional methods of fault diagnosis. For example, timedomain analysis supplies limited information, so it can be only applied to simple machinery. Frequencydomain analysis is based on an assumption of stationary signals, and it analyzes the frequencydomain information of signals but it neglects timedomain characteristics of signals. However, actual signals of mechanical faults are usually nonstationary. In other words, the frequency of signals always changes with time, so timedomain and frequencydomain analysis are no longer suitable to handle complex fault diagnosis. Timefrequency analysis method is an effective method to cope with nonstationary and nonlinear signals. It includes short time Fourier transform [1], Hilbert timefrequency spectrum analysis [2], WignerVille distribution [3], Wavelet analysis [4], HHT timefrequency analysis, and high order spectral analysis [5, 6]. Timefrequency methods described above can solve problems of fault signal processing and features extraction successfully, but it is hard to analyze potential weak fault feature as vibration spectrum images are too subjective and the observation ability of operators varies. Only to use the combination of image features extraction and model recognition can avoid the subjective mistakes. Furthermore, the intelligent mechanical fault diagnosis is demanded urgently by the market, and researches to new feature extractions and recognition methods become a common trend.
With the development of information theories and machinevision, fault diagnosis technologies based on image processing has improved a lot. The focus on intelligent image processing is how to exact information to characterize the image features from images. The methods of features extraction which are related to the vibration spectral images usually can be divided into three categories: spectral feature, texture feature, and shape feature. The spectral feature recognition is on the basis of grey level histogram with single spectral band or color histogram with multispectral bands. It is simple, intuitive, and robustness. Zhou et al. [7] and others identified the voice signals successfully according to combinations of spectral features and prosodic features. Combined with spectroscopy, structural, and semantic features, the classification for highresolution remote sensing images is performed in the literature [8], and it resulted significantly with the support vector machine (SVM). However, different objects may have the same spectral feature, so these algorithms cannot describe specific object accurately. The texture features describe the spatial changes situation of the color, the gradation or the fine structure, and shape of the images [9]. Texture feature extraction proposed in literature [10] considered global nature of images, so it can describe regional features steadily. In addition, it has favorable rotation invariance and antiinterference capability. Haralick et al. [11] firstly proposed fourteen feature parameters to show the gray information of the images based on the texture features. The literature [12] effectively realized the image classification by utilizing the wavelet transform to extract the texture feature of the SAR images, selecting wavelet energy feature kurtosis which constitutes the feature vector for the SVM classifier. In literature [13], threedimensional discrete wavelet texture was utilized to realize the image classification and achieved better result than the spectralspatial classification method widely used. The disadvantages are that the more the changes of image pixel resolution, the more the deviation of texture feature. Shape features are not only an inherent feature of the rigid body but also an important feature of human object recognition. Shape features can be expressed by form factors and profile factors. Literature [14] applied a feature extraction based on anomaly detection and the local shape features to the high resolution video images of the vessel detection. In addition, some researchers have proposed some image shape feature detection operators. Gool et al. [15] constructed a line moment feature, which is not effected by translation and rotation; Belongie et al. [16] proposed SC (Shape Context) features for detecting edge points; Fergus et al. [17] used the Canny algorithm to select the image edge; Kadir et al. [18] used to calculate the color histogram of the local maximum entropy, and circular the feature extraction area. Shape feature extraction is researched on the necessary part of the image and the integrality of the target is grasped well. The disadvantage is that when the image’s target is deformed, the stability of the characterization would decline [19]. However, there is no target deformation on the mechanical vibration’s timefrequency images that we researched on; consequently, it is appropriate for the field of mechanical fault diagnosis. The “,” “,” and “” channels of color image can be processed independently as three gray images. However, the relationships between trichannels are ignored, and the processed image will get distortion. Quaternion is a branch of nonEuclidean geometry and is firstly put forward by mathematician Hamilton. In literature [20], some researchers raised up a color images processing method based on the quaternion. Compared with traditional methods, it can sufficiently describe the color relations of the images. Literature [21] rose up a quaternion geometric moment calculation method of the medical images to realize the registration of the skull images.
Images recognition technology is one of the hot spots in the area of artificial intelligence, which has the intuitive and extensive statistical significances. When this method is applied to the faults diagnosis, not only human experience evaluation can be liberated, but also accurate and timely fault warning can be achieved to prevent faults. Timefrequency analysis of the vibration signals can get 2dimensional or 3dimensional timefrequency distribution images which reflect the features, and image recognition technology can be applied to the color timefrequency images [22]. Image recognition methods are usually divided into statistical method [23], syntactic recognition method [24], neural networks method [25], template matching method [26], and geometrical transformation method [27]. Furthermore, multiagent system is used in image recognition widely [28–30]. Multiweight neural network is utilized to identify the faults in the paper. This method especially has unique advantages on solving the pattern recognition problems for the few samples, and it successfully solves a lot of few samples, nonlinearity, and high dimensional and other problems. It has the practical significance on the solutions of faults diagnosis.
This paper utilized image recognition technology to process the mechanical vibration spectral images. The main idea is to analyze mechanical fault signals by using pseudo WignerVille distribution method and explore the feature information of different mechanical faults in timefrequency domain. In order to prevent the color image information losing and ensure the integrality and relevance of the image’s triphosphor, the images’ invariant moment in quaternion domain has been analyzed; at last quaternion invariable moment features were sent to multiweights neural network for training and testing, and recognition results were got by using unknown samples’ test. The experimental result proved that modeling for color images by using quaternion invariable moment created a new method for color images processing, and this method has some practical values as well.
2. Generation of Mechanical Vibration TimeFrequency Spectral Image
Methodology of fault diagnosis based on mechanical vibration timefrequency color image recognition mainly consists of three steps: (1) generating of vibration spectral color image reflecting mechanical vibration fault’s feature; (2) extracting features of vibration spectral image accurately; (3) recognizing the extracted vibration spectral feature effectively. WignerVille distribution can accurately define timefrequency structure of signals [31]. For a certain signal denoted here as , the WignerVille distribution can be defined asThere is interference of cross terms in the transformation expressed as (1), and the cross terms will not disappear with the changes of timefrequency distance. In order to eliminate the interference of these cross terms, windowing processing is applied to (1), and a novel transformation, denoted here as pseudo WignerVille distribution, can consequently be obtained and expressed as follows:Wherein, of (2) should meet some requirements listed as follows [5]:(I)Symmetry is .(II)Normalization is .(III)The length of window is limited, which can be set as , and then for any , is right.(IV)Fourier Transform for is a lowpass function.
Here, is selected as rectangular window. Because actual collected signal is discrete, should be discredited. Considering the length of denoted as is limited, that is, for any , it has , and then (2) can be adapted asEquation (3) is discredited in time domain; that is, , , where represents sampling frequency.
Then (3) can be adapted asThe length of discrete window is .
Equation (4) is discretized in frequency domain; that is, , , and ; then (4) can be adapted asEquation (5) is a formulation of discrete pseudo WignerVille distribution. Pseudo WignerVille distribution inherits good performance in timefrequency concentration of WignerVille distribution. In addition, it can suppress crossterms interference and has higher timefrequency resolution. So it can reflect timefrequency feature information of vibration signals of rotation machinery.
3. Vibration Spectrum Feature Extraction
From the perspective of image pattern recognition analysis, it is a new way to extract effective vibration spectrum characteristics and combine with intelligent recognition method. To solve the key technology of image features extract and reflect accurately in the fault diagnosis of vibration spectra, this paper used quaternion invariant moment and gray levelgradient cooccurrence matrix for image feature extraction.
3.1. Quaternion Geometric Moment of Image
When different faults happen, their vibration signals have different vibration information. In terms of the same fault, timefrequency distributions of the vibration signal in different cycle period are similar. However, under different fault conditions, timefrequency distribution of vibration varies significantly, especially in the features extraction method of quaternion invariable moment. Taking advantage of this nature, effective features can be extracted from the timefrequency spectral images of rotation machinery.
Quaternion theory [32] is an extension of plurality basic theory, which consists of onedimensional real part together with threedimensional imaginary part, and can be expressed aswhere , and are real numbers and , and are quaternion’s imaginary part element. If real part , is a pure imaginary quaternion and meets Hamilton rules; that is,Conjugate of quaternion is named as and can be expressed asThe mold of quaternion is denoted as and can be expressed asColor images can commonly be described as combinations of three primary colors, that is, Red, Green, and Blue (RGB) [33]. In the present work, three primary values of color image pixel values are represented by pure imaginary quaternion; therefore, any pixel denoted here as can be expressed aswhere , , and represent the channels’ values of , , and of image coordinates ; therefore, a color image can be expressed as a purely imaginary quaternion matrix. Compared with traditional methods, quaternion model can better reflect the integrity and relevance of the color image.
Regarding a timefrequency distribution of certain vibration for rotation machinery, its essence is a twodimensional color image, expressed here as . From the uniqueness theorem proposed by Hu in literature [34], it can be deduced that if a twodimensional image is piecewise continuous and has nonzero value in limited area, then any order moment will exist. Set of quaternion moment can be uniquely identified by ; otherwise, twodimensional image can also be determined by set of quaternion moment ; that is, the set of quaternion moment that we obtained can describe information of this timefrequency image partly. Moreover, based on Hu, and combining quaternion as well, taking full advantage of color information of images, quaternion invariant moment reflects the integrity and relevance of the image.
Given a twodimensional color image , then the corresponding twodimensional geometric quaternion moment denoted here as can be indicated aswhere indicates the order of , and it exists by , , , and .
The image quaternion geometric moment is an expansion of grayscale images geometric moment concept [35]. is a quaternion vector and contains images’ color information. It not only describes global features of the image shape but also provides plenty of geometric feature information of images. Each order moment has different physical meanings.
3.1.1. Quaternion ZerothOrder Moment
In (11), if , is quaternion zerothorder moment and represents color sum of entire image. We can consider it as the “mass” of the image; moreover, this “mass” takes information of both size and direction.
3.1.2. Quaternion FirstOrder Moment
, are firstorder moment of the color image regarding axis and axis, and color moment center denoted here as can be expressed as where indicates module of the quaternion moment, which expresses color intensity, and coordinate is a physical centroid of the image.
3.1.3. Quaternion Central Moment
The geometric moment calculated with respect to centroid of the image is called central moment and expressed aswhere is the centroid of the image, and , . The secondorder moments are a measure of variance of the image intensity distribution about the origin. The central moments , represent the variance about the mean (centroid ). is the covariance measure of them.
In order to make with scale and rotation invariance [36], (13) should be transferred to In 1962, Hu firstly proposed invariant moment in literature. Seven invariant moments are also a nonlinear combination of secondorder and thirdorder central moment given in (14). They have performance of invariance in transformation of translation and rotation and can be expressed as follows:where matrix is a matrix extracted from a characteristic matrix of matrix with respect to each timefrequency color image.
3.2. Gray LevelGradient Cooccurrence Matrix Feature Extraction
The element of gray levelgradient cooccurrence matrix is defined as gray value and grade value pixels between normalization gray image and normalization grade image. In the process of normalization gray levelgradient cooccurrence matrix, the sum of each element is expressed as follows:And ; the above formula can be expressed as follows: where is the element of gray levelgradient cooccurrence matrix, is gray value, and is grade value.
Exacted 15 feature statistics, small gradient advantage, large gradient advantage, gray distribution heterogeneous, gradient distribution heterogeneous, energy, the average of gray, the average of gradient, gray mean square error, gradient mean square error, correlation, gray entropy, gradient entropy, mixed entropy, inertia, and inverse difference moment [37], are recorded as , and , respectively, and expressed in Table 1.

Gray image texture feature vector is composed of the 15 numerical characteristics in Table 1:
4. Construction of Multiweights Neural Network
Main difficulty in intelligent diagnosis of mechanical faults is lack of typical fault samples which consequently result in quantity imperfectness of training samples. Other difficulties are shortage of faults characteristic knowledge. These two difficulties restrict the development of mechanical faults diagnosis. Multiweights neural network takes the best coverage of similar samples in specific space as a target. All samples with respect to similar objects exist in an irregular connected area, and a highdimensional complex geometry in corresponding space is built to cover the area according to certain rules. If a sample is in the enclosed area recognition is correct, on the contrary, the recognition is rejected. In other words, the subspace structure of each type sample depends only on the sample type. Relationships among the samples of each type are analyzed based on highdimensional space complex geometry coverage recognition method [38]. The method described above is completely different with traditional method of pattern recognition, which depends on a subspace constructed according to relationships of various types of samples. Therefore, multiweights neural network is selected to identify the quaternion invariant moment features of rotation machine vibration spectrum images, and a multiweights neuron model in the multiweights neural network is illustrated in Figure 1.
A basic model of the general neurons is established as where are weight vectors. In present work, a matrix denoted here as extracted from color image according to the methodology in above chapter is selected as one weight vector. are input vectors, that is, feature vectors of samples to be identified. In this paper, test samples are also taken from . is quantity of test samples, is a function determined by multiweights neuron, and is a step function. This model describes not only relationships amongst sample points in feature space, but also relationships between testing samples and training samples.
Threeweight neuron denoted here as is utilized in present work; that is, value of the parameter is equal to 3. It can be expressed aswhere , and represents a triangle area constructed as three points, denoted, respectively, as , and . Then describes a coverage that is determined by a topological product of the and a hypersphere whose radius is denoted as .
Certain training sample set extracted from specific type of fault is supposed as , is total number of sample points, represents one of the samples, represents the dimensions of the sample point, represents the Euclidean distance from point to space , , and . Construction steps of the multiweights neural network are described as follow.
Step 1. Calculations of Euclidean distances between all the points in the training samples and denoted two points the Euclidean distance of them is minimum as and . Among the rest of the points in the sample set , a point to which the sum of distance from and is shortest is selected and denoted as . Three points described above can construct a closed region, denoted here as . A neuron is applied to cover this region, and the coverage is
Step 2. It is judged that whether the reset of points are covered by the geometry constructed in the previous step. Every point should be eliminated if it lies within the coverage of geometry. With respect to those points lying out of the geometry, we should find a point to which the sum of distance from , , and is shortest according to Step 1 and mark it as point . Besides, we calculate the Euclidean distances from to , from to , and from to , respectively, to facilitate description, then choose two of them with smaller and smallest values, and mark them as and . Thus, , , and construct the second closed region, denoted here as . Therefore, another neuron is applied to cover this region, and the coverage iswhere represents the distance from point to the closed region .
Step 3. Eliminate all points lying in the coverage of region and then calculate the minimum value of the sum of Euclidean distance from each of the rest of points to the group of three points that make up . We mark this point as . Similarly, we calculate the Euclidean distance from to the points that make up , choose two of the points with smaller and smallest values, and mark them as and . Points , , and construct the th closed region which is marked as . At last, a neuron is applied to cover this region, and the coverage is
Step 4. Repeat the operation described in Step 3 until the processing of all sample points is over. Ultimately, we can construct geometry and denote their combo as :Schematic diagram of fault recognition for rotating machinery is illustrated as Figure 2. Firstly, vibration signals that can reflect characteristic of mechanical faults are collected and preprocessed by pseudo WignerVille distribution, and a twodimensional color timefrequency image can subsequently be obtained. Secondly, methodology of quaternion invariable moment is utilized to extract character of the image obtained in previous step and a feature Matrix can consequently be obtained. Thirdly, different image feature matrixes are randomly selected as the training sets and test sets. Then, multiweights neural network is built and trained. At last, Euclidean distances from each point in testing set to multiweights neuron are calculated. According to methodology described above, faults can be recognized.
5. Case Study of Two Methods of Features Extraction
Timefrequency images reflect vibration of rotating machinery under the condition that the frequency changes at any time. In present work, bearing faults of rotation machines are analyzed by utilizing pseudo WignerVille distribution. Experimental data come from bearing fault research center of Case Western Reserve University of United States. We chose 3 hp motor, rotation speed is 1725 rpm, and bearing to be detected at the drive end is SKF 6205. And the installed acceleration sensor at the bearing seat of drive end is used to collect the vibration acceleration signal. Figure 3, respectively, gives timefrequency images with respect to inner race fault, outer race fault, normal bearing, and ball bearing fault.
In the present work, image processing technology is utilized to find characteristic information from the image. Figure 3 is illustrated for distinguishing the faults. Traditional images processing technology of color image is a continuation of grayscale images. In which, color image is divided into three components, which are processed, respectively. Unfortunately, relationships amongst three components are ignored, so some information of the color image will be lost. In the present work, quaternion invariant moment theory is firstly utilized to the extraction of mechanical faults features. A pure imaginary quaternion is utilized to indicate three colors (Red, Green, and Blue) of an image and loss of images’ color information can consequently be avoided. Furthermore, invariant moment is an effective description for twodimensional image shape features, and it describes images’ constant feature when the image changes on the plane. Quaternion invariant moment makes effective description on color and shape features of color images. Different samples with respect to the same fault can be obtained after the transformations of translation, flex, and rotation to one image. However geometric moments are an effective method to deduce translation, flex, and rotation invariant. We selected three samples of four types of bearing and got their pseudo WignerVille timefrequency image and subsequently calculated their quaternion invariant moment. Result is listed in Table 2. It can be found from this table that seven quaternion invariant moments of the same faults close to each other; however, quaternion invariant moments of different faults (~) change significantly and can be selected as the features of the fault recognition. Texture feature data from four types of bearing are shown in Table 3.


6. Fault Recognition
Every sample’s timefrequency image can be obtained by pseudo WignerVille distribution, and quaternion invariant moment of all the samples can be calculated as fault characteristics of the late recognition. Twenty samples are selected randomly as training set in every kinds of fault, the remaining twenty samples as test set. A diagnostic test was carried out on the test sample, and the classification result is shown in Table 4. At the same time, texture features of timefrequency gray image were carried out to recognize faults based on geometrical study and we got test results in Table 5.


We can know that normal bearing and rolling bearing fault are all correctly identified in Table 4, and there are two inner race faults, and an outer race fault was identified incorrectly as rolling ball fault. The correct recognition rate is 96.25% in entire identification process. We can know that outer race fault and rolling bearing fault are all correctly identified in Table 5, and the correct recognition rate is 95% in entire identification process. The experimental results show that recognition rate of quaternion invariant moment based on geometric learning is higher than gray levelgradient cooccurrence matrix based on geometric learning.
In additional, we also used the same training and test samples to carry out the probabilistic neural network (PNN) identification test. PNN is a feedforward neural network; when it classified through pattern, the structure of PNN generally consists of input layer, hidden layer, and output layer. It is classified by using bias minimum risk rule to carry out parallel processing and classification estimation operation. The recognition result of using quaternion invariant moment feature extraction method was shown in Figure 4. “1” means “normal bearing,” “2” means “inner race fault,” “3” means “outer race fault,” and “4” means “rolling ball fault.” Seven samples did not identify in Figure 4, and the recognition rate is 91.25%. The recognition results of using gray levelgradient cooccurrence matrix feature extraction method were shown in Figure 5. “1” means “inner race fault,” “2” means “outer race fault,” “3” means “normal bearing,” and “4” means “rolling ball fault.” Four inner race faults were identified incorrectly, sixnormal bearing was identified as others in Figure 5, and the recognition rate is 87.5%.
In order to test the stability and practicality of geometric learning recognition method with few samples, the training samples in turn reduced 20%, and test sample still took the original 20 samples to identify. Four kinds of recognition results are shown in Table 6.

With the increase of the number of training samples, the known fault recognition rates also grow in Figure 6. The rates of geometric learning algorithm are higher than probabilistic neural network algorithm in the same feature extraction method. The recognition rates of quaternion invariant moment are higher than gray levelgradient cooccurrence matrix in the same recognition. Quaternion invariant moment combined geometric learning, which improved recognition rates.
7. Conclusion
Mechanical fault diagnosis is a complex process and it mainly includes vibration signal of the fault acquisition, the fault information of feature extraction, and pattern recognition. In order to prevent losses of color image information, considering the relevance of triphosphor, the timefrequency image information would be processed in the quaternion field, and invariant moment features are also extracted in the quaternion field. From this perspective, the paper took fault diagnosis of rolling bearing as an example and used pseudo WignerVille distribution timefrequency to get timefrequency spectrum of fault signals. Color image feature extraction method based on quaternion invariant moment was presented creatively and expressed the information of color image and shape characteristics of the image effectively. In addition, the paper compared quaternion invariant moment with gray levelgradient cooccurrence matrix. In order to avoid identifying the error of vibration spectrum through experience, we applied geometric learning method to intelligent mechanical diagnosis. The superiority of geometric learning algorithm was verified by comparing with probabilistic neural network algorithm. It was useful to combine the feature extraction method of quaternion invariant moment with geometric learning algorithm through the simulation results. And it could distinguish normal bearing and fault types. It had a good recognition result in reducing the number of training samples.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (61273024, 61305031, and 51307089), the Natural Science Foundation of Jiangsu Province (KB2012227), and the “333” and “226” Highlevel Personnel Training Project, a Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).
References
 R. De Frein and S. T. Rickard, “The synchronized shorttimefouriertransform: properties and definitions for multichannel source separation,” IEEE Transactions on Signal Processing, vol. 59, no. 1, pp. 91–103, 2011. View at: Publisher Site  Google Scholar  MathSciNet
 H. Li, S. Kwong, L. Yang, D. Huang, and D. Xiao, “Hilberthuang transform for analysis of heart rate variability in cardiac health,” IEEE/ACM Transactions on Computational Biology and Bioinformatics, vol. 8, no. 6, pp. 1557–1567, 2011. View at: Publisher Site  Google Scholar
 V. ClimenteAlarcon, J. A. AntoninoDaviu, M. RieraGuasp, and M. Vlcek, “Induction motor diagnosis by advanced notch FIR filters and the wignerville distribution,” IEEE Transactions on Industrial Electronics, vol. 61, no. 8, pp. 4217–4227, 2014. View at: Publisher Site  Google Scholar
 K. Huang and S. Aviyente, “Wavelet feature selection for image classification,” IEEE Transactions on Image Processing, vol. 17, no. 9, pp. 1709–1720, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 J. Wang and Y. Pi, “SAR tomography imaging via higherorder spectrum analysis,” Journal of Systems Engineering and Electronics, vol. 20, no. 4, pp. 748–754, 2009. View at: Google Scholar
 A. Kachenoura, L. Albera, J.J. Bellanger, and L. Senhadji, “Nonminimum phase identification based on higher order spectrum slices,” IEEE Transactions on Signal Processing, vol. 56, no. 5, pp. 1821–1829, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Zhou, Y. Sun, J. Zhang, and Y. Yan, “Speech emotion recognition using both spectral and prosodic features,” in Proceedings of the International Conference on Information Engineering and Computer Science (ICIECS '09), pp. 1–4, December 2009. View at: Publisher Site  Google Scholar
 X. Huang and L. Zhang, “An SVM ensemble approach combining spectral, structural, and semantic features for the classification of highresolution remotely sensed imagery,” IEEE Transactions on Geoscience and Remote Sensing, vol. 51, no. 1, pp. 257–272, 2013. View at: Publisher Site  Google Scholar
 H. Tamura, S. Mori, and T. Yamawaki, “Textural features corresponding to visual perception,” IEEE Transactions on Systems, Man and Cybernetics, vol. 8, no. 6, pp. 460–473, 1978. View at: Publisher Site  Google Scholar
 H. Lategahn, S. Gross, T. Stehle, and T. Aach, “Texture classification by modeling joint distributions of local patterns with Gaussian mixtures,” IEEE Transactions on Image Processing, vol. 19, no. 6, pp. 1548–1557, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 R. M. Haralick, K. Shanmugam, and I. Dinstein, “Textural features for image classification,” IEEE Transactions on Systems, Man and Cybernetics, vol. 3, no. 6, pp. 610–621, 1973. View at: Publisher Site  Google Scholar
 G. Akbarizadeh, “A new statisticalbased kurtosis wavelet energy feature for texture recognition of SAR images,” IEEE Transactions on Geoscience and Remote Sensing, vol. 50, no. 11, pp. 4358–4368, 2012. View at: Publisher Site  Google Scholar
 X. Guo, X. Huang, and L. Zhang, “Threedimensional wavelet texture feature extraction and classification for multi/hyperspectral imagery,” IEEE Geoscience and Remote Sensing Letters, vol. 11, no. 12, pp. 2183–2187, 2014. View at: Publisher Site  Google Scholar
 Z. Shi, X. Yu, Z. Jiang, and B. Li, “Ship detection in highresolution optical imagery based on anomaly detector and local shape feature,” IEEE Transactions on Geoscience and Remote Sensing, vol. 52, no. 8, pp. 4511–4523, 2014. View at: Publisher Site  Google Scholar
 L. Gool, T. Moons, and D. Ungureanu, “Affine/photometric invariants for planar intensity patterns,” in Computer Vision—ECCV '96: 4th European Conference on Computer Vision Cambridge, UK, April 15–18, 1996 Proceedings, Volume I, vol. 1064 of Lecture Notes in Computer Science, pp. 642–651, Springer, Berlin, Germany, 1996. View at: Publisher Site  Google Scholar
 S. Belongie, J. Malik, and J. Puzicha, “Shape context: a new descriptor for shape matching and object recognition,” in Proceedings of the Neural Information Processing Systems (NIPS '00), pp. 831–837, 2000. View at: Google Scholar
 R. Fergus, P. Perona, and A. Zisserman, “A sparse object category model for efficient learning and exhaustive recognition,” in Proceedings of the IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR '05), vol. 1, pp. 380–387, IEEE, June 2005. View at: Publisher Site  Google Scholar
 T. Kadir, A. Zisserman, and M. Brady, “An affine invariant salient region detector,” in Computer Vision—ECCV 2004: 8th European Conference on Computer Vision, Prague, Czech Republic, May 11–14, 2004. Proceedings, Part I, vol. 3021 of Lecture Notes in Computer Science, pp. 228–241, Springer, Berlin, Germany, 2004. View at: Publisher Site  Google Scholar
 G. Kumar and P. K. Bhatia, “A detailed review of feature extraction in image processing systems,” in Proceedings of the 4th International Conference on Advanced Computing and Communication Technologies (ACCT '14), pp. 5–12, IEEE, Rohtak, India, February 2014. View at: Publisher Site  Google Scholar
 G. K. Rohde, A. Aldroubi, and D. M. Healy Jr., “Interpolation artifacts in subpixel image registration,” IEEE Transactions on Image Processing, vol. 18, no. 2, pp. 333–345, 2009. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Huang, L. Hua, H. Feng, L. Ding, and Y. Chen, “Color medical image registration based on quaternion moment theory,” OptoElectronic Engineering, vol. 40, no. 3, pp. 102–107, 2013. View at: Publisher Site  Google Scholar
 S. Uhlmann and S. Kiranyaz, “Integrating color features in polarimetric SAR image classification,” IEEE Transactions on Geoscience and Remote Sensing, vol. 52, no. 4, pp. 2197–2206, 2014. View at: Publisher Site  Google Scholar
 R. S. H. Istepanian, L. J. Hadjileontiadis, and S. M. Panas, “ECG data compression using wavelets and higher order statistics methods,” IEEE Transactions on Information Technology in Biomedicine, vol. 5, no. 2, pp. 108–115, 2001. View at: Publisher Site  Google Scholar
 G. Badr and B. J. Oommen, “On optimizing syntactic pattern recognition using tries and AIbased heuristicsearch strategies,” IEEE Transactions on Systems, Man, and Cybernetics B: Cybernetics, vol. 36, no. 3, pp. 611–622, 2005. View at: Publisher Site  Google Scholar
 S. Razavi and B. A. Tolson, “A new formulation for feedforward neural networks,” IEEE Transactions on Neural Networks, vol. 22, no. 10, pp. 1588–1598, 2011. View at: Publisher Site  Google Scholar
 S. Omachi and M. Omachi, “Fast template matching with polynomials,” IEEE Transactions on Image Processing, vol. 16, no. 8, pp. 2139–2149, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 M. Holden, “A review of geometric transformations for nonrigid body registration,” IEEE Transactions on Medical Imaging, vol. 27, no. 1, pp. 111–128, 2008. View at: Publisher Site  Google Scholar
 Y. Tang, H. Gao, W. Zhang, and J. Kurths, “Leaderfollowing consensus of a class of stochastic delayed multiagent systems with partial mixed impulses,” Automatica, vol. 53, pp. 346–354, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Tang, Z. Wang, H. Gao, H. Qiao, and J. Kurths, “On controllability of neuronal networks with constraints on the average of control gains,” IEEE Transactions on Cybernetics, vol. 44, no. 12, pp. 2670–2681, 2014. View at: Publisher Site  Google Scholar
 W. Zhang, Y. Tang, Q. Miao, and J.A. Fang, “Synchronization of stochastic dynamical networks under impulsive control with time delays,” IEEE Transactions on Neural Networks and Learning Systems, vol. 25, no. 10, pp. 1758–1768, 2014. View at: Publisher Site  Google Scholar
 X.G. Xia and V. C. Chen, “A quantitative SNR analysis for the pseudo WignerVille distribution,” IEEE Transactions on Signal Processing, vol. 47, no. 10, pp. 2891–2894, 1999. View at: Publisher Site  Google Scholar
 D. Xu and D. P. Mandic, “The theory of quaternion matrix derivatives,” IEEE Transactions on Signal Processing, vol. 63, no. 6, pp. 1543–1556, 2015. View at: Publisher Site  Google Scholar  MathSciNet
 E. G. Karakasis, G. A. Papakostas, D. E. Koulouriotis, and V. D. Tourassis, “A unified methodology for computing accurate quaternion color moments and moment invariants,” IEEE Transactions on Image Processing, vol. 23, no. 2, pp. 596–611, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 M. K. Hu, “Visual pattern recognition by moment invariants,” IEEE Transactions on Information Theory, vol. 8, no. 2, pp. 179–187, 1962. View at: Publisher Site  Google Scholar
 Y. M. Wang, W. Yan, and S. Q. Yu, “Moment feature extraction of image based on Radon transform and its application in image recognition,” Computer Technology, vol. 27, no. 2, pp. 82–84, 2001. View at: Google Scholar
 Y. Wang, “The geometric moment and its invariants,” Journal of Shanghai Dianji University, vol. 9, no. 2, pp. 7–10, 2006. View at: Google Scholar
 J. Hong, “Gray levelgradient cooccurrence matrix texture analysis method,” Acta Automatica Sinica, vol. 10, no. 1, pp. 22–25, 1984. View at: Google Scholar
 S. Wang, Y. Liu, J. Lai, and X. Liu, Biomimetic Pattern Recognition and MultiWeight Neuron, National Defense Industry Press, Beijing, China, 2012.
Copyright
Copyright © 2015 Liang Hua et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.