Next Article in Journal
Catalytic Hydrolysis of Tricresyl Phosphate by Ruthenium (III) Hydroxide and Iron (III) Hydroxide towards Sensing Application
Next Article in Special Issue
Precise Loran-C Signal Acquisition Based on Envelope Delay Correlation Method
Previous Article in Journal
Selective Detection of Folic Acid Using 3D Polymeric Structures of 3-Carboxylic Polypyrrole
Previous Article in Special Issue
Double-Constraint Inpainting Model of a Single-Depth Image
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Small Foreign Object Debris Detection for Millimeter-Wave Radar Based on Power Spectrum Features

Ministerial Key Laboratory of JGMT, Nanjing University of Science and Technology, Xiao Ling Wei200#, Nanjing 210094, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(8), 2316; https://doi.org/10.3390/s20082316
Submission received: 10 March 2020 / Revised: 15 April 2020 / Accepted: 15 April 2020 / Published: 18 April 2020
(This article belongs to the Special Issue Data, Signal and Image Processing and Applications in Sensors)

Abstract

:
Foreign object debris (FOD) detection can be considered a kind of classification that distinguishes the measured signal as either containing FOD targets or only corresponding to ground clutter. In this paper, we propose a support vector domain description (SVDD) classifier with the particle swarm optimization (PSO) algorithm for FOD detection. The echo features of FOD and ground clutter received by the millimeter-wave radar are first extracted in the power spectrum domain as input eigenvectors of the classifier, followed with the parameters optimized by the PSO algorithm, and lastly, a PSO-SVDD classifier is established. However, since only ground clutter samples are utilized to train the SVDD classifier, overfitting inevitably occurs. Thus, a small number of samples with FOD are added in the training stage to further construct a PSO-NSVDD (NSVDD: SVDD with negative examples) classifier to achieve better classification performance. Experimental results based on measured data showed that the proposed methods could not only achieve a good detection performance but also significantly reduce the false alarm rate.

1. Introduction

Foreign object debris (FOD) refers to any object located in the airport environment that can injure or harm the aircraft. Typical examples of FOD include aircraft parts, twisted metal strips, ramp garbage, plastic products, and so on [1]. FOD poses a safety risk to aircraft and a significant economic loss to airlines. The estimated value of FOD damage reaches several billion USD per year for airlines in terms of structural damage and maintenance costs [2,3]. After the accident of a Concorde airplane due to a metal strip on the runway, FOD detection has become an exciting research area and gained significant attention [4].
At present, the outstanding FOD detection systems are the Tarsier system by QinetiQ, FODFinder system by Trex, iFerret system by Stratech, and FODetect system by Xsight [5,6,7,8]. These systems mainly realize FOD detection through radar or electro-optical technologies [9]. They follow the standards of an advisory circular issued by the Federal Aviation Administration (FAA), and the performance of these systems has been tested. For example, as a monitoring system relying only on video image processing, the iFerret system can detect targets with a height of 4 cm, and a target with a height of 2 cm can be detected under good weather conditions. However, this system is highly vulnerable to haze weather [5]. The FODFinder system adopts millimeter-wave radar and video image processing to detect FOD targets, but it can only detect targets as small as 2.5 cm × 2.5 cm (diameter × height) [7]. Thus, these systems only satisfy the demands of civil airports. For military airports, it is highly desirable to detect much smaller FOD targets, such as a 2 cm diameter metal ball, though this is extremely difficult. For FOD detection systems that rely on millimeter-wave radar, they are generally able to detect targets with a radar cross-section (RCS) of not less than −30 dBsm. However, for a 2 cm diameter metal ball, the RCS is as low as −35 dBsm, which is very difficult to detect. As the FOD target becomes smaller, its RCS will decrease to a very small value, resulting in the target echo being close to the ground clutter. Therefore, the signal-to-clutter ratio (SCR) is greatly reduced, which makes the detection difficulty increase sharply. This paper mainly aims to solve the difficult problem of small FOD detection in military airports.
There are several studies have been published in recent years, which mainly focus on image processing and radar signal processing. Different from the traditional image processing methods, novel algorithms based on a convolutional neural network (CNN) are proposed in References [10,11] to detect FOD based on optical images, which achieve a high recall rate for small FOD, such as small screws. However, the images obtained at night or in hazy weather are not ideal due to the influence of light and brightness. Yigit et al. detected a 1.2 cm × 3.5 cm (diameter × height) screw at a distance of 2 m from the radar by using millimeter-wave, ground-based synthetic aperture radar imaging, while the detection capability of remote small FODs has not been verified [12]. In radar signal processing methods, the constant false alarm rate (CFAR) detection algorithm is commonly utilized for FOD detection. Cell average (CA) and ordered statistics (OS) CFAR methods have good detection performance with a homogeneous background, while the detection effect for small FOD targets is reduced with non-homogeneous backgrounds [13,14]. The clutter map (CM) CFAR algorithm can estimate the average clutter power through multiple iterations but the “self-shielding” is inevitable, leading to a reduction in detection performance [15]. Leonardi et al. proposed the fusion of raw signal level data for the echo data received by multiple radars, which improved the detection rate of small FOD targets [16]. However, for a small target whose SCR is equal to 4 dB, the detection rate of the method only reaches 35%, which is still a very low value. Therefore, the principle of radar automatic target recognition was applied in this study to transform the target detection into classification, and then implement FOD detection through a classifier [17]. To verify the practical effect of this method, we collected data at a military airport in Jiangsu, China. Based on the measured data, we compared the methods proposed in this paper with the traditional CA-CFAR and CM-CFAR algorithms. For a 2 cm diameter metal ball, the results showed that the detection rate of the CM-CFAR algorithm was slightly higher than that of the CA-CFAR algorithm, but the false alarm rate of the CA-CFAR algorithm was lower. However, the methods proposed in this paper not only achieved a higher detection rate than the CM-CFAR algorithm, but also obtained a lower false alarm rate than the CA-CFAR algorithm. Therefore, it is helpful to realize the difficulty of detecting small FOD targets in military airfields.
This paper is arranged as follows: Section 2 introduces the signal model of the linear frequency modulated continuous wave (LFMCW) radar used in this study and outlines the principle of the traditional CFAR algorithms for FOD detection. In Section 3, the principle of the classifier is introduced, and it is applied for FOD detection. The feature extraction is carried out on the data first, then the theory of the feature-based SVDD classifier is expounded. Finally, we describe the process of optimizing the classifier parameters using the particle swarm optimization (PSO) algorithm. Section 4 gives the detection results and analysis of comparing traditional FOD detection with classifier-based FOD detection. Some conclusions are collected in Section 5.
Notations: In this paper, the bold uppercase letters are matrices, the bold lowercase letters are vectors, and the italic letters show scalar variables.

2. Traditional CA-CFAR for FOD Detection

2.1. Signal Model of LFMCW Radar

The typical LFMCW scheme is used for millimeter-wave FOD detection radar systems. The application of LFMCW has resulted in reduced costs, power consumption, and hardware and software complexity. For these reasons, LFMCW radar is attractive in near-field and high-resolution applications. For a frequency sweep cycle, the saw-tooth LFMCW signal emitted by the radar can be expressed as:
s t ( t ) = A exp [ j 2 π ( f 0 t + 1 2 μ t 2 ) ] ,
where A and f 0 are the amplitude and center frequency of the transmission signal, respectively; μ = B r / T r is the sweeping slope; B r is the bandwidth of the sweep; and T r is the duration of a single frequency sweep. For a stationary FOD target on the airport runway, the echo signal received by the radar after the target reflection can be written as:
s r ( t ) = K a A exp { j 2 π [ f 0 ( t τ ) + 1 2 μ ( t τ 2 ) ] } ,
where K a is the reflection coefficient and τ = 2 R / c is the two-way time delay for a FOD target located at range R . Substituting τ into Equation (2), the echo signal is down-converted with the transmitted signal, then the beat frequency signal can be obtained as follows:
s b ( t ) = K a A 2 2 exp [ j 2 π ( μ 2 R c t μ 2 R 2 c 2 + 2 R f 0 c ) ] .
It can be seen that the signal s b ( t ) is a single-frequency signal with the frequency f = 2 μ R / c . Then, performing an N-point fast Fourier transform (FFT) on it, and the sampling frequency is f s . The signal after the FFT can be expressed as:
S b ( n ) = K a A 2 N 2 sin c ( n 2 μ R c N f s ) exp { j [ N 1 N ( n 2 μ R c f s N N 1 ) π ] } exp [ j 2 π ( 2 R f 0 c μ 2 R 2 c 2 ) ]
where sin c ( x ) = sin ( π x ) / ( π x ) and n = 1 ,   2 ,   ,   N 1 . After the FFT processing, the range of the FOD target can be estimated from the spectrum of the S b ( n ) signal.

2.2. Traditional CFAR Algorithm

The constant false alarm rate (CFAR) algorithm is a typical method for FOD detection. Consider a CA-CFAR detector, the basic principle of which is shown in Figure 1 [18], where s b ( t ) is the beat frequency signal and S b ( n ) is the signal after the FFT, which is fed into the CA-CFAR detector. The shadow parts on both sides of the cell under test (CUT) are guard cells, which do not participate in the estimation of the ground clutter power to avoid missing detection. The detection threshold can be expressed as shown in Equation (5) where p i ,   q i ( i ,   j = 1 ,   ,   n ) are the reference cells on both sides of the CUT; the power summation of the front and back edge reference cells are P and Q , respectively; Z is the average value of all the reference cells; and γ is the threshold coefficient:
T = γ Z = 1 2 n γ ( P + Q ) = 1 2 n γ ( i = 1 n p i + j = 1 n q j ) .
The CA-CFAR detector outputs the detection result as either FOD present or FOD absent according to whether the power of the CUT is greater than the threshold. For the CM-CFAR detector, it mainly includes two steps: updating the clutter map and detecting the target [18]. First, the clutter map is obtained through multiple iterations of ground clutter, and then the detection threshold is determined. Finally, similar to the CA-CFAR detector, the power of the CUT is compared with the threshold to obtain the detection result.

3. FOD Detection with the PSO-SVDD Classifier

FOD detection involves binary hypothesis testing and can also be regarded as a binary classification [19]. There are only two different decision results, namely FOD present or FOD absent. The small FOD targets detection model can be simply represented as:
{ FOD absent : S b ( n ) = S c ( n ) FOD present : S b ( n ) = S c ( n ) + S F ( n )
where S c ( n ) and S F ( n ) are the signals of ground clutter and the FOD target, respectively, and they are independent of each other. Therefore, this study aimed to use the PSO-SVDD classifier to replace the CFAR detectors for FOD detection. The data of signal S b ( n ) was categorized into two kinds, one was the data corresponding to ground clutter (with FOD absent), and the other was the data that contained the FOD target (with FOD present). As shown in Figure 2, the FOD absent data were taken as the training data in the training stage. First, the power spectrum features of each training sample were extracted and combined as an eigenvector to describe it. Then, all the eigenvectors were combined into an eigenvector matrix to train the SVDD classifier. Meanwhile, the PSO algorithm was used to optimize the parameters of the classifier. In the testing stage, the FOD absent and FOD present data were utilized to validate the classification performance. For each test sample, the power spectrum features were also extracted, and the eigenvector was input into the trained PSO-SVDD classifier to obtain the final classification result: FOD present or FOD absent.

3.1. Feature Extraction

As shown in Figure 3, for the cell S b ( u ) under test, k cells before and after it were combined into a small block of data, and feature extraction was carried out on the block of data.
The data in this block can be written as S = [ S b ( u k ) , , S b ( u 1 ) , S b ( u ) , S b ( u + 1 ) , , S b ( u + k ) ] , mapping it into the power spectrum domain; then, the power spectrum of it can be expressed as D = S 2 = [ D ( u k ) , , D ( u 1 ) , D ( u ) , D ( u + 1 ) , , D ( u + k ) ] . The power spectrum features can be extracted using the following equation:
f c = m = u k u + k ( m m ˜ ) 2 v ( m ) ,
which represents the second-order central moment of the power spectrum, where m ˜ = m = u k u + k m v ( m ) denotes the first-order central moment of the power spectrum and v ( m ) = D ( m ) / m = u k u + k D ( m ) is the normalized power spectrum.
f a = 1 2 k + 1 m = u k u + k D ( m )
represents the average power spectrum. Dimensionality reduction was achieved by extracting the two power spectrum features from the data, and the data with FOD present and FOD absent was converted into the power spectrum domain with higher discrimination. The two features could be combined into an eigenvector of x = ( f c ,   f a ) to describe the block of data.

3.2. Principle of SVDD

SVDD is a one-class classifier proposed by Tax and Duin, which is widely used in anomaly detection [20,21]. The SVDD models a class of data by fitting a hypersphere with center a and radius r around all or most of the samples. Assume that we are given a set of training samples { x i ,   i = 1 ,   2 ,   ,   l } , where l is the number of samples. The SVDD aims to minimize the volume of the hypersphere by minimizing r 2 . However, the training data will not be spherical in many cases; therefore, the kernel function is introduced to map it into the kernel space. The kernel function can be written as:
K ( x i , x j ) = Φ ( x i ) , Φ ( x j ) = Φ ( x i ) T Φ ( x j ) ,
where { Φ ( x i ) ,   i = 1 ,   2 ,   ,   l } is the data set mapped into the kernel space. The solution of the optimal decision surface of SVDD can be expressed as the following optimization problem:
{ min r , a , ξ i F ( r , a ) = r 2 + C i = 1 l ξ i s . t . Φ ( x i ) a 2 r 2 + ξ i ,   ξ i 0 ,   i = 1 ,   2 ,   ,   l .
The C parameter controls the tradeoff between the volume of the hypersphere and the number of target objects rejected. Since the training data may contain outliers, ξ i are slack variables that relax the constraints. A lower classification error rate and a smoother decision boundary containing the training samples can be obtained by adjusting C . The constraints can be incorporated into the objective function by using Lagrange multipliers:
L ( r , a , ξ i , α i , β i ) = r 2 + C i = 1 l ξ i i = 1 l β i ξ i i = 1 l α i { r 2 + ξ i Φ ( x i ) a } ,
where α i and β i are Lagrange multipliers, where α i , β i 0 . The above optimization is a minimax problem, which can be written as:
max α i , β i min r , a , ξ i L ( r , a , ξ i , α i , β i ) .
L should be minimized with respect to r , a , and ξ i and maximized with respect to α i and β i . Setting the partial derivatives to zero gives the constraints:
{ L r = 0 : i = 1 l α i = 1 L a = 0 : a = i = 1 l α i x i L ξ i = 0 : C α i β i = 0 .
Substituting Equation (13) into Equation (11) and introducing the kernel function, the final optimization problem can be obtained as follows:
{ max α L = i = 1 l α i K ( x i , x i ) i , j = 1 l α i α j K ( x i , x j ) s . t . i = 1 l α i = 1 ,   0 α i C ,   i = 1 ,   2 ,   ,   l .
According to the Karush–Kuhn–Tucher (KKT) conditions, only part of α i satisfies 0 < α i C and its corresponding x i is the support vector. The radius and center of the hypersphere are both determined by the support vectors. Suppose the number of support vectors is n S V , and S V s denotes the set of support vectors; then, the radius and center of the hypersphere are:
{ r 2 = K ( x i , x i ) 2 i = 1 n S V α i K ( x i , x i ) + i , j = 1 n S V α i α j K ( x i , x j ) a = i = 1 l α i Φ ( x i ) .  
In the SVDD model, to determine whether a new sample x * lies in the description, the distance from the center of the sphere to x * must be less than r 2 . Hence, the decision function for x * is:
d ( x * ) = Φ ( x * ) Φ ( a ) 2 r 2 .
Consequently, the discriminant result for the FOD detection is:
{ d ( x * ) 0 , FOD absent d ( x * ) > 0 , FOD present .
Several different choices of kernel functions exist. In this study, we used the well-known Gaussian kernel function. The Gaussian kernel has only one free parameter to be tuned and is shown to yield tighter boundaries than other kernel choices [22]. The Gaussian kernel is given by the following:
K ( x , y ) = exp ( x y 2 / 2 σ 2 ) ,
where σ is a free parameter that is adjusted to control the tightness of the boundary and is typically optimized through cross-validation [23]. In the kernel-based SVDD, the selection of free parameters C and σ affects the classification performance; therefore, they need to be optimized to achieve the best detection performance.

3.3. Parameter Optimization Based on the PSO

The grid method is adopted in traditional SVDD parameter optimization, which is time-consuming and may not find the optimal parameters. This study applied the PSO algorithm for parameter optimization, which is a global search evolutionary algorithm with a simple structure and low complexity [24]. The purpose of the PSO algorithm is to obtain the optimal parameters C and σ to minimize the classification error rate. For a testing set with N samples, its classification error rate is:
e r r = F P + F N N ,
where F P is the number of samples where FOD is judged as FOD absent and F N is the number of samples where FOD is judged as FOD present. The steps of the PSO can be described in Algorithm 1.
Algorithm 1. The optimization steps of parameters using PSO.
Step 1: Input the training and testing data, initialize the parameters C and σ of SVDD model, and set the searchable range of the parameters.
Step 2: Initialize the particle swarm, including the population size W , acceleration constants c 1 and c 2 , inertia weight ω , maximum number of iterations I t , and the particle speed and position.
Step 3: Determine the individual extremum of the initial position and the optimal position of the particle swarm.
Step 4: Calculate the fitness value of the new position of each particle in the swarm.
Step 5: Compare the current optimal position of each particle with the optimal position of the particle swarm and update the optimal solution to the current optimal position of particle swarm.
Step 6: Update the speed and position of the particle.
Step 7: Determine whether the SVDD model with the current parameters can minimize the error rate or reach the maximum number of iterations. If one of them is satisfied, the optimal parameters C and σ are obtained. Otherwise, return to step 4 to recalculate the particle fitness value.

4. Detection Results and Analysis

In this section, the simulations that were conducted to compare the PSO-SVDD classification method with the classical CA-CFAR and CM-CFAR algorithms are described. To obtain the measured data, a millimeter-wave radar was used to collect data at a military airport in Jiangsu, China, where the airport runways were paved with asphalt. The operation and the signal processing parameters of the radar system were as follows: the operation frequency was f 0 = 94.5   GHz , the bandwidth of the sweep was B r = 600   MHz , the duration of a single frequency sweep was T r = 1   ms , the transmitting power was 20 dBm, and the vertical and horizontal beamwidth of the antenna was 8° and 2°, respectively. The sampling frequency was f s = 1   MHz and the number of points for the FFT was N = 1024 . The radar sensor was placed at 1.5 m above the ground. Additionally, during the feature extraction, the number of the intercepted cells was set to k = 6 . For the CA-CFAR detector, there were two guard cells and four reference cells on either side. According to the principle of CA-CFAR, the threshold coefficient was set to γ = 3.5 in this paper.
Before collecting the data, the experimental scene was manually checked first to ensure that there were no FOD targets. Then, the millimeter-wave radar scanned the experimental scene and recorded the echo data of the ground clutter. After that, FOD targets were manually placed on the airport runway. Taking the detection of a 2 cm diameter metal ball as an example, we placed the target at a distance of 42 m from the radar. Then, the echoes of it were acquired using radar scanning. Figure 4a shows the amplitude spectrum when the metal ball was present and Figure 4b depicts the amplitude spectrum of the ground clutter (without the metal ball). It can be seen from Figure 4 that in the area closer to the radar, since the energy of the echo entered the radar from the sidelobe, the energy of the ground clutter echo was relatively strong. However, in the area far from the radar, the target echo received by the radar was weaker, which further increased the difficulty of detecting small FOD targets. According to the local magnified image near 42 m in Figure 4a, the SCR of the target was only 4 dB, which was almost submerged in the ground clutter. Since the ground clutter echo and the target echo were difficult to distinguish, it was very difficult to detect the target via the CA-CFAR and CM-CFAR algorithms. Therefore, we extracted the implicit internal features of the ground clutter echo and target echo to further enlarge the distinction between them, and then the classification method could be applied for the FOD detection.
Based on the method proposed in this paper, the two power spectrum features were extracted. Figure 5a,b shows the distribution of the two features when the FOD was present and absent. The second-order central moment describes the distribution of the power spectrum energy relative to its geometric centroid. The more dispersed the energy distribution, the smaller the second-order central moment. From Figure 5a, when the FOD was present, the second-order central moment was smaller than that when it was absent. For the average power spectrum, the greater the energy, the greater the average. From Figure 5b, it was obvious that the average power spectrum in the presence of the FOD was greater than that in the absence of FOD. Compared with Figure 4, by extracting these two features, the discriminability between FOD present and FOD absent was significantly improved. Then, the extracted features could be used to train the classifier or test the performance of it.
To verify the FOD detection performance of the PSO-SVDD detector, we counted 300 samples of ground clutter data and 150 samples of data containing FOD. In the training stage, 80% of the FOD absent samples were randomly selected to train the PSO-SVDD classifier, and the rest were taken as the testing data. Based on the analysis in Section 3.3, the parameters were first set to W = 30 , c 1 = 2 , c 2 = 2 , ω = 0.5 , I t = 200 , 0 < C < 1 , and 0 < σ < 2 8 . Then, the PSO algorithm was used to optimize the parameters of the SVDD model. The optimization process is shown in Figure 6, where the minimum classification error rate was 25.71% after 11 iterations, and the optimal decision boundary could be obtained with parameters C = 0.63 and σ = 63.77 . In the classification process, the test samples inside the decision boundary were judged FOD absent; otherwise, they were judged as FOD present.
The classification results of the PSO-SVDD are depicted in Figure 7. In Figure 7a, the red circle markers denote the feature vectors of the samples without FOD, which were utilized to train the PSO-SVDD classifier. In Figure 7a,b, the green circle markers represent the support vectors, and the decision boundary is denoted by the purple line. The decision boundary distinguished the FOD present samples from FOD absent samples in the testing procedure. It can be seen from Figure 7b that there were two FOD absent samples that were classified outside the decision boundary, that is, they were judged as FOD present, resulting in false alarms. Furthermore, the remaining samples with FOD absent were all divided into the interior of the decision boundary, indicating that the judgment was correct. Meanwhile, most of the FOD present samples were correctly divided outside the decision boundary, while a small part of them were wrongly classified inside as FOD absent.
Since only the FOD absent samples were utilized to train the classifier in the training procedure, overfitting inevitably occurred. To improve the performance further, 10% of the samples containing FOD were randomly added to train a PSO-NSVDD (NSVDD: SVDD with negative examples) classifier. The principle of the NSVDD model is similar to that of the SVDD model [25], while the parameters that need to be optimized are C + , C , and σ , where C + and C are the penalty parameters for the FOD absent and FOD present samples, respectively, and σ is the kernel width of the NSVDD model. After 23 iterations, the optimal parameters of the model were C + = 0.76 , C = 0.58 , and σ = 12.87 . The final detection results of the PSO-NSVDD classifier are shown in Figure 8, where the yellow circles in Figure 8a indicate the newly added FOD present samples. Furthermore, it can be seen from Figure 8b that, due to the addition of samples with FOD present in the training stage, the decision boundary was described more accurately and the overfitting was avoided. Therefore, compared with the results of the PSO-SVDD detector in Figure 7b, more samples with FOD present were classified outside the decision boundary and the FOD detection performance was further improved. Meanwhile, almost all the FOD absent samples were correctly classified within the decision boundary, which further reduced the false alarm rate.
To detect a target that exists at any distance from the radar, according to the method proposed in this paper, we first divided the collected ground clutter data into several small blocks in the range dimension. Similarly, we also divided the data containing the minimum FOD target to be detected into small blocks. Then, the feature extraction was performed on each block, and finally, the optimal PSO-SVDD and PSO-NSVDD classifiers were trained in these blocks. Within each block, there was little variation of ground clutter echo and the PSO-NSVDD classifier trained by the data containing the minimum FOD target could also be used to detect larger targets. Therefore, the parameters of these classifiers could be fixed. In the later detection process, there was no need to retrain new classifiers, which meant that these trained classifiers could be directly utilized to detect the FOD targets that existed at any distance. Therefore, we placed a 2 cm diameter metal ball and a 4.3 cm diameter golf ball at a distance of 10 m from the radar, and placed another golf ball with the same size at a distance of 42 m to verify the detection effect. We also counted 150 samples containing these FOD targets for testing. Table 1 shows the detection results of these targets using four detection methods. For the metal ball and golf ball, SCR = 19 dB and 30 dB at 10m and SCR = 4 dB and 13 dB at 42 m, respectively. According to the detection rate P d and the false alarm rate P f of the metal ball at different ranges, compared with the CA-CFAR algorithm, the detection rate of the CM-CFAR algorithm was slightly higher, but the false alarm rate was also higher, which was not what we expected. For the CA-CFAR and CM-CFAR algorithms, the golf ball was detected much easier than the metal ball, suggesting that it was difficult for the conventional CFAR methods to detect such small-volume FOD targets. However, for the metal ball at two ranges, the P d of the PSO-SVDD classifier was significantly higher, reaching 87.33% at 10 m and 64.67% at 42 m, and the P f was significantly lower for both. Additionally, after adding the samples with FOD present, the P d of the PSO-NSVDD was further increased to 92.61% at 10 m and 78.18% at 42 m. Meanwhile, for both the golf ball and metal ball, the PSO-NSVDD classifier produced a much lower false alarm rate. Based on the analysis above, we have reason to believe that the methods proposed in this paper not only showed a higher detection rate, but also a lower false alarm rate for the target with the low SCR. For the target with the high SCR, the detection rate was the same and the false alarm rate was even lower. Thus, these proposed methods were effective for a wide SCR range.
Finally, the two methods proposed in this paper were comprehensively analyzed as follows. For the PSO-SVDD classifier, only the ground clutter (with FOD absent) data needed to be utilized in the training stage. In an actual airport environment, a mass of ground clutter data can be collected to train the classifier; therefore, this method is easily implemented. For the PSO-NSVDD classifier, a small number of samples with FOD present needed to be added during the training stage, meaning that the data containing the minimum FOD target needed to be collected in advance, while it is unrealistic to collect such data in practical applications. However, the detection performance of the latter was further improved. Furthermore, compared with the traditional CA-CFAR and CM-CFAR methods, the detection performance of these two methods for small FOD targets was much higher.

5. Conclusions

Small FOD target detection under strong ground clutter is a key issue for millimeter-wave radar FOD detection, which is difficult to achieve using traditional CFAR algorithms. This paper transformed FOD detection into a kind of classification, and the PSO-SVDD and PSO-NSVDD classifiers were used to classify the FOD target and ground clutter to detect FOD from the background. The experimental results based on the measured data showed that the proposed methods achieved a higher detection rate and lower false alarm rate for small FOD targets than the classical CFAR algorithms.

Author Contributions

Conceptualization, C.M. and P.N.; methodology, C.M. and P.N.; software, P.N.; validation, C.M., P.N., and H.T.; formal analysis, C.M.; investigation, C.M.; resources, P.N.; data curation, P.N.; writing—original draft preparation, P.N.; writing—review and editing, C.M., P.N., and M.J.; visualization, P.N.; supervision, W.W.; project administration, W.W.; funding acquisition, W.W. All authors have read and agreed to the published version of the manuscript.

Funding

The research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xu, H.; Han, Z.; Feng, S.; Zhou, H.; Fang, Y. Foreign object debris material recognition based on convolutional neural networks. EURASIP J. Image Video Proc. 2018, 21. [Google Scholar] [CrossRef] [Green Version]
  2. Wang, X.; Lan, Z.; Li, Z.; Lin, Y.; Xiao, Q. A practical application of runway Foreign Object Debris detection system at the airport. In Proceedings of the 2018 International Conference on Microwave and Millimeter Wave Technology, Chengdu, China, 7–11 May 2018. [Google Scholar]
  3. Mollo, G.; Napoli, R.D.; Naviglio, G.; Chiara, C.D.; Capasso, E.; Alli, G. Multifrequency Experimental Analysis (10 to 77 GHz) on the Asphalt Reflectivity and RCS of FOD Targets. IEEE Geosci. Remote Sens. Lett. 2017, 14, 1441–1443. [Google Scholar] [CrossRef]
  4. Mehdi, G.; Miao, J. Millimeter wave FMCW radar for Foreign Object Debris (FOD) detection at airport runways. In Proceedings of the IEEE International Bhurban Conference on Applied Sciences & Technology, Islamabad, Pakistan, 9–12 January 2012; pp. 407–412. [Google Scholar]
  5. iFerret on Scratch. Available online: https://www.sourcesecurity.com/news/stratech-iferret-assess-ground-surveillance-systems-co-7811-ga.19664.html/ (accessed on 24 January 2019).
  6. Tarsier®: Automatic Runway FOD Detection System. Available online: https://www.tarsierfod.com/ (accessed on 19 December 2018).
  7. FOD Finder™. Available online: https://www.xsightsys.com/index.php/fodetect/ (accessed on 4 January 2016).
  8. What is FODetect? Available online: https://www.xsigthsys.com/fodetect.html (accessed on 24 January 2018).
  9. Zhang, H. The current status and inspiration of FOD industry in international civil aviation. Civ. Aviat. Manag. 2015, 295, 58–61. [Google Scholar]
  10. Cao, X.; Wang, P.; Meng, C.; Bai, X.; Gong, G.; Liu, M.; Qi, J. Region Based CNN for Foreign Object Debris Detection on Airfield Pavement. Sensors 2018, 18, 737. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Liu, Y.; Li, Y.; Liu, J.; Peng, X.; Zhou, Y.; Murphey, Y. FOD detection using DenseNet with focal loss of object samples for airport runway. In Proceedings of the IEEE Symposium Series on Computational Intelligence, Bangalore, India, 18–21 November 2018; pp. 547–554. [Google Scholar]
  12. Yigit, E.; Demirci, S.; Unal, A.; Ozdemir, C.; Vertiy, A. Millimeter-wave ground-based synthetic aperture radar imaging for Foreign Object Debris detection: Experimental studies at short ranges. J. Infrared Milli. Terahertz Waves 2012, 33, 1227–1238. [Google Scholar] [CrossRef]
  13. Tom, A.; Viswanathan, R. Switched order statistics CFAR test for target detection. In Proceedings of the IEEE Radar Conference, Rome, Italy, 26–30 May 2008. [Google Scholar]
  14. Zattouta, B.; Farrouki, A.; Barkat, M. Automatic censoring detection using binary clutter-map estimation for non-gaussian environments. In Proceedings of the IEEE International Conference on Signal Processing and Communications, Dubai, UAE, 24–27 November 2007; pp. 205–208. [Google Scholar]
  15. Nitzberg, R. Clutter map CFAR analysis. IEEE Trans. Aerosp. Electron. Syst. 1986, 22, 419–421. [Google Scholar] [CrossRef]
  16. Leonardi, M.; Piracci, E.G.; Fastella, V. W-Band multiple-radar processing for airport Foreign Object Debris and humans detection. In Proceedings of the 1st European Microwave Conference in Central Europe, Prague, Czech Republic, 13–15 May 2019; pp. 285–288. [Google Scholar]
  17. Du, L.; Wang, B.; Wang, P.; Ma, Y.; Liu, H. Noise reduction method based on principal component analysis with beta process for micro-doppler radar signatures. IEEE J. Sel. Top. Appl. Earth Observ. Remote Sens. 2015, 8, 4028–4040. [Google Scholar] [CrossRef]
  18. Li, H.; Zhang, Z.; Wang, Y.; Yu, X. Performance analysis and comparison of CFAR methods for FOD detection in airport runway environment. Radio Eng. 2015, 45, 53–57. [Google Scholar]
  19. Wang, L.; Tang, J.; Liao, Q. A study on radar target detection based on deep neural networks. IEEE Sens. Lett. 2019, 3. [Google Scholar] [CrossRef]
  20. Tax, D.M.J.; Duin, R.P.W. Support vector domain description. Pattern Recognit. Lett. 1999, 20, 1191–1199. [Google Scholar] [CrossRef]
  21. Tax, D.M.J.; Duin, R.P.W. Support vector data description. Mach. Learn. 2004, 54, 45–66. [Google Scholar] [CrossRef] [Green Version]
  22. Munoz-Marf, J.; Bruzzone, L.; Camps-Vails, G. A support vector domain description approach to supervised classification of remote sensing images. IEEE Trans. Geosci. Remote Sens. 2007, 45, 2683–2692. [Google Scholar] [CrossRef]
  23. Banerjee, A.; Burlina, P.; Diehl, C. A support vector method for anomaly detection in hyperspectral imagery. IEEE Trans. Geosci. Remote Sens. 2006, 44, 2282–2291. [Google Scholar] [CrossRef]
  24. Liu, Y.; Qin, Z.; Shi, Z.; Lu, J. Center particle swarm optimization. Neurocomputing 2007, 70, 672–679. [Google Scholar] [CrossRef]
  25. Mu, T.; Nandi, A.K. Multiclass classification based on extended support vector data description. IEEE Trans. Syst. Man Cybern. Part B Cybern. 2009, 39, 1206–1212. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Basic principle of the cell average constant false alarm rate (CA-CFAR) detector. FFT: Fast Fourier Transform, FOD: Foreign Object Debris.
Figure 1. Basic principle of the cell average constant false alarm rate (CA-CFAR) detector. FFT: Fast Fourier Transform, FOD: Foreign Object Debris.
Sensors 20 02316 g001
Figure 2. Flowchart of the FOD classification. PSO: Particle Swarm Optimization.
Figure 2. Flowchart of the FOD classification. PSO: Particle Swarm Optimization.
Sensors 20 02316 g002
Figure 3. Flowchart of the feature extraction.
Figure 3. Flowchart of the feature extraction.
Sensors 20 02316 g003
Figure 4. (a) Amplitude spectrum with FOD present. (b) Amplitude spectrum with FOD absent.
Figure 4. (a) Amplitude spectrum with FOD present. (b) Amplitude spectrum with FOD absent.
Sensors 20 02316 g004
Figure 5. (a) Feature 1: the second-order central moment of the power spectrum. (b) Feature 2: the average power spectrum.
Figure 5. (a) Feature 1: the second-order central moment of the power spectrum. (b) Feature 2: the average power spectrum.
Sensors 20 02316 g005
Figure 6. The optimization process of the PSO algorithm.
Figure 6. The optimization process of the PSO algorithm.
Sensors 20 02316 g006
Figure 7. Detection results with the PSO support vector domain description (PSO-SVDD): (a) training procedure and (b) testing procedure.
Figure 7. Detection results with the PSO support vector domain description (PSO-SVDD): (a) training procedure and (b) testing procedure.
Sensors 20 02316 g007
Figure 8. Detection results with the PSO-NSVDD: (a) training procedure and (b) testing procedure.
Figure 8. Detection results with the PSO-NSVDD: (a) training procedure and (b) testing procedure.
Sensors 20 02316 g008
Table 1. The detection results using different methods.
Table 1. The detection results using different methods.
MethodsMetal Ball d = 2   cm , R = 10   m Golf Ball d = 4.3   cm , R = 10   m Metal Ball d = 2   cm , R = 42   m Golf Ball d = 4.3   cm , R = 42   m
P d ( % ) P f ( % ) P d ( % ) P f ( % ) P d ( % ) P f ( % ) P d ( % ) P f ( % )
CA-CFAR46.673.21001.3513.334.671001.11
CM-CFAR53.147.411004.5821.577.681003.14
PSO-SVDD87.331.251000.6564.670.951000.83
PSO-NSVDD92.610.391000.178.180.511000.12

Share and Cite

MDPI and ACS Style

Ni, P.; Miao, C.; Tang, H.; Jiang, M.; Wu, W. Small Foreign Object Debris Detection for Millimeter-Wave Radar Based on Power Spectrum Features. Sensors 2020, 20, 2316. https://doi.org/10.3390/s20082316

AMA Style

Ni P, Miao C, Tang H, Jiang M, Wu W. Small Foreign Object Debris Detection for Millimeter-Wave Radar Based on Power Spectrum Features. Sensors. 2020; 20(8):2316. https://doi.org/10.3390/s20082316

Chicago/Turabian Style

Ni, Peishuang, Chen Miao, Hui Tang, Mengjie Jiang, and Wen Wu. 2020. "Small Foreign Object Debris Detection for Millimeter-Wave Radar Based on Power Spectrum Features" Sensors 20, no. 8: 2316. https://doi.org/10.3390/s20082316

APA Style

Ni, P., Miao, C., Tang, H., Jiang, M., & Wu, W. (2020). Small Foreign Object Debris Detection for Millimeter-Wave Radar Based on Power Spectrum Features. Sensors, 20(8), 2316. https://doi.org/10.3390/s20082316

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop