## Abstract

We propose a polarimetric microwave imaging technique that exploits recent advances in computational imaging. We utilize a frequency-diverse cavity-backed metasurface, allowing us to demonstrate high-resolution polarimetric imaging using a single transceiver and frequency sweep over the operational microwave bandwidth. The frequency-diverse metasurface imager greatly simplifies the system architecture compared with active arrays and other conventional microwave imaging approaches. We further develop the theoretical framework for computational polarimetric imaging and validate the approach experimentally using a multi-modal leaky cavity. The scalar approximation for the interaction between the radiated waves and the target— often applied in microwave computational imaging schemes—is thus extended to retrieve the susceptibility tensors, and hence provides additional information about the targets. Computational polarimetry has relevance for existing systems in the field that extract polarimetric imagery, and particular for ground observation. A growing number of short-range microwave imaging applications can also notably benefit from computational polarimetry, particularly for imaging objects that are difficult to reconstruct when assuming scalar estimations.

© 2017 Optical Society of America

## 1. Introduction

Recent advances by the microwave community have resulted in the development of innovative imaging modalities for medical diagnosis [1–6], concealed threat detection [7–13], through-wall imaging [14–17], and non-destructive testing [18–21]. Radio-frequency (RF) bands are particularly suited to these applications, since RF waves can penetrate through many materials that are opaque at optical frequencies. Furthermore, since RF waves are non-ionizing, they are considered safe for human exposure at suitably low power levels.

The potential for efficient, cost-effective, and high-resolution systems that can achieve fast acquisition rates have recently been demonstrated in computational imaging systems based on cavity-backed [22–24] and metasurface [25–30] apertures. These systems radiate pseudo-orthogonal field distributions in transmission and—by exploitation of the reciprocity principle—in reception, to multiplex information and reconstruct an image. The frequency-diverse aperture limits the complexity of the hardware architecture required for real-time high-resolution imaging, obviating the actively controlled components or the need for mechanical motion that is typically required in conventional systems. Recently, frequency-diverse metasurface apertures were adapted for intensity-only measurements, demonstrating that the phase information can be coded in the dispersive radiation of these complex, but passive, antennas and solved for within the inverse problem [31]. This allows further simplification of the associated measurement electronics. However, the application of such computational systems to polarimetric imaging remains unexplored.

Polarization is a relevant source of information in radar imaging, allowing determination of the nature of the interaction between in scene targets and the spatial components of the radiated electromagnetic fields [32]. Polarimetric radars were initially and widely implemented in geoscience applications, allowing for the interrogation of relevant parameters unavailable to single-polarized radars (soil classification, soil moisture, snow distribution, etc.) [33–36]. Following the growing interest in short-range imagery, polarization sensitive radars have been adapted for exploiting the vectorial nature of electromagnetic waves. These approaches thus improve the quality of the reconstructed images, enhancing the accuracy of the estimated contrasts and revealing target features invisible to the approaches based on single-polarized measurements [37–39]. The elimination of artifacts due to multipath in concealed weapon detection has also been demonstrated by illuminating targets using a combinations of left and right-hand circular polarizations [40, 41]. Finally, full-polarimetric imaging has found recent applications in medical diagnosis [42]. However, we note that the technical constraints associated with the development of polarimetric radars are exacerbated when striving for high-resolution and real-time systems, therefore making full reconstruction of polarimetric images extremely burdensome.

In this paper, we propose to extend the framework of microwave computational imaging to the measurement of polarization information in short-range applications by encoding the susceptibility of the target in the physical layer of the antenna, i.e. as a single frequency-dependent signal. In the Section 2 of this article, the theoretical principle of polarimetric computational imaging is first introduced. A generic polarization-sensitive metasurface is considered for simplifying the architecture required for retrieving vectorial information in the target space. The cavity-backed metasurface developed to this end is then presented in the Section 3. A proof of concept of short-range imaging is then presented for validating this novel principle in the Section 4. Finally Section 5 provides the concluding remarks.

## 2. Computational polarimetry: the interrogation of a susceptibility tensor from a signal coded in the physical layer

We first describe the underlying principle of polarimetric microwave imaging, using the illustrative geometry shown in Fig. 1.

In this conceptual example, two arrays of dual-polarized antennas are used to measure the interaction between radiated waves and a target. In the configuration shown, the field is generated by the transmit array and the scattered fields are detected by the receive array. The electric susceptibility tensor *χ̄* (*r*) accounts for the interaction between the radiated waves and the illuminated scene. In a conventional computational imaging approach, the field radiated and measured by a pair of transmit and receive antennas would interact with the target susceptibility, generating the frequency-dependent signal *S̄* (*r** _{t}*,

*r**,*

_{r}*ν*). This 2 × 2 matrix contains the interaction between the excited

*x*or

*z*-oriented element in transmission at location

*r**with the*

_{t}*x*or

*z*-oriented element in reception at location

*r**, where*

_{r}*x*and

*z*are coordinates in the plane of the arrays (both arrays are assumed to be in the same plane). In the first Born approximation, the interaction between the dyadic fields radiated by the

*x*-polarized antenna in transmission, with the

*z*-polarized antenna in reception is:

The vector fields *E** _{x}*(

*r**,*

_{t}*r*,

*ν*) and

*E**(*

_{z}

*r**,*

_{r}*r*,

*ν*) are radiated by each of the transmit and receive antennas from their respective locations

*r**and*

_{t}

*r**. The total measured frequency-dependent signal can be expressed in matrix form for each dual-polarized transmit and receive pair as:*

_{r}*Ē*is a 2 × 3 matrix representing the three-dimensional vector field radiated by two co-localized transverse polarized antennas.

Assuming that these fields are determined either by an experimental measurement or from a numerical analysis, the susceptibility tensor can then be interrogated by solving the inverse problem corresponding to the forward model specified by Eq. (2). Controlling spatially the radiated polarized fields over two apertures is particularly burdensome in applications where real-time and high-resolution constraints are imposed. The RF architectures associated with each array would thus have to include as many active chains as dual-polarized antennas, working simultaneously over an ultra-wide bandwidth. Fast switches would also need to be implemented for selecting the excited polarization of each antenna leading to three-dimensional vector fields.

In this paper, the polarimetric imaging approach is simplified by implementing a computational technique (Fig. 2).

The two arrays are replaced by metasurfaces, which radiate pseudo-orthogonal patterns in space and frequency. In this way, the large number of active chains implemented for the independent transmission and reception of waves (multiplied by two for a polarization sensitive system) is substantially decreased. Instead a single compressed signal can be measured through the input and output ports and the radiating antennas. The amount of information included in a set of these measurements directly depends on the spatial diversity of the fields emitted and received by the antennas [28], and is now studied considering the potential contribution of polarization diversity.

The expression of the fields ** E**(

*r*_{t,r},

**,**

*r**ν*) radiated respectively by the transmit and receive arrays is defined by the expression of the equivalent vector potentials

*A*_{t,r}(

*r*_{t,r},

*ν*), according to the corresponding current distributions

*J*_{t,r}(

*r*_{t,r},

*ν*) in the radiating apertures that

*A*_{t,r}(

*r*_{t,r},

*ν*) satisfies the Helmoltz equation [32]:

*A*_{t,r}(

**,**

*r**ν*) is thus obtained by computing the convolution product between the current distribution

*J*_{t,r}(

*r*_{t,r},

*ν*) and a scalar free-space Green’s function

*g*(

**,**

*r**ν*) in the antenna space

*r*_{t,r}:

*k*= 2

*πν/c*is the wavenumber and

*Ḡ*(

*r**,*

_{a}

*r**,*

_{b}*ν*) is the dyadic Green’s function modeling the propagation from

*r**to*

_{a}

*r**, written as:*

_{b}The susceptibility tensor *χ̄* (** r**) can be estimated from the interaction between the fields radiated by two metasurfaces [28], leading to the following measured frequency signal:

*Ḡ*(

**,**

*r*

*r*_{t},

*ν*) =

*Ḡ*(

**,**

*r*

*r*_{t},

*ν*)

*. As demonstrated in [28], the outer product*

^{T}

*E*_{r}(

**,**

*r**ν*)

*E*_{t}(

**,**

*r**ν*)

*forms a basis interrogating the reflectivity of the scene. In previous studies of computational imaging using a signal coded by the physical layer, a scalar approximation was used for the susceptibility tensor, allowing for an estimation of a target reflectivity illuminated by several metasurfaces. Here, we extend this framework to a tensorial reconstruction of the susceptibility in the target space, passively encoding the three-dimensional polarization information into a unique frequency signal measured between two ports. This signal is simplified when restricting the study to anisotropic media defined by a diagonalized susceptibility tensor. From a fundamental perspective, we consider electrons within the target medium interacting with the incident electric field as being represented by a simple mechanical model (Fig. 3).*

^{T}In the depicted model, the principal axes of the bound electron (*x′*, *y′*, *z′*) are not aligned with the coordinate axes (*x*, *y*, *z*), the latter defined as references for the transmitted and received electromagnetic fields. Thus, the non-diagonal terms are created by projections of the radiated field over the main axis of the tensor at each location, leading to the following eigendecomposition:

*R̄*(

**) is a rotation matrix and**

*r**ξ*(

**) is the diagonalized susceptibility tensor. Assuming the first Born approximation, a forward model is defined that links the measurements to the susceptibility tensor defined at each location of the target space. The principal axes of each tensor can thus be estimated from a diagonalization determining the orientation of the currents on the targets, which is homogenized over the size of the reconstructed voxels. Finally, These estimations can be represented using the weighted eigenvectors or equivalent ellipsoids, showing the anisotropy of the scattering (Fig. 4). Similar studies are notably developed for the estimation of anisotropy in the the field of diffusion tensor imaging [43–45].**

*r*Having introduced the essential aspects of the computational imaging, we take a closer look at the radiation of polarized fields by metasurfaces. Indeed, simplifications can be applied when considering antennas radiating in only one direction of the space. The radiated field can thus be deduced by the measurement of the tangential electric field on a synthetic aperture plane with vanishing boundary conditions, denoted ${\mathit{E}}_{\text{t},\text{r}}^{\text{tan}}(\mathit{r},\nu )$ [46]. From this field, we can extract a set of equivalent dipoles, as

*p*_{t,r}(

*r*_{t,r}) are the equivalent electric dipoles and

*m*_{t,r}(

*r*_{t,r}) the magnetic ones, computed in transmission and reception over small 2D domains

*r′*_{t,r}, effectively defining a grid where each facet is centered at

*r*_{t,r}. The vector

**is defined as the normal to the measurement plane and**

*n***is a unit vector perpendicular to the plane (optical axis). For the sake of simplicity, the transmit and receive metasurfaces can be defined as co-located and oriented such that**

*ŷ***=**

*n***for each radiating plane. The electric dipoles**

*ŷ*

*p*_{t,r}are thus canceled out by the dot product $\widehat{\mathit{y}}\cdot {\mathit{E}}_{\text{t},\text{r}}^{\text{tan}}=0$, allowing for a decomposition of the field in the target space as the superposition of contributions radiated by magnetic dipoles only, as depicted in Fig. 5.

In this example, different distributions of random magnetic dipoles are radiating according to the excitation port and to the frequency. The expression of the radiated electric field in the target space is finally expressed from the contributions of all the magnetic dipoles, summing over the discrete dipole space *r*_{t,r} [46,47]:

The projection property of the cross product allows for a random distribution of magnetic sources contained in a plane to be sufficient for radiating 3D polarized electric near fields that are able to interrogate the full susceptibility tensor. A pseudo-random distribution of sources in the flat metasurface aperture thus radiates a complex field able to probe the medium of interest, encoding the susceptibility variation of the target space into a unique frequency-dependent signal. This system is able to passively multiplex the target space information in transmission and reception, extending earlier approaches based on a scalar approximation of the susceptibility tensor in [24,48].

The susceptibility of the target space can be estimated by a pseudo-inversion of the radiated fields in transmission and reception of the form

*E*

_{t,r}(

**,**

*r**ν*)

^{+}is the pseudo-inverse computed for each independent polarization state on the discrete representation of the radiated fields in transmission and reception [49–51], determined from Eqs. (12) and (13) and assuming that the tangential field ${\mathit{E}}_{\text{t},\text{r}}^{\text{tan}}(\mathit{r},\nu )$ can be determined from an analytical model (or measured).

We stress that these estimations can only be achieved efficiently if these fields present a low level of correlation, spatially, and for each polarization state. Conveniently, we can exploit frequency diversity to encode the measured information onto the ports connected to the metasurface antenna.

To be more specific, we propose to develop the expression of the sub-element ${\widehat{\overline{\chi}}}_{(m,n)}$ of the retrieved susceptibility tensor at the location *r*_{0}:

*R*

_{(m,n,p,q)}(

*r*

_{0}) is a function of the degree of correlation in space and polarization of the radiated fields, and is defined as:

An ideal estimation of the full susceptibility tensor at this specific location would thus impose to ensure that

In practice, the sum of non-zero cross-terms can thus introduce a speckle noise impacting the accuracy of the estimation. In general, the reconstructions will thus highly depend on the technique implemented for solving the inverse problem formulated in Eq. (15), using direct or iterative strategies of reconstruction. In this article, all the results were computed by implementing an iterative least-square algorithm.

## 3. Cavity-backed polarization-sensitive metasurface

We validate the computational polarimetric imaging technique using the frequency-diverse, chaotic leaky microwave cavity illustrated in Fig. 6.

The radiated patterns must exhibit a low level of correlation among the radiated field patterns in frequency, space, and polarization—ensuring an appropriate multiplexing of the polarimetric information gathered from the target space and encoding it into a single frequency vector. These properties are obtained by exploiting the modal diversity in a leaky cavity designed to be several times larger than the operating wavelength. The frequency bandwidth is thus set between 17.5 and 26.5 GHz, corresponding to wavelengths of 1.15 cm to 1.67 cm, for a cavity with inner dimensions of 15.2 × 28.5 × 28.5 cm^{3}. This type of over-sized reverberating cavity finds applications in various domains such as the evaluation of the gain of antennas or for electromagnetic vulnerability measurements [52,53]. A metallic ball is inserted into this air-filled cavity to optimize the modal diversity, to create irregular and convex boundaries and prevent the development of degenerate states [54,55]. As opposed to previous microwave computational imaging demonstrations [23,31] also based on air-filled leaky cavities, this work introduces a means of estimating the vector response of the target space instead of its scalar approximation. In addition to the low level of correlation of the fields in space and frequency that is required in the previous experiments, we also require an additional constraint of pseudo-orthogonality in polarization. The various correlations in frequency, space and polarization can be probed by analyzing the near field scans of ${\mathit{E}}_{1,2}^{\text{tan}}({\mathit{r}}_{\text{t},\text{r}},\nu )$ obtained by exciting the Ports 1 and 2 connected to the cavity successively (Fig. 7).

For the measurements, a commercial near-field scanner was used, NSI 200V-3×3 [23, 27]. The measured field distribution is equivalent to the one depicted in Fig. 5, since the tangential fields and the magnetic dipoles are directly related (Eq. 12). One can observe that both the location of the feeding port and the frequency, even when the shift represents a tiny fraction of the operating frequency, have a noticeable impact on the near field distribution in the aperture of the metasurface. Before the introduction of any imaging experiment, the correlation between the radiated fields is studied. This analysis is initiated by considering only the field radiated on the *x*-polarization when feeding port 1. Near field scans are recorded over 4001 frequency points between 17.5 GHz and 26.5 GHz, leading to a frequency sampling interval of 2.25 MHz. Pearson’s formula is used for the computation of the correlation, computed between two scalar zero-mean fields *E*_{1}(*r** _{t}*,

*ν*) and

*E*

_{2}(

*r**,*

_{t}*ν*) as:

*E*] and

_{i}*σ*are the expectation and the standard deviation of

_{i}*E*, respectively, computed over the frequency dimension for each spatial coordinate

_{i}**. The autocorrelation of the first scan is depicted in Fig. 8.**

*r*The pre-processed fields, initially back-propagated to a 18.65 × 18.65 cm^{2} area and sampled using 63 points over the *x* and *z* axes, are downsampled to 16 × 16 points per scans for obtaining a spatial grid of 12 mm, close to the central wavelength of the operating bandwidth *λ _{mean}* = 13.6 mm. The two-dimensional coordinate

*r*is vectorized on each axis of this figure. The depicted result, corresponding mainly to a diagonal line linking the components sharing the same index, denotes the low level of spatial correlation in this scan obtained by exploiting the frequency diversity in the same way as in a computational imaging experiment. Two low level parallel line are obtained on both sides of the central diagonal, corresponding to the spatial correlation of adjacent points. This phenomenon has already been partially mitigated by resampling the field on a grid closer to the operating wavelength. This study is then extended to both ports and both polarizations, leading to 16 combinations of correlations given in Fig. 9.

_{t}This global representation of the spatial correlation of the near fields for each couple of ports and polarizations presents the properties expected for interrogating the susceptibility tensor in the target space: the four autocorrelations are close to diagonal matrices, while all the other cross correlations are low. In this way, the polarization information is encoded by the radiation of the metasurface in transmission and reception, and can be retrieved from the inverse formulation by exploiting the frequency diversity.

## 4. Practical implementation

We investigate the polarimetric imaging capabilities of the system using a set of 2 mm thick copper wires as 8 × 12 cm^{2} targets, bent to form the four letters of “DUKE” shown in Fig. 10. A target made of thin conductive elements is expected to have an anisotropic scattering according to the orientation of the wires, revealed by the proposed computational imaging system. The letters are positioned at a distance of 30 cm in front of the radiating aperture and sequential measurements are performed using a vector network analyzer connected to ports 1 and 2 of the cavity. For each measured frequency signal, the inverse problem stated in Eq. 15 is computed using an iterative generalized minimal residual method.

The estimation of the first susceptibility tensor $\widehat{\overline{\chi}}$ is shown in Fig. 11 where, first, only the magnitude is displayed. The computation is performed for the 3 × 3 interactions between all the polarization components in a target space discretized in 51200 voxels. The transverse axes ** x̂** and

**are both sampled on a 2.5 mm grid of 80 elements and the optical axis**

*ẑ***is sampled on a 5 mm grid of 8 elements. The susceptibility tensor is thus a complex matrix of dimensions $\widehat{\overline{\chi}}\in {\u2102}^{3\times 3\times 80\times 8\times 80}$. The cross-polarized symmetric elements of the tensor are averaged to facilitate numerical processing and to reduce the speckle level.**

*ŷ*The two transverse co-polarized terms “xx” and “zz” have higher amplitudes and allow for the reconstruction of a continuous target. Interestingly, the contribution of the co-polarized “yy” term is quite significant despite a higher level of speckle. To the best of our knowledge, this feature has never been exploited in conventional polarimetric microwave imaging where most applications are restricted to far-field synthetic aperture radars. In a short-range scenario, the paraxial approximation is no longer valid—leading to an increase amount of measurable information projected along the optical axis. The remaining elements of the susceptibility tensor have a smaller amplitude since they correspond to polarization conversions. This first representation only allows for a study of the magnitude of the tensor in the 3D target space. The resolution of the reconstructed images is directly defined, as in conventional radar imaging, by the size of the radiating aperture and the operating bandwidth [56].This analysis is thus continued by extracting the elements in the target plane to consider the impact of both the magnitude and the phase distributions (Fig. 12).

#### 4.1. Correlation of the co-polarized transversed terms

The results of the two transverse co-polarized images “xx” and “zz” seem similar, but noticeable phase differences are obtained on the horizontal and vertical parts of the target due to the anisotropic nature of the scattering on these thin copper wires. A simple visualization of the anisotropy is thus proposed by computing the correlation *χ̂ _{corr}* between the two co-polarized transverse terms

*χ̂*

_{x,x}and

*χ̂*

_{z,z}:

The result is presented in Fig. 13, controlling the opacity and the color of each pixel respectively with the magnitude and the phase of *χ̂ _{corr}*. Similar methods are applied in the field of synthetic aperture radar polarimetry by computing covariance matrices for reducing the speckle noise and extracting relevant soil parameters [35,57].

With this simple processing, the orientation of the wires can easily be retrieved according to the phase value, i.e. close to −*π*/2 for the vertical direction and *π*/2 for the horizontal ones. The elements oriented diagonally are represented by the white color standing for a 0 radian phase. However, this approach does not allow for the discrimination of the direction of the diagonal orientations of the wires. A comparable analysis can be performed for a three-dimensional visualization of *χ̂ _{corr}*, displaying a −20 dB isosurface of the reconstructed data and using a color coded representation of the phase (Fig. 14).

Once again, the interference between the two transverse co-polarized terms allows for the reduction of the speckle noise and the visualization of the anisotropic scattering. However, such an approach does not make it possible to take into account all the information contained in the reconstructed susceptibility tensor. A more complex processing is studied based on the eigendecomposition of the susceptibility tensor to directly identify the main axis of each voxel.

#### 4.2. Eigendecomposition of the susceptibility tensor

The ratio between the small distance of the targets (30 cm) and the size of the radiating aperture (15 cm) makes it possible to obtain polarimetric components oriented along the optical axis. An eigendecomposition is performed on the real and imaginary parts of the retrieved susceptibility tensor. Each voxel of the target space is represented by a 3D ellipsoid oriented and scaled according to the associated eigenvalues and eigenvectors. The ellipsoids are color-coded according to the orientation of the main axis. The results extracted from the target plane are depicted in Fig. 15.

In contrast to the previous section where the phase information was used to represent the anisotropic nature of the scattering, this additional analysis relies on the amplitude of the set of components composing the susceptibility tensor, where the real and imaginary parts are processed independently. Some ellipsoids reconstructed on the targets clearly correspond to anisotropic signatures, with an elongated shape corresponding to a eigenvalue larger than the other two, while some voxels have a signature closer to an isotropic response, represented by more spherical ellipsoids. Although the principal axes of the ellipsoids are not necessarily all oriented in the directions of the copper wires, it is possible to discern a direct correlation between the shape of the targets and the colors of the ellipsoids. In particular, the red elements oriented along the optical axis mainly occur for diffraction on the edges of the targets.

## 5. Conclusion

A polarimetric microwave imaging technique has been presented in this paper by extending the computational principle of the scalar approaches that have been previously developed in the literature. The theoretical principle proposed in the introduction suggests a key simplification of the architectures for RF systems that seek to achieve high-resolution images from a single or few source system. The leaky multi-modal cavity used to produce pseudo-orthogonal field patterns in frequency and polarization can be generalized to many metasurface aperture paradigms—all of which are capable of generating field patterns with low correlation as a function of frequency or other parameters. The multiplexing of an increased quantity of information compared to the scalar case requires the development of radiating structures with high quality factors. It is possible in this context to imagine the application of such a technique using a set of metasurfaces operating in a cooperative manner in order to reduce the radiation efficiency and modal diversity constraints imposed on the unique element presented in this proof of concept. In the present example, measurements were taken between two ports of a multi-modal cavity over the band 17.5–26.5 GHz. From the frequency-index measurements, it was possible to reconstruct targets made of copper wires forming the word “DUKE”, demonstrating the different modes of operation permitted by the proposed polarimetric approach. Through the two last sections, the anisotropic behavior of the scattering has been highlighted by computing the correlation between the tranverse co-polarized components of the susceptibility tensor, which also allows the reduction of speckle noise, and by computing independent eigendecompositions of the real and imaginary parts of the full susceptibility tensor..

## Funding

Air Force Office of Scientific Research (AFOSR, Grant No. FA9550-12-1-0491).

## References and links

**1. **E. C. Fear, X. Li, S. C. Hagness, and M. A. Stuchly, “Confocal microwave imaging for breast cancer detection: Localization of tumors in three dimensions,” IEEE Trans. Biomed. Eng. **49**, 812–822 (2002). [CrossRef] [PubMed]

**2. **T. Rubæk, P. M. Meaney, P. Meincke, and K. D. Paulsen, “Nonlinear microwave imaging for breast-cancer screening using gauss–newton’s method and the cgls inversion algorithm,” IEEE Trans. Antennas Propag. **55**, 2320–2331 (2007). [CrossRef]

**3. **T. C. Williams, J. M. Sill, and E. C. Fear, “Breast surface estimation for radar-based breast imaging systems,” IEEE Trans. Biomed. Eng. **55**, 1678–1686 (2008). [CrossRef] [PubMed]

**4. **N. K. Nikolova, “Microwave imaging for breast cancer,” IEEE Microw. Mag. **12**, 78 (2011).

**5. **G. Bellizzi, O. M. Bucci, and I. Catapano, “Microwave cancer imaging exploiting magnetic nanoparticles as contrast agent,” IEEE Trans. Biomed. Eng. **58**, 2528–2536 (2011). [CrossRef] [PubMed]

**6. **A. Abbosh, B. Mohammed, and K. Bialkowski, “Differential microwave imaging of the breast pair,” IEEE Antennas Wireless Propag. Lett **15**, 1434–1437 (2016). [CrossRef]

**7. **D. M. Sheen, D. L. McMakin, and T. E. Hall, “Three-dimensional millimeter-wave imaging for concealed weapon detection,” IEEE Trans. Microw. Theory Techn. **49**, 1581–1592 (2001). [CrossRef]

**8. **D. Sheen, D. McMakin, and T. Hall, “Near-field three-dimensional radar imaging techniques and applications,” Appl. Opt. **49**, E83–E93 (2010). [CrossRef] [PubMed]

**9. **X. Zhuge and A. G. Yarovoy, “A sparse aperture mimo-sar-based uwb imaging system for concealed weapon detection,” IEEE Trans. Geosci. Remote Sens. **49**, 509–518 (2011). [CrossRef]

**10. **S. S. Ahmed, A. Schiessl, and L.-P. Schmidt, “A novel fully electronic active real-time imager based on a planar multistatic sparse array,” IEEE Trans. Microw. Theory Techn. **59**, 3567–3576 (2011). [CrossRef]

**11. **S. S. Ahmed, A. Schiessl, F. Gumbmann, M. Tiebout, S. Methfessel, and L. Schmidt, “Advanced microwave imaging,” IEEE Microw. Mag. **13**, 26–43 (2012). [CrossRef]

**12. **Y. Rodriguez-Vaqueiro, Y. A. Lopez, B. Gonzalez-Valdes, J. A. Martinez, F. Las-Heras, and C. M. Rappaport, “On the use of compressed sensing techniques for improving multistatic millimeter-wave portal-based personnel screening,” IEEE Trans. Antennas Propag. **62**, 494–499 (2014). [CrossRef]

**13. **B. Gonzalez-Valdes, Y. Álvarez, Y. Rodriguez-Vaqueiro, A. Arboleya-Arboleya, A. García-Pino, C. M. Rappaport, F. Las-Heras, and J. A. Martinez-Lorenzo, “Millimeter wave imaging architecture for on-the-move whole body imaging,” IEEE Trans. Antennas Propag. **64**, 2328–2338 (2016). [CrossRef]

**14. **E. J. Baranoski, “Through-wall imaging: Historical perspective and future directions,” J. Franklin Inst. **345**, 556–569 (2008). [CrossRef]

**15. **M. Dehmollaian, M. Thiel, and K. Sarabandi, “Through-the-wall imaging using differential sar,” IEEE Trans. Geosci. Remote Sens. **47**, 1289–1296 (2009). [CrossRef]

**16. **T. S. Ralston, G. L. Charvat, and J. E. Peabody, “Real-time through-wall imaging using an ultrawideband multiple-input multiple-output (mimo) phased array radar system,” in “Phased Array Systems and Technology (ARRAY), 2010 IEEE International Symposium on,” (IEEE, 2010), pp. 551–558.

**17. **Y. Wang and A. E. Fathy, “Advanced system level simulation platform for three-dimensional uwb through-wall imaging sar using time-domain approach,” IEEE Trans. Geosci. Remote Sens. **50**, 1986–2000 (2012). [CrossRef]

**18. **S. Caorsi, A. Massa, M. Pastorino, and M. Donelli, “Improved microwave imaging procedure for nondestructive evaluations of two-dimensional structures,” IEEE Trans. Antennas Propag. **52**, 1386–1397 (2004). [CrossRef]

**19. **M. Benedetti, M. Donelli, A. Martini, M. Pastorino, A. Rosani, and A. Massa, “An innovative microwave-imaging technique for nondestructive evaluation: applications to civil structures monitoring and biological bodies inspection,” IEEE Trans. Instrum. Meas. **55**, 1878–1884 (2006). [CrossRef]

**20. **S. Kharkovsky and R. Zoughi, “Microwave and millimeter wave nondestructive testing and evaluation-overview and recent advances,” IEEE Instrum. Meas. Mag. **10**, 26–38 (2007). [CrossRef]

**21. **U. Hasar, “Non-destructive testing of hardened cement specimens at microwave frequencies using a simple free-space method,” NDT and E Int. **42**, 550–557 (2009). [CrossRef]

**22. **D. Carsenat and C. Decroze, “Uwb antennas beamforming using passive time-reversal device,” IEEE Antennas Wireless Propag. Lett **11**, 779–782 (2012). [CrossRef]

**23. **T. Fromenteze, O. Yurduseven, M. F. Imani, J. Gollub, C. Decroze, D. Carsenat, and D. R. Smith, “Computational imaging using a mode-mixing cavity at microwave frequencies,” Appl. Phys. Lett. **106**, 194104 (2015). [CrossRef]

**24. **T. Fromenteze, E. L. Kpré, D. Carsenat, C. Decroze, and T. Sakamoto, “Single-shot compressive multiple-inputs multiple-outputs radar imaging using a two-port passive device,” IEEE Access **4**, 1050–1060 (2016). [CrossRef]

**25. **J. Hunt, T. Driscoll, A. Mrozack, G. Lipworth, M. Reynolds, D. Brady, and D. R. Smith, “Metamaterial apertures for computational imaging,” Science **339**, 310–313 (2013). [CrossRef] [PubMed]

**26. **G. Lipworth, A. Mrozack, J. Hunt, D. L. Marks, T. Driscoll, D. Brady, and D. R. Smith, “Metamaterial apertures for coherent computational imaging on the physical layer,” J. Opt. Soc. Am. A **30**, 1603–1612 (2013). [CrossRef]

**27. **O. Yurduseven, J. N. Gollub, D. L. Marks, and D. R. Smith, “Frequency-diverse microwave imaging using planar mills-cross cavity apertures,” Opt. Express **24**, 8907–8925 (2016). [CrossRef] [PubMed]

**28. **D. L. Marks, J. Gollub, and D. R. Smith, “Spatially resolving antenna arrays using frequency diversity,” J. Opt. Soc. Am. A **33**, 899–912 (2016). [CrossRef]

**29. **J. Gollub, O. Yurduseven, K. Trofatter, D. Arnitz, M. Imani, T. Sleasman, M. Boyarsky, A. Rose, A. Pedross-Engel, H. Odabasi, T. Zvolensky, G. Lipworth, D. Brady, D. L. Marks, M. S. Reynolds, and D. R. Smith, “Large metasurface aperture for millimeter wave computational imaging at the human-scale,” Sci. Rep. **7**, 42650 (2017). [CrossRef] [PubMed]

**30. **O. Yurduseven, D. L. Marks, T. Fromenteze, J. N. Gollub, and D. R. Smith, “Millimeter-wave spotlight imager using dynamic holographic metasurface antennas,” Opt. Express **25**, 18230–18249 (2017). [CrossRef] [PubMed]

**31. **T. Fromenteze, X. Liu, M. Boyarsky, J. Gollub, and D. R. Smith, “Phaseless computational imaging with a radiating metasurface,” Opt. Express **24**, 16760–16776 (2016). [CrossRef] [PubMed]

**32. **W. C. Chew, *Waves and fields in inhomogeneous media*, vol. 522 (IEEE press, 1995).

**33. **F. T. Ulaby and C. Elachi, *Radar polarimetry for geoscience applications*, vol. 1 (Artech House, 1990).

**34. **R. N. Treuhaft and P. R. Siqueira, “Vertical structure of vegetated land surfaces from interferometric and polarimetric radar,” Radio Sci. **35**, 141–177 (2000). [CrossRef]

**35. **J.-S. Lee and E. Pottier, *Polarimetric radar imaging: from basics to applications* (CRC, 2009). [CrossRef]

**36. **D. Velotto, F. Nunziata, M. Migliaccio, and S. Lehner, “Dual-polarimetric terrasar-x sar data for target at sea observation,” IEEE Geosci. Remote Sens. Lett. **10**, 1114–1118 (2013). [CrossRef]

**37. **R. K. Amineh, A. Khalatpour, and N. K. Nikolova, “Three-dimensional microwave holographic imaging using co-and cross-polarized data,” IEEE Trans. Antennas Propag. **60**, 3526–3531 (2012). [CrossRef]

**38. **R. Salman, I. Willms, L. Reichardt, T. Zwick, and W. Wiesbeck, “On polarization diversity gain in short range uwb-radar object imaging,” in “IEEE Int. Conf. UWB,” (IEEE, 2012), pp. 402–406.

**39. **J. Wang, P. Aubry, and A. Yarovoy, “A novel approach to full-polarimetric short-range imaging with copolarized data,” IEEE Trans. Antennas Propag. **64**, 4733–4744 (2016). [CrossRef]

**40. **D. M. Sheen, D. L. McMakin, W. M. Lechelt, and J. W. Griffin, “Circularly polarized millimeter-wave imaging for personnel screening,” in “Defense and Security,” (International Society for Optics and Photonics, 2005), pp. 117–126.

**41. **D. M. Sheen, D. L. McMakin, W. M. Lechelt, and J. W. Griffin, “Circularly polarized millimeter-wave imaging for personnel screening,” in “Proc. SPIE,”, vol. 5789 (2005), vol. 5789, pp. 117–126.

**42. **A. Abubakar, P. M. Van den Berg, and J. J. Mallorqui, “Imaging of biomedical data using a multiplicative regularized contrast source inversion method,” IEEE Trans. Microw. Theory Techn. **50**, 1761–1771 (2002). [CrossRef]

**43. **P. J. Basser, J. Mattiello, and D. LeBihan, “Mr diffusion tensor spectroscopy and imaging,” Biophys. J. **66**, 259–267 (1994). [CrossRef] [PubMed]

**44. **D. Le Bihan, J.-F. Mangin, C. Poupon, C. A. Clark, S. Pappata, N. Molko, and H. Chabriat, “Diffusion tensor imaging: concepts and applications,” J. Mag. Res. Imag. **13**, 534–546 (2001). [CrossRef]

**45. **S. Mori and J. Zhang, “Principles of diffusion tensor imaging and its applications to basic neuroscience research,” Neuron **51**, 527–539 (2006). [CrossRef] [PubMed]

**46. **J. D. Jackson, *Classical Electrodynamics* (Wiley, 1999).

**47. **G. Lipworth, A. Rose, O. Yurduseven, V. R. Gowda, M. F. Imani, H. Odabasi, P. Trofatter, J. Gollub, and D. R. Smith, “Comprehensive simulation platform for a metamaterial imaging system,” Appl. Opt. **54**, 9343–9353 (2015). [CrossRef] [PubMed]

**48. **T. Fromenteze, C. Decroze, and D. Carsenat, “Passive uwb multiplexing device for a single-port mimo radar,” in “IEEE Int. Microw. Symp.”, (2014), pp. 1–4.

**49. **R. Barrett, M. Berry, T. F. Chan, J. Demmel, J. Donato, J. Dongarra, V. Eijkhout, R. Pozo, C. Romine, and H. Van der Vorst, *Templates for the solution of linear systems: building blocks for iterative methods* (SIAM, 1994).

**50. **O. Yurduseven, M. F. Imani, H. Odabasi, J. Gollub, G. Lipworth, A. Rose, and D. R. Smith, “Resolution of the frequency diverse metamaterial aperture imager,” Prog. Electromag. Res. **150**, 97–107 (2015). [CrossRef]

**51. **T. Fromenteze, E. L. Kpré, C. Decroze, D. Carsenat, O. Yurduseven, M. Imani, J. Gollub, and D. R. Smith, “Unification of compressed imaging techniques in the microwave range and deconvolution strategy,” in “Eur. Conf. Radar,” (IEEE, 2015), pp. 161–164.

**52. **P.-S. Kildal, K. Rosengren, J. Byun, and J. Lee, “Definition of effective diversity gain and how to measure it in a reverberation chamber,” Microw. Opt. Technol. Lett **34**, 56–59 (2002). [CrossRef]

**53. **D. A. Hill and J. M. Ladbury, “Spatial-correlation functions of fields and energy density in a reverberation chamber,” IEEE Trans. Electromagn. Compat **44**, 95–101 (2002). [CrossRef]

**54. **C. Draeger and M. Fink, “One-channel time reversal of elastic waves in a chaotic 2d-silicon cavity,” Phys. Rev. Lett. **79**, 407 (1997). [CrossRef]

**55. **G. Montaldo, D. Palacio, M. Tanter, and M. Fink, “Time reversal kaleidoscope: A smart transducer for three-dimensional ultrasonic imaging,” Appl. Phys. Lett. **84**, 3879–3881 (2004). [CrossRef]

**56. **O. Yurduseven, T. Fromenteze, D. L. Marks, J. N. Gollub, and D. R. Smith, “Frequency-diverse computational microwave phaseless imaging,” IEEE Antennas Wireless Propag. Lett (2017). [CrossRef]

**57. **J. J. van Zyl, M. Arii, and Y. Kim, “Model-based decomposition of polarimetric sar covariance matrices constrained for nonnegative eigenvalues,” IEEE Trans. Geosci. Remote Sens. **49**, 3452–3459 (2011). [CrossRef]