We propose a three-dimensional microwave holographic imaging method based on the forward-scattered waves only. In the proposed method, one transmitter and multiple receivers perform together a two-dimensional scan on two planar apertures on opposite sides of the inspected domain. The ability to achieve three-dimensional imaging without back-scattered waves enables the imaging of high-loss objects, for example, tissues, where the back-scattered waves may not be available due to low signal-to-noise ratio or nonreciprocal measurement setup. The simulation and experimental results demonstrate the satisfactory performance of the proposed method in providing three-dimensional images. Resolution limits are derived and confirmed with simulation examples.
1. Introduction
Better penetration of microwaves inside lossy dielectric objects, compared to visible light, provides a means for inspection of biological tissues, wood, concrete, ceramics, soil, and so forth. Various techniques have been proposed to harness the ability of microwaves for two-dimensional (2D) and three-dimensional (3D) imaging in a wide range of applications such as biomedical imaging [1], concealed weapon detection [2], through-the-wall imaging [3], and nondestructive testing and evaluation [4]. All these methods operate on the scattered signals at the antenna terminals but they differ significantly in the data acquisition schemes and the processing algorithms. A summary of major microwave imaging techniques is available in [5].
In this paper, the focus is on microwave holographic imaging. Microwave holography was originally developed by Farhat et al. [2, 6, 7] in a form similar to acoustic holography [8, 9], which can be viewed as a long-wavelength implementation of the original optical holography developed by Gabor [10]. These microwave/acoustic holography techniques operate by recording a “hologram” which contains information about the magnitude and the phase of the wave reflected by a target. The hologram is captured on a film by interfering the scattered wave due to the examined object with a reference wave. This interference pattern or “hologram” can then be used to reconstruct an image of the object by illuminating the hologram with the same reference wave. In modern microwave/acoustic holography, scanned transceivers capture the scattered wavefront’s phase and amplitude distributions directly. Once recorded, the image can be reconstructed digitally using a Fourier-optics image reconstruction that essentially emulates the hologram-based reconstruction process.
In [11], a 3D microwave holographic imaging technique has been proposed. The authors described the technique as the merging of the single-frequency 2D holographic imaging methods with wideband 2D synthetic aperture radar (SAR). The SAR imaging systems are commonly used to obtain high-resolution 2D images of a terrain at long ranges. These systems operate by scanning a wide bandwidth radar transceiver along a long linear aperture, typically using an airplane. Then they reconstruct a 2D image along the scanned line (e.g., see [12–15]). In [11], wideband back-scattered data is collected over a rectangular aperture. The system is quasi-monostatic that is, the transmitting and the receiving antennas are separate but colocated during the 2D scan. The processing relies on an assumed analytical (exponential) form of the incident field and Green’s function in order to cast the inversion in the form of a 3D inverse Fourier transform (FT). This limits the technique to far-zone imaging. Resampling of the data in kz-space is also necessary. Also, in [16–18], similar 3D holographic imaging techniques have been proposed where the wideband data is acquired on a cylindrical aperture to extend the “view angles” for the object. This leads to better imaging results.
In [19, 20], we extended the 3D holography technique in [11] to near-field imaging. The proposed method allows for incorporating forward-scattered signals in addition to the back-scattered signals. The method also allows for incident-field and Green’s function distributions in numeric forms. These are obtained through simulations of the particular background medium with the particular Tx/Rx antennas (Tx and Rx stand for transmitter and receiver, resp.). The accurate representations of the incident field and Green’s function for the particular acquisition setup and antennas (not assuming a plane-wave exponential form as in [11]) are crucial in near-field imaging. The procedure in [11] is inapplicable with numeric representations of the incident field and the Green’s function. Thus, the numerical form of the incident field and the Green’s function required a new inversion procedure. In [19, 20], a linear system of equations is solved for each spatial frequency pair (kx, ky) and 2D inverse FT is applied to the solution on planes (slices) at all desired range locations. These linear systems of equations have much smaller dimensions and they are less ill-conditioned compared to the systems of equations in regular optimization-based microwave imaging techniques.
Also, the algorithm is not limited to propagating waves (where kz is assumed real as in [11]) and is thus capable of processing the evanescent waves available in near-field imaging. Processing larger values for (kx, ky) leads to better spatial resolution of the reconstructed images. Furthermore, in the approach proposed in [11], resampling of the data in kz-space is necessary which may lead to additional errors. Algorithms such as the Stolt interpolation [21] and nonuniform fast Fourier transform (NUFFT) [22] have been proposed to alleviate or avoid this problem. Even if the resampling (interpolation) errors can be alleviated, assuming kx, ky, and kz as independent variables still leads to errors in the image reconstruction process. These issues are irrelevant in our technique since the use of the parameter kz is avoided.
In [20], we presented experimental holographic imaging results when employing the TEM horn antennas proposed in [23]. Thus, the main focus was on applications for breast-cancer diagnostics and artificial glycerin-based phantoms emulating the breast tissues were imaged. Only 2D holographic imaging results were presented employing forward-scattered waves. This is due to the fact that the back-scattered (or reflected) signals were measured through the complex reflection coefficient of the Tx antenna. Such signals are so weak that they are masked by the measurement noise and uncertainties. On the other hand, the back-scattered signals are crucial for implementing 3D holographic imaging with planar scanning.
In this paper, we propose a configuration which allows for 3D holographic imaging where only forward-scattered data is available. We develop expressions for the resolution limits of bistatic holography using an approach developed for bistatic SAR [24]. We confirm that the availability of forward-scattered data in a planar surface acquisitions similar to those in [19, 20] cannot provide enough information to perform 3D holographic imaging. This is consistent with our previous observations. The images obtained when we exclude the back-scattered data do not have depth resolution. Moreover, the image quality in the cross range is also compromised when attempting 3D reconstruction without back-scattered signals, although a 2D holographic reconstruction in the cross-range plane of the object is successful. Overall, without the back-scattered data, the holographic reconstructed images contain strong artifacts along both range and cross range. The solution to this problem as proposed in this work is a scanning setup comprising one transmitter and five receivers which move together during the 2D scan on two rectangular apertures on opposite sides of the inspected region. We show that 3D image reconstruction is possible with this setup when employing only forward-scattered data. We also derive the range and cross-range resolution limits.
2. Resolution in Two-Antenna Configuration
We consider a two-antenna microwave holography setup where the antennas scan simultaneously two rectangular planar surfaces in a raster pattern. This setup is illustrated in Figure 1 where antenna 1 and antenna 2 perform the scan together on aperture 1 at z=0 and aperture 2 at z=z-. At each sampling step, the wideband transmission coefficient of this two-port system is acquired, which represents the forward-scattered wave. We ignore the back-scattered waves since in many applications in lossy media, for example, tissue imaging, these waves are too weak and are likely to be masked by noise and uncertainties in the measurements. Also, the measurement system may be nonreciprocal due to the use of amplifiers (or other nonreciprocal components) at the transmitting and/or receiving sides. This makes the measurement of the reflected signals impossible unless separate Tx and Rx channels are provided.
Microwave holography setup. The grey thick arrows show the directions of 2D scan.
In [19, 20], both antennas are assumed to be always in the same (x′,y′) position during the scan. Here, we consider the general case where antenna 2 may have an offset of ±Dx and ±Dy with respect to antenna 1 along the x- and y-axes, respectively.
In the following, we study the spatial resolution limits of the two-antenna imaging setup employing an approach previously proposed for bistatic SAR [24]. Note that the resolution limits derived with this method (which assumes propagating waves) provide good estimates in the case of the two-antenna holographic setup where the object is in the far-field region of the antennas. However, in the near-field imaging it is possible to obtain better resolution due to the availability of evanescent waves. The closer the antennas are to the object, the better the resolution is in a near-field imaging scenario.
With reference to Figure 1, at each sampling position (x′,y′,0) for antenna 1—this implies that antenna 2 is at (x′+Dx, y′+Dy,z-)—the total path traveled by the wave through the object is denoted by R and is written as
(1)R=R1+R2,
where
(2)R1=(x-x′)2+(y-y′)2+z2,R2=[(x′+Dx)-x]2+[(y′+Dy)-y]2+(z--z)2.
The differential change in R, denoted by ΔR, due to an incremental change (Δρss^) in the position of the object in any arbitrary direction s^, is written as [24]
(3)ΔR=(s^T·∇R)Δρs,
where ∇ is the gradient operator and the superscript T represents the transpose operator. Here, s^ is a unit vector.
On the other hand, the minimum “measureable” ΔR is related to the velocity of the wave in the medium c and the bandwidth of the imaging system B as [24]
(4)ΔRmin≈cB.
From (3) and (4), the resolution limit in the direction s^, denoted by ρ-s, is estimated as [24]
(5)ρ-s=cB·1|s^T∇R|.
Please note that the resolution in this paper is defined as the half-power (or half-intensity) width of a point-like object.
To estimate the resolution of the imaging system along the x, y, and z directions, first the gradient of R in a rectangular coordinate system is written as
(6)∇R=[∂R∂x∂R∂y∂R∂z]=[(x-x′)R1-[(x′+Dx)-x]R2(y-y′)R1-[(y′+Dy)-y]R2zR1-(z--z)R2].
Then, from (5) and (6), the cross-range resolution ρ-u(u=xory) and the range resolution ρ-z are obtained as
(7)ρ-u=cB·1|((u-u′)/R1)-([(u′+Du)-u]/R2)|,u=x,y,u′=x′,y′,ρ-z=cB·1|(z/R1)-((z--z)/R2)|.
In the previous equations, ρ-u and ρ-z are obtained when s^ in (5) is substituted with the unit vector in the particular direction, x, y, or z.
We assume that the object is very small and it is positioned on the z-axis that is, the object is at (0,0,z). Then, (7) are simplified as
(8)ρ-u=cB·1|(-u′R1′)-(u′+DuR2′)|,u′=x′,y′,Du=Dx,Dy,ρ-z=cB·1|(z/R1′)-((z--z)/R2′)|,
where
(9)R1′=x′2+y′2+z2,R2′=(x′+Dx)2+(y′+Dy)2+(z--z)2.
As observed from the previous equations, the range and cross-range resolution limits depend on the positions of the transceivers and the object.
To derive the resolution limits, we consider a 2D cut in the 3D setup in Figure 1 as shown in Figure 2, where the axis u represents either x or y. Du and u′ denote the offsets of the receiver with respect to the transmitter and the position of the transmitter, respectively.
Simplifying the 3D setup in Figure 1 to a 2D setup for deriving resolution limits in terms of θ1 and θ2. When Du=0, the range resolution is poor in the shaded region. The grey thick arrows show the directions of the scan.
With reference to Figure 2, we first rewrite the resolution limits in (8) in terms of the angles θ1 and θ2 as
(10)ρ-u=cB·1|sinθ1+sinθ2|,(11)ρ-z=cB·1|cosθ1-cosθ2|.
These are the angles of incidence and scattering, respectively, associated with the imaged point-like object.
From (11), it follows that when the transmitter and the receiver are at the same position (x′,y′) during the 2D scan (as in [19, 20]), that is, when Du=0, there are two main regions in the inspected volume where the range resolution ρ-z is poor. These two regions are (1) the middle of the inspected domain, that is, z≈z-/2 and (2) the direct path connecting the transceivers when it passes right through the scatterer; that is, θ1≈θ2≈0. In Figure 2, the regions of poor range resolution are shaded in gray.
On the other hand, when the receiver has an offset with respect to the transmitter, that is, when Du≠0, we obtain finite values for both the cross-range and the range resolution limits everywhere in the inspected volume. In this configuration, as per (10), the best cross-range resolution is obtained when the antennas are at the edges of the apertures; that is, the angles θ1 and θ2 are close to π/2:
(12)ρ-umin|θ1=θ2=±(π/2)=c2B.
Regarding the cross-range resolution, as per (10), it is best when one antenna is on the z-axis (θ1=0) while the other antenna is at u→∞ (Du→∞, θ2→π/2):
(13)ρ-u=cB.
Under these conditions, the range resolution is
(14)ρ-zmin|θ1=0,θ2=±π/2θ1=±π/2,θ2=0=cB.
It is worth noting that in reflection holography where θ2=π-θ1, the cross-range and range resolution limits are
(15)ρ-u=cB·1|2sinθ1|,ρ-z=cB·1|2cosθ1|.
Assuming B≈2fc with fc being the center frequency of the band, the cross-range resolution can be written as follows:
(16)ρ-u=λc4sinθ1,
where λc is the wavelength at fc. Thus, in this case, the cross-range resolution improves if the size of the aperture is large (θ1→π/2). On the other hand, the best range resolution limit is obtained when θ1→0 as
(17)ρ-zmin=c2B.
These results are consistent with the expressions for the resolution limits of reflection holography in [11].
3. Multiple Receiver Setup
In the previous section, we showed that having a nonzero offset distance for the receiver along the x- or y-axis (Du≠0, u=x,y) leads to improving the range and cross-range resolutions in a scenario where forward-scattered signals are acquired. Thus, here we propose a multiple receiver setup in a star distribution to achieve satisfactory range and cross-range resolutions in 3D microwave imaging with planar raster scanning.
Figure 3 illustrates the setup where one transmitter illuminates the object while five receivers measure the forward-scattered waves. Antenna 2 is aligned with the transmitter (antenna 1) while the other receivers (antennas 3 to 6) have offset distances of ±Dx and ±Dy along the ±x and ±y directions, respectively. The transmitter and the five receivers move together during the 2D scan on the two planar apertures. From the results in Section 2, it follows that larger offset distances lead to improved resolution. However, at large offset distances the scattered wave travels along longer paths from the object to the receiver. This weakens the signal due to two factors: spatial spread and attenuation if the medium is lossy. Besides, signal strength can also weaken due to the antenna pattern. These factors impose upper limits on the offset distances. In addition, increasing the offset distances increases the size of aperture 2. If the size of aperture 1 is 2Lx×2Ly, then the size of aperture 2 is 2(Lx+Dx)×2(Ly+Dy).
Proposed 3D microwave holography setup using forward-scattered waves only. Antenna 1 illuminates the object while antennas 2 to 6 receive the scattered waves. The grey thick arrows show the directions of the scan.
4. 3D Holographic Image Reconstruction
In this section, we extend our previously proposed 3D holographic microwave imaging algorithm [19] to process the data obtained with the multiple receiver configuration proposed here.
With reference to Figure 3, the transmitting antenna (antenna 1) and the receiving antenna array (antennas 2 to 6) perform a 2D scan at the planes z=0 and z=z-, respectively. For simplicity, we assume that the antennas are x-polarized (e.g., dipoles oriented along the x-axis). Thus, the field can be reasonably approximated by a TMx polarization and the theory is scalar in nature (the algorithm for a full-vector 3D holography is available in [20]). The approach here is directly applicable to acquisition of co- and cross-polarized data. In the scalar case, Green’s function G can be viewed as the Gxx element of the full dyadic while the E-field is represented by its x-component only. From now on we omit the subscript x for brevity.
Assume that at any measurement frequency ωl(l=1,2,…,Nω) we know the incident field Einc(0,0,0;x,y,z;ωl) at any point (x,y,z) in the inspected volume due to antenna 1 when it is at (0,0,0). In addition, Green’s function for antenna 2 G2(x,y,z;0,0,z-;ωl) is known for an x-polarized scattering point source at (x,y,z) and the x-polarized response at (0,0,z-). This information can be obtained via simulations. For brevity, we introduce the notations
(18)Einc(x,y,z,ωl)≡Einc(0,0,0;x,y,z;ωl),G2(x,y,z,ωl)≡G2(x,y,z;0,0,z-;ωl).
In [20], we show that while Einc can be obtained from simulations, Green’s function can be obtained from Einc using the reciprocity principle. Thus, assuming that antennas 1 to 6 are identical, only one simulation suffices to obtain Einc(x,y,z,ωl), from which Green’s functions for antennas 2 to 6 are obtained as
(19)G2(x,y,z,ωl)=Einc(x,y,z¯-z,ωl),G3(x,y,z,ωl)=Einc(x,y-Dy,z¯-z,ωl),G4(x,y,z,ωl)=Einc(x,y+Dy,z¯-z,ωl),G5(x,y,z,ωl)=Einc(x-Dx,y,z¯-z,ωl),G6(x,y,z,ωl)=Einc(x+Dx,y,z¯-z,ωl).
Let Eksc(x′,y′,z¯,ωl) be the scattered E-field received by the kth antenna (k=2,…,6) when the transmitting antenna is at (x′,y′,0). Following the approach presented in [19], Eksc(x′,y′,z¯,ωl) is written as
(20)Eksc(x′,y′,z¯,ωl)=∫z∫y∫xf(x,y,z)·ak(x′-x,y′-y,z,ωl)dxdydz,fork=2,…,6,
where f(x,y,z)=ks2(x,y,z)-kb2 is the contrast function, ks and kb are the wavenumbers of the object and the background mediums, respectively, and
(21)ak(x,y,z,ωl)=Einc(-x,-y,z,ωl)Gk(-x,-y,z,ωl).
Notice that in (20), the integration over x and y can be interpreted as a 2D convolution integral. Thus, the 2D FT of Eksc is written as
(22)E~ksc(kx,ky,z¯,ωl)=∫zF~(kx,ky,z)A~k(kx,ky,z,ωl)dzfork=2,…,6,
where F~ and A~k are the 2D FTs of f(x,y,z) and ak(x,y,z,ωl), respectively. By approximating the integral in (22) by a discrete sum and employing the data collected at all frequencies and at all receivers, we construct a system of equations at each (kx,ky) as
(23)E~_=A~__F~_,
where
(24)E~_=[[E~_2sc]⋮[E~_6sc]],A~__=[[A~__2]⋮[A~__6]],F~_=[F~(z1)⋮F~(zNz)],
and, for each k=2,…,6,
(25)[E~_ksc]=[E~ksc(z-,ω1)⋮E~ksc(z-,ωNω)],[A~__k]=[A~k(z1,ω1)Δz⋯A~k(zNz,ω1)Δz⋮⋮⋮A~k(z1,ωNω)Δz⋯A~k(zNz,ωNω)Δz].
Here, Δz is the discretization step along the z-axis. In (24) and (25), the arguments kx and ky have been omitted for brevity.
The system in (23) is solved at each (kx,ky) pair for F~(zn), n=1,…,Nz, in a least-square sense. To reduce the ill-conditioning of the least-square solution we employ the Tikhonov regularization [25]. Thus, the solution is obtained from
(26)F~_=(A~__*A~__+Γ__*Γ__)-1A~__*E~_,
where (·)* is the complex conjugate operator, Γ__=αΙ is the Tikhonov matrix, α is the Tikhonov regularization parameter, and I is the identity matrix. The optimal value of the regularization parameter is usually determined by an adhoc method in practical problems. In this study, we use the plot of the Lagrange function to determine α. The MATLAB code for this method is available in MATLAB central [26].
Once the systems of equations for all (kx,ky) are solved, the inverse 2D FT is applied to F~(kx,ky,zn), n=1,2,…,Nz, to reconstruct a 2D slice of the function f(x,y,zn) at each z=zn plane. Then, the normalized modulus of f(x,y,zn), |f(x,y,zn)|/M, where M is the maximum of |f(x,y,zn)| for all zn, is plotted versus the spatial coordinates x and y to obtain 2D images of the object at all Nz planes. By putting together all 2D slice images, a 3D image of the object is obtained.
5. Image Reconstruction Results
The performance of the proposed multiple receiver planar-scan imaging method is studied through various simulation examples using FEKO [27] and an experimental example. First, we present an example demonstrating the improvement achieved in the multiple receiver configuration compared to the single receiver configuration where Dx=Dy=0. Then, we study the resolution limits of the multiple receiver setup.
As a quick note, the choice of frequency range in microwave imaging of lossy objects is primarily driven by the trade-off between resolution, which improves with increasing frequency and penetration, which usually improves with decreasing frequency. Thus, it has to be chosen properly and in accordance with the properties of the inspected dielectric medium.
5.1. Two-Antenna System versus Six-Antenna System
In the first example, six λ/2 (at 6.5 GHz) x-polarized dipole antennas are employed to scan two cuboids as shown in Figure 4. The objects have a side of 2 mm and center-to-center distance of 8 mm. They are at the range position of z=30 mm. The frequency-independent properties of the background medium are εrb=16 and σb=0.5 S/m while those of the two objects are εro=32 and σo=1S/m. Antenna 2 is aligned with Antenna 1 (transmitter). Antennas 3 to 6 have offsets of ±Dx and ±Dy along the ±x and ±y directions, respectively, as illustrated in Figure 4. Antenna 1 scans an area of 120 mm × 120 mm. Thus, with reference to Figure 3, aperture 1 has a size of 120 mm × 120 mm while aperture 2 has a size of (120mm+2Dx)×(120mm+2Dy). The number of samples along both x and y directions is 61.
Simulation example in FEKO. All dimensions are in mm. The background has properties of εrb=16 and σb=0.5 S/m. Imaging is performed for two objects along the x axis with properties εro=32 and σo=1 S/m.
The transmission S-parameters Sk1 (k = 2 to 6) are acquired over the frequency range from 3 GHz to 10 GHz. The sampling rate in frequency is determined based on the Nyquist criterion as suggested in [11, 19]. The maximum phase shift that results from a change in the wavenumber Δk is ΔkRmax, where Rmax is the maximum distance traveled by the wave from the transmitter to the receiver through the object. Requiring that this phase shift be less than π rad yields
(27)Δf<c2Rmax.
Here, a sampling step of Δf=0.25 GHz is chosen which is approximately c/2Rmax with Rmax=150 mm. In general, based on our studies, decreasing Δf further improves the robustness to noise to some extent.
When calibrated (as described in [19]), the Sk1(fork=2,…,6) parameters provide good approximations of the x-component of the forward-scattered waves at the center of the antennas.
We first show that using the forward-scattered waves acquired in a single receiver configuration when Dx=Dy=0 is not sufficient to implement 3D holography as we discussed in Section 2. Figure 5(a) shows the reconstructed images when only S21 is processed using the 3D holography technique presented in [19]. The quality of the reconstructed images is not satisfactory. Not only the two objects are not resolved well at z=30 mm (their true range location) but also the images at other range locations show the presence of strong artifacts.
Reconstructed images for the objects shown in Figure 4 when applying 3D microwave holography to (a) only S21, (b) all five transmission S-parameters when Dx=Dy=10 mm, and (c) all five transmission S-parameters when Dx=Dy=20 mm. The horizontal and vertical axes show the positions along the x and y directions in terms of mm.
Then, we process the forward-scattered waves sampled by the five offset receivers and apply the reconstruction technique proposed here. Figure 5(b) shows the images obtained when all transmission S-parameters (Sk1, k=2,…,6) are employed and when Dx=Dy=10 mm. As the figures show, the two objects are more discernible in the image at z = 30 mm and the artifacts in the images at other range locations are much weaker compared to those in Figure 5(a). Figure 5(c) shows the images obtained when Dx=Dy=20 mm. It is obvious that increasing the offset distances Dx and Dy improves the quality of the images as expected from the discussion in Section 2.
5.2. Simulation Examples for the Six-Antenna System
We conduct FEKO simulations for the setups illustrated in Figures 6(a)–6(d).
Simulation examples in FEKO. The background has properties of εrb=16 and σb=0.5 S/m. Imaging is performed for (a) a single object with properties εro=32 and σo=1 S/m, (b) two objects with properties εro=32 and σo=1 S/m along the z-axis, (c) two objects with properties εro=48 and σo=2 S/m along the x-axis embedded in a lower contrast medium with properties εrm=32 and σm=1 S/m, and (d) two X-shape objects with properties εro=32 and σo=1 S/m. All dimensions are in mm.
In the first example (Figure 6(a)), a single small object is present, which allows the study of the point-spread function (PSF) of the method. The distance between the two apertures is 70 mm and one cuboid object is placed midway between them. The size and properties of the object and the properties of the background medium are the same as in the previous example. The offset distances are Dx = Dy = 20 mm. The data acquisition process and sampling in space and frequency is the same as in the previous example. To keep the simulations of the 2D scan affordable in terms of computational time, we have used a relatively coarse mesh size of λ/10, where λ is the wavelength at 6.5 GHz in the background medium. To estimate the signal-to-noise ratio (SNR) of the simulations, we also perform a single simulation with a mesh size of λ/30 and evaluate the change in the magnitude of the transmission S-parameter between Antenna 1 and Antenna 2, that is, S21. If we denote S21 obtained with the coarse and fine meshes as S21c and S21f, respectively, then we evaluate the numerical SNR at each frequency ω as follows:
(28)SNR(ω)=20log[|S21f(ω)||S21f(ω)-S21c(ω)|].
Figure 7 shows the variation of SNR with frequency. The minimum SNR is as low as −5 dB at the lower frequencies. The maximum SNR is still low and around 12 dB in the middle of the frequency band. For image reconstruction, we intentionally use the set of data acquired with the coarse mesh to show the performance of the algorithm in the presence of high levels of noise and uncertainties. We note that the variation of the SNR versus frequency is similar in all other simulation examples presented here.
Variation of the numerical SNR versus frequency in the simulations.
Figure 8(a) shows the images reconstructed with Δz = 5 mm and with the Tikhonov regularization parameter α being 2.5×10-5 (the same value is used in the other simulation examples). For brevity we only show images every 10 mm. It is observed that the object is recovered well at z = 35 mm and the artifacts at other range positions are weak. Figures 8(b) and 8(c) show the 1D images of this object along the z-axis and also along the x- and y-axes in the z = 35 mm plane. From the half-power width of these functions, we estimate the range and cross-range resolutions to be ρ-z=12 mm and ρ-u=4 mm. These are close to the values computed from (12) and (14), that is, ρ-zmin=10.8 mm and ρ-umin=5.4 mm. One may argue that the offsets are not very large in this example (Dx=Dy=20 mm) and the resolution limits should not be close to the resolution limits computed using (12) and (14). Similar conclusions are drawn here from another example, where z-=50 mm. We believe that the slight improvement in the resolution is due to the combination of data obtained from various look angles (we have five receivers and for each receiver the data is collected over the aperture). It is wellknown that in the diffraction tomographic imaging techniques based on many viewing angles, the resolution is improved compared to the diffraction-limited resolution [28]).
(a) Reconstructed images (shown every 10 mm along z and only at three range locations) when the distance between the reconstruction planes is Δz = 5 mm for the example in Figure 6(a). Color bars are the same as in Figure 5. Normalized 1D magnitude of the reconstructed images along (b) range direction (z direction) and (c) cross-range directions (x and y directions) in the plane of the object. The red dashed line shows the 3-dB level below the peak (0.7 level on linear scale).
Also, it is observed that the resolution is improved compared to the diffraction-limited resolution in free space. This is due to the wavelength contraction resulting from the high relative dielectric constant of the background medium. For example, the diffraction-limited cross-range resolution in free space is approximately 24 mm while the obtained resolution is 4 mm.
Figure 6(b) shows the setup for another simulation example in which two small objects are along the z-axis positioned at z = 25 mm and z = 45 mm. The properties of the background medium and objects and the dimensions and sampling rates are the same as in the previous examples.
The reconstruction in this case is more challenging due to the multiple scattering between the two objects and because they are at the same x and y coordinates. Figure 9(a) shows the reconstructed images for this example. It is observed that the two objects are reconstructed at their true locations with weak artifacts at other planes. Figure 9(b) shows the magnitude of the image values along the z-axis. It is again observed that the two objects are detected well with the dip between the two peaks dropping below 0.7 level. Figure 9(c) also shows the magnitude of the image values along the x- and y-axes on the z = 25 mm plane (similar variation is observed at z = 45 mm plane).
(a) Reconstructed images (shown every 10 mm along z and only at three range locations) when the distance between the reconstruction planes is Δz = 5 mm for the example in Figure 6(b). Color bars are the same as in Figure 5. Normalized 1D magnitude of the reconstructed images along (b) range direction (z direction) and (c) cross-range directions (x and y directions) in the plane of the object at z = 45 mm.
Figure 6(c) shows the next simulation example in which two cuboids with size of 2 mm × 2 mm × 2 mm are along the x-axis positioned at x=-4 mm and x=4 mm. The properties of the background medium as well as the dimensions, and the spatial and frequency sampling rates are the same as in the previous examples. However, the objects have properties of εro=48 and σo=2 S/m. They are embedded in a lower contrast medium which is shaped like a box of size 14 mm × 8 mm × 2 mm and has properties of εrm=32 and σm=1 S/m. Figure 10(a) shows the reconstructed images in this example. It is observed that the two objects are reconstructed at their true locations with a shadow around them which represents the lower contrast medium. The images at other range locations show some weak artifacts. Figure 10(b) shows the magnitude of the image values along the z-axis. It is observed that the objects are detected well at their true location of z = 35 mm. Figure 10(c) also shows the magnitude of the image values along the x- and y-axes on the z = 35 mm plane. It is observed again that the two objects are detected well at their true locations.
(a) Reconstructed images (shown every 10 mm along z and only at three range locations) when the distance between the reconstruction planes is Δz = 5 mm for the example in Figure 6(c). Color bars are the same as in Figure 5. Normalized 1D magnitude of the reconstructed images along (b) range direction (z direction) and (c) cross-range directions (x and y directions) in the plane of the objects.
Figure 6(d) shows the last simulation example in which two X-shape objects are at the range positions of 25 mm and 45 mm. The properties of the background medium as well as the dimensions and the spatial and frequency sampling rates are the same as in the previous examples. The objects have electrical properties of εro=32 and σo=1S/m. Figure 11 shows the reconstructed images in this example. It is observed that the two objects are reconstructed at their true range locations with a weaker artifact appearing at z = 35 mm. The images at other range locations show some weak artifacts.
Reconstructed images (shown every 10 mm along z and only at three range locations) when the distance between the reconstruction planes is Δz = 5 mm for the example in Figure 6(d). Color bars are the same as those in Figure 5.
5.3. Experimental Results
We also demonstrate the capabilities of the proposed multiple receiver approach through an experiment in free space. We emphasize that confirming the capability of 3D imaging in free space would also imply the capability of 3D imaging in lossy mediums as long as the SNR is sufficient for reliable reconstruction using transmission coefficients. Besides, 3D imaging in free space itself is important when the setup is nonreciprocal and reflection coefficients cannot be acquired (e.g., due to the use of power amplifiers at the transmitter side or low-noise amplifiers at the receiver side).
Here, six X-band open-ended waveguides are employed, one as a transmitter and five as receivers (as shown in Figure 12). The waveguides have a flange of size 40 mm × 40 mm. Thus, for the receivers we have Dx=Dy=40 mm. The distance between the two apertures is z-=110 mm. The object is a very thin copper sheet with size of 40 mm × 20 mm centered at (0,0,58) mm with its edges along the y- and x-axes. The antennas are stationary while the object is scanned over a 2D region of size 160 mm × 200 mm with steps of 5 mm along x and y. The motion of the object is controlled by a computer through two stepper motors. At each scanning step, the five transmission S-parameters (Sk1, k=2,…,6) are recorded by an Advantest R3770 vector network analyzer via an RF switch. The measurements are performed in anechoic chamber to minimize the electromagnetic interference.
(a) Illustration of the experimental imaging setup consisting of six X-band open-ended waveguides where one acts as a transmitter and five act as receivers. A rectangular metallic sheet is being imaged. (b) Image of five open-ended waveguides attached together to construct the multiple receiver array.
The 2D scan is performed two times, once with the object and a second time without the object, in order to calibrate the acquired transmission S-parameters as described in [19].
The structure of a commercially available X-band waveguide, which also contains a flange, is modeled in FEKO to obtain the incident field/Green’s functions for the transmitter/receivers. Also, to better align the strengths of the simulated fields and the ones produced by the waveguide in practice, we calibrate the simulated incident field as
(29)Einc,kcal(r,ω)=1-|Skkmeas(ω)|21-|Skksim(ω)|2Einc,ksim(r,ω),fork=2,…,6,
where Einc,kcal(r,ω) is the calibrated incident field, Einc,ksim(r,ω) is the simulated incident field, r is the position at which the field is evaluated, Skkmeas(ω) is the measured reflection S-parameter, and Skksim(ω) is the simulated reflection S-parameter. In effect, (29) aligns the amplitude of the radiated field in the simulations with that in the measurements.
We perform the measurement over a wide band from 3 GHz to 20 GHz. Although the antennas are not well matched in the whole band, the signature of the object can be observed and processed from the data collected over the entire band. Figure 13 shows the images obtained with the proposed 3D imaging algorithm. Again, reconstruction is performed with Δz=5 mm but for brevity we only show the images every 10 mm. The value of the Tikhonov regularization parameter α is 5×10-2. The object is detected well at z = 60 mm plane (close to its true location of z = 58 mm) while the artifacts on the other planes are relatively low. The strongest artifact is on the plane z = 50 mm since it is closer to the object than the calculated range resolution limit of 17 mm. Figure 14 shows the magnitude of the image values along the z axis. It is observed again that the object is detected well at its true location. Figure 14 also shows the magnitude of the image values along the x- and y-axes on the z = 60 mm plane. It is observed again that the object is detected well at its true location.
Reconstructed images (shown every 10 mm along the z axis) when the distance between the reconstruction planes is Δz = 5 mm using the data in the frequency range from 3 GHz to 20 GHz in the experimental example. Color bars are the same as those in Figure 5. The horizontal and vertical axes show the positions along the x and y directions in terms of mm.
Normalized 1D magnitude of the reconstructed images along the range direction (z direction) and cross-range directions (x and y directions) on the plane of the object in the experiment.
We note that here it is possible to use the waveguides in multimode operation and even below their cut-off frequency. This is because their exact incident field/Green’s function is obtained through simulations and then calibrated through (29). In general, acquisition of data over a wider frequency band leads to better images. For example, Figure 15 shows the images when we use the data collected from 8.5 GHz to 20 GHz. Although this time the lowest frequency is higher than the cut-off frequency of the X-band waveguide (8.2 GHz), the images suffer higher levels of artifacts.
Reconstructed images (shown every 10 mm along z) when the distance between the reconstruction planes is Δz = 5 mm using the data in the frequency range from 8.5 GHz to 20 GHz in the experimental example. Color bars are the same as those in Figure 5. The horizontal and vertical axes show the positions along the x and y directions in terms of mm.
5.4. Computational Complexity of the Inversion Technique
Here, we provide an estimate of the computational complexity of our 3D image reconstruction process. First, we denote the number of samples along x′ and y′ on the acquisition apertures by Nx′ and Ny′. The number of samples along x and y in each reconstructed cross-range plane is denoted by Nx and Ny. We also denote the number of samples of kx and ky by Nkx and Nky, respectively. Table 1 summarizes the computational complexity of our approach. The flops for the FT and the inverse FT operations are provided based on the implementation of fast FT and inverse fast FT operations. The computational complexity of solving the systems of equations has been provided with the assumption that they are solved with QR factorization. The total number of flops for the image reconstruction process is the sum of all the flops in Table 1.
Details of computational complexity of the proposed image reconstruction process.
Operation
Number of flops
FT of the scattered fields
5Nx′Ny′log(Nx′Ny′)
Computing functions ak
5NxNy
FT of functions ak
5NxNylog(NxNy)
Solving the systems of equations in all combinations of kx and ky
NkxNky(10NωNz2+10NωNz+Nz2)
Inverse FT of the contrast function
5NxNylog(NxNy)
The execution time for the examples provided in this paper is typically less than 1 minute on a computer with Intel Xeon CPU 2.93 GHz with 48 GB of RAM.
6. Conclusion and Discussion
We proposed a multiple receiver microwave holography setup to perform 3D imaging with forward-scattered waves only that are acquired with planar raster scanning. This new setup paves the way toward microwave imaging of tissues and other materials of high loss where the back-scattered waves are too weak or are not available due to the nonreciprocal measurement system. The resolution limits for such setup are derived assuming far-field approximation.
The resolution limits proposed here are approximate values derived based on the following assumptions. (1) The background medium is homogeneous, lossless, and nondispersive. In the scenarios where the background medium is not homogenous or dispersive, these expressions can be applied with averaged parameter values for approximate estimation of the resolution limits. If the background medium is lossy, then the maximum Dx and Dy offset values would be limited due to the attenuation of the signal for longer paths. (2) The object is in the far zone of the antennas; that is, we only take into account the propagating waves. In near-field imaging, however, the contribution of evanescent waves leads to better resolution beyond what is computed from the expressions derived here.
Through simulation results, we showed that for dimensions typically used in microwave imaging, the spatial resolution limits provide acceptable image quality. The achieved range and cross-range resolution limits for the proposed setup are comparable or below a centimeter in the frequency range from 3 GHz to 10 GHz.
We also evaluated the performance of the proposed 3D imaging method via an experiment in free space with a 5-element array of X-band open-ended waveguides. The satisfactory results justify further experimental work toward tissue imaging with microwave holography where the sensors must have significantly reduced size and improved sensitivity. In tissue imaging, the SNR values can vary significantly based on (1) the utilized antennas, (2) the acquisition setup and measurement instrumentation, and (3) the properties of the tissue. All of these factors would affect the choice of the proper Δz value as well as the sensitivity of the system. However, these can be determined empirically using simple phantom measurements with the particular hardware setup.
Although the objects are symmetric in the cross-range plane, the reconstructed images are slightly asymmetric. The first factor which affects both the simulation and measurement results is the noise in the scattered field. The simulated scattered field suffers from the numerical noise as discussed in Section 5 while the measured scattered field suffers from mechanical positioning errors during raster scanning, from electronic noise in the measurement equipment and from electromagnetic interference from the environment. In the case of the experimental results, the slight errors in the numerical evaluation of the kernel in (20) intensify this issue further.
As we discussed earlier, increasing the offset values Dx and Dy leads to improvement in the resolution. However, in lossy background media the offsets can be increased as much as the received signals remain reliable. The optimal offsets are determined based on the amount of loss in the background medium, size and contrast of the objects to be detected, radiation pattern of the antennas, and the maximum allowable dimension of the imaging setup. Thus, the optimal offset values would be highly casedependent.
We reiterate that the proposed reconstruction method is based on the linear Born approximation with the entailing limitations. The inclusion of electrically large objects in the imaged region will inevitably lead to artifacts in the reconstructed images. Yet, the large object and other possible inclusions in this object will be detected as demonstrated by a simulation example. In the experimental example, the object is again a large copper sheet and our algorithm is again capable of producing satisfactory images where the object is clearly visible. In general, it is well known that the accuracy of the imaging results with methods based on the Born approximation degrades for high-contrast and large objects, the degradation usually being in the form of increased image artifacts. On the other hand, the algorithms which are not based on the Born approximation (e.g., model-based optimization approaches) suffer from ill-conditioning especially in 3D imaging. The ill-conditioning also makes them prone to errors in the forward models they use, which may be significant in full-wave simulation models. Such algorithms also suffer from nonuniqueness of the solution and often fail to converge to a meaningful solution. From a practical standpoint, such algorithms are also often in a disadvantage because they require extensive computational time and memory. Our proposed 3D imaging is fast, does not require a full-wave forward model, and is significantly less ill-conditioned. Future work is now focused on an approach to tissue imaging, which builds on the result of the presented holography method and overcomes the limitations of the Born approximation.
NikolovaN. K.Microwave imaging for breast cancer20111277894FarhatN. H.GuardW. R.Millimeter wave holographic imaging of concealed weapons1971599138313842-s2.0-001512936910.1109/PROC.1971.8441BaranoskiE. J.Through-wall imaging: historical perspective and future directions200834565565692-s2.0-4584913377110.1016/j.jfranklin.2008.01.005ZoughiR.2000Kluwer Academic PublishersPastorinoM.2010Hoboken, NJ, USAWileyFarhatN. H.Microwave holography and its applications in modern aviationEngineering Applications of Holography Symposium1972295314Proceedings of SPIETricolesG.FarhatN. H.Microwave holography: applications and techniques19776511081212-s2.0-001741798210.1109/PROC.1977.10435HildebrandB. P.HainesK. A.Holography by scanning1969591162-s2.0-0014442238HildebrandB. P.BrendenB. B.1972New York, NY, USAPlenumGaborD.A new microscopic principle194816140987777782-s2.0-34250769340SheenD. M.McMakinD. L.HallT. E.Three-dimensional millimeter-wave imaging for concealed weapon detection2001499158115922-s2.0-003544603510.1109/22.942570SoumekhM.Bistatic synthetic aperture radar inversion with application in dynamic object imaging1991399204420552-s2.0-002622135810.1109/78.134436SoumekhM.A system model and inversion for synthetic aperture radar imaging19921164762-s2.0-0026624352SoumekhM.1994Englewood Cliffs, NJ, USAPrentice-HallMunsonD. C.Jr.O'BrienJ. D.JenkinsW. K.A tomographic formulation of spotlight-mode synthetic aperture radar19831789179252-s2.0-0020802427SheenD.McMakinD.HallT.Near-field three-dimensional radar imaging techniques and applications20104919E83E932-s2.0-7795597709210.1364/AO.49.000E83DallingerA.SchelkshornS.DetlefsenJ.Efficient ω-k-algorithm for circular SAR and cylindrical reconstruction areas2006485912-s2.0-33748414030TanW. X.HongW.WangY. P.WuY. R.A novel spherical-wave three-dimensional imaging algorithm for microwave cylindrical scanning geometries201111143702-s2.0-78751484683AminehR. K.RavanM.KhalatpourA.NikolovaN. K.Three-dimensional near-field microwave holography using reflected and transmitted signals20115912477747892-s2.0-8145515878910.1109/TAP.2011.2165496AminehR. K.KhalatpourA.XuH.BaskharounY.NikolovaN.Three-dimensional near-field microwave holography for tissue imaging201220121110.1155/2012/291494291494LiA.Algorithms for the implementation of Stolt interpolation in SAR processing1Proceedings of the IEEE International Geoscience and Remote Sensing Symposium1992360362SubizaB.Gimeno-NievesE.Lopez-SanchezJ. M.Fortuny-GuaschJ.An approach to SAR imaging by means of non-uniform FFT's6Proceedings of the IEEE International Geoscience and Remote Sensing SymposiumJuly 2003408940912-s2.0-0242456861AminehR. K.RavanM.TrehanA.NikolovaN. K.Near-field microwave imaging based on aperture raster scanning with TEM horn antennas20115939289402-s2.0-7995214142510.1109/TAP.2010.2103009WalterscheidI.BrennerA. R.EnderJ. H. G.Results on bistatic synthetic aperture radar20044019122412252-s2.0-494424152910.1049/el:20045466TikhonovA. N.ArseninV. Y.1977New York, NY, USAWileyhttp://www.mathworks.com/matlabcentral/fileexchange/52-regtoolsEM softwares & systems-S.A. (Pty) Ltd., http://www.feko.info/WehnerD. R.1995Artech House