# On the Use of Passive Reflecting Surfaces and Compressive Sensing Techniques for Detecting Security Threats at Standoff Distances

Authors: Yolanda Rodriguez-Vaqueiro; Jose Angel Martinez-Lorenzo
Journal: International Journal of Antennas and Propagation
 (2014)

---

## Abstract

This work presents a new radar system concept, capable of detecting explosive related threats at standoff distances. The system consists of a two-dimensional aperture of randomly distributed transmitting/receiving antenna elements and a set of passive reflecting surfaces (PRS) positioned in the vicinity of the target. The PRS act as a mirror that enhances the field of view of the radar system, thus increasing its resolution. A 3D imaging algorithm, based on novel compressive sensing techniques, is used in this work. This system configuration provides a resolution of 7.1 mm in cross-range and 25 mm in range, when the target is at 10 m range, and the radar works at 60 GHz center frequency and has 6 GHz bandwidth.


---

## Body

## 1. Introduction

During the last decade, new systems based on Millimeter-Wave-Radar technology have been deployed in airport checkpoints all around the world [1]. Millimeter-wave imaging systems are preferred to X-ray imaging systems [2–4] for this particular application, because the former do not use ionizing radiation. Millimeter-wave imaging systems have proved to be successful in finding explosives concealed underclothing [5–8]; the success of this technology is mainly due to the short range existing between the sensing components of the system and the person under test. A new important challenge arises when this technology is used for detecting security threats at standoff distances [9–12], which include ranges running between ten to fifty meters. For this particular concept of operation, the resolution of the imaging system is reduced when the range is increased.Previous work [13] showed a successful 2D reconstruction at 40 m standoff distance, when a passive array of scatters (PAS) was located between the radar and the person under test. The PAS can be thought of as a millimeter wave lens that enhances the field of view of the radar system, thus increasing its resolution. The same paper also demonstrated that compressive sensing techniques [14–21], which have been successfully used in near-field imaging problems [22, 23], can be used for imaging security threats at standoff distances [13].In this paper, a new type of passive surfaces is used by the imaging system. In particular, the PAS, which works in transmission mode, is replaced by a set of passive reflecting surfaces (PRS), which works in reflection mode. The PRS add additional functionality to the imaging system. Additionally, the imaging algorithm presented in [13] used a single frequency to create 2D images; in this paper, multiple frequencies are used to generate 3D images. Finally, this paper improves the accuracy of the electromagnetic modeling of the passive surfaces by using a physical optics algorithm instead of the ideal point scatter algorithm used in [13].This paper is structured as follows. Section2 presents the configuration of the system, Section 3 describes the mathematical formulation for the imaging problem, and Section 4 shows a numerical example. Finally, the last section presents the conclusions.


## 2. System Configuration

### 2.1. System Concept of Operation

The proposed system configuration is shown schematically in Figure1. It is composed of an inexpensive, high-resolution radar system that can still fit in or on a van. The field of view of the radar system is increased by using a PRS, which can be thought of as a millimeter wave mirror that enhances the resolution of the imaging system. The PRS are a passive structure that backscatters the electromagnetic waves impinging on its surface; this concept of operation is different than that of the passive array of scatter (PAS) used in [13], which works in transmission mode, acting as a millimeter wave lens, rather than in reflection mode.Figure 1
General sketch of our van-based, high resolution radar system for standoff detection of potential suicide bombers. The circles in the vicinity of the person under test represent the passive reflecting surfaces used to enhance the field of view and resolution of the imaging system.


### 2.2.  System Parameters

Figure2 represents a top view of the configuration and the parameters of the system. The blue dots, on the left, represent the positions of the transmitting and receiving antennas. The radar is located on a square aperture of width L1, and the total number of transmitting/receiving antennas is na. The orange dots represent the positions of the reflecting elements composing the PRS. The PRS is also located on a square aperture of width L2, and the total number of elements on the PRS is ns. The person under test is represented by the red silhouette on the top right, and the image reconstruction is performed in np pixels located on multiple two-dimensional planes, represented by green lines in Figure 2. The distance between the geometrical center of the radar and the geometrical center of the PAS is z0, and the distance between the latter point and the person under test is x0. The cross-range resolution l of the imaging system can be approximated by the following equation:
(1)l=λL2cos⁡αx0,
where λ is the wavelength at the center frequency of the imaging system, and α is the angle existing between the normal vector of the plane that contains the PRS and the vector that goes from the center of the PRS to the center of the imaging region. The range resolution lr of the imaging system is given by the following equation [9]:
(2)lr=c2BW,
where c is the propagation speed of electromagnetic waves in free space, and BW is the bandwidth of the imaging system. The system uses nf frequencies.Figure 2
Top view of the radar configuration. The blue circles on the left represent a thinned array of transmitter/receiver antennas; the orange dots on the right represent the passive reflecting surfaces, which randomly redirect the energy of the radar towards the target; the person under test (target) is represented by the red silhouette on the top, and the two dimensional planes over which the reconstruction is implemented are represented by the green lines.


## 2.1. System Concept of Operation

The proposed system configuration is shown schematically in Figure1. It is composed of an inexpensive, high-resolution radar system that can still fit in or on a van. The field of view of the radar system is increased by using a PRS, which can be thought of as a millimeter wave mirror that enhances the resolution of the imaging system. The PRS are a passive structure that backscatters the electromagnetic waves impinging on its surface; this concept of operation is different than that of the passive array of scatter (PAS) used in [13], which works in transmission mode, acting as a millimeter wave lens, rather than in reflection mode.Figure 1
General sketch of our van-based, high resolution radar system for standoff detection of potential suicide bombers. The circles in the vicinity of the person under test represent the passive reflecting surfaces used to enhance the field of view and resolution of the imaging system.


## 2.2.  System Parameters

Figure2 represents a top view of the configuration and the parameters of the system. The blue dots, on the left, represent the positions of the transmitting and receiving antennas. The radar is located on a square aperture of width L1, and the total number of transmitting/receiving antennas is na. The orange dots represent the positions of the reflecting elements composing the PRS. The PRS is also located on a square aperture of width L2, and the total number of elements on the PRS is ns. The person under test is represented by the red silhouette on the top right, and the image reconstruction is performed in np pixels located on multiple two-dimensional planes, represented by green lines in Figure 2. The distance between the geometrical center of the radar and the geometrical center of the PAS is z0, and the distance between the latter point and the person under test is x0. The cross-range resolution l of the imaging system can be approximated by the following equation:
(1)l=λL2cos⁡αx0,
where λ is the wavelength at the center frequency of the imaging system, and α is the angle existing between the normal vector of the plane that contains the PRS and the vector that goes from the center of the PRS to the center of the imaging region. The range resolution lr of the imaging system is given by the following equation [9]:
(2)lr=c2BW,
where c is the propagation speed of electromagnetic waves in free space, and BW is the bandwidth of the imaging system. The system uses nf frequencies.Figure 2
Top view of the radar configuration. The blue circles on the left represent a thinned array of transmitter/receiver antennas; the orange dots on the right represent the passive reflecting surfaces, which randomly redirect the energy of the radar towards the target; the person under test (target) is represented by the red silhouette on the top, and the two dimensional planes over which the reconstruction is implemented are represented by the green lines.


## 3. Mathematical Formulation for the Imaging Problem

### 3.1. Sensing Matrix

The system works on a multiple monostatic configuration, in which each element of the array transmits and receives in different slots of time without interacting with the radiation of other elements in the array.Under this configuration, the sensing matrixA∈Cnt×np establishes a linear relationship between the unknown complex vector x∈Cnp, with np the number of unknowns, and the measured complex field data y∈Cnt, with nt=na·nf, the number of antennas times the number of frequencies. This relationship can be expressed in a matrix form as follows:
(3)A·x+n=y,
where n∈Cnt represents the noise collected by each receiving antenna for a given frequency.Each elementaij of the sensing matrix is computed as follows: (1) the field radiated by the ith antenna, modeled as an infinitesimal aperture [24], is used to illuminate the PRS, for i=1,2,…,nt; (2) the scattered field produced by the PRS on the imaging region is computed by a physical optics algorithm [25]; (3) the field in the imaging region is backscattered by an infinitesimal sphere located at the center of the jth pixel, for j=1,2,…,np; (4) this backscattered field is used to illuminate the PRS again; and (5) the backscattered field produced by the PRS on the ith antenna is stored as the coefficient of the sensing matrix aij. This first order approximation neglects the mutual interaction amongst scatters located in the imaging region.The PRS are made of metal, and they have fixed reflection coefficient equal to −1. Having PRS made of different materials, therefore having different reflection coefficient, does not affect to the quality of the reconstructed image if this reflection coefficient is taken into consideration by the sensing matrix.


### 3.2.  Imaging Algorithm Using Compressive Sensing Approach

The proposed radar system is designed in accordance with the compressive sensing theory [14–21]. In order to apply such principles for standoff detection of explosive related threats, certain mathematical conditions must be satisfied by the sensing matrix A and the reconstructed reflectivity image x. These conditions can be summarized as follows [18]: (1) the sensing matrix must satisfy the restricted-isometry-property condition, which is related to the independence of the columns of the matrix and (2) the unknown reflectivity vector must accept a sparse representation as a solution, which is related to the number of nonzero entries in the solution vector. The parameters of the system can be modified until these two conditions are satisfied; the optimized parameters include the following: aperture length of the radar, aperture length of the PRS, resolution on the reconstruction plane, number of antennas on the radar aperture, number of scatters of the PRS, working frequency, separation between the radar and the PRS, and the separation between the PRS and the target. In this work, this optimization is done manually, but it is expected that in further research contributions such optimization process will be automated.If the two aforementioned conditions are satisfied, then the reconstruction of the unknown vector can be performed with a small number of measurements (transmitting/receiving antennas) by solving the following convex problem [20]:
(4)min⁡∥x∥1,s.t.ii∥Ax-y∥2<δ,
where ∥x∥1 represents the norm-one of the vector x, ∥Ax-y∥2 represents the norm-two of the residual error Ax-y, and δ is an upper bound for the residual error.Additionally, it is possible to perform the reconstruction using CS even when the unknown vector is not sparse by using a transformation in a dictionaryW. The analysis-based problem is expressed as
(5)min⁡∥W*x∥1,s.t.ii∥Ax-y∥2<δ,
where W is a dictionary that allows for a sparse representation of x, and (·)* is the conjugate operator.In this particular problem, the magnitude of the gradient of the unknown reflectivity function, which describes the reflectivity properties of the object under test, has been minimized. This transformation of the unknown reflectivity vectorx[k,m,n] accepts a sparser representation more than the vector x[k,m,n], where k and m are the cross-range indices and n is the range index of a given pixel in the reconstruction domain. For this reason, a total variation (TV) transformation is used to reconstruct the image [20]. The TV problem can be written as
(6)min⁡∥x∥TV,s.t.ii∥Ax-y∥2<δ,
where ∥x∥TV represents the TV norm of the unknown vector x, and it is computed as follows:
(7)∥x∥TV=∑k,m‍∥∇x[k,m,n]∥2,
where ∇ is given by
(8)∇x[k,m,n]=[x[k+1,m,n]-x[k,m,n]x[k,m+1,n]-x[k,m,n]].
It should be noted that ∇ only considers the partial differences in the two cross-range directions. In this way, the reconstruction for each range plane can be performed independently, thus allowing the use of the TV algorithm described in [20].


## 3.1. Sensing Matrix

The system works on a multiple monostatic configuration, in which each element of the array transmits and receives in different slots of time without interacting with the radiation of other elements in the array.Under this configuration, the sensing matrixA∈Cnt×np establishes a linear relationship between the unknown complex vector x∈Cnp, with np the number of unknowns, and the measured complex field data y∈Cnt, with nt=na·nf, the number of antennas times the number of frequencies. This relationship can be expressed in a matrix form as follows:
(3)A·x+n=y,
where n∈Cnt represents the noise collected by each receiving antenna for a given frequency.Each elementaij of the sensing matrix is computed as follows: (1) the field radiated by the ith antenna, modeled as an infinitesimal aperture [24], is used to illuminate the PRS, for i=1,2,…,nt; (2) the scattered field produced by the PRS on the imaging region is computed by a physical optics algorithm [25]; (3) the field in the imaging region is backscattered by an infinitesimal sphere located at the center of the jth pixel, for j=1,2,…,np; (4) this backscattered field is used to illuminate the PRS again; and (5) the backscattered field produced by the PRS on the ith antenna is stored as the coefficient of the sensing matrix aij. This first order approximation neglects the mutual interaction amongst scatters located in the imaging region.The PRS are made of metal, and they have fixed reflection coefficient equal to −1. Having PRS made of different materials, therefore having different reflection coefficient, does not affect to the quality of the reconstructed image if this reflection coefficient is taken into consideration by the sensing matrix.


## 3.2.  Imaging Algorithm Using Compressive Sensing Approach

The proposed radar system is designed in accordance with the compressive sensing theory [14–21]. In order to apply such principles for standoff detection of explosive related threats, certain mathematical conditions must be satisfied by the sensing matrix A and the reconstructed reflectivity image x. These conditions can be summarized as follows [18]: (1) the sensing matrix must satisfy the restricted-isometry-property condition, which is related to the independence of the columns of the matrix and (2) the unknown reflectivity vector must accept a sparse representation as a solution, which is related to the number of nonzero entries in the solution vector. The parameters of the system can be modified until these two conditions are satisfied; the optimized parameters include the following: aperture length of the radar, aperture length of the PRS, resolution on the reconstruction plane, number of antennas on the radar aperture, number of scatters of the PRS, working frequency, separation between the radar and the PRS, and the separation between the PRS and the target. In this work, this optimization is done manually, but it is expected that in further research contributions such optimization process will be automated.If the two aforementioned conditions are satisfied, then the reconstruction of the unknown vector can be performed with a small number of measurements (transmitting/receiving antennas) by solving the following convex problem [20]:
(4)min⁡∥x∥1,s.t.ii∥Ax-y∥2<δ,
where ∥x∥1 represents the norm-one of the vector x, ∥Ax-y∥2 represents the norm-two of the residual error Ax-y, and δ is an upper bound for the residual error.Additionally, it is possible to perform the reconstruction using CS even when the unknown vector is not sparse by using a transformation in a dictionaryW. The analysis-based problem is expressed as
(5)min⁡∥W*x∥1,s.t.ii∥Ax-y∥2<δ,
where W is a dictionary that allows for a sparse representation of x, and (·)* is the conjugate operator.In this particular problem, the magnitude of the gradient of the unknown reflectivity function, which describes the reflectivity properties of the object under test, has been minimized. This transformation of the unknown reflectivity vectorx[k,m,n] accepts a sparser representation more than the vector x[k,m,n], where k and m are the cross-range indices and n is the range index of a given pixel in the reconstruction domain. For this reason, a total variation (TV) transformation is used to reconstruct the image [20]. The TV problem can be written as
(6)min⁡∥x∥TV,s.t.ii∥Ax-y∥2<δ,
where ∥x∥TV represents the TV norm of the unknown vector x, and it is computed as follows:
(7)∥x∥TV=∑k,m‍∥∇x[k,m,n]∥2,
where ∇ is given by
(8)∇x[k,m,n]=[x[k+1,m,n]-x[k,m,n]x[k,m+1,n]-x[k,m,n]].
It should be noted that ∇ only considers the partial differences in the two cross-range directions. In this way, the reconstruction for each range plane can be performed independently, thus allowing the use of the TV algorithm described in [20].


## 4. Results and Discussion

The performance of the millimeter-wave imaging system is evaluated using a numerical example that simulates security threats at a standoff distance. Table1 shows the parameters used in the numerical simulation. The total number of measurements used for the reconstruction is given by the number of antennas times the number of frequencies, which is equal to 2340. The system works in a step frequency configuration, using thirteen uniformly spaced frequencies over a 6 GHz bandwidth centered at 60 GHz. The frequency step is chosen so that no aliasing replicas of the image appear in the imaging domain [9]. The PRS is depicted in Figure 3; it is formed by small perfect electric conductor (PEC) facets, which have a flat square shape of side length 3λ=15 mm. Uniform white noise is considered in the simulation, producing a signal to noise ratio of 25 dB.Table 1
Parameters for the numerical example.
Parameter
Configuration
z
0
2000λ = 8.75 m
x
0
250λ = 1.25 m
L
1
80λ = 0.4 m
L
2
250λ = 1.25 m
n
a
180
n
s
720
n
f
13
l
7.1 mm
l
r
25 mmFigure 3
Overview of the passive reflecting surfaces (PRS). The variablesX′ and Y′, expressed in meters, are relative to the center of the PRS.In this simulation, the person under test carries two types of explosive threats on the chest: (1) two high-reflectivity metallic pipes and (2) one low-reflectivity dielectric box made of TNT. In order to accelerate the speed of the numerical simulation, the reflectivity function of the person under test is projected onto four parallel planes, each one defined at a different range. These projections of the reflectivity function are considered to be the ground truth for the imaging algorithm, and they have been depicted in Figure4. The colorbar in the image indicates the absolute value of the reflectivity function divided by the average reflectivity of the whole image.Projected images of the person under test, used as ground truth by the imaging algorithm, at four different ranges: (a)x0+0 mm, (b) x0+30 mm, (c) x0+60 mm, and (d) x0+90 mm.
(a)
(b)
(c)
(d)Figure5 presents the reconstruction obtained using a traditional synthetic aperture radar (SAR) algorithm [9]. Due to the limited number of measurements taken, the recovered image does not provide information about the target.Reconstructed images of the person under test using a traditional SAR algorithm at four different ranges: (a)x0+0 mm, (b) x0+30 mm, (c) x0+60 mm, and (d) x0+90 mm.
(a)
(b)
(c)
(d)The reconstructed image of the person under test is shown in Figure6. The compressive sensing imaging algorithm successfully detects and localizes the security threats on the person's torso. Three dimensional representations of the ground truth and the reconstructed images are shown in Figure 7, which are also successful in detecting and localizing the security threats.Reconstructed images of the person under test using CS methods at four different ranges: (a)x0+0 mm, (b) x0+30 mm, (c) x0+60 mm, and (d) x0+90 mm.
(a)
(b)
(c)
(d)Three dimensional images created using voxels: (a) ground truth and (b) reconstruction.
(a)
(b)


## 5. Conclusions

This paper describes a millimeter wave imaging system, which is capable of detecting security threats at standoff distances. Unlike traditional imaging systems, in which the radar system directly interrogates the target under test, this system illuminates passive reflecting surfaces that redirect the energy of the radar towards the person under test. The passive reflecting surfaces act as a millimeter wave mirror that increases the field of view of the imaging system, thus increasing its resolution.Novel compressive sensing techniques have been used in order to reduce the amount of data collected by the system. The imaging algorithm is based on minimizing the total variation norm, which provides a quick iterative reconstruction of the person under test, thus improving traditional, noniterative, Fourier-based imaging algorithms.This system configuration provides a resolution of 7.1 mm in cross-range and 25 mm in range, and it is capable of producing 3D images.



---

*Source: 248351-2014-07-02.xml*