Electromagnetic Imaging of Two-Dimensional Geometries by Multipulse Interference Using the Inverse FDTD Method

The size, shape, and location of unknown objects in the ground and in the body can be estimated by an electromagnetic imaging technique. An imaging approach to clear detection of two-dimensional geometries is proposed in this paper. Based on the inverse finite-difference time-domain (FDTD) method, a phase interference technique using multidirectional pulses is employed. The advantage of the proposed method is that it can clearly reconstruct the geometry in a simple calculation. Sample imaging results are demonstrated. The analysis of the FDTD results shows that the detectable object size is limited by the incident wavelength and the measurement spacing and illustrates the detectability of multiple objects.


Introduction
Electromagnetic imaging is an inverse scattering technique to estimate unknown objects.The scattered wave from an object carries information about its geometry, size, and location.The original object shape can be reconstructed by numerically time reversing the scattering process.The inverse scattering technique has been widely studied for various applications [1][2][3][4][5] such as computed tomography, nondestructive evaluation, and geophysical remote sensing.
A variety of electromagnetic imaging techniques [6][7][8][9] have been proposed based on the finite-difference timedomain (FDTD) method.FDTD works in the time domain with arbitrary structures [10,11] and thus is suitable to calculate inverse scattering problems.
One of the problems of electromagnetic imaging is to reconstruct the detailed geometry of original objects.Previous works successfully detected the position and size of the two-and three-dimensional objects [5][6][7][8][9], but the reconstructed geometries were unclear even for an object immersed in a uniform medium [12][13][14][15][16][17][18].A comparatively successful imaging of the detailed geometry was demonstrated using phase interference with inverse FDTD simulations, in which the scattered field was calculated, and then using time reversal; the original shape was extracted from the interference patterns [19].However this requires both the scattered electric and magnetic fields in the entire space, and it is difficult to clearly discern the original boundaries due to the continuities of electromagnetic fields unless the original shape is used in the visualization.
In this paper, we propose a multipulse interference technique using the inverse FDTD method to improve the detection of geometrical details.In Section 2, we introduce the fundamentals of the proposed approach and demonstrate the imaging for some basic geometries.In Section 3, the sensitivity of our method to small and multiple objects is verified.In Section 4, the results are summarized.

Multipulse Interference Technique
In Section 2.1, the FDTD configuration for the electromagnetic imaging is defined.In Section 2.2, the imaging procedure of our multipulse interference technique is described.In Section 2.3, an example imaging for a square is demonstrated and our imaging quality is compared with the conventional interference technique.In Section 2.4, numerical results for some geometries, refractive indices, and polarizations are verified.

FDTD Configuration.
The conventional phase interference technique [19] must detect both the scattered electric and magnetic fields in the entire space (all grid points), but this is the impractical measurement actually.To more closely mimic what might be experimentally observed, we employ a circular array of measurement points in common use [5,[15][16][17][18]. Figure 1 shows our FDTD configuration.In order to illuminate the object over a wide range, a pair of opposing Gaussian pulses propagating along the ±-axis is employed (the plane wave is not suitable because it makes the temporal discontinuity when the detection is terminated).The center wavelength is   and the beam width of the pulse is 3  , where   is the arbitrary constant for scalability.The polarization is the transverse electric (TE) mode, in which the electromagnetic fields reduce to three nonzero components   ,   , and   .Scattered fields are observed at measurement points spaced at 5 ∘ intervals on the circumference of a circle of radius 6  .In the FDTD calculation, the grid spacing is   /20 according to the accuracy considerations given in [20].The computational domain of 15  × 15  is terminated by perfectly matched layer (PML) absorbing boundaries [11].
(1) The forward propagating scattered wave for the pulse pair is calculated using the FDTD method.The scattered electric fields are observed on the measurement points for 20 wave periods (enough time to detect primary scattered waves) after introducing the initial pulses.
(2) The object is removed from the computational domain in Figure 1, and the scattered fields sampled at the measurement points are time-reversely input as point sources (on the measurement points).Visualizing the intensity distribution, the inverse propagation is calculated using FDTD until the primary scattered pulse pair impinges on each other.The original shape emerges in the interference pattern of intensity, because the equiphase plane caused by reflection delay is parallel to the object surface.
(3) Iterating steps 1 and 2 with a pulse pair incident along a different direction, the reconstructed images are obtained.To get a clear image independent of pulse direction, the separate intensity distributions are summed.
The crucial difference between the proposed method and the conventional interference technique [19] is the use of the pulse pair.The incident pulse should illuminate the object surface to detect the reflected wave from the surface, and thus the single directional pulse in the conventional method is not enough to accurately reconstruct the shape.The multipulse interference improves the reproducibility of original boundaries dramatically as discussed in Section 2.3.Additionally, our approach without iterative processes is performed in a much simpler calculation than statistical optimization [21].

Example
Imaging for a Square.We illustrate the details of our approach by applying it to a square object.Figure 2 shows the imaging results for a square of side 3  .The refractive indices inside and outside the object are 3.5 (silicon in microwave region) and 1 (air), respectively.Figure 2(a) depicts the original shape (scale bar is 3  ).A pair of opposing pulses is input along the ±-axis and the scattered   fields are detected on the probe array.Using the sampled scattered fields, time-reversed propagation is carried out.Figure 2(b) shows the intensity distribution of the scattered   field due to the interference of the primary scattered pulses.The original shape is successfully reconstructed in the interference patterns, but the boundaries parallel to the axis are obscure compared with those parallel to the -axis.This is because the incident pulse pair propagating along the ±-axis mainly interacts with the boundaries parallel to the -axis.Hence using a pulse pair propagating along the ±axis clarifies the boundaries parallel to the -axis as shown in Figure 2(c).In order to obtain a clear image independent of pulse direction, we use pulse pairs propagating along the and -axes.However, the interference of electromagnetic waves propagating perpendicularly to each other gives rise to diagonal fringe patterns, and such intermittent boundaries are obscure.To circumvent this problem, we add the two separately calculated intensity distributions shown in Figures 2(b) and 2(c).Figure 2(d) shows the combined image (sum of the two intensity distributions).All boundary shapes  To compare our imaging quality with the conventional interference technique [19], in Figure 3, we plot the square reconstruction with the original shape line using (a) our method and (b) the conventional method.As shown in Figure 3, our result accurately captures the size and boundary, whereas the conventional method merely reconstructs the rough shape.Thus the crucial difference between the proposed method and the conventional approach is the reproducibility of the detailed geometry.

Numerical Verifications.
We verify the imaging qualities of our interference technique with respect to changes of the object shape, refractive index, and incident polarization.First we verify our approach for other basic geometries.Figure 4 shows (a) original shapes and (b) reconstructed images for triangle, cylinder, heart, z-like, and flower shapes.The refractive indexes inside and outside the object are 3.5 and 1, respectively.Successful images can be obtained not only for the straight-edged shape (triangle), but also for the curvilinear objects (cylinder and heart).Although the zlike shape is not axisymmetric, its geometry is successfully reconstructed (the small distortion of the shape is caused by the multiple scattering discussed in Section 3.2).In contrast, the dips in the flower image are comparatively obscure.This is because the incident pulses cannot effectively arrive in the dips smaller than   .Details of imaging such as small features are discussed in Section 3.1.We test our approach for other refractive indices.Figure 5 shows the imaging results for a square of refractive index (a) 2 (silicon nitride in microwave region) and (b) 1.5 (silicon dioxide in microwave region).The side of the square is 3  .The original geometry is correctly detected for both refractive indices.Because the equiphase plane in interference patterns is formed by the object shape, the geometry can be reconstructed unless the refractive index is too low to make a detectable reflection.
The proposed approach can be easily applied to the transverse magnetic (TM) mode, in which the electromagnetic fields reduce to   ,   , and   .The   fields propagating along the and -axes in the TE mode are replaced by the   and   fields, respectively, in the TM mode.In the observation step, both the scattered   and   fields are detected on the measurement points.In the imaging step, we visualize the   intensity after backwards propagation, calculated from the input   and   fields at the observation points.Figure 6 shows the scattered   intensity distribution for (a) square and (b) heart shape reconstructions in the TM mode (original shapes are given in Figures 2 and 4).Since both   and   fields are sampled on the probes, the reconstructions of the square (with boundaries parallel to the axes) and heart shape (with boundaries that are not parallel to the axes) are not inferior to the TE one (see Figure 4).

Sensitivity to Small and Multiple Objects
In this section, we examine the sensitivity of the proposed multipulse interference technique to small and multiple objects.We show that the minimum size of detectable objects is limited by the incident wavelength and the spacing of the measurement points.Adaptability for practical problems is tested by trying to detect multiple objects.

Minimum Size of Detectable Objects.
In principle, the minimum size of objects detected by electromagnetic imaging is restricted by the incident wavelength, because information about subwavelength features is carried by evanescent waves [22].But since we use a Gaussian pulse which has a wideband spectrum, features smaller than the center wavelength are detectable.In Figure 7, we show the imaging results for a cylinder with small bumps using different center wavelengths,   = 2 0 ,  0 , and  0 /2 ( 0 is a constant).The width of bumps is  0 /4.The refractive indices inside and outside the object are 3.5 and 1, respectively.Figure 7(a) depicts the original object (scale bar is 3 0 ).Figures 7(b), 7(c), and 7(d) show the results for   = 2 0 ,  0 , and  0 /2, respectively.When the short wavelength (shorter than the bump size) components are very small in the incident pulse, the bumps cannot be accurately reconstructed as shown in Figure 7(b).As the incident wavelength is decreased, as shown in Figures 7(c) and 7(d), the bumps become clearly detectable.
When the object is much larger than the incident wavelength, the minimum size of detectable objects is constrained by the spacing of the measurement probe array.As shown in Figure 8, a large object produces a scattered wave of similar width to its size without the strong diffraction.In order to capture the scattered wave passing through the gap between measurement points, the spacing must be no greater than the object size.We show the imaging results for a square using different probe spacing in Figure 9. Figure 9(a) depicts the original square of side  = 4  and of refractive index 3.5.Figures 9(b), 9(c), and 9(d) show the reconstructed images for the probe spacing of , 1.5 s, and 2 s, respectively.Using a spacing of , we successfully obtain the original shape as shown in Figure 9(b).Using a spacing of 1.5 s, the imaging quality is greatly degraded as shown in Figure 9(c).In Figure 9(d) with a spacing of 2 s, we can no longer discriminate the boundaries parallel to the -axis because a large part of scattered field propagating along -axis misses the probe array.

Multiple
Objects.We verify our approach for detecting multiple objects.Figure 10 shows the imaging for multiple objects.In Figure 10(a), we depict two cylinders of radius 1.5  and of refractive index 3.5.Figure 10(b) shows the imaging result for these two cylinders.We can discern the two cylinders, but the area between the objects is comparatively unclear.This is because the primary waves are multiply scattered.When the wave scattered by a cylinder propagates along -axis, it impinges on the other cylinder.The longer distance between objects has the smaller effect with each other in reconstruction, but basically in the inverse calculation such multiple scatterings cannot be captured without information about the original object distributions.
In general, however, the objects to be detected are buried with noisy structures which induce multiple scattering.In order to simulate such conditions, we configure a square surrounded by randomly arrayed particles in Figure 10(c).The side and refractive index of the square are 3  and 3.5, respectively.The radius and refractive index of the particles are   /10 and 1.5, respectively.The total area of 30 particles is 10% of the square.Figure 10(d) shows the imaging result for the square surrounded by noise particles.Although the centrally placed square is detectable, multiple scatterings by noise particles degrade the reconstruction compared with the result for the same square in free space (Figure 2).Especially in upper right corner of the square, the closely This result for noisy particles suggests that the multipulse interference technique is not suitable for nested structures.Figure 10(e) depicts a cylinder surrounded by a square (both of refractive indices are 3.5).The radius of the cylinder is   .The side and thickness of the square are 2  and   /2, respectively.In such a nested structure, only the outmost geometry is detectable because the scattered fields due to the inner structures are trapped by the outer structures.As shown in Figure 10(f), the inner structure cannot be detected, but the outmost square is clearly detected.

Conclusion
A multipulse interference technique for two-dimensional geometries has been proposed based on the inverse FDTD method.The proposed approach, using the sum of separately calculated intensity distributions, improves the imaging quality compared with the conventional method.The interference technique is based on the principle that the equiphase plane caused by reflection delay is parallel to the object surface and shows the robustness for changes of the object shape, refractive index, and incident polarization.The minimum size of detectable objects is similar in size to the incident wavelength due to the diffraction limit and is larger than the probe spacing needed to capture a scattered wave that is not greater than the object size.Thus, the shorter the wavelength and the closer the probe spacing, the better the imaging quality in our method.Multiple objects are detectable unless the multiple scattering is dominant.Noise particles that account up to about 10% of the object area within one object width of its surface are acceptable in our reconstruction.For nested structures, only the outmost geometry is detectable because the scattered fields due to the inner structures are trapped by the outer structures.
In future work we will try to extend the proposed method into the three-dimensional problems and the pseudo-spectral Advances in Optical Technologies time-domain method [23,24] that can improve the computational cost.

2 AdvancesFigure 1 :
Figure 1: FDTD configuration.A pair of Gaussian pulses of scalable center wavelength   propagates along ±-axis in the TE polarization and impinges on the object.Dots spaced at 5 ∘ intervals on the circumference of a circle of radius 6  are measurement points.Grid spacing and computational domain terminated by PML absorbing boundaries are depicted.

Figure 2 :Figure 3 :
Figure 2: Electromagnetic imaging for a square of side 3  .(a) Original shape (scale bar is 3  ).(b) Reconstructed intensity due to a pulse pair on -axis.(c) Reconstructed intensity due to a pulse pair on -axis.(d) Sum of intensity distributions of (b) and (c) yields an image reconstruction.

Figure 4 :Figure 5 :
Figure 4: Reconstruction of some basic geometries.(a) and (b) depict original shapes (scale bar is 3  ) and reconstructed images.From left to right, the triangle, cylinder, heart, z-like, and flower shapes, respectively, are depicted.

Figure 6 :
Figure 6: Reconstruction of (a) square and (b) heart shape in the TM mode (original shapes are given in Figures2 and 4).
Figure 6: Reconstruction of (a) square and (b) heart shape in the TM mode (original shapes are given in Figures2 and 4).

Figure 8 :
Figure 8: Scattering by an object larger than the incident wavelength.The object size and the measurement probe spacing are indicated.

Figure 9 :
Figure 9: Reconstruction of a square using different probe spacing.(a) Original square of side  = 4  (scale bar is 6  ).(b), (c), and (d) depict the reconstructed images for probe spacing of , 1.5 s, and 2 s, respectively.

Figure 10 :Figure 11 :
Figure 10: Reconstruction of multiple objects.(a) and (b) are original and reconstructed shapes of two cylinders, respectively.(c) and (d) are original and reconstructed shapes of a square surrounded by randomly arrayed particles, respectively.(e) and (f) are original and reconstructed shapes of a cylinder surrounded by a square.Scale bars are 3  .