# 3D Imaging of Dielectric Objects Buried under a Rough Surface by Using CSI

Authors: Evrim Tetik; Ibrahim Akduman
Journal: International Journal of Antennas and Propagation
 (2015)

---

## Abstract

A 3D scalar electromagnetic imaging of dielectric objects buried under a rough surface is presented. The problem has been treated as a 3D scalar problem for computational simplicity as a first step to the 3D vector problem. The complexity of the background in which the object is buried is simplified by obtaining Green’s function of its background, which consists of two homogeneous half-spaces, and a rough interface between them, by using Buried Object Approach (BOA). Green’s function of the two-part space with planar interface is obtained to be used in the process. Reconstruction of the location, shape, and constitutive parameters of the objects is achieved by Contrast Source Inversion (CSI) method with conjugate gradient. The scattered field data that is used in the inverse problem is obtained via both Method of Moments (MoM) and Comsol Multiphysics pressure acoustics model.


---

## Body

## 1. Introduction

Imaging of objects embedded in a layer or beneath a rough surface has been a popular subject of remote sensing due to its wide range of application areas such as nondestructive testing, mine detection, medical imaging, and geophysical or archaeological exploration. In most of the cases, the interface separating two layers is rough, and the object of interest is embedded in the inaccessible layer beneath the rough surface. In order to image the object, the area in which the object is presumed to be located is illuminated by a number of antennas, and the field scattered by the object and the roughness is measured by receiving antennas above the surface. The measured data is then used to obtain the unknown properties of the object, such as shape, location, dielectric permittivity, and conductivity.There have been a wide range of approaches introduced to investigate and reconstruct the object and the roughness. Interactions between the object and the rough surface are calculated in [1, 2]. The forward-backward (FB) method and propagation inside layer expansion (PILE) are used to calculate the scattering from an object above and below rough surface, respectively. In order to reconstruct the rough surface and the object buried beneath it simultaneously, an iterative method presented in [3] decreases the cost function by optimizing the boundary control parameters using least-squares and the Semianalytic Mode Matching (SAMM) forward model.In [4], a two-step method is used to detect and characterize a target buried beneath a rough surface. First, the target is detected by analyzing the frequency-averaged Wigner-Ville function in order to filter out rough surface scattering, and then it is characterized by applying an iterative solution derived from the Newton-Kantorovich algorithm to the Wigner-Ville function. In [5], surface and buried mines are detected by using mid-wavelength infrared (MWIR) images and image segmentation based on a wavelet thresholding algorithm. A hybrid numerical method based on the parallel genetic algorithm (GA) and the finite difference time domain (FDTD) is used to obtain the location and dimensions of 2D inhomogeneous objects buried in a lossy ground in [6].Surface impedance is obtained from the impedance boundary condition, by using the electric field and its normal derivative on the surface. By observing the surface impedance, the location of the dielectric objects buried beneath the surface is obtained in [7, 8]. Several methods also use standard, modified, or distorted Born iterative methods along with methods such as conjugate gradient, fast Fourier transform, and singular value decomposition (SVD) for imaging objects buried beneath a rough surface [9–12].In this study, objects buried beneath a rough surface in a 3D geometry are reconstructed by using CSI [13]. To this aim, Green’s function of the background is obtained by using BOA [14] in order to simplify the complexity of the geometry. A similar approach is presented in [15] for 2D case.Formulation of the problem is given in Section2. In Section 3, Green’s function of the background is obtained. The inverse problem is explained in Section 4, and numerical results of simulations and conclusion are given in Sections 5 and 6, respectively.Time dependenceexp-iωt has been used and suppressed throughout this problem.


## 2. Formulation of the Problem

In this paper, the problem of the imaging of an unknown object buried under a rough surface has been treated as a 3D scalar problem for computational simplicity as a first step to the 3D vector problem. The configuration of the problem can be seen in Figure1. In this configuration, the interface between the two half-spaces is a locally rough interface Γ=fx1,x2. The half-spaces x3>fx1,x2 and x3<fx1,x2 are composed of simple materials with constitutive parameters ε1,σ1 and ε2,σ2, respectively. In the lower half-space, an arbitrary shaped object with parameters εobj,σobj is located. Both half-spaces are assumed to be filled with nonmagnetic materials, whose magnetic permeabilities are equal to free-space permeability μ0.Figure 1
Cross-sectional geometry of the problem.The region in which the object is assumed to be located is illuminated by microwave point sources at pointszj, j=1,2,…,J, on surface T, and total electric field ujx=Enx1,x2,zj is measured on surface S for each illumination. The incident scalar electric field function for each illumination, which is time harmonic with angular frequency ω, is Eix1,x2,zj=eikr/4πr, where r represents the distance between the source and observation points. The inverse scattering problem here is reconstructing the object’s location, shape, and constitutive parameters by using the total electric field ujx=Enx1,x2,zj, n∈S, measured on the surface S.Letuj0x be the background field, the total field for the jth illumination in the absence of the object. It can be computed by using BOA [14], which treats the roughness as a series of objects embedded in upper and lower half-spaces that can be seen in Figure 2. Therefore, if the field scattered by the object is represented by ujSx=uj0x-ujx, the inverse problem can be expressed as the reconstruction of χx in the following system of integral equations:(1)ujx=uj0x+kb2∭DGbx,yχyujydy,x∈D,ujSx=kb2∭DGbx,yχyujydy,x∈S,which are theobject anddata equations [13], respectively. In these equations, Gbx,y is Green’s function of the background system which consists of the two half-spaces and the rough interface that separates them. The function χx is the object or contrast function which is defined by (2)χx=k2xkb2x-1,where kx and kbx are the wavenumbers defined by(3)kb2x=k12=ω2ε1μ0+iωσ1μ0,x3>fx1,x2k22=ω2ε2μ0+iωσ2μ0,x3<fx1,x2,k2x=ω2εobjμ0+iωσobjμ0,x3∈Ckb2x,x3∉C.It can be seen in (2) that the contrast function χx is zero outside the object, and its real and imaginary parts are related to εobj and σobj, respectively.Figure 2
Cross-sectional representation of the series of objects which construct the roughness according to BOA.


## 3. Green’s Function of the Background Medium

BOA is used to obtainGbx,y, which is Green’s function of the background, for solving the inverse scattering problem given by (1). This approach considers the roughness as a series of scatterers whose parameters differ according to their position. The parameters of the sections B1,B3,…,B2N-1 that lie in the x3<0 half-space and the sections B2,B4,…,B2N that lie in the x3>0 region are ε1,σ1 and ε2,σ2, respectively.G
b
x
,
y, Green’s function of the background, can be expressed as the sum of the contribution of its two components: (4)Gbx,y=GSx,y+G0x,y,where G0x,y is Green’s function of the two-part space with a planar interface and GSx,y is the contribution of the roughness. GSx,y can also be considered as the field scattered from the objects that form the roughness due to a point source with unit strength which is located at y, and it satisfies the integral equation(5)I-LGSx,y=LG0x,y,where the integral operator L is defined by(6)LGx,y=∭DG0x,zkb2z-k02zGx,zdz.k0x is the wavenumber of the two-part space with a planar interface, defined by(7)k02x=k12,x3>0k22,x3<0.Equation (5) can be solved for GSx,y by applying the forward solution method given in [16]. G0x,y, which is needed to solve this equation, can be given as(8)G0x,y=eik1r-r′4πr-r′+14π∫0∞R12γ1kρJ0kρρ-ρ′e-γ1z+z′dkρ,x3>0,y3>014π∫0∞T12γ1kρJ0kρρ-ρ′e-γ1z′+γ2zdkρ,x3<0,y3>014π∫0∞T21γ2kρJ0kρρ-ρ′e-γ1z+γ2z′dkρ,x3>0,y3<0eik2r-r′4πr-r′+14π∫0∞R21γ2kρJ0kρρ-ρ′eγ2z+z′dkρ,x3<0,y3<0,where(9)R12=γ1-γ2γ1+γ2,T12=2γ1γ1+γ2,R21=γ2-γ1γ1+γ2,T21=2γ2γ1+γ2,(10)γ1kρ=kρ2-k12,γ2kρ=kρ2-k22.The square root functions in (10) are defined on complex kρ plane under the conditions(11)γn0=-ikn,n=1,2.


## 4. Inverse Problem

If we rewrite (1) in a symbolic form, we obtain the object equation (12)uj=uj0+GjDχjuj,x∈Dand the data equation(13)fj=GjSχjuj,x∈S,where the integral operators GD and GS are defined by(14)GD,SΨx=∭DGbx,ykb2yΨydy,x∈D,S.Because the product of the contrast function and the fields are seen in the integral equation system of the problem, the quantity(15)wj=χjuj,which is known as a contrast source, can be replaced in (12) and (13) to obtain(16)uj=uj0+GjDwj,x∈D,fj=GjSwj,x∈S.Using CSI method, the inverse problem becomes the minimization of the cost functional defined as(17)F=∑jfj-GjSwjS2∑jfjS2+∑jχuj0-wj+χjGjDwjD2∑jχjuj0D2,where ·S2 and ·D2 denote the norms on L2S and L2D, respectively. During the iteration process, the contrast sources wj and the object function χj are alternately constructed by using conjugate gradient method such that they minimize the whole cost functional until it reaches a predetermined value. The explicit expressions and the iteration steps for the CSI method can be found in [13].


## 5. Numerical Results

Throughout the simulations, the regionx>0 is considered to be free-space, and x<0 dry soil with relative dielectric permittivity εr2=3.6 and conductivity σ2=10-4. 25 equally spaced point sources are located on a square surface, which is defined by -0.5λ2<x1<0.5λ2, -0.5λ2<x2<0.5λ2, and x3=0.57λ2, and the scattered field is measured on equally spaced positions on a surface which is defined by -0.57λ2<x1<0.57λ2, -0.57λ2<x2<0.57λ2, and x3=0.38λ2, for each illumination, where λ2 is the wavelength of the soil. The reconstruction domain is divided by cubical cells of size λ2/20 on each side. The scattered field data are obtained synthetically by using Comsol Multiphysics pressure acoustics with attenuation model, or MoM, and 5% random noise; ηuS,jei2πrd is added, where η is the noise level and rd is a uniformly distributed random variable between 0 and 1, and therefore the corresponding signal to noise ratio is SNR=-20log10η. All simulations were carried out at a frequency of 300 MHz.According to Comsol Multiphysics pressure acoustics with attenuation model, the total acoustic pressure functionpt created by a point source with strength Q satisfies the differential equation(18)∇·1ρ∇pt-q-k2ptρ=Q,k=ωc-iα,where k is the wavenumber, and ρ, c, and α are the density, speed of sound, and attenuation, respectively.In this model, it can be seen from (18) that pt corresponds to scalar electric field of the problem at hand. If ρ is matched to the relative magnetic permeability, and c and α are selected so that the real and imaginary parts of the wavenumbers match, this pressure acoustics model can be used to simulate scalar electromagnetic model.MoM has been used to obtain Green’s function of the background in (4) and (5). A cell size of λ2/20×λ2/20×λ2/20 is selected.For the first simulation, a simple roughness, which is located in the domain defined by-λ2<x1<λ2, -λ2<x2<λ2, and -0.28λ2<x3<0.28λ2, is selected. The profile of the roughness can be seen in Figure 3(a). A spherical object with relative dielectric permittivity εr=5 and conductivity σ=0.1 and a radius of 0.1λ2 is centered at (0,0,0.28λ2). The scattered field was measured at 400 points. The reconstruction was carried out in a domain defined by -0.57λ2<x1<0.57λ2, -0.57λ2<x2<0.57λ2, and -0.76λ2<x3<0. The iteration is stopped at 500. The original and reconstructed relative dielectric permittivity are shown in Figure 3(b), and the original and reconstructed conductivity are shown in Figure 3(c). It has been observed that the shape, location, and the relative dielectric permittivity of the object are obtained quite well. However, the conductivity of the object is underestimated.Figure 3
(a) Rough surface; (b) original and reconstructed relative dielectric permittivity; (c) original and reconstructed conductivity.
(a)
(b)
(c)A second simulation has been performed with the synthetic data obtained by solving the forward problem using MoM for the same configuration. The original and reconstructed values of relative dielectric permittivity and conductivity of the object can be seen in Figures4(a) and 4(b), respectively. This simulation shows that, using MoM, instead of Comsol Multiphysics, although the object appears to be slightly above its actual position, its shape and location are still clear, and the reconstructed value of the relative permittivity improves; however, the reconstructed conductivity still does not improve significantly.Figure 4
(a) Original and reconstructed relative dielectric permittivity; (b) original and reconstructed conductivity.
(a)
(b)A larger and more complex roughness, which is located in the domain defined by-2λ2<x1<2λ2, -2λ2<x2<2λ2, and -0.38λ2<x3<0.38λ2, and that can be seen in Figure 5(a), is selected for the third simulation. The same object used in the previous examples is centered at the same location in order to investigate the effect of the roughness. The scattered field was measured at 625 points. After 500 iterations, the original and reconstructed relative dielectric permittivity and conductivity can be seen in Figures 5(b) and 5(c), respectively.Figure 5
(a) Rough surface; (b) original and reconstructed relative dielectric permittivity; (c) original and reconstructed conductivity.
(a)
(b)
(c)It is observed that although the location and the value of the relative permittivity of the object are still well determined, the shape is not as clear, indicating that the method is sensitive to the complexity of the structure, and the higher the complexity is, the lower the resolution is with the given configuration of the reconstruction.For the fourth simulation, two spherical objects with a radius of0.1λ2 are assumed to be buried under the rough surface that can be seen in Figure 5(a). The center of the object with relative dielectric permittivity εr1=5 and conductivity σ1=0.1 is located at 0.3λ2,0,-0.32λ2, and the center of the object with relative dielectric permittivity εr2=7 and conductivity σ2=0.2 is located at 0.3λ2,0,-0.2λ2. The scattered field was measured at 900 points. After 500 iterations, the original and reconstructed relative dielectric permittivity and conductivity can be seen in Figures 6(a) and 6(b), respectively.Figure 6
(a) Original and reconstructed relative dielectric permittivity; (b) original and reconstructed conductivity.
(a)
(b)It can be seen in Figure6(a) that the value of the relative dielectric permittivity is determined, and although the value of the conductivity in Figure 6(b) is underestimated, the shape and location of the objects are quite clear. This simulation suggests that when the complexity of the geometry rises, the clarity and the precision of the reconstructions decrease.As it can be seen in some of the simulations, the position of the object may appear slightly different from its actual position although its shape and constructive parameters remain in an acceptable range. In such cases, measuring the error based on pixel-by-pixel comparison may be misleading. Therefore, in order to prevent this, instead of a measure of the errors, the exact and the reconstructed profiles are given for visual interpretation.


## 6. Conclusion

A 3D scalar electromagnetic imaging of dielectric objects buried under a rough surface is presented. BOA is used to obtain Green’s function of the background, which includes the contribution of the roughness in order to be able to separate the contribution of the object from its background. Green’s function of the two-part space with a planar interface is obtained to be used in the process. CSI method with conjugate gradient is used to reconstruct the location, shape, and constitutive parameters of the object.The numerical simulations with data obtained by using both MoM and Comsol Multiphysics pressure acoustics model show that the method shows satisfactory results although, due to the computational weight of the 3D structure, larger cell size, relative to the complexity of the problem, has been used.



---

*Source: 179304-2015-10-01.xml*