# Through-Wall Detection with LS-SVM under Unknown Wall Characteristics

Authors: Fangfang Wang; Yerong Zhang; Huamei Zhang
Journal: International Journal of Antennas and Propagation
 (2016)

---

## Abstract

One of the main challenges in through-wall imaging (TWI) is the presence of the walls, whose returns tend to obscure the target behind the walls and must be considered and computed in the imaging procedure. In this paper, a two-step procedure for the through-wall detection is proposed. Firstly, an effective clutter mitigation method based on singular value decomposition (SVD) is used. It does not require knowledge of the background scene or rely on accurate modeling and estimation of wall parameters. Then, TWI problem is cast as a regression one and solved by means of least-squares support vector machine (LS-SVM). The complex scattering process due to the presence of the walls is automatically included in the nonlinear relationship between the feature vector extracted from the target scattered fields and the position of the target. The relationship is obtained through a training phase using LS-SVM. Simulated results show that the proposed approach is effective. We also analyze the impacts of training samples and signal-to-noise ratio (SNR) on test detection accuracy. Simulated results reveal that the proposed LS-SVM based approach can provide comparative performances in terms of accuracy, convergence, robustness, and generalization in comparison with the support vector machine (SVM) based approach.


---

## Body

## 1. Introduction

Sensing through obstacles such as walls, doors, and other visually opaque materials using microwave signals is emerging as a powerful tool which supports a range of civilian and military applications [1]. Through-wall imaging (TWI) is a research field of increasing interest to detect and locate scatterers which are hidden beyond an opaque obstacle, such as a building wall [2].The presence of the walls makes the TWI problem different and more challenging than free space imaging. Compensation of wall effects such as attenuation and dispersion must be accounted for in imaging procedure; otherwise it will lead to the degradation of imaging such as the shift of true positions and the smearing and blurring of the image [3, 4]. This has been done by using some modified synthetic aperture radar (SAR) processing methods relying on an approximate model of the walls [5] or casting the original problem into an inverse scattering one governed by wave equations [6]. The problems of detecting the target behind the walls are usually formulated and solved as optimization problems which are solved using iterative minimization procedures. Iterative methods involve solving a well-behaved forward problem in a feedback loop. Thus, these methods are computationally expensive.All of the imaging approaches mentioned above for TWI problem are based on physical modeling. However, the lack of required data and the expense of data acquisition can limit the practical application of physics-based models. In this situation, data-driven models are good alternative, or a complement, to physics-based models for TWI problem [7]. They are characterized by their ability to capture the underlying physics of the system simply by examination of the inputs and outputs of the system. Taking advantage of data-driven model, the TWI problem is reformulated as a regression one, where the data (i.e., scatter-fields measurement) and the unknowns (i.e., the position of the object) are related by means of an approximated function to be estimated through an offline data fitting process. In this way, the effects caused by the walls can be included in the mapping obtained after the training phase. Therefore, this technology is not sensitive to a priori knowledge on the walls or the room geometry.Recently, the support vector machine (SVM) has been successfully applied in solving regression problems and has shown excellent performances. SVM is developed on the basis of statistic learning theory and exhibits good generalization capabilities thanks to an intrinsic tradeoff between fitting of the training data and model complexity [8, 9]. Nonlinear SVM provides robust solutions by mapping the input space into a higher dimensional feature space using kernel functions. However, the SVM is obtained using quadratic programming methods, which can be time-consuming and difficult to implement [10]. Least-squares support vector machine (LS-SVM) is a reformulation of the standard SVM [11]. LS-SVM adopts equality constraints instead of the inequality constraints implemented in SVM and a quadratic error term to obtain a linear set of equations in a dual space. LS-SVM performs well for its simple calculation, faster convergence, and higher precision than traditional SVM [12].The above imaging methods are based on the assumption that the scattered fields from the target have been already extracted from the total fields. Generally, the strong clutter induced by the walls and crosstalk between antennas often obscure the target electromagnetic (EM) returns. Therefore, a clutter mitigation technique is also required before applying imaging algorithms.Recently, clutter mitigation based on statistical techniques such as singular value decomposition (SVD) has been introduced in TWI applications successfully [13, 14]. The clutter mitigation technique based on SVD does not require a priori information wall parameters. The target signal is reconstructed using singular values and vectors corresponding to the target based on the different property compared to the clutter.Inspired by the advantage of LS-SVM and SVD algorithms, we introduce in this paper a two-stage procedure for detection of the target behind the walls under wall parameters ambiguities. Firstly, the total received data is dealt with to mitigate the clutter including the wall reflection and crosstalk between antennas using SVD. Then, some features are extracted from the target scattered fields and a correspondence between the feature vector and target information such as positions is established by LS-SVM using the training data. For a new detection task, the target positions are regressed by the correspondence.The remainder of this paper is organized as follows. In Section2, we present the through-wall detection problem briefly and SVD is introduced to mitigate the clutter from the B-scan radar data. Then, we recast the problem as a regression one. In Section 3, LS-SVM is introduced and used for solving the regression problem. In Section 4, numerical results for synthetic data are shown, and the performances of the proposed through-wall detection technique are assessed. Finally, some conclusions and final remarks are made.


## 2. Problem Description

TWI scenario in this paper is referring to a three-layer media model as depicted in Figure1. The first and third layers are assumed to be the free space with dielectric permittivity and magnetic permeability denoted as ε
0 and μ
0, respectively. The second layer is the nonmagnetic wall with thickness d and permittivity ε
b and conductivity is σ
b. The target is located in the investigation domain D. The antenna works in a monostatic mode; the reflected fields are collected at the same positions as the transmitting antenna, while the latter moves over a set of N different uniform positions in order to synthesize a measurement line ∑.Figure 1
Geometry of the scattering problem.When the wave propagates through the wall, it will reflect at the air-to-wall and wall-to-air interfaces causing multiple reverberations within the wall [15]. Returns also involve multiple bounces between the wall and target, so the whole receiving signal is a nonlinear system. This gives rise to ghost targets and hence increases probability of false alarm. It is another challenge in through-wall imaging and is not considered in this work. A common way to model a through-wall scene is to consider it as a linear system, so the received echo represents the total contributions of wall reflection and target scattering fields and noise written as(1)
e
n
t
=
e
n
a
t
+
e
n
w
t
+
e
n
t
t
,where e
n
a
(
t
) accounts for the antenna crosstalk and e
n
w
(
t
) is the wall reflection collectively called the clutter. e
n
t
(
t
) is the fields scattered from the target behind the walls. Upon discretizing time, data can be arranged in a M
×
N B-scan matrix as(2)
B
=
b
1
,
b
2
,
…
,
b
N
,where b
j
=
(
y
j
1
,
y
j
2
,
…
,
y
j
M
)
T represents the A-scan data collected at the jth measurement position. N and M are the number of sensor positions and time samples, respectively.In TWI scenario, the signal scattered from the target is overwhelmed by the clutter. The clutter mitigation is becoming the key to detecting the target without any wall parameters information. Therefore, the first step is to extract the target scattered fieldse
n
t
(
t
) from the total fields e
n
(
t
). The SVD of B is given as(3)
B
=
U
S
V
H
,where U
=
u
1
,
u
2
,
…
,
u
M and V
=
v
1
,
v
2
,
…
,
v
N are unitary matrices, respectively. u
i
=
u
1
i
,
u
2
i
,
…
,
u
M
i
H and v
i
=
v
1
i
,
v
2
i
,
…
,
v
N
i
H are the eigenvectors of B
B
H and B
H
B. S
=
diag
⁡
(
σ
1
,
σ
2
,
…
,
σ
r
) and σ are the singular values with the order σ
1
≥
σ
2
≥
⋯
≥
σ
r
≥
0. Therefore, B can be expressed as a linear combination of all eigencomponents:(4)
B
=
∑
n
=
1
r
E
n
=
σ
1
u
1
v
1
H
+
⋯
+
σ
r
u
r
v
r
H
,where E
n stands for the nth eigencomponent. B can be decomposed into two subspaces, target B
C and clutter B
T, respectively:(5)
B
=
B
C
+
B
T
=
∑
n
=
1
p
σ
n
u
n
v
n
H
+
∑
n
=
p
+
1
r
σ
n
u
n
v
n
H
.
p is the boundary of the two subspaces.In monostatic radar scheme, the amplitude and phase of the clutter at each antenna location across the synthetic aperture are the same. As a result, the clutter measurement vectors corresponding to various antenna locations are linearly dependent leading to the rank of the clutter data matrix being one. When applying SVD to the data matrix, the clutter, due to its strength relative to the target, resides in the first singular vector corresponding to the maximal singular value. Therefore, after removing the singular vector corresponding to the maximal singular value, the residual singular vectors are used to construct the target subspace. The target subspace can be written as(6)
B
T
=
∑
n
=
2
r
σ
n
u
n
v
n
H
.TWI aims at retrieving the location, the shape, and the dielectric property of the target from the measured fields when illuminated by known incident fields. The target scattered fields can be thought of as a function of their parameters. Therefore, TWI problem can be considered as a regression one, which involves estimating or approximating the underlying relationship from detected signal data set. Generally, one target is centered at(
p
_
x
,
p
_
y
) with a diameter of ρ
t, and the conductivity and dielectric permittivity of it are characterized by σ
t, ε
t, respectively. x is the feature vector extracted from the target scattered fields B
T such as a set of amplitude at each observation location, and y
=
(
p
_
x
,
p
_
y
,
ρ
t
,
ε
t
,
σ
t
) represents the target parameters. Once the mapping from x to y is established, we can predict the vector y from the data vector x.


## 3. LS-SVM for Function Regression

The training dataset{
(
x
1
,
y
1
)
,
…
,
(
x
N
,
y
N
)
}
⊂
R
N
×
R is given, where R
N denotes the space of the input patterns; x
i represents the scattered fields due to the target, and y
i represents the position of the target such as x-coordinate or y-coordinate (i.e., y is an element of the vector y) for the through-wall detection. The LS-SVM model for nonlinear function regression is represented in feature space as(7)
y
x
=
w
T
φ
x
+
b
.Here, the nonlinear functionφ
(
·
)
:
R
N
→
R
N
k maps the input space to a higher-dimension feature space. b is a bias term and w
∈
R
N is the weight vector. The optimization problem consists in minimizing(8)
J
w
,
e
=
1
2
w
T
w
+
γ
1
2
∑
i
=
1
N
e
i
2subject to the equality constraints(9)
y
i
=
w
T
φ
x
i
+
b
+
e
i
i
=
1
,
…
,
N
,where the fitting error is denoted by e
i. The hyperparameter γ controls the tradeoff between the smoothness of the function y and the accuracy of the fitting. The solution is obtained after constructing the Lagrangian(10)
L
w
,
b
,
e
,
α
=
J
w
,
e
-
∑
i
=
1
N
α
i
w
T
φ
x
i
+
b
+
e
i
-
y
i
,where α
i are Lagrangian multipliers. The optimality of upper function is as the following sets of linear equations:(11)
∂
L
∂
w
=
0
⟶
w
=
∑
i
=
1
N
α
i
φ
x
i
∂
L
∂
b
=
0
⟶
∑
i
=
1
N
α
i
=
0
∂
L
∂
e
i
=
0
⟶
α
i
=
γ
e
i
∂
L
∂
α
i
=
0
⟶
w
T
φ
x
i
+
b
+
e
i
-
y
i
=
0
.After eliminating variables w and e, we get matrix equations:(12)
0
1
T
1
Ω
+
γ
-
1
I
b
α
=
0
y
,where y
=
y
1
;
…
;
y
N, 1
=
1
;
…
;
1, α
=
α
1
;
…
;
α
N, and Ω
k
l
=
K
(
x
k
,
x
l
) is the kernel matrix.According to Mercer’s condition, there is mappingφ and kernel function(13)
K
x
k
,
x
l
=
φ
x
k
T
φ
x
l
.The resulting LS-SVM model for nonlinear function regression becomes(14)
y
x
=
∑
i
=
1
N
α
i
K
x
,
x
i
+
b
,where α
i and b are obtained by solving (12). The most usual kernel functions are linear, polynomial, and radial basis function (RBF). In this paper, RBF kernel is used: (15)
K
x
k
,
x
l
=
exp
⁡
-
x
k
-
x
l
2
2
σ
2
.Here σ is a constant defining the kernel width. The efficient performance of the LS-SVM model is determined by an optimal selection of the hyperparameters γ and σ, which are tuned using cross-validation in the training phase. The whole procedure of the proposed method is shown in Figure 2.Figure 2
Procedure of the through-wall detection.


## 4. Results and Discussion

In this section, the direct electromagnetic scattering problem, which is exploited to collect data, is formulated. Figure3 shows the general scenario of an arbitrary dielectric scatterer residing in a simple room. The thickness of walls is 0.2 m. The investigation domain is D
=
[
-
1.08
,
1.08
]
×
[
0.25
,
3.64
]
m
2. The conductivity σ
b and relative permittivity ε
r of the walls are set as a value of 0.003 S/m and 2.75, respectively.Figure 3
Schematic of through-wall geometry.The inside and outside regions of room (i.e., walls are excluded) are free space with dielectric permittivity and permeability denoted byε
0 and μ
0. The geometry of the room is discretized with the FDTD square cells of 1 cm length. The time resolution is 19.25
×
1
0
-
12s. The antenna (TX/RX) transmits electromagnetic radiation and measures the radiation that is scattered back to the antenna. Then it moves and synthesizes a measurement aperture L
=
2.4 m. To model the EM illumination of the room and the object in it with an UWB short pulse, the transmitter dipole antenna is fed by a 4.5 ns Gaussian pulse modulated by a 0.5 GHz sine wave.When the observation points areN
=
9, the time domain responses are sampled and arranged into a B-scan matrix B as shown in Figure 4(a). It is shown that the target scattered fields are overwhelmed by the clutter signals. Using SVD, matrix B can be decomposed into several singular elements, and the singular values of the matrix are depicted in Figure 4(b). In monostatic radar model, a homogenous wall subspace can be assumed to be spanned by the first singular vector associated with the dominant singular value. Therefore, a SVD-based method is proposed to separate the clutter and the target subspace since the first singular element is to account for the clutter returns. After the clutter mitigation based on the proposed method, the target scattered fields are reconstructed as shown in Figure 5(a). The results based on background subtraction are depicted in Figure 5(b). From Figures 5(a) and 5(b), it shows that the proposed method is as effective as background subtraction in removing the clutter and revealing the target behind the wall, while it does not know a prior knowledge of the background scene. Figure 5(a) shows that there are also signals caused by mutual interaction between the wall and target after the clutter mitigation process.Figure 4
B-scan data for the total received fields and its singular values: (a) B-scan data and (b) singular values.
(a)
(b)Figure 5
B-scan data for the target scattering fields based on (a) the proposed method and (b) background subtraction.
(a)
(b)Once obtaining the scattered fields from the target, some features can be extracted from these fields. In this paper, we consider singular values except for the first one, time of arrival, and amplitude of the scattered fields at individual observation location as the features. Ultimately, a vector containing these features is fed into LS-SVM predictor as an input and the position of the scatterer is the output.The LS-SVM must be trained by exploiting some input-output patterns suitably chosen so as to give an exhaustive representation of the through-wall detection problem. To do this, the training dataset is achieved by repeated simulations with varying the position of the target as follows:(16)
p
_
x
n
=
-
1
+
n
Δ
p
_
x
,
n
=
0,1
,
…
,
N
x
,
Δ
p
_
x
=
0.1
m
,
p
_
y
n
=
0.3
+
n
Δ
p
_
y
,
n
=
0,1
,
…
,
N
y
,
Δ
p
_
y
=
0.1
m
,where N
x
=
21, N
y
=
33. The total 693 samples are divided into two sets, where 643 samples are selected as the training set to construct LS-SVM model and the remaining is selected as the testing set.Since the assessment is focused on the target detection through a wall, the degree of matching between the actual and predicted position comprised of horizon and depth is firstly evaluated using the constructed LS-SVM model. The results are shown in Figure6. The more the points are concentrated around the diagonal line, the better the prediction of the model. The results show that most of the points locate near the diagonal line.Figure 6
Predicted values of (a) horizon and (b) depth of the target center versus the actual ones.
(a)
(b)We define the mean square error (MSE) as(17)
M
S
E
=
1
N
t
e
s
t
∑
i
=
1
N
t
e
s
t
p
i
-
q
i
2
,where N
t
e
s
t is the number of the samples in testing dataset and p
i and q
i are the ith actual and predicted position of the target, respectively. Firstly, the MSE against the number of training samples using LS-SVM compared with the results obtained by SVM are illustrated in Figure 7.Figure 7
Impact of the number of training samples on MSE for (a) horizon and (b) depth of the target center.
(a)
(b)The number of training samples changes from 200 to 643; it can be observed that the MSE decreases rapidly as a function of the number of training examples. Results in Figure7 declare that it is necessary to increase the training samples to get higher detection accuracy if possible. However, the CPU time increases with the growth of the training set size. As a result, the tradeoff between two factors should be considered intensively.The above-mentioned scene is for noiseless circumstance. In reality, it is important that the reconstruction algorithm is robust under noise contamination. To do this, after the feature vectors are generated, they are successively corrupted by an additive Gaussian noise with mean value equal to zero and variance fixed according to the desired signal-to-noise ratio (SNR). The MSE as a function of SNR is displayed in Figure8. It can be shown that the added noise almost has no effect on through-wall detection when using SVM model. However, the predicted error decreases rapidly as a function of the SNR when using LS-SVM model.Figure 8
MSE as a function of SNR for (a) horizon and (b) depth of the target center.
(a)
(b)


## 5. Conclusion

In this paper, we have presented a LS-SVM based method for target detection in through-wall imaging. SVD is firstly applied to remove or at least significantly suppress the clutter. Then, a flexible mapping between the scattered fields and the position of the target behind the walls is build using LS-SVM model. Simulated results reveal that the proposed approach does not assume prior knowledge of the scene nor the wall EM characteristics. The proposed approach can achieve high detection accuracy and keep low computation complexity. Moreover, the results confirm that the LS-SVM predictor can closely match the SVM predictor performances in terms of accuracy, convergence, robustness, and generalization.



---

*Source: 2790782-2016-04-28.xml*