# An Approach for Predicting the Shape and Size of a Buried Basic Object on Surface Ground Penetrating Radar System

Authors: Nana Rachmana Syambas
Journal: International Journal of Antennas and Propagation
 (2012)

---

## Abstract

Surface ground-penetrating radar (GPR) is one of the radar technology that is widely used in many applications. It is nondestructive remote sensing method to detect underground buried objects. However, the output target is only hyperbolic representation. This research develops a system to identify a buried object on surface GPR based on decision tree method. GPR data of many basic objects (with circular, triangular, and rectangular cross-section) are classified and extracted to generate data training model as a unique template for each type of basic object. The pattern of object under test will be known by comparing its data with the training data using a decision tree method. A simple powerful algorithm to extract feature parameters of object which is based on linear extrapolation is proposed. The result showed that tested buried basic objects can be correctly predicted and the developed system works properly.


---

## Body

## 1. Introduction

Surface ground-penetrating radar (GPR) which has evolved as popular technology for nondestructive testing method has been widely used for the detection of underground buried objects since 1980. The interesting applications are for the detection of underground infrastructure [1] and the landmine [2]. GPR transmits short pulses of high frequency electromagnetic into the ground. These waves propagate with a velocity that depends on the dielectric property of medium. If the waves encounter a buried object with different reactive indices, some of the waves are reflected back and the receiver processes them to create a hyperbolic image as the object representation. In order to create an image, the GPR transmits electromagnetic pulses at a certain frequency for a certain time slot and samples the responses as an A-scan corresponding to a single position of the GPR. By moving the GPR in x direction, a collection of these A-scans called a B-scan which represent different GPR positions are constructed. Then it is processed to become a hyperbolic image. Ideally, the imaging process of GPR data is described in Figure 1.Figure 1
Imaging process of GPR data.Automatic generation of hyperbolic image from B-scan data has been widely reported using various methods and techniques such as Hough transform [3, 4], statistical and probability [5, 6], genetic algorithm [7], neural networks [4, 8, 9], and the latest technique is compressive sensing [10].Figure2 shows the B-scan of GPR data for basic circular and rectangular objects. From the image, it is hard to recognize the shape of the object. Image of an object generated by GPR is not corresponding to its geometrical representation.B-scan of GPR data for basic circular and rectangular objects.Regeneration of a geometric model is a much more complex procedure and is not usually attempted [11]. However, because the output target is only a hyperbolic representation, interpretation of a geometric buried object is depends on operator’s expertise and experience.Recent developments have shown that significant efforts have been focused towards the extraction of meaningful physical interpretation from GPR data. Researches on generating geometrical representation have been reported in [12–14]. Unfortunately, these reports are based on diffraction tomography or cross-hole GPR method, which is not practical to drill holes and is time consuming [15].This research tries to interpret geometrical representation from a B-scan of GPR data based on template matching. The pattern of object under test will be interpreted by comparing its data with the training data using a simple decision tree method.


## 2. The Proposed System

Proposed system for identification [16] is shown in Figure 3. Firstly an output GPR (a file*a.out) is extracted to generate reflected signal E(z) (V/m) as function of time (ns).   The difference of the output GPR format and the input format of the C.45 [17] software tool, needs a program to convert those formats.  The program is based on Java called DataCreator, NamesCreator, TstCreator, and AhCreatorwhich generates a file*.data, a file*.names, a file*.tst, and file*.ah, respectively. The first two files are needed to generate the training data and the last twofiles are neededto generate and to test the test data.Figure 3
Proposed system and algorithm number 1.The C.45 software tool needsabuilt-tree of training data to predict a decision. Generation of built-tree can be done after file*.data and file*.names have been generated. Values of training data are calculated to determine entropy and gain ratio for each attribute. Entropy is used to measure how informative is a node in a tree, and gain ranks attribute to build a decision tree of the training data.Interpretation of a test data can be done after the file*.tst is generated. Using the built-tree, a tested data record that has unknown attribute values can be classified by estimating the probability of the various possible results.Finally the ResultShower.jar program converts the output decision of C.45 in text form to visual pattern form to show the visual of the tested unknown object. These steps of the process are summarized in algorithm 1.As example on this experiment, training data consist of 30 samples of three types (circular, triangular, and rectangular) of object data. For each type object data has 10 different size samples (r: 7, 10, 12, 13, 15, 17, 18, 19, 20, and 23). Each sample has 10180 attributes of the file*.data, so the training data has totally 305400 attributes. Figure 4 shows the built-tree of the training data.Figure 4
Built-tree of the training data with algorithm number 1.


### 2.1. GPR Test Range Model

The model GPR test range is a box with 3 meter length, 3 meter wide, and 2 meter high (3×3×2 meter), containing a half full of sand. A 1 GHz transmitter and receiver of GPR, which has 10 centimetre fix distance and located 20 centimetre above the sand, moves to x direction in 5 centimetre intervals until 1 meter length (20 scan positions). In each interval, the pulse signal with 12 ns width is transmitted. The sand’s parameter relative permittivity is 3.0, conductivity is 10−2 Siemens/meter, and relative permeability is 1. A simple buried object (with three types of patterns: circular, triangular, or rectangular cross-section) is made from perfect conductor material and is lied in the sand at 50 cm depth. Figure 5 shows the model of simulation.Figure 5
GPR test range model.For training purpose, 30 samples of three types (circular, triangular, and rectangular) of object data are generated. Each type of object data has 10 different size samples. Test data consist of 12 samples of the three types of objects in which each type has 4 different sample sizes.Successful GPR modelling based on the finite-difference time-domain (FDTD) theory [18] has been reported by many authors [19–22].Giannopoulos [19]has developed GPR modelling facilities called GprMax2D/3D ver. 2.0 which are used to generate the GPR data. The output simulator can be a file*b.out and file*a.out to represent hyperbolic image in binary file and ASCII file, respectively. Sampling experiment by using real GPR tool will be done in GPR test range ICTR-ITB laboratory to validate and verify the GPR data simulator. An example of input and output GPR data simulator is shown in Figure 6. The research takes files*a.out as input of the proposed system to generate training data and tests data files.Figure 6
An input and an output GPR data simulator.


### 2.2. Testing Scenarios

There are scenarios to cover both the same environment test condition and the different environment test condition. The test data for each scenario consist of 12 samples of the three types of object in which each type has 4 different sample sizes (r: 8, 9, 21, and 22). The total samples test data for all scenarios are 96 samples. These scenarios are as follows.(1)
Tested and training objects are in the same 50 centimetre depth condition.
(i)
Scenario number 0. The sand and the tested object have the same parameters as in the training data (the sand:εr=3, μr=1, and the object: perfect conductor). Tested and training objects have the same environmental condition.(ii)
Scenario number 1. The sand has the same parameter as in training data and the tested object has different material (the sand:εr=3, μr=1, and the object: fresh water, ρr=81).(iii)
Scenario number 2. The tested object has the same parameter as training data and the sand has different permittivity (sand:εr=4, μr=1, and the object: perfect conductor).(iv)
Scenario number 3. The sand medium and the tested object have different parameters compared to training data (sand:εr=4, μr=1, and the object: fresh water, ρr=81).(2)
Tested and training objects are in the different depth.
(i)
Scenario number 4. The sand and the tested object have the same parameters as in the training data (the sand:εr=3, μr=1, and the object: perfect conductor), and depth of the tested object is 40 cm.(ii)
Scenario number 5. The sand and the tested object have the same parameters as in the training data (the sand:εr=3, μr=1, and the object: perfect conductor), and depth of the tested object is 55 cm.(iii)
Scenario number 6. The sand has the same parameter as in training data and the object to be tested has different material (the sand:εr=3, μr=1, and the object: fresh water, ρr=81), and the depth of the tested object is 40 cm.(iv)
Scenario number 7. The sand has the same parameter as in training data and the tested object has different material (the sand:εr=3, μr=1, and the object: fresh water, ρr=8), and depth of the tested object is 55 cm.


### 2.3. The Test Result Based on Algorithm Number 1

An example of interpretation of the test result for scenario number 6 is shown in Figure7. This scenario in which tested and training objects are in the different environmental condition, generates the extreme result, that is, the most of interpretation is false. On the other hand, scenario number 0, in which tested andtraining objects are in the same environmental condition, generates the best result, that, is all interpretation is true.Figure 7
An example of the partial test result for scenario number 6.The complete test results are summarized in Table1. From 96 tested data samples, the system can truly well interpret only for 50 samples (52%). Algorithm 1, which represent direct comparison between test and training data with template matching method, has inaccurate interpretation for condition where the tested and the training object have different environmental condition.Table 1
The test result based on algorithm number 1.
DUT
Sc no.
Sc no. 0
Sc no. 1
Sc no. 2
Sc no. 3
Sc no. 4
Sc no. 5
Sc no. 6
Sc no. 7
TRI_8
True
True
True
True
True
True
False
True
TRI_9
True
True
True
True
True
True
False
True
TRI_21
True
True
True
True
True
True
False
True
TRI_22
True
True
True
True
True
True
False
True
CIR_8
True
True
False
False
False
False
False
False
CIR_9
True
True
False
False
False
False
False
False
CIR_21
True
True
False
False
False
False
False
False
CIR_22
True
True
False
False
False
False
False
False
REC_8
True
True
False
False
True
False
False
False
REC_9
True
True
False
False
True
False
False
False
REC_21
True
True
False
False
True
False
True
False
REC_22
True
True
False
False
True
False
True
False


## 2.1. GPR Test Range Model

The model GPR test range is a box with 3 meter length, 3 meter wide, and 2 meter high (3×3×2 meter), containing a half full of sand. A 1 GHz transmitter and receiver of GPR, which has 10 centimetre fix distance and located 20 centimetre above the sand, moves to x direction in 5 centimetre intervals until 1 meter length (20 scan positions). In each interval, the pulse signal with 12 ns width is transmitted. The sand’s parameter relative permittivity is 3.0, conductivity is 10−2 Siemens/meter, and relative permeability is 1. A simple buried object (with three types of patterns: circular, triangular, or rectangular cross-section) is made from perfect conductor material and is lied in the sand at 50 cm depth. Figure 5 shows the model of simulation.Figure 5
GPR test range model.For training purpose, 30 samples of three types (circular, triangular, and rectangular) of object data are generated. Each type of object data has 10 different size samples. Test data consist of 12 samples of the three types of objects in which each type has 4 different sample sizes.Successful GPR modelling based on the finite-difference time-domain (FDTD) theory [18] has been reported by many authors [19–22].Giannopoulos [19]has developed GPR modelling facilities called GprMax2D/3D ver. 2.0 which are used to generate the GPR data. The output simulator can be a file*b.out and file*a.out to represent hyperbolic image in binary file and ASCII file, respectively. Sampling experiment by using real GPR tool will be done in GPR test range ICTR-ITB laboratory to validate and verify the GPR data simulator. An example of input and output GPR data simulator is shown in Figure 6. The research takes files*a.out as input of the proposed system to generate training data and tests data files.Figure 6
An input and an output GPR data simulator.


## 2.2. Testing Scenarios

There are scenarios to cover both the same environment test condition and the different environment test condition. The test data for each scenario consist of 12 samples of the three types of object in which each type has 4 different sample sizes (r: 8, 9, 21, and 22). The total samples test data for all scenarios are 96 samples. These scenarios are as follows.(1)
Tested and training objects are in the same 50 centimetre depth condition.
(i)
Scenario number 0. The sand and the tested object have the same parameters as in the training data (the sand:εr=3, μr=1, and the object: perfect conductor). Tested and training objects have the same environmental condition.(ii)
Scenario number 1. The sand has the same parameter as in training data and the tested object has different material (the sand:εr=3, μr=1, and the object: fresh water, ρr=81).(iii)
Scenario number 2. The tested object has the same parameter as training data and the sand has different permittivity (sand:εr=4, μr=1, and the object: perfect conductor).(iv)
Scenario number 3. The sand medium and the tested object have different parameters compared to training data (sand:εr=4, μr=1, and the object: fresh water, ρr=81).(2)
Tested and training objects are in the different depth.
(i)
Scenario number 4. The sand and the tested object have the same parameters as in the training data (the sand:εr=3, μr=1, and the object: perfect conductor), and depth of the tested object is 40 cm.(ii)
Scenario number 5. The sand and the tested object have the same parameters as in the training data (the sand:εr=3, μr=1, and the object: perfect conductor), and depth of the tested object is 55 cm.(iii)
Scenario number 6. The sand has the same parameter as in training data and the object to be tested has different material (the sand:εr=3, μr=1, and the object: fresh water, ρr=81), and the depth of the tested object is 40 cm.(iv)
Scenario number 7. The sand has the same parameter as in training data and the tested object has different material (the sand:εr=3, μr=1, and the object: fresh water, ρr=8), and depth of the tested object is 55 cm.


## 2.3. The Test Result Based on Algorithm Number 1

An example of interpretation of the test result for scenario number 6 is shown in Figure7. This scenario in which tested and training objects are in the different environmental condition, generates the extreme result, that is, the most of interpretation is false. On the other hand, scenario number 0, in which tested andtraining objects are in the same environmental condition, generates the best result, that, is all interpretation is true.Figure 7
An example of the partial test result for scenario number 6.The complete test results are summarized in Table1. From 96 tested data samples, the system can truly well interpret only for 50 samples (52%). Algorithm 1, which represent direct comparison between test and training data with template matching method, has inaccurate interpretation for condition where the tested and the training object have different environmental condition.Table 1
The test result based on algorithm number 1.
DUT
Sc no.
Sc no. 0
Sc no. 1
Sc no. 2
Sc no. 3
Sc no. 4
Sc no. 5
Sc no. 6
Sc no. 7
TRI_8
True
True
True
True
True
True
False
True
TRI_9
True
True
True
True
True
True
False
True
TRI_21
True
True
True
True
True
True
False
True
TRI_22
True
True
True
True
True
True
False
True
CIR_8
True
True
False
False
False
False
False
False
CIR_9
True
True
False
False
False
False
False
False
CIR_21
True
True
False
False
False
False
False
False
CIR_22
True
True
False
False
False
False
False
False
REC_8
True
True
False
False
True
False
False
False
REC_9
True
True
False
False
True
False
False
False
REC_21
True
True
False
False
True
False
True
False
REC_22
True
True
False
False
True
False
True
False


## 3. The Feature Extraction

An effort to enhance the performance of the system is done by adding a feature extraction of the object to algorithm 1. This improved the algorithm called as algorithm 2.Investigation of reflected signalE(z) received by receiver for 30 samples of three types of (circular, triangular, and rectangular) object of training data carries out an important information, that is, each type of object has a unique pattern signal. In the same type of object, different size will generate different amplitude signal with the same pattern signal. This phenomenon can be seen in Figure 8.Received signal pattern for three types of the basic object.
(a)
Rectangular
(b)
Circular
(c)
TriangularAfterward, the unique pattern signal of object is extracted as the objectfeature, that is, arepresentation of the object. To simplify the feature and to reduce the number attribute of complex extracted signal pattern, the time variable is excluded by sorting values of theE(z) from the smallest to the largest, and then a linear extrapolation is applied. Changing the testing condition (εr,μr,ρr, and depth) will shift the position in time axis and change the amplitude signal without changing the pattern of object.By this process, each reflected signal pattern of an object, which in algorithm 1 has 10180 attributes, will be transformed to a straight line that only has 3 attributes, called slope of the line (a), intersection point of the line with the E(z) axis (b), and standard deviation (s). Furthermore, these attributes are used as the object representations are processed using template matching method by comparing attribute of tested data to training data, the same as in algorithm 1.Figure9 shows algorithm 2 illustrated in each step, starting from extracting the reflected signal, extracting the feature, and matching the tested and training attribute data using decision tree to generate the interpretation result.Figure 9
Proposed algorithm number 2 and its illustration.


### 3.1. The Test Result Based on Algorithm Number 2

The system based on algorithm 2 is tested by the test data and scenarios the same as in the system based on algorithm number 1. An example of interpretation of the partial test result for scenario number 5 is shown in Figure10 and all the test result is summarized in Table 2. For 96 test samples in different scenarios, the result showed that the system correctly interprets all the samples. Algorithm number 2 has accurate interpretations for all environmental testing conditions.Table 2
The test result based on algorithm number 2.
DUT
Sc no.
Sc no. 0
Sc no. 1
Sc no. 2
Sc no. 3
Sc  no. 4
Sc no. 5
Sc no. 6
Sc no. 7
TRI_8
True
True
True
True
True
True
True
True
TRI_9
True
True
True
True
True
True
True
True
TRI_21
True
True
True
True
True
True
True
True
TRI_22
True
True
True
True
True
True
True
True
CIR_8
True
True
True
True
True
True
True
True
CIR_9
True
True
True
True
True
True
True
True
CIR_21
True
True
True
True
True
True
True
True
CIR_22
True
True
True
True
True
True
True
True
REC_8
True
True
True
True
True
True
True
True
REC_9
True
True
True
True
True
True
True
True
REC_21
True
True
True
True
True
True
True
True
REC_22
True
True
True
True
True
True
True
TrueFigure 10
An example of the partial test result for scenario number 5.


## 3.1. The Test Result Based on Algorithm Number 2

The system based on algorithm 2 is tested by the test data and scenarios the same as in the system based on algorithm number 1. An example of interpretation of the partial test result for scenario number 5 is shown in Figure10 and all the test result is summarized in Table 2. For 96 test samples in different scenarios, the result showed that the system correctly interprets all the samples. Algorithm number 2 has accurate interpretations for all environmental testing conditions.Table 2
The test result based on algorithm number 2.
DUT
Sc no.
Sc no. 0
Sc no. 1
Sc no. 2
Sc no. 3
Sc  no. 4
Sc no. 5
Sc no. 6
Sc no. 7
TRI_8
True
True
True
True
True
True
True
True
TRI_9
True
True
True
True
True
True
True
True
TRI_21
True
True
True
True
True
True
True
True
TRI_22
True
True
True
True
True
True
True
True
CIR_8
True
True
True
True
True
True
True
True
CIR_9
True
True
True
True
True
True
True
True
CIR_21
True
True
True
True
True
True
True
True
CIR_22
True
True
True
True
True
True
True
True
REC_8
True
True
True
True
True
True
True
True
REC_9
True
True
True
True
True
True
True
True
REC_21
True
True
True
True
True
True
True
True
REC_22
True
True
True
True
True
True
True
TrueFigure 10
An example of the partial test result for scenario number 5.


## 4. Calculation of the Object Size

Based on reflected signalE(z) for 30 samples training data of three types (circular, triangular, and rectangular), the maximum and minimum signals amplitude can be determined. In the same type of object, different size will generate different amplitude signal with the same pattern signal. By linear or logarithmic extrapolation of signals amplitude, a formula to calculate another size of objects in the same object type with similar environmental condition (d, εr, μr, and σr) can be generated as in Figure 11. Unfortunately for the triangular object, the generated formula is constant. Constant relationship between amplitude and size of triangular shape can be due to the backscattered rays which are not directed towards the antenna.Figure 11
Relationship of amplitude and size for same type of the object.The formula for a tested object has different condition (d, εr, μr and σr) with the training object, which is need to be corrected by a correction factor as can be derived from Figure 12.Figure 12
Reflected signal of nonperfect conductor object.Intrinsic impedance, transmission, and reflection factors are as follows:(1)η=μεη=με,T=2η1η1+η0,Γ=η2-η1η1+η2,Tup=21/εr11/εr0+1/εr1=21+εr1/εr0=21+εr1,Tpu=21/εr01/εr0+1/εr1=2εr1/εr01+εr1/εr0=2εr11+εr1,Tpb=21/εrb1/εrb+1/εr1=2εr1/εrb1+εr1/εrb=21+εrb/εrp,Tbp=21/εr11/εr1+1/εrb=2εrb/εr11+εrb/εr1=21+εrp/εrb,Γpb=ηb-ηpηb+ηp=1/εrb-1/εrp1/εrb+1/εrp=εrb-εrpεrb+εrp,Γbp=ηp-ηbηp+ηb=1/εrp-1/εrb1/εrp+1/εrb=εrp-εrbεrp+εrb,Ez=Ez0×e-2αpd×Tup×Γpb×Tpu+Ez0e-2αpde-2αbdb×Tup×Tpb×Γbp×Tbp×Tpu,αp=ωμrpεrp2[1+(σrpωεrp)2-1]1/2Npm,αb=ωμrbεrb2[1+(σrbωεrb)2-1]1/2Npm.
The correction factor for different environmental condition is(2)k=|Ez1||Ez2|=(Ez0×e-2αpd×Tup×Γpb×Tpu+Ez0e-2αpde-2αbdb×Tup×Tpb×Γbp×Tbp×Tpu)1(Ez0×e-2αpd×Tup×Γpb×Tpu+Ez0e-2αpde-2αbdb×Tup×Tpb×Γbp×Tbp×Tpu)2.As example for scenario number 1, correction factor k became
(3)|Ez1||Ez2|=1(Γpb+e-2αbdb×Tpb×Γbp×Tbp).
For εrp=3, d=50 cm and εrb=81, the factor correction is
(4)Tpb=21/εrb1/εrb+1/εr1=2εr1/εrb1+εr1/εrb=21+εrb/εrp=21+81/3=0.323,Tbp=21/εr11/εr1+1/εrb=2εrb/εr11+εrb/εr1=21+εrp/εrb=21+3/81=1.5,Γpb=εrb-εrpεrb+εrp=81-381+3=0.677,Γpb=εrp-εrbεrp+εrb=3-8181+3=-0.677,αb=ωμε2[1+(σωε)2-1]1/2=2π×109812×3×108[1+(0.0022)2-1]1/2=0.2,|Ez1||Ez2|=k=1(0.677+e-2(0.2)0.5×0.323×0.677×1.5)=1.06.Table3 shows that the correction factor for all scenarios has been investigated in this research.Table 3
Correction factor for all scenarios.
Scenario no. 0
k
=
1
Scenario no. 1
k
=
1.06
Scenario no. 2
k
=
0.9
Scenario no. 3
k
=
0.95
Scenario no. 4
k
=
0.8
Scenario no. 5
k
=
1.13
Scenario no. 6
k
=
0.85
Scenario no. 7
k
=
1.2Using the formula and the correction factor, size of tested object can be predicted. Table4 shows the prediction result for rectangular and circular tested object. Unfortunately, due to the constant relationship between signal amplitude and object size in generated formula for triangular object, this approach did not success in predicting the triangular tested object. The overestimate or underestimate size prediction for the circular and rectangular objects occur due to the linear or logarithmic approximation of the experimental result graph as shown in Figure 11.Table 4
Size of prediction of tested object.
Size of tested objectr (cm)
Prediction result (cm)
Rectangular
Circular object
Scenario
Sizer
Sizer
8
7.5
8.0
Sc no. 0
9
8.4
8.9
21
21.0
21.8
22
21.2
21.8
8
8.1
10.6
Sc no. 4
9
9.2
11.7
21
21.0
26.8
22
21.7
27.1
8
7.5
9.6
Sc no. 5
9
8.3
10.5
21
20.6
22.7
22
21.0
22.7
8
7.9
9.6
Sc no. 2
9
8.9
11.2
21
22.2
23.5
22
22.6
23.9
8
7.5
5.1
Sc no. 1
9
8.5
5.6
21
21.2
19.6
22
21.4
22.8
8
7.6
6.7
Sc no. 7
9
8.4
7.7
21
20.1
20.6
22
20.7
23.1
8
8.7
6.6
Sc no. 6
9
10.0
8.8
21
23.1
21.9
22
24.1
25.2
8
8.4
5.5
Sc no. 3
9
9.5
7.1
21
21.1
20.6
22
21.6
23.0


## 5. Testing Measurement in the Laboratory

Testing data based on sampling experiment using real GPR tool has been done in GPR test range at ICTR-ITB laboratory. The tool consists of test range, GPR scanner, spectrum analyser, and personal computer with the software tools for controlling scanner and processing data. Figure13 shows the GPR test range in ICTR-ITB laboratory.GPR test range in ICTR-ITB laboratory.Three types objects are tested. These objects are buried objects (circular, triangular, and rectangular cross-section with diameter 9, 20, and 8 cm, resp.). The object made from perfect conductor material is buried in the sand at 35, 25, and 25 cm depth, respectively. Figure14 shows the measurement result for the circular object.A-scan and B-scan of the tested circular object.
(a)
Files A-scan of 20 trace points
(b)
A-scan for the first trace
(c)
The received signal (Ez)
(d)
The raw B-scan
(e)
The processed B-scanFrom Figure15 which represents the zooming and adding cursor on B-scan, the depth of buried object as the top of hyperbolic line in y axis is shown by the cursor on 517, surface of the sand on 491, and position of the floor on 566. The real height of the sand is 100 cm represented by value 566-491=75 point and the position of the object is represented by value 566-517=26 point. So, the real position of the buried object is 26/75×100cm=34.7 cm depth.Figure 15
Zooming and marking the B-scan.Algorithm number 2 has been applied to the measured tested data. The reflected signal received by GPR receiver that is represented by 20 A-scan of 20 point trace experiment result is extracted to get only the reflected signal of buried object by capturing the reflected signal in the axis range where the hyperbolic line is located and then converting the signal to the binary data in ASCII file for further manipulation.By sorting values of the binary data from the smallest to the largest and then applying a linear extrapolation, the feature parameter of the buried object can be obtained, that is,a=0.00007, b=0.0307, and c=0.8829.Matching the feature parameter with the built-tree that has been generated from training data using the C.45 tool as illustrated in Figure9 will generate a decision of buried object interpretation. The extract feature parameter and the interpretation result of this testing are shown in Figure 16.Figure 16
Feature extraction and test result of the circular object.The same process is applied to the measured tested data for triangular and rectangular objects, generating feature parameter with values ofa=0.0001, b=0.048, and c=0.9524, and values of  a=0.0001, b=0.055, and c=0.8592, respectively.The amplitude of extracted signalEz for circular object is 119 mV/m. Entering this value to the equation y=2.33x+43.24 results in the size value, x of 32.5 cm. The testing experiment condition in laboratory is different from training condition, so the size value should be corrected using correction factor derives from (2) as follows:
(5)k=Ez01×e-2αp1d1×Tup1×Tpu1Ez02×e-2αp2d2×Tup2×Tpu2,d1=50,d2=35,ε0=8.86×10-12,εr1=3,σ=0.01,εr2=5,σ=0.01,k=1×e-2(1.09)(0.5)×0.732×1.2685×e-2(0.84)(0.35)×0.618×1.382=0.29.
The size prediction of the test circular object can be derived by multiplying x by factor k, that is, 0.29×32.5=9.4 cm.Similar process can be applied to rectangular and triangular tested object. All the test result of measured test data in the ICTR-ITB laboratory can be summarized in Figure17.Figure 17
The test result of measured data in the ICTR_ITB Laboratory.


## 6. Conclusions

This paper proposed a system for identifying and predicting target shape and size of buried basic object on surface ground-penetrating radar system. Feature parameters of object are extracted by using simple proposed algorithm which is based on linear extrapolation. 96 samples of basic objects which are laid underground and have both same and different condition with the data training are tested. Results show that all basic objects can be correctly interpreted and the proposed algorithm number 2 can improve the performance of the shape object interpretation from 52% to 100%. Three sampling experiments of buried object have been done in GPR test range ICTR-ITB laboratory and the result of the test showed that the developed system works properly.The proposed method is appropriate to be implemented in the industry to enhance the capability of GPR tools. The next output GPR target is not only in a hyperbolic representation but also can directly interpret the shape and the size of the detected object.  The main application is for the underground utility infrastructure detection and interpretation such as the telecommunication and electrical cable and duct which have regular shape and size of the rectangular or circular cross-section object.However, there is a weakness in predicting the size of triangular object due to the constant relationship between amplitude and size of the triangular shape that is generated from B-scan GPR data. For further research, instead of B-scan GPR data, using C-scan GPR data as the system input might solve this problem, even though the generation of C-scan GPR data needs more time.



---

*Source: 919741-2012-08-09.xml*