|
1.IntroductionGround-penetrating radar (GPR) is a nondestructive method using electromagnetic radiation to locate shallow geological subsurface features and underground utilities buried in the ground. It has become a valuable tool in several applications,1,2 such as archaeological explorations, environmental engineering, and geological problems. The effective imaging of buried objects is a key part of GPR, and the efficiency and resolution of the imaging results are the measure of the imaging algorithm.3 The theories of the present imaging methods are based on diffraction tomography (DT), reverse time migration (RTM), range migration (RM), and back projection (BP). The principle of the DT algorithm is based on the first-order Born approximation which assumes that the buried object of interest is a weak scatterer.4 A few additional assumptions are also invoked during the process of DT derivation to simplify and linearize the nonlinear electric field integral equation. These assumptions incur a trade-off to the reconstruction of the buried objects especially for the practical usage when noise is present in the collected field data. Taking advantage of the multiple reflections in the propagation medium, the RTM algorithm allows high-resolution focusing.5 However, the number of transmitting and receiving antennas must be more than the number of scatterers in the medium. The RM algorithm can work well only when the imaging scene can be modeled as a single homogeneous medium.6 When the GPR antennas and buried objects are distributed in different media, the imaging result of the RM algorithm will be blurred or possibly not focused at all. The standard two-dimensional (2-D) depth migrations7 can recover the location and shape of buried objects with arbitrary precision, depending on the accuracy of the velocity model used. The BP algorithm is one of the most practical imaging methods because of its convenience and robustness, particularly when the imaging scene can be modeled as layered media.8 Based on the aforementioned theories, some of the improvements and optimization imaging methods have been advanced to distinguish the shape of buried objects in GPR imaging. A modified split-step method9 was applied to extract structural information from a complex synthetic data set as accurately as possible, based on the standard 2-D depth migrations. Furthermore, a synthetic aperture radar technique10 was implemented for GPR image reconstruction, which can recover the shape of buried objects. However, the present imaging methods depend too much on the application environment or prior knowledge of the medium being imaged, which limits the popularization and application of GPR technology. The BP theory can accurately compensate the distortion of the wave path caused by GPR pulse signal when GPR pulse signal passes an interface of two media. It has become one of the potential GPR imaging algorithms. In order to recover the shape of buried objects in GPR imaging, the self-correlation back-projection (SBP) algorithm is proposed in this paper. The following improvements are made in the proposed algorithm. First, the valid echo information sequence of the imaging points is adaptively chosen by setting up a correlation threshold. Then, the imaging result is postprocessed by the depth energy compensation algorithm. Finally, the performance of the proposed algorithm is verified through serial contrast experiment. 2.Self-Correlation Back-Projection AlgorithmFor the convenience of discussion, the 2-D imaging configuration of the GPR system is shown in Fig. 1. The scene is divided into two regions by . The upper region is air with relative permittivity and conductivity . The lower region is a homogeneous medium with relative permittivity and conductivity . The GPR system works in a monostatic way. The antennas transmit and receive signals in each of the positions on a survey line . 2.1.Two-Way TraveltimeIn the 2-D imaging configuration shown in Fig. 1, the current concerned antenna position is represented by a black cube with sequence number , the coordinates of which are . For a given point with coordinates in the imaging area, the transmitting signal travels from to , with a turning at the inflection point , and returns along the same path in reverse direction. According to Snell’s law, the geometry between the incidence angle and the refraction angle satisfies where is the velocity of the electromagnetic wave in free space. Equation (1) can be turned into Eq. (2) based on the coordinates in Fig. 1Then, can be obtained by solving Eq. (2). If is known, the two-way traveltime from the imaging point to the ’th antenna position can be given by 2.2.Self-Correlation Back-Projection Algorithm Developing StepsFor a given point with coordinates in the imaging area shown in Fig. 1, its projection point on the survey line is defined as with coordinates . The ’th antenna position with coordinates is the nearest to on the survey line . At the ’th antenna position, the A-scan of GPR data can be given by where is the number of sampling points; is the sampling instant of . The correlation coefficient between the vector and the vector is defined as where is the length of and , and and are the mean values of and , respectively. In the following, the SBP algorithm will be described for GPR imaging.
If the received GPR data have been amplified, the depth energy compensation algorithm will not be necessary. Otherwise, step 5 should be performed. The aforementioned steps will be repeated until all the points in the imaging scene are considered. 3.Experimental Results and AnalysisTo verify the performance of the proposed SBP algorithm, we design a synthetic and a practical field experiments for the objects imaging of complex shape. In addition, a serial contrast experiments are performed with the classic BP, the fast back projection (FBP)8, the SBP, and the 2-D depth migration algorithm with straight ray assumption.12 To solve Eq. (2), the binary search algorithm is used by both the classic BP and SBP algorithms; the approximate algorithm proposed in Ref. 8 is used by the FBP algorithm. We evaluate the imaging error by the absolute cumulative error (ACE), which is defined as where and are the normalized matrices of the objects position matrix and imaging results matrix, and and are the elements in the matrices and , respectively. In this paper, all the programs are run by MATLAB 7.1 on the same hardware.3.1.Tests with Synthetic DataAccording to Fig. 1, the geoelectric model of synthetic GPR data is formulated. The antennas are located 1 m above the surface of the surveyed structure. The scene is divided into two regions by . The upper region is air with relative permittivity 1 and conductivity . The lower region is homogeneous medium with relative permittivity 16 and conductivity . The buried objects are shown in Fig. 2 by the blue curve with relative permittivity 9. The GPR wave fields of the geoelectric model are simulated by using a finite-difference time-domain13 solution of the Maxwell’s equations. In the synthetic experiment, a Ricker wavelet is taken as a pulse source with a center frequency of 100 MHz. The B-scans of the synthetic GPR data are shown in Fig. 3. By using the mean removal algorithm, we can preprocess the synthetic GPR data to suppress the interference partially from the terrain echo and the coupled wave between receiving and transmitting antennas. The normalized imaging results of the classic BP, the FBP, and the SBP are shown in Figs. 4, 5, and 6, respectively. The aforementioned imaging results are shown in Table 1. The velocity model of the 2-D depth migration algorithm is shown in Fig. 7, where the electromagnetic wave travels at in air and in earth, respectively. The normalized imaging results of the 2-D depth migration algorithm are shown in Fig. 8. The theoretical distributions of buried objects are labeled with the white circle curve, which is convenient to analyze the imaging results. Table 1Performance comparisons table about different imaging methods.
The experimental results show that the objects shape cannot be distinguished from Fig. 3. It may be recovered by the GPR signal-processing techniques. Due to the fact that buried objects have complex shapes in this experiment, each A-scan signal is the stack of multiple reflection echo signals. For a given imaging point, the echo information extracted by the BP algorithm from each A-scan signal may be the stack of the echo information of many imaging points. It reduces the imaging accuracy of the classic BP and FBP algorithms. For the same B-scan in GPR imaging, at the same condition of imaging resolution, the larger the number of the sensing locations, the longer the running time of both the BP and FBP algorithms. The SBP algorithm can adaptively choose the valid echo information sequence, with its running time related to the specific B-scan. The larger the number of the valid echo information sequence, the longer the running time of the SBP algorithm. However, the SBP algorithm usually has faster calculation speed because the valid echo information sequence is generally located in the multiple nearest neighbors of the imaging point. The 2-D depth migration algorithm can distinguish the objects shape effectively. However, the locations of the objects are moved upwards compared with their theoretical distribution. The GPR antennas and objects are distributed in different media. The precision of velocity model is reduced by the air between the GPR antennas and the earth. It is the main reason for the aforementioned imaging error. The experimental results illustrate that the proposed SBP algorithm is superior to the existing BP algorithms in terms of computing speed and imaging accuracy. Compared with the 2-D depth migration algorithm, the proposed SBP algorithm has a significant advantage in providing a rough outline of buried objects without prior knowledge of the velocity distribution. 3.2.Field Test CaseThe experimental data are collected at the Georgia Institute of Technology14 and is publicly available in Ref. 15 in MATLAB format files, and the field map in the detecting location is shown in Fig. 9. The data consist of different burial and no-object scenarios, and are taken with a multistatic stepped-frequency continuous-wave GPR. The GPR is scanned over a region at a constant height above the surface of the ground. The scan region is discretized into a grid of 91 points. At each scan position, GPR takes 401 equally spaced frequency measurements from 60 MHz to 8.06 GHz with 20-MHz increments. The pulse source is a differentiated Gaussian pulse with a center frequency of 2.5 GHz. The GPR B-scan is shown in Fig. 10, which is selected from the aforementioned experimental data. There are five buried objects at , , 0, 20, and 45 in cross-range () dimension and . The distance between the transmitting and receiving antennas is 12 cm. The height of antennas is 27.8 cm above the surface of the surveyed structure. Each A-scan is acquired every 2 cm. The electromagnetic wave travels at in air and in sand.16 The field GPR data are preprocessed by the removing mean value method to suppress the interference partially from the terrain echo and the coupled wave between receiving and transmitting antennas. The normalized imaging result of the SBP algorithm is shown in Fig. 11. The theoretical distributions of buried objects are labeled with the black curve, which is convenient to analyze the imaging results. The experimental result shows that the proposed SBP algorithm can effectively recover the multiple objects’ position information. The calculated ACE for Fig. 11 is 426.5, and the running time of the SBP algorithm is 3.5 s. It further validates the effectiveness of the proposed SBP algorithm. 4.ConclusionBased on the existing BP algorithms, the SBP algorithm is proposed in this paper, which can reconstruct the shape of buried objects in GPR imaging. By setting up the correlation threshold, the SBP algorithm can adaptively choose the valid echo information sequence of the imaging points, which improves the image resolution. Because the valid echo information sequence is generally located in the multiple nearest neighbors of the imaging point, the SBP algorithm usually has a faster calculation speed. In addition, the imaging result is postprocessed by the depth energy compensation algorithm to improve the imaging resolution. The experimental results show that the SBP algorithm is superior to the classic BP and FBP algorithms in terms of computing speed and imaging accuracy. It has a significant advantage in providing a rough outline of buried objects without prior knowledge of the velocity distribution. The application of the proposed SBP algorithm is valuable in the imaging of underground objects with fast speed and high quality. AcknowledgmentsThis work was funded by the National Natural Science Foundation of China (Nos. 61102115, 61362020, and 61371186), the Natural Science Foundation of Guangxi Province (Nos. 2012GXNSFBA053177, 2013GXNSFAA019327, and 2013GXNSFFA019004), and the Guangxi Experiment Center of Information Science, Guilin University of Electronic Technology (No. 20130112). ReferencesA. Tzanis,
“Detection and extraction of orientation-and-scale-dependent information from two-dimensional GPR data with tuneable directional wavelet filters,”
J. Appl. Geophys., 89 48
–67
(2013). http://dx.doi.org/10.1016/j.jappgeo.2012.11.007 JAGPEA 0926-9851 Google Scholar
H. Zhang et al.,
“Matched filtering algorithm based on phase-shifting pursuit for ground-penetrating radar signal enhancement,”
J. Appl. Remote Sens., 8
(1), 083593
(2014). http://dx.doi.org/10.1117/1.JRS.8.083593 Google Scholar
J. Yang et al.,
“Sparse MIMO array forward-looking GPR imaging based on compressed sensing in clutter environment,”
IEEE Trans. Geosci. Remote Sens., 52
(7), 4480
–4494
(2014). http://dx.doi.org/10.1109/TGRS.2013.2282308 IGRSD2 0196-2892 Google Scholar
S. H. H. Lum, T. Tang and J. Lyall,
“Improved GPR image quality by truncating the bandwidth of the Fourier transform in 2D diffraction tomographic inversion,”
in IEEE Int. RF and Microwave Conf.,
391
–394
(2008). Google Scholar
M. Razavian, M. H. Hosseini and R. Safian,
“Time-reversal imaging using one transmitting antenna based on independent component analysis,”
IEEE Geosci. Remote Sens. Lett., 11
(9), 1574
–1578
(2014). http://dx.doi.org/10.1109/LGRS.2014.2301724 Google Scholar
L. Zhou and Y. Su,
“GPR imaging with RM algorithm in layered mediums,”
IEEE Geosci. Remote Sens. Lett., 8
(5), 934
–938
(2011). http://dx.doi.org/10.1109/LGRS.2011.2138116 Google Scholar
B. Hua and G. A. McMechan,
“Parsimonious 2D prestack Kirchhoff depth migration,”
Geophysics, 68
(3), 1043
–1051
(2003). http://dx.doi.org/10.1190/1.1581075 GPYSA7 0016-8033 Google Scholar
L. Zhou, C. Huang and Y. Su,
“A fast back-projection algorithm based on cross correlation for GPR imaging,”
IEEE Geosci. Remote Sens. Lett., 9
(2), 228
–232
(2012). http://dx.doi.org/10.1109/LGRS.2011.2165523 Google Scholar
A. Tzanis,
“The curvelet transform in the analysis of 2-D GPR data: signal enhancement and extraction of orientation-and-scale-dependent information,”
J. Appl. Geophys., 115 145
–170
(2015). http://dx.doi.org/10.1016/j.jappgeo.2015.02.015 JAGPEA 0926-9851 Google Scholar
Y. Zhang and T. Xia,
“OFDM and compressive sensing based GPR imaging using SAR focusing algorithm,”
Proc. SPIE, 9437 943725
(2015). http://dx.doi.org/10.1117/12.2083857 PSISDG 0277-786X Google Scholar
Y. Zhang et al.,
“Data analysis technique to leverage ground penetrating radar ballast inspection performance,”
in 2014 IEEE Radar Conf.,
463
–468
(2014). Google Scholar
B. Biondi and T. Tisserant,
“3D angle-domain common-image gathers for migration velocity analysis,”
Geophys. Prospect., 52
(6), 575
–591
(2004). http://dx.doi.org/10.1111/j.1365-2478.2004.00444.x GPPRAR 0016-8025 Google Scholar
G. R. Werner, C. A. Bauer and J. R. Cary,
“A more accurate, stable, FDTD algorithm for electromagnetics in anisotropic dielectrics,”
J. Comput. Phys., 255
(12), 436
–455
(2013). http://dx.doi.org/10.1016/j.jcp.2013.08.009 JCTPAH 0021-9991 Google Scholar
A. C. Gurbuz,
“Determination of background distribution for ground-penetrating radar data,”
IEEE Geosci. Remote Sens. Lett., 9
(4), 544
–548
(2012). http://dx.doi.org/10.1109/LGRS.2011.2174137 Google Scholar
W. R. Scott,
“Multistatic ground-penetrating radar experiments,”
(2007) http://users.ece.gatech.edu/~wrscott/ Google Scholar
T. Counts et al.,
“Multistatic ground-penetrating radar experiments,”
IEEE Trans. Geosci. Remote Sens., 45
(8), 2544
–2553
(2007). http://dx.doi.org/10.1109/TGRS.2007.900677 IGRSD2 0196-2892 Google Scholar
BiographyHairu Zhang received his BEng degree from Wuhan University of Science and Technology, Wuhan, China, in 2009 and obtained his MEng degree from Guilin University of Electronic Technology, Guilin, China, in 2012. Currently, he is a PhD student at Xidian University, Xi’an, China. His fields of interest include signal processing for ground-penetrating radar, computational intelligence, and adaptive signal processing. Shan Ouyang received his BEng degree from Guilin University of Electronic Technology, Guilin, China, in 1986 and obtained his MEng and PhD degrees in electronic engineering from Xidian University, Xi’an, China, in 1992 and 2000, respectively. He received the National Excellent Doctoral Dissertation of China in 2002. Currently, he is a professor and PhD supervisor at Guilin University of Electronic Technology. His fields of interest include signal processing for communications and radar. Guofu Wang received his BEng degree from Wuhan University of Technology, Wuhan, China, in 2002 and obtained his MEng and PhD degrees from Xi’an Institute Optics and Precision Mechanics of CAS, Xi’an, China, in 2005 and 2007, respectively. Currently, he is a professor at Guilin University of Electronic Technology. His fields of interest include signal processing for ground-penetrating radar and adaptive signal processing. Jingjing Li received his BEng degree from Henan Polytechnic University, Jiaozuo, China, in 2010 and obtained his MEng degree from Guilin University of Electronic Technology, Guilin, China, in 2014. Currently, she is a PhD student at Guilin University of Electronic Technology, Guilin, China. Her fields of interest include signal processing for radar. |