# Imaging Internal Defects with Synthetic and Experimental Data

Authors: Hongwei Zhou; Guanghui Hu; Ling Ma
Journal: International Journal of Antennas and Propagation
 (2019)

---

## Abstract

This work concerns an inverse time-dependent electromagnetic scattering problem of imaging internal defects in a homogeneous isotropic medium. The position and cross section of the defects are detected by transient electromagnetic pulses in the case of TE polarization. We apply the Kirchhoff migration scheme to locate the position of small objects from both synthetic and experimental data. The multiple-input-multiple-out scheme is used to recover extended scatterers from the data generated by the software GprMax. Numerical experiments show that the Kirchhoff migration method is not only efficient but also robust with respect to polluted data at high noise levels. Experimental results show good quantitative agreement with numerical simulations.


---

## Body

## 1. Introduction

Inverse scattering is to recover physical and geometrical information of inaccessible objects from scattered fields measured outside. It has been one of the most challenging problems with considerable practical applications in many areas of technology such as nondestructive evaluation, subsurface and ground-penetrating radar, geophysical remote sensing, medical imaging, seismology, and target identification; see [1–5]. One of the most prominent features of inverse scattering problem is its noninvasiveness, along with the affordability due to cheap nonionizing sensors. However, solving inverse scattering problems is difficult due to the inherent ill-posedness and nonlinearity. Small variation in the measured data can lead to large errors in the reconstruction of the scatterer, unless regularization methods are used. Extensive investigation has been carried out and a variety of inversion algorithms have been proposed.In the time-harmonic regime, Kirsch and Kress [6, 7], Colton and Monk [8, 9], and Angrll, Kleiman, and Roach [10] theoretically separate the nonlinearity and the ill-posedness of inverse obstacle problems, giving rise to the decomposition method. Muhammad [11], Liang Yang [12], and Hui Het al. [13] employed iterative methods to solve inverse obstacle problems. In these optimization-based iterative approaches, an efficient forward solver is needed for each iteration, and good a priori information might be required in order to choose an initial guess that ensures numerical convergence. Noniterative sampling methods were intensively studied over the last twenty years, for instance, linear sampling method [14], factorization method [15], enclosure method [16], and singular point source method [17]. The key ingredient of the sampling approach is to design an appropriate indicator function from measured data for characterizing the region occupied by the scatterer. Forward solvers are not needed in the process of inversion. The above-mentioned approaches are mostly applicable to measured data irradiated at a fixed energy with many incoming directions. A recursive linearization method with multifrequency data was investigated in [18–22].There exist a number of reconstruction algorithms in the time domain with dynamic measurement data. Here we mention the time reversal techniques [23, 24], the reversed time migration [25], and the boundary control method [26]. For time-dependent sampling type methods, we refer to the point source method [27], the enclosure method [28], and the total focusing method (TFM) arising from nondestructive evaluation [29, 30]. In one of the author's previous works [31], the TFM (which is also known as the Kirchhoff migration approach used in geophysics; see, e.g., [32]) was examined for imaging acoustically extended scatterers in two dimensions and the indicator behavior was mathematically analyzed. This work concerns imaging internal defects in a homogeneous isotropic medium. It is shown that the TFM is not only efficient but also robust to measurement noise. The aim of this paper is to test the TFM for inverse electromagnetic scattering problem of imaging perfectly conducting cylinders buried in a half-space homogeneous isotropic background medium. Such kind of inverse problems has many applications in searching for internal defects with radar techniques like ground-penetrating radar (GPR); see, e.g., [5, 33, 34]. We note that the light speed is much faster than the sound speed, leading to additional difficulties in computational simulation [35]. In this paper, we make use of the open software GprMax [35, 36] to generate the forward scattering data. The GprMax is an electromagnetic simulation tool based on the Finite-Difference Time-Domain (FDTD) approach coupled with the perfectly matched layer (PML) technique. In our experiments, the real data are gained by GPR.In this paper, we apply the TFM (or Kirchhoff migration scheme) to inverse electromagnetic scattering problems. In comparison to other algorithms, the implementation of TFM is quite simple, since the image is formed through a superposition of the scattered signals irradiated by each transducer. In fact, our imaging functions are explicit and involve only integral calculations on the measurement surface. Hence, our inversion algorithm is totally “direct”. This also explains why this method is very robust with respect to measurement noise at high levels. Although the inverse scattering problem is difficult due to the inherent ill-posedness and nonlinearity, we do not need to solve the ill-posed and nonlinear problems, even without approximation and iteration. In this work, experiments with synthetic and experimental data are conducted to demonstrate the feasibility and applicability of this approach in the TE polarization case of inverse electromagnetic scattering problems.The remaining part of this paper is organized as follows. In Section2 we rigorously formulate the forward and inverse electromagnetic scattering problems. Section 3 is devoted to the description and implementation of the TFM for imaging small scatterers from both synthetic and real data. Extended scatterers will be reconstructed in Section 4 by using GprMax simulated data. Conclusion follows in Section 5.


## 2. Problem Setting

Consider the electromagnetic wave propagation in a homogeneous isotropic medium. This background medium can be characterized by the electric permittivityε, magnetic permeability μ, and electric conductivity σ, all of which are assumed to be constant. The propagation of the electric and magnetic fields (E,H) is governed by the Maxwell system.(1)curlE+μ∂H∂t=0(2)curlH-ε∂E∂t=σEAssume that D=Ω×R is an infinitely long perfectly conducting cylinder buried in the lower half-space x2<0, which remains invariant in x3-direction; see Figure 1(a). The medium in the exterior De:=R3∖D- ofD is assumed to be dielectric, which means that σ=0. The cross section of D in the ox1x2-plane is denoted by Ω⊂R-2≔x=x1,x2:x2<0. Incident cylindrical source waves emitted from the surface Γ⊂R3∖D-∩x3=0 will be utilized for the purpose of detecting the position and shape of Ω. Throughout the paper, we assume that the incident wave Eχix,t;y with x=x,x3, y=y,y3∈R3, t>0 is a filamentary current generated at the source location y∈Γ with a temporal pulse signal χ:R→R. The temporal function is supposed to have a compact support χt∈0,T0 for some T0>0. Moreover, we consider the Transverse Electric (TE) case of the monopole with the polarization direction e3=0,0,1. Then, the incident electric and magnetic waves Ei,Hi are governed by the system(3)curlEix,t;y+μ∂Hix,t;y∂t=0curlHix,t;y-ε∂Eix,t;y∂t=-δx-yχte3for all x≠y, x∈De, and t>0, where δ is the Dirac distribution. Denote by Es,Hs the scattered fields and write the total field as E=Ei+Es, H=Hi+Hs. Then, we have(4)curlEx,t;y+μ∂Hx,t;y∂t=0inDe∖y×R+,curlHx,t;y-ε∂Ex,t;y∂t=-δx-yχte3inDe∖y×R+,v×Ex,t;y=0on∂D×R+,together with the zero initial condition at t=0 in R3. Here v=v1,v2,v3 is the unit outward normal to the boundary ∂D. Eliminating the magnetic field, we arrive at(5)∂2Ex,t;y∂t2=-1εμcurlcurlEx,t;y+1εδx-y∂tχte3,x≠y.In the TM polarization case, the electric field takes the form Ex=0,0,ux. Hence, the perfectly conducting boundary condition can be written as(6)v×E=ijkv1v2000u=v2u,-v1u,0on∂D×R+,implying that u=0 on ∂D×R+. Using(7)curlcurlEx=-ΔEx+graddivEx,we deduce the reduced wave equation from the Maxwell system (5) as follows:(8)∂2ux,t;y∂t2=1εμΔux,t;y+1εδx-y∂tχt,x∈R2∖Ω-∪y×R+.In the particular case that χt=sin⁡ωt is a 2π-periodic function, where ω denotes the frequency, the incident field uin can be explicitly represented as(9)uinx,t;y=μωRei4H01ωcx-yeiωt,x≠y.Here H01 is the Hankel function of first kind of order zero; see [37, Chapter 3.4].Figure 1
Geometrical schematic.
(a)
Geometrical model (b)
Sample grid subdivisionIn this paper we shall consider the following inverse problem:(IP): determine the position and shape of∂Ω from knowledge of us(x,t;y) irradiated by incident dipoles located at y∈Γ.Throughout this paper, we use the temporal function of the form(10)χt=-A2π2f2t-2f2-1e-π2f2t-2/f2,t>0,which is known as the Ricker pulse. Here f>0 is the center frequency and A>0 0 is the amplitude. The incident pulse function χ(t) with f=6 GHz is depicted in Figure 2.Figure 2
The Ricker pulse withf = 6 GHz and A = 1.We shall consider two cases:(i)Ω∈R-2 is a small object (compared to the incident wavelength), Γ=(x1,b):a1<x1<a2 is a finite line segment, and the data are given by us(x,t;x) for all (x,t)∈Γ×(0,T).(ii)Ω is an extended scatterer, Γ=x∈R2:x=R for some R>0 such that Ω⊂BR≔x∈R2:x<R, and the data are given by us(x,t;y) for all (x,t)∈Γ×(0,T) irradiated by many dipoles emitted from Γ.Note that in case (i), the positions of the receivers and transmitters are identical, leading to the so-called inverse back scattering problem, whereas the case (ii) corresponds to a multiple-input-multiple-output (MIMO) radar system. The aim of this paper is to apply the total focusing method to the above inverse electromagnetic scattering problems with synthetic and experimental data. The TFM is sometimes described as the “gold standard” in the classical beamforming. Both TFM and the Kirchhoff migration form the image with the superposition of the scattered signals irradiated by each transmitter.


## 3. Imaging Small Scatterers

In this section we suppose that a single metal stick is buried in sand below the ground planex2=b. The cross section Ω of the metal is so small compared to the wavelength such that Ω can be regarded as a point-like object. The electric dipoles are generated on part of the ground plane Γ=x1,b:a1<x1<a2 for some aj>0j=1,2. The Kirchhoff migration imaging scheme for inverse back scattering consists of the following three steps:Firstly, we take the searching area as a rectangular domain in the lower half-spaceR-2 and divide it into an M×N grid. The grid points, as shown in Figure 1(b), will be referred to as sampling points z∈R-2. We shall design an imaging function defined on these grid points.Secondly, the filamentary current is transmitted from the source position{x(k)=(k(a2-a1)/N,0)} and the scattered signals are also recorded at the same position x(k), which we denote by uks(t)≔us(x(k),t;x(k)) for k=1,2,…,N. Note that us(x,t;y)=u(x,t;y)-ui(x,t;y) stands for the scattered wave field. Calculate the flight time of the signal emitted from x(k) to a sampling point z and then returned back to the receiver x(k), given by(11)tkz=2xk-zc,k=1,2,…,N.Here c=1/εμ>0 denotes the light speed in the background medium. In particular, we have c=c0ε in terms of the light speed c0 in the vacuum.Thirdly, collect all signalsuks for k=1,2,…,N and plot the indicator function(12)Iz=∑k=1Nukstkz2>0.Here we have assumed that the terminal T is large enough. By the mathematical analysis performed in [21], the local maximizers of I can represent the location of the small scatterer.


### 3.1. Imaging with Synthetic Data

To implement the Kirchhoff migration scheme, we assume that the interface between the air and the sand is given byx1,b:x1∈R with =-0.1. The dielectric properties of the sand are considered as 4. We choose N=41 receivers and transmitters uniformly lying between 0.1,-0.1 and 0.1,-0.3. The cross section of the metal stick is supposed to be the circle of radius 0.009 centered at 0.2,-0.167; see Figure 3 for the locations of the metal stick in the ox1x2-plane and the first transmitter (receiver). The center frequency of Ricker waves is uniformly taken as f=2.6GHz and the amplitude is set to be A=1. The terminal time is taken as(13)T=2hc+1f,where h>0 denotes the travel distance of the radar. In our simulations we take T=7×10-9s.Figure 3
Locations of the first transmitter/receiver and the target.The scattered data were generated by GprMax based on the FDTD method (see [20, 38]). The unbounded exterior domain is truncated by an absorbing boundary condition. The mesh of the forward solver is successively refined until the relative error of the successive measured scattered data is below 0.1%.For thek-th transmitter, the scattered data uks are recorded at 1063 time points given by(14)tk=k1063T,k=1,2,…,1063.Linear interpolation is adopted to obtain the value of ukst for an arbitrary t∈0,T. In Figure 4, we show the signals of the incident electric field and the total fields irradiated by the first transmitter. The incident field obviously has a launch waveform, as shown in Figure 4(a). The significant changes before 0.2ns are due to the transmitter itself rather than the reflected wave field from the metal stick. In Figure 4(b), one can observe the changes of the signal at about 0.4ns caused by the presence of the object. Subtracting the incident field from the total field, one can obtain the scattered signal, which will be utilized to locate the position of the stick.Figure 4
Electric signals generated by the first transmitter.
(a)
Incident field (b)
Total fieldForx=x1,x2∈R2 with x1∈0.1+kh,0.1+k+1h, h=a2-a1/N=0.2/268, we first plot the function(15)Bx≔uks2x2-0.1c,-0.4<x2<-0.1,which forms the B-scan image of the object. Note that only the information of the signal of uks is involved in B-scan for x1∈0.1+kh,0.1+k+1h, whereas all signals are required in our indicator equation (15). Figure 5(a) shows the B-scan image of the small object located at (0.2, 0.167). The x1-coordinate of the location corresponds to the trace number 100. In Figure 5(b), one can see an approximate location of the target using our indicator (15). Numerical experiments show that the Kirchhoff migration method is not only efficient but also robust with respect to polluted data at high noise levels (Figure 6). Note that in this paper the scattered near-field data is polluted by(16)uδscx,t=uscx,t+δζuscx,t,where δ is the noise ratio and ζ is the standard normal distribution.Figure 5
Imaging a small scatterer using synthetic data without noise.
(a)
B-scan image formed by the indicator equation (15) (b)
The Kirchhoff migration imageFigure 6
Imaging a small scatterer using polluted synthetic data at the noise levelδ=100%.
(a)
B-scan image formed by the indicator equation (15) (b)
The Kirchhoff migration image


### 3.2. Imaging with Experimental Data

In our experiment settings, a metal stick is detected in a sand box by a GSSI Mini HR hand-held ground-penetrating radar; see Figure7 for the experiment site. The size of the sand box is 0.45 m×0.3m×0.2m, the thickness of the wooden shell is 0.015m, the dielectric properties of the sand are 4, and the dielectric properties of the wooden box are considered approximately the same as the sand. For electric dipoles with TE polarization, the three-dimensional Maxwell equation can be reduced to a two-dimensional model for scalar wave equations. We set up a two-dimensional coordinate system and take the initial position of the source and receiver of the radar at (0.1, 0.1). The metal stick is centered at (0.2, 0.167) with the radius 0.009. The radar rolls along the straight line x2 = 0.1 on the wooden shell and moves forward in the positive x1-direction at an even pace. In this process there are totally 268 electric dipoles excited by the radar before it reaches the ending point (0.25, 0.1). The temporal function of the incident wave is a Ricker wavelet excited at the frequency 2.6GHz with the amplitude 1.Figure 7
Experiment site.The B-scan image formed by our experimental data is plotted with Matlab R2013b, as shown in Figure8(a). Figure 8(b) shows the image using the indicator (12), where the position of the target is precisely located. In comparison with the counterpart from synthetic data (see Figure 5(b)), we find that the inversion of the depth of the stick contains relatively large errors.Figure 8
Imaging with experimental data.
(a)
B-scan image (b)
Imaging by using the Kirchhoff migrationNext, we want to find two metal sticks with different positions. The centers of the cross section in theox1x2-plane are given by (0.15, 0.15) and (0.25, 0.17), respectively. The other parameters in our experiment are the same as before. See Figure 9(b) for the experiment site and Figure 9(a) for the location of the two sticks.Figure 9
Locating two small scatterers buried in sand.
(a)
Locations of the first transmitter/receiver and the two small targets (b)
Experiment siteThe B-scan images with synthetic and experimental data are shown in Figures10(a) and 10(b), respectively. The final inversion images based on the Kirchhoff migration are illustrated in Figures 11(a) and 11(b), respectively.Figure 10
Comparison of the B-scan images with the synthetic data at noise level 20% (a) and experimental data (b).
(a)
Synthetic data at the noise level 20% (b)
Experimental dataFigure 11
Comparison of the images formed by the Kirchhoff migration scheme.
(a)
Synthetic data at the noise level 20% (b)
Experimental dataAs shown in Figure10, the B-scan results in both simulation and experiment can indicate a local peak in the position of the metal sticks, verifying the correctness of the inversion scheme. In Figure 11, the local maximum values of the indicator function are obtained approximately at (0.15, 0.15) and (0.25, 0.17), which agree well with the positions of the targets in our model. From the experiment image, it can be concluded that the contour image and the position of the left metal stick agree well with the true scatterers, but the size of the right one is not well reconstructed. Further, there is an untrustworthy imaging area between the two metal sticks, maybe due to the offset of the two maximizes nearby.


## 3.1. Imaging with Synthetic Data

To implement the Kirchhoff migration scheme, we assume that the interface between the air and the sand is given byx1,b:x1∈R with =-0.1. The dielectric properties of the sand are considered as 4. We choose N=41 receivers and transmitters uniformly lying between 0.1,-0.1 and 0.1,-0.3. The cross section of the metal stick is supposed to be the circle of radius 0.009 centered at 0.2,-0.167; see Figure 3 for the locations of the metal stick in the ox1x2-plane and the first transmitter (receiver). The center frequency of Ricker waves is uniformly taken as f=2.6GHz and the amplitude is set to be A=1. The terminal time is taken as(13)T=2hc+1f,where h>0 denotes the travel distance of the radar. In our simulations we take T=7×10-9s.Figure 3
Locations of the first transmitter/receiver and the target.The scattered data were generated by GprMax based on the FDTD method (see [20, 38]). The unbounded exterior domain is truncated by an absorbing boundary condition. The mesh of the forward solver is successively refined until the relative error of the successive measured scattered data is below 0.1%.For thek-th transmitter, the scattered data uks are recorded at 1063 time points given by(14)tk=k1063T,k=1,2,…,1063.Linear interpolation is adopted to obtain the value of ukst for an arbitrary t∈0,T. In Figure 4, we show the signals of the incident electric field and the total fields irradiated by the first transmitter. The incident field obviously has a launch waveform, as shown in Figure 4(a). The significant changes before 0.2ns are due to the transmitter itself rather than the reflected wave field from the metal stick. In Figure 4(b), one can observe the changes of the signal at about 0.4ns caused by the presence of the object. Subtracting the incident field from the total field, one can obtain the scattered signal, which will be utilized to locate the position of the stick.Figure 4
Electric signals generated by the first transmitter.
(a)
Incident field (b)
Total fieldForx=x1,x2∈R2 with x1∈0.1+kh,0.1+k+1h, h=a2-a1/N=0.2/268, we first plot the function(15)Bx≔uks2x2-0.1c,-0.4<x2<-0.1,which forms the B-scan image of the object. Note that only the information of the signal of uks is involved in B-scan for x1∈0.1+kh,0.1+k+1h, whereas all signals are required in our indicator equation (15). Figure 5(a) shows the B-scan image of the small object located at (0.2, 0.167). The x1-coordinate of the location corresponds to the trace number 100. In Figure 5(b), one can see an approximate location of the target using our indicator (15). Numerical experiments show that the Kirchhoff migration method is not only efficient but also robust with respect to polluted data at high noise levels (Figure 6). Note that in this paper the scattered near-field data is polluted by(16)uδscx,t=uscx,t+δζuscx,t,where δ is the noise ratio and ζ is the standard normal distribution.Figure 5
Imaging a small scatterer using synthetic data without noise.
(a)
B-scan image formed by the indicator equation (15) (b)
The Kirchhoff migration imageFigure 6
Imaging a small scatterer using polluted synthetic data at the noise levelδ=100%.
(a)
B-scan image formed by the indicator equation (15) (b)
The Kirchhoff migration image


## 3.2. Imaging with Experimental Data

In our experiment settings, a metal stick is detected in a sand box by a GSSI Mini HR hand-held ground-penetrating radar; see Figure7 for the experiment site. The size of the sand box is 0.45 m×0.3m×0.2m, the thickness of the wooden shell is 0.015m, the dielectric properties of the sand are 4, and the dielectric properties of the wooden box are considered approximately the same as the sand. For electric dipoles with TE polarization, the three-dimensional Maxwell equation can be reduced to a two-dimensional model for scalar wave equations. We set up a two-dimensional coordinate system and take the initial position of the source and receiver of the radar at (0.1, 0.1). The metal stick is centered at (0.2, 0.167) with the radius 0.009. The radar rolls along the straight line x2 = 0.1 on the wooden shell and moves forward in the positive x1-direction at an even pace. In this process there are totally 268 electric dipoles excited by the radar before it reaches the ending point (0.25, 0.1). The temporal function of the incident wave is a Ricker wavelet excited at the frequency 2.6GHz with the amplitude 1.Figure 7
Experiment site.The B-scan image formed by our experimental data is plotted with Matlab R2013b, as shown in Figure8(a). Figure 8(b) shows the image using the indicator (12), where the position of the target is precisely located. In comparison with the counterpart from synthetic data (see Figure 5(b)), we find that the inversion of the depth of the stick contains relatively large errors.Figure 8
Imaging with experimental data.
(a)
B-scan image (b)
Imaging by using the Kirchhoff migrationNext, we want to find two metal sticks with different positions. The centers of the cross section in theox1x2-plane are given by (0.15, 0.15) and (0.25, 0.17), respectively. The other parameters in our experiment are the same as before. See Figure 9(b) for the experiment site and Figure 9(a) for the location of the two sticks.Figure 9
Locating two small scatterers buried in sand.
(a)
Locations of the first transmitter/receiver and the two small targets (b)
Experiment siteThe B-scan images with synthetic and experimental data are shown in Figures10(a) and 10(b), respectively. The final inversion images based on the Kirchhoff migration are illustrated in Figures 11(a) and 11(b), respectively.Figure 10
Comparison of the B-scan images with the synthetic data at noise level 20% (a) and experimental data (b).
(a)
Synthetic data at the noise level 20% (b)
Experimental dataFigure 11
Comparison of the images formed by the Kirchhoff migration scheme.
(a)
Synthetic data at the noise level 20% (b)
Experimental dataAs shown in Figure10, the B-scan results in both simulation and experiment can indicate a local peak in the position of the metal sticks, verifying the correctness of the inversion scheme. In Figure 11, the local maximum values of the indicator function are obtained approximately at (0.15, 0.15) and (0.25, 0.17), which agree well with the positions of the targets in our model. From the experiment image, it can be concluded that the contour image and the position of the left metal stick agree well with the true scatterers, but the size of the right one is not well reconstructed. Further, there is an untrustworthy imaging area between the two metal sticks, maybe due to the offset of the two maximizes nearby.


## 4. Imaging Extended Scatterers

The main purpose of this section is to test the Kirchhoff migration for imaging extended scatterers from synthetic data of several incident electric dipoles. As in the previous section, the forward data are obtained by the software GprMax. We consider the second case of our inverse electromagnetic scattering problem formulated at the end of Section2, where the source positions are uniformly distributed on the circle x-O=R for some R>0, where O=(o1,o2)∈R2. Denote by M>0 the number of incident dipoles. The Cartesian coordinates of the m-th dipole can be formulated as(17)ym=O+Rcos⁡θm,sin⁡θm,θm≔2mπM,m=1,2,…,M.In all examples we use N=16 receivers x(n)(n=1,2,…,16) equally spaced on x-O=R. For notational convenience, we use un.ms(t)=un,ms(x(m),t;x(n)) to represent the scattered signal excited by the n-th transmitter recorded at the m-th receiver. The flight time of an impulsive signal from y(n) to a sampling point z=(z1,z2) and then back to xm is then given by(18)tn,mz=yn-z+z-xmc,n,m=1,2,…,16.In this case the imaging function based on Kirchhoff migration scheme can be formulated as(19)Iz=∑n=1N∑m=1Mun,mstn,mz2>0.In all of our numerical examples to be reported below, we set O=0.35,0.35. The sampling region is set to be the rectangular domain 0,0.7×0,0.7. The homogeneous background medium is set to be air, with the dielectric constant ε=1 and the electric conductivity σ=0. Unless otherwise stated, the excitation frequency of a Ricker wave is taken as f=10GHz. The amplitude is set to be A=1 and the time window is T=4ns.In the first example, we excite only one electric dipole located aty4 and measure the scattered field at xn for n=1,2,…,16. Three types of scatterers will be recovered: two small circular metals, a square-like metal, and an L-shaped metal; see Figures 12(a), 12(c), and 12(e). The last two targets will be referred to as extended scatterers, since their size is much bigger than the incident wavelength. The circular metals are centered at A0.25,0.35 and B0.45,0.35 with the radius 0.01 m. They can be treated as point-like scatterers, because their radiuses are both much smaller than the wavelength λ=0.3m corresponding to f=10GHz. The coordinates of the corners of the square-like metal are given by A0.275,0.425, B0.425,0.425, C0.425,0.275, and B0.275,0.275, and those for the L-shaped metal are A0.275,0.475, B0.35,0.475, C0.35,0.35, D0.475,0.35, E0.475,0.275, and F0.275,0.275. Figures 12(b), 12(d), and 12(f) show the reconstructions from the data of one dipole only. It can be seen that a single source can be used to capture the position of point scatterers and the illuminated boundary of the square-like and L-shaped scatterers. However, the entire shape of extended objects cannot be well reconstructed.Figure 12
Imaging small and extended metals with a single incoming dipole.
(a)
Two small objects (b)
Location of the two small objects (c)
A square-like metal (d)
Imaging of the square-like metal (e)
An L-shaped metal (f)
Imaging of the L-shaped metalIn the second example, we increase the number of incident dipoles for reconstructing the entire shape; see Figure13. The data of four incident waves without noise can generally produce a full-range image of an extended scatterer. However, there are bigger background disturbances and the scatterer cannot be accurately recovered. Using eight and sixteen incoming sources can image the outline of the scatterer. In Figures 13(e) and 13(f), the maximum values of the indicator are almost distributed on the contour of the object and the imaging is better than the case of eight dipoles. This means that sixteen sources are capable of imaging the scatterer, and therefore it is not necessary to use additional sources. Hence, the data of multiple point sources in different directions need to be measured in order to better capture the shape of the target.Figure 13
Imaging a square-like object (left) and an L-shape object (right) byM transmitters and N = 16 receivers.
(a)
M = 4 (b)
M = 4 (c)
M = 8 (d)
M = 8 (e)
M = 16 (f)
M = 16In the third example, we discuss the sensitivity of the imaging quality to the center frequency of incident pulse functions. We useM=16 incident dipoles andN=64 receivers. As shown in Figure 14, the shape of an extended object can be imaged at various frequencies from 1GHz to 10GHz. However, at lower frequencies the target cannot be accurately recovered. Increasing the excitation frequency of the source (which means a shorter wavelength) may lead to an image with higher resolutions.Figure 14
Imaging an L-shape object byM=16 transmitters and N=64 receivers (f is the excitation frequency of the sources).
(a)
f=1GHz (b)
f=2GHz (c)
f=3GHz (d)
f=4GHz (e)
f=5GHz (f)
f=10GHz


## 5. Conclusions

We apply the Kirchhoff migration approach to solve inverse scattering problems for time-dependent electromagnetic waves. The inversion scheme involves only integral calculations and is robust to polluted data at high noisy levels. We use both synthetic and experimental data to examine the performance of such a direct sampling scheme for locating point-like scatterers. The multiple-input-multiple-out scheme is used for imaging extended scatterers from the data generated by the software GprMax. Our experiments show that the data of one electric dipole can correctly reconstruct the location of two small scatterers, and that sixteen incoming sources can perfectly reconstruct the shape of an extended target. In addition, the source irradiated at 10GHz can be used to more clearly reconstruct the geometrical shape than the data of lower frequencies ranging from 1GHz to 5GHz. Our future work consists of imaging extended obstacles from experimental data buried in a random background medium, which is more challenging than the problem in a stationary homogeneous medium.



---

*Source: 6956809-2019-07-14.xml*