Next Article in Journal
High Temporal Resolution Monitoring of Suspended Matter Changes from GOCI Measurements in Lake Taihu
Previous Article in Journal
Remote Sensing of Pigment Content at a Leaf Scale: Comparison among Some Specular Removal and Specular Resistance Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GPR Antipersonnel Mine Detection Based on Tensor Robust Principal Analysis

1
College of Electronic Science and Technology, National University of Defense Technology, Changsha 410073, China
2
National Innovation Institute of Technology, Academy of Science, Beijing 100166, China
*
Author to whom correspondence should be addressed.
Submission received: 21 February 2019 / Revised: 15 April 2019 / Accepted: 17 April 2019 / Published: 24 April 2019

Abstract

:
The ground Penetrating Radar (GPR) is a promising remote sensing modality for Antipersonnel Mine (APM) detection. However, detection of the buried APMs are impaired by strong clutter, especially the reflection caused by rough ground surfaces. In this paper, we propose a novel clutter suppression method taking advantage of the low-rank and sparse structure in multidimensional data, based on which an efficient target detection can be accomplished. We firstly created a multidimensional image tensor using sub-band GPR images that are computed from the band-pass filtered GPR signals, such that differences of the target response between sub-bands can be captured. Then, exploiting the low-rank and sparse property of the image tensor, we use the recently proposed Tensor Robust Principal Analysis to remove clutter by decomposing the image tensor into three components: a low-rank component containing clutter, a sparse component capturing target response, and noise. Finally, target detection is accomplished by applying thresholds to the extracted target image. Numerical simulations and experiments with different GPR systems are conducted. The results show that the proposed method effectively improves signal-to-clutter ratio by more than 20 dB and yields satisfactory results with high probability of detection and low false alarm rates.

1. Introduction

Being capable of sensing discontinuity of electromagnetic properties in the ground, Ground Penetrating Radar (GPR) has been used as nondestructive remote sensing modality in many applications, including tree root detection [1], frozen layer monitoring [2], structural assessment [3], infrastructure inspection [4], and landmine detection [5,6,7,8,9]. However, GPR detection of Antipersonnel Mines (APMs) is always challenging because of clutter contamination, such as antenna coupling, ground reflection, and friendly objects (rocks, voids, etc.) [10,11]. As a major source of clutter, the ground reflection is of large magnitude, and usually varies with position due to surface roughness [11]. Since many APMs are of low-metal content and shallow burial, the responses of APMs are often obscured by ground reflections, rendering target visualization and detection very difficult.
Many efforts have been devoted to developing clutter mitigation methods, and these methods can be briefly categorized into three types: model-based methods, filtering, and subspace projection [11,12,13].
In the model-base methods, clutter models are defined, and objects whose responses significantly deviate from the clutter models are detected as potential targets. These methods vary depending on the clutter model chosen to capture the characteristics of the clutter. In [14], a Wide-Sense Stationary (WSS) random process is adopted to describe the clutter whereas in [15], the background clutter is estimated using the least squares technique. A linear prediction model with adaptive coefficients is proposed to characterize the clutter in [5]. However, the model-based methods need “clean” clutter (received signal that contains clutter only) to estimate model parameters which tends to be contaminated by other factors. In addition, the performances of these methods depend strongly on the clutter models and degrade when clutter models deviate from the real situations.
The filtering methods adopt filers in a certain domain to filter out clutter while preserving the target response. An intuitive technique termed time-gating is to mute the received signal within a time window [16]. The problem is how to set a proper time window that covers clutter without compromising target reflections. In [13], the entropy-based criterion for time window selection is proposed under the assumption that clutter has similar behavior for each trace. Another popular method is mean subtraction, which can be viewed as a high-pass filter in the frequency domain, but the calculation of average background clutter trace can be affected by the target responses [16]. A physically-based filtering method is presented in [11] which uses frequency-dependent linear transfer functions to eliminate the strong ground clutter as well as the antenna coupling by subtraction of a simulated Green function. To obtain the linear transfer functions, calibration under controlled environments is required.
The subspace projection method assumes that the ground clutter and the target responses are in different subspaces. Therefore, clutter and target responses can be separated by projecting the recorded GPR data into clutter and target subspaces, respectively. The Principal Component Analysis (PCA) is a representative subspace projection method [17,18]. Usually, the subspace spanned by the first dominant eigencomponent is considered as a clutter subspace to be removed. In practice, the uneven ground surface results in nonstationary clutter whose rank is greater than one and the PCA method cannot remove the clutter completely. Moreover, the PCA method is so sensitive to strong noise that even a single grossly corrupted observation can lead to error results [19]. As an improvement of the traditional PCA, the Robust Principle Component Analysis (RPCA) is proposed for recovering a low-rank subspace from noisy measurements [19]. For GPR data traces, the ground clutter as well as antenna coupling are considered as low-rank, while the target responses are assumed to be contained in the sparse component [20]. Hence, the clutter suppression problem is cast as a low-rank and sparse representation problem and can be solved by the RPCA-based method [12,21]. The recently proposed Tensor Robust Principal Component Analysis (TRPCA) [22,23] provides an alternative perspective of solving the low-rank and sparse decomposition problem. Extended from the RPCA, the TRPCA operates on tensors to seek a low-rank tensor embedded in sparsely corrupted multidimensional data. Avoiding the information loss caused by restructuring high order tensor into matrix, TRPCA makes full use of the abundant information contained in multidimensional data to guarantee flexibility as well as accuracy in video processing [24], seismic data denoising [25], and infrared target detection [26].
Based on TRPCA, this paper presents a GPR APM detection method utilizing the low-rank and sparse properties of multidimensional data for the separation of target response and background clutter. To solve the clutter suppression and target detection problem under tensor algebra framework, we first propose a tensor formation scheme using a sub-band technique. Inspired by the observation that targets give distinct responses at different frequency bands [5,27,28], we exploit sub-band images to characterize the frequency dependent features of APMs. Specifically, the whole frequency band of the recorded GPR signal is divided into several sub-bands to compute individual GPR image for each sub-band and then an image tensor can be formed with these sub-band images inserted as frontal slices (definition in Table 1). By this means, the spatial-spectral information is incorporated in the images tensor, which can be simultaneously handled by a tensor based method to facilitate the detection. Then, TRPCA is used to decompose the obtained image tensor into low-rank component and sparse component to separate the target response (sparse) from the clutter (low-rank). Since most unwanted factors are eliminated, the final detection result can be acquired by applying thresholds on the target image. The proposed method is an extension of the RPCA-based clutter suppression method [21] and the contributions are threefold: (1) a tensor formation scheme is proposed to enable the spatial-spectral processing with TRPCA, and as far as we know, the application of TRPCA in GPR APM detection has not yet been reported; (2) TRPCA is adapted to the APM detection scenario and the problem is solved by a tractable optimization—the Alternating Direction Method of Multiplier (ADMM); (3) our proposed method achieves high Probability of Detection (PD) while effectively reducing the False Alarm Rate (FAR). Numerical simulations and laboratory experiments with impulse GPR and Stepped Frequency Continuous Wave (SFCW) GPR were conducted to validate our proposed method. Both of low-metal APMs and nonmetallic APMs were used as targets and were successfully detected in the experiments. It is noteworthy that the goal of detection is to flag areas that may contain anomalies [29]. Therefore, the subsurface (and surface) anomalies, including APMs and friendly objects, are considered as suspected targets, and the clutter refers to antenna coupling and ground reflection in this paper. The rejection of friendly objects is to be done by a classifier, which is topic of our future work.
The paper is organized as follows. Section 2 introduces the notations and briefly reviews the TRPCA method. The mathematical formulation and the proposed TRPCA-based clutter suppression and target detection method are illustrated in Section 3. Section 4 presents the experiment results as well as discussions of the results to demonstrate the effectiveness of our proposed method. Conclusions is drawn in Section 5.

2. Notations and Preliminaries

2.1. Notations

The notations used in this paper are listed in the Table 1. Below are some necessary definitions.
For a tensor A R n 1 × n 2 × n 3 , we define a block diagonal matrix A ¯ R n 1 n 3 × n 2 n 3 whose ith block on diagonal is the ith frontal slice A ¯ ( i ) of A ¯
A ¯ = A ¯ ( 1 ) A ¯ ( 2 ) A ¯ ( n 3 ) .
A block circular matrix l bcirc A of size n 1 n 3 × n 2 n 3 is denoted as
bcirc ( A ) = A ( 1 ) A ( 2 ) A ( n 3 ) A ( 2 ) A ( 1 ) A ( n 3 ) A ( n 3 ) A ( n 3 1 ) A ( 1 ) .
Two operators, unfold and fold are respectively defined as
unfold ( A ) = A ( 1 ) A ( 2 ) A ( n 3 )
fold ( unfold ( A ) ) = A .
Then, we have a new tensor-tensor production termed as t-product.
Definition 1
(t-product [22]). Let A R n 1 × n 2 × n 3 and B R n 2 × n 4 × n 3 . Then the t-product C = A B is a tensor of size n 1 × n 4 × n 3 defined as
C = A B = fold bcirc A · unfold B .
Using the t-product, a tensor factorization technique named tensor-SVD (t-SVD) can be defined.
Definition 2
(t-SVD [30]). For A R n 1 × n 2 × n 3 , the t-SVD of A is given by
A = U S V T
where U R n 1 × n 1 × n 3 and V R n 2 × n 2 × n 3 are orthogonal tensors, and S is a f-diagonal tensor, T denotes the transpose of a tensor. The detail descriptions of orthogonal tensor, f-diagonal tensor, and tensor transpose are given in [22,30].
Based on t-SVD, the tensor tubal rank as well as the Tensor Nuclear Norm (TNN) is defined.
Definition 3
(tensor tubal rank [22]). The tensor tubal rank of A R n 1 × n 2 × n 3 , denoted as rank t A is defined as the number of nonzero singular tubes of S , where S is from the t-SVD of A = U S V T
r a n k t A = # { i : S ( i , i , : ) 0 } .
Definition 4
(TNN [23]). The tensor nuclear norm of tensor A R n 1 × n 2 × n 3 , denoted as A , is defined as the average of the nuclear norm of all the frontal slices of A ¯
A = 1 n 3 i = 1 n 3 A ¯ ( i ) .
We can further define the tensor spectral norm as
Definition 5
(tensor spectral norm [23]). The tensor spectral norm of tensor A R n 1 × n 2 × n 3 , denoted as A 2 , is defined as the spectral norm of the block diagonal matrix A ¯
A 2 = A ¯ 2 .

2.2. Principle of the TRPCA

Extended from the matrix RPCA, TRPCA considers a similar problem that is to recover a low-rank component corrupted by sparse errors. Different from the RPCA, which handles with matrices, TRPCA manipulates tensor data to take the advantage of abundant information contained in the multidimensional data structure [23]. More specifically, TRPCA aims to separate the low tubal rank tensor L and the sparse component S from the observation X
X = L + S .
Based on t-SVD, the above mentioned problem can be addressed by solving the following convex optimization, where the tensor tubal rank is replaced by the TNN [23]
min L , S L + λ S 1 , s . t . X = L + S .
However, the solution adopted in [23] is aiming at the exactly recovery of tensor. TRPCA need to be adapted to the GPR APM detection application and take the ubiquitous noise in practical measurements into consideration.
Note that when the tensor reduces to matrix, i.e., the third dimension is one, the t-product reduces to standard matrix product and the tensor nuclear norm reduces to matrix nuclear norm. Therefore, RPCA is a special case of TRPCA [23]. The consistency between TRPCA and RPCA helps a lot in our GPR APM detection application, which is to be shown in the later sections.

3. Proposed Method

The proposed ground clutter suppression and target detection method includes three procedures: tensor construction, clutter suppression, and target detection. To begin with, the mathematical formulation of GPR APM detection problem under the tensor framework will be introduced.

3.1. Problem Formulation and Tensor Construction

The scenario of GPR detection is depicted in Figure 1. Emitting and receiving signals at each location, the down-looking GPR scans along a line over the ground surface. The recorded data traces are collected as columns to form a two dimensional matrix (B-scan), ldenoted by X
X = X 1 + X 2 + X 3
where X 1 , X 2 , and X 3 are the clutter, the target response, and noise, respectively. As mentioned above, the clutter that consist of ground reflection and antenna corss-talk is considered as low-rank component in the B-scan, whereas the target response are contained in sparse component [12,21,31]. Therefore, the low-rank and sparse property of GPR image data is utilized to discriminate the clutter X 1 and the target response X 2 by RPCA decomposition.
To further improve the performance of GPR APM detection, we consider incorporating more information contained in the spectrum by constructing an image tensor with multiple sub-band images. As reported in [5,27,28], the sub-band processing takes the distinct features of APMs at different frequency bands into account, thereby increasing detection accuracy. To capture the frequency-dependent characteristics of target responses, sub-band images are formed. Specifically, as shown in Figure 2, the recorded GPR data X is first filtered by two band-pass filters, which divide the whole frequency band into the Low-Frequency (LF) band and the High-Frequency (HF) band while permitting overlap. Two sub-band GPR images, i.e., X L F and X H F , are formed by focusing the LF and the HF GPR data, respectively. Then, we stack the sub-band images as frontal slices, that is, X ( : , : , 1 ) = X L F and X ( : , : , 2 ) = X H F , to create an image tensor resembling a hyperspectral image cube [32], where the spatial characteristics lies in the first two dimensions and the spectral information is contained in the third dimension. According to the low-rank and sparse prior of GPR data [12,21], the obtained GPR image tensor X can be represented as
X = X 1 + X 2 + X 3
where the X 1 , X 2 , and X 3 represent the low-rank clutter tensor, the sparse target image tensor, and noise, respectively. Therefore, the clutter suppression and target detection is cast as a TRPCA problem to separate the sparse target from the low-rank background clutter.

3.2. Clutter Suppression via TRPCA

In this section, clutter is removed by solving the aforementioned TRPCA problem Equation (13) via a convex optimization program. Mathematically,
min X 1 , X 2 , X 3 X 1 1 + λ 2 X 2 + λ 3 X 3 F 2 s . t . X = X 1 + X 2 + X 3 .
The constrained problem Equation (14) can be addressed by the augmented Lagrangian function [33]
L ( X 1 , X 2 , X 3 , U ; ρ , λ 2 , λ 3 ) = i = 1 N ψ i ( X i ) + ρ 2 X i = 1 N X i + U F 2
where i = 1 , , N , N = 3 , ρ > 0 , and functions ψ i , i = 1 , , 3 are defined as
ψ 1 ( X ) = X
ψ 2 ( X ) = λ 2 X 1
ψ 3 ( X ) = λ 3 X F 2
where λ 2 > 0 and λ 3 > 0 are tuning parameters.
We use the ADMM [33] to solve the problem Equation (15) by alternately solving the ( X 1 , X 2 , X 3 , U )
X i k + 1 = argmin X i ψ i ( X i ) + ρ 2 X i D i k F 2
U k + 1 = U k + X ˜ k + 1 1 N X
where i = 1 , , N , N = 3 , the average of tensors X ˜ k = ( 1 / N ) i = 1 N X i k , and D i k is defined as
D i k = X i k X ˜ k + 1 N X U k .
To update the X 1 k + 1 , we use the tensor Singular Value Thresholding (t-SVT) [34]
X 1 k + 1 = F 1 ρ ( D 1 k )
where the t-SVT operator F 1 / ρ is defined as follows: Let X = U S V T be the t-SVD of X R n 1 × n 2 × n 3 . For each 1 / ρ > 0 , the t-SVT on X is
F 1 ρ ( X ) = U S 1 ρ V T
where
S 1 ρ = ifft ( S ¯ 1 ρ ) + , [ ] , 3 .
The subscript ( x ) + denotes the positive part of x.
The X 2 k + 1 is updated using the elementwise soft-shrinkage operator [33]
X 2 k + 1 = shrink D 2 k , λ 2 ρ
where the shrink ( x , λ 2 / ρ ) is
shrink x , λ 2 ρ = x λ 2 ρ x > λ 2 ρ 0 | x | λ 2 ρ x + λ 2 ρ x < λ 2 ρ .
For X 3 k + 1 , we have the closed form solution as [33]
X 3 k + 1 = 1 1 + λ 3 / ρ D 3 k .
Finally, the solving procedures are described in Algorithm 1.
Algorithm 1 The ADMM for TRPCA
Input:
   Image tensor X
Output:
   Target image tensor X 2
Initialize:
   X 1 0 = X 2 0 = X 3 0 = 0 , U 0 = 0 , ρ = 1 , ϵ = 1 e 6 , m = 200 , λ 2 = 0.15 X 2 , λ 3 = 0.15 X
1:
while not converged do
2:
    X 1 k + 1 F 1 / ρ ( D 1 k )
3:
    X 2 k + 1 shrink D 2 k , λ 2 / ρ
4:
    X 3 k + 1 D 3 k / ( 1 + λ 3 / ρ )
5:
    U k + 1 U k + X ˜ k + 1 X / N
6:
     c 1 X 2 k + 1 X 2 k
     c 2 X 3 k + 1 X 3 k
     c 3 U k + 1 U k
7:
   if max ( c 1 , c 2 , c 3 ) < ϵ or k > m then
8:
     Break
9:
   else
10:
      k k + 1
11:
   end if
12:
end while

3.3. Target Detection

After the TRPCA decomposition, a clear target image tensor with most of the clutter removed is obtained, yet, there is residual clutter that may raise false alarms. To facilitate the detection, we sum up the frontal slices of the target image tensor to obtain a target image, that is, X 2 = i = 1 2 X 2 ( : , : , i ) .
With a clear target image, detection can be accomplished as follows. First, an energy detector is applied to the target image. Top α 1 percent of pixels are retained. Then, remove isolated detections whose areas are smaller than α 2 (pixels). At last, cluster the detections and recalculate the center of mass. Those closely spaced detections with distances between each other that are smaller than α 3 (centimeters) are clustered.
As shown in Figure 2, the procedures of our method can be summarized as follows:
  • Tensor construction: Equally divide the bandwidth into two sub-bands and form a sub-band GPR image on each sub-band data using the Layered Range Migration (LRM) method [35]. Create the three-dimensional image tensor by inserting the sub-band images as frontal slices.
  • TRPCA: Decompose the formed image tensor to separate the low-rank component X 1 and the sparse component X 2 . The target response is contained in the sparse component X 2 .
  • Target detection: Calculate the target image X 2 by summing up the frontal slices of X 2 . Remove clutter residue with thresholds and finally locate the target. The center coordinates of the remaining nonzero regions in the target image are output as the detection results.

4. Results and Discussion

To coroborate the proposed clutter suppression and target detection method, both numerical simulations and experiments using different GPR systems were carried out. The simulation data were synthesized using the gprMax [36]. In the experiments, two datasets were collected and tested: one was acquired with an impulse GPR developed by our research group; the other was collected by a SFCW GPR and published by the Georgia Technology Institute [37,38].

4.1. Metrics and Baselines

The Signal to Clutter Ratio (SCR) is adopted to assess the effectiveness of clutter suppression processing. The SCR is calculated as
SCR = N c N t p R t | I ( p ) | 2 p R c | I ( p ) | 2
where I ( p ) is the p-th pixel in image. N t and N c denotes the number of pixels in the target region R t and the number of pixels in the clutter region R c , respectively. The target region is indicated by a box (or boxes in a scenario with multiple targets) enclosing the location of target. The size of the box is manually set to cover the target image. The entire image excluding the target region is defined as the clutter region.
The metrics of evaluating detection performance are PD P d and FAR F a , which are defined as
P d = number   of   correct   detections number   of   actual   targets
F a = number   of   false   alarms total   length   of   the   surveyed   scenario
Note that F a is calculated on a meter basis. For example, a 1% FAR means that, on average, there is one false alarm per 100 m.
For comparison purposes, two representative clutter suppression and target detection methods, the PCA-CFAR method and the RPCA-based method, are adopted as baseline methods. The PCA-CFAR method utilizes the PCA [17] to remove clutter from a GPR image by discarding the largest principal component, and then detects target by Cell Averaging Constant FAR (CA-CFAR) technique [39]. The RPCA-based detection method [40] first extracts a target image from clutter using the Go-Decomposition (GoDec) [41] and then applies thresholds to reveal the target location. Same thresholds α 1 , α 2 , and α 3 are used here.

4.2. Numerical Simulations

Shown as Figure 1, point source and receiver spaced 2 cm apart were 15 cm distant from the ground and scanned with 1 cm step length. A Ricker waveform with 2 GHz center frequency and 1.6 GHz bandwidth were adopted. The time window is set to 3.11 ns and the sampling interval is 2.59 ps. The soil is considered homogeneous, and the relative permittivity and conductivity of soil are ϵ = 5 and σ = 0.01 S/m, respectively. As shown in Figure 1, the APM was modeled as a hollow cylinder containing another cylinder inside. The outer cylinder was of radius r 1 = 3.5 cm and height h 1 = 4.5 cm, while the inner cylinder was of r 2 = 2.5 cm radius and h 2 = 3 cm height. The relative permittivity of the outer and the inner cylinders were ϵ 1 = 2.9 and ϵ 2 = 3.6 , respectively.
One thousand B-scans were synthesized with the parameters listed above. Each B-scan is of the same size with 60 columns and 1198 rows, that is, it covers a length of 60 cm. The rough ground surfaces in these simulations were randomly generated using the fractal correlated noise model [42]. The fractal dimension was set to 1.5 and variation in height was 6 cm. In addition, Gaussian noise was added to the data and thus the Signal to Noise Ratio (SNR) was reduced to 0 dB.
Figure 3a displays one of the simulated noise-free B-scans where an APM is flush buried. Therefore, a part of ground reflections coincides with the target responses in time domain and results in waveform distortion. After noise is added, the target responses in Figure 3b are overwhelmed and hardly any hyperbola signature of the target are extracted. Using the LRM to focus the GPR data yields a full band imaging result (Figure 3c), where the target size and location are indicated by a dash-dotted box. Due to obscuration of the dominant ground clutter, the target image cannot be observed. Similar to [27], we equally divide the frequency bandwidth into two sub-bands while allowing 50% overlap: the LF band of 0.1–3.8 GHz and the HF band of 1.9–5.6 GHz. The optimal selection of the overlapping range will be discussed at the end of this section. As shown Figure 3d,e, two sub-band images, i.e., the LF image and the HF image are computed using these sub-band GPR data. Note that the uneven surface has similar shapes in these images, whereas the obscuration caused by the ground clutter renders target detection difficult. Then, we create an image tensor with these sub-band images and decompose the image tensor with TRPCA to obtain the target image tensor. In order to display the decomposition results, frontal slices of the obtained target image tensor are summed up to form a two-dimensional image, as shown in Figure 3f. The target image in Figure 3f is clear enough to recognize the APM. As listed in Table 2, the SCR of Figure 3f is improved by 24 dB compared with the full band image shown in Figure 3c. Since the major part of the clutter is removed, we can use thresholds to remove clutter residues and get the final detection result. As shown in Figure 3g, the target is indicated by a point.
The reasons why TRPCA performs well in the environment with heavy clutter are due to two aspects: (1) the low-rank and sparse properties of the GPR images, and (2) the frequency dependency of target response. Since the model and theory of TRPCA are consistent with RPCA [23], the low-rank and sparse structure of GPR image can be utilized by the TRPCA for the separation of target image and clutter in each sub-band image (frontal slice), which is similar to the RPCA-based methods [12,21]. The proposed method uses imaging to enhance the sparsity and intensity of a target response, which generally leads to a better decomposition [19]. Moreover, our method takes it one step further than the RPCA-based method by making use of the spectral information of target response lies in the third dimension. Though the differences caused by frequency dependency of target response are hard to discern by naked eyes in the sub-band images (Figure 3d,e), they can be captured and revealed by the TRPCA, which is intended for extraction of the variations in frontal slices [22].
A different situation is shown in Figure 4a where an APM is buried deeper (7 cm beneath the ground) due to the raised surface. The ground clutter and target responses are separated in time, but two pits in the ground at x = 14 cm and x = 42 cm causes strong reflections with hyperbolic curve signatures. After imaging, the target images can be found in the centers of sub-band images (Figure 4d,e). However, the clutter, especially images of the pits are of large magnitude that can raise false alarms. To remove the ground clutter and avoid false alarms, our proposed method is applied to improve SCR of the target image in Figure 4f by 24.8 dB, compared to the full band image Figure 4c. The target location is marked by a point in Figure 4g after the threshold-based detection procedure.
Processing all of the one thousand simulated B-scans with our proposed method, the detection performance in the low SNR (0 dB) environment is of 99.5% PD and 10% FAR, as shown in Table 3.
The PCA-CFAR method and the recently proposed RPCA-based method were tested on the same synthesized dataset. The full band images (Figure 3c and Figure 4c) were used as the inputs for these two methods. Figure 3h shows the detection result of the PCA-CFAR method in the shallow-burying scenario, where the target and false alarms are indicated with point and triangles, respectively. Due to the uneven ground surface, the rank of clutter is bigger than one and thus cannot be fully removed because only the first principal component is removed. In addition, the traditional PCA is so sensitive to strong noise that the processed image is contaminated by remaining noise. Therefore, together with the target marked by point, there is residual clutter and noise which degrades SCR and causes false alarms. This is similar to the case with deeply buried targets; there is clutter and noise residue in Figure 4h due to the deficiencies of PCA-CFAR. The improvements of SCR of the PCA-CFAR method in these two scenarios are listed in Table 2, which are 6.8 dB and 4.7 dB, respectively. The detection performance on the whole dataset is low PD (87.0%) with a high FAR of 67.0%. It shows that the PCA-CFAR method cannot handle the clutter caused by an uneven ground surface very well and its performance degrades severely under low SNR.
Thanks to the low-rank and sparse signal model, the RPCA-based method outperforms the PCA-CFAR method in the simulation by yielding higher SCRs (20.3 dB in Figure 3i, 21 dB in Figure 4i) and better detection performance (93.6% PD and 27.7% FAR). However, the variation of surface height results in nonstationary ground clutter whose low-rank property is degraded and thus, the RPCA-based method performs inadequately. This is the reason why both bulge and cavity in the ground cause false alarms in detection results Figure 3i and Figure 4i. Thus, in the numerical simulation, the detection performance of RPCA-based method is inferior to that of the proposed TRPCA-based method.
The time consumptions of these three methods are listed in Table 4. The solving programs of these methods are running on a same computer with 2.4 GHz CPU and 8 GB memory. It takes 0.24 s for the proposed method to complete the detection on a single B-scan, whereas the PCA-CFAR method takes 3.7 s. By using a accelerated solution program [41], time consumption of the RPCA-based method is 0.14 s. Note that the aforementioned optimization program ADMM can be further improved by some techniques [32,43,44]. However, these discussions are out of the scope of this paper.

4.3. Laboratory Experiments

Laboratory experiments using different GPR systems were conducted to demonstrate the effectiveness of our proposed method on real data.
In the first experiment, an impulse GPR of 1.6 GHz center frequency and 1.2 GHz bandwidth was used to scan a sand pit. The horn antennas were 15 cm above the surface. To deal successfully with the strong ground reflections, a dynamic range of over 69 dB is expected [45]. The space sampling step length was 5 mm and the time sampling interval was 9.77 ps. As Figure 5 shows, the pit was filled with dry sand, whose relative permittivity was 4. The surface of sand was flat. Seventy-two B-scans containing 90 APMs were collected. Each B-scan covered a length of 0.9 m, where only one APM or at most two APMs were buried. Both of metal APMs and nonmetallic APMs were used as targets in the experiment. We take two B-scans out from the dataset as examples to show the feasibility of our method.
As shown in Figure 6a, one of the B-scans contains a plastic APM buried 5 cm under the ground. Indicated by a dash-dotted box, the APM is of 4 cm height and 3.5 cm radius. We can hardly notice the subtle target responses around 2.0 ns because the direct waves and ground reflections (overlapped between 0.2 ns and 1.8 ns) are dominant in the B-scan (Figure 6a). The imaging result of the B-scan is displayed in Figure 6b where an ambiguous dot appears at 15 cm location. Though focused and enhanced by the imaging processing, the images of target in sub-band images Figure 6c,d are still obscured by the clutter. Processing the sub-band images with our method yields much more clearer target image (Figure 6e) that are brighter and larger than the remaining clutter. Comparing Figure 6e and Figure 6b, the SCR is improved by 26.3 dB. After applying thresholds, clutter residues are eliminated and the target is detected, as shown in Figure 6f. The results of competitive methods are shown in Figure 6g,h. Since the target is small and contains no metal, the responses are too weak to be detected by the PCA-CFAR method and the RPCA-based method. In addition, clutter residues in Figure 6g cause multiple false alarm in the detection results.
In the other B-scan, there are two shallowly buried (3 cm) APMs containing both metal and nonmetal. These APMs are identical and their heights and radiuses are 5.4 cm and 6 cm, respectively. The size and location of the APMs are indicated by dash-dotted boxes. The reflections from metal are strong enough to be observed in Figure 7a. However, the target image in Figure 7b is overlapped with strong ground clutter and the SCR is insufficient to confirm the target. Similar are the sub-band images shown in Figure 7c,d, where the target images are of small magnitude compared with the ground reflections. Note that more details of the target can be found in the HF image Figure 7d than in LF image Figure 7c, which demonstrates the frequency dependency of target responses in different sub-bands. Applying our method to this GPR data, we obtain the prominent target image with a much higher SCR (30 dB) in Figure 7e. Then, the targets’ locations are successfully detected and marked by points in Figure 7f. Figure 7g,h display the results of the PCA-CFAR method and the RPCA-based method, and the SCR of these images are 12.6 dB and 19.8 dB, respectively.
As listed in Table 3, the performance of our method on this dataset outperforms the PCA-CFAR with more accurate detection (100% PD) and less false alarms (12.3% FAR), whereas the PCA-CFAR method yields a higher FAR (54.5%) to achieve 89.4% PD. The performance of RPCA-based method is of 93.3% PD and 10.8% FAR.
The second experiment was conducted on the dataset published by the Georgia Technology Institute. Seven types of mines were buried in a 1.8 m × 1.8 m region. The test site was different from the actual minefield because the targets were close to each other, but this dataset can still be used to demonstrate the applicability of our method to the SFCW GPR data. The experimental setup is detailed in [37,38]. Fifty B-scans acquired along the x axis were selected to test our method. As shown in the burial map Figure 8a, the coverage area of these B-scans is enclosed by a rectangle, and a dash line across the center indicates location of the example B-scan. The locations and sizes of these four target are indicated by four dash-dotted boxes. As we can see in the imaging result (Figure 8b), the horizontal stripe-like clutter is overlapped with the images of targets. Two deep buried targets can be found due to their big sizes and strong reflections while the other two small targets buried shallowly near the center are hidden in the clutter. In the LF image Figure 8c, two targets of the bigger sizes can be noticed, while the smallest simulant mine is prominent in the HF image Figure 8d. Removing clutter using our method, the result is shown in Figure 8e, where the SCR is increased to 26.8 dB. Finally, all of the four targets are correctly detected and indicated with points in Figure 8f after the thresholds are applied. After processing the data with the PCA-CFAR method and the RPCA-based method, SCRs are improved by 11 dB and 27.9 dB, respectively. Also, the targets are detected by these two methods, as shown in Figure 8g,h. However, there is a false alarm indicated by a triangle in Figure 8h, which may be caused by the residual ground clutter.
After processing this dataset, our proposed method detects all of the targets shown in Figure 8a with an FAR of 7.0%, while the PCA-CFAR miss two small targets and the FAR is 19.2%. The RPCA-based method achieves 100% PD and 8.4% FAR. The performances are summarized in Table 3.

4.4. Discussion

  • Reconstruction of the target image. As a part of our proposed method, GPR data is migrated using the LRM algorithm. The resolution of the signal is not enough for an accurate shape reconstruction of our mine-like targets but improves the SCR of the image, making the targets look like points with high energy.
  • Overlap of sub-band. The idea of improving performance by allowing some overlap comes form [27]. According to our investigation on the synthesized dataset, the best performance is obtained with an overlap close to 50–60% (see Figure 9).

5. Conclusions

In this paper, we propose a novel TRPCA-based method for clutter suppression and target detection. A tensor formation scheme capturing the target sub-band features is firstly presented. Then, taking advantage of the low-rank and sparse property of the multidimensional tensor, our method formulates the clutter removal and target detection as a TRPCA problem and then solves the problem by decomposing the obtained tensor with a tractable method—ADMM. The feasibility and effectiveness of our proposed method for APM detection applications are verified by both numerical simulations and laboratory experiments with different GPR systems. Our proposed method is distinguished from traditional ones by the advantage of multidimensional data structure, which enables exploiting of the information contained in different sub-bands. Since the sub-band technique is a rough description of the frequency-dependent features, our future works will focus on designing a precise characterization of target features to guide the tensor construction and applying classification to discard false positives. Considering the conditions of laboratory experiments are different from those of the actual minefields, the analysis should be extended to field experiments under realistic clutter conditions.

Author Contributions

X.S. proposed the original idea and wrote the manuscript; T.L. designed and performed the experiments; X.S., T.L. and D.X. analyzed the data; Y.S. contributed to the direction, content, and revised the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhu, S.; Huang, C.; Su, Y.; Sato, M. 3D Ground Penetrating Radar to Detect Tree Roots and Estimate Root Biomass in the Field. Remote Sens. 2014, 6, 5754–5773. [Google Scholar] [CrossRef] [Green Version]
  2. Jadoon, K.; Weihermüller, L.; McCabe, M.; Moghadas, D.; Vereecken, H.; Lambot, S. Temporal Monitoring of the Soil Freeze-Thaw Cycles over a Snow-Covered Surface by Using Air-Launched Ground-Penetrating Radar. Remote Sens. 2015, 7, 12041–12056. [Google Scholar] [CrossRef] [Green Version]
  3. Catapano, I.; Ludeno, G.; Soldovieri, F.; Tosti, F.; Padeletti, G. Structural Assessment via Ground Penetrating Radar at the Consoli Palace of Gubbio (Italy). Remote Sens. 2017, 10, 45. [Google Scholar] [CrossRef]
  4. Dong, Z.; Ye, S.; Gao, Y.; Fang, G.; Zhang, X.; Xue, Z.; Zhang, T. Rapid Detection Methods for Asphalt Pavement Thicknesses and Defects by a Vehicle-Mounted Ground Penetrating Radar (GPR) System. Sensors 2016, 16, 2067. [Google Scholar] [CrossRef]
  5. Ho, K.C.; Gader, P.D. A linear prediction land mine detection algorithm for hand held ground penetrating radar. IEEE Trans. Geosci. Remote Sens. 2002, 40, 1374–1384. [Google Scholar] [CrossRef]
  6. Slob, E.; Sato, M.; Olhoeft, G. Surface and borehole ground-penetrating-radar developments. Geophysics 2010, 75, A103–A120. [Google Scholar] [CrossRef]
  7. Feng, X.; Sato, M.; Liu, C.; Takahashi, K.; Zhang, Y. Topographic Correction of Elevated GPR. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 799–804. [Google Scholar] [CrossRef]
  8. Núñez-Nieto, X.; Solla, M.; Gómez-Pérez, P.; Lorenzo, H. GPR Signal Characterization for Automated Landmine and UXO Detection Based on Machine Learning Techniques. Remote Sens. 2014, 6, 9729–9748. [Google Scholar] [CrossRef] [Green Version]
  9. Gao, X.; Podd, F.; van Verre, W.; Daniels, D.; Peyton, A. Investigating the Performance of Bi-Static GPR Antennas for Near-Surface Object Detection. Sensors 2019, 19, 170. [Google Scholar] [CrossRef]
  10. Kovalenko, V.; Yarovoy, A.G.; Ligthart, L.P. A novel clutter suppression algorithm for landmine detection with GPR. IEEE Trans. Geosci. Remote Sens. 2007, 45, 3740–3751. [Google Scholar] [CrossRef]
  11. Lopera, O.; Slob, E.C.; Milisavljevic, N.; Lambot, S. Filtering Soil Surface and Antenna Effects From GPR Data to Enhance Landmine Detection. IEEE Trans. Geosci. Remote Sens. 2007, 45, 707–717. [Google Scholar] [CrossRef]
  12. Tivive, F.H.C.; Bouzerdoum, A.; Abeynayake, C. GPR Target Detection by Joint Sparse and Low-Rank Matrix Decomposition. IEEE Trans. Geosci. Remote Sens. 2018. [Google Scholar] [CrossRef]
  13. Solimene, R.; Cuccaro, A.; Dell’Aversano, A.; Catapano, I.; Soldovieri, F. Ground Clutter Removal in GPR Surveys. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 792–798. [Google Scholar] [CrossRef]
  14. Dogaru, T.; Carin, L. Time-Domain Sensing of Targets Buried Under a Rough Air–Ground Interface. IEEE Trans. Antennas Propag. 1998, 46, 360–372. [Google Scholar] [CrossRef]
  15. Brunzell, H. Detection of shallowly buried objects using impulse radar. IEEE Trans. Geosci. Remote Sens. 1999, 37, 875–886. [Google Scholar] [CrossRef]
  16. Daniels, D.J. Ground Penetrating Radar, 2nd ed.; IEE: London, UK, 2004. [Google Scholar]
  17. Karlsen, B.; Larsen, J.; Sorensen, H.B.D.; Jakobsen, K.B. Comparison of PCA and ICA based clutter reduction in GPR systems for anti-personal landmine detection. In Proceedings of the 11th IEEE Signal Processing Workshop Statisttical Signal Processing, Singapore, 8 August 2001; pp. 146–149. [Google Scholar]
  18. Abujarad, F.; Omar, A. GPR Data Processing Using the Component-Separation Methods PCA and ICA. In Proceedings of the 2006 IEEE International Workshop on Imagining Systems and Techniques, Minori, Italy, 29 April 2006; pp. 60–64. [Google Scholar]
  19. Candès, E.J.; Li, X.; Ma, Y.; Wright, J. Robust Principal Component Analysis? J. ACM 2011, 58, 11. [Google Scholar] [CrossRef]
  20. Kalika, D.; Knox, M.T.; Collins, L.M.; Torrione, P.A.; Morton, K.D. Leveraging Robust Principal Component Analysis To Detect Buried Explosive Threats In Handheld Ground-Penetrating Radar Data. In Detection and Sensing of Mines, Explosive Objects, and Obscured Targets XX; Bishop, S.S., Isaacs, J.C., Eds.; Spie-Int Soc Optical Engineering: Bellingham, WA, USA, 2015; Volume 9454, p. 94541D. [Google Scholar]
  21. Song, X.; Xiang, D.; Zhou, K.; Su, Y. Improving RPCA-based clutter suppression in GPR detection of antipersonnel mines. IEEE Geosci. Remote. Sens. 2017, 17, 1338–1442. [Google Scholar] [CrossRef]
  22. Zhang, Z.; Ely, G.; Aeron, S.; Hao, N.; Kilmer, M. Novel Methods for Multilinear Data Completion and De-noising Based on Tensor-SVD. In Proceedings of the 2014 IEEE Conference on Computer Vision and Pattern Recognition, Columbus, OH, USA, 23–28 June 2014; pp. 3842–3849. [Google Scholar]
  23. Lu, C.; Feng, J.; Chen, Y.; Liu, W.; Lin, Z.; Yan, S. Tensor Robust Principal Component Analysis: Exact Recovery of Corrupted Low-Rank Tensors via Convex Optimization. In Proceedings of the 2016 IEEE Conference on Computer Vision and Pattern Recognition, Las Vegas, NV, USA, 27–30 June 2016; pp. 5249–5257. [Google Scholar]
  24. Wei, D.; Wang, A.; Feng, X.; Wang, B.; Wang, B. Tensor Completion Based on Triple Tubal Nuclear Norm. Algorithms 2018, 11, 94. [Google Scholar] [CrossRef]
  25. Ely, G.; Aeron, S.; Hao, N.; Kilmer, M.E. 5D seismic data completion and denoising using a novel class of tensor decompositions. Geophysics 2015, 80, V83–V95. [Google Scholar] [CrossRef] [Green Version]
  26. Zhang, L.; Peng, Z. Infrared Small Target Detection Based on Partial Sum of the Tensor Nuclear Norm. Remote Sens. 2019, 11, 382. [Google Scholar] [CrossRef]
  27. Wang, T.; Keller, J.M.; Gader, P.D.; Sjahputera, O. Frequency sub-band Processing and Feature Analysis of Forward-Looking Ground-Penetrating Radar Signals for Land-Mine Detection. IEEE Trans. Geosci. Remote Sens. 2007, 45, 718–729. [Google Scholar] [CrossRef]
  28. Webb, A.; Havens, T.C.; Schulz, T.J. Spectral diversity for ground clutter mitigation in forward-looking GPR. In Detection and Sensing of Mines, Explosive Objects, and Obscured Targets XXI; Bishop, S.S., Isaacs, J.C., Eds.; SPIE: Bellingham, WA, USA, 2016; Volume 9823, p. 98231M. [Google Scholar]
  29. Torrione, P.A.; Throckmorton, C.S.; Collins, L.M. Performance of an adaptive feature-based processor for a wideband ground penetrating radar system. IEEE Trans. Aerosp. Electron. Syst. 2006, 42, 644–658. [Google Scholar] [CrossRef]
  30. Kilmer, M.E.; Martin, C.D. Factorization strategies for third-order tensors. Linear Algebra Appl. 2011, 435, 641–658. [Google Scholar] [CrossRef] [Green Version]
  31. Masarik, M.P.; Burns, J.; Thelen, B.T.; Kelly, J.; Havens, T.C. GPR anomaly detection with robust principal component analysis. In Detection and Sensing of Mines, Explosive Objects, and Obscured Targets XX; SPIE: Bellingham, WA, USA, 2015; Volume 9454, p. 945414. [Google Scholar]
  32. Fan, H.; Chen, Y.; Guo, Y.; Zhang, H.; Kuang, G. Hyperspectral Image Restoration Using Low-Rank Tensor Recovery. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2017, 10, 4589–4604. [Google Scholar] [CrossRef]
  33. Boyd, S.; Parikh, N.; Chu, E.; Peleato, B.; Eckstein, J. Distributed Optimization and Statistical Learning via the Alternating Direction Method of Multipliers. Found. Trends Mach. Learn. 2011, 3, 1–122. [Google Scholar] [CrossRef]
  34. Lu, C.; Feng, J.; Chen, Y.; Liu, W.; Lin, Z.; Yan, S. Tensor Robust Principal Component Analysis with a New Tensor Nuclear Norm. 2018. Available online: http://arxiv.org/abs/1804.03728v1 (accessed on 1 November 2018).
  35. Zhou, L.; Su, Y. GPR Imaging With RM Algorithm in Layered Mediums. IEEE Geosci. Remote Sens. Lett. 2011, 8, 934–938. [Google Scholar] [CrossRef]
  36. Giannakis, I.; Giannopoulos, A.; Warren, C. A Realistic FDTD Numerical Modeling Framework of Ground Penetrating Radar for Landmine Detection. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 37–51. [Google Scholar] [CrossRef]
  37. Counts, T.; Gurbuz, A.C.; Scott, W.R.; McClellan, J.H.; Kim, K. Multistatic Ground-Penetrating Radar Experiments. IEEE Trans. Geosci. Remote Sens. 2007, 45, 2544–2553. [Google Scholar] [CrossRef]
  38. Scott, W.R. Multistatic Beamforming Data. 2017. Available online: http://waymond-scott.ece.gatech.edu/ multistaticbeamformingdata/multistaticbeamformingdata/ (accessed on 1 December 2017).
  39. Eldarymli, K.; Mcguire, P.; Power, D.; Moloney, C. Target detection in synthetic aperture radar imagery: A state-of-the-art survey. J. Appl. Remote Sens. 2013, 7, 071598. [Google Scholar] [CrossRef]
  40. Song, X.; Xiang, D.; Zhou, K.; Su, Y. Fast Prescreening for GPR Antipersonnel Mine Detection via Go Decomposition. IEEE Geosci. Remote Sens. Lett. 2019, 16, 15–19. [Google Scholar] [CrossRef]
  41. Zhou, T.; Tao, D.; Zhou, T.; Tao, D. GoDec: Randomized Lowrank & Sparse Matrix Decomposition in Noisy Case. In Proceedings of the 28th International Conference on Machine Learning, Bellevue, WA, USA, 28 June–2 July 2011; ACM: Bellevue, WA, USA, 2011; pp. 33–40. [Google Scholar]
  42. Giannakis, I.; Giannopoulos, A.; Yarovoy, A. Model-Based Evaluation of Signal-to-Clutter Ratio for Landmine Detection Using Ground-Penetrating Radar. IEEE Trans. Geosci. Remote Sens. 2016, 54, 3564–3573. [Google Scholar] [CrossRef]
  43. Ghadimi, E.; Teixeira, A.; Shames, I.; Johansson, M. Optimal Parameter Selection for the Alternating Direction Method of Multipliers (ADMM): Quadratic Problems. IEEE Trans. Autom. Control 2015, 60, 644–658. [Google Scholar] [CrossRef] [Green Version]
  44. Chang, T.H.; Hong, M.; Liao, W.C.; Wang, X. Asynchronous Distributed ADMM for Large-Scale Optimization—Part I: Algorithm and Convergence Analysis. IEEE Trans. Signal Process. 2016, 64, 3118–3130. [Google Scholar] [CrossRef]
  45. Yarovoy, A.G.; Ligthart, L.P.; Schukin, A.; Kaploun, I. Polarimetric Video Impulse Radar for Landmine Detection. Subsurf. Sens. Technol. Appl. 2002, 3, 271–293. [Google Scholar] [CrossRef]
Figure 1. Measurement Geometry.
Figure 1. Measurement Geometry.
Remotesensing 11 00984 g001
Figure 2. Flow chart of the TRPCA-based clutter suppression and target detection method.
Figure 2. Flow chart of the TRPCA-based clutter suppression and target detection method.
Remotesensing 11 00984 g002
Figure 3. Results of the first simulation where a target is shallowly buried. (a,b) show the synthesized data without and with noise, respectively. (c) Full band imaging result. (d) LF imaging result. (e) HF imaging result. (f) Extracted target image. (g) Detection result. (h) Result of PCA-CFAR. (i) Result of RPCA-based method.
Figure 3. Results of the first simulation where a target is shallowly buried. (a,b) show the synthesized data without and with noise, respectively. (c) Full band imaging result. (d) LF imaging result. (e) HF imaging result. (f) Extracted target image. (g) Detection result. (h) Result of PCA-CFAR. (i) Result of RPCA-based method.
Remotesensing 11 00984 g003
Figure 4. Results of the second simulation where a target is deep buried. (a,b) show the synthesized data without and with noise, respectively. (c) Full band imaging result. (d) LF imaging result. (e) HF imaging result. (f) Extracted target image. (g) Detection result. (h) Result of PCA-CFAR. (i) Result of RPCA-based method.
Figure 4. Results of the second simulation where a target is deep buried. (a,b) show the synthesized data without and with noise, respectively. (c) Full band imaging result. (d) LF imaging result. (e) HF imaging result. (f) Extracted target image. (g) Detection result. (h) Result of PCA-CFAR. (i) Result of RPCA-based method.
Remotesensing 11 00984 g004
Figure 5. Photograph of measurement setup in the impulse GPR experiment.
Figure 5. Photograph of measurement setup in the impulse GPR experiment.
Remotesensing 11 00984 g005
Figure 6. Results of the first impulse GPR experiment. A plastic APM is buried in the scenario. (a) GPR data. (b) Full band imaging result. (c) LF imaging result. (d) HF imaging result. (e) Extracted target image. (f) Detection result. (g) Result of PCA-CFAR. (h) Result of RPCA-based method.
Figure 6. Results of the first impulse GPR experiment. A plastic APM is buried in the scenario. (a) GPR data. (b) Full band imaging result. (c) LF imaging result. (d) HF imaging result. (e) Extracted target image. (f) Detection result. (g) Result of PCA-CFAR. (h) Result of RPCA-based method.
Remotesensing 11 00984 g006
Figure 7. Results of the second impulse GPR experiment. Two APMs containing both metal and nonmetal are buried in the scenario. (a) GPR data. (b) Full band imaging result. (c) LF imaging result. (d) HF imaging result. (e) Extracted target image. (f) Detection result. (g) Result of PCA-CFAR. (h) Result of RPCA-based method.
Figure 7. Results of the second impulse GPR experiment. Two APMs containing both metal and nonmetal are buried in the scenario. (a) GPR data. (b) Full band imaging result. (c) LF imaging result. (d) HF imaging result. (e) Extracted target image. (f) Detection result. (g) Result of PCA-CFAR. (h) Result of RPCA-based method.
Remotesensing 11 00984 g007aRemotesensing 11 00984 g007b
Figure 8. Results of the SFCW GPR experiment. Four mines of different sizes are buried in the scenario. (a) Burial map (axis in cm) [37]. (b) Full band imaging result. (c) LF imaging result. (d) HF imaging result. (e) Extracted target image. (f) Detection result. (g) Result of PCA-CFAR. (h) Result of RPCA-based method.
Figure 8. Results of the SFCW GPR experiment. Four mines of different sizes are buried in the scenario. (a) Burial map (axis in cm) [37]. (b) Full band imaging result. (c) LF imaging result. (d) HF imaging result. (e) Extracted target image. (f) Detection result. (g) Result of PCA-CFAR. (h) Result of RPCA-based method.
Remotesensing 11 00984 g008
Figure 9. PD/FAR vs overlap.
Figure 9. PD/FAR vs overlap.
Remotesensing 11 00984 g009
Table 1. Metamathematical notations.
Table 1. Metamathematical notations.
NotationsDefinition
A , A , a , a tensor, matrix, vector, scalar
R , C real numbers and complex numbers
A i , j , k or a i , j , k the ( i , j , k ) th entry of A
A ( i , j , : ) the tube of A
A ( i , : , : ) , A ( : , i , : ) the ith horizontal or lateral slice of A
A ( : , : , i ) or A ( i ) the ith frontal slice of A
A ¯ = fft A , [ ] , 3 Fast Fourier Transformation (FFT) on A along the third dimension
A = ifft A ¯ , [ ] , 3 Inverse FFT on A along the third dimension
A 1 1 norm of A , which is computed as i , j , k | a i , j , k |
A infinity norm of A , which is defined as max | a i , j , k |
A F Frobenius norm of A , which is computed as i , j , k | a i , j , k | 2
A nuclear norm of A , which is computed as sum of singular values
A 2 spectral norm of A , which is defined as the largest singular value
Table 2. SCR of different clutter suppression methods.
Table 2. SCR of different clutter suppression methods.
ImaingThe Proposed MethodPCA-CFARRPCA-Based Method
Simulated data 15.629.712.420.3
Simulated data 26.731.511.421.0
Impulse GPR data 1−4.022.30.1
Impulse GPR data 2−0.530.012.619.8
SFCW GPR data−10.626.80.417.3
Table 3. Comparison of detection performance (PD/FAR).
Table 3. Comparison of detection performance (PD/FAR).
The Proposed MethodPCA-CFARRPCA-Based Method
Simulated dataset99.5%/10.0%87.0%/67.0%93.6%/27.7%
Impulse GPR dataset100%/12.3%89.4%/54.5%93.3%/10.8%
SFCW GPR dataset100%/7%84.6%/19.2%100%/8.4%
Table 4. Time consumption of different methods.
Table 4. Time consumption of different methods.
The Proposed MethodPCA-CFARRPCA-Based Method
Time consumption (s)0.243.80.14

Share and Cite

MDPI and ACS Style

Song, X.; Liu, T.; Xiang, D.; Su, Y. GPR Antipersonnel Mine Detection Based on Tensor Robust Principal Analysis. Remote Sens. 2019, 11, 984. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11080984

AMA Style

Song X, Liu T, Xiang D, Su Y. GPR Antipersonnel Mine Detection Based on Tensor Robust Principal Analysis. Remote Sensing. 2019; 11(8):984. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11080984

Chicago/Turabian Style

Song, Xiaoji, Tao Liu, Deliang Xiang, and Yi Su. 2019. "GPR Antipersonnel Mine Detection Based on Tensor Robust Principal Analysis" Remote Sensing 11, no. 8: 984. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11080984

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