This paper presents a method to identify landmines in various burial conditions. A ground penetration radar is used to generate data set, which is then processed to reduce the ground effect and noise to obtain landmine signals. Principal components and Fourier coefficients of the landmine signals are computed, which are used as features of each landmine for detection and identification. A database is constructed based on the features of various types of landmines and the ground conditions, including the different levels of moisture and types of ground and the burial depths of the landmines. Detection and identification is performed by searching for features in the database. For a robust decision, the counting method and the Mahalanobis distance-based likelihood ratio test method are employed. Four landmines, different in size and material, are considered as examples that demonstrate the efficiency of the proposed method for detecting and identifying landmines.
1. Introduction
Landmine removal is a critical problem faced by many countries around the world, and the situation can be compounded by natural disasters or land development. Therefore, it is an urgent issue to detect landmines in the ground and remove them safely. The process of landmine removal starts with the detection of landmines in the ground. For safe detection, non-touch-based detection methods are required. These methods involve the detection of landmines in the signals obtained by non-touch-based sensors, such as metal detectors and radars. Among those sensors, ground penetrating radars, or GPRs, are an attractive choice for landmine detection due to their advantages over other sensors. The GPR can be used as a stand-alone sensor or as a complementary sensor to a metal detector [1, 2]. It can detect both metal and nonmetal landmines [3]. Moreover, its weight can be made light, so that it can be installed in a hand-held system or in a vehicle-mounted system in the form of an array of multiple antenna elements [4–6].
In the landmine detection step, the key factor is to obtain unique signatures of a landmine from the signal, which are used as tags of each landmine. The signal may be, however, contaminated by noise, surface reflections, and so forth. Therefore, minimizing such influences and extracting the unique signatures for each landmine is an active research topic. Many studies show various methods of computing signatures from GPR data [7–12]. Tantum et al. [13] compare different algorithms for landmine detection. Each approach, however, utilizes one feature for detection, which may not be enough for robust detection because the discernibility of a single feature could be compromised by disturbances to the signal. This problem can be avoided with multiple features extracted by different methods from a signal using by sensor with no parameter changes. These features reflect different aspects of a landmine signal and would enhance their discernibility to a large extent. One paper [14] presents feature fusion, which includes either features obtained by different sensors or features obtained by one sensor with different parameter settings. This concept, however, is not useful because it requires different sensors or one sensor with different parameter settings, which is not a practical configuration.
The landmine detection step may be followed by landmine identification, which is a process to identify the types of landmines and their burial depths. This additional information can help detect and remove landmines more efficiently. Identification involves more steps than detection and has not been significantly investigated in the related field.
In this study, the problem of detection and identification of landmines using GPRs is addressed, and a novel method of detection and identification using multiple features is proposed. Principal component analysis and Fourier Transform are employed to obtain eigenvalues and Fourier coefficients, respectively, which are used as features associated with each landmine. Algorithms for detection and identification using multiple features are presented and tested with various examples.
2. Overall Process
The overall process of the proposed pipeline of landmine detection and identification is illustrated in Figure 1. The input of the process is the GPR signal, and part of the signal containing landmine effects is extracted after the noise and ground effects are reduced. The principal component analysis (PCA) and Fourier Transform (FT) methods compute features from the landmine signal. A database is constructed by storing the features with the associated landmine information. This database construction uses experimental data with known landmine information. For detection and identification, the input signal is processed to generate features, and similar features are identified in the database. The decision for detection and identification employs two methods: the counting value method first performs detection, followed by identification by the likelihood ratio test method. Once a landmine is successfully identified, the buried depth of the landmine is estimated. In the subsequent sections, the theories and algorithms of each step are introduced.
The overall process of landmine detection and identification.
2.1. Hardware and GPR Signals
The GPR system used in this work has a GPR antenna at the end of an arm that moves horizontally, covering the scanning area, as shown in Figure 2. The specification of the GPR is summarized in Table 1.
The specification of the GPR used for the test.
Property
Value
Internal pulse duration
240 ps (−6 dB width) unipolar
Internal pulse amplitude
15 V Peak
Pulse repetition frequencies
1 MHz
Bandwidth of operation
250 MHz to 2.5 GHz
EIRP
<−41 dBm/MHz
Radar sampling
512 samples per scan
Receiver time range
19.2 ns
Receiver sampling interval
37.5 ps
Output scan rate
61 Hz
Sweep speed
<1.5 ms^{−l}
The GPR system used in the experiment.
Targets are buried at different depths in a 2m×1m×1m wooden box. The types of ground, the moisture levels of the ground, and the antenna heights and the target depths are changed for each scan. The antenna acquires 256 depth data samples at each location and scans the survey area to yield a dense data set, as shown in Figure 3. Next, all the data are regulated into 150 by 253 matrices with 5 mm intervals. The rows and the columns of the matrices indicate sampling locations and depths, respectively.
An example of raw GPR data.
2.2. Reduction of Noise and Ground Effects
The GPR data (D) are assumed to be the sum of separate signals consisting of three components [7]: ground (G), landmine (L), and noise (N). Therefore, it is represented by D = G + L + N. From this model, the signal L is obtained by subtracting G + N from D. Scanning a region without a landmine would generate G + N. Although this method is conceptually straightforward, there are two problems: one is that the signal should be normalized to eliminate the radar power difference for each scan, and the other is that maintaining the same ground condition with and without a landmine is nearly impossible. Therefore, the approach cannot yield L consistently. An alternative approach is to select the reference signal, a signal obtained from D with the least amount of landmine influences, and to subtract it from D. The reference signal is computed as the average of two columns of data far away from the column of the largest strength. This approach is advantageous in that there is no need to consider the ground condition and the radar signal normalization to obtain the G + N signal since it is obtained in the same data set, eliminating the need for normalization. One example of GPR data is shown in Figure 3. The reference signal is subtracted to yield L, which is shown in Figure 4.
The GPR data after reducing the ground effects and noise.
2.3. Normalization
Once the landmine signal (L) is obtained, the intensity of the signal is normalized to reduce the influence caused by the difference of signal strength. The maximum and minimum intensity values, Imax and Imin, are selected from the signal L, which are then mapped to the values of Nmax and Nmin, the limiting values for normalization. Then the normalized signal Inormalized is computed by
(1)Inormalized=(I-IminImax-Imin)×(Nmax-Nmin)+Nmin.
Here, I is an intensity value of L.
3. Feature Extraction and Database Construction
Landmines are made of various materials, such as plastics and metals, and are different in size and shape. Therefore, they respond differently and bear unique signatures. In this work, two methods are used to extract such features from the signal L: principal component analysis (PCA) and Fourier Transform (FT).
3.1. Principal Component Analysis
PCA is a statistical method that can provide patterns of scattered data, N, in K-dimension. It constructs a K×K covariance matrix from the data. The eigenvectors and eigenvalues of the matrix are Ec=(e1,e2,…,eK) and Ev=(e1,e2,…,eK), respectively. A value of ei indicates how strong a pattern the data shows in the ei direction. Using this property, a signal of a landmine would contain eigenvalues and eigenvectors that are unique to the landmine. In this work, each GPR data point in 3D space as shown in Figure 4 is projected onto the XZ plane to generate a set of points in the 2D space as shown in Figure 5. Here, the X axis indicates the column number and Z is the intensity. Then the data in the XZ plane are processed to compute two principal components (N=2). Therefore two eigenvalues and eigenvectors are obtained, which can be considered as features of each landmine. Nondimensional eigenvectors could be considered for GPR landmine discrimination in the similar way as [15]. As the GPR data show, however, the eigenvectors would not be used for GPR landmine discrimination in this work. Most of the GPR data are clustered near the zero intensity and one principal component should follow that distribution pattern as e_{1} in Figure 5. Therefore, the eigenvectors for different GPR data show little difference and cannot be used as a feature. The difference between GPR landmine signals, however, can be observed by comparing the distribution away from the zero intensity, Z = 0, which is reflected in the second eigenvalue. Namely, as shown in Figure 5, the eigenvalue in the direction of e_{2} would reflect such a unique pattern, which is used as a feature for a landmine.
GPR data projected onto the XZ plane and eigenvectors of the data set using the PCA method.
3.2. Fourier Transform
Fourier Transform is a way to represent a signal in the time domain as a combination of frequency components. Consider a column of the GPR data as shown in Figure 6. The signal represented in the frequency domain is in Figure 7. The components C1 and C2 at two frequency values, f1 and f2, respectively, are used as the features of the given signal. The choice of these two frequencies can be made as follows. Consider la landmines and lb burial conditions. Then, la·lb GPR signals can be obtained, each of which is transformed to yield la·lb signals in the frequency domain. Then, the frequencies of the coefficients which can differentiate each landmine as much as possible are chosen after all the signals in the frequency domain are investigated. This approach can be easily extended to include more than two frequency components as features of each landmine.
A column of the signal of maximum intensity.
The FT result of the signal.
3.3. Database Construction
A series of experiments for various landmines buried in diverse ground conditions and depths is performed and a database is constructed for storing the features, associated landmine information, and ground conditions.
4. Landmine Detection and Identification
Landmine detection is a process to determine whether a signal from GPR contains a landmine, whereas identification includes not only detecting a landmine in the signal but also determining the type of landmine. Therefore, identification is more involved than detection, providing more information about a landmine.
Given a feature extracted from the GPR data, detection and identification is performed by searching for the same feature in the database and retrieving the corresponding landmine information. The possibility that exactly the same feature is found in the database is slim, however, due to various noise and uncertainty contained in the acquired signal. Instead, it is highly likely that a feature close to the given one is found, or multiple landmines with features close to it are returned. Therefore, for robust detection and identification, a systematic approach should be used to eliminate confusion in the decision. In this work, two methods, the counting and Mahalanobis distance-based likelihood ratio test methods, are proposed.
4.1. Counting Method
The counting method counts the number of matches and chooses the one with the largest number of counts as the identified landmine. Each landmine is given a count variable. A given feature is compared with those in the database. If a match is found, the count variable corresponding to the landmine is increased by one. Here, the matching criterion is to check if the absolute difference between the feature and that in the database is within a given tolerance. There may be multiple matches for one feature. In this case, all of the count variables that match are increased by one. This counting process is performed for all of the features considered. The one with the largest counting value is reported as the detected landmine with its identification information.
4.2. Mahalanobis Distance-Based Likelihood Ratio Test Method
Consider an n-dimensional vector x = (x1,x2,…,xn). The Gaussian probability distribution function for x is given by
(2)g(μ,Σ)(x)=12πn|Σ|e-(1/2)(x-μ)TΣ-1(x-μ).
Here, μ is the n-dimensional mean vector (E[x1],E[x2],…,E[xn]) and Σ the n×n covariance matrix Σ = (Cov[xi,xj]) and (i,j=1,2,…,n). Assume that there are two data groups G1 and G2 with Gaussian distributions g1(x) and g2(x), respectively. Suppose that P(x) is the probability that an input x is a landmine and P(ωl) is the probability that ωl is a member of Gl. If the probability of x belonging to G1, denoted as P(ω1∣x), is greater than that of x belonging to G2, P(ω2∣x), x is decided to be a member of G1. Applying Bayes’ theorem, we have
(3)P(ω1∣x)>P(ω2∣x)⇒P(x∣ω1)P(ω1)P(x)>P(x∣ω2)P(ω2)P(x),
which is reduced to
(4)Λ(x)=P(x∣ω1)P(x∣ω2)>P(ω2)P(ω1).
The ratio Λ(x) is called the likelihood ratio, and the test using this ratio is the likelihood ratio test (LRT). Therefore, an input vector x is decided to belong to G1 if the ratio Λ(x) satisfies the condition (4). Otherwise, it belongs to G2. In the LRT, the Mahalanobis distances from x to Gk (k=1,2,…,r) are used as input, which is a useful scheme for evaluating similarity between data sets [16]. The distance DM(x) of a vector x from a data set with mean μ in n-dimension is given by [16]
(5)DM(x)=(x-μ)TS-1(x-μ).
Here, S is the pooled covariance matrix of the two data sets, which is the weighted average of the covariance matrices involved.
5. Experiments
In this section, a series of experiments is performed to demonstrate the performance of the proposed method. Four types of landmines (two for antihuman, M14 and M16, and two for antitank, M15 and M19) are used for the experiments. The relative sizes of the landmines are M14<M16<M15<M19. M15 and M16 are made of metal, while M14 and M19 are made of plastic. Various ground conditions are prepared in the experiments that consider three moisture levels (dry, moderate, and damp), three materials (sand, soil, and gravel) and seven burial depths (from 0 cm to 30 cm with a 5 cm interval). Each GPR data set is obtained by repeating the scanning 20 times to ensure the accuracy of the scanned data.
5.1. Preparation of GPR Data
The GPR data sets are processed to reduce ground effects from the landmine signals, which are then normalized based on the maximum and minimum signal values.
5.2. Features
In the experiments, the second eigenvalues from the principal component analysis and the coefficients at two frequencies, 0.71 GHz and 0.95 GHz, are selected as features for each landmine. These three features are chosen after checking the differentiability of the landmines considered in the experiment.
For each landmine in all the ground conditions, those three features are computed and stored in the database. Figure 8 shows the variation of the second eigenvalues of each landmine with respect to the depths. Figure 9 shows the variation of the Fourier coefficients at 0.71 and 0.95 GHz with respect to the burial depths for the dry sand condition. Note that the two curves for the M15 and M19 landmines in Figure 9(b) cross. However, if the eigenvalues and the coefficients at 0.71 GHz are used first, those two landmines can be discriminated. The coefficients at 0.95 GHz are able to differentiate the upper two curves from the rest. The matching tolerance for the counting method is ±5%.
The curves of the second eigenvalues of each landmine with respect to the depth for the dry sand condition.
The curves of the Fourier coefficients for (a) 0.71 GHz and (b) 0.95 GHz of the data for each landmine with respect to the depth for the dry sand condition.
To use the LRT, the second eigenvalues and the coefficients at 0.71 GHz are plotted in a 2D plane as shown in Figure 10, which shows the plot of the second eigenvalues and the coefficients at 0.71 GHz of the four landmines buried at different depths in the dry sand. As illustrated in the figure, each landmine is clustered irrespective of burial depths and the LRT can be successfully applied for landmine identification. If a set of test input data is given as in Figure 10, they can be identified as M15 or M19 depending on the positions of the input data in the plot.
The plots of the second eigenvalues and the Fourier coefficients for 0.71 GHz of each landmine of variable burial depths for the dry sand condition.
With the arbitrary noise added to the landmine signals, the proposed method generates satisfactory detection and identification results. Among all of the experiments, one case with landmines buried in sand is presented. Depending on the noise levels added to the signal, from ±1.5% to ±9.6%, the detection rates are 100% for all cases, and the identification ratio changes from 100% to 79.2%, as shown in Figure 11, which illustrates the variation of the identification rates with respect to the noise levels. Additionally, the higher the moisture levels are, the lower the identification rates become. The detection rates for the soil and gravel conditions are also 100% in the experiments. However, the identification rates vary from approximately 30% to 100% depending on the ground condition and the added noise levels. If both methods are compared for each case, the LRT method works better than the counting method for the soil and gravel conditions, whereas the counting method performs better for the cases with signals contaminated by a large amount of noise and landmines buried in sand.
The identification rate variation in sand.
6. Conclusions
In this paper, we propose a procedure using data from a GPR to identify landmines buried in the ground of diverse conditions. The proposed method uses multiple features that reflect different aspects of landmines, which enables more stable detection and identification of a landmine. Moreover, using a series of experiments with respect to burial depths, the depth at which an identified landmine exists could also be estimated. Such additional information about a landmine could enhance the landmine removal process.
The proposed method is limited in that it has been tested with only a single landmine case. Therefore, if there are multiple landmines and other foreign objects, the proposed method needs to be refined to handle such cases. The proposed method requires thorough evaluation with field data before it is used in practice. In addition, extending the work to the cases of multiple landmines with diverse objects is recommended for future research.
Acknowledgments
This work was supported by Unmanned Technology Research Center, Defense Acquisition Program Administration, and Agency for Defense Development.
AmazeenC. A.LockeM. C.Developmental status of the U.S. Army's new handheld standoff mine detection system (HSTAMIDS)Proceedings of the 2nd International Conference on the Detection of Abandoned Land MinesOctober 19984581931972-s2.0-0032306014AmazeenC. A.LockeM. C.US Army's new handheld standoff mine detection system (HSTAMIDS)Proceedings of the EUREL International Conference on The Detection of Abandoned Land Mines: A Humanitarian Imperative Seeking a Technical SolutionOctober 1996431172176MacDonaldJ.LockwoodJ. R.McFeeJ.AltshulerT.BroachT.CarinL.HarmonR.RappaportC.ScottW.WeaverR.Alternatives for landmine detectionSatoM.KobayashiT.TakahashiK.FujiwaraJ.FengX.Vehicle mounted SAR-GPR and its evaluation6217Detection and Remediation Technologies for Mines and Minelike Targets XIApril 2006Proceedings of SPIE2-s2.0-3374734287710.1117/12.665416HarmonR. S.HollowayJ. H.BroachJ. T.Processing of GPR data from NIITEK landmine detection system5089The International Society for Optical Engineering: Detection and Remediation Technologies for Mines and Minelike Targets VIIIApril 200313751382Proceedings of SPIE2-s2.0-024271246010.1117/12.497133YarovoyA.SavelyevT.ZhugeX.AubryP.John SchavemakerL. L.TettelaarP.Den BreejenE.Performance of UWB array-based radar sensor in a multi-sensor vehicle-based suit for landmine detectionProceedings of the 5th European Radar Conference (EuRAD '08)October 20082882912-s2.0-63149149484BrunzellH.Detection of shallowly buried objects using impulse radarGaderP. D.MystkowskiM.ZhaoY.Landmine detection with ground penetrating radar using hidden Markov modelsYuS. H.MehraR. K.WittenT. R.Automatic mine detection based on ground penetrating radarDetection and Remediation Technologies for Mines and Minelike Targets IVApril 1999Orlando, Fla, USA961972Proceedings of SPIE2-s2.0-0032656782CarevicD.Clutter reduction and target detection in Ground Penetrating Radar data using waveletsDetection and Remediation Technologies for Mines and Minelike Targets IVApril 1999Orlando, Fla, USA973978Proceedings of SPIE2-s2.0-0032637839RattoC. R.TorrioneP. A.CollinsL. M.Context-dependent feature selection for landmine detection with ground-penetrating radar7303Detection and Sensing of Mines, Explosive Objects, and Obscured Targets XIVApril 2009112Proceedings of SPIE2-s2.0-6964908494110.1117/12.817946HaskettH. T.BroachJ. T.Automatic mine detection algorithm using ground penetration radar signaturesDetection and Remediation Technologies for Mines and Minelike Targets IVApril 1999Orlando, Fla, USA942952Proceedings of SPIE2-s2.0-0032676772TantumS. L.WeiY.MunshiV. S.CollinsL. M.A comparison of algorithms for landmine detection and discrimination using ground penetrating radar4742Detection and Remediation Technologies for Mines and Minelike Targets IVApril 2002Orlando, Fla, USA728735Proceedings of SPIE2-s2.0-003644004910.1117/12.479146KovalenkoV.SavelyevT. G.Van KempenL.SahliH.SachsJ.SatoM.Investigation of time-frequency features for GPR landmine discriminationDe MaesschalckR.Jouan-RimbaudD.MassartD. L.The Mahalanobis distance