Three-Dimensional Microwave Holographic Imaging Employing Forward-Scattered Waves Only

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 nonreciprocalmeasurement 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.


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][13][14][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).

International Journal of Antennas and Propagation
This limits the technique to far-zone imaging.Resampling of the data in   -space is also necessary.Also, in [16][17][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 incidentfield 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 nearfield 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 (  ,   ) 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   is assumed real as in [11]) and is thus capable of processing the evanescent waves available in near-field imaging.Processing larger values for (  ,   ) leads to better spatial resolution of the reconstructed images.Furthermore, in the approach proposed in [11], resampling of the data in  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   ,   , and   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   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 breastcancer diagnostics and artificial glycerin-based phantoms emulating the breast tissues were imaged.Only 2D holographic imaging results were presented employing forwardscattered 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 x L (x  , y  , 0) bistatic SAR [24].We confirm that the availability of forwardscattered 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.

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  = 0 and aperture 2 at  = .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.
In [19,20], both antennas are assumed to be always in the same (  ,   ) position during the scan.Here, we consider the general case where antenna 2 may have an offset of ±  and ±  with respect to antenna 1 along the and -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 nearfield 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 (  ,   , 0) for antenna 1-this implies that antenna 2 is at (  +   ,   +   , )-the total path traveled by the wave through the object is denoted by  and is written as where ( The differential change in , denoted by Δ, due to an incremental change (Δ  ŝ) in the position of the object in any arbitrary direction ŝ, is written as [24] Δ = (ŝ  ⋅ ∇) Δ  , where ∇ is the gradient operator and the superscript  represents the transpose operator.Here, ŝ is a unit vector.
On the other hand, the minimum "measureable" Δ is related to the velocity of the wave in the medium  and the bandwidth of the imaging system  as [24] Δ min ≈   .
From ( 3) and (4), the resolution limit in the direction ŝ, denoted by   , is estimated as [24] 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 , , and  directions, first the gradient of  in a rectangular coordinate system is written as Then, from ( 5) and ( 6), the cross-range resolution   ( =  or ) and the range resolution   are obtained as .
In the previous equations,   and   are obtained when ŝ in ( 5) is substituted with the unit vector in the particular direction, , , or .
We assume that the object is very small and it is positioned on the -axis that is, the object is at (0, 0, ).Then, (7) are simplified as where 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  represents either  or .  and   denote the offsets of the receiver with respect to the transmitter and the position of the transmitter, respectively.
With reference to Figure 2, we first rewrite the resolution limits in (8) in terms of the angles  1 and  2 as 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 (  ,   ) during the 2D scan (as in [19,20]), that is, when   = 0, there are two main regions in the inspected volume where the range resolution   is poor.These two regions are (1) the middle of the inspected domain, that is,  ≈ /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   ̸ = 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:

International Journal of Antennas and Propagation
Regarding the cross-range resolution, as per (10), it is best when one antenna is on the -axis ( 1 = 0) while the other antenna is at  → ∞ (  → ∞,  2 → /2): Under these conditions, the range resolution is It is worth noting that in reflection holography where  2 =  −  1 , the cross-range and range resolution limits are Assuming  ≈ 2  with   being the center frequency of the band, the cross-range resolution can be written as follows: where   is the wavelength at   .Thus, in this case, the crossrange 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 These results are consistent with the expressions for the resolution limits of reflection holography in [11].

Multiple Receiver Setup
In the previous section, we showed that having a nonzero offset distance for the receiver along the or -axis (  ̸ = 0,  = ,) 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 crossrange 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 forwardscattered waves.Antenna 2 is aligned with the transmitter (antenna 1) while the other receivers (antennas 3 to 6) have offset distances of ±  and ±  along the ± and ± 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 2  × 2  , then the size of aperture 2 is 2(  +   ) × 2(  +   ).

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  = 0 and  = , respectively.For simplicity, we assume that the antennas are -polarized (e.g., dipoles oriented along the -axis).Thus, the field can be reasonably approximated by a TM  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  can be viewed as the    element of the full dyadic while the -field is represented by its -component only.From now on we omit the subscript  for brevity.
The system in ( 23) is solved at each (  ,   ) pair for F(  ),  = 1, . . .,   , 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 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  ℎ 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 (  ,   ) are solved, the inverse 2D FT is applied to F(  ,   ,   ),  = 1, 2, . . .,   , to reconstruct a 2D slice of the function (, ,   ) at each  =   plane.Then, the normalized modulus of (, ,   ), |(, ,   )|/, where  is the maximum of |(, ,   )| for all   , is plotted versus the spatial coordinates  and  to obtain 2D images of the object at all   planes.By putting together all 2D slice images, a 3D image of the object is obtained.

Image Reconstruction Results
The performance of the proposed multiple receiver planarscan 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   =   = 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.

Two-Antenna System versus Six-Antenna
System.In the first example, six /2 (at 6.5 GHz) -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  = 30mm.The frequency-independent properties of the background medium are  b  = 16 and  b = 0.5 S/m while those of the two objects are  o  = 32 and  o = 1 S/m.Antenna 2 is aligned with Antenna 1 (transmitter).Antennas 3 to 6 have offsets of ±  and ±  along the ± and ± 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 (120 mm + 2  ) × (120 mm + 2  ).The number of samples along both  and  directions is 61.
The transmission -parameters  1 ( = 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 Δ is Δ max , where  max 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 Here, a sampling step of Δ = 0.25 GHz is chosen which is approximately /2 max with  max = 150 mm.In general, based on our studies, decreasing Δ further improves the robustness to noise to some extent.When calibrated (as described in [19]), the  1 (for  = 2, . . ., 6) parameters provide good approximations of the 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   =   = 0 is not sufficient to implement 3D holography as we discussed in Section 2. Figure 5(a) shows the reconstructed images when only  21 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  = 30 mm (their true range location) but also the images at other range locations show the presence of strong artifacts.
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 -parameters ( 1 ,  = 2, . . ., 6) are employed and when   =   = 10 mm.As the figures show, the two objects are more discernible in the image at  = 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   =   = 20 mm.It is obvious that increasing the offset distances   and   improves the quality of the images as expected from the discussion in Section 2.

Simulation Examples for the Six-Antenna System.
We conduct FEKO simulations for the setups illustrated in Figures 6(a)-6(d).
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   =   = 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 -parameter between Antenna 1 and Antenna 2, that is,  21 .If we denote  21 obtained with the coarse and fine meshes as  c 21 and  f 21 , respectively, then we evaluate the numerical SNR at each frequency  as follows: 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.
Figure 8(a) shows the images reconstructed with Δ = 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  = 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 -axis and also along the and -axes in the  = 35 mm plane.From the half-power width of these functions, we estimate the range and cross-range resolutions to be   = 12 mm and   = 4 mm.These are close to the values computed from ( 12) and ( 14), that is,  min  = 10.8 mm and  min  = 5.4 mm.One may argue that the offsets are not very large in this example (  =   = 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  = 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 International Journal of Antennas and Propagation 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]).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.The reconstruction in this case is more challenging due to the multiple scattering between the two objects and because they are at the same  and  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 -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 and -axes on the  = 25 mm plane (similar variation is observed at  = 45 mm plane).

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   =   = 40 mm.The distance between the two apertures is  = 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 and 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  and .The motion of the object is controlled by a computer through two stepper motors.At each scanning step, the five transmission -parameters ( 1 ,  = 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.
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 -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 where  cal inc, (r, ) is the calibrated incident field,  sim inc, (r, ) is the simulated incident field, r is the position at which the field is evaluated,  meas  () is the measured reflection -parameter, and  sim  () is the simulated reflection -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 Δ = 5 mm but for brevity we only

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   and   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 nearfield 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.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., modelbased 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 illconditioned.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.

Figure 1 :
Figure 1: Microwave holography setup.The grey thick arrows show the directions of 2D scan.

Figure 2 :
Figure2: Simplifying the 3D setup in Figure1to a 2D setup for deriving resolution limits in terms of  1 and  2 .When   = 0, the range resolution is poor in the shaded region.The grey thick arrows show the directions of the scan.

Figure 3 :
Figure 3: Proposed 3D microwave holography setup using forwardscattered 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.

Figure 4 :
Figure 4: Simulation example in FEKO.All dimensions are in mm.The background has properties of  b  = 16 and  b = 0.5 S/m.Imaging is performed for two objects along the  axis with properties  o  = 32 and  o = 1 S/m.

Figure 6 :
Figure 6: Simulation examples in FEKO.The background has properties of  b  = 16 and  b = 0.5 S/m.Imaging is performed for (a) a single object with properties  o  = 32 and  o = 1 S/m, (b) two objects with properties  o  = 32 and  o = 1 S/m along the -axis, (c) two objects with properties  o  = 48 and  o = 2 S/m along the -axis embedded in a lower contrast medium with properties  m  = 32 and  m = 1 S/m, and (d) two X-shape objects with properties  o  = 32 and  o = 1 S/m.All dimensions are in mm.

Figure 6 (
b) shows the setup for another simulation example in which two small objects are along the -axis positioned at  = 25 mm and  = 45 mm.The properties of the background medium and objects and the dimensions and sampling rates are the same as in the previous examples.

10 Figure 7 :
Figure 7: Variation of the numerical SNR versus frequency in the simulations.

Figure 6 (
c) shows the next simulation example in which two cuboids with size of 2 mm × 2 mm × 2 mm are along the -axis positioned at  = −4mm and  = 4mm.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  o  = 48 and  o = 2S/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  m  = 32 and  m = 1S/m.Figure10(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.Figure10(b)shows the magnitude of the image values along the -axis.It is observed that the objects are detected well at their true location of  = 35 mm.Figure10(c) also shows the magnitude of the image values along the and -axes on the  = 35 mm plane.It is observed again that the two objects are detected well at their true locations.

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  o  = 32 and  o = 1 S/m.Figure11shows 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  = 35 mm.The images at other range locations show some weak artifacts.

Figure 9 :
Figure 9: (a) Reconstructed images (shown every 10 mm along  and only at three range locations) when the distance between the reconstruction planes is Δ = 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 ( direction) and (c) cross-range directions ( and  directions) in the plane of the object at  = 45 mm.

Table 1 :
Details of computational complexity of the proposed image reconstruction process.5      log (      ) Computing functions   5    FT of functions   5    log (    ) Solving the systems of equations in all combinations of   and         (10   2  + 10    +  2  ) Inverse FT of the contrast function 5    log (    )