Abstract

A fast and efficient microwave tomographic algorithm is proposed for 2-D and 3-D real-time intrawall imaging. The exploding reflection model is utilized to simplify the imaging formulation, and the half-space Green’s function is expanded in the spectral domain to facilitate the easy implementation of the imaging algorithm with the fast Fourier transform (FFT) and inverse fast Fourier transform (IFFT). The linearization of the inversion scheme and employment of FFT/IFFT in the imaging formula make the algorithm suitable for various applications pertaining to the inspection of a large probed region and allow real-time processing. Representative numerical and experimental results are presented to show the effectiveness and efficiency of the proposed algorithm for real-time intrawall characterization.

1. Introduction

Ground penetrating radar (GPR) provides unique noninvasive means for retrieving information on the internal construction of the building walls. The ability to determine fidelity information of the wall, such as hidden wires and tubes, cracks, and voids, and/or the existence of reinforcements is beneficial for both military and civilian applications, that is, building safety and durability assessment, bridge/dam defect characterization, and through-the-wall radar imaging (TWRI) [18].

The detection and imaging of concealed targets inside a visually opaque building wall can be addressed by means of the inverse scattering technique. In the past decade, various linear and nonlinear inverse scattering algorithms have been developed for nondestructive testing of abnormalities in concrete structures. In [9], a nonlinear inversion scheme based on the Gauss-Newton method is presented for the quantitative imaging of reinforcement bars in the concrete. The nonlinear inverse scattering algorithm is capable of quantitative reconstruction of the targets with super-resolution. However, it is based on an optimization scheme and thus is CPU intensive and very time-consuming. To accelerate the imaging speed, linear inverse scattering algorithms based on Born or Rytov approximations are proposed for radar imaging [1013]. Linear inversion schemes ignore multiple scattering effects and thus are only capable of reconstruction of weak scatterers. For strong scattering objects, it yields a qualitative reconstruction; that is, only the geometric properties, such as location and shape, are reconstructed [12, 13]. In [14], a qualitative imaging approach using the half-space Green’s function in combination with truncated singular value decomposition is proposed for reinforced concrete diagnostics. The multiple signal classification algorithm and linear sampling method are applied to nondestructive imaging of concrete structures in [1517]. Backprojection algorithms based on ray tracing technique taking into account the wave propagation velocity in the concrete have been widely used for intrawall nondestructive inspection in [18, 19].

Although successful imaging results can be achieved with the aforementioned algorithms, these algorithms mainly deal with 2-D intrawall imaging, which only reconstruct the cross sections of the concealed objects in the wall. In [20], 3-D imaging is applied to damage the assessment of the concrete structure using the delay-and-sum beamforming algorithm. A strategy for 3-D imaging through the fusion of 2-D images for concrete and pavement imaging is presented in [21]. 3-D imaging provides valuable information about the target’s extent in length, height, and width for enhanced target characterization and feature extraction. Most existing 3-D imaging techniques are computationally intensive and are usually processed off-line. For many applications, it is desirable that the intrawall imaging algorithms are computationally efficient such that the localization of the concealed objects inside the wall can be efficiently determined, preferably within a few seconds for 2-D imaging and a few minutes for 3-D imaging with a portable computer.

In this paper, a fast and efficient microwave tomographic algorithm is proposed for 2-D and 3-D real-time intrawall imaging. The imaging algorithm is based on the first-order Born approximation and exploits the half-space Green’s function. The exploding reflection model is employed, and then, the Green’s function is expanded in the spectrum domain to facilitate the implementation of the imaging algorithm with the fast Fourier transform (FFT) and inverse fast Fourier transform (IFFT). The linearization of the inversion scheme and employment of FFT/IFFT in the imaging formula make the imaging algorithm suitable for various applications pertaining to the inspection of large probed domain and allow near real-time processing.

2. Problem Formulation

2.1. Three-Dimensional Intrawall Imaging

Consider a typical scenario for 3-D intrawall inspection using monostatic synthetic aperture radar (SAR), as shown in Figure 1. The transmitter transmits an ultra-wideband (UWB) electromagnetic signal and collects the reflected signal over a rectangular planar aperture that may have one or more targets within its field of view. The reflected signal is recorded coherently by the receiver, digitized, and stored in the computer. Then, the region between the front and back walls can be efficiently reconstructed to form a focused image of the targets inside the wall using a subsurface imaging algorithm.

As shown in Figure 1, the radar operates in region I, which is the free-air in front of the wall. The targets of interest are embedded in region II, which is the building wall whose dielectric constant and thickness are denoted as εb and d, respectively. As opposed to TWRI, which aims at imaging of targets behind the wall (region III) [58], intrawall imaging is interested in the imaging of targets inside the wall (region II). The radar is located at and operates over the frequency range from fmin to fmax. We assume the investigated domain is a cuboid region denoted as Dinv = [−xmin, xmax] × [−ymin, ymax] × [−zmin, zmax]. The complex permittivity inside Dinv is defined as where is the position vector of the target.

The scattered electric field observed at the receiver location can be expressed as where k is the wavenumber in the free space, Et (r, rR, k) is the total electrical field inside the target, and O(r, k) is the contrast function defined as where is the relative complex permittivity profile of the target and G (rR, r, k) is the dyadic half-space Green’s function which can be expressed in spectrum domain as [12, 22] where

As the total electrical field inside the target is also a function of the contrast function, (2) is a nonlinear equation. To linearize the inversion scheme, the first-order Born approximation can be used, which assumes that the total field inside the target is approximated by the incident field [23]. Then, the mathematical relationship between the scattered field and the contrast function can be written as where the incident field is given by wherein is the antenna polarization direction.

Substituting (7) into (6) gives where is the free space wave impedance and the reciprocity relationship has been used in the above equation.

For simplicity, (8) can be expressed in the following operator form: where the operator L maps the contrast function from the object space to the data space,

Then the image can be formulated as where La is the adjoint operator of L and maps the scattered field from the data space to the target space [24, 25], where the subscript is the complex conjugation. From the concept of time reversal imaging, the complex conjugation in the frequency domain is equivalent to the time reversal in the time domain. The back propagation of the time-reversed field will focus on the target.

It is found from (12) that an efficient evaluation of the subsurface Green’s function is critical to the imaging. However, the exact calculation of the dyadic Green’s function requires the evaluation of the Sommerfeld integral, which is generally complicated and computationally expensive. To simplify and accelerate the imaging speed, the assumption of exploding reflection model is used [24, 26], where the squared term, representing the two-way propagation, is replaced by a one-way propagation with half the propagation speed. Then (12) can be further written as

It is noticed that the radar system both transmits and receives in the y direction, . By substituting the 3-D spectrum form of the half-space Green’s function into (14), the image can be reconstructed as where and where k and kb are the wavenumbers in the free-space and, in the wall, is the 2-D FFT of with respect to and ,

The 3-D intrawall imaging formulation in (15) can now be efficiently implemented in the following steps: (1) perform 2-D FFT to compute the spatial Fourier transform of the scattered field ; (2) multiply with all the functions inside the inner integral except the exponential term ; and (3) at each down range pixel z applies 2-D IFFT to evaluate the inner integral and then sums over all the frequencies to calculate the outer integral.

2.2. Two-Dimensional Intrawall Imaging

The above derivation of a 3-D imaging algorithm can be simply reduced to 2-D intrawall imaging by replacing the half-space dyadic Green’s function with 2-D half-space Green’s function [22]

Substitute the 2-D half-space Green’s function into (14); the 2-D imaging formula can be derived as where , , , and are the FFT of with respect to :

3. Numerical Results

Representative numerical examples are presented in this section to verify and show the efficiency of the proposed algorithms for 2-D and 3-D intrawall imaging. In all the following examples, the imaging is performed on a four-core P4 2.66 GHz desktop computer with 16 GB of RAM.

In the first example, we present the 2-D imaging of a reinforced wall with a sinusoidal distribution of steel bars inside the wall as shown in Figure 2(a). The diameter of each steel bar is 1 cm with a spacing of 20 cm in the horizontal direction. The distribution of the steel bars is an approximate sinusoid function. The received reflected signal was simulated using a full wave electromagnetic solver based on the finite difference time domain (FDTD) method. The monostatic radar measures the return signal at a standoff distance of 0.3 m over a synthetic aperture extending from −2 m to 2 m with a step of 0.025 m. The dielectric constant, conductivity, and thickness of the wall are εb = 6, σ = 0.01 S/m, and d = 0.2 m, respectively. The operating frequency is 0.8 GHz to 3.2 GHz with a step Δf = 22.2 MHz. Figure 2(b) is the time domain B-scan raw data of the simulated received signal. Frequency domain data is obtained, by performing an FFT on the B-scan time domain signal, and is used in the imaging.

In the following simulations, unless otherwise specified, we assume that the wall dielectric constant is known and is used in the imaging. The estimation of wall parameters associated with single or multilayered building walls has been extensively studied in the past decade within the framework of a one-dimensional inverse scattering problem [2734]. For a single-layer wall, radar-based closed-form estimation of the wall dielectric constant, conductivity, and thickness is given in [30, 31].

Figure 2(c) shows the imaging result of the steel bars inside the wall using the proposed algorithm with the known wall dielectric constant εb = 6, in which the back wall is highlighted with a dashed black line. From this figure, we find that the 20 sinusoid distributed steel bars and the back boundary of the wall are clearly identified and correctly localized at their true locations. For comparison, the imaging results using the ray tracking-based back-projection algorithm from [19] are provided in Figure 2(d). It takes about 57.34 s to reconstruct the image using the ray tracing-based back-projection algorithm. The proposed imaging algorithm, however, requires only 0.64 s to reconstruct the image in Figure 2(c) on the same computer. Due to the utilization of FFT/IFFT in the proposed imaging algorithm, the imaging process can be achieved in a short computation time, which is desirable for on-site applications. From the comparison of Figures 2(c) and 2(d), we also find less cluttered and a cleaner image can be achieved using the proposed algorithm due to its more accurate modeling of the wave propagation phenomenon through incorporation of the half-space Green’s function.

Incorrect estimation of the wall parameters could result in the shift of target position and defocussing of the target image. Figures 2(e) and 2(f) show the imaging result with an underestimated wall parameter (εb = 3) and an overestimated wall parameter (εb = 10). In the two figures, the back boundary of the wall is highlighted with a dashed black line. It is observed from the two figures that underestimation of the wall dielectric constant causes a shift of the target position toward the radar in the downrange. Overestimation of the wall dielectric constant causes a shift of targets position away from the radar system in downrange.

In the above simulation, the radar scans the wall from −2 m to 2 m at a step size of 2.5 cm. To investigate the effect of step size on the imaging, we present imaging results with 2.5 cm, 5 cm, and 10 cm step sizes in Figures 3(a)3(c), respectively. From the three figures, we find degradation in the azimuth resolution by increasing the step size of the synthesized aperture. The minimum step size should be no less than the Nyquist sampling interval , where λ is the wavelength of the transmitted signal and εb is the dielectric constant of the wall.

The downrange resolution is primarily determined by the bandwidth of the transmitted signal. Figures 4(a)4(c) are the imaging results with 2.4 GHz, 1.2 GHz, and 0.6 GHz bandwidth, respectively. As expected, the downrange resolution degrades as the bandwidth is reduced. At the narrowest bandwidth of 0.6 GHz, we can only resolve the front and back boundaries of the wall but cannot resolve the steel bars inside the wall, as shown in Figure 4(c).

To quantitatively characterize the results, the signal fidelity factor in [3537] is calculated to quantify the similarity of the received and transmitted signals. In order to compare only the shape of the signals but not their magnitudes, the received and the transmitted signals are normalized according to

The signal fidelity factor (SFF) is then defined as the peak of the cross-correlation of the normalized transmitted and received signals and is given as [37]

We note that a 180-degree phase change occurs as the signal is reflected from the wall; thus, the minimum of the correlation is used and a minus sign is included in (22). Figure 5 shows the SFF of the reinforcement wall, shown in Figure 2(a), at different antenna locations. For a perfect metal wall, the SFF is expected to be 1 as no distortion occurs in the shape of the transmitted and received signal. Due to the scattering from the steel bars and the multiple scattering effects inside the wall, signal distortion occurs as it is reflected from a reinforced wall.

In the second example, we present the 3-D imaging of a reinforcement grid inside the wall. Figure 6(a) shows the dimensions of the reinforcement grid and the interelement spacing between the reinforcement bars. The radius of each steel bar is 1 cm. The spacing between the fifth and sixth reinforcement bars is 0.3 m, and the spacing between the sixth and seventh reinforcement bars is 0.1 m in the vertical direction. The interelement spacing between other steel bars is 0.2 m. The permittivity, conductivity, and thickness of the wall are εb = 6, σ = 0.01 S/m, and d = 0.3 m, respectively. The simulation was performed over a 2 m × 2 m square aperture parallel to the wall at a standoff distance of 0.1 m. The operating frequency covers a range from 0.8 GHz to 2.5 GHz with a step of 26 MHz. The imaging region is 2 m × 2 m × 0.6 m and is discretized into 160 × 160 × 80 pixels. Figure 6(b) shows the 3-D imaging result of the reinforcement grid inside the wall. From this figure, it is clear that the proposed algorithm can successfully identify and correctly localize not only the back boundary of the wall but also all the steel bars inside the wall. The large spacing between the fifth and the sixth steel bars and the small spacing between the sixth and seventh steel bars can be clearly observed in the image. The 3-D imaging of the reinforcement grid in Figure 6(b) took only 28.6 s using the proposed imaging algorithm.

In the last example, we present the 3-D imaging of two perpendicular tubes hidden inside the wall. The dimension and orientations of the tubes are shown in Figure 7(a). The lengths of the long and short tubes are 0.8 m and 0.2 m, respectively. Figure 7(b) shows the 3-D imaging result of the two tubes inside wall. We find that the long tube and short tube can be clearly observed in the 3-D image and the orientations of the tubes are also clearly identified. The successful retrieval of the width, length, and orientation of the concealed objects inside the wall provides valuable information for building maintenance and building quality and durability assessment as well as in detection of weapons and other related objects in various homeland security and law-enforcement scenarios. It took only 28.3 s to reconstruct the 3-D image. The fast imaging speed shows the promise of the proposed 3-D intrawall imaging for on-site applications.

4. Experimental Results

Some experimental examples are presented in this section to verify the proposed algorithm for 2-D and 3-D intrawall imaging in a controlled lab environment. The 2-D imaging algorithm was applied to three experimental cases of a solid concrete wall with (1) embedded rebar, (2) the wall with no backing, and (3) a mirror placed on the back face of the wall. The last test of the mirror behind the concrete wall was also processed using the 3-D imaging algorithm.

The 2-D intrawall imaging algorithm was applied to the case of metallic rebar embedded within a solid concrete wall. Due to the limitations of fabricating this type of scenario in full scale, a 1/3 scale model was created which utilized multiple layers of a cement board readily available in home improvement stores. Multiple layers of the board were required as it was sourced at 1/4 thickness, and thus to make a 1/3 scale replica of a full thickness wall, multiple pieces were layered to create a single wall. This layering process allowed for the inclusion of metallic strips which were analogous to metallic rebar placed within a full size solid concrete wall. Additionally, because of the layering, the metallic strips were composed of copper tape in 3-layer thickness. Thicker strips could not be used while still maintaining a continuous contact between cement board layers. Due to the 1/3 scale of the wall, the frequency range for the test was increased by a factor of 3 to cover 3 GHz to 9 GHz which corresponds to a standard TWRI frequency range of 1 GHz to 3 GHz without scaling. An Agilent PNA with an associated horn antenna which covered 3–9 GHz was used to measure the wall analogue. The measurement setup along with the diagram of the wall with metallic strip locations is shown in Figure 8. The resulting image using the proposed 2-D imaging algorithm is shown in Figure 8 where the metallic inclusions are well localized within the center of the wall. Additionally, striations are visible in the image corresponding to the interfaces between the various layers of cement board.

In the next set of experiments, radar measurements were performed for a concrete wall, with and without backing. The wall under test was constructed using solid concrete bricks with a known dielectric constant of 7.66 and a loss tangent of 0.158 S/m and downrange thickness of 6 inches. The measurement array was fixed at 3 m standoff distance, measured from the front of the wall to the feed point on the antenna. The backing, which was a mirror, was placed at the rear surface of the wall in a configuration which allowed for the most direct contact between the backing and the rear surface of the wall. A frequency range of 0.7 GHz to 3.1 GHz was chosen for imaging using a step size of 3 MHz. An Agilent ENA 5071B was used for signal transmission and data collection. A dual-polarized horn antenna, with an operational bandwidth from 0.7 GHz to 6 GHz, was used as the transceiver. Figure 9 shows the experimental setup of the solid concrete wall with no backing, a selected reflection signal from the wall, and its associated image resulting from the proposed 2-D microwave tomographic technique. Figure 10 shows the results of imaging the same solid concrete wall of Figure 9 with the addition of a vertically oriented mirror placed directly behind the wall and in direct contact with the wall. In Figure 10(c), the imaging result clearly shows the reflections from the front of the wall, at 0.05 m, and the mirror on the back of the wall, at 0.2 m. The reflection from the rear face of the wall is low compared to that of the mirror but is still visible at 0.2 m downrange to the left and right of the mirror reflection. The third ringing of the wall and mirror can also be seen in this image at 0.37 m. In all cases, the computation time was less than 1 second for image generation. In Figure 11, the comparison of the SFF calculations for the solid concrete wall using simulated and measured results is presented. A relatively close agreement between the SFF values of simulated and measured received signals versus the antenna positions is obtained.

The 3-D imaging algorithm, discussed above, was also applied to the measured scenario of a mirror on the rear surface of the wall in both vertical and horizontal configurations. Measured data was collected in a grid with 97 positions in the x-direction and 23 positions in the y-direction and at a 3 m standoff distance from the wall. At each location, step frequency data was collected from 1 GHz to 6 GHz with 1601 frequency steps. Both measured data sets were then processed using the 3-D microwave tomography algorithm given in (15). Figures 12(a) and 12(b) show a side by side comparison of the vertical and horizontal mirror imaging results. In both cases, the mirror orientation is readily identifiable. We note that the above and other cases of detection and imaging of the wall’s conductor backings are of interest in the planning of various search and rescue mission scenarios. Additionally, it is of note that the experimental results shown in Figures 811 are individual experiments which are representative examples from the multiple data collections performed during this research. These results are not an average over the many experiments performed but rather are selected examples which are indicative of the results.

5. Conclusion

The fast and nondestructive inspection of wall interior structures has promising applications in building safety and durability assessment as well as in various defense and law-enforcement scenarios. In this paper, a fast and efficient microwave tomographic algorithm is proposed for near real-time 2-D and 3-D intrawall imaging. The linearization of the inversion scheme and employment of the FFT/IFFT in the imaging formulation make the proposed imaging algorithm suitable for various applications pertaining to the inspection of a large probed region and allow real-time processing. Numerical and experimental 2-D and 3-D examples are presented with each showing the successful retrieval and localization of the wall interior or wall backing information using the proposed imaging algorithm. It takes less than one second/minute for 2-D/3-D imaging using the proposed algorithm, a feature that is very attractive for many applications which require on-site processing. Although this paper is concentrated on the intra-wall inspection, the presented algorithm is applicable to other concrete structure nondestructive testing, such as bridge/dam defect characterization, pavement inspection, and archeology. We note that the proposed radar imaging algorithm utilizes a linearization scheme to significantly speed up the imaging of intrawall objects and results in the object’s accurate geolocation and shape reconstruction but not its material composition, such as permittivity and loss tangent. The later requires nonlinear inverse-profiling schemes which usually necessitate a numerically computed forward model and a global optimization inversion scheme thus making the inversion very time-consuming and not suitable for on-site processing.

Conflicts of Interest

The authors declare that they have no conflicts of interest.