Progress in Microwave Imaging: From Theoretical Developments to CuttingEdge Applications
View this Special IssueResearch Article  Open Access
Cylindrical ThreeDimensional MillimeterWave Imaging via Compressive Sensing
Abstract
Millimeterwave (MMW) imaging techniques have been used for the detection of concealed weapons and contraband carried by personnel. However, the future application of the new technology may be limited by its large number of antennas. In order to reduce the complexity of the hardware, a novel MMW imaging method based on compressive sensing (CS) is proposed in this paper. The MMW images can be reconstructed from the significantly undersampled backscattered data via the CS approach. Thus the number of antennas and the cost of system can be further reduced than those based on the traditional imaging methods that obey the Nyquist sampling theorem. The effectiveness of the proposed method is validated by numerical simulations as well as by real measured data of objects.
1. Introduction
Millimeterwave (MMW) imaging techniques have been applied for the detection of concealed weapons and contraband carried by personnel at airports and other secure locations [1–5]. Millimeter waves can penetrate common clothing barriers to form an image of a person as well as any concealed items with different reflectivity or emissivity. Also, millimeter waves are nonionizing and, therefore, pose no known health hazard at moderate power levels, not like Xray systems.
However, a large number of antennas and switches (maybe from hundreds to ten thousands) are needed to construct one or twodimensional (1D or 2D) array to obtain threedimensional (3D) images, which increases the cost of the imaging system and then limits its wide application. In order to reduce the complexity and cost of the system, we propose a compressive sensing (CS) based method for 3D MMW imaging. The area of CS was initiated in 2006 by Candès et al. [6] and by Donoho [7].
CS has been widely explored in the domain of radar imaging [8–11], in which the image was usually converted into a vector. And the 1D optimization algorithm was used to recover the image. However, if we convert the 3D image into a vector, a huge sensing matrix will be introduced. Then the computational complexity will be largely increased. In this paper, we introduce an operator as the sensing matrix to represent the traditional imaging process based on wavenumber domain algorithm (WDA) [12]. Thus, largescale matrix computations are avoided. The image resolution of the WDA is determined by signal bandwidth, signal frequency, and the length of synthetic aperture. And the resolution of the CSbased imaging algorithm is determined mainly by the signaltonoise ratio and the matching degree of signal model to the requirements of CS, such as the restricted isometry property (RIP). The simulation and experiment results show that the CSbased method offers a much more accurate image than the traditional WDA just using a small subset of the full samples.
The rest of the paper is organized as follows. In Section 2, the CS imaging method based on the traditional WDA is introduced. The results demonstrating the efficiency and accuracy of the method are shown in Section 3. Section 4 summarizes the conclusions.
2. CSBased Method for 3D MMW Imaging
For the 3D imaging geometry is illustrated in Figure 1. An object formed by a distribution of independent scatterers of reflectivity is located at the center of the imaging scene. The backscattered field after demodulation is given by where is the spatial frequency, represents the light speed, is the distance from each antenna element at to each scattering center at , and is the radius of the cylindrical aperture, as shown in Figure 1.
The scheme of the traditional WDA [12] for nearfield cylindrical 3D imaging is illustrated in Figure 2. We represent the imaging process as an operator. Then, the inverse process can be represented as the inversion of the operator. Thus the imaging process is given by where the operator denotes the imaging process as shown in Figure 2; that is, where represents 2D interpolation.
Then can be expressed as the inversion of (2) as follows: where denotes the inverse process of (3).
Due to the sparsity of the image scene , only a small subset of scattered data is enough to be used to reconstruct the image exactly by solving the optimization problem [13]: where is simply the number of nonzero components of , and represents a sparse transform, making a sparse vector.
However, solving (5) is an NPhard problem and thus practically not feasible. Instead, its convex relaxation is usually considered as where and denote the norm and norm, respectively. In this paper, we use the totalvariation (TV) as a sparse transform. The parameter is used to balance the twin objectives of minimizing both error and sparsity.
The operators and correspond to the Fourier matrices if we express the signal model directly. Thus, and satisfy the RIP required by CS, and both (5) and (6) have the same unique solution [6]. There are many stateoftheart optimization algorithms designed to solve the convex optimization problem of (6) [14–17]. Here we use a nonlinear conjugate gradient (CG) descent algorithm described in [18]. The computational complexities of and both are an order of , , and are the sizes of the scattered data, being the interpolator length. And the computational complexity of the key part of the CG algorithm during an iteration is an order of [18].
3. Results
The aforementioned algorithm is tested with numerical simulations and experimental data. The simulation parameters are shown in Table 1. Figure 3 shows the target model with 8 pointlike scatterers with equal radar cross section. The 3D imaging result (isosurface of −13 dB relative to the maximum value) of WDA using the full sampling data and the result of the proposed CSbased method using only 25% of the full data are shown in Figures 4(a) and 4(b), respectively. Clearly, even with a small subset of the full data, the CSbased method can provide a more accurate imaging result with higher resolution than WDA.

(a)
(b)
The experiments on real measured data are performed. A steppedfrequency Wband radar with two circular polarized horn antennas is used for the measurement, which is shown in Figure 5. Three screws and a card are measured, respectively, as shown in Figure 6. The diameter of the cross section of the screw is about 1 cm.
The corresponding MMW imaging results of WDA using the full measured data are shown in Figures 7(a) and 8(a), respectively. And the results of the CSbased method using only 25% of the full data are demonstrated in Figures 7(b) and 8(b), respectively. Due to the fact that the card is a distributed target, the imaging result in Figure 8(b) is not so good as the pointlike target in Figure 7(b). However, the main characteristics remained; for example, the chip inside the card can be easily seen in Figure 8(b).
(a)
(b)
(a)
(b)
4. Conclusion
A CSbased method has been proposed for cylindrical nearfield 3D MMW imaging. The imaging process and its inversion are represented by an operator and its inversion, respectively, to avoid largescale matrix computation. Thus the CS approach can be applied to the cylindrical 3D imaging by using the operator as the sensing matrix. Simulation and experiment results show that the CSbased method can offer highly accurate and reliable results, even using only a small subset of the full samples compared with the traditional wavenumber domain algorithm.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 D. M. Sheen, D. L. McMakin, and T. E. Hall, “Threedimensional millimeterwave imaging for concealed weapon detection,” IEEE Transactions on Microwave Theory and Techniques, vol. 49, no. 9, pp. 1581–1592, 2001. View at: Publisher Site  Google Scholar
 R. Appleby and R. N. Anderton, “Millimeterwave and submillimeterwave imaging for security and surveillance,” Proceedings of the IEEE, vol. 95, no. 8, pp. 1683–1690, 2007. View at: Publisher Site  Google Scholar
 D. M. Sheen, D. McMakin, and T. E. Hall, “Nearfield threedimensional radar imaging techniques and applications,” Applied Optics, vol. 49, no. 19, pp. E83–E93, 2010. View at: Publisher Site  Google Scholar
 S. S. Ahmed, A. Genghammer, A. Schiessl, and L.P. Schmidt, “Fully electronic Eband personnel imager of 2 m2 aperture based on a multistatic architecture,” IEEE Transactions on Microwave Theory and Techniques, vol. 61, no. 1, pp. 651–657, 2013. View at: Publisher Site  Google Scholar
 D. M. Sheen and T. E. Hall, “Reconstruction techniques for sparse multistatic linear array microwave imaging,” in Passive and Active MillimeterWave Imaging XVII, vol. 9078 of Proceedings of SPIE, p. 12, June 2014. View at: Publisher Site  Google Scholar
 E. J. Candès, J. Romberg, and T. Tao, “Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information,” IEEE Transactions on Information Theory, vol. 52, no. 2, pp. 489–509, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 D. L. Donoho, “Compressed sensing,” IEEE Transactions on Information Theory, vol. 52, no. 4, pp. 1289–1306, 2006. View at: Publisher Site  Google Scholar  MathSciNet
 R. G. Baraniuk and P. Steeghs, “Compressive radar imaging,” in Proceedings of the IEEE Radar Conference, pp. 128–133, April 2007. View at: Publisher Site  Google Scholar
 L. C. Potter, E. Ertin, J. T. Parker, and M. Çetin, “Sparsity and compressed sensing in radar imaging,” Proceedings of the IEEE, vol. 98, no. 6, pp. 1006–1020, 2010. View at: Publisher Site  Google Scholar
 X. Zhang, P. Huang, and X. Li, “2D SAR imaging scheme based on compressive sensing,” Electronics Letters, vol. 50, no. 2, pp. 114–116, 2014. View at: Publisher Site  Google Scholar
 W. J. Zhang and A. Hoorfar, “A generalized approach for SAR and MIMO radar imaging of building interior targets with compressive sensing,” IEEE Antennas and Wireless Propagation Letters, vol. 14, pp. 1052–1055, 2015. View at: Publisher Site  Google Scholar
 D. M. Sheen, D. L. McMakin, T. E. Hall, and R. H. Severtsen, “Realtime wideband cylindrical holographic surveillance system,” US Patent 5 859 609, 1999. View at: Google Scholar
 S. Sun, G. Zhu, and T. Jin, “Novel methods to accelerate CS radar imaging by NUFFT,” IEEE Transactions on Geoscience and Remote Sensing, vol. 53, no. 1, pp. 557–566, 2015. View at: Publisher Site  Google Scholar
 Y. Nesterov, Introductory Lectures on Convex Optimization. A Basic Course, Kluwer Academic Publishers, Boston, Mass, USA, 2004. View at: Publisher Site  MathSciNet
 J. A. Tropp and A. C. Gilbert, “Signal recovery from random measurements via orthogonal matching pursuit,” IEEE Transactions on Information Theory, vol. 53, no. 12, pp. 4655–4666, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. A. T. Figueiredo, R. D. Nowak, and S. J. Wright, “Gradient projection for sparse reconstruction: application to compressed sensing and other inverse problems,” IEEE Journal on Selected Topics in Signal Processing, vol. 1, no. 4, pp. 586–597, 2007. View at: Publisher Site  Google Scholar
 A. Beck and M. Teboulle, “A fast iterative shrinkagethresholding algorithm for linear inverse problems,” SIAM Journal on Imaging Sciences, vol. 2, no. 1, pp. 183–202, 2009. View at: Publisher Site  Google Scholar
 M. Lustig, D. Donoho, and J. M. Pauly, “Sparse MRI: the application of compressed sensing for rapid MR imaging,” Magnetic Resonance in Medicine, vol. 58, no. 6, pp. 1182–1195, 2007. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Guoqiang Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.