# A New Cooperative Anomaly Detection Method for Stacker Running Track of Automated Storage and Retrieval System in Industrial Environment

Authors: Darong Huang; Mengting Lin; Lanyan Ke; Zhenping Deng
Journal: Journal of Control Science and Engineering
 (2018)

---

## Abstract

Considering the complexity and the criticality of the stacker equipment, in order to solve the problem that the stop accuracy of the stacker reduces or even fails to work due to abrasion of the running rail, this paper proposes a cooperative detection method based on Pulse Coupling Neural Network (PCNN) and wavelet transform theory to detect the abnormal points of the stacker running rail in industrial environment by analyzing the variation signals. First of all, considering the fact that the data is mixed up with noises because of the environment at the site and the possibility of the data acquisition equipment breaking down, a noise reduction method for the vibration signal data of stacker is constructed based on PCNN. Then, the basic theory of wavelet transform is introduced and then the rules of judging anomaly points on stackers’ running tracks are discussed based on wavelet transform. In addition, a cooperative detection method based on PCNN and wavelet transform theory is carried out based on the space-time distribution feature of the vibration of the stacker orbits in the industrial environment. Then the rationality of the proposed algorithm is verified by simulation through data provided by State Grid Measuring Center of China. This paper constructs a model of the abnormal point detection of the stackers in an industrial environment. The experimental simulation and example simulation show that the cooperative detection method based on PCNN and wavelet transform theory can effectively detect and locate the anomaly points of the stacker running tracks. The expansibility in engineering applications is promising. Lastly, some conclusions are discussed.


---

## Body

## 1. Introduction

As is known to all, the stacker is a key equipment of automated storage and retrieval system (ASRS) in industrial environment. In practice scene of ASRS, the main function of the stacker is to grab, move, and stack goods from one shelf to another. Thus, the stability of stacker running track will influence the accuracy of grabbing and moving goods in practical engineering. As one of the key components of the stacker equipment, the track of the stacker is divided into upper and lower track. In fact, the stacking machine would wear, crack, sag, and bulge by long time running. And these defects will reduce the accuracy of the stacking machine. If the field engineer cannot detect and repair the defect as soon as possible, they may even cause the whole equipment to be worn and shut down. So, how to detect the anomaly point of stacker running rail is very important for safety performance of the whole system.In anomaly detection, domestic and foreign researchers have made a lot of progress in all aspects. Since the 1980s, the anomaly detection problem has been widely researched in the field of statistics. For examples, [1] presented an anomaly detection algorithm based on statistical theory. Unfortunately, practice has proved that the computing results of the statistical method are not ideal, because the distribution characteristics of the data must be known in advance. Obviously, the shortage limits the application scope of the method. To solve the defect, some improved detection algorithms were put forward combining with distance in [2, 3]. The obvious characteristic of the improved algorithm is required to have a lot of domain knowledge of real systems or equipment. However, it is known to us all that the running performance of real systems or equipment is affected by various factors. Therefore, it is difficult to determine the input parameters of the algorithm. In other words, it has some disadvantages such as instability and high algorithm complexity. Thus, once the testing data is incomplete, it is difficult to determine the degree of abnormal [4]. In particular, when the data has features such as high-dimensional and sparseness, the performance of the algorithm is very poor. In order to solve the problem caused by sparse data, Rastogi and Ramaswamy proposed an improved algorithm based on density in [5]. The method can avoid the limitation of the dependence of distance to a certain context. To better handle testing data to find the abnormal data point, some detecting ideas based on depth continue to appear. For instance, some researchers presented the anomaly detection algorithm based on depth to mark each record as a point in the k dimension space. In fact, the typical DEEPLOC algorithm was proposed by Struyf and Rousseeuw in [6]. One big advantage of the algorithm is that the depth value of each point may be given according to the definition of depth. The detection rule is that a record with a smaller depth is more likely to be an anomaly point than a larger one. Thus, once the data set is organized by the depth value, the algorithm only needs to carry out the outlier detection on the lower layer of the depth value. Meanwhile, the outlier detection is not required in the record on a layer with a large depth value.Furthermore, Zhou S and Xu W have constructed the local anomaly detection algorithm based on the deviation in [7]. The algorithm and rule may solve the problem that the existing local anomaly detection algorithms do not perform data object partition. But in real running process of the systems, the high dimensional test data would lead to increased computational complexity [8]. So, the trick to find out the abnormal point from high-dimension data set is to reduce its dimensionality. Aiming at the high dimensionality of the data set, some researchers introduced the concept of local projection scoring (LPS) and proposed an efficient abnormal point detection method of high-dimensional data in [9]. Of course, some new algorithms (for instance, neural network [10] and clustering algorithm [11]) continue to be introduced and presented to detect and monitor the abnormal points of test data. The simulation results show that these models and methods are effective in real engineering and application.In addition, some scholars had discussed the anomaly detection in the frequency domain. For instance, to overcome the disadvantages that the window size does not change with frequency, the wavelet transform theory is usually introduced to compensate for the localization defects of the short-time Fourier transform [12]. The characteristics of multiresolution (also called multiscale) can be used to observe signal gradually from coarse to fine. The detection and observational modes are favorable for detecting singularities of the signal step by step. So, the wavelet theorem is also an ideal signal time-frequency analysis and processing tools. Among them, the most representative algorithm is the time series anomaly detection method based on wavelet transform proposed in [12].Although the engineers and scholars have made a lot of achievements in the field of the detection of anomaly spot, the identifying and detection of abnormal data is still in its infancy of research for the stacker of ASRS. The practice application results show that there are still some problems in the methods of safety supervisory of the stacker of ASRS. As is known to all, the running performance of stacker is affected by various factors such as running environment, structural characteristics, and the optimal goal of the whole system. Hence, how to implement the detection of the abnormal points is the kernel problem of safety maintenance of stacker in ASRS. At present, there are some research achievements about the detection of a stacker’s running performance. But the research is mainly focused on the structure analysis or the design of the system, and few people have done their researches from the perspective of abnormal point’s detection [13]. In recent years, many new algorithms are presented and proposed based on Internet, OPC, and the fault tree to analyze the data set acquired and measured from the running state of the stacker in [14]. In engineering application, the detecting ability of the remote fault of stacker may meet the desired purpose using the proposed algorithm presented in [15, 16]. Meanwhile, Kai Zhang and his coauthors have discussed and analyzed the monitoring method based on the multimode and multivariate statistics to monitor the running state of stacker crane in [17, 18]. Experimental results show that the anomaly detection method is effective.However, because the running state of stacker is influenced by the state of the stacker’s component, the scene environment, the data acquisition equipment, and so on, the data set acquired from the system will contain a lot of noise. Obviously, the noise of data will reduce the accuracy of abnormal point’s judgment. Therefore, deleting and cleaning the noise from testing data is necessary to implement and accomplish the monitoring of running state of whole stacker system. To restrain the interferences of strong background noise, Huang D.R. and his coauthors have constructed a cooperated denoising algorithm for rolling bearing of stacker in [19]. The simulation results have verified the effectiveness of health monitoring of ASRS. Notice that, in the previous scenario, the actual monitoring data obtained from the real stacker is the aliasing vibration signal. So, how to accomplish the separation of multivibration signal is a difficult problem in vibration process. If the method is reasonable, the detecting accuracy of stacker’s running abnormal points will be greatly improved.Based on the analysis and the thesis above, to ensure the effectiveness of the incomplete data processing of real system, it is necessary to construct and design a cooperative anomaly detection algorithm so that the abnormal spot can be detected and located as quickly as possible in on-site industrial environment. Notice that the timing of the real time data process is vital in the industrial environment, and then the data signal needs to be treated from coarse to fine as soon as possible. On the basis, the Pulse Coupling Neural Network (PCNN) presented in [20, 21] is introduced to denoise, because it has the advantage that the industrial data process does not depend on precise mathematics model. Meanwhile, due to the fact that the wavelet transform may complete the itemization of the data collection, it is introduced to construct and establish the abnormal point detecting algorithm to locate the stacker’s defect through the pure data processed by PCNN.Hence, the rest of this paper will discuss the details of the algorithm and thesis. The layout of the rest of the paper is organized as follows: Section2 will introduce the basic concepts of PCNN and the modified PCNN is introduced to construct the data denoising model. In Section 3, in order to locate the anomaly points, the wavelet transform is introduced as the anomaly detection and location algorithm. Also, the cooperative anomaly judgment algorithm and rule for anomaly detection will be discussed in detail combined with modified PCNN and wavelet transform. In addition, the algorithm flow chart is also drawn. Later, the contrast experiments and numerical simulation to detect the abnormal point of stacker running track are carried out to verify the effect of the algorithm and rule using data provided by State Grid Measuring Center of China. Finally, some conclusions and the directions for future engineering application are discussed according to the real stacker running track of ASRS in industrial environment.


## 2. Improved PCNN Denoising Model and Algorithm for Vibration Signals of Stacker Running Track

In a real working condition, the actual vibration signals measured from this system will be unavoidably affected by many complicated environmental factors. Obviously, the data package usually includes strong noises. So, to guarantee the effectiveness of anomaly detection for stacker running track, a reasonable data preprocessing procedure is very crucial to eliminate the noises that are contained in the dataset. In this context, constructing an effective denoising model and algorithm to process the original signal is of great theoretical and practical significance for the condition monitoring of the stacker running track.However, most industrial monitoring and control applications require high performance, timeliness, and reliability. Then, most administrators and engineers hope to effectively operate the system without knowing the accurate model. Based on this thesis, the PCNN will be introduced later. On this basis, an improved PCNN denoising model and algorithm are analyzed and designed according to the actual situation to ensure the timeliness and stability of the performance of the stacker running track.


### 2.1. Basic Theory of PCNN

As we all know, PCNN is presented by Eckhorn based on the observed synchronous pulse transmission after the experiments of the cerebral cortex of the animals [22]. In fact, due to its scale invariance, rotation invariance, intensity invariance, distortion invariant, and other characteristics, PCNN is widely used in image smoothing, image segmentation, image edge detection, image fusion, optimal solution, and so on. Moreover, in theory, the PCNN model has similar group neurons synchronization release pulse characteristic and the accurate model is not needed to parse the structure of dataset. Thus, the engineers not only reveal the inherent ability of PCNN, but also explore the application of PCNN in the signal denoising. So, the basic concept of PCNN will be introduced in the next context.According to [22], the PCNN model can be expressed by the following equation:(1)Fijn=e-τFijFijn-1+VF∑k,lMijklYkln-1+Sij(2)Lijn=e-τLijLijn-1+VL∑k,lWijklYkln-1(3)Uijn=Fijn1+βLijn(4)θijn=e-τθijθijn-1+VθYijn-1(5)Yijn=1Uijn>θijn0Uijn≤θijnwhere Fij[n] denotes the nth feedback input for the (i,j) neuron; Sij, Uij, and Yij represent the external stimulation, internal behavior, and output of neuron Nij, respectively; Lij, Fij are the two input channels for the link domain and the feedback domain of the neuron Nij; M and W are the connection weight coefficient matrix of feedback domain and link domain; θij and Vθ are the output and threshold amplification factor and the variable threshold function; β denotes synaptic link coefficient; τL, τF, τθ are the time attenuation constants of link domain, feedback domain, and variable threshold function, respectively.From the perspective of simulation, the PCNN neuron consists of three parts: receiving domain, modulation domain, and pulse generation domain. In real application, PCNN has the advantage that the data processing does not depend on precise mathematics model. That is to say, in the pretreatment of denoising, once the network interface of PCNN receives the input signal, the receiving field transmits it through two channels,F and L. The impulse response function of F channel changes slowly with time compared to L channel. The modulation part combines the signal Lij from channel L with an offset and multiplies the signal Fij from channel F to generate the internal signal Uij. Then, θij and Uij will be compared to control the firing of signal neurons. If θij>Uij, the neurons will be activated. Otherwise, the neurons may be deleted from the structure of PCNN.According to the link coefficient ofL channel and F channel, PCNN can be divided into two cases: coupled and uncoupled. When β=0, each neuron was separately operated and unaffected by the pulse output of other neurons around it. In addition, considering the friendly interface for end users, the basic structure of PCNN is shown as Figure 1.Figure 1
Basic structure of PCNN model diagram.Obviously, the original dataset may be used to identify the character of noises by training the PCNN network. Of course, the noise-polluted data may be cleaned through analyzing the PCNN firing matrix and then the noises may also be filtered from the original dataset. In general, the engineers may test the abnormal point in a long data sequence to better locate the defect of stacker running track.Unfortunately, although the basic structure of PCNN can clean the real dataset mixed with noise, the complexity of this topological structure can cause various issues that negatively impacts the engineer’s operation play experience. In real scenario of ASRS, the difficulty in condition monitoring is to ensure the timeliness of locating the abnormal position of stacker running track. In particular, because the running process of stacker is complex and changeable, the original dataset signified the running state of the stacker including lots of noise. From an operational perspective, the traditional denoising method based on nominal model may hardly achieve the expected denoising performance. So, it is necessary to simplify network structure to meet the actual demand and it is a nontrivial problem to construct an improved denoising algorithm based on the existing PCNN model to meet the actual needs.


### 2.2. Modified PCNN Denoising Model for Vibration Signals of Stacker Running Track

In the industrial scenario, due to very poor measurement environment, the quality of the measurement signals may be affected by all kinds of factors. Meanwhile, the influential factors are connected with each other in vibration signals of stacker running track and they also have indirect effects. How to solve the coupling relations within the limits of real time control is very important to denoise the original dataset. In fact, considering that the basic topological structure of PCNN has the excessive parameters and the implementation is more complicated, it is not widely used in practice compared to the modified PCNN model. For a handy operation or statistics, the modified uncoupled model proposed in [21] is adopted to simplify the denoising process.According to [21], the simplified PCNN model can be expressed by the following equation:(6)Fijn=Sijn(7)Yijn=1Fijn>θijn0Fijn≤θijn(8)θijn=VθYijn=1e-τθθijn-1Yijn=0(9)Bijn=1Yijn=10Yijn=0(10)Tijn=nYijn=1Tijn-1Yijn=0whereFij[n] denotes the nth feedback input for the (i,j) neuron; Sij denotes the external stimulation of neurons Nij; τθ denotes time attenuation constants of variable threshold function.In this model, the signal from the channelL is not accepted by Fij. As a result, Fij is only affected by external input. The B matrix records whether the neuron fires at the nth input. And when the neuron fires, it will be assigned a larger threshold, so that it will not be fired again in a short time. Conversely, the threshold of the neuron that is not fired will decay over time. Therefore, the T matrix is used to record the time when each neuron is firstly fired. The noise point in the data can be filtered out by the value of T matrix.The modified topological structure of simplified PCNN is shown as Figure2.Figure 2
Modified PCNN model diagram.Obviously, the modified model has advantages of less parameters and simple implementation. In other words, the simplified denoising system has operational simplicity and high efficiency for the denoising processing of the raw incoming data. For the sake of applying in industry, the simplified denoising process is more applicable. Thus, we can make use of the modified denoising model to denoise the vibration signals of stacker running track and finally realize the recognition of anomaly points.


## 2.1. Basic Theory of PCNN

As we all know, PCNN is presented by Eckhorn based on the observed synchronous pulse transmission after the experiments of the cerebral cortex of the animals [22]. In fact, due to its scale invariance, rotation invariance, intensity invariance, distortion invariant, and other characteristics, PCNN is widely used in image smoothing, image segmentation, image edge detection, image fusion, optimal solution, and so on. Moreover, in theory, the PCNN model has similar group neurons synchronization release pulse characteristic and the accurate model is not needed to parse the structure of dataset. Thus, the engineers not only reveal the inherent ability of PCNN, but also explore the application of PCNN in the signal denoising. So, the basic concept of PCNN will be introduced in the next context.According to [22], the PCNN model can be expressed by the following equation:(1)Fijn=e-τFijFijn-1+VF∑k,lMijklYkln-1+Sij(2)Lijn=e-τLijLijn-1+VL∑k,lWijklYkln-1(3)Uijn=Fijn1+βLijn(4)θijn=e-τθijθijn-1+VθYijn-1(5)Yijn=1Uijn>θijn0Uijn≤θijnwhere Fij[n] denotes the nth feedback input for the (i,j) neuron; Sij, Uij, and Yij represent the external stimulation, internal behavior, and output of neuron Nij, respectively; Lij, Fij are the two input channels for the link domain and the feedback domain of the neuron Nij; M and W are the connection weight coefficient matrix of feedback domain and link domain; θij and Vθ are the output and threshold amplification factor and the variable threshold function; β denotes synaptic link coefficient; τL, τF, τθ are the time attenuation constants of link domain, feedback domain, and variable threshold function, respectively.From the perspective of simulation, the PCNN neuron consists of three parts: receiving domain, modulation domain, and pulse generation domain. In real application, PCNN has the advantage that the data processing does not depend on precise mathematics model. That is to say, in the pretreatment of denoising, once the network interface of PCNN receives the input signal, the receiving field transmits it through two channels,F and L. The impulse response function of F channel changes slowly with time compared to L channel. The modulation part combines the signal Lij from channel L with an offset and multiplies the signal Fij from channel F to generate the internal signal Uij. Then, θij and Uij will be compared to control the firing of signal neurons. If θij>Uij, the neurons will be activated. Otherwise, the neurons may be deleted from the structure of PCNN.According to the link coefficient ofL channel and F channel, PCNN can be divided into two cases: coupled and uncoupled. When β=0, each neuron was separately operated and unaffected by the pulse output of other neurons around it. In addition, considering the friendly interface for end users, the basic structure of PCNN is shown as Figure 1.Figure 1
Basic structure of PCNN model diagram.Obviously, the original dataset may be used to identify the character of noises by training the PCNN network. Of course, the noise-polluted data may be cleaned through analyzing the PCNN firing matrix and then the noises may also be filtered from the original dataset. In general, the engineers may test the abnormal point in a long data sequence to better locate the defect of stacker running track.Unfortunately, although the basic structure of PCNN can clean the real dataset mixed with noise, the complexity of this topological structure can cause various issues that negatively impacts the engineer’s operation play experience. In real scenario of ASRS, the difficulty in condition monitoring is to ensure the timeliness of locating the abnormal position of stacker running track. In particular, because the running process of stacker is complex and changeable, the original dataset signified the running state of the stacker including lots of noise. From an operational perspective, the traditional denoising method based on nominal model may hardly achieve the expected denoising performance. So, it is necessary to simplify network structure to meet the actual demand and it is a nontrivial problem to construct an improved denoising algorithm based on the existing PCNN model to meet the actual needs.


## 2.2. Modified PCNN Denoising Model for Vibration Signals of Stacker Running Track

In the industrial scenario, due to very poor measurement environment, the quality of the measurement signals may be affected by all kinds of factors. Meanwhile, the influential factors are connected with each other in vibration signals of stacker running track and they also have indirect effects. How to solve the coupling relations within the limits of real time control is very important to denoise the original dataset. In fact, considering that the basic topological structure of PCNN has the excessive parameters and the implementation is more complicated, it is not widely used in practice compared to the modified PCNN model. For a handy operation or statistics, the modified uncoupled model proposed in [21] is adopted to simplify the denoising process.According to [21], the simplified PCNN model can be expressed by the following equation:(6)Fijn=Sijn(7)Yijn=1Fijn>θijn0Fijn≤θijn(8)θijn=VθYijn=1e-τθθijn-1Yijn=0(9)Bijn=1Yijn=10Yijn=0(10)Tijn=nYijn=1Tijn-1Yijn=0whereFij[n] denotes the nth feedback input for the (i,j) neuron; Sij denotes the external stimulation of neurons Nij; τθ denotes time attenuation constants of variable threshold function.In this model, the signal from the channelL is not accepted by Fij. As a result, Fij is only affected by external input. The B matrix records whether the neuron fires at the nth input. And when the neuron fires, it will be assigned a larger threshold, so that it will not be fired again in a short time. Conversely, the threshold of the neuron that is not fired will decay over time. Therefore, the T matrix is used to record the time when each neuron is firstly fired. The noise point in the data can be filtered out by the value of T matrix.The modified topological structure of simplified PCNN is shown as Figure2.Figure 2
Modified PCNN model diagram.Obviously, the modified model has advantages of less parameters and simple implementation. In other words, the simplified denoising system has operational simplicity and high efficiency for the denoising processing of the raw incoming data. For the sake of applying in industry, the simplified denoising process is more applicable. Thus, we can make use of the modified denoising model to denoise the vibration signals of stacker running track and finally realize the recognition of anomaly points.


## 3. Design and Analysis of Cooperative Anomaly Detection Method and Algorithm for Stacker Running Track of ASRS

In the above section, one valid dimensional data cleaning model and framework was discussed and analyzed. Once the original dataset measured from the track running system was processed by the modified PCNN, the errors, noises, or missing data that are contained in the dataset should be eliminated and removed. That is to say, the output dataset may be used to accurately locate the anomalous point for the stacker running track of ASRS. Obviously, the key stage of locating the abnormal defect of stacker running track is to find a reasonable model. The engineering practice shows that the energy of induction signal induced by the abnormal point of the stacker running track concentrates on the high frequency band. Of course, the fact has also provided a basis for determining anomaly point of stacker running stacker. Based on this, it is imperative to find a reasonable way to decompose the output dataset of modified PCNN into two parts: low frequency part and high frequency part.Simultaneously, because the wavelet transform developing since 80’s of last century has specific property of time-frequency localization, the method is especially fit for analyzing and detecting local signal mutation for the stacker running track of ASRS. Next step, the cooperative anomaly detection model of stacker running track will be constructed and analyzed based on wavelet transform.


### 3.1. Cooperative Anomaly Detection Method Based Wavelet Transform

In our experiments, the main feature of signal mutation is that the signal has local changes in time and space. To guarantee the accuracy of locating and identifying the defects of stacker running track, we have tried to design and construct a cooperative detection method combined with wavelet theory. The principle of using wavelet transform to detect the anomalous points is to decompose the signals in different resolutions. When the signal is abrupt, the coefficient gotten by the wavelet transform has a modulus maximum value. Therefore, the location of the outliers can be filtered out through the detection of the maximum modulus point.Without loss of generality, the basic principle of wavelet transform is clearly defined: Supposeϕ(t)∈L2(R) is a square integral function, if its Fourier transform satisfies the following condition.(11)∫RΨω2ω-1dω<+∞whereψ(ω) is wavelet transform of ϕt and ϕt is a wavelet function or base wavelet.In fact, when the output dataset of modified PCNN was decomposed into high frequency and low frequency, the key factor is to carry the expansion or contraction, and the translation based on the base wavelet. If the expansion or contraction factor is supposed asa and the translation factor is supposed asb, then the wavelet function can be rewritten as follows.(12)ϕa,bt=a-1/2ϕt-baa>0,b∈Rwherea denotes the scaling factor (or scale factor) and b is translation factor. And ϕa,b(t) is wavelet base function which is concerned with a and b. It can be obtained by scaling or moving the parent function ϕt.So, for an arbitrary signalx(t)∈L2(R), the wavelet transform can be expressed by the inner product of the signal and wavelet base function (13)Wxa,b=x,ϕa,bt=a-1/2∫Rxtϕa,bt-badtNotice that the actual measured signals indicating the changes of the running state of the stacker running track are discrete in the stacker running track in ASRS, and the wavelet transform should be rewritten in a discrete form to treat the output dataset provided by modified PCNN. In general, the scaling factora and the translation factor b are put as follows.(14)a=2-j,b=2-jkThus, Wx(a,b) became a discrete binary wavelet function as follows.(15)ϕj,k=2j/2ϕj,k2jt-k(16)Wxj,k=2j/2∫Rxtϕ2-jt-kdtwhere j,k∈Z.Then, the output dataset denoising by modified PCNN may be decomposed into high frequency and low frequency by the following formula.(17)Wxj,k=∑j,k=-∞∞xtϕj,ktObviously, the wavelet series information may be used to show the changes of the running state of the stacker track. Moreover, to ensure that the energy of induction signal induced by the abnormal point of the stacker running track concentrates on the high frequency band, we have selected the orthogonal wavelet transform to decompose the output dataset into high frequency and low frequency. The basic decomposing mode is defined as follows.(18)Cj+1=H∗CjDj+1=D∗CjwhereC is the high frequency part and D is the low frequency part.In our experiments, if the raw output signalx(t) is described as C0, the wavelet decomposing process may be shown as Figure 3.Figure 3
The wavelet decomposing process of output dataset.Obviously, according to the decomposing process in Figure3, the high frequency band of the dataset, which indicates the energy of induction signal induced by the abnormal point of the stacker running track, can be gotten. To add extra flexibility, the high frequency band is described as Cj(j=0,1,⋯).In practical project, the low-frequency coefficients reflect the contour of the original signal and the high-frequency coefficients describe the details of the signal. In particular, the singularity in the signals is often caused by a sudden change in frequency domain. For the engineers in health monitoring of stacker running state, it means that the high-frequency coefficients of wavelet transform can highlight singularity and can be used for detection and localization of the defect of the stacker running track. From the angle of engineering application, the damage of stacker running track will cause perturbations of measured signal at damage sites. Moreover, the measurement dataset is an aggregate of the running information, and the singularity of high frequency may depict the character of damage position on stacker running track. In other words, the damage concussion will cause the signal saltation. So, in practical engineering applications, it is very important to design a reasonable judgment rule of anomaly points.To better locate the defects of stacker running track, the Lipschitz index is used to describe and design the judgment rule of singularity of the dataset. For the simplicity of analysis, the corresponding concept of Lipschitz index was defined as follows.Definition 1.
If there exists a constantc>0 andh, such that f(t) satisfies the following condition(19)ft0+h-ft0≤chαthus, α is the Lipschitz index of f(t) at t=t0. And then if α<1, f(t) is singular at t=t0.In general, the size of Lipschitz index is related to the value of the singularity in actual project; i.e., the more severe the degree of mutation, the steeper the peak of the catastrophic point, and the smaller the singularity index, and vice versa. Based on this thesis, we can define the local singularity of high frequency signals decomposed by wavelet transform as follows.Definition 2.
∀ x ∈ δ x 0, if there exists a constant K>0 ands, the wavelet function ϕ(x) is continuous and differentiable and the wavelet transform has N vanishing moment (N is a positive integer); i.e., (20)Wfs,x≤KsαThus, α is the local singularity at x0.Notice that wavelet transform is applied to vibration signal analysis of stacker running track to detect the meshing abnormality of track with local defects, and the amplitudes of the vibration single will decrease or increase to some extent. So, by Definition1, once there exists α<1 at one point, the engineers may judge that the position of the stacker running track may be damaged. Thus, the judgment rule of anomaly point is designed as follows.Rule 1.
∀ x ∈ δ x 0, if the wavelet transform coefficient satisfies (21)Wfs,x≤Wfx,x0thus, x0 may be judged as the abnormal point of the stacker running stack.Notice that the mutation information of flaws may be accumulated in running process of stacker’s track, and the ability of a single high-frequency detail signal to reflect an abnormal point has limitations to locate the defect to stacker running track. To solve and overcome the problem, the product of the high-frequency coefficient, which can amplify the detail signal, is selected as the basis for the final judgment of the abnormal point in our experiments. To satisfy the need of engineering design, the accumulated information of track’s flaws can be described by using the following formula.(22)C=Cj×Cjn-1×Cjn-2where C is the product of high-frequency signal. Cj, Cjn-1, Cjn-2 are the n-order, n-1-order, and n-2-order high-frequency detail signals at the decomposition scalen, respectively.Then, the basic rule is shown as follows.Rule 2.
If a mutation appears in the product of the high-frequency coefficient, it is an abnormal point.Through the analysis and rules above, once the output dataset processed by modified PCNN has been decomposed into the low frequency and high frequency components of the signals, the high frequency components can be used to detect and locate the anomaly points of stacker running track.


### 3.2. Design and Analysis of Cooperative Anomaly Detection Algorithm Based on Modified PCNN and Wavelet Transform

Based on the above analysis and discussion, combining with the modified PCNN and wavelet transform, the cooperative anomaly detection algorithm of the stacker running track may be designed in detail as below.Step 1.
Normalize the data. Calculate the initial mean square error and initialize the network parameters by (1)-(9).Step 2.
Use the normalized data as input data. End the loop if matrixB=0. Get matrix T when the loop ends.Step 3.
The noise points are determined according to the elements in the matrix, and then filter each data point.Step 4.
Calculate mean square error (MSE), and compare it with the mean square error before, if the mean square error is smaller than the mean square error before, back to Step2. Input data into that model to continue to reduce noise. Otherwise, output the signal last time as the result.Step 5.
Select wavelet base function and determine the number of layers of wavelet decomposition.Step 6.
The wavelet transform coefficients of each layer are obtained by wavelet decomposition.Step 7.
Multiply the detail signal to amplify the mutation signal.Step 8.
Obtain the location of the outliers according to the location of the signal mutation.The algorithm flow chart is shown as Figure4.Figure 4
Flow chart of the algorithm.


## 3.1. Cooperative Anomaly Detection Method Based Wavelet Transform

In our experiments, the main feature of signal mutation is that the signal has local changes in time and space. To guarantee the accuracy of locating and identifying the defects of stacker running track, we have tried to design and construct a cooperative detection method combined with wavelet theory. The principle of using wavelet transform to detect the anomalous points is to decompose the signals in different resolutions. When the signal is abrupt, the coefficient gotten by the wavelet transform has a modulus maximum value. Therefore, the location of the outliers can be filtered out through the detection of the maximum modulus point.Without loss of generality, the basic principle of wavelet transform is clearly defined: Supposeϕ(t)∈L2(R) is a square integral function, if its Fourier transform satisfies the following condition.(11)∫RΨω2ω-1dω<+∞whereψ(ω) is wavelet transform of ϕt and ϕt is a wavelet function or base wavelet.In fact, when the output dataset of modified PCNN was decomposed into high frequency and low frequency, the key factor is to carry the expansion or contraction, and the translation based on the base wavelet. If the expansion or contraction factor is supposed asa and the translation factor is supposed asb, then the wavelet function can be rewritten as follows.(12)ϕa,bt=a-1/2ϕt-baa>0,b∈Rwherea denotes the scaling factor (or scale factor) and b is translation factor. And ϕa,b(t) is wavelet base function which is concerned with a and b. It can be obtained by scaling or moving the parent function ϕt.So, for an arbitrary signalx(t)∈L2(R), the wavelet transform can be expressed by the inner product of the signal and wavelet base function (13)Wxa,b=x,ϕa,bt=a-1/2∫Rxtϕa,bt-badtNotice that the actual measured signals indicating the changes of the running state of the stacker running track are discrete in the stacker running track in ASRS, and the wavelet transform should be rewritten in a discrete form to treat the output dataset provided by modified PCNN. In general, the scaling factora and the translation factor b are put as follows.(14)a=2-j,b=2-jkThus, Wx(a,b) became a discrete binary wavelet function as follows.(15)ϕj,k=2j/2ϕj,k2jt-k(16)Wxj,k=2j/2∫Rxtϕ2-jt-kdtwhere j,k∈Z.Then, the output dataset denoising by modified PCNN may be decomposed into high frequency and low frequency by the following formula.(17)Wxj,k=∑j,k=-∞∞xtϕj,ktObviously, the wavelet series information may be used to show the changes of the running state of the stacker track. Moreover, to ensure that the energy of induction signal induced by the abnormal point of the stacker running track concentrates on the high frequency band, we have selected the orthogonal wavelet transform to decompose the output dataset into high frequency and low frequency. The basic decomposing mode is defined as follows.(18)Cj+1=H∗CjDj+1=D∗CjwhereC is the high frequency part and D is the low frequency part.In our experiments, if the raw output signalx(t) is described as C0, the wavelet decomposing process may be shown as Figure 3.Figure 3
The wavelet decomposing process of output dataset.Obviously, according to the decomposing process in Figure3, the high frequency band of the dataset, which indicates the energy of induction signal induced by the abnormal point of the stacker running track, can be gotten. To add extra flexibility, the high frequency band is described as Cj(j=0,1,⋯).In practical project, the low-frequency coefficients reflect the contour of the original signal and the high-frequency coefficients describe the details of the signal. In particular, the singularity in the signals is often caused by a sudden change in frequency domain. For the engineers in health monitoring of stacker running state, it means that the high-frequency coefficients of wavelet transform can highlight singularity and can be used for detection and localization of the defect of the stacker running track. From the angle of engineering application, the damage of stacker running track will cause perturbations of measured signal at damage sites. Moreover, the measurement dataset is an aggregate of the running information, and the singularity of high frequency may depict the character of damage position on stacker running track. In other words, the damage concussion will cause the signal saltation. So, in practical engineering applications, it is very important to design a reasonable judgment rule of anomaly points.To better locate the defects of stacker running track, the Lipschitz index is used to describe and design the judgment rule of singularity of the dataset. For the simplicity of analysis, the corresponding concept of Lipschitz index was defined as follows.Definition 1.
If there exists a constantc>0 andh, such that f(t) satisfies the following condition(19)ft0+h-ft0≤chαthus, α is the Lipschitz index of f(t) at t=t0. And then if α<1, f(t) is singular at t=t0.In general, the size of Lipschitz index is related to the value of the singularity in actual project; i.e., the more severe the degree of mutation, the steeper the peak of the catastrophic point, and the smaller the singularity index, and vice versa. Based on this thesis, we can define the local singularity of high frequency signals decomposed by wavelet transform as follows.Definition 2.
∀ x ∈ δ x 0, if there exists a constant K>0 ands, the wavelet function ϕ(x) is continuous and differentiable and the wavelet transform has N vanishing moment (N is a positive integer); i.e., (20)Wfs,x≤KsαThus, α is the local singularity at x0.Notice that wavelet transform is applied to vibration signal analysis of stacker running track to detect the meshing abnormality of track with local defects, and the amplitudes of the vibration single will decrease or increase to some extent. So, by Definition1, once there exists α<1 at one point, the engineers may judge that the position of the stacker running track may be damaged. Thus, the judgment rule of anomaly point is designed as follows.Rule 1.
∀ x ∈ δ x 0, if the wavelet transform coefficient satisfies (21)Wfs,x≤Wfx,x0thus, x0 may be judged as the abnormal point of the stacker running stack.Notice that the mutation information of flaws may be accumulated in running process of stacker’s track, and the ability of a single high-frequency detail signal to reflect an abnormal point has limitations to locate the defect to stacker running track. To solve and overcome the problem, the product of the high-frequency coefficient, which can amplify the detail signal, is selected as the basis for the final judgment of the abnormal point in our experiments. To satisfy the need of engineering design, the accumulated information of track’s flaws can be described by using the following formula.(22)C=Cj×Cjn-1×Cjn-2where C is the product of high-frequency signal. Cj, Cjn-1, Cjn-2 are the n-order, n-1-order, and n-2-order high-frequency detail signals at the decomposition scalen, respectively.Then, the basic rule is shown as follows.Rule 2.
If a mutation appears in the product of the high-frequency coefficient, it is an abnormal point.Through the analysis and rules above, once the output dataset processed by modified PCNN has been decomposed into the low frequency and high frequency components of the signals, the high frequency components can be used to detect and locate the anomaly points of stacker running track.


## 3.2. Design and Analysis of Cooperative Anomaly Detection Algorithm Based on Modified PCNN and Wavelet Transform

Based on the above analysis and discussion, combining with the modified PCNN and wavelet transform, the cooperative anomaly detection algorithm of the stacker running track may be designed in detail as below.Step 1.
Normalize the data. Calculate the initial mean square error and initialize the network parameters by (1)-(9).Step 2.
Use the normalized data as input data. End the loop if matrixB=0. Get matrix T when the loop ends.Step 3.
The noise points are determined according to the elements in the matrix, and then filter each data point.Step 4.
Calculate mean square error (MSE), and compare it with the mean square error before, if the mean square error is smaller than the mean square error before, back to Step2. Input data into that model to continue to reduce noise. Otherwise, output the signal last time as the result.Step 5.
Select wavelet base function and determine the number of layers of wavelet decomposition.Step 6.
The wavelet transform coefficients of each layer are obtained by wavelet decomposition.Step 7.
Multiply the detail signal to amplify the mutation signal.Step 8.
Obtain the location of the outliers according to the location of the signal mutation.The algorithm flow chart is shown as Figure4.Figure 4
Flow chart of the algorithm.


## 4. Examples and Analysis

### 4.1. Simulation Examples

#### 4.1.1. Simulation Results of Noise Reduction Effect of PCNN

To verify the noise reduction effect of modified PCNN, the simulation examples were first used to test the denoising ability of PCNN. In our simulation experiments, the sinusoidal signal was selected and adopted as the input data to train and test PCNN rule and the Gaussian white noise with a mean of zero and a variance of 1 is added to the sinusoidal signal as the experimental signal. The testing signalx(t) was described as follows.(23)xt=sin⁡t+vwherev is the Gaussian white noise, and v~N(0,1).In actual imitative preparation,v,n,and t were set up as below, respectively.(24)vt=normrnd0,0.1,nt=0:0.1:100n=lengthsin⁡tBy training the neural network of PCNN, the denoising results are shown as Figure 5.Figure 5
PCNN noise reduction results.Simultaneously, to explain the validity of PCNN denoising algorithm, the wavelet threshold denoising algorithm was used to compare with the presented algorithm in our simulation examples. The effect of wavelet threshold denoising algorithm is shown as Figure6.Figure 6
Wavelet threshold denoising results.Obviously, although the comparison results between Figures5 and 6 indicate the differences of the different denoising algorithm, the details of the differences are not obvious and significant. So, we must have a valid index selected to show the differences. To go through this process, the signal-to-noise ratio (SNR) and minimum mean square error (MSE) were used to compare the advantages of proposed denoising algorithm. The basic formulas were shown as follows.(1) Signal-to-Noise Ratio (25) S N R = S ^ t v   (2) Mean Square Error (26) M S E = 1 n ∑ t = 1 n S t - S ^ t 2whereS^t is denosing signal, St is the original signal, and v is the noise signal.By (25) and (26), the corresponding comparison results are shown in Table 1.Table 1
Comparison of wavelet threshold denoising and PCNN noise reduction.
Wavelet denoising PCNN denoising SNR 15.15 16.47 MSE 0.37 0.36It can be seen in Table1 that PCNN can retain the characteristics of the original signal while getting good noise reduction with higher SNR and lower MSE. So, it means that the denoising effect of modified PCNN is better than wavelet threshold denoising method.


#### 4.1.2. Simulation Results of Anomaly Detection

In order to further verify the validity of cooperative anomaly detection method for singular point in signal dataset, a set of vibration signal with abnormal points is selected to test the effect of presented cooperative anomaly detection algorithm; i.e.,(27)xt=sin⁡0.02πt0≤t≤400sin⁡0.2πt401≤t≤600In our simulation experiments, the simulation results are shown as Figure7.Figure 7
Simulation result of cooperative anomaly extraction.
(a)
The original signal containing the exception point (b)
Wavelet decomposition result (c)
Abnormal point detection resultsIn fact, it can be seen from the simulation results that the sudden change of the signal appears at t=400 which can be observed from detail coefficients clearly. That indicates that the anomalies appear at t=400. It can be concluded that the presented cooperative anomaly detection can mark the location of the abnormal point of the vibration signal clearly. On the other hand, after testing the function, the presented cooperative algorithm may be applied to actual project.


### 4.2. Real Experiment Examples

Moreover, to further verify the validity of the proposed algorithm, the actual stacking machine vibration signal of ASRS, which has been developed and applied in State Grid Measuring Center of China, was selected to simulate the performance of the anomaly detection algorithm. The prototype systems of the stackers of State Grid Measuring Center of China is shown in Figure8.Figure 8
Stackers in State Grid Measuring Center of China.As can be seen in Figure8, the main function is to grab, move, and stack goods from one goods shift to another. In actual project, the degree of unevenness of flaws (such as rusts, crack, the concavo-convex fluctuation, or any welding defects) will directly influence the safety and effectiveness of whole ASRS. In real engineering, the test rig of stacker running track is shown as Figure 9.Figure 9
Basic structure of stacker running track.Figures8 and 9 demonstrate that the crucial equipment of ASRS is the stacker running track. And then the operation track systems of the stacker contain the suspension transporting system (i.e., upper track) and road tracks (i.e., lower tracks). Thus, if outliers are presented and detected within the dataset, the defects of the stacker running track may be located and identified according to the precise geospatial data.In real application, the sampling time is from March to April of 2018. The size of the sampling data source measured from stacker running track is 7370. The data may contain noise because of the interference of the site environment and collecting data equipment and other reasons. To the Management Department of State Grid Measuring Center of China, the engineers especially want to find the anomaly points and locate the defects of stacker running track from the experimental data. If the engineers can do this, the whole security of ASRS may be better ensured. So, the simulation results between the original signal and compounded signal with noise caused by the defects of stacker’s running track were simulated in ASRS. In order to verify the effectiveness of the proposed model, wavelet transform and box-plot are chosen to compare with the result of PCNN-wavelet model. The results were shown as Figure10.Figure 10
Comparison of abnormal point detection of three methods.
(a)
Data anomaly detection results based on wavelet (b)
Data anomaly detection results based on box-plot (c)
Results of data anomaly detection based on PCNN-waveletAs seen from Figure10(a), due to the interference of noise points, it is impossible to accurately judge the abnormal points of the stacker by using only the wavelet transform. As can be seen from Figure 10(b), the abnormal point determination interval is 1 m-3.6m, 5.4m-6.4m, 11.8m-13.2m, 14.8m-16.1m, and 22.3m-23.2m. Correspondingly, the anomaly extensions may be at 1.5m-2.0m, 6.7m-7.6m, 10 m, 17.8m, 20 m, and 22.4m by the judgement rules in Figure 10(c). The results can be shown as Table 2.Table 2
Comparison of abnormal point detection results.
Wavelet Box-plot PCNN-wavelet ∖ 1m-3.6m 1.5m-2.0m ∖ 5.4m-6.4m 6.7m-7.6m ∖ 11.8m-13.2m 10m ∖ 14.8m-16.1m 17.8m ∖ ∖ 20m ∖ 22.3m-23.2m 22.4mObviously, after the noise reduction of the input signal, the location of the abrupt change in the vibration signal may be depicted and described accurately by PCNN-wavelet. To ensure the operability in real engineering practice, some actual damage scenes of stacker running track were used to verify the effect of the presented algorithm. For instance, the actual testing and detecting sample photos at 1.5m, 7.2m, and 18 m are shown in Figure11, respectively.Figure 11
Testing and detecting sample photos from scene.
(a)
Depression of lower track at 1.5m (b)
Cracks appearing on upper track at 7.2m (c)
Interface deformation on upper track at 18 mIn addition, there are rail welds at 10 m, 20 m, and 22 m, respectively.From the experimental simulation results and the actual situation verification, although the box diagram can roughly locate the abnormal point interval, the PCNN-wavelet is more accurate and closer to the actual situation. In Figure10(c), the simulative curve appears to be changed suddenly in these places. In other words, the cooperative anomaly detection method presented in this paper is effective.In fact, the emergence of data mutation points may be caused by the depressions, cracks, and deformation at the interface in the tracks. On the other hand, the response information of the damage position may be developed to make different damage positions in different information range. And the administrators of the industrial system may locate damage positions according to numerical characteristics in the damage unit. Furtherly, the relevant departments may design the good maintenance strategies according to the monitoring curves. That means our algorithm may be applied to the real project.


## 4.1. Simulation Examples

### 4.1.1. Simulation Results of Noise Reduction Effect of PCNN

To verify the noise reduction effect of modified PCNN, the simulation examples were first used to test the denoising ability of PCNN. In our simulation experiments, the sinusoidal signal was selected and adopted as the input data to train and test PCNN rule and the Gaussian white noise with a mean of zero and a variance of 1 is added to the sinusoidal signal as the experimental signal. The testing signalx(t) was described as follows.(23)xt=sin⁡t+vwherev is the Gaussian white noise, and v~N(0,1).In actual imitative preparation,v,n,and t were set up as below, respectively.(24)vt=normrnd0,0.1,nt=0:0.1:100n=lengthsin⁡tBy training the neural network of PCNN, the denoising results are shown as Figure 5.Figure 5
PCNN noise reduction results.Simultaneously, to explain the validity of PCNN denoising algorithm, the wavelet threshold denoising algorithm was used to compare with the presented algorithm in our simulation examples. The effect of wavelet threshold denoising algorithm is shown as Figure6.Figure 6
Wavelet threshold denoising results.Obviously, although the comparison results between Figures5 and 6 indicate the differences of the different denoising algorithm, the details of the differences are not obvious and significant. So, we must have a valid index selected to show the differences. To go through this process, the signal-to-noise ratio (SNR) and minimum mean square error (MSE) were used to compare the advantages of proposed denoising algorithm. The basic formulas were shown as follows.(1) Signal-to-Noise Ratio (25) S N R = S ^ t v   (2) Mean Square Error (26) M S E = 1 n ∑ t = 1 n S t - S ^ t 2whereS^t is denosing signal, St is the original signal, and v is the noise signal.By (25) and (26), the corresponding comparison results are shown in Table 1.Table 1
Comparison of wavelet threshold denoising and PCNN noise reduction.
Wavelet denoising PCNN denoising SNR 15.15 16.47 MSE 0.37 0.36It can be seen in Table1 that PCNN can retain the characteristics of the original signal while getting good noise reduction with higher SNR and lower MSE. So, it means that the denoising effect of modified PCNN is better than wavelet threshold denoising method.


### 4.1.2. Simulation Results of Anomaly Detection

In order to further verify the validity of cooperative anomaly detection method for singular point in signal dataset, a set of vibration signal with abnormal points is selected to test the effect of presented cooperative anomaly detection algorithm; i.e.,(27)xt=sin⁡0.02πt0≤t≤400sin⁡0.2πt401≤t≤600In our simulation experiments, the simulation results are shown as Figure7.Figure 7
Simulation result of cooperative anomaly extraction.
(a)
The original signal containing the exception point (b)
Wavelet decomposition result (c)
Abnormal point detection resultsIn fact, it can be seen from the simulation results that the sudden change of the signal appears at t=400 which can be observed from detail coefficients clearly. That indicates that the anomalies appear at t=400. It can be concluded that the presented cooperative anomaly detection can mark the location of the abnormal point of the vibration signal clearly. On the other hand, after testing the function, the presented cooperative algorithm may be applied to actual project.


## 4.1.1. Simulation Results of Noise Reduction Effect of PCNN

To verify the noise reduction effect of modified PCNN, the simulation examples were first used to test the denoising ability of PCNN. In our simulation experiments, the sinusoidal signal was selected and adopted as the input data to train and test PCNN rule and the Gaussian white noise with a mean of zero and a variance of 1 is added to the sinusoidal signal as the experimental signal. The testing signalx(t) was described as follows.(23)xt=sin⁡t+vwherev is the Gaussian white noise, and v~N(0,1).In actual imitative preparation,v,n,and t were set up as below, respectively.(24)vt=normrnd0,0.1,nt=0:0.1:100n=lengthsin⁡tBy training the neural network of PCNN, the denoising results are shown as Figure 5.Figure 5
PCNN noise reduction results.Simultaneously, to explain the validity of PCNN denoising algorithm, the wavelet threshold denoising algorithm was used to compare with the presented algorithm in our simulation examples. The effect of wavelet threshold denoising algorithm is shown as Figure6.Figure 6
Wavelet threshold denoising results.Obviously, although the comparison results between Figures5 and 6 indicate the differences of the different denoising algorithm, the details of the differences are not obvious and significant. So, we must have a valid index selected to show the differences. To go through this process, the signal-to-noise ratio (SNR) and minimum mean square error (MSE) were used to compare the advantages of proposed denoising algorithm. The basic formulas were shown as follows.(1) Signal-to-Noise Ratio (25) S N R = S ^ t v   (2) Mean Square Error (26) M S E = 1 n ∑ t = 1 n S t - S ^ t 2whereS^t is denosing signal, St is the original signal, and v is the noise signal.By (25) and (26), the corresponding comparison results are shown in Table 1.Table 1
Comparison of wavelet threshold denoising and PCNN noise reduction.
Wavelet denoising PCNN denoising SNR 15.15 16.47 MSE 0.37 0.36It can be seen in Table1 that PCNN can retain the characteristics of the original signal while getting good noise reduction with higher SNR and lower MSE. So, it means that the denoising effect of modified PCNN is better than wavelet threshold denoising method.


## 4.1.2. Simulation Results of Anomaly Detection

In order to further verify the validity of cooperative anomaly detection method for singular point in signal dataset, a set of vibration signal with abnormal points is selected to test the effect of presented cooperative anomaly detection algorithm; i.e.,(27)xt=sin⁡0.02πt0≤t≤400sin⁡0.2πt401≤t≤600In our simulation experiments, the simulation results are shown as Figure7.Figure 7
Simulation result of cooperative anomaly extraction.
(a)
The original signal containing the exception point (b)
Wavelet decomposition result (c)
Abnormal point detection resultsIn fact, it can be seen from the simulation results that the sudden change of the signal appears at t=400 which can be observed from detail coefficients clearly. That indicates that the anomalies appear at t=400. It can be concluded that the presented cooperative anomaly detection can mark the location of the abnormal point of the vibration signal clearly. On the other hand, after testing the function, the presented cooperative algorithm may be applied to actual project.


## 4.2. Real Experiment Examples

Moreover, to further verify the validity of the proposed algorithm, the actual stacking machine vibration signal of ASRS, which has been developed and applied in State Grid Measuring Center of China, was selected to simulate the performance of the anomaly detection algorithm. The prototype systems of the stackers of State Grid Measuring Center of China is shown in Figure8.Figure 8
Stackers in State Grid Measuring Center of China.As can be seen in Figure8, the main function is to grab, move, and stack goods from one goods shift to another. In actual project, the degree of unevenness of flaws (such as rusts, crack, the concavo-convex fluctuation, or any welding defects) will directly influence the safety and effectiveness of whole ASRS. In real engineering, the test rig of stacker running track is shown as Figure 9.Figure 9
Basic structure of stacker running track.Figures8 and 9 demonstrate that the crucial equipment of ASRS is the stacker running track. And then the operation track systems of the stacker contain the suspension transporting system (i.e., upper track) and road tracks (i.e., lower tracks). Thus, if outliers are presented and detected within the dataset, the defects of the stacker running track may be located and identified according to the precise geospatial data.In real application, the sampling time is from March to April of 2018. The size of the sampling data source measured from stacker running track is 7370. The data may contain noise because of the interference of the site environment and collecting data equipment and other reasons. To the Management Department of State Grid Measuring Center of China, the engineers especially want to find the anomaly points and locate the defects of stacker running track from the experimental data. If the engineers can do this, the whole security of ASRS may be better ensured. So, the simulation results between the original signal and compounded signal with noise caused by the defects of stacker’s running track were simulated in ASRS. In order to verify the effectiveness of the proposed model, wavelet transform and box-plot are chosen to compare with the result of PCNN-wavelet model. The results were shown as Figure10.Figure 10
Comparison of abnormal point detection of three methods.
(a)
Data anomaly detection results based on wavelet (b)
Data anomaly detection results based on box-plot (c)
Results of data anomaly detection based on PCNN-waveletAs seen from Figure10(a), due to the interference of noise points, it is impossible to accurately judge the abnormal points of the stacker by using only the wavelet transform. As can be seen from Figure 10(b), the abnormal point determination interval is 1 m-3.6m, 5.4m-6.4m, 11.8m-13.2m, 14.8m-16.1m, and 22.3m-23.2m. Correspondingly, the anomaly extensions may be at 1.5m-2.0m, 6.7m-7.6m, 10 m, 17.8m, 20 m, and 22.4m by the judgement rules in Figure 10(c). The results can be shown as Table 2.Table 2
Comparison of abnormal point detection results.
Wavelet Box-plot PCNN-wavelet ∖ 1m-3.6m 1.5m-2.0m ∖ 5.4m-6.4m 6.7m-7.6m ∖ 11.8m-13.2m 10m ∖ 14.8m-16.1m 17.8m ∖ ∖ 20m ∖ 22.3m-23.2m 22.4mObviously, after the noise reduction of the input signal, the location of the abrupt change in the vibration signal may be depicted and described accurately by PCNN-wavelet. To ensure the operability in real engineering practice, some actual damage scenes of stacker running track were used to verify the effect of the presented algorithm. For instance, the actual testing and detecting sample photos at 1.5m, 7.2m, and 18 m are shown in Figure11, respectively.Figure 11
Testing and detecting sample photos from scene.
(a)
Depression of lower track at 1.5m (b)
Cracks appearing on upper track at 7.2m (c)
Interface deformation on upper track at 18 mIn addition, there are rail welds at 10 m, 20 m, and 22 m, respectively.From the experimental simulation results and the actual situation verification, although the box diagram can roughly locate the abnormal point interval, the PCNN-wavelet is more accurate and closer to the actual situation. In Figure10(c), the simulative curve appears to be changed suddenly in these places. In other words, the cooperative anomaly detection method presented in this paper is effective.In fact, the emergence of data mutation points may be caused by the depressions, cracks, and deformation at the interface in the tracks. On the other hand, the response information of the damage position may be developed to make different damage positions in different information range. And the administrators of the industrial system may locate damage positions according to numerical characteristics in the damage unit. Furtherly, the relevant departments may design the good maintenance strategies according to the monitoring curves. That means our algorithm may be applied to the real project.


## 5. Conclusion

In this paper, a cooperative anomaly detection method for the stacker running track in the industrial environment is presented, which is based on PCNN and wavelet transform. Firstly, the data denoising model is built based on PCNN. Then, the data is detected by wavelet transform. Finally, the rationality and validity of the proposed method are verified by example analysis and simulation. The main conclusions are as follows:( 1 ) In ASRS, under the influence of complex industrial environment, the data collected by the sensors will inevitably be interfered by noise, which will bring difficulties in data analysis. Simulation and real experiment examples prove that PCNN performs well in noise reduction. Moreover, it can eliminate noise interference as well as retain the original signal characteristics.( 2 ) After getting clean data, the wavelet transform can effectively locate the location of the abnormal point and has good applicability.( 3 ) After using PCNN to eliminate noise, the interference of noise-to-signal data is reduced, which is more conducive to the detection of signal data anomaly spot by wavelet transform. Through the comparison of the field and experimental data, it is shown that using PCNN and wavelet transform can realize the detection of the orbital anomaly spot of stacking machine.



---

*Source: 1938490-2018-11-11.xml*