#### Abstract

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 -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 (, ) 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 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 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 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 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] 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 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 . 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 and 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 , there are two main regions in the inspected volume where the range resolution is poor. These two regions are the middle of the inspected domain, that is, and the direct path connecting the transceivers when it passes right through the scatterer; that is, . 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 , 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 and are close to : Regarding the cross-range resolution, as per (10), it is best when one antenna is on the -axis while the other antenna is at (, ): Under these conditions, the range resolution is

It is worth noting that in reflection holography where , the cross-range and range resolution limits are Assuming 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 cross-range resolution improves if the size of the aperture is large (). On the other hand, the best range resolution limit is obtained when as 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 - or -axis (, ) 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 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 , then the size of aperture 2 is .

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

Assume that at any measurement frequency we know the incident field at any point in the inspected volume due to antenna 1 when it is at . In addition, Green’s function for antenna 2 is known for an -polarized scattering point source at and the -polarized response at . This information can be obtained via simulations. For brevity, we introduce the notations In [20], we show that while can be obtained from simulations, Green’s function can be obtained from using the reciprocity principle. Thus, assuming that antennas 1 to 6 are identical, only one simulation suffices to obtain , from which Green’s functions for antennas 2 to 6 are obtained as Let be the scattered -field received by the th antenna when the transmitting antenna is at . Following the approach presented in [19], is written as where is the contrast function, and are the wavenumbers of the object and the background mediums, respectively, and Notice that in (20), the integration over and can be interpreted as a 2D convolution integral. Thus, the 2D FT of is written as where and are the 2D FTs of and , 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 as where and, for each , Here, is the discretization step along the -axis. In (24) and (25), the arguments and have been omitted for brevity.

The system in (23) is solved at each pair for , , 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 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 , , 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.

#### 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 . 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 (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 mm. The frequency-independent properties of the background medium are and S/m while those of the two objects are and 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 (. The number of samples along both and directions is 61.

The transmission -parameters ( = 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 , where 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 GHz is chosen which is approximately with mm. In general, based on our studies, decreasing further improves the robustness to noise to some extent.

When calibrated (as described in [19]), the 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 is not sufficient to implement 3D holography as we discussed in Section 2. Figure 5(a) shows the reconstructed images when only 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 mm (their true range location) but also the images at other range locations show the presence of strong artifacts.

**(a)**

**(b)**

**(c)**

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 (, ) are employed and when 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 mm. It is obvious that increasing the offset distances and 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).

**(a)**

**(b)**

**(c)**

**(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 , 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 and evaluate the change in the magnitude of the transmission -parameter between Antenna 1 and Antenna 2, that is, . If we denote obtained with the coarse and fine meshes as and , 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 (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 mm and mm. These are close to the values computed from (12) and (14), that is, mm and mm. One may argue that the offsets are not very large in this example ( 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 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)**

**(b)**

**(c)**

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

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

**(a)**

**(b)**

**(c)**

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 mm and 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 and 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 and 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 -axis. It is observed that the objects are detected well at their true location of = 35 mm. Figure 10(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.

**(a)**

**(b)**

**(c)**

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 and . 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 = 35 mm. The images at other range locations show some weak artifacts.

**(a)**

**(b)**

**(c)**

##### 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 mm. The distance between the two apertures is 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 (, ) 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)**

**(b)**

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 is the calibrated incident field, is the simulated incident field, is the position at which the field is evaluated, is the measured reflection -parameter, and 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 mm but for brevity we only show the images every 10 mm. The value of the Tikhonov regularization parameter is . The object is detected well at = 60 mm plane (close to its true location of = 58 mm) while the artifacts on the other planes are relatively low. The strongest artifact is on the plane = 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 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 - and -axes on the = 60 mm plane. It is observed again that the object is detected well at its true location.

**(a)**

**(b)**

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.

##### 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 and on the acquisition apertures by and . The number of samples along and in each reconstructed cross-range plane is denoted by and . We also denote the number of samples of and by and , 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.

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. 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. 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 the utilized antennas, the acquisition setup and measurement instrumentation, and the properties of the tissue. All of these factors would affect the choice of the proper 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 and 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.