Journal of Electrical and Computer Engineering

Volume 2013 (2013), Article ID 636972, 11 pages

http://dx.doi.org/10.1155/2013/636972

Research Article

Performance of 2D Compressive Sensing on Wide-Beam Through-the-Wall Imaging

Edison Cristofani, Mathias Becquaert, and Marijke Vandewal

CISS Department, Royal Military Academy, 30 Avenue de la Renaissance, 1000 Brussels, Belgium

Received 27 May 2013; Revised 18 July 2013; Accepted 24 July 2013

Academic Editor: Sandra Costanzo

Copyright © 2013 Edison Cristofani 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.

Abstract.

Compressive sensing has become an accepted and powerful alternative to conventional data sampling schemes. Hardware simplicity, data, and measurement time reduction and simplified imagery are some of its most attractive strengths. This work aims at exploring the possibilities of using sparse vector recovery theory for actual engineering and defense- and security-oriented applications. Conventional through-the-wall imaging using a synthetic aperture configuration can also take advantage of compressive sensing by reducing data acquisition rates and omitting certain azimuth scanning positions. An ultra-wideband stepped frequency system carrying wide beam antennas performs through-the-wall imaging of a real scene, including a hollow concrete block wall and a corner reflector behind it. Random downsampling rates lower than those announced by Nyquist’s theorem both in the fast-time and azimuth domains are studied, as well as downsampling limitations for accurate imaging. Separate dictionaries are considered and modeled depending on the objects to be reconstructed: walls or point targets. Results show that an easy interpretation of through-the-wall scenes using the -norm and orthogonal matching pursuit algorithms is possible thanks to the simplification of the reconstructed scene, for which only as low as 25% of the conventional SAR data are needed.

1. Introduction

Previous research studies have explored numerous applications of through-the-wall imaging (TWI) in domains which can exploit the convenient penetration capabilities of the 1–5 GHz frequency band in most nonmetal building construction materials [1, 2]. Indeed, being able to perform nondestructive inspection of objects behind a solid wall is a noteworthy application. The literature on TWI widely copes with object or human detection—either moving [3–5] or static [6, 7]. Among others, military, law enforcement, and security operations have already taken profit of the satisfactory performance of TWI. Although most of the studies suggest using an active TWI approach wherever needed, some authors opt for adequately processing opportunistic signals in the surroundings of the scene such as widespread GSM [8] or wireless networks [9]. Despite presenting interesting wall penetration properties, the noncooperative—thus passive—nature of these approaches cannot assure complete reliability and availability, and consequently, an active approach is preferred.

Synthetic Aperture Radar (SAR) allows competitive TWI using rather common antennas with wide beams and simple out-of-focus data acquisition. Thanks to coherent signal integration along a known measurement path, these data are focused on postprocessing, yielding an improved lateral resolution which remains constant at all ranges [10]. Given a measurement path, sufficient data both in the azimuth and fast-time directions must be obtained in order to perform SAR focusing and to satisfy the Nyquist-Shannon anti-aliasing criterion. This, generally, implies using complex hardware which generates amounts of data which can become unbearable when fine resolution imagery is needed.

Compressive Sensing (CS) introduces the possibility of reducing the amount of sensed signals for further digital reconstruction and, likewise, to perform accurate TWI. First used to solve sparse vector recovery problems within mathematical applications, CS has been recently used in a growing list of fields for sparse signal reconstruction, including TWI applications [11–14]. CS is able to reconstruct signals with high probability provided random sampling and sufficient sparsity levels are satisfied [15]. Since CS announces a convex function optimization problem, solutions can be found by using the -norm minimization or, in a less demanding approach, the Orthogonal Matching Pursuit (OMP) greedy algorithm, among an endless list of possibilities.

Actual TWI experiments were conducted for this study. Hollow concrete stretcher blocks were used to build up a solid wall, replicating load-bearing foundation walls used in most permanent constructions. A corner reflector, serving as an object of interest, was placed behind the wall. Scene reconstruction using an ultra-wideband moving sensor was achieved, both in conventional SAR and CS approaches. This study will explore the possibilities of random sampling, with a focus on minimizing the amount of data required to produce accurate enough reconstructions of simple TWI scenes. This action will lead to a significant simplification in hardware, as well as reducing processing complexity and data handling. The remaining of this document is divided into four further sections: Section 2 describes the used experimental setup, as well as signal processing concepts that will be used throughout the work. Section 3 formulates the CS problem to sparse data recovery, and strategies for creating dictionaries are described. Section 4 presents conventional and CS results obtained from actual TWI measurements. Finally, Section 5 details conclusions drawn and future work.

2. Conventional TWI Measurements

2.1. TWI and the Measurement Site

The sensor used for the presented TWI comprises two Schwarzbeck BBHA 9120 D horn antennas [16] (for transmitting and receiving, providing gain stability in the 1–5 GHz band) on a stepped frequency (SF) SAR moving scanner, which translates for 1500 mm following an equispaced azimuth sampling scheme. The sensor used is an Anritsu 37169A Vector Network Analyzer (VNA). A concrete wall was built up for this study in order to obtain results comparable to those of a real TWI situation (see Figure 1). The wall was located at 950 mm from the antennas, in the far-field region for the entire sensor’s frequency band, and a corner reflector was placed at 1240 mm behind the wall.

Figure 1: Schematic representation of the proposed experimental setup.

Hollow concrete stretcher blocks were used (see Figure 2 for shape and dimensions), commonly found in masonry and especially in load-bearing foundation walls.

Figure 2: Dimensions of a hollow concrete stretcher block are 385 mm × 185 mm × 135 mm (length, height, and width, resp.). The outer structure is 35 mm thick.

Range estimation is greatly affected by the change of relative permittivity occurring in the air-material interfaces. Consequently, imaged targets and walls will appear further away and wider than in reality. Given a material with relative permittivity , any transmitted electromagnetic wave through it will propagate at a speed , being the speed of light in vacuum. For an air-concrete-air interface of 35 mm, and considering approximately [17], the second interface would be located at  mm, or  mm further away. As depicted in Figure 3, any object behind a concrete wall will suffer from this phenomenon, whereas its actual location is closer than shown by TWI. Waves traveling through hollow concrete blocks, as those used in the proposed experiments, will encounter two 35-mm air-concrete-air transitions, and therefore the expected range of the corner reflector would be shifted to ca.  mm.

Figure 3: The ranging information is altered due to the different transmission media in the wall. Objects will appear further away.

2.2. The SF Sensor and Signal Model

The SAR scanner operated in this study is based on an SF ultra-wideband sensor in the 1–5 GHz band. Using such SF sensor guarantees low hardware requirements, since the instantaneous bandwidth is kept narrow, while achieving a remarkable range resolution and power output. Instead of transmitting a linear frequency modulated chirp using a continuous or pulsed approach, the SF sensor produces signals by sequentially generating carrier frequencies covering an effective system bandwidth . These frequencies pulses are described by

where is the th frequency generated, is the starting frequency (i.e., 1 GHz for this setup), and is the frequency step. Figure 4 depicts how the transmitted signals are generated.

Figure 4: Waveform generation scheme of an SF sensor.

For any SF systems, the maximum unambiguous range and range resolution parameters are respectively obtained as follows:

In Table 1, the most important parameters of the available experimental setup are shown. During the waveform generation phase, each frequency pulse is transmitted for a given time . Averaging is performed in detection in order to reduce noise and spurious contributions to the signal of interest.

Table 1: Descriptive parameters of the TWI experiment.

VNA model
Anritsu 37169A

Start frequency
1 GHz

End frequency
5 GHz

101

Averaging
10

Antenna opening angle
14°

3750 mm

40 MHz

4 GHz

Tx power
+10 dBm

Tx/RX antennas
Schwarzbeck BBHA 9120 D

Typ. ant. gain
11 dBi

Scene size
1500 mm 2500 mm

37.5 mm

The signal model of the transmitted waveform   for a given can be described as

where corresponds to the fast-time dimension. The sensor will actually transmit a series of individual frequencies which will cover the whole bandwidth. The detected signal model can be described as follows:

where corresponds to the delay caused by the round-trip propagation time of a given th scatterer in the scene located at a distance from the sensor. Distance is simply obtained as the sensor-to-scatterer Euclidean distance. The VNA performs a homogeneous demodulation of the detected signal, which is known as beat signal and is modeled as

2.3. SAR Imaging

SAR measurements use relatively simple antennas with a small aperture length and, therefore, poor azimuth resolution. By translating the sensor as shown in Figure 5 appropriately and coherently combining received signals from the different azimuth positions, the equivalent physical length of the aperture can be increased. An antenna of a theoretically unlimited aperture length could be synthesized yielding very fine azimuth resolutions. Building such antenna would not be feasible in practice; therefore, SAR processing results in a useful tool.

Figure 5: Experimental scene description. The SAR sensor moves in azimuth illuminating the scene with a wide-beam antenna.

Collected unfocused SAR data must be combined coherently using an SAR focusing algorithm. These algorithms perform range/azimuth compression in order to compact received energy, whether in the time or frequency domains, and are mostly based on the matched filter [18–20] assuming the prior knowledge of the transmitted signal model. Figure 6 shows a simplified representation of a range-compressed SAR measurement based on the experimental setup described in Section 2.1.

Figure 6: Expected results after an SAR measurement. Energy is spread in azimuth following a hyperbolic function. Certain elements will be shifted in range (dashed lines).

The distance between the sensor’s positions along the azimuth measurement track and any point in the scene is computed using the Euclidean distance and will follow a hyperbolic function. The collected energy from any given point target will present such hyperbolic shape and an amplitude variation in function of the antenna radiation diagram (both transmitting and receiving antennas).

Produced SAR images achieve interesting improved azimuth resolutions, which can be theoretically computed as

where is the carrier or center frequency of the sensor and is the opening angle of the antenna. According to the used sensor described in Table 1, expected resolutions range and azimuth are and , respectively.

3. Compressive Sensing

3.1. Conventional Sensing

In signal acquisition and processing, the Nyquist-Shannon theorem [21] specifies that any signal can be perfectly reconstructed only after being uniformly sampled at, at least, twice the maximum frequency or bandwidth describing it. Consider

where and are the sampling and maximum frequencies, respectively. This condition can be a remarkable burden for certain systems, where challenging amounts of data are needed when large areas must be measured, high bandwidth is needed, or a combination of both. Handling, storing, or transmitting large amounts of data are onerous tasks, and SAR imaging systems are affected. Complexity of analog-to-digital converters would increase if fine-resolution SAR imagery was needed.

3.2. Problem Formulation

Data recovery of sparse vectors is extensively described in the literature [15, 22, 23]. Let be an arbitrary -dimensional vector or complete measurement, and let be an incomplete observation of , being . Vectors and are related to the linear expression

where is a measurement matrix containing a plausible observation scheme. Illustratively, the th element of would be stored in the th element of if . A perfect recovery of vector can be obtained from if two conditions are fulfilled: vector must be sufficiently sparse and matrix must satisfy the Restricted Isometry Property (RIP) [15, 22]. Vector   is said to be -sparse if the -norm of satisfies that is, the maximum number of nonzero elements describing must be . Nevertheless, having such a vector satisfying the -sparsity condition or demonstrating the RIP for all the possible combinations can be cumbersome tasks, if not implausible. This latter requirement can be relaxed by defining matrix following a random distribution which satisfies the RIP with high probability [24].

In radar imaging, vector is defined as the full observation of the returns from an illuminated scene. Provided that the signal model is known (refer to Section 2.2), vector can be expressed as

where for all positions in the scene, describes the reflectivity coefficients, being , whereas matrix , an over complete dictionary, contains fully sampled received signals as described by the theoretical signal model. Therefore, is populated column wise as

Let then be a matrix containing either zeros or ones, in perfect analogy to matrix , describing a certain sampling scheme. Let a sensor obey such scheme; the sensed signal for a given is now defined as

For the application proposed in this study, only the scene reflectivity coefficients in are needed and not the original measurement vector . The recovery problem to be solved can then be formulated as

being the -norm of vector . However, the -norm is an NP-hard problem and solutions may not exist. Since matrix   sufficiently satisfies the RIP condition, the -norm can be used as a viable alternative to noise-free convex optimization problems

Optimization of noisy measurements must be approximated since a perfect reconstruction will not be possible [22]. A tolerance threshold related to the data noise level or accepted reconstruction error must be imposed as follows:

Several categories different from convex optimization or basis pursuit can be found, each one offering multiple algorithms and variations. One of the most popular in terms of simplicity and efficiency is Orthogonal Matching Pursuit (OMP) [25] and its variations [26]. Being a greedy and simple algorithm (see Algorithm 1), OMP requires a prior knowledge or estimation of the sparsity in order to recover a sparse vector. Assessing the sparsity of a scene is not always possible especially in TWI. A relatively performing solution is implementing a stop condition to avoid the algorithm iterating more than what required by the scene. -norm and OPM will be used in Section 4.

Inputs:

        (i) measured incomplete signal and matrix , populated with simulated wall or point-target model data

         (ii) maximum number of iterations

Output: reconstructed scene vector

() Initialize residual to the measured signal vector and secondary matrix to zero

() Increase . Find the column-vector of that provides highest correlation with the residual

() Add column of to and zero from

() Solve

      and update residual with

() Go to step 2 until reaches stop condition (

Algorithm 1: OMP algorithm description. Adapted from [26].

3.3. CS Random Sampling and Dictionary Structure

Random downsampling in the fast-time and azimuth dimensions implies a substantial reduction in data collection, in measurement time (less azimuth scans), and sensor complexity (fast-time sampling is less demanding). Random sampling schemes must be cautiously imitated when constructing the CS dictionary with the known signal model. Figure 7 describes how azimuth and fast-time random samplings are generated using the conventional SAR approach as a reference. Sampling intervals and azimuth steps follow adjustable random statistical distributions.

Figure 7: Sample selection scheme on full SAR data simulating a real CS sensor. Fast-time samples in white and azimuth measurements in grey are skipped.

3.4. Modeling CS Dictionaries

CS dictionaries can be modeled according to the type of scene that must be reconstructed. A possible solution for the problem here addressed is considering two regions in the measurement (Figure 8), which also add complexity scalability to CS on TWI scenes as follows.

(i)A region closer to the sensor (wall region) where only walls parallel to the sensors’ motion path are expected to appear. Since walls present very strong azimuth redundancies in SAR measurements, the wall region can be reconstructed forcing coarse azimuth and range resolutions in the final imagery. For every possible range in the scene, the wall dictionary will contain in its columns simulated received SAR data based on a simple infinite wall model (Figure 9(a)).(ii)A second region is defined further away from the sensor, where scatterers other than walls are expected. This region is known as the wall-free region. Sufficient azimuth information is necessary to characterize point targets in the scene, and reducing the number of azimuth positions is not a sound option. For each position of a point target in the scene, received SAR signals are simulated and included in the columns of the point-target dictionary (Figure 9(b)).

Figure 8: Two regions are described and will be treated differently. Walls may appear closer to the sensor. Further away only point targets are expected.

(a)

(b)

Figure 9: Range-compressed images of a simulated stepped frequency SAR system. (a) a scene containing scatterers modeling a simple wall model located at 0.99 m in range, used to populate the wall dictionary. (b) a scene containing one-single point target located at range 1.9 m and azimuth 0.75 m, used to populate the point-target dictionary.

3.5. Random Sampling Limitations

Estimating the number of required random samples for CS reconstruction is a complex task since scene sparsity cannot be easily assessed prior to the measurement. However, an approximate expertise can be gained by systematically testing out several random downsampling possibilities. As discussed in Section 3.3, the main restriction will be established by the wall-dictionary reconstruction. Figure 10 shows Peak Signal-to-Noise Ratio (PSNR) reconstruction values, in which from 10% to 90% of the azimuth and fast-time samples of a complete SAR TWI measurement were randomly selected. The underlying purpose of the experiment is determining how accurate the image of a wall can be reconstructed if low, medium, or severe downsampling happens and selecting an adequate working point. Analysis of results suggests, as expected, that using severe random downsampling in any dimension impedes an accurate data reconstruction. Using higher values produces an evident improvement in the reconstructions, but those would not yield hardware or computational time reductions. Based on the compromise between accuracy and simplicity, random sampling should range between 20% and 50% in the fast time and between 20% and 50% for azimuth.

Figure 10: PSNR values in dB for different random downsampling possibilities in fast time and azimuth. A time-domain signal model was used to create the dictionary.

Given the presented results, this study will use a random downsampling of 50% of the usual samples for both fast-time and azimuth dimensions in the wall-free region, which implies keeping only 25% of the necessary samples for conventional SAR processing. Additionally, decimation by two of the azimuth measurements can be safely applied for the wall region, yielding 12.5% of the original data.

3.6. Practical Experimental Details

The experimental TWI setup uses a conventional SAR system, and original CS measurements cannot be obtained. A full SAR measurement is performed with the presented SF SAR system in order to obtain sufficient fast-time and azimuth information for further random downsampling. The SAR simulator used to model the expected scatterers operates likewise.

4. Results and Discussion

4.1. Conventional SAR TWI Results

The experimental setup depicted in Figures 1 and 5 was measured using the sensor described in Table 1, and conventional time-domain SAR processing was applied. In Figure 11(a), in a range-compressed image of the acquired SAR signals, the wall is clearly visible in the leftmost side due to the strong reflections generated. In the center image, the received energy from corner reflector behind the wall is greatly spread in azimuth following the hyperbolic structure depicted in Figure 6.

(a)

(b)

(c)

Figure 11: (a) A range-compressed SAR image. The wall is clearly seen. (b) The energy of the corner reflector is spread in azimuth at approximate range 2.25 m. (c) A fully processed SAR image where the presence of a corner reflector is clearly spotted at range 2.25 m.

After having taken into account azimuth information, the fully compressed SAR image is reconstructed as shown in Figure 11(c). Indeed, the wall is clearly seen as well as several artifacts in adjacent ranges further from the sensor. Such artifacts correspond to low returns from the second layer of the hollow concrete blocks as well as multiple internal reflections which appear in object-free ranges, as predicted in Section 2.1.

4.2. CS Experimental Results

A CS head and processor are imitated in software following the considerations explained in Sections 2 and 3. Raw SAR data shown in Section 4.1 are sampled at different random sampling rates and reconstructed assuming scene sparsity. In this section, presented results were obtained according to the methodology described in Section 3.4 for efficient wall reconstruction and ranging and behind-the-wall accurate reconstruction. The reconstruction dictionaries are built up using simple wall and point-target models and a stepped-frequency SAR simulator.

Figure 12(a) describes how the signal simulator is fed with coordinates produced by the models, which follows to a CS sampler simulator which downsamples the simulated data in range and azimuth according to the provided sampling parameters. The resulting downsampled data are vectorized and stored in the columns of the wall or point-target dictionary, or , respectively. Scene reconstruction of actual measurements is depicted in Figure 12(b). A CS head is partially simulated in software, downsampling the measured signals following a given random sampling scheme identical to that used during the dictionary generation. The CS reconstruction block is fed with the incomplete measurements and the previously generated dictionaries and operational parameters for reconstruction, which may vary depending on the reconstruction algorithms (OMP or the -norm).

(a)

(b)

Figure 12: Block diagrams of the (a) wall and point-target dictionary generation and (b) scene reconstruction from actual incomplete measurements.

Figure 13 shows the reconstruction of the wall reflectivity coefficients vector using OMP and a wall-only dictionary coarsely modeled for efficient processing purposes. The range of maximum reflectivity, 0.95 m, corresponds to the imaged wall, which satisfies expected results.

Figure 13: CS reconstruction of the scene using a coarse wall dictionary for wall ranging. The actual wall was located at 0.95 m.

A more accurate CS dictionary is then used for point-target reflectivity coefficients reconstruction, or vector . Based on , it is assumed that one-single wall was present in the scene. A guard range is kept to avoid the clutter region in ranges immediately after the wall (refer to Figures 3, 6, and 8). Figures 14–16 present reconstructed scenes, where a corner reflector is expected to appear at range 2.25 m, for several range downsampling options (60%, 50%, and 40%). OMP generates simplified but correct reconstructions of two or three orders of magnitude faster, on average, than the -norm algorithm.

(a)

(b)

Figure 14: CS reconstructions of the behind-the-wall scene using a point-target dictionary. Downsampling factors were set to 60% for fast-time and 50% for azimuth. (a) -norm; (b) OMP.

(a)

(b)

Figure 15: CS reconstructions of the behind-the-wall scene using a point-target dictionary. Downsampling factors were set to 50% for fast-time and 50% for azimuth. (a) -norm; (b) OMP.

(a)

(b)

Figure 16: CS reconstructions of the behind-the-wall scene using a point-target dictionary. Downsampling factors were set to 40% for fast-time and 50% for azimuth. (a) -norm; (b) OMP.

Figure 14 depicts results using downsampling factors of 50% and 60% for azimuth and fast time, respectively. Accurate ranging of the corner reflector is obtained (2.25 m) and comparable to that of the conventional SAR imagery shown in Figure 10 right. Results generated by the -norm clearly indicate where the corner reflector is located, while OPM extracts the most significant reflectivity coefficients. Given the theoretical range resolution of 37.5 mm, the two targets extracted using OMP can be assumed to be in the same resolution cell. As the number of fast-time samples kept is reduced (Figures 15 and 16), inaccuracies in the corner reflector’s positioning and false alarms appear, although a general agreement can be reached on the presence of a strong reflector in the scene.

5. Conclusions and Future Work

In this work, CS has been successfully applied to TWI of a real concrete wall. SAR measurements were preprocessed in software to simulate a CS sensor which would perform random sampling in the fast-time and azimuth dimensions. The particular considerations regarding the application of CS onto SAR measurements for TWI were addressed. By using CS, any kind of SAR processing was avoided. A conceptually simple procedure was designed to automatically distinguish between wall and wall-free regions, which exploits the possibilities that CS can offer and reduces computational complexity and avoids having to treat wall clutter. A systematic investigation was performed to determine the minimum random downsampling with respect to conventional SAR TWI without compromising reconstruction accuracy. Data volume was reduced to as low as 25% of the usual required samples, enabling measurement time and hardware complexity reduction.

Future work is focused on finding alternatives to further exploit the CS capabilities for TWI. Some of the lines of investigation the authors are following related to signal processing and CS are as follows.

(i)Reducing the effects of strong reflections from the wall introduces the reconstruction of point-target regions and improves the accuracy of detection. The authors currently are working on strategies to overcome these negative effects.(ii)Multiple echoes and patterns due to the wall structure appear in further adjacent ranges. Using a complete SAR model of a wall will help identify which scatterers can be considered as wall clutter or targets of interest behind the wall.(iii)Range inaccuracies due to inhomogeneous transmission media will be corrected to provide a correct range estimation. A complete wall model will be developed based on existing publications.(iv)The SF system used in the experiments as well as the antennas has a frequency-dependent performance. A complete SAR signal model must be developed to create CS dictionaries as accurate as possible. Moreover, the mutual interaction between both adjacent antennas modifies their radiation patterns. The antennas will be measured in an anechoic chamber to include this effect in the complete signal model. (v)CS applied to 3D SAR imaging proposes challenging open questions regarding sparsity assessment and strategies to guarantee the viability the of proposed methodologies in this work.(vi)Data transformation bases will be studied to enhance data sparsity and facilitate accurate scene reconstructions with even less random samples and simpler hardware.(vii)Other efficient sparse data reconstruction algorithms will be explored and modified to the TWI needs where applicable.

Acknowledgments

The authors of this study would like to show their gratitude to the OMRA and LEMA laboratories of the Royal Military Academy, Belgium, for providing both technical support and logistic support during construction and measurement of the presented experimental scenes.

References

O. Landron, M. J. Feuerstein, and T. S. Rappaport, “A comparison of theoretical and empirical reflection coefficients for typical exterior wall surfaces in a mobile radio environment,” IEEE Transactions on Antennas and Propagation, vol. 44, no. 3, pp. 341–351, 1996.

A. Muqaibel, A. Safaai-Jazi, A. Bayram, A. M. Attiya, and S. M. Riad, “Ultra-wideband through-the-wall propagation,” IEE Proceedings, pp. 581–588, 2005.

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 Proceedings of the 4th IEEE International Symposium on Phased Array Systems and Technology (Array '10), pp. 551–558, October 2010.

F. Ahmad, M. G. Amin, and P. Setlur, “Through-the-wall target localization using dual-frequency CW radars,” in Proceedings of the Sensors, and Command, Control, Communications, and Intelligence (C3I) Technologies for Homeland Security and Homeland Defense V, vol. 6201 of Proceedings of SPIE, 62010H, May 2006.

Y. Yoon and M. G. Amin, “Spatial filtering for wall-clutter mitigation in through-the-wall radar imaging,” IEEE Transactions on Geoscience and Remote Sensing, vol. 47, no. 9, pp. 3192–3208, 2009.

R. Solimene, F. Soldovieri, G. Prisco, and R. Pierri, “Three-dimensional through-wall imaging under ambiguous wall parameters,” IEEE Transactions on Geoscience and Remote Sensing, vol. 47, no. 5, pp. 1310–1317, 2009.

W. Zhang, L. Li, and F. Li, “Autofocusing imaging through the unknown building walls,” in Proceedings of the Asia Pacific Microwave Conference (APMC '08), pp. 1–5, December 2008.

D. K. P. Tan, M. Lesturgie, H. Sun, W. Li, and Y. Lu, “Gsm based through-the-wall passive radar demonstrator for motion sensing,” in Proceedings of the New Trends for Environmental Monitoring Using Passive Systems (PASSIVE '08), pp. 1–4, Hyères, France, October 2008.

K. Chetty, G. E. Smith, and K. Woodbridge, “Through-the-wall sensing of personnel using passive bistatic wifi radar at standoff distances,” IEEE Transactions on Geoscience and Remote Sensing, vol. 50, no. 4, pp. 1218–1226, 2012.

I. G. Cumming and F. H. Wong, Digital Processing of Synthetic Aperture Radar Data: Algorithms and Implementation, Artech House, Norwood, Mass, USA, 2005.

M. Leigsnering, C. Debes, and A. M. Zoubir, “Compressive sensing in through-the-wall radar imaging,” in Proceedings of the 36th IEEE International Conference on Acoustics, Speech, and Signal Processing (ICASSP '11), pp. 4008–4011, May 2011.

F. Ahmad and M. G. Amin, “Partially sparse reconstruction of behind-the-wall scenes,” in Proceedings of the SPIE Symposium on Defense, Security, and Sensing, Compressive Sensing Conference, Baltimore, Md, USA, April 2012.

Q. Huang, L. Qu, B. Wu, and G. Fang, “UWB through-wall imaging based on compressive sensing,” IEEE Transactions on Geoscience and Remote Sensing, vol. 48, no. 3, pp. 1408–1415, 2010.

E. Lagunas, M. G. Amin, F. Ahmad, and M. Najar, “Compressive sensing for through wall radar imaging of stationary scenes using arbitrary data measurements,” in Proceedings of the 11th International Conference on Information Science, Signal Processing and their Applications (ISSPA '12), pp. 1347–1352, July 2012.

D. L. Donoho, “Compressed sensing,” IEEE Transactions on Information Theory, vol. 52, no. 4, pp. 1289–1306, 2006.

(n. d.), Schwarzbeck Mess—Elektronik, Double Ridged Broadband Horn BBHA, 9120 D, http://www.schwarzbeck.de/Datenblatt/k9120d.pdf

C. Thajudeen, A. Hoorfar, F. Ahmad, and T. Dogaru, “Measured complex permittivity of walls with different hydration levels and the effect on power estimation of twri target returns,” Progress in Electromagnetics Research B, vol. 30, pp. 177–199, 2011.

R. Bamler, “A comparison of range-Doppler and wavenumber domain SAR focusing algorithms,” IEEE Transactions on Geoscience and Remote Sensing, vol. 30, no. 4, pp. 706–713, 1992.

I. G. Cumming, Y. L. Neo, and F. H. Wong, “Interpretations of the omega-K algorithm and comparisons with other algorithms,” in Proceedings of the IEEE International Geoscience and Remote Sensing Symposium (IGARSS '03), pp. 1455–1458, July 2003.

H. L. Van Trees, “Detection of signals—estimation of signal parameters,” in Detection, Estimation, and Modulation Theory: Detection, Estimation, and Linear Modulation Theory, John Wiley & Sons, New York, NY, USA, 2002.

C. E. Shannon, “A mathematical theory of communication,” The Bell Systems Technical Journal, vol. 27, pp. 379–423, 623–656, 1948.

E. J. Candès, J. K. Romberg, and T. Tao, “Stable signal recovery from incomplete and inaccurate measurements,” Communications on Pure and Applied Mathematics, vol. 59, no. 8, pp. 1207–1223, 2006.

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.

H. Rauhut, “Compressive sensing and structured random matrices,” in Theoretical Foundations and Numerical Methods for Sparse Recovery, M. Fornasier, Ed., vol. 9 of Radon Series on Computational and Applied Mathematics, pp. 1–92, deGruyter, Berlin, Germany, 2010.

T. Blumensath and M. E. Davies, “In greedy pursuit of new directions: (nearly) orthogonal matching pursuit by directional optimisation,” in Proceedings of the 15th European Signal Processing Conference (EUSIPCO '07), Poznan, Poland, September 2007.

J. A. Tropp and D. Needell, “CoSaMP: iterative signal recovery from incomplete and inaccurate samples,” presented at CoRR, 2008