# Buried Object Detection by an Inexact Newton Method Applied to Nonlinear Inverse Scattering

Authors: Matteo Pastorino; Andrea Randazzo
Journal: International Journal of Microwave Science and Technology
 (2012)

---

## Abstract

An approach to reconstruct buried objects is proposed. It is based on the integral equations of the electromagnetic inverse scattering problem, written in terms of the Green’s function for half-space geometries. The full nonlinearity of the problem is exploited in order to inspect strong scatterers. After discretization of the continuous model, the resulting equations are solved in a regularization sense by means of a two-step inexact Newton algorithm. The capabilities and limitations of the method are evaluated by means of some numerical simulations.


---

## Body

## 1. Introduction

Subsurface imaging is an important subject in several applicative areas, including seismic and geophysical prospecting, nondestructive, testing and medical diagnostics [1–39]. Electromagnetic techniques are widely applied to face this problem. Ground penetrating radar [40] is the basic instrumentation able to retrieve discontinuities in the lower subspace (the ground). Processing of GPR data can allow an improvement in the detection and localization capabilities of the system. In some cases, neural network-based approaches have been adopted [41]. Moreover, inverse-scattering-based techniques have been proposed, too. By using these methods, the measured values of the field scattered by the discontinuities present in the lower half space are “inverted” in order to obtain the spatial distributions of dielectric parameters (e.g., the dielectric permittivity and the electric conductivity) inside a fixed inspection domain. Approaches of this kind have been proposed, for example, for the detection of water leaking from pipes [42]. If an accurate model of the buried scatterers is available, it would be possible to retrieve information about the buried objects in terms of their dielectric parameters.The present authors developed in [43] an inverse-scattering-based method for buried object imaging under the second-order Born approximation (SOBA) [21, 44–48]. By using the SOBA, the reconstruction of buried objects has been obtained with a better accuracy with respect to linearized approximations. In particular, due to the nonlinear nature of the scattering equations, in which the electric field inside the inspected region is an unknown quantity as well as the spatial distributions of the dielectric properties of the medium, the SOBA allows to approximately retrieve the field distribution, too.However, in order to detect strong discontinuities, the exact nonlinear inverse scattering formulation must be considered. This approach has been faced in [49, 50], in which strong scatterers have been reconstructed in free space under tomographic imaging conditions. The approach proposed in [49], which is based on the use an inexact Newton method, has been experimentally tested in [51]. The same method has been successfully applied in [52] within the contrast source formulation (under free-space conditions).In the present paper, the approach developed in [43] for buried object detection under the SOBA is extended to treat strong scatterers. In this case, the exact nonlinear inverse scattering formulation, in terms of the half-space Green’s function, is taken into account. In particular, starting from the integral equations for the scattered field for cylindrical objects, slides of the spatial distributions of the dielectric parameters are reconstructed by numerically solving these equations in a regularization sense.The paper is organized as follows. In Section2 the mathematical formulation of the approach is reported. In Section 3, the application of the inexact Newton method is described. Numerical results are reported in Section 4, and conclusions are drawn in Section 5.


## 2. Mathematical Formulation

Let us consider a cylindrical object buried in a homogeneous-half-space medium (Figure1). A 2D transverse-magnetic (TM) configuration is assumed. The dielectric properties of the buried object are assumed to be independent of the z coordinate. The transmitting and receiving antennas are located in the transverse plane and generate a z-polarized electric field in both the half spaces (this field is independent of the axial coordinate, too).Figure 1
Geometrical configuration of the problem.The dielectric characteristics of the investigation domainS,in the lower half space, are described by the so-called object function, which is defined as
(1)τ(r)=jω[ϵ(r)-ϵb],
where ϵ(r) is the complex dielectric permittivity at point r, and ϵb denotes the dielectric permittivity of the soil.The investigation area is illuminated byV line-current sources, with angular frequency ω (in the following, the ejωt time dependence is assumed and omitted). A mixed illumination configuration (partially borehole and partially from the top), which allows direct comparisons with the results reported in [43], is exploited. Under such assumptions, also the scattered electric field is z-polarized and independent of z. Consequently, only the transverse plane can be considered for the imaging procedure.The scattered electric field,uscatt(v), v=1,…,V, is collected for any incident wave in M measurement points sm(v), which constitute the so-called observation domain O. The dielectric properties of S can be related to the scattered field by means of the following set of nonlinear integral equations:
(2)uscatt(v)(r)=Glowerupperτutot(v)(r),r∈O,utot(v)(r)=uinc(v)(r)+Glowerlowerτutot(v)(r),r∈S,
where Glowerupper(·)(r)=jωμ0∫S(·)gHS(r,r′)dr′, r∈O and Glowerlower(·)(r)=jωμ0∫S(·)gHS(r,r′)dr′, r∈S are the data and state operators, where gHS is the half-space Green’s function [43] (superscripts in the operators recall where they are calculated, whereas subscripts recall where they are integrated).Equations (2) can be rewritten in a compact operator form as
(3)HS(x)=y,
where x=[τutot(1)⋯utot(V)]t is an array containing the values of the object function and of the total electric field inside the investigation area for any illumination. In (3), y=[uscatt(1)⋯uscatt(V)uinc(1)⋯uinc(V)]t is an array containing the values of the scattered electric field acquired by the measurement probes in the upper half space and the values of the incident electric field inside the investigation area. The half-space operator HS can be expressed, according to (2), as
(4)HS(x)=[Glowerupperτutot(1)⋮Glowerupperτutot(V)utot(1)-Glowerlowerτutot(1)⋮utot(V)-Glowerlowerτutot(V)].
As it is well known, the inverse problem (4) is non-linear and ill-posed. Consequently, nonlinear regularizing methods must be used.


## 3. Regularized Solution of (3)

Due to the mentioned ill posedness, (3) must be solved in a regularized sense. The inexact Newton regularizing scheme developed in [49, 50] is applied here. This method includes two nested loops, in which a linearization (with respect to the current solution guess xk) of the operator HS is constructed (outer loop) and solved by means of the truncated Landweber method (inner loop). Exactly as in [49], the algorithm performs the following steps.(1)
Fork=0 (k denotes the outer iteration index), a starting guess x0 is chosen (initialization phase).(2)
The following linearized version of equation is constructed(5)HSxk′h=yk,where h=[hτhu(1)⋯hu(V)]t, yk=y-HS(xk), and HSxk′ is the Fréchet derivative of HS at point xk. A regularized solution h^k to (5) is obtained by applying an iterative truncated Landweber algorithm, for which we have
(6)hk,0=0,hk,l+1=hk,l-γHSxk′*(HSxk′hk,l-yk),where HSxk′* is the adjoint operator of HSxk′.(3)
The current solution is updated by using the following relation:(7)xk+1=xk+h^k.(4)
A predefined stopping rule is finally applied. If this criterion is not satisfied, a fixed maximum number of outer iterations,kmax, is performed. In both cases, at any iteration, the procedure can be stopped, otherwise steps (2)–(4) are repeated.It should be noted that step(2) is actually the inner loop of the algorithm, whereas steps (1)–(3) constitute the outer loop.


## 4. Numerical Results

The developed method has been tested by means of numerical simulations. The working frequency has been set equal to 300 MHz (i.e., similar to that used in ground penetrating radars). A homogeneous lossy soil characterized by a relative dielectric permittivityϵrb=4 and an electric conductivity σb=0.01S/m has been considered (in order to simulate a dry sandy soil [53]). Similar to [43], a mixed measurement configuration is assumed, where reflection and cross-borehole operating modes are mixed together in order to acquire as much information as possible. Correspondingly, the antennas are located at points (Cartesian coordinates)
(8)rnant={(-23λ0,-43λ0+(n-1)2λ021),n=1,…,15,(-23λ0+(n-15)2λ021,0),n=16,…,28,(23λ0,-(n-29)2λ021),n=29,…,43,
where λ0 is the free-space wavelength in the upper medium. A subset of these antennas act both as transmitter (TX) and receiver (RX). When an antenna operates in TX mode, all the remaining ones collect the scattered electric field (i.e., they work in RX mode). Consequently, for every illumination, M=42 field samples are collected. In particular, in the cases reported in the following, V=7 antennas are used. They are located at positions denoted by the indexes given by n=7(v-1)+1,v=1,…,V.The investigation area is a square domain of dimensionsλ0×λ0, whose upper side coincides with the air-ground interface. In order to numerically solve the equations involved in the inverse-scattering problem, such domain has been discretized into N=20×20 square subdomains.The scattered field values used as input data for the considered simulations have been numerically obtained by using a computational code based on the method of moments [54]. In order to avoid the so-called inverse crime, a finer discretization (such as the ratio of the cell dimensions in the direct and inverse problems is not an integer number) has been used for the computation of the synthetic data. Moreover, in order to simulate a more realistic scenario, a Gaussian noise with zero mean value and variance corresponding to a signal-to-noise ratio SNR=25dB has also been added to the computed data.In all the simulations reported in this section, the initialization phase has been always performed by assuming an empty investigation domain in the lower half space. Consequently, in this phase, the object function is equal to zero, and the total electric field is assumed equal to the one transmitted in the lower half space.The first reported case concerns the reconstruction of a single circular cylinder of radiusa=0.2λ0 centered at robj=(0.2λ0,-0.5λ0). The dielectric properties of the cylinder are ϵrobj=5.6 and σobj=0.01S/m. The parameters of the inexact Newton algorithm are the following: maximum number of inner iterations, kLWmax=30; maximum number of outer iterations, kINmax=100. Figure 2 reports the behavior of two error parameters for different values of the number of outer iterations of the iterative procedure. In particular, the two parameters are defined as
(9)rdata=‖uactual-ucalc‖2‖uactual‖2,rsolution=‖τactual-τrec‖2‖τactual‖2,
with τactual and τrec being the actual and reconstructed permittivity distributions (at the current iteration) and uactual and ucalc the measured and computed values of the electric field data (at the current iteration, too). As can be seen from Figure 2, the plot of the “residual” on the data (the error parameter rdata) decreases monotonically as expected. Moreover, the plot of the solution “residual” (the error parameter rsolution) shows a typical semiconvergence behavior, that is, after a certain value of the iteration number the reconstruction quality degrades [55, 56].Figure 2
Data and solution residuals versus the outer iteration number (error parameters of (9)). Single circular cylinder buried in a lossy half space.The final image of the distribution of the relative dielectric permittivity inside the investigation region (corresponding to the optimum iteration number) is shown in Figure3. As can be seen, the method is able to correctly retrieve the dielectric profile of the investigation area. In particular, both the shape of the buried inclusion and the value of its dielectric permittivity are correctly identified.Figure 3
Reconstructed distribution of the relative dielectric permittivity. Single circular cylinder buried in a lossy half space.For the present case, an analysis of the effects of the noise has also been performed. In particular, the mean relative error on the reconstruction, defined as theL1-norm
(10)e=1N∑n=1N|ϵnactual-ϵnrec||ϵnactual|
has been calculated. In (10), ϵnactual and ϵnrec denote the values of the original and reconstructed complex permittivity inside the nth discretization subdomain of the investigation region. Moreover, the optimal value of the outer iteration have been evaluated, too. The obtained results are reported in Table 1. As expected, as the noise level increases, the reconstruction “quality” decreases. However, even in presence of strong noise levels, the method is able to yield rather good reconstructions (e.g., the error parameter e is equal to about 8% error for SNR=15dB).Table 1
Mean relative errors and optimal iteration numbers versus the signal-to-noise ratio. Single circular cylinder buried in a lossy half space.
SNR
25
20
15
10
Mean relative error
0.065
0.078
0.085
0.12
Optimum outer iteration number
53
34
14
5In the second set of simulations, the reconstruction of a more complex configuration has been considered. In particular, the presence of a double-layer circular cylinder has been assumed. The outer layer has radiusa1=0.2λ0 and is characterized by a relative dielectric permittivity ϵr1=5.6 and an electric conductivity σ1=0.01 S/m, whereas the inner one has a radius a2=0.1λ0 and is characterized by a relative dielectric permittivity ϵr2=8 and an electric conductivity σ2=0.01 S/m. The object is centered at robj=(0.2λ0,-0.5λ0). The parameters of the inexact Newton method are the same as the ones adopted in the first case. The reconstructed distribution of the relative dielectric permittivity (final image) is shown in Figure 4. As can be seen, the object is correctly located, and the values of the dielectric properties are identified with a quite good accuracy.Figure 4
Reconstructed distribution of the relative dielectric permittivity. Double-layer circular cylinder buried in a lossy half space.Finally, the reconstruction of two separate circular cylinders have been considered. The radii of the two cylinders area1=0.2λ0 and a2=0.15λ0. They are buried with centers located at points r1obj=(0.2λ0,-0.5λ0) and r2obj=(-0.25λ0,-0.75λ0). The dielectric properties of the two scatterers are the following: ϵr1=5.6, ϵr2=8, and σ1=σ2=0.01 S/m. The signal-to-noise ratio is SNR=25dB. The parameters of the iterative inversion algorithm are the same as in the previous cases. The reconstructed distribution of the relative dielectric permittivity is shown in Figure 5. This image corresponds to the optimal solution. As can be seen, also in this quite complex case, the method is able to correctly reconstruct the investigated scene. In particular, the two objects are clearly identified and separated. Moreover, both shapes and dielectric properties of the two targets are estimated with quite good accuracies (although for the second cylinder the dielectric permittivity is slightly underestimated). In order to better assess the reconstruction accuracy, Figure 6 gives the profiles of the retrieved distribution of Figure 5 along some cuts parallel to the coordinate axes.Figure 5
Reconstructed distribution of the relative dielectric permittivity. Two separate circular cylinders buried in a lossy half space.Horizontal cuts ((a)y=-0.525λ0 and (b) y=-0.775λ0)) and vertical cuts ((c) x=0.325λ0 and (d) x=-0.325λ0)) of the reconstructed distribution of the relative dielectric permittivity (Figure 4). Two circular cylinders buried in a lossy half space.
(a)
(b)
(c)
(d)


## 5. Conclusion

In this paper, a previously proposed approach for the imaging of buried objects has been extended to inspect strong scatterers, for which the exact nonlinear equations of the inverse scattering problem must be considered. The developed method is based on a two-step inexact Newton method, which allows us to obtain a regularized solution of the (discretized) integral equations relating the samples of scattered electric field collected in the upper half space to the distributions of the dielectric parameters inside the test domain located in the ground. The reported preliminary results of numerical simulations seem to indicate the validity of the proposed approach, which needs further evaluation assessments, including an experimental validation, which will be the subject of future research activity in this area.



---

*Source: 637301-2012-07-29.xml*