# A Through-the-Wall Radar Imaging Method Based on a Realistic Model

Authors: Tian Jin; Alexander Yarovoy
Journal: International Journal of Antennas and Propagation
 (2015)

---

## Abstract

An image focusing method based on a realistic model for a wall is proposed for through-the-wall radar imaging using a multiple-input multiple-output array. A technique to estimate the wall parameters (i.e., position, thickness, and permittivity) from the radar returns is developed and tested. The estimated wall properties are used in the developed penetrating image formation to form images. The penetrating image formation developed is computationally efficient to realize real-time imaging, which does not depend on refraction points. The through-the-wall imaging method is validated on simulated and real data. It is shown that the proposed method provides high localization accuracy of targets concealed behind walls.


---

## Body

## 1. Introduction

Through-the-wall imaging (TWI) and, in general, imaging of buildings internal structure have received much interest in recent years [1–4]. The crucial point in TWI is the presence of the wall. Traditional image formations are mostly based on the assumption of homogenous propagation medium, which is no longer valid in TWI. Refraction of electromagnetic (EM) wave may occur at the front and rear surfaces of the wall, where the propagation path is not a straight line any more. Although the traditional image formations, that is, the back-projection (BP) algorithm, still work if the influence of the wall is not taken into account, the formed image of behind-the-wall targets will be blurred and displaced from their true locations, which degrades radar performance.In order to obtain a good focusing quality of behind-the-wall targets, the effect of wall should be compensated during the imaging procedure, where the wall parameters need to be known or estimated. Fortunately, radar receives backscattered signals from targets of interest and their surroundings, that is, the wall, which potentially allows extracting the surrounding information from the received echo.In TWI, the wall parameters, such as width and permittivity, are required for good quality imaging. There are some methods to estimate the wall parameters, such as using different array structures or multiple standoff distances [5], but they involve additional measurement procedures. Some approaches use iterative optimization schemes [6], but they are computationally intensive. Other authors retrieved the unknown wall parameters by using an inverse scattering algorithm, which minimizes the mean square error between the measured and the predicted reflection coefficient calculated from a forward model [7]. The inverse approaches require precise measurement with a high dynamic range and a time-consuming calibration procedure. The permittivity of a wall is commonly estimated from the amplitude of the echoes backscattered from the front surface of the wall [8]. However, this approach requires precise calibration similar to inverse methods. Therefore, estimation of wall parameters using radar return needs further study in TWI. There are several through-the-wall image formations that can be used in TWI [9–12], where the refraction point needs to be determined to calculate the propagation path between the transmitting (and/or receiving) antenna and the behind-the-wall target. A precise but time-consuming method is proposed to calculate the refraction point by solving several goniometrical functions [9]. In order to simplify the computation, an equivalent two-layer model is introduced to represent the three-layer model (air-wall-air) without introducing any error in calculation of propagation path, where only one refraction point needs to be determined [13]. However, the position of the refraction point changes with the positions of transmitting/receiving antennas. Therefore, the calculation of refraction point has to be done for each pair of transmitting/receiving antennas with respect to every imaging grid point, which is still a huge computational burden even with the equivalent two-layer model. In order to reduce the computational complexity, a novel penetrating image formation is proposed, which does not need to determine the refraction point and thus is referred to as refraction-point-free (RPF) penetrating image formation.In this paper, we propose a novel through-the-wall image formation for virtual aperture radar (VAR). Compared with synthetic aperture radar (SAR), which obtains the imaging aperture by the platform movement, a VAR system forms the imaging aperture by a multiple-input multiple-output (MIMO) array. For a MIMO array withM transmitting antennas and N receiving antennas, a virtual aperture with M×N virtual transceivers can be obtained. The data acquisition time of a virtual aperture is usually in the order of ms, which is much shorter than that of a synthetic aperture. Therefore, a moving target, defocusing in SAR images, can be well focused in VAR images because of the fact that a moving target, that is, a moving person, can be considered stationary during such a short time. Therefore, VAR systems are widely adopted in TWI. This paper is organized as follows. In Section 2, the realistic model for through-the-wall VAR is developed, which includes three wall parameters: position, width, and permittivity. In Section 3, a method to estimate the three parameters using the radar return is proposed. Based on the estimated wall parameters, the RPF penetrating image formation is proposed to form through-the-wall images. In Section 4, simulated and real data are used to validate the efficiency of our proposed method. Discussions and conclusions are given in Section 5.


## 2. Realistic Model for Through-the-Wall VAR

In a typical scenario for through-the-wall radar (TWR), a wall separates radar and targets. Figure1 illustrates the 2D imaging geometry of a linear MIMO array parallel to the wall. Using the linear MIMO array, we can obtain a linear virtual aperture to form the image in the x-y imaging plane, where x and y denote the cross-range and the down-range, respectively.Figure 1
Imaging geometry of TWR with MIMO array.The wall and behind-the-wall targets will both cause backscattered signals, where the propagation path of EM wave with respect to themth transmitting (Tx) and the nth receiving (Rx) antennas is depicted in Figure 1. Therefore, the received echo of TWR is composed of two parts: one is from the wall and the other is from behind-the-wall targets. The wall also serves as the propagation channel for the echo model of behind-the-wall targets, as depicted in Figure 2. In TWR imaging, the realistic model is characterized by three wall parameters: position, width, and permittivity. When the linear MIMO array is parallel to the wall, the wall position can be defined by the distance between the wall’s front surface and the linear MIMO array.Figure 2
Developed realistic model of TWR.For the echo from the wall, the transmitted signal will be reflected at the front and rear surfaces, which is denoted as the 1st and 2nd reflections, respectively. The echo phase history of the 1st and 2nd reflections can be defined by the electrical length of the two-way path traveled by a spherical wave from the transmitting antenna to the reflection point and back to the receiving antenna.Given themth transmitting element and the nth receiving element located at (xTm,0) and (xRn,0), respectively, the propagation path lW1 of the 1st reflection is(1)lW1=2rTXP1,where rTXP1 is the distance from the transmitting antenna and the reflection point at the front surface P1.The propagation medium of the 2nd reflection is air-wall, which refracts atP3 on the front surface and reflects at P2 on the rear surface, whose equivalent propagation path lW2 in air is(2)lW2=2rTXP3+εrrP3P2,where εr is the relative permittivity of the wall and rTXP3 and rP3P2 are the distances from the transmitting antenna to P3 and from P3 to P2, respectively.When the transmitting signal refracts at the pointP3, the incident angle θWi and the refraction angle θWt obey Snell’s law as(3)sin⁡θWisin⁡θWt=εr.According to the geometry relationship depicted in Figure 1, we have(4)rTXP3≈rTXP1-dtan⁡θWtsin⁡θWi,rP3P2=dcos⁡θWt,where d is the wall width.Substituting (4) into (2) and considering the relationship of θWi and θWt shown in (3), one derives(5)lW2=2rTXP1+dεr-sin2⁡θWi=lW1+2dεr-sin2⁡θWi.


## 3. Proposed Imaging Method

### 3.1. Estimation of Wall Parameters

The realistic model is characterized by three wall parameters: position, width, and permittivity, which should be estimated from the received echo in practice. For the MIMO array, the time delay between the 1st reflection and the 2nd reflection with respect to themth transmitting and the nth receiving elements is (6)tdm,n=2dεr-sin2⁡θWim,ncwith (7)θWim,n≈arctan⁡xRn-xTm2R,where R is the distance between the MIMO array and the wall’s front surface and can be estimated as(8)R^=12MN∑m=1M∑n=1Nc2tf2m,n-xRn-xTm2,where tf(m,n) is the time delay of the 1st reflection with respect to the mth transmitting and the nth receiving antennas.The range resolutionρr=c/(2B) is determined by the system bandwidth B, where c is the speed of EM wave in the free space. According to the imaging geometry, when the system bandwidth fulfils the condition (9)B>ccos⁡θWt2dεr,the reflections of the wall’s front and rear surfaces can be separated in the down-range (or fast-time) domain. The relative permittivity of concrete with no reinforcement is in the range of 5 to 9 under different hydration conditions [14]. Assuming εr=6 and d=0.2 m, the condition in (9) gives B>307 MHz for all refraction angles, which is usually fulfilled for most TWI radars.Based on (6), (7), and (8), we have(10)0.25c2td2m,n=d2εr-d2xRn-xTm2xRn-xTm2+4R^2,which is written in the matrix form as(11)Ap=bwith(12)A=1-xR1-xT12xR1-xT12+4R^21-xR2-xT12xR2-xT12+4R^2⋮⋮1-xRN-xTM2xRN-xTM2+4R^2MN×2,p=d2εrd22×1=p1p2,b=0.25c2td21,10.25c2td21,2⋮0.25c2td2M,NMN×1.The solution of (11) is(13)p=ATA-1ATb,where the superscript (·)T is the transpose operator.The estimations of the width and the relative permittivity are obtained accordingly as(14)d^=p2,ε^r=p1p2.


### 3.2. RPF Penetrating Image Formation

According to the equivalent two-layer model [13], the equivalent electrical length is independent of the distance between the radar and the air-wall interface. Therefore, the equivalent electrical length can be calculated on the geometry depicted in Figure 3 instead, where the MIMO array is on the front surface and P4 is the equivalent refraction point on the rear surface. θTi and θTt are the incident and the refraction angles, respectively, with the relationship of (15)sin⁡θTisin⁡θTt=εr.Therefore, the equivalent electrical length from the transmitting element to the target in air can be expressed as (16)lTxTar=rTxP4εr+rP4Tar,where rTxP4 and rP4Tar are the distances from the transmitting antenna to the refraction point P4 and from the refraction point P4 to the target, respectively.Figure 3
Equivalent two-layer propagation model.According to the geometry relationship depicted in Figure5, rTxP4 and rP4Tar are calculated as(17)rTxP4=dcos⁡θTt,rP4Tar=rTxTar-rP4Q2=rTxTar-dcos⁡θTi-θTtcos⁡θTt=rTxTar-dcos⁡θTi+sin⁡θTitan⁡θTt,where rTxTar is the distance from the transmitting antenna to the target.Substituting (15), (17) into (16), lTxTar can be rewritten as (18)lTxTar=rTxTar+dεr-sin2⁡θTi-cos⁡θTi.Similarly, the equivalent electrical length from the target to the receiving element lTarRx in air can be calculated as(19)lTarRx=rTarRx+dεr-sin2⁡θRi-cos⁡θRi.When the target located at (x,y) is behind the wall, the time delay of its return with respect to the mth transmitting and the nth receiving antennas can be estimated as (20)tT1m,n=x-xTm2+y-yTm2+x-xRn2+y-yRn2+d^ε^r-sin2⁡θTim,n-cos⁡θTim,n+ε^r-sin2⁡θRim,n-cos⁡θRim,ncwith (21)θTim,n≈arctan⁡x-xTmy-yTm,θRim,n≈arctan⁡x-xRny-yRn.When the target is in front of the wall, the time delay of its return can be computed as(22)tT2m,n=x-xTm2+y-yTm2+x-xRn2+y-yRn2c.According to (20) and (22), based on the realistic model of TWR, the RPF penetrating image formation is (23)Ix,y=∑m=1M∑n=1N∫st,m,nδt-tT1m,ndt,y≥R^+0.5d^∑m=1M∑n=1N∫st,m,nδt-tT2m,ndt,y<R^+0.5d^,where I(x,y) is the formed image, t is the fast-time, and s(t,m,n) is the received signal. It is noticed that when the transmitting signal is not the impulse signal but the stepped-frequency signal, s(t,m,n) should be the range compressed signal and a phase compensation term should be involved in (23) to ensure coherent accumulation [15].


## 3.1. Estimation of Wall Parameters

The realistic model is characterized by three wall parameters: position, width, and permittivity, which should be estimated from the received echo in practice. For the MIMO array, the time delay between the 1st reflection and the 2nd reflection with respect to themth transmitting and the nth receiving elements is (6)tdm,n=2dεr-sin2⁡θWim,ncwith (7)θWim,n≈arctan⁡xRn-xTm2R,where R is the distance between the MIMO array and the wall’s front surface and can be estimated as(8)R^=12MN∑m=1M∑n=1Nc2tf2m,n-xRn-xTm2,where tf(m,n) is the time delay of the 1st reflection with respect to the mth transmitting and the nth receiving antennas.The range resolutionρr=c/(2B) is determined by the system bandwidth B, where c is the speed of EM wave in the free space. According to the imaging geometry, when the system bandwidth fulfils the condition (9)B>ccos⁡θWt2dεr,the reflections of the wall’s front and rear surfaces can be separated in the down-range (or fast-time) domain. The relative permittivity of concrete with no reinforcement is in the range of 5 to 9 under different hydration conditions [14]. Assuming εr=6 and d=0.2 m, the condition in (9) gives B>307 MHz for all refraction angles, which is usually fulfilled for most TWI radars.Based on (6), (7), and (8), we have(10)0.25c2td2m,n=d2εr-d2xRn-xTm2xRn-xTm2+4R^2,which is written in the matrix form as(11)Ap=bwith(12)A=1-xR1-xT12xR1-xT12+4R^21-xR2-xT12xR2-xT12+4R^2⋮⋮1-xRN-xTM2xRN-xTM2+4R^2MN×2,p=d2εrd22×1=p1p2,b=0.25c2td21,10.25c2td21,2⋮0.25c2td2M,NMN×1.The solution of (11) is(13)p=ATA-1ATb,where the superscript (·)T is the transpose operator.The estimations of the width and the relative permittivity are obtained accordingly as(14)d^=p2,ε^r=p1p2.


## 3.2. RPF Penetrating Image Formation

According to the equivalent two-layer model [13], the equivalent electrical length is independent of the distance between the radar and the air-wall interface. Therefore, the equivalent electrical length can be calculated on the geometry depicted in Figure 3 instead, where the MIMO array is on the front surface and P4 is the equivalent refraction point on the rear surface. θTi and θTt are the incident and the refraction angles, respectively, with the relationship of (15)sin⁡θTisin⁡θTt=εr.Therefore, the equivalent electrical length from the transmitting element to the target in air can be expressed as (16)lTxTar=rTxP4εr+rP4Tar,where rTxP4 and rP4Tar are the distances from the transmitting antenna to the refraction point P4 and from the refraction point P4 to the target, respectively.Figure 3
Equivalent two-layer propagation model.According to the geometry relationship depicted in Figure5, rTxP4 and rP4Tar are calculated as(17)rTxP4=dcos⁡θTt,rP4Tar=rTxTar-rP4Q2=rTxTar-dcos⁡θTi-θTtcos⁡θTt=rTxTar-dcos⁡θTi+sin⁡θTitan⁡θTt,where rTxTar is the distance from the transmitting antenna to the target.Substituting (15), (17) into (16), lTxTar can be rewritten as (18)lTxTar=rTxTar+dεr-sin2⁡θTi-cos⁡θTi.Similarly, the equivalent electrical length from the target to the receiving element lTarRx in air can be calculated as(19)lTarRx=rTarRx+dεr-sin2⁡θRi-cos⁡θRi.When the target located at (x,y) is behind the wall, the time delay of its return with respect to the mth transmitting and the nth receiving antennas can be estimated as (20)tT1m,n=x-xTm2+y-yTm2+x-xRn2+y-yRn2+d^ε^r-sin2⁡θTim,n-cos⁡θTim,n+ε^r-sin2⁡θRim,n-cos⁡θRim,ncwith (21)θTim,n≈arctan⁡x-xTmy-yTm,θRim,n≈arctan⁡x-xRny-yRn.When the target is in front of the wall, the time delay of its return can be computed as(22)tT2m,n=x-xTm2+y-yTm2+x-xRn2+y-yRn2c.According to (20) and (22), based on the realistic model of TWR, the RPF penetrating image formation is (23)Ix,y=∑m=1M∑n=1N∫st,m,nδt-tT1m,ndt,y≥R^+0.5d^∑m=1M∑n=1N∫st,m,nδt-tT2m,ndt,y<R^+0.5d^,where I(x,y) is the formed image, t is the fast-time, and s(t,m,n) is the received signal. It is noticed that when the transmitting signal is not the impulse signal but the stepped-frequency signal, s(t,m,n) should be the range compressed signal and a phase compensation term should be involved in (23) to ensure coherent accumulation [15].


## 4. Verification

### 4.1. Simulation Results

In the simulation, the finite-difference time-domain (FDTD) method is used to simulate the echo. The transmitting signal is the first derivative Gaussian impulse with the lowest and highest frequencies of the transmitted signal being 0.5 GHz and 2 GHz, respectively. The antenna array is composed of 1 transmitting antenna and 21 receiving antennas. The receiving antennas form a linear receiving array of 2 m length. The transmitting antenna is placed at the center of the receiving array with its coordinates (0, 0). The array is parallel to the wall with the distance to the front surface of 1.2 m. The thickness of the wall is 0.2 m and its relative permittivity is 6.25. One target, with its coordinates of (0, 1.7), is behind the wall.The envelope amplitude image of the BP imaging result is depicted in Figure4, where the effect of wall is not considered and thus the image of the target behind the wall is blurred and displaced from its true position with the error of about 0.3 m in the down-range direction.Figure 4
Imaging result using traditional imaging method.Figure 5
Imaging result of proposed method with SNR = 10 dB.The additive white Gaussian noise (AWGN) with zero-mean is added to the above FDTD simulated data to evaluate the performance of proposed estimation method under different signal-to-noise ratios (SNRs). The SNR is defined as the power of transmitted signal over the variance of AWGN. The estimated wall parameters under different SNRs are shown in Table1, where 100 Monte Carlo experiments are done for each SNR.Table 1
Mean value of the estimated wall parameters in different SNRs.
SNR (dB)
R
^ (m)
d
^ (m)
ε
^
r
5
1.19
0.220
5.355
10
1.20
0.211
6.623
15
1.20
0.205
6.368
20
1.20
0.203
6.291The imaging result of proposed method in SNR = 10 dB is depicted in Figure5. The cross-range resolution is improved from 0.20 m in Figure 4 to 0.13 m in Figure 5, and the behind-the-wall target and the rear surface are both corrected to their true positions.


### 4.2. Real Data Results

We use the real data, collected by our built through-the-wall VAR system (Figure6), to verify the proposed method further. The VAR system transmits the stepped-frequency signal from 500 MHz to 2.5 GHz. The MIMO array consists of 2 transmitting and 11 receiving antennas. The receiving antennas, spaced 0.25 m center to center, form a linear receiving array. The 2 transmitting antennas are located at the two ends of the receiving array and work sequentially. Therefore, there are 22 transmitting-receiving channels in one virtual aperture. The Archimedean spiral antennas are adopted for transmitting and receiving antennas, whose dispersive characteristic should be compensated to ensure the high resolution in down-range. The compensation procedure can be performed via the predistortion technique [16] or the calibration function method [17], where the latter is used in our system.Figure 6
Photograph of the through-the-wall VAR system.The experimental scene is depicted in Figure7. The MIMO array is 8.5 m away from the exterior wall and a trihedral reflector is in the building. The estimation results of the wall parameters are shown in Table 2.Table 2
Estimation of wall parameters for real data.
R
d
ε
r
Estimated value
8.51 m
0.27 m
5.86
True value
8.50 m
0.28 m
5.56Figure 7
Data collection scenario.The imaging results of the BP algorithm and the method proposed are compared in Figure8, where the dashed circle and the dashed lines indicate the location of the target and the building walls, respectively. If the refraction of electromagnetic wave is not considered, the locating error is about 0.4 m in the down-range direction for the target. The result of proposed RPF penetrating image formation with the estimated wall parameters provides almost correct target localization.Figure 8
Comparison of imaging results of real data with 20 dB dynamic range: (a) BP algorithm; (b) proposed method.
(a)
(b)


## 4.1. Simulation Results

In the simulation, the finite-difference time-domain (FDTD) method is used to simulate the echo. The transmitting signal is the first derivative Gaussian impulse with the lowest and highest frequencies of the transmitted signal being 0.5 GHz and 2 GHz, respectively. The antenna array is composed of 1 transmitting antenna and 21 receiving antennas. The receiving antennas form a linear receiving array of 2 m length. The transmitting antenna is placed at the center of the receiving array with its coordinates (0, 0). The array is parallel to the wall with the distance to the front surface of 1.2 m. The thickness of the wall is 0.2 m and its relative permittivity is 6.25. One target, with its coordinates of (0, 1.7), is behind the wall.The envelope amplitude image of the BP imaging result is depicted in Figure4, where the effect of wall is not considered and thus the image of the target behind the wall is blurred and displaced from its true position with the error of about 0.3 m in the down-range direction.Figure 4
Imaging result using traditional imaging method.Figure 5
Imaging result of proposed method with SNR = 10 dB.The additive white Gaussian noise (AWGN) with zero-mean is added to the above FDTD simulated data to evaluate the performance of proposed estimation method under different signal-to-noise ratios (SNRs). The SNR is defined as the power of transmitted signal over the variance of AWGN. The estimated wall parameters under different SNRs are shown in Table1, where 100 Monte Carlo experiments are done for each SNR.Table 1
Mean value of the estimated wall parameters in different SNRs.
SNR (dB)
R
^ (m)
d
^ (m)
ε
^
r
5
1.19
0.220
5.355
10
1.20
0.211
6.623
15
1.20
0.205
6.368
20
1.20
0.203
6.291The imaging result of proposed method in SNR = 10 dB is depicted in Figure5. The cross-range resolution is improved from 0.20 m in Figure 4 to 0.13 m in Figure 5, and the behind-the-wall target and the rear surface are both corrected to their true positions.


## 4.2. Real Data Results

We use the real data, collected by our built through-the-wall VAR system (Figure6), to verify the proposed method further. The VAR system transmits the stepped-frequency signal from 500 MHz to 2.5 GHz. The MIMO array consists of 2 transmitting and 11 receiving antennas. The receiving antennas, spaced 0.25 m center to center, form a linear receiving array. The 2 transmitting antennas are located at the two ends of the receiving array and work sequentially. Therefore, there are 22 transmitting-receiving channels in one virtual aperture. The Archimedean spiral antennas are adopted for transmitting and receiving antennas, whose dispersive characteristic should be compensated to ensure the high resolution in down-range. The compensation procedure can be performed via the predistortion technique [16] or the calibration function method [17], where the latter is used in our system.Figure 6
Photograph of the through-the-wall VAR system.The experimental scene is depicted in Figure7. The MIMO array is 8.5 m away from the exterior wall and a trihedral reflector is in the building. The estimation results of the wall parameters are shown in Table 2.Table 2
Estimation of wall parameters for real data.
R
d
ε
r
Estimated value
8.51 m
0.27 m
5.86
True value
8.50 m
0.28 m
5.56Figure 7
Data collection scenario.The imaging results of the BP algorithm and the method proposed are compared in Figure8, where the dashed circle and the dashed lines indicate the location of the target and the building walls, respectively. If the refraction of electromagnetic wave is not considered, the locating error is about 0.4 m in the down-range direction for the target. The result of proposed RPF penetrating image formation with the estimated wall parameters provides almost correct target localization.Figure 8
Comparison of imaging results of real data with 20 dB dynamic range: (a) BP algorithm; (b) proposed method.
(a)
(b)


## 5. Discussions and Conclusions

In this paper, the realistic model based TWI method is proposed. The way to estimate the wall parameters (i.e., position, width, and permittivity) from the radar return is suggested. A RPF penetrating image formation algorithm is proposed to focus behind-the-wall targets based on the estimated wall parameters, which is computationally efficient to realize the real-time imaging. The imaging results of simulated and real data have shown that the proposed method provides high localization accuracy of behind-the-wall targets.The realistic model with a real-valued permittivity is assumed in this paper. Although the focusing quality can be further improved when considering a more complex realistic model with a frequency-dependent and/or complex-valued permittivity with sacrifice of the computational time, it is not necessary for TWR application in practice. When TWRs are used for detection of behind-the-wall persons, a radar video is more useful than a single radar image in detection, where a computationally efficient imaging method is preferred to increase the frame-rate rather than to obtain the perfect focusing quality. For some walls, their dispersive effect will change the shape of signal and make it hard to discern the 1st and 2nd reflections from the wall. In such cases, the higher down-range resolution is required. An alternative approach without increasing the hardware burden is to correct the dispersive effect with data processing methods, that is, deconvolution methods. Nevertheless, further research is required.The proposed realistic model is possible to be extended to describe a multilayered homogenous wall or multiple homogenous walls, whose derivation is analogous to that of the equivalent two-layer model, but more restriction requirement on the down-range resolution is needed. For an inhomogeneous wall, the more complex realistic model is required to characterize the scattering of the wall and its effects on behind-the-wall targets. Fortunately, inhomogeneous walls in practice usually have the periodic structure, such as reinforced concrete walls and cinder block walls, which might simplify the modeling work based on their approximate scattering solutions [18].



---

*Source: 539510-2015-10-04.xml*