Research article

Improved FCM algorithm for fisheye image cluster analysis for tree height calculation

  • Received: 01 July 2021 Accepted: 25 August 2021 Published: 09 September 2021
  • The height of standing trees is an important index in forestry research. This index is not only hard to measure directly but also the environmental factors increase the measurement difficulty. Therefore, the measurement of the height of standing trees is always a problem that experts and scholars are trying to improve. In this study, improve fuzzy c-means algorithm to reduce the calculation time and improve the clustering effect, used on this image segmentation technology, a highly robust non-contact measuring method for the height of standing trees was proposed which is based on a smartphone with a fisheye lens. While ensuring the measurement accuracy, the measurement stability is improved. This method is simple to operate, just need to take a picture of the standing tree and determine the shooting distance to complete the measurement. The purpose of the fisheye lens is to ensure that the tree remains intact in the photograph and to reduce the shooting distance. The results of different stability experiments show that the measurement error ranged from -0.196m to 0.195m, and the highest relative error of tree measurement was 3.05%, and the average relative error was 1.45%. Analysis shows that compared with previous research, this method performs better at all stages. The proposed approach can provide a new way to obtain tree height, which can be used to analyze growing status and change in contrast height because of high accuracy and permanent preservation of images.

    Citation: Jiayin Song, Yue Zhao, Zhixiang Chi, Qiang Ma, Tianrui Yin, Xiaopeng Zhang. Improved FCM algorithm for fisheye image cluster analysis for tree height calculation[J]. Mathematical Biosciences and Engineering, 2021, 18(6): 7806-7836. doi: 10.3934/mbe.2021388

    Related Papers:

    [1] Ming Chen, Menglin Gong, Jimin Zhang, Lale Asik . Comparison of dynamic behavior between continuous- and discrete-time models of intraguild predation. Mathematical Biosciences and Engineering, 2023, 20(7): 12750-12771. doi: 10.3934/mbe.2023569
    [2] Renji Han, Binxiang Dai, Lin Wang . Delay induced spatiotemporal patterns in a diffusive intraguild predation model with Beddington-DeAngelis functional response. Mathematical Biosciences and Engineering, 2018, 15(3): 595-627. doi: 10.3934/mbe.2018027
    [3] Adnan Sami, Amir Ali, Ramsha Shafqat, Nuttapol Pakkaranang, Mati ur Rahmamn . Analysis of food chain mathematical model under fractal fractional Caputo derivative. Mathematical Biosciences and Engineering, 2023, 20(2): 2094-2109. doi: 10.3934/mbe.2023097
    [4] Peter A. Braza . Predator-Prey Dynamics with Disease in the Prey. Mathematical Biosciences and Engineering, 2005, 2(4): 703-717. doi: 10.3934/mbe.2005.2.703
    [5] Jie Hu, Juan Liu, Peter Yuen, Fuzhong Li, Linqiang Deng . Modelling of a seasonally perturbed competitive three species impulsive system. Mathematical Biosciences and Engineering, 2022, 19(3): 3223-3241. doi: 10.3934/mbe.2022149
    [6] Jian Zu, Wendi Wang, Bo Zu . Evolutionary dynamics of prey-predator systems with Holling type II functional response. Mathematical Biosciences and Engineering, 2007, 4(2): 221-237. doi: 10.3934/mbe.2007.4.221
    [7] Moitri Sen, Malay Banerjee, Yasuhiro Takeuchi . Influence of Allee effect in prey populations on the dynamics of two-prey-one-predator model. Mathematical Biosciences and Engineering, 2018, 15(4): 883-904. doi: 10.3934/mbe.2018040
    [8] Yazhi Wu, Guangyao Tang, Changcheng Xiang . Dynamic analysis of a predator-prey state-dependent impulsive model with fear effect in which action threshold depending on the prey density and its changing rate. Mathematical Biosciences and Engineering, 2022, 19(12): 13152-13171. doi: 10.3934/mbe.2022615
    [9] Gianni Gilioli, Sara Pasquali, Fabrizio Ruggeri . Nonlinear functional response parameter estimation in a stochastic predator-prey model. Mathematical Biosciences and Engineering, 2012, 9(1): 75-96. doi: 10.3934/mbe.2012.9.75
    [10] Saheb Pal, Nikhil Pal, Sudip Samanta, Joydev Chattopadhyay . Fear effect in prey and hunting cooperation among predators in a Leslie-Gower model. Mathematical Biosciences and Engineering, 2019, 16(5): 5146-5179. doi: 10.3934/mbe.2019258
  • The height of standing trees is an important index in forestry research. This index is not only hard to measure directly but also the environmental factors increase the measurement difficulty. Therefore, the measurement of the height of standing trees is always a problem that experts and scholars are trying to improve. In this study, improve fuzzy c-means algorithm to reduce the calculation time and improve the clustering effect, used on this image segmentation technology, a highly robust non-contact measuring method for the height of standing trees was proposed which is based on a smartphone with a fisheye lens. While ensuring the measurement accuracy, the measurement stability is improved. This method is simple to operate, just need to take a picture of the standing tree and determine the shooting distance to complete the measurement. The purpose of the fisheye lens is to ensure that the tree remains intact in the photograph and to reduce the shooting distance. The results of different stability experiments show that the measurement error ranged from -0.196m to 0.195m, and the highest relative error of tree measurement was 3.05%, and the average relative error was 1.45%. Analysis shows that compared with previous research, this method performs better at all stages. The proposed approach can provide a new way to obtain tree height, which can be used to analyze growing status and change in contrast height because of high accuracy and permanent preservation of images.



    Bearing is an important component of rotating machinery and has been widely used in various fields such as automation and medical treatment. Due to the harsh and complex working environment, it has become one of the most vulnerable parts of mechanical equipment. Its running state directly affects the safety, reliability, and service life of the whole mechanical system ADDIN EN.CITE.DATA [1]. Therefore, it is extremely important to conduct fault diagnosis research on bearings. Bearing fault diagnosis is very complicated, especially since its early fault signal is extremely weak, which is easy to be submerged in the noise signal of other components and cannot be detected ADDIN EN.CITE.DATA [2]. However, the above-mentioned traditional vibration analysis method and current analysis method mainly classify the bearing faults at high speed, and it is difficult to accurately classify the vibration signals buried in noise and the current signals in the low-speed bearing faults. Compared with vibration analysis ADDIN EN.CITE.DATA [3,4,5] and current signal analysis ADDIN EN.CITE.DATA [6], AE signal has the following advantages: a) It is not affected by mechanical background noise; b) It is more sensitive to early and low-speed bearing faults; c) It is sensitive to the location of the fault; d) It can provide good trend parameters; e) It has the characteristics of high frequency and obvious frequency ADDIN EN.CITE.DATA [7]. Due to the above advantages, bearing fault diagnosis based on AE signals has been widely used, and some research results have been obtained. However, AE technology has the disadvantages of a large amount of time-domain data and difficulty in storage and processing. At the same time, high-frequency AE signals are accompanied by high sampling rates, and the combination of high sampling rates and massive data will bring greater challenges to the cost of acquisition hardware. This time-domain signal obtained based on the Nyquist sampling theorem has a large amount of redundant information, and then the workload of feature extraction and analysis for such a huge amount of data has increased exponentially. Therefore, fundamentally reducing the amount of data becomes an urgent problem to be solved by the AE signal fault diagnosis method.

    The CS theory proposed in recent years uses the transform space projection method to realize the compressed sampling of the original signal ADDIN EN.CITE.DATA [8,9,10]. Most of the information contained in the original signal is obtained with very few measured data, thereby greatly improve the data transmission efficiency and reduce the data storage space, which provides a new idea for bearing AE signal fault diagnosis. At present, the bearing fault diagnosis method based on CS has received widespread attention.

    There are many application research of compressed sensing fault diagnosis focus on the optimized feature extraction method. Literature [11] proposed a bearing fault diagnosis method based on CS and matching pursuit (MP) reconstruction algorithm. Literature [12] proposed a CS framework for characteristic harmonics to detect bearing faults. This method uses a compressed MP strategy to detect characteristic harmonics from sparse measurements under the condition of incomplete signal reconstruction. Literature [13] proposed a bearing fault diagnosis method based on CS. This method trains multiple over-complete dictionaries through the dictionary learning method. These dictionaries are effective for the sparse decomposition of signals in each specific state, while signals in other states cannot be sparsely decomposed. This method uses this characteristic to determine the fault state of the bearing. Literature [14] proposed a sparse representation classification strategy, which combined sparse representation with random dimensionality reduction to extract and classify the fault features of rotating machinery. Literature [15] proposed a fault diagnosis method based on sparse representation of time-frequency features, which can reconstruct the time-frequency features of fault signals from a small amount of compressed sampled data containing noise. Literature [16] proposed a method to directly extract the acoustic emission signal compression feature (AECF) from the CS data, and the AECF trend is used to evaluate the running state of the bearing. Literature [22] proposed a fault diagnosis method based on SAE, which combines CS and wavelet packet energy entropy.

    Another hot-spot about compressed sensing fault diagnosis focus on the design classification method. Literature ADDIN EN.CITE [17] uses three methods to process the compressed data and the compressed measured value is directly used as the input of the classifier to diagnose bearing faults. Literature [18] proposed a new intelligent classification method, which uses sparse over-complete features and a deep neural network (DNN) with an unsupervised feature learning algorithm based on sparse autoencoders (SAE) to classify bearing faults in compression measurement. Literature [19] proposed a method for bearing fault diagnosis in a fluctuating environment based on the theory of compressed sensing. The proposed method can effectively reduce the amount of data required for bearing diagnosis and maintain similar accuracy to the current method. Moreover, the reconstructed signal can be used for other fault diagnosis methods. Literature [20] proposed a bearing fault diagnosis method based on CS and heuristic neural networks. Literature ADDIN EN.CITE [21] proposed a new type of intelligent diagnosis method based on CS and deep learning. It uses the advantages of CS and deep learning to obtain high recognition accuracy with a small amount of measurement data.

    The characteristics of bearing fault diagnosis technology based on compression perception above were comprehensively analyzed. Based on the research in reference [16], this paper optimizes and improves the signal decomposition method and feature extraction method in the compressed domain. This method obtains the compressed signal by projecting the signal from the time domain to the compressed domain and is regarded as the original measurement signal; the wavelet packet decomposition of the measured signal is carried out using the wavelet packet change matrix in the compression domain. The characteristic parameters of each frequency band information obtained are taken as feature vectors; the PSO algorithm is used to optimize the weighted coefficient of the feature vector to obtain the fault feature with enhanced feature; the SVM is used to identify and classify bearing fault types. The accuracy and practicability of the method are verified by fault experiments. The contribution of this method is as follows:

    a. A fault diagnosis method combining CS technology and AE signals is proposed. Compressed sampling greatly reduces the amount of AE data and retains most of the effective information of the signal for fault diagnosis. The problem of high hardware requirements for the AE signal due to large data volume is solved.

    b. The transform matrix of wavelet packet decomposition in the compression domain is derived from the time domain wavelet packet decomposition matrix, which is used to perform wavelet packet decomposition of the compression signal. The fault features of the compression domain of different frequency bands are extracted.

    c. Feature weighting is adopted to deal with fault features in the compression domain. The PSO algorithm is used to optimize the weighting coefficient of the fault features. This method realizes the enhancement of weak fault features and avoids the uncertainty of diagnosis results caused by the subjective and one-sided selection of feature information of different frequency bands.

    d. Experimental verification shows that the proposed method has a more satisfactory diagnostic performance for bearing fault diagnosis under low-speed conditions compared with the traditional method.

    CS was proposed by Candes et al. ADDIN EN.CITE.DATA [8,9]. It makes use of the sparse characteristic of signals in a certain transformation domain to compress data while sampling the signal. The compressed data retains most of the information of the original signal, and the reconstruction algorithm uses the observation data to reconstruct the original signal. Most signals in nature are not sparse signals, but a specific transform domain Ψ can be found to make the signal sparse. CS theory points out that if the measured signal xRN is sparse in a certain transform domain Ψ, a measurement matrix ΦRM×N (MN) that is not related to the transform domain can be used to linearly project the measured signal to obtain the compressed signal yRN. Then the original signal ˆx is reconstructed with high probability by solving the optimization problem. The sparse representation of a signal x on the orthogonal basis ΨRN×N is

    x=Ψs (1)

    where the transform coefficient s is sparse and contains only k(kN) non-zero elements. We select the measurement matrix Φ that is not related to Ψ. The signal x is projected to the measurement matrix Φ. Then the reduced-dimensional projection data y is expressed as

    y=Φx=As (2)

    where A=ΦΨ is the perception matrix. To carry out effective and unique signal reconstruction, the expression is often solved by optimizing l1 norm [23] and the greedy algorithm

    ˆx=arg minx1 s.t. As=y (3)

    The signal xRN uses a random measurement matrix Φ of size M×N for compression measurement to obtain compressed data y of size M×1. Then the relationship between the original signal and the compressed data can be expressed as

    (1ε)xy2Φ(xy)2(1+ε)xy2 (4)

    where ε(0,1), the original signal x has approximate distance preserving in Euclidean space before and after projection. Thus, sufficient effective information of the original signal is preserved in the compressed data, which provides a theoretical basis for the feature extraction method of the compressed domain in this paper.

    When the measurement matrix meets the requirement of random projection distance preserving, the signal can be accurately reconstructed. In this paper, the Gaussian random matrix that satisfies the random projection distance preserving property with a high probability is used as the measurement matrix to compress and sample the signal. The measurement matrix is constructed by generating a matrix Φ with the size of M×N, so that each element Φ independently obeys a Gaussian distribution with a mean value of 0 and a variance of 1/1MM [24,25]. This measurement matrix is a random measurement matrix, which is not related to most orthogonal bases or orthogonal dictionaries, and the number of measurements required for accurate reconstruction is small.

    The reduction of signal feature information after original signal compression will affect the effect of subsequent fault diagnosis. This paper proposes a week feature enhancement method for improving the accuracy of fault classification.

    Wavelet packet decomposition can decompose the signal into different frequency bands adaptively, without leakage, and without overlapping according to the feature of the signal. It not only improves the time-frequency resolution of the signal but also obtains more detailed information about the signal. In this paper, wavelet packet transform is used to decompose the compressed signal to extract signal features. The time-domain wavelet packet transform matrix of the original signal xRN is TRN×N, then the decomposition process of the original signal x can be expressed as

    f=Tx (5)

    where fRN is the signal component of the signal x after wavelet packet decomposition in time domain space. Similarly, the wavelet packet transform matrix of the compressed signal yRM is ˆTRM×M, then the decomposition process of the compressed signal is

    ˆf=ˆTy (6)

    where ˆfRM is the signal component of the compressed signal y after the wavelet packet decomposition in the compressed space. According to the distance preserving property of random projection, the dimension reduction observation vector ˆf is obtained by projecting the transformation vector f on the random measurement matrix Φ. Similar to the original signal transformation process, the dimensionality reduction observation vector ˆf is regarded as the compressed data y obtained through signal transformation, then the compressed domain signal transformation can be expressed as ˆf=ˆTy. The premise of analyzing the dimensionality reduction observation vector is to obtain the compression domain transformation matrix ˆT, then the solution formula is

    ˆf=Φf=ΦTx=ˆTΦx (7)
    ˆT=ΦTΦ1 (8)

    In Eq (8), the random measurement matrix Φ is not a square matrix, and the least-squares approximation solution needs to be solved by pseudo-inverse Φ1 so as to solve the transformation matrix ˆT in the compression domain. In this paper, the compression domain transformation matrix is firstly used to decompose the compressed signal, and then the signal features are extracted. As shown in Figure 1, the schematic diagram of the compression domain transformation matrix based on random projection distance preservation is presented. Therefore, the method in this paper can extract the characteristic information of the compressed signal for diagnostic analysis on the premise of ensuring that the compressed signal retains enough effective information of the original signal.

    Figure 1.  Schematic diagram of the compressed domain transformation matrix based on random projection distance preservation.

    The compressed signal is decomposed by wavelet packet to obtain signal features of different frequency bands. Since different frequency bands all contain the local feature information of the signal, if this feature information can be fully utilized, the accuracy of fault classification can be further improved. Therefore, this paper proposes to weight the decomposed signal features, and the PSO method is used to optimize the weight coefficient so as to realize the adaptive selection of the signal features of different frequency bands.

    The PSO algorithm is an intelligent optimization algorithm that tracks the optimal particles in the solution space for searching by simulating the predation behavior of birds. Assuming that there are n particles in the D -dimensional search space, where the position of the i -th particle is Xi=(xi1,xi2,,xiD), and the velocity of the corresponding i -th particle is Vi=(vi1,vi2,,viD). The best position experienced by the i-th particle individual is pbesti=(pi1,pi2,,piD), the best position experienced by the population is gbest=(g1,g2,,gD), then the velocity and position update formula of each particle is

    vk+1id=wvkid+c1r1(pbestidxkid)+c2r2(pbestdxkid) (9)
    xk+1id=xkid+vk+1id (10)

    where w is the inertia weight; c1, c2 are acceleration factors; r1, r2 are random numbers between [0, 1]; d=1,2,,D; i=1,2,,D; k is the current iteration frequency.

    In the optimization process of the particle swarm optimization algorithm, each particle has a fitness function to determine whether the position it has experienced so far is optimal, and then the position is updated after comparison. The distribution of signal feature weight coefficients directly affects the fault classification effect. At this time, it is necessary to ensure that the distance between compressed signal feature classes is as large as possible. Therefore, this paper takes the average distance between compressed signal feature classes as the fitness function, and the expression is as follow

    Di=1n1ni,j=1ijdij (11)

    where Di is the average distance between the i -th type and other classes, dij is the distance between the i -th type feature and the j -th type feature, and n is the number of fault types.

    The particle swarm optimization algorithm steps are as follows:

    Step 1: Initialize the particle swarm and set various parameters.

    Step 2: Evaluate the fitness of each particle according to the fitness function.

    Step 3: Compare the fitness value and update the individual historical best position pbest and the global best position gbest.

    Step 4: Update the position x and velocity v of each particle according to Eqs (9) and (10).

    Step 5: Turn to step 2 to loop iteratively until the number of iterations is satisfied or within the allowable error range, and output the best fitness value and particle position.

    The flow chart of the feature enhancement algorithm based on PSO is shown in Figure 2. When setting the parameters of the particle swarm, the main parameter values are shown in Table 1.

    Figure 2.  Flow chart of feature enhancement algorithm based on PSO.
    Table 1.  PSO algorithm parameter settings.
    Allowable Error The Maximum Number of Iterations Particle Swarm Size n Acceleration Factor c1 Acceleration Factor c2 Inertia Weight
    0.0001 400 700 1.5 1.5 0.8

     | Show Table
    DownLoad: CSV

    By combining CS theory with the PSO algorithm, the problem of AE signal data redundancy can be effectively solved, and the frequency band features containing the main feature information can be adaptively selected to avoid the uncertain influence of subjective parameter setting on experimental results. After the original signal is compressed and sampled, the amount of data is reduced. The compressed signal is decomposed by the wavelet packet in the compressed domain. The PSO algorithm is used to adaptively select the obtained frequency band features, and the fault features are enhanced. Finally, the faults are classified by the SVM. The basic process of bearing AE signal feature enhancement fault diagnosis is shown in Figure 3. The specific steps of the method for enhancing the diagnosis of bearing AE signal weak fault features based on CS are as follows:

    Figure 3.  Flow chart of fault diagnosis for bearing AE signal feature enhancement.

    a. AE signals of different fault states of the rolling bearing are collected by the AE acquisition system to form a data sample set.

    b. The Gaussian random matrix is selected as the measurement matrix Φ. Using Eq (2), the N -dimensional AE signal x is projected to obtain the M -dimensional measurement value y.

    c. The wavelet packet decomposition matrix ˆT in the compressed domain is obtained by Eq (7). The transformation matrix ˆT is used for wavelet packet decomposition of the measured value y. The time domain and frequency domain features of different frequency bands are extracted to form the feature set of the compression domain.

    d. The PSO algorithm is used to calculate the weight of signal features of each frequency band so as to obtain feature sets with sufficient distance between classes.

    e. Input the enhanced feature set into the SVM classifier for training and testing, and finally realize the fault diagnosis of the rolling bearing.

    In order to verify the effectiveness of the method in this paper, thrust ball bearings are used for fault simulation experiments, which simulate three states of normal bearing, raceway failure, and rolling element failure. The ball bearing is installed on the bearing failure simulation test bench, which is composed of six parts: drive device, transmission device, test device, loading device, guide rail device, and anti-overload protection device. AE sensors are used to collect bearing AE signals. The AE acquisition system consists of sensors, preamplifiers, data acquisition cards, hosts, and display. The layout of the test bench and sensors is shown in Figure 4. The experiment uses electrical discharge machining technology to process single-point pits on the central position of the bearing raceway and rolling elements to simulate the weak bearing fault. The corresponding fault diameters are 0.5 mm, 1 mm, 1.5 mm, and the depth is 0.65 mm, as shown in Figure 5. The spindle speed of the test bench is 400 rpm, the preamplifier is set to 60 dB, the sampling frequency is 1 MHz, and the number of data points is 1, 638, 400.

    Figure 4.  Bearing failure test bench.
    Figure 5.  Failure diagram of bearing parts.

    This paper mainly studies the AE signal with reduced data volume after compression. The AE signal is obtained by constructing a Gaussian random measurement matrix and performing random projection compression on it. Figure 6 shows the time domain waveform diagram and frequency spectrum of the AE signal before and after compression of a sample under the three states of normal, race failure, and rolling element failure. It can be seen from Figure 6 that compared with the original AE signal, the compressed AE signals in the three states exhibit relatively similar random characteristics in both the time domain and the frequency domain. The obvious periodic fault features of the original signal cannot be well-reflected. This is due to the loss of some time-domain features of the signal during the compression measurement process. At this time, conventional AE signal analysis methods cannot accurately and effectively perform feature extraction and classification.

    Figure 6.  Time-domain waveforms and spectrograms of three states before and after AE signal compression. (a) Original AE signal waveform; (b) Original AE signal spectrum; (c) Compressed measurement data waveform; (d) Compressed measurement data spectrum.

    Since compressed data shows randomness in amplitude, in order to obtain enough feature information of compressed domain signals, this paper extracts features from compressed data in the time domain and frequency domain. This paper extracts twenty-eight fault features including thirteen-time domain parameters (T1T13) and fifteen-frequency domain parameters (T14T28) from the wavelet packet decomposition components of the compressed data. The specific feature parameter names and calculation formulas are shown in Table 2. T19 and T20 indicate the change of the position of the main frequency band and T21T25 indicate the degree of dispersion or concentration of the spectrum. In Table 2, Xk is the spectrum of signal x, where k=1,2,,N/2.56 is the number of spectral lines, and N is the number of sampling points of signal x. PS[k]=2|XK|2/2|XK|2N2N2, f=2.56fsk/2.56fskNN, M=1+N/N22, P=PS[0]+N/2.56k=1PS[k],

    P0=(N/2.56k=0PS[k]f)/P.
    Table 2.  Characteristic parameter table.
    Serial Number Feature Name Serial Number Feature Name
    T1 Absolute mean T15 Spectral mean
    T2 Mean T16 Spectral skewness index
    T3 Peak index T17 Spectral kurtosis index
    T4 Margin index T18 Root mean square frequency
    T5 Kurtosis index T19 P1
    T6 Skewness index T20 P2
    T7 Zero peak T21 P3
    T8 Impulse indicator T22 P4
    T9 Peak-to-peak T23 P5
    T10 Effective value T24 P6
    T11 Waveform index T25 P7
    T12 Variance T26 Frequency center
    T13 Square root amplitude T27 Spectral variance
    T14 Total spectrum T28 Amplitude spectral entropy
    P1=N/2.56k=0(f2PS[k])P, P2=N/2.56k=0(f2PS[k])PN/2.56k=0(f4PS[k]), P3=1/1MMN/2.56k=0(PS[k]P/PMM)3(1/1MMN/2.56k=0(PS[k]P/PMM)2)3, P4=1/1MMN/2.56k=0(PS[k]P/PMM)4(1/1MMN/2.56k=0(PS[k]P/PMM)2)2, P5=2.56N/2.56k=0((fP0)2PS[k])N, P6=P5P0, P7=N/2.56k=0((fp0)3PS[k])P5

     | Show Table
    DownLoad: CSV

    According to the feature enhancement method in this article, firstly, the AE data of different states are obtained. Two hundred sets of data are taken from each state, and 8192 sample points are intercepted from each set of data. Secondly, the AE data is compressed and measured, and a total of 600 sets of data are compressed and projected to obtain compressed data with a length of 4096 (downsampling rate R=2). Then feature extraction is carried out for the compressed domain data. According to experience, three-layer wavelet packet decomposition is selected, and each group of sample signals will generate 8 component signals of different frequency bands so that there are 4800 component signals of 600 × 8. Twenty-eight eigenvalues of each group of sample wavelet packet decomposition data are calculated to obtain an eigenmatrix of the sample signal compression domain with a size of 4800 × 28. Finally, most of the multi-component signals of these compressed data contain low-frequency information and high-frequency noise signals of bearing operation, and the fault information contained therein is masked, which has a strong interference effect on fault identification. Therefore, in order to improve the classification accuracy of compressed data, weaken the interference of high-frequency noise, and enhance the weak fault characteristic signal, the PSO method can be used to optimize the characteristic weight coefficients of the eight components of the wavelet packet decomposition of each sample to realize the characteristic enhancement and obtain a 600 × 28 eigenvalue matrix.

    One feature that is relatively obvious in comparison among the three sample data is selected. As shown in Figure 7, the signal feature is extracted directly after compression and the Characteristic parameters of the signal after the fault feature is enhanced. It can be seen from Figure 7(a) that the compressed feature parameters are in disorder and difficult to classify. In Figure 7(b), the compressed signal is decomposed by wavelet packet by the method of this paper, and the features of each compon-ent signal are extracted for weighted summation. It can be seen that the characteristic parameters of different types of faults obviously have certain laws. Normal and faulty bearings are effectively separated. The two faults of the raceway and rolling element are relatively poor, but they can still be clearly identified. This paper takes different types of fault features as input and selects a suitable proportion of training set and test set. The fault type and classification accuracy are determined by the matching degree of training set prediction classification and test set classification.

    Figure 7.  Comparison of features extracted directly from compressed domain and feature parameters after feature enhancement. (a) Compressed domain feature parameters; (b) Feature parameters after fault feature enhancement.

    The feature set is formed by feature extraction and enhancement of the compressed data, and the feature vector matrix is divided into the training set and the test set by 4:1. The kernel function of SVM selects the Gaussian radial basis function, and the results of classification using this method are shown in Table 3. It can be seen from Table 3 that the data samples are divided into three categories. The test sample numbers 1–40, 41–80, 81–120 are a group of samples corresponding to normal, race failure, and rolling element failure. The 120 samples in the test set use the weak fault feature enhancement algorithm to accurately identify various fault samples with a classification accuracy of 100%.

    Table 3.  Classification results of different types of faults.
    Fault Type Speed/rpm Fault Diameter/mm Sample Number Classification Accuracy Rate/%
    Normal 400 0 1-200 100
    Race failure 400 0.5 201-400 100
    rolling element failure 400 0.5 401-600 100

     | Show Table
    DownLoad: CSV

    In order to further verify the effectiveness of the method in this paper, the different faults of bearing races and rolling elements are diagnosed. When the spindle speed is 400 rpm, the diameter of the race failure and the rolling element failure are respectively 0.5 mm, 1 mm and 1.5 mm. There are three types of failure samples. The specific parameter settings are shown in Table 4. After the same compression measurement, three-layer wavelet packet decomposition, feature enhancement, and classification process, the ratio of the training set to the test set is 4:1. The experimental characteristic parameters of different degrees of failure before and after the compression domain feature enhancement are shown in Figure 8. It can be seen from Figure 8(a) and (c) that feature extraction is conducted on the compressed data obtained from the race fault degree experiment and the rolling element fault degree experiment, and the characteristic parameters of different fault types are all mixed, irregular and difficult to separate. From Figure 8(b) and (d), it can be seen that the feature parameters of different fault levels obtained by the feature enhancement method of the PSO algorithm to optimize the weight coefficients show more regularity, better classification and recognition, and larger class spacing. So as to provide more accurate classification information for subsequent bearing fault diagnosis and greatly improve the accuracy of diagnosis. The classification results of the fault feature enhancement diagnosis method are shown in Table 5. The weak fault feature enhancement method based on the proposed method can accurately identify faults of different degrees with 100% classification accuracy for different experimental types. Therefore, the four different degrees of race faults and rolling element faults can be accurately identified with 100% classification accuracy for fault diagnosis using this method.

    Table 4.  Different degrees of seat ring/rolling element failure parameter table.
    Fault State Speed/rpm Fault Diameter/mm Sample Number
    Race/rolling element failure degree 400 0 1–200
    400 0.5 201–400
    400 1 401–600
    400 1.5 601–800

     | Show Table
    DownLoad: CSV
    Figure 8.  Compressed domain features of different experimental types enhance the front and back feature parameters. (a) Feature parameters directly extracted from race faults; (b) Feature parameters after enhanced race fault features; (c) Feature parameters directly extracted from rolling element faults; (d) Feature parameters after enhanced rolling element fault features.
    Table 5.  Classification results of different degrees of seat/rolling element failure.
    Experiment Type Training Samples Test Sample Identify the Sample Diagnostic Accuracy Rate/%
    Race failure degree test 640 160 160 100
    Rolling element failure degree experiment 640 160 160 100

     | Show Table
    DownLoad: CSV

    Through the verification of experimental data, we can see that the fault diagnosis method in this paper has obtained perfect results. It can not only correctly classify various faults but also identify different degrees of faults, and the recognition accuracy of each fault can reach 100%. It shows that the method proposed in this paper is effective, which can accurately diagnose different types of bearing faults while reducing the amount of data.

    The compression rate reflects the degree of compression of the original signal. Restricted by the constrained equidistance condition, an excessively large compression rate under the condition of a limited number of sampling points will cause serious information loss of the compressed signal. The signal reconstruction error of the normal bearing AE signal under different compression ratios is shown in Figure 9. It can be seen that as the compression rate increases, the reconstruction error of the signal gradually increases, and the reconstruction error is basically in the range of 30 to 50%. When the compression ratio is in the range of 2–4, the OMP method can obtain a smaller reconstruction error, and the reconstruction effect has become better. When the compression ratio is greater than 10, the reconstruction error is large, and the feature information in the compressed data is greatly reduced.

    Figure 9.  Signal reconstruction errors under different compression rates.

    Discrete Cosine Transform (DCT) is an orthogonal projection transformation method that transforms a signal from the time domain to the frequency domain. The key information of the signal is mainly concentrated on a few low-frequency coefficients. DCT is suitable for sparse representation of local singular parts and has a good effect on sparse representation of low-frequency parts. It has strong energy concentration characteristics and decorrelation performance and has good robustness to noise. The Orthogonal Match Pursuit (OMP) algorithm is an iterative greedy algorithm and a classic algorithm for solving the l1 norm. The algorithm calculates the current optimal solution in each iteration, updates the residuals, and finds the optimal global solution through continuous iterations. The calculation speed is fast, and it is easy to implement. It is very suitable for signal reconstruction.

    According to the above analysis, the DCT method is used to sparsely represent the original AE signal, the OMP method is used to reconstruct the signal, and the compression rate R = 2 is selected to reconstruct the original signal. The signal reconstruction error is 0.31, and the signal-to-noise ratio is 5, indicating that the reconstruction method using CS can restore the original signal more accurately. Therefore, the effective key information in the original signal is retained in the compressed measurement data, and the original signal can be reconstructed with a lower reconstruction error. This provides sufficient feature information for feature extraction in the compressed domain.

    In order to further verify the effect of the bearing AE signal fault feature enhancement method based on CS, we selected several bearing fault diagnosis methods to compare with the method in this paper. Mainly include: a) Input the compressed data directly into the neural network and use the Softmax classifier to classify; b) Extract the features from the compressed data and input the support vector machine; c) The original AE signal input the support vector machine after extracting the features. Under the conditions of the same parameters, the classification accuracy of the four bearing AE signal fault diagnosis methods is shown in Figure 10.

    Figure 10.  Classification results of different methods under three experimental types.

    It can be seen from Figure 10 that four-fault diagnosis methods are used for different degrees of seat ring failure, different degrees of rolling element failure, and different types of failures. The method in this paper has high diagnostic efficiency in the three types of experiments. The accuracy rate reaches 100%, which can realize accurate diagnosis of different fault states. The accuracy of the compressed data input neural network method is low under the three experimental conditions; the method of direct classification of compressed domain feature extraction has the lowest diagnostic accuracy, and the lowest is only 57.5%; the diagnostic accuracy of the method of extracting features from the original AE signals for classification is relatively high, up to 96.67%, but fluctuates greatly under different experiments. In summary, compared with the traditional methods, the method in this paper has obvious advantages by enhancing the weak fault characteristics after CS. No matter in terms of diagnostic accuracy and the wide applicability and stability of the diagnostic method, the method in this paper shows good results.

    The reason why the method in this paper is superior to other methods is that the processing object of this paper is to compress the signal, and the less signal contains most of the characteristic information of the original time-domain signal. Then the weak fault features of the compressed signal can be enhanced by some adaptive selection of frequency band information, and the classification accuracy can be higher with fewer data; however, the traditional signal classification method based on time-frequency statistical feature parameters cannot accurately obtain all the information of the signal, so the diagnosis accuracy rate is low when combined with the SVM classifier; in addition, the compression signal classification method based on neural network is limited by the number of hidden layers, the number of training samples and the training method, resulting in unstable network performance and low diagnostic recognition rate.

    Aiming at the problems of high cost and difficulty in data processing caused by a large amount of data transmission in the fault diagnosis of bearing AE signal, this paper proposes an enhanced diagnosis method for weak fault features of bearing AE signal based on CS. The method presented in this paper is successfully applied to the fault diagnosis of rolling bearings under the condition of ensuring high diagnostic accuracy. The method directly deals with the AE signal in the compression domain of the rolling bearing fault. The wavelet packet decomposition and the PSO method are used to realize the intelligent selection of frequency bands with more fault feature information. It not only solves the problem of the large data volume of AE signal data processing in fault diagnosis but also solves the problem of weak fault characteristics of compressed data. The signal appears in a random and disordered state in the compressed domain. The components decomposed by wavelet packet contain a large amount of feature information. The method of adaptive feature weighting is used to extract effective and enhanced fault information, which avoids the blindness of component selection. The effectiveness of the proposed fault feature enhancement diagnosis method is verified by bearing fault experiments. The method in this paper can accurately identify different types of faults and different degrees of faults. The appropriate compression ratio is selected through the reconstruction error curve of the compression ratio so as to prove that the compression domain signal has enough feature information, and ensure the effectiveness of feature enhancement and the accuracy of fault diagnosis. Finally, compared with the other three traditional fault diagnosis methods, the method in this paper has better performance in terms of the stability of diagnosis and the accuracy of recognition. The method proposed in this paper provides a new idea for the fault diagnosis of compressed signal processing. At the same time, it provides a new diagnosis method for equipment fault diagnosis under big modern data, which has a high engineering application value.

    This work was supported in part by the National Key R & D Program of China under Awards 2018YFB1306100 and in part by the National Natural Science Foundation of China under Awards 51875272 and in part by Kunming University of Science and Technology Scientific Research Start-up Fund for Talents Introduction under Awards KK23201801048, and in part by Yunnan Provincial School Education Cooperation Key Project under Awards KKDA202001003.

    The authors declare no conflict of interest.



    [1] Y. Wang, M. Lehtomäki, X. Liang, J. Pyörälä, A. Kukko, A. Jaakkola, et al., Is field-measured tree height as reliable as believed a comparison study of tree height estimates from field measurement, airborne laser scanning and terrestrial laser scanning in a boreal forest, ISPRS J. Photogramm., 147 (2019), 132-145. doi: 10.1016/j.isprsjprs.2018.11.008
    [2] I.S. Saliu, B. Satyanarayana, M.A. B. Fisol, G. Wolswijk, C. Decannière, R. Lucas, et al., An accuracy analysis of mangrove tree height mensuration using forestry techniques, hypsometers and UAVs, Estuar. Coast. Shelf Sci., 248 (2020), 106971.
    [3] D.W. Wanik, J.R. Parent, E.N. Anagnostou, B.M. Hartman, Using vegetation management and LiDAR-derived tree height data to improve outage predictions for electric utilities, Electr. Pow. Syst. Res., 146 (2017), 236-245. doi: 10.1016/j.epsr.2017.01.039
    [4] Y. D. Huang, M. Z. Li, S. Q. Ren, M. J. Wang, P. Y. Cui, Impacts of tree-planting pattern and trunk height on the airflow and pollutant dispersion inside a street canyon, Build. Environ., 165 (2019), 106385. doi: 10.1016/j.buildenv.2019.106385
    [5] Y. Xu, C. Li, Z. Sun, L. Jiang, J. Fang, Tree height explains stand volume of closed-canopy stands: evidence from forest inventory data of China, Forest Ecol. Manag., 438 (2019), 51-56. doi: 10.1016/j.foreco.2019.01.054
    [6] C. Cabo, C. Ordóñez, C. A. López-Sánchez, J. Armesto, Automatic dendrometry: tree detection, tree height and diameter estimation using terrestrial laser scanning, Int. J. Appl. Earth Obs., 69 (2018), 164-174. doi: 10.1016/j.jag.2018.01.011
    [7] J. S. Cui, J. Huo, M. Yang, Y. K. Wang, Research on the rigid body posture measurement using monocular vision by coplanar feature points, Optics, 126 (2015), 5423-5429.
    [8] H. Du, M. G. Li, The study for particle image velocimetry system based on binocular vision, Measurement, 42 (2009), 619-627. doi: 10.1016/j.measurement.2008.10.011
    [9] X. Hu, H. Zheng, Y. Chen, L. Chen, Dense crowd counting based on perspective weight model using a fisheye camera, Optik, 126 (2015), 123-130. doi: 10.1016/j.ijleo.2014.08.132
    [10] H. Kim, J. Jung, L. Paik, Fisheye lens camera based surveillance system for wide field of view monitoring, Optik, 127 (2016), 5636-5646. doi: 10.1016/j.ijleo.2016.03.069
    [11] J. Zhu, J. Zhu, X. Wan, C. Wu, C. Xu, Object detection and localization in 3D environment by fusing raw fisheye image and attitude data, J. Vis. Commun. Image R., 59 (2019), 128-139. doi: 10.1016/j.jvcir.2019.01.005
    [12] W. C. Wang, C. H. Hwang, C. I. Chu, Y. H. Chen, Displacement measurement of interior wall of hollow cylinder by digital image correlation method using fisheye lens, Procedia Eng., 79 (2014), 437-446. doi: 10.1016/j.proeng.2014.06.365
    [13] Z. Y. Liu, F. Ding, Y. Xu, X. Han, Background dominant colors extraction method based on color image quick fuzzy c-means clustering algorithm, Def. Technol., (2020).
    [14] J. C. Dunn, A fuzzy relative of the ISODATA process and its use in detecting compact well-separated clusters, Taylor Francis, 3 (1973), 32-57.
    [15] Y. H. Chen, L. J. Zhang, H. T. Lang, Acceleration of fuzzy c-means clustering based target detection for large size image, J. Ocean U. China, 47 (2017), 94-100.
    [16] S. Sandhya, B. Chidambararajan, M. S. Kumar, FCM-based segmentation and neural network classification of tumor in brain mri images, in Intelligence in Big Data Technologies-Beyond the Hype, Springer, (2021), 371-378.
    [17] B. J. Shankar, K. Murugan, A. Obulesu, D. Finney, R. Anitha, MRI image segmentation using bat optimization algorithm with fuzzy c means (BOA-FCM) clustering, J. Med. Imag. Health Inform., 11 (2021), 661-666. doi: 10.1166/jmihi.2021.3365
    [18] M. Li, Y. Wang, Q. Sun, Y. Liu, Research of ASW-FCM-based algorithm for clustered wind turbine group equivalent modeling, J. Electr. Eng. Technol., 15 (2020), 1555-1566. doi: 10.1007/s42835-020-00439-0
    [19] J. Fan, J. Wang, A two-phase fuzzy clustering algorithm based on neurodynamic optimization with its application for polsar image segmentation, IEEE Trans. Fuzzy Syst., 26 (2016), 72-83.
    [20] R. Shang, P. Tian, L. Jiao, R. Stolkin, J. Feng, B. Hou, et al., A spatial fuzzy clustering algorithm with kernel metric based on immune clone for SAR image segmentation, IEEE J. Stars., 9 (2016), 1640-1652.
    [21] Y. Zhang, X. Bai, R. Fan, Z. Wang, Deviation-sparse fuzzy c-means with neighbor information constraint, IEEE Trans. Fuzzy Syst., 27 (2019), 185-199. doi: 10.1109/TFUZZ.2018.2883033
    [22] S. Krinidis, V. Chatzis, A robust fuzzy local information C-means clustering algorithm, IEEE Trans. Image Process., 19 (2010), 1328-1337. doi: 10.1109/TIP.2010.2040763
    [23] X. Bai, Z. Chen, Y. Zhang, Z. Liu, Y. Lu, Infrared ship target segmentation based on spatial information improved FCM, IEEE Trans. Cyber., 46 (2015), 3259-3271.
    [24] H. Zhang, Q. Wang, W. Shi, M. Hao, A novel adaptive fuzzy local information C-Means clustering algorithm for remotely sensed imagery classification, IEEE Trans. Geosci. Remote Sen., 55 (2017), 5057-5068. doi: 10.1109/TGRS.2017.2702061
    [25] P. K. Mishro, S. Agrawal, L. Dora, R. Panda, A fuzzy C-means clustering approach to HMRF-EM model for MRI brain tissue segmentation, in 2017 6th International Conference on Computer Applications in Electrical Engineering-Recent Advances (CERA), Springer, (2017), 371-376.
    [26] H. Zhang, Q. Wu, Y. Zheng, T. M. Nguyen, D. Wang, Effective fuzzy clustering algorithm with Bayesian model and mean template for image segmentation, IET Image Process, 8 (2014), 571-581. doi: 10.1049/iet-ipr.2013.0178
    [27] L. Szilagyi, Z. Benyo, S. M. Szilagyi, H. S. Adam, MR brain image segmentation using an enhanced fuzzy C-means algorithm, in Proceedings of International Conference of the IEEE Engineering in Medicine & Biology Society, Springer, (2003), 724-726.
    [28] D. Comaniciu, P. Meer, Mean shift: a robust approach toward feature space analysis, IEEE Transactions Pattern Anal. Mach. Intell., 24 (2002), 603-619. doi: 10.1109/34.1000236
    [29] A. Y. Ng, M. I. Jordan, Y. Weiss, On spectral clustering: analysis and an algorithm, in Advances in Neural Information Processing Systems, Springer, (2002), 849-856.
    [30] T. Lei, X. Jia, Y. Zhang, S. Liu, H. Meng, A.K. Nandi, Superpixel-based fast fuzzy c-means clustering for color image segmentation, IEEE Trans. Fuzzy Syst., (2018), 1753-1766.
    [31] P. Arbeláez, M. Maire, C. Fowlkes, J. Malik, Contour detection and hierarchical image segmentation, IEEE Transactions Pattern Anal. Mach. Intell., 33 (2011), 898-916. doi: 10.1109/TPAMI.2010.161
    [32] Q. Lin, Q. Zhu, P. Huang, J. Chen, Z. Ming, J. Yu, Computers, A novel hybrid multi-objective immune algorithm with adaptive differential evolution, Comput. Oper. Res., 62 (2015), 95-111. doi: 10.1016/j.cor.2015.04.003
    [33] L. Zuo, C. Luo, Y. Zuo, Paralleled segmentation cluster algorithm based on En FCM for large-scale image, Microcomput Its Appl., 34 (2015), 55-58.
    [34] J. Song, Z. Zhang, A modified robust fcm model with spatial constraints for brain MR image segmentation, Information, 10 (2019), 74. doi: 10.3390/info10020074
    [35] H. X. Pei, Z. R. Zheng, C. Wang, C. N. Li, Y. H. Yao, D-fcm: density based fuzzy C-means clustering algorithm with application in medical image segmentation, Procedia Comput. Sci., 122 (2017), 407-414. doi: 10.1016/j.procs.2017.11.387
  • This article has been cited by:

    1. Ariane Verdy, Priyanga Amarasekare, Alternative stable states in communities with intraguild predation, 2010, 262, 00225193, 116, 10.1016/j.jtbi.2009.09.011
    2. Vardit Makler-Pick, Matthew R. Hipsey, Tamar Zohary, Yohay Carmel, Gideon Gal, Intraguild Predation Dynamics in a Lake Ecosystem Based on a Coupled Hydrodynamic-Ecological Model: The Example of Lake Kinneret (Israel), 2017, 6, 2079-7737, 22, 10.3390/biology6020022
    3. Yuanshi Wang, Hong Wu, Shikun Wang, Wen Shi, Dynamics of Intraguild Predation Systems with Intraspecific Competition, 2018, 80, 0092-8240, 2408, 10.1007/s11538-018-0467-6
    4. Richard C. Bruce, Intraguild Interactions and Population Regulation in Plethodontid Salamanders, 2008, 22, 0733-1347, 31, 10.1655/07-015.1
    5. K. Blue Pahl, David J. Yurkowski, Kirsty J. Lees, Nigel E. Hussey, Measuring the occurrence and strength of intraguild predation in modern food webs, 2020, 25, 23522496, e00165, 10.1016/j.fooweb.2020.e00165
    6. Melanie J. Hatcher, Jaimie T. A. Dick, Alison M. Dunn, How parasites affect interactions between competitors and predators, 2006, 9, 1461-023X, 1253, 10.1111/j.1461-0248.2006.00964.x
    7. Yuanshi Wang, Hong Wu, Population dynamics of intraguild predation in a lattice gas system, 2015, 259, 00255564, 1, 10.1016/j.mbs.2014.11.001
    8. Zhiguang Liu, Fengpan Zhang, Species coexistence of communities with intraguild predation: The role of refuges used by the resource and the intraguild prey, 2013, 114, 03032647, 25, 10.1016/j.biosystems.2013.07.010
    9. Yun Kang, Lauren Wedekin, Dynamics of a intraguild predation model with generalist or specialist predator, 2013, 67, 0303-6812, 1227, 10.1007/s00285-012-0584-z
    10. Yuanshi Wang, Donald L. DeAngelis, Stability of an intraguild predation system with mutual predation, 2016, 33, 10075704, 141, 10.1016/j.cnsns.2015.09.004
    11. Gaku Takimoto, Takeshi Miki, Maiko Kagami, Intraguild predation promotes complex alternative states along a productivity gradient, 2007, 72, 00405809, 264, 10.1016/j.tpb.2007.04.005
    12. Richard J. Hall, Intraguild predation in the presence of a shared natural enemy, 2011, 92, 0012-9658, 352, 10.1890/09-2314.1
    13. Kejun Zhuang, Hongjun Yuan, Spatiotemporal dynamics of a predator–prey system with prey-taxis and intraguild predation, 2019, 2019, 1687-1847, 10.1186/s13662-019-1945-3
    14. Manuel Falconi, Yrina Vera-Damián, Claudio Vidal, Predator interference in a Leslie–Gower intraguild predation model, 2020, 51, 14681218, 102974, 10.1016/j.nonrwa.2019.102974
    15. Yiding Yang, Zhilan Feng, Dashun Xu, Gregory J. Sandland, Dennis J. Minchella, Evolution of host resistance to parasite infection in the snail–schistosome–human system, 2012, 65, 0303-6812, 201, 10.1007/s00285-011-0457-x
    16. Suparna Dash, Subhas Khajanchi, Dynamics of intraguild predation with intraspecies competition, 2023, 1598-5865, 10.1007/s12190-023-01956-7
    17. Emanuel A. Fronhofer, Lynn Govaert, Mary I. O'Connor, Sebastian J. Schreiber, Florian Altermatt, The shape of density dependence and the relationship between population growth, intraspecific competition and equilibrium population density, 2023, 0030-1299, 10.1111/oik.09824
    18. Purnedu Mishra, Dariusz Wrzosek, Schoener-Polis-Holt's model of the intraguild predation with predator taxis and repulsive chemotaxis, 2024, 0, 1531-3492, 0, 10.3934/dcdsb.2024062
    19. Debjit Pal, Dipak Kesh, Debasis Mukherjee, Study of a diffusive intraguild predation model, 2025, 2, 2997-6006, 10.20935/AcadEnvSci7508
    20. Jiang Li, Xianning Liu, Yangjiang Wei, Multistable switches induced by prudent predation in three-species food web models with omnivory, 2025, 198, 09600779, 116577, 10.1016/j.chaos.2025.116577
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3222) PDF downloads(71) Cited by(3)

Figures and Tables

Figures(18)  /  Tables(7)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog