# Moving Target Imaging and Ghost Mitigation in Through-the-Wall Sensing Application

Authors: Xin Sun; BiYing Lu; ZhiMin Zhou
Journal: International Journal of Antennas and Propagation
 (2014)

---

## Abstract

Human is one kind of the most interesting targets in through-the-wall imaging. In high-resolution imaging applications, human is no longer a point target. Therefore, the previous signal models constructed by point target assumption cannot accurately describe real characteristics of EM propagation. We construct the signal models based on extended target theory in this paper. Compared with previous works, the main contributions are as follows. Firstly, human is considered as an extended target. The expressions of target scattering and wall reflections are derived. Secondly, target scattering is no longer isotropic in new model. A new kind of ghost problem which is caused by target obscuring in EM propagation is discovered and exploited. Thirdly, to improve image quality in moving target imaging, an efficient approach which adopts CFAR, clustering method, and spatial geometry relationship is proposed to remove the ghosts. The derived models are shown to agree with synthetic and experimental results. And the efficiency of proposed method is also validated, which illuminates that the ghosts are efficiently mitigated and the image quality is significantly improved.


---

## Body

## 1. Introduction

Ultra-wideband (UWB) through-wall-imaging (TWI) approaches that can detect objects through obstacles, such as walls, doors, and other opaque materials, are considered as a powerful tool for a variety of civilian and military applications [1–5]. In real TWI application, human is one kind of the most interesting targets [6–10]. Studying human imaging, especially moving human imaging, is quite significant. Therefore, the corresponding moving target imaging (MTI) technique is widely researched in TWI application [11–15].For exploiting the MTI, an accurate signal model is necessary and urgently required. Many efforts have been made to analyze the signal models to improve image quality. In [16, 17], an inverse scattering signal model is derived. Based on the model, tomography imaging method and change detection technique are adopted in MTI processing. Besides inverse scattering signal model, a back projection signal model is discussed in [18–20]. Electromagnetic synthetic experiments are given in this paper. In order to analyze the EM propagation in a building, ray multipath is computed with ATrace for a building structure, and all paths linking a transmitter-receiver are discovered in [21]. And in [22, 23], a multipath model for sensing through walls is proposed. To sum up, it is noted that all these models consider the target as point one and assume its scattering to be isotropic. Since the target is point one, effects of the target itself on EM propagation are not paid enough attention in these models.However, in high-resolution imaging radar applications [24–27], most of the targets, such as human and furniture, are no longer point targets. For example, in a common high-resolution through-the-wall radar system, the frequency range is usually 1–3 GHz or even larger. If the frequency band is 2 GHz, the range resolution is 7.5 cm. However, the size of a common human is 1.8 m × 0.5 m × 0.35 m [28, 29]. In this case, the human will occupy several resolution units. Considering the relationship between human size and resolution, human must be considered as an extended target. So, in order to describe the echo and analyze its effects more accurately, the signal models need to be reconstructed by extended target theory.In this paper, we construct the TWI signal models using the extended target theory in two typical applications: the MIMO through-the-wall mode and monostatic through-the-wall mode. Based on the models, a new kind of ghost is discovered and exploited. To improve the image quality, an efficient ghost mitigation approach is proposed. The paper is organized as follows. In Section2, compared with the conventional point target signal models, two new signal models, for MIMO and monostatic modes, are constructed by extended target theory. Then moving target imaging for extended target is discussed in Section 3. The ghost caused by target obscuring is discovered and analyzed. In Section 4, to improve the imaging quality, a ghost mitigation method which applies CFAR, clustering method, and spatial geometry relationship is proposed. Simulations and experiments are given in Section 5. Conclusions end this paper.


## 2. TWI Signal Model

After comparing with the conventional point target signal models, two new signal models for extended target, in MIMO and monostatic modes, are constructed in this section.


### 2.1. Conventional Point Target Signal Model

In TWI applications, transmissivity coefficient and the reflectivity coefficient can be measured or computed by EM models. Besides, building walls, such as brick, adobe, and poured concrete walls, can be modeled by homogeneous dielectric slabs. Its transmission coefficient,Γ, and the reflection coefficient, γ, are given in [30]. Their expressions are
(1)
Γ
=
4
e
j
(
k
1
x
-
k
0
x
)
d
(
1
+
p
01
)
(
1
+
p
10
)
(
1
+
R
01
R
10
e
2
i
k
1
x
d
)
,
(2)
γ
=
R
01
+
R
10
e
i
2
k
1
x
d
(
1
+
R
01
R
10
e
i
2
k
1
x
d
)
,
where d is the thickness of the wall and k
0
x
=
2
π
f
m
cos
⁡
φ
i
/
c and k
1
x
=
2
π
ε
f
m
cos
⁡
φ
r
/
c are the normal components of the propagation constants in the air and in the wall, respectively. Here f
m is the mth frequency point. c is the velocity of EM in free space, and ε is the relative permittivity of the wall. φ
i and φ
r are the incident angle and refraction angle in the air-wall interface. p
10 is the ratio of k
0
x to k
1
x; p
10
=
k
0
x
/
k
1
x for horizontal polarization, and p
10
=
ε
k
0
x
/
k
1
x for vertical polarization. p
01 is the reciprocal of p
10; that is, p
01
=
1
/
p
10. The expressions of variables R
10 and R
01 are
(3)
R
10
=
1
-
p
10
1
+
p
10
,
R
01
=
p
10
-
1
1
+
p
10
.Using the Taylor series method, the expanded form of (2) is
(4)
γ
=
-
R
10
+
[
1
-
(
R
10
)
2
]
R
10
·
∑
n
=
1
∞
(
R
10
)
2
n
exp
⁡
[
j
2
n
k
1
x
·
d
]
,
where n denotes the reflection times when EM wave propagates inside the wall. The first term of (4) is the reflection coefficient of wall exterior. The second term corresponds to the reflection coefficient of the first reflection from the interior and multireflection inside the wall.When EM wave propagates in the wall, the energy of multireflections is rather weak because of amplitude attenuation. Therefore, wall echo is mainly composed of the first reflections from the exterior and interior. For convenience hereinafter, the multi-reflections are neglected in the following derivation. After considering the penetration losses in the wall, the reflections coefficients for the exterior and interior can be obtained from (5). One has
(5)
Λ
exterior
=
|
R
10
|
,
Λ
interior
=
α
|
[
R
10
-
(
R
10
)
3
]
exp
⁡
(
j
2
k
1
x
·
d
)
|
.Here, the variableα denote EM penetration losses of the interior reflection.(A) MIMO Through-the-Wall Mode. Figure 1(a) is the conventional point target signal model in the MIMO mode. In the figure, we take an antenna array with two transmitters and seven receivers as an example. EM illuminated by T
1 is marked with red lines and that transmitted by T
2 is marked with blue lines.The conventional TWI point target signal models: (a) MIMO mode; (b) monostatic mode.
(a)
(b)Therefore, taking stepped frequency signal, for example, the wall echo received by antennaR is the sum of the returns corresponding to all the transmitters. One has
(6)
s
w
f
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
w
·
Z
i
(
k
m
)
·
Λ
exterior
(
x
T
i
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
1
f
(
x
T
i
,
x
R
)
)
+
A
w
·
Z
i
(
k
m
)
·
Λ
interior
(
x
T
i
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
2
f
(
x
T
i
,
x
R
)
)
]
∑
i
=
1
2
}
,wheres
w
f
(
t
,
x
R
) denotes the wall echo received by receiver R. A
w is the total factor involving EM wave propagation attenuation and the influence of the antenna footprint. k
m is the wave number with k
m
=
2
π
f
m
/
c. FT
t
↔
k
-
1
(
·
) is the inverse Fourier transform operator with respect to t. Z
i
(
·
) is the spectrum of the signal transmitted by antenna T
i. x
T
i and x
R are the positions of transmitter T
i and receiver R, with x
T
i
=
(
x
T
i
,
y
T
i
) and x
R
=
(
x
R
,
y
R
), respectively. The variables Δ
l
1
f
(
·
) and Δ
l
2
f
(
·
) are signal distance of the exterior reflection and interior reflection, respectively. Taking T
1 in Figure 1(a), for example, the expressions of Δ
l
1
f
(
·
) and Δ
l
2
f
(
·
) are(7a)
Δ
l
1
f
(
x
T
1
,
x
R
)
=
∥
x
T
1
x
W
∥
+
∥
x
W
x
R
∥
,
(7b)
Δ
l
2
f
(
x
T
1
,
x
R
)
=
∥
x
T
1
x
W
1
∥
+
(
∥
x
W
1
x
W
1
′
∥
+
∥
x
W
1
′
x
W
1
′′
∥
)
·
ε
+
∥
x
W
1
′′
x
R
∥
.In (7a), x
W is a reflection position on air-wall interface. In (7b), x
W
1 and x
W
1
′′ are refraction positions on air-wall interface and wall-air interface, respectively. x
W
1
′ is the reflection position on wall-air interface of the wall. The operation ∥
·
∥ denotes the Descartes range.Based on the model, returns of the point target scattering are(8)
s
g
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
g
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
P
,
k
m
)
·
Γ
(
x
P
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
P
,
x
R
)
)
]
∑
i
=
1
2
}
,
where s
g
(
t
,
x
R
) is the target echo received by antenna R. A
g is the factor involving target scattering amplitude, EM wave propagation attenuation, and the influence of the antenna footprint. Δ
l
g
(
·
) is the signal distance of target scattering. Taking T
1 in Figure 1(a), for example, the expression of Δ
l
g
(
·
) is
(9)
Δ
l
g
(
x
T
1
,
x
P
,
x
R
)
=
∥
x
T
1
x
A
∥
+
∥
x
A
x
A
′
∥
·
ε
+
∥
x
A
′
x
P
∥
+
∥
x
P
x
B
′
∥
+
∥
x
B
′
x
B
∥
·
ε
+
∥
x
B
x
R
∥
,
where x
A and x
B are refraction positions on air-wall interface and wall-air interface of the exterior, respectively, and x
A
′ and x
B
′ are refraction positions on wall-air interface and air-wall interface of the interior, respectively.(B) Monostatic Through-the-Wall Mode. Figure 1(b) shows the conventional point target signal model for monostatic through-wall-imaging applications. In this mode, signal is emitted and received by the same antenna. The antennas work alternately. Based on the signal model, wall reflection is computed by
(10)
s
w
f
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
(
-
j
k
m
Δ
l
2
f
(
x
T
R
)
)
A
w
·
Z
(
k
m
)
·
Λ
exterior
(
x
T
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
1
f
(
x
T
R
)
)
+
A
w
·
Z
(
k
m
)
·
Λ
interior
(
x
T
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
2
f
(
x
T
R
)
)
}
.
Here, the signal distances of the wall exterior and interior, that is, Δ
l
1
f
(
·
) and Δ
l
1
f
(
·
), are(11a)
Δ
l
1
f
(
·
)
=
2
d
1
,
(11b)
Δ
l
2
f
(
·
)
=
2
(
d
1
+
d
·
ε
)
,whered
1 is the distance from the antenna to the wall.Anonymously, the point target scattering is(12)
s
g
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
A
g
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
P
,
k
m
)
h
h
h
h
h
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
P
)
)
}
,where the signal distance of target scattering is(13)
Δ
l
g
(
x
T
R
,
x
P
)
=
2
(
∥
x
T
R
x
A
∥
+
∥
x
A
x
A
′
∥
·
ε
+
∥
x
A
′
x
P
∥
)
.


### 2.2. TWI Signal Model for Extended Target

As mentioned above, target is considered as an isotropic point target in conventional models. However, in real TWI applications, most of the targets, such as the human and furniture, are extended targets rather than point targets. And these scatterings are no longer isotropic. Therefore, the signal models for TWI applications need to be reconstructed.In the new models, we assume a simple scene where a human is placed in a four-wall room (Figure2). The human is located at x
p
=
(
x
p
,
y
p
). Considering amplitude attenuation, reflections from the interior on the back wall and the multipath components are neglected in the derivation.Reflections from the back wall in extended target signal model for the MIMO and monostatic TWI applications. (a) The MIMO mode. In this figure, the antenna array includes two transmitters and seven receivers. For convenience, only the signal channels transmitted byT
1 are marked. (b) The monostatic mode. In this figure, EM waves emitted by T
R
5 are mainly obscured by the target.
(a)
(b)(A) MIMO Through-the-Wall Mode. In the new signal model, reflection from the front wall is the same as that in point target signal model, that is, (6). But the target scattering is changed. As an extended one, the target scattering is consequently the sum of the returns from its scattering points. Therefore, the extended target returns received by antenna x
R can be expressed as
(14)
s
g
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
n
=
1
N
∑
i
=
1
2
[
A
g
(
x
U
n
)
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
U
n
,
k
m
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
·
Γ
(
x
U
n
,
x
R
,
k
m
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
U
n
,
x
R
)
)
]
∑
n
=
1
N
∑
i
=
1
2
}
,
where x
U
n denotes the nth scattering point of the target. N is the count of scatter points. A
g
(
·
) is the factor of each scattering point.Because the target is an extended one, it will partially obscure the EM waves. In Figure2, H and H′ denote the left end and right ends of the target. When EM propagates through the front wall, the waves get partially reflected by the back wall and partially obscured by the target (shown in Figure 2(a)). Now, we analyze the range of the obscured area on the back wall.Taking antennaT
1, for example, when EM is illuminated by transmitter T
1, there are two cases about obscuring (shown in Figure 2(a)). One occurs before the reflection from the back wall (marked in red lines); the other occurs after the reflection from the back wall (marked in black lines). In the first case, the obscured area, Π
11
=
(
x
ob
,
y
ob
), on the back wall conforms to (15a), (15b), and (15c). Here, variable L denotes the length of the target along the cross range axis,(15a)
Π
11
=
{
x
ob
∈
[
x
P
H
x
P
H
′
]
y
ob
=
d
1
+
d
+
d
2
+
y
T
1
x
p
>
x
T
1
+
L
2
,
(15b)
x
P
H
=
x
p
-
L
2
+
(
(
x
p
-
L
2
-
x
T
1
)
(
ε
y
p
-
ε
y
T
1
-
ε
d
-
d
)
×
(
y
T
1
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
1
-
d
)
(
y
p
-
y
T
1
-
d
1
)
-
d
d
1
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
>
x
T
1
+
L
2
,
(15c)
x
P
H
′
=
x
p
+
L
2
+
(
(
x
p
+
L
2
-
x
T
1
)
(
ε
y
p
-
ε
y
T
1
-
ε
d
-
d
)
×
(
y
T
1
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
1
-
d
)
(
y
p
-
y
T
1
-
d
1
)
-
d
d
1
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
>
x
T
1
+
L
2
.In the second case, the obscured area,Π
12, on the back wall conforms to (16a), (16b), and (16c). The derivations about (15a), (15b), and (15c) and (16a), (16b), and (16c) are given in the Appendix. One has(16a)
Π
12
=
{
x
ob
∈
[
x
Q
H
x
Q
H
′
]
y
ob
=
d
1
+
d
+
d
2
+
y
T
1
x
p
>
x
T
1
+
L
2
,
(16b)
x
Q
H
=
(
y
p
-
d
1
-
d
-
d
2
)
x
T
1
-
d
2
(
x
p
-
L
/
2
)
y
p
-
d
1
-
d
-
2
d
2
+
(
(
x
p
-
L
2
-
x
T
1
)
·
(
y
p
-
d
1
-
d
-
d
2
)
·
[
(
d
1
2
+
2
d
d
1
+
2
d
2
d
1
)
ε
+
d
2
+
2
d
2
d
-
(
d
1
ε
+
d
)
y
p
]
L
2
)
×
(
[
ε
(
2
d
1
+
d
+
2
d
2
-
y
p
)
×
(
d
1
+
2
d
+
2
d
2
-
y
p
)
-
d
d
1
]
×
(
y
p
-
d
1
-
d
-
2
d
2
)
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
x
p
>
x
T
1
+
L
2
,
(16c)
x
Q
H
′
=
(
y
p
-
d
1
-
d
-
d
2
)
x
T
1
-
d
2
(
x
p
+
L
/
2
)
y
p
-
d
1
-
d
-
2
d
2
+
(
(
x
p
+
L
2
-
x
T
1
)
·
(
y
p
-
d
1
-
d
-
d
2
)
·
[
(
d
1
2
+
2
d
d
1
+
2
d
2
d
1
)
ε
+
d
2
+
2
d
2
d
-
(
d
1
ε
+
d
)
y
p
]
L
2
)
×
(
[
ε
(
2
d
1
+
d
+
2
d
2
-
y
p
)
×
(
d
1
+
2
d
+
2
d
2
-
y
p
)
-
d
d
1
]
×
(
y
p
-
d
1
-
d
-
2
d
2
)
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
>
x
T
1
+
L
2
.Anonymously, when EM is illuminated by transmitterT
2, there are also two cases about obscuring. In the first case, the obscured area, Π
21, on the back wall conforms to(17a)
Π
21
=
{
x
ob
∈
[
x
P
H
′
x
P
H
′
′
]
y
ob
=
d
1
+
d
+
d
2
+
y
T
1
x
p
<
x
T
2
-
L
2
,
(17b)
x
P
H
′
=
x
p
-
L
2
-
(
(
x
T
2
-
x
p
+
L
2
)
(
ε
y
p
-
ε
y
T
2
-
ε
d
-
d
)
×
(
y
T
2
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
2
-
d
)
(
y
p
-
y
T
2
-
d
1
)
-
d
d
1
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
<
x
T
2
-
L
2
,
(17c)
x
P
H
′
′
=
x
p
+
L
2
-
(
(
x
T
2
-
x
p
-
L
2
)
(
ε
y
p
-
ε
y
T
2
-
ε
d
-
d
)
×
(
y
T
2
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
2
-
d
)
(
y
p
-
y
T
2
-
d
1
)
-
d
d
1
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
<
x
T
2
-
L
2
.In the second case, the obscured area, Π
22, on the back wall conforms to(18a)
Π
22
=
{
x
ob
∈
[
x
Q
H
′
x
Q
H
′
′
]
,
y
ob
=
d
1
+
d
+
d
2
+
y
T
1
,
(18b)
x
Q
H
′
=
(
y
p
-
d
1
-
d
-
d
2
)
x
T
2
-
d
2
(
x
p
-
L
/
2
)
y
p
-
d
1
-
d
-
2
d
2
+
(
(
x
T
2
-
x
p
+
L
2
)
·
(
y
p
-
d
1
-
d
-
d
2
)
·
[
(
d
1
2
+
2
d
d
1
+
2
d
2
d
1
)
ε
+
d
2
+
2
d
2
d
-
(
d
1
ε
+
d
)
y
p
]
L
2
)
×
(
[
ε
(
2
d
1
+
d
+
2
d
2
-
y
p
)
×
(
d
1
+
2
d
+
2
d
2
-
y
p
)
-
d
d
1
]
×
(
y
p
-
d
1
-
d
-
2
d
2
)
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
<
x
T
2
-
L
2
,
(18c)
x
Q
H
′
′
=
(
y
p
-
d
1
-
d
-
d
2
)
x
T
2
-
d
2
(
x
p
+
L
/
2
)
y
p
-
d
1
-
d
-
2
d
2
+
(
(
x
T
2
-
x
p
-
L
2
)
·
(
y
p
-
d
1
-
d
-
d
2
)
·
[
(
d
1
2
+
2
d
d
1
+
2
d
2
d
1
)
ε
+
d
2
+
2
d
2
d
-
(
d
1
ε
+
d
)
y
p
]
L
2
)
×
(
[
ε
(
2
d
1
+
d
+
2
d
2
-
y
p
)
×
(
d
1
+
2
d
+
2
d
2
-
y
p
)
-
d
d
1
]
×
(
y
p
-
d
1
-
d
-
2
d
2
)
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
<
x
T
2
-
L
2
.Using setΠ
MIMO to denote the total obscured area, as in the above analysis, set Π
MIMO is
(19)
Π
MIMO
=
Π
11
+
Π
12
+
Π
21
+
Π
22
.It is noted that setΠ
MIMO is a function of transmitter positions and target position. Besides, signal distance in the first case is smaller than that in the second case. Furthermore, reflections from the obscured area in Π
12 and Π
22 will be scattered by the target and then received by other antennas. So in real application, EM is mainly obscured by the target in the first case, that is, Π
11 and Π
21.Therefore, reflection from the back wall is(20)
s
w
b
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
w
′
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
i
,
x
A
3
′′
,
x
R
,
k
m
)
·
Γ
(
x
A
3
′′
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
i
,
x
A
3
′′
,
x
R
)
)
]
∑
i
=
1
2
}
,
hhhhhhhhhhhhhhhhhhhhhhh
x
A
3
′′
∉
Π
MIMO
,
where s
w
b
(
t
,
x
R
) is reflections from the back wall received by antenna R. x
A
3
′′ is the reflection position on the back wall. Δ
l
b
(
·
) denotes signal distance of the back wall reflections. Taking T
1 in Figure 2(a), for example, the expression of Δ
l
b
(
·
) is
(21)
Δ
l
b
(
x
T
1
,
x
A
3
′′
,
x
R
7
)
=
∥
x
T
1
x
A
3
∥
+
∥
x
A
3
x
A
3
′
∥
·
ε
+
∥
x
A
3
′
x
A
3
′′
∥
+
∥
x
A
3
′′
x
A
4
′
∥
+
∥
x
A
4
′
x
A
4
∥
·
ε
+
∥
x
A
4
x
R
7
∥
,
where x
A
3 and x
A
3
′ are refraction positions on air-wall interface and wall-air interface of the front wall after signal transmitted by antenna T
1, respectively. x
A
4
′ and x
A
4 are refraction positions on air-wall interface and wall-air interface of the front wall before signal received by the receiver, respectively.(B) Monostatic Through-the-Wall Mode. Figure 2(b) demonstrates the signal model for monostatic through-the-wall mode when an extended target is in a four-wall room. In this mode, reflection from the front wall is the same as (10). According to the extended target theory, the extended target scattering is
(22)
s
g
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
∑
n
=
1
N
A
g
(
x
U
n
)
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
U
n
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
U
n
)
)
∑
n
=
1
N
}
.As shown in the figure, the obscured area is different from that in MIMO mode. Obscured areaΠ
monos
. on the back wall conforms to(23a)
Π
monos
.
=
{
x
ob
∈
[
x
P
H
x
P
H
′
]
x
p
-
L
2
<
x
T
R
<
x
p
+
L
2
,
y
ob
=
d
1
+
d
+
d
2
+
y
T
1
(23b)
x
P
H
=
x
p
-
L
2
+
(
(
x
p
-
L
2
-
x
T
R
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
×
(
ε
y
p
-
ε
y
T
R
-
ε
d
-
d
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
×
(
y
T
R
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
R
-
d
)
(
y
p
-
y
T
R
-
d
1
)
-
d
d
1
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
-
L
2
<
x
T
R
<
x
p
+
L
2
,
(23c)
x
P
H
′
=
x
p
+
L
2
+
(
(
x
p
+
L
2
-
x
T
R
)
iiiiiiiiiiiiiiiiiiiiiiiiiii
×
(
ε
y
p
-
ε
y
T
R
-
ε
d
-
d
)
iiiiiiiiiiiiiiiiiiiiiiiiiii
×
(
y
T
R
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
R
-
d
)
(
y
p
-
y
T
R
-
d
1
)
-
d
d
1
)
-
1
,
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
x
p
-
L
2
<
x
T
R
<
x
p
+
L
2
.Therefore, the back wall reflection is(24)
s
w
b
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
A
w
′
·
Z
(
k
m
)
·
Γ
2
(
x
T
,
x
A
3
′′
,
k
m
)
hhhhhh
·
Λ
exterior
(
x
T
R
,
x
A
3
′′
,
k
m
)
h
h
h
h
h
h
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
,
x
A
3
′′
)
)
A
w
′
}
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
A
3
′′
∉
Π
monos
.
,
where the signal distance of back wall reflection is
(25)
Δ
l
b
(
x
T
,
x
A
3
′′
)
=
2
(
d
1
+
d
ε
+
d
2
)
.


### 2.3. TWI Signal Model for Moving Extended Target

We assume that an extended target moves with constant speedv
→. At time t
0
=
0, the target is located at x
p
(
t
0
), with x
p
(
t
0
)
=
(
x
p
(
t
0
)
,
y
p
(
t
0
)
). At time t, the target moves to a new position, namely, x
p
(
t
)
=
(
x
p
(
t
)
,
y
p
(
t
)
). For convenience in the following derivation, we define v
x and v
y as the velocity components along the cross range axis and the range axis, respectively:
(26)
x
p
(
t
)
=
x
p
(
t
0
)
+
v
x
t
,
y
p
(
t
)
=
y
p
(
t
0
)
+
v
y
t
.(A) MIMO Through-the-Wall Mode. When the extended target moves in MIMO mode, target scattering is computed as
(27)
s
g
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
n
=
1
N
∑
i
=
1
2
[
A
g
(
x
U
n
(
t
)
)
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
U
n
(
t
)
,
k
m
)
·
Γ
(
x
U
n
(
t
)
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
U
n
(
t
)
,
x
R
)
)
]
∑
n
=
1
N
∑
i
=
1
2
}
,
where variable x
U
n
(
t
) is the position of the nth scattering point of the extended target at time t.When target moves in this room, reflection from the first wall is the same as (6). However, the back wall reflection is different because the obscured area varies with target movement. In this model, the obscured area on the back wall at time t is
(28)
Π
MIMO
(
t
)
=
Π
11
(
t
)
+
Π
12
(
t
)
+
Π
21
(
t
)
+
Π
22
(
t
)
.Therefore, reflection from the back wall can be derived:(29)
s
w
b
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
w
′
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
i
,
x
A
3
′′
,
x
R
,
k
m
)
·
Γ
(
x
A
3
′′
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
i
,
x
A
3
′′
,
x
R
)
)
]
∑
i
=
1
2
}
,
hhhhhhhhhhhhhhhhhhhhhhhhh
x
A
3
′′
∉
Π
MIMO
(
t
)
.So, when an extended target moves in a four-wall room, the total echo received by antennaR at time t is
(30)
S
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
w
·
Z
i
(
k
m
)
·
Λ
exterior
(
x
T
i
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
1
f
(
x
T
i
,
x
R
)
)
+
A
w
·
Z
i
(
k
m
)
·
Λ
interior
(
x
T
i
,
x
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
2
f
(
x
T
i
,
x
R
)
)
]
∑
i
=
1
2
}
+
FT
t
↔
k
-
1
{
∑
n
=
1
N
∑
i
=
1
2
[
A
g
(
x
U
n
(
t
)
)
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
U
n
(
t
)
,
k
m
)
·
Γ
(
x
U
n
(
t
)
,
x
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
U
n
(
t
)
,
x
R
)
)
]
∑
i
=
1
2
}
+
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
w
′
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
i
,
x
A
3
′′
,
x
R
,
k
m
)
·
Γ
(
x
A
3
′′
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
i
,
x
A
3
′′
,
x
R
)
)
]
∑
i
=
1
2
}
,
hhhhhhhhhhhhhhhhhhhhhhhhhhhh
A
3
′′
∉
Π
MIMO
(
t
)
.(B) Monostatic Through-the-Wall Mode. When the extended target moves in monostatic mode, target scattering is computed as
(31)
s
g
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
∑
n
=
1
N
A
g
(
x
U
n
(
t
)
)
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
U
n
(
t
)
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
U
n
(
t
)
)
)
∑
i
=
1
2
}
.Target movement does not affect the reflection from the first wall, which is denoted by (10). However, the obscured area varies with time t. So reflections from the back wall are
(32)
s
w
b
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
A
w
′
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
R
,
x
A
3
′′
)
)
}
,
hhhhhhhhhhhhhhhh
x
A
3
′′
∉
Π
monos
.
(
t
)
.As a result, in monostatic through-the-wall mode, the total echo received by antennaR at time t is
(33)
S
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
exp
⁡
(
-
j
k
m
Δ
l
2
f
(
x
T
R
)
)
A
w
·
Z
(
k
m
)
·
Λ
exterior
(
x
T
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
1
f
(
x
T
R
)
)
+
A
w
·
Z
(
k
m
)
·
Λ
interior
(
x
T
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
2
f
(
x
T
R
)
)
}
+
FT
t
↔
k
-
1
{
∑
n
=
1
N
A
g
(
x
U
n
(
t
)
)
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
U
n
(
t
)
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
U
n
(
t
)
)
)
∑
i
=
1
2
}
+
FT
t
↔
k
-
1
{
A
w
′
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
R
,
x
A
3
′′
)
)
∑
i
=
1
2
}
,
hhhhhhhhhhhhhhhhhh
A
3
′′
∉
Π
monos
.
(
t
)
.


## 2.1. Conventional Point Target Signal Model

In TWI applications, transmissivity coefficient and the reflectivity coefficient can be measured or computed by EM models. Besides, building walls, such as brick, adobe, and poured concrete walls, can be modeled by homogeneous dielectric slabs. Its transmission coefficient,Γ, and the reflection coefficient, γ, are given in [30]. Their expressions are
(1)
Γ
=
4
e
j
(
k
1
x
-
k
0
x
)
d
(
1
+
p
01
)
(
1
+
p
10
)
(
1
+
R
01
R
10
e
2
i
k
1
x
d
)
,
(2)
γ
=
R
01
+
R
10
e
i
2
k
1
x
d
(
1
+
R
01
R
10
e
i
2
k
1
x
d
)
,
where d is the thickness of the wall and k
0
x
=
2
π
f
m
cos
⁡
φ
i
/
c and k
1
x
=
2
π
ε
f
m
cos
⁡
φ
r
/
c are the normal components of the propagation constants in the air and in the wall, respectively. Here f
m is the mth frequency point. c is the velocity of EM in free space, and ε is the relative permittivity of the wall. φ
i and φ
r are the incident angle and refraction angle in the air-wall interface. p
10 is the ratio of k
0
x to k
1
x; p
10
=
k
0
x
/
k
1
x for horizontal polarization, and p
10
=
ε
k
0
x
/
k
1
x for vertical polarization. p
01 is the reciprocal of p
10; that is, p
01
=
1
/
p
10. The expressions of variables R
10 and R
01 are
(3)
R
10
=
1
-
p
10
1
+
p
10
,
R
01
=
p
10
-
1
1
+
p
10
.Using the Taylor series method, the expanded form of (2) is
(4)
γ
=
-
R
10
+
[
1
-
(
R
10
)
2
]
R
10
·
∑
n
=
1
∞
(
R
10
)
2
n
exp
⁡
[
j
2
n
k
1
x
·
d
]
,
where n denotes the reflection times when EM wave propagates inside the wall. The first term of (4) is the reflection coefficient of wall exterior. The second term corresponds to the reflection coefficient of the first reflection from the interior and multireflection inside the wall.When EM wave propagates in the wall, the energy of multireflections is rather weak because of amplitude attenuation. Therefore, wall echo is mainly composed of the first reflections from the exterior and interior. For convenience hereinafter, the multi-reflections are neglected in the following derivation. After considering the penetration losses in the wall, the reflections coefficients for the exterior and interior can be obtained from (5). One has
(5)
Λ
exterior
=
|
R
10
|
,
Λ
interior
=
α
|
[
R
10
-
(
R
10
)
3
]
exp
⁡
(
j
2
k
1
x
·
d
)
|
.Here, the variableα denote EM penetration losses of the interior reflection.(A) MIMO Through-the-Wall Mode. Figure 1(a) is the conventional point target signal model in the MIMO mode. In the figure, we take an antenna array with two transmitters and seven receivers as an example. EM illuminated by T
1 is marked with red lines and that transmitted by T
2 is marked with blue lines.The conventional TWI point target signal models: (a) MIMO mode; (b) monostatic mode.
(a)
(b)Therefore, taking stepped frequency signal, for example, the wall echo received by antennaR is the sum of the returns corresponding to all the transmitters. One has
(6)
s
w
f
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
w
·
Z
i
(
k
m
)
·
Λ
exterior
(
x
T
i
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
1
f
(
x
T
i
,
x
R
)
)
+
A
w
·
Z
i
(
k
m
)
·
Λ
interior
(
x
T
i
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
2
f
(
x
T
i
,
x
R
)
)
]
∑
i
=
1
2
}
,wheres
w
f
(
t
,
x
R
) denotes the wall echo received by receiver R. A
w is the total factor involving EM wave propagation attenuation and the influence of the antenna footprint. k
m is the wave number with k
m
=
2
π
f
m
/
c. FT
t
↔
k
-
1
(
·
) is the inverse Fourier transform operator with respect to t. Z
i
(
·
) is the spectrum of the signal transmitted by antenna T
i. x
T
i and x
R are the positions of transmitter T
i and receiver R, with x
T
i
=
(
x
T
i
,
y
T
i
) and x
R
=
(
x
R
,
y
R
), respectively. The variables Δ
l
1
f
(
·
) and Δ
l
2
f
(
·
) are signal distance of the exterior reflection and interior reflection, respectively. Taking T
1 in Figure 1(a), for example, the expressions of Δ
l
1
f
(
·
) and Δ
l
2
f
(
·
) are(7a)
Δ
l
1
f
(
x
T
1
,
x
R
)
=
∥
x
T
1
x
W
∥
+
∥
x
W
x
R
∥
,
(7b)
Δ
l
2
f
(
x
T
1
,
x
R
)
=
∥
x
T
1
x
W
1
∥
+
(
∥
x
W
1
x
W
1
′
∥
+
∥
x
W
1
′
x
W
1
′′
∥
)
·
ε
+
∥
x
W
1
′′
x
R
∥
.In (7a), x
W is a reflection position on air-wall interface. In (7b), x
W
1 and x
W
1
′′ are refraction positions on air-wall interface and wall-air interface, respectively. x
W
1
′ is the reflection position on wall-air interface of the wall. The operation ∥
·
∥ denotes the Descartes range.Based on the model, returns of the point target scattering are(8)
s
g
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
g
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
P
,
k
m
)
·
Γ
(
x
P
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
P
,
x
R
)
)
]
∑
i
=
1
2
}
,
where s
g
(
t
,
x
R
) is the target echo received by antenna R. A
g is the factor involving target scattering amplitude, EM wave propagation attenuation, and the influence of the antenna footprint. Δ
l
g
(
·
) is the signal distance of target scattering. Taking T
1 in Figure 1(a), for example, the expression of Δ
l
g
(
·
) is
(9)
Δ
l
g
(
x
T
1
,
x
P
,
x
R
)
=
∥
x
T
1
x
A
∥
+
∥
x
A
x
A
′
∥
·
ε
+
∥
x
A
′
x
P
∥
+
∥
x
P
x
B
′
∥
+
∥
x
B
′
x
B
∥
·
ε
+
∥
x
B
x
R
∥
,
where x
A and x
B are refraction positions on air-wall interface and wall-air interface of the exterior, respectively, and x
A
′ and x
B
′ are refraction positions on wall-air interface and air-wall interface of the interior, respectively.(B) Monostatic Through-the-Wall Mode. Figure 1(b) shows the conventional point target signal model for monostatic through-wall-imaging applications. In this mode, signal is emitted and received by the same antenna. The antennas work alternately. Based on the signal model, wall reflection is computed by
(10)
s
w
f
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
(
-
j
k
m
Δ
l
2
f
(
x
T
R
)
)
A
w
·
Z
(
k
m
)
·
Λ
exterior
(
x
T
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
1
f
(
x
T
R
)
)
+
A
w
·
Z
(
k
m
)
·
Λ
interior
(
x
T
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
2
f
(
x
T
R
)
)
}
.
Here, the signal distances of the wall exterior and interior, that is, Δ
l
1
f
(
·
) and Δ
l
1
f
(
·
), are(11a)
Δ
l
1
f
(
·
)
=
2
d
1
,
(11b)
Δ
l
2
f
(
·
)
=
2
(
d
1
+
d
·
ε
)
,whered
1 is the distance from the antenna to the wall.Anonymously, the point target scattering is(12)
s
g
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
A
g
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
P
,
k
m
)
h
h
h
h
h
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
P
)
)
}
,where the signal distance of target scattering is(13)
Δ
l
g
(
x
T
R
,
x
P
)
=
2
(
∥
x
T
R
x
A
∥
+
∥
x
A
x
A
′
∥
·
ε
+
∥
x
A
′
x
P
∥
)
.


## 2.2. TWI Signal Model for Extended Target

As mentioned above, target is considered as an isotropic point target in conventional models. However, in real TWI applications, most of the targets, such as the human and furniture, are extended targets rather than point targets. And these scatterings are no longer isotropic. Therefore, the signal models for TWI applications need to be reconstructed.In the new models, we assume a simple scene where a human is placed in a four-wall room (Figure2). The human is located at x
p
=
(
x
p
,
y
p
). Considering amplitude attenuation, reflections from the interior on the back wall and the multipath components are neglected in the derivation.Reflections from the back wall in extended target signal model for the MIMO and monostatic TWI applications. (a) The MIMO mode. In this figure, the antenna array includes two transmitters and seven receivers. For convenience, only the signal channels transmitted byT
1 are marked. (b) The monostatic mode. In this figure, EM waves emitted by T
R
5 are mainly obscured by the target.
(a)
(b)(A) MIMO Through-the-Wall Mode. In the new signal model, reflection from the front wall is the same as that in point target signal model, that is, (6). But the target scattering is changed. As an extended one, the target scattering is consequently the sum of the returns from its scattering points. Therefore, the extended target returns received by antenna x
R can be expressed as
(14)
s
g
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
n
=
1
N
∑
i
=
1
2
[
A
g
(
x
U
n
)
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
U
n
,
k
m
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
·
Γ
(
x
U
n
,
x
R
,
k
m
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
U
n
,
x
R
)
)
]
∑
n
=
1
N
∑
i
=
1
2
}
,
where x
U
n denotes the nth scattering point of the target. N is the count of scatter points. A
g
(
·
) is the factor of each scattering point.Because the target is an extended one, it will partially obscure the EM waves. In Figure2, H and H′ denote the left end and right ends of the target. When EM propagates through the front wall, the waves get partially reflected by the back wall and partially obscured by the target (shown in Figure 2(a)). Now, we analyze the range of the obscured area on the back wall.Taking antennaT
1, for example, when EM is illuminated by transmitter T
1, there are two cases about obscuring (shown in Figure 2(a)). One occurs before the reflection from the back wall (marked in red lines); the other occurs after the reflection from the back wall (marked in black lines). In the first case, the obscured area, Π
11
=
(
x
ob
,
y
ob
), on the back wall conforms to (15a), (15b), and (15c). Here, variable L denotes the length of the target along the cross range axis,(15a)
Π
11
=
{
x
ob
∈
[
x
P
H
x
P
H
′
]
y
ob
=
d
1
+
d
+
d
2
+
y
T
1
x
p
>
x
T
1
+
L
2
,
(15b)
x
P
H
=
x
p
-
L
2
+
(
(
x
p
-
L
2
-
x
T
1
)
(
ε
y
p
-
ε
y
T
1
-
ε
d
-
d
)
×
(
y
T
1
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
1
-
d
)
(
y
p
-
y
T
1
-
d
1
)
-
d
d
1
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
>
x
T
1
+
L
2
,
(15c)
x
P
H
′
=
x
p
+
L
2
+
(
(
x
p
+
L
2
-
x
T
1
)
(
ε
y
p
-
ε
y
T
1
-
ε
d
-
d
)
×
(
y
T
1
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
1
-
d
)
(
y
p
-
y
T
1
-
d
1
)
-
d
d
1
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
>
x
T
1
+
L
2
.In the second case, the obscured area,Π
12, on the back wall conforms to (16a), (16b), and (16c). The derivations about (15a), (15b), and (15c) and (16a), (16b), and (16c) are given in the Appendix. One has(16a)
Π
12
=
{
x
ob
∈
[
x
Q
H
x
Q
H
′
]
y
ob
=
d
1
+
d
+
d
2
+
y
T
1
x
p
>
x
T
1
+
L
2
,
(16b)
x
Q
H
=
(
y
p
-
d
1
-
d
-
d
2
)
x
T
1
-
d
2
(
x
p
-
L
/
2
)
y
p
-
d
1
-
d
-
2
d
2
+
(
(
x
p
-
L
2
-
x
T
1
)
·
(
y
p
-
d
1
-
d
-
d
2
)
·
[
(
d
1
2
+
2
d
d
1
+
2
d
2
d
1
)
ε
+
d
2
+
2
d
2
d
-
(
d
1
ε
+
d
)
y
p
]
L
2
)
×
(
[
ε
(
2
d
1
+
d
+
2
d
2
-
y
p
)
×
(
d
1
+
2
d
+
2
d
2
-
y
p
)
-
d
d
1
]
×
(
y
p
-
d
1
-
d
-
2
d
2
)
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
x
p
>
x
T
1
+
L
2
,
(16c)
x
Q
H
′
=
(
y
p
-
d
1
-
d
-
d
2
)
x
T
1
-
d
2
(
x
p
+
L
/
2
)
y
p
-
d
1
-
d
-
2
d
2
+
(
(
x
p
+
L
2
-
x
T
1
)
·
(
y
p
-
d
1
-
d
-
d
2
)
·
[
(
d
1
2
+
2
d
d
1
+
2
d
2
d
1
)
ε
+
d
2
+
2
d
2
d
-
(
d
1
ε
+
d
)
y
p
]
L
2
)
×
(
[
ε
(
2
d
1
+
d
+
2
d
2
-
y
p
)
×
(
d
1
+
2
d
+
2
d
2
-
y
p
)
-
d
d
1
]
×
(
y
p
-
d
1
-
d
-
2
d
2
)
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
>
x
T
1
+
L
2
.Anonymously, when EM is illuminated by transmitterT
2, there are also two cases about obscuring. In the first case, the obscured area, Π
21, on the back wall conforms to(17a)
Π
21
=
{
x
ob
∈
[
x
P
H
′
x
P
H
′
′
]
y
ob
=
d
1
+
d
+
d
2
+
y
T
1
x
p
<
x
T
2
-
L
2
,
(17b)
x
P
H
′
=
x
p
-
L
2
-
(
(
x
T
2
-
x
p
+
L
2
)
(
ε
y
p
-
ε
y
T
2
-
ε
d
-
d
)
×
(
y
T
2
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
2
-
d
)
(
y
p
-
y
T
2
-
d
1
)
-
d
d
1
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
<
x
T
2
-
L
2
,
(17c)
x
P
H
′
′
=
x
p
+
L
2
-
(
(
x
T
2
-
x
p
-
L
2
)
(
ε
y
p
-
ε
y
T
2
-
ε
d
-
d
)
×
(
y
T
2
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
2
-
d
)
(
y
p
-
y
T
2
-
d
1
)
-
d
d
1
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
<
x
T
2
-
L
2
.In the second case, the obscured area, Π
22, on the back wall conforms to(18a)
Π
22
=
{
x
ob
∈
[
x
Q
H
′
x
Q
H
′
′
]
,
y
ob
=
d
1
+
d
+
d
2
+
y
T
1
,
(18b)
x
Q
H
′
=
(
y
p
-
d
1
-
d
-
d
2
)
x
T
2
-
d
2
(
x
p
-
L
/
2
)
y
p
-
d
1
-
d
-
2
d
2
+
(
(
x
T
2
-
x
p
+
L
2
)
·
(
y
p
-
d
1
-
d
-
d
2
)
·
[
(
d
1
2
+
2
d
d
1
+
2
d
2
d
1
)
ε
+
d
2
+
2
d
2
d
-
(
d
1
ε
+
d
)
y
p
]
L
2
)
×
(
[
ε
(
2
d
1
+
d
+
2
d
2
-
y
p
)
×
(
d
1
+
2
d
+
2
d
2
-
y
p
)
-
d
d
1
]
×
(
y
p
-
d
1
-
d
-
2
d
2
)
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
<
x
T
2
-
L
2
,
(18c)
x
Q
H
′
′
=
(
y
p
-
d
1
-
d
-
d
2
)
x
T
2
-
d
2
(
x
p
+
L
/
2
)
y
p
-
d
1
-
d
-
2
d
2
+
(
(
x
T
2
-
x
p
-
L
2
)
·
(
y
p
-
d
1
-
d
-
d
2
)
·
[
(
d
1
2
+
2
d
d
1
+
2
d
2
d
1
)
ε
+
d
2
+
2
d
2
d
-
(
d
1
ε
+
d
)
y
p
]
L
2
)
×
(
[
ε
(
2
d
1
+
d
+
2
d
2
-
y
p
)
×
(
d
1
+
2
d
+
2
d
2
-
y
p
)
-
d
d
1
]
×
(
y
p
-
d
1
-
d
-
2
d
2
)
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
<
x
T
2
-
L
2
.Using setΠ
MIMO to denote the total obscured area, as in the above analysis, set Π
MIMO is
(19)
Π
MIMO
=
Π
11
+
Π
12
+
Π
21
+
Π
22
.It is noted that setΠ
MIMO is a function of transmitter positions and target position. Besides, signal distance in the first case is smaller than that in the second case. Furthermore, reflections from the obscured area in Π
12 and Π
22 will be scattered by the target and then received by other antennas. So in real application, EM is mainly obscured by the target in the first case, that is, Π
11 and Π
21.Therefore, reflection from the back wall is(20)
s
w
b
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
w
′
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
i
,
x
A
3
′′
,
x
R
,
k
m
)
·
Γ
(
x
A
3
′′
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
i
,
x
A
3
′′
,
x
R
)
)
]
∑
i
=
1
2
}
,
hhhhhhhhhhhhhhhhhhhhhhh
x
A
3
′′
∉
Π
MIMO
,
where s
w
b
(
t
,
x
R
) is reflections from the back wall received by antenna R. x
A
3
′′ is the reflection position on the back wall. Δ
l
b
(
·
) denotes signal distance of the back wall reflections. Taking T
1 in Figure 2(a), for example, the expression of Δ
l
b
(
·
) is
(21)
Δ
l
b
(
x
T
1
,
x
A
3
′′
,
x
R
7
)
=
∥
x
T
1
x
A
3
∥
+
∥
x
A
3
x
A
3
′
∥
·
ε
+
∥
x
A
3
′
x
A
3
′′
∥
+
∥
x
A
3
′′
x
A
4
′
∥
+
∥
x
A
4
′
x
A
4
∥
·
ε
+
∥
x
A
4
x
R
7
∥
,
where x
A
3 and x
A
3
′ are refraction positions on air-wall interface and wall-air interface of the front wall after signal transmitted by antenna T
1, respectively. x
A
4
′ and x
A
4 are refraction positions on air-wall interface and wall-air interface of the front wall before signal received by the receiver, respectively.(B) Monostatic Through-the-Wall Mode. Figure 2(b) demonstrates the signal model for monostatic through-the-wall mode when an extended target is in a four-wall room. In this mode, reflection from the front wall is the same as (10). According to the extended target theory, the extended target scattering is
(22)
s
g
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
∑
n
=
1
N
A
g
(
x
U
n
)
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
U
n
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
U
n
)
)
∑
n
=
1
N
}
.As shown in the figure, the obscured area is different from that in MIMO mode. Obscured areaΠ
monos
. on the back wall conforms to(23a)
Π
monos
.
=
{
x
ob
∈
[
x
P
H
x
P
H
′
]
x
p
-
L
2
<
x
T
R
<
x
p
+
L
2
,
y
ob
=
d
1
+
d
+
d
2
+
y
T
1
(23b)
x
P
H
=
x
p
-
L
2
+
(
(
x
p
-
L
2
-
x
T
R
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
×
(
ε
y
p
-
ε
y
T
R
-
ε
d
-
d
)
h
h
h
h
h
h
h
h
h
h
h
h
h
h
×
(
y
T
R
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
R
-
d
)
(
y
p
-
y
T
R
-
d
1
)
-
d
d
1
)
-
1
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
p
-
L
2
<
x
T
R
<
x
p
+
L
2
,
(23c)
x
P
H
′
=
x
p
+
L
2
+
(
(
x
p
+
L
2
-
x
T
R
)
iiiiiiiiiiiiiiiiiiiiiiiiiii
×
(
ε
y
p
-
ε
y
T
R
-
ε
d
-
d
)
iiiiiiiiiiiiiiiiiiiiiiiiiii
×
(
y
T
R
+
d
1
+
d
+
d
2
-
y
p
)
L
2
)
×
(
ε
(
y
p
-
y
T
R
-
d
)
(
y
p
-
y
T
R
-
d
1
)
-
d
d
1
)
-
1
,
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
i
x
p
-
L
2
<
x
T
R
<
x
p
+
L
2
.Therefore, the back wall reflection is(24)
s
w
b
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
A
w
′
·
Z
(
k
m
)
·
Γ
2
(
x
T
,
x
A
3
′′
,
k
m
)
hhhhhh
·
Λ
exterior
(
x
T
R
,
x
A
3
′′
,
k
m
)
h
h
h
h
h
h
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
,
x
A
3
′′
)
)
A
w
′
}
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
x
A
3
′′
∉
Π
monos
.
,
where the signal distance of back wall reflection is
(25)
Δ
l
b
(
x
T
,
x
A
3
′′
)
=
2
(
d
1
+
d
ε
+
d
2
)
.


## 2.3. TWI Signal Model for Moving Extended Target

We assume that an extended target moves with constant speedv
→. At time t
0
=
0, the target is located at x
p
(
t
0
), with x
p
(
t
0
)
=
(
x
p
(
t
0
)
,
y
p
(
t
0
)
). At time t, the target moves to a new position, namely, x
p
(
t
)
=
(
x
p
(
t
)
,
y
p
(
t
)
). For convenience in the following derivation, we define v
x and v
y as the velocity components along the cross range axis and the range axis, respectively:
(26)
x
p
(
t
)
=
x
p
(
t
0
)
+
v
x
t
,
y
p
(
t
)
=
y
p
(
t
0
)
+
v
y
t
.(A) MIMO Through-the-Wall Mode. When the extended target moves in MIMO mode, target scattering is computed as
(27)
s
g
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
n
=
1
N
∑
i
=
1
2
[
A
g
(
x
U
n
(
t
)
)
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
U
n
(
t
)
,
k
m
)
·
Γ
(
x
U
n
(
t
)
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
U
n
(
t
)
,
x
R
)
)
]
∑
n
=
1
N
∑
i
=
1
2
}
,
where variable x
U
n
(
t
) is the position of the nth scattering point of the extended target at time t.When target moves in this room, reflection from the first wall is the same as (6). However, the back wall reflection is different because the obscured area varies with target movement. In this model, the obscured area on the back wall at time t is
(28)
Π
MIMO
(
t
)
=
Π
11
(
t
)
+
Π
12
(
t
)
+
Π
21
(
t
)
+
Π
22
(
t
)
.Therefore, reflection from the back wall can be derived:(29)
s
w
b
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
w
′
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
i
,
x
A
3
′′
,
x
R
,
k
m
)
·
Γ
(
x
A
3
′′
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
i
,
x
A
3
′′
,
x
R
)
)
]
∑
i
=
1
2
}
,
hhhhhhhhhhhhhhhhhhhhhhhhh
x
A
3
′′
∉
Π
MIMO
(
t
)
.So, when an extended target moves in a four-wall room, the total echo received by antennaR at time t is
(30)
S
(
t
,
x
R
)
=
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
w
·
Z
i
(
k
m
)
·
Λ
exterior
(
x
T
i
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
1
f
(
x
T
i
,
x
R
)
)
+
A
w
·
Z
i
(
k
m
)
·
Λ
interior
(
x
T
i
,
x
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
2
f
(
x
T
i
,
x
R
)
)
]
∑
i
=
1
2
}
+
FT
t
↔
k
-
1
{
∑
n
=
1
N
∑
i
=
1
2
[
A
g
(
x
U
n
(
t
)
)
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
U
n
(
t
)
,
k
m
)
·
Γ
(
x
U
n
(
t
)
,
x
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
U
n
(
t
)
,
x
R
)
)
]
∑
i
=
1
2
}
+
FT
t
↔
k
-
1
{
∑
i
=
1
2
[
A
w
′
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
i
,
x
A
3
′′
,
x
R
,
k
m
)
·
Γ
(
x
A
3
′′
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
i
,
x
A
3
′′
,
x
R
)
)
]
∑
i
=
1
2
}
,
hhhhhhhhhhhhhhhhhhhhhhhhhhhh
A
3
′′
∉
Π
MIMO
(
t
)
.(B) Monostatic Through-the-Wall Mode. When the extended target moves in monostatic mode, target scattering is computed as
(31)
s
g
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
∑
n
=
1
N
A
g
(
x
U
n
(
t
)
)
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
U
n
(
t
)
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
U
n
(
t
)
)
)
∑
i
=
1
2
}
.Target movement does not affect the reflection from the first wall, which is denoted by (10). However, the obscured area varies with time t. So reflections from the back wall are
(32)
s
w
b
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
A
w
′
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
R
,
x
A
3
′′
)
)
}
,
hhhhhhhhhhhhhhhh
x
A
3
′′
∉
Π
monos
.
(
t
)
.As a result, in monostatic through-the-wall mode, the total echo received by antennaR at time t is
(33)
S
(
t
,
x
T
R
)
=
FT
t
↔
k
-
1
{
exp
⁡
(
-
j
k
m
Δ
l
2
f
(
x
T
R
)
)
A
w
·
Z
(
k
m
)
·
Λ
exterior
(
x
T
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
1
f
(
x
T
R
)
)
+
A
w
·
Z
(
k
m
)
·
Λ
interior
(
x
T
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
2
f
(
x
T
R
)
)
}
+
FT
t
↔
k
-
1
{
∑
n
=
1
N
A
g
(
x
U
n
(
t
)
)
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
U
n
(
t
)
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
U
n
(
t
)
)
)
∑
i
=
1
2
}
+
FT
t
↔
k
-
1
{
A
w
′
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
R
,
x
A
3
′′
)
)
∑
i
=
1
2
}
,
hhhhhhhhhhhhhhhhhh
A
3
′′
∉
Π
monos
.
(
t
)
.


## 3. Moving Target Imaging and Ghost Analysis

Compared with the conventional point target signal models, new characteristics can be concluded when the target is an extended one. (i) The extended target, such as human, obscures the EM propagation. As a result, reflection from the back wall changes. (ii) The obscured area varies with target movement. Therefore, the moving target imaging results will be affected. In this section, we will analyze the effect of target obscuring on moving target imaging and propose a method to remove or minimize these effects.In moving target imaging, change detection method is the first step because of its efficiency in removing the strong clutter and noise. The change detection operation is(34)
Δ
S
(
t
,
x
R
)
=
S
(
t
+
Δ
t
,
x
R
)
-
S
(
t
,
x
R
)
,
where Δ
S
(
·
) is the subtraction result after change detection.Then we make imaging processing for the subtract result,Δ
S
(
t
,
x
R
). Here we adopt the back projection (BP) imaging method [18–20]. In the imaging operation, the region of interest is divided into a finite number of pixels in range and cross range directions. The complex amplitude image value for the pixel located at x
q
=
(
x
q
,
y
q
) is obtained by applying frequency-dependent phase and weights to all the received data:
(35)
I
(
x
q
,
y
q
)
=
∬
w
·
Δ
S
(
t
,
x
R
)
·
exp
⁡
(
j
2
π
f
m
r
(
x
T
,
x
q
,
x
R
)
c
)
d
x
T
d
x
R
,
where w is the weighting function to shape the beam. r
(
x
T
,
x
q
,
x
R
) is the compensation signal distance for the pixel x
q. It is computed as
(36)
r
(
x
T
,
x
q
,
x
R
)
=
(
x
T
-
x
q
)
2
+
(
y
T
-
y
q
)
2
+
(
x
q
-
x
R
)
2
+
(
y
q
-
y
R
)
2
.According to change detection and BP method, moving target imaging for extended target in MIMO and monostatic modes can be computed.(A) MIMO Through-the-Wall Mode. By using (30), the subtract result, Δ
S
(
t
,
x
R
), after change detection is (37)
Δ
S
(
t
,
x
R
)
=
S
(
t
+
Δ
t
,
x
R
)
-
S
(
t
,
x
R
)
=
{
FT
t
↔
k
-
1
{
∑
n
=
1
N
∑
i
=
1
2
[
A
g
(
x
U
n
(
t
+
Δ
t
)
)
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
U
n
(
t
+
Δ
t
)
,
k
m
)
·
Γ
(
x
U
n
(
t
+
Δ
t
)
,
x
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
U
n
(
t
+
Δ
t
)
,
x
R
)
)
]
}
-
FT
t
↔
k
-
1
{
∑
n
=
1
N
∑
i
=
1
2
[
A
g
(
x
U
n
(
t
)
)
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
U
n
(
t
)
,
k
m
)
·
Γ
(
x
U
n
(
t
)
,
x
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
U
n
(
t
)
,
x
R
)
)
]
}
}
+
{
FT
t
↔
k
-
1
{
A
w
′
·
Z
1
(
k
m
)
·
Γ
(
x
T
1
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
1
,
x
A
3
′′
,
x
R
,
k
m
)
·
Γ
(
x
A
3
′′
,
x
R
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
b
(
x
T
1
,
x
A
3
′′
,
x
R
)
)
}
}
︸
A
3
′′
∈
(
Π
-
1
(
t
+
Δ
t
)
-
Π
-
1
(
t
)
)
∪
(
Π
-
1
(
t
)
-
Π
-
1
(
t
+
Δ
t
)
)
+
{
FT
t
↔
k
-
1
{
A
w
′
·
Z
2
(
k
m
)
·
Γ
(
x
T
2
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
2
,
x
A
3
′′
,
x
R
,
k
m
)
·
Γ
(
x
A
3
′′
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
2
,
x
A
3
′′
,
x
R
)
)
}
}
︸
A
3
′′
∈
(
Π
-
2
(
t
+
Δ
t
)
-
Π
-
2
(
t
)
)
∪
(
Π
-
2
(
t
)
-
Π
-
2
(
t
+
Δ
t
)
)
,where the first term in (37) is the subtraction result for target scattering. The second and the third terms denote the subtraction result for the back wall reflection when T
1 and T
2 work, respectively. Two sets, Π
1
(
·
) and Π
2
(
·
), denote the obscured areas corresponding to transmitter T
1 and transmitter T
2, respectively. Referring to (15a), (15b), (15c), (16a), (16b), (16c), (17a), (17b), (17c), (18a), (18b), and (18c), expressions of the two sets are(38a)
Π
1
(
t
)
=
Π
11
(
t
)
+
Π
12
(
t
)
,
(38b)
Π
2
(
t
)
=
Π
21
(
t
)
+
Π
22
(
t
)
.After BP imaging, there are two types of targets in the result; namely,(39)
I
(
x
q
,
y
q
)
=
I
target
(
x
q
,
y
q
)
+
G
T
1
(
x
q
,
y
q
)
+
G
T
2
(
x
q
,
y
q
)
=
∬
{
FT
t
↔
k
-
1
{
∑
n
=
1
N
∑
i
=
1
2
[
A
g
(
x
U
n
(
t
+
Δ
t
)
)
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
U
n
(
t
+
Δ
t
)
,
k
m
)
·
Γ
(
x
U
n
(
t
+
Δ
t
)
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
U
n
(
t
+
Δ
t
)
,
x
R
)
)
]
}
-
FT
t
↔
k
-
1
{
∑
n
=
1
N
∑
i
=
1
2
[
A
g
(
x
U
n
(
t
)
)
·
Z
i
(
k
m
)
·
Γ
(
x
T
i
,
x
U
n
(
t
)
,
k
m
)
·
Γ
(
x
U
n
(
t
)
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
i
,
x
U
n
(
t
)
,
x
R
)
)
]
}
}
×
exp
⁡
(
j
2
π
f
m
r
(
x
T
i
,
x
q
,
x
R
)
c
)
d
x
T
i
d
x
R
+
∬
{
FT
t
↔
k
-
1
{
A
w
′
·
Z
1
(
k
m
)
·
Γ
(
x
T
1
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
1
,
x
A
3
′′
,
x
R
,
k
m
)
·
Γ
(
x
A
3
′′
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
1
,
x
A
3
′′
,
x
R
)
)
}
}
︸
A
3
′′
∈
(
Π
-
1
(
t
+
Δ
t
)
-
Π
-
1
(
t
)
)
∪
(
Π
-
1
(
t
)
-
Π
-
1
(
t
+
Δ
t
)
)
×
exp
⁡
(
j
2
π
f
m
r
(
x
T
1
,
x
q
,
x
R
)
c
)
︸
A
3
′′
∈
(
Π
-
1
(
t
+
Δ
t
)
-
Π
-
1
(
t
)
)
∪
(
Π
-
1
(
t
)
-
Π
-
1
(
t
+
Δ
t
)
)
d
x
R
+
∬
{
FT
t
↔
k
-
1
{
A
w
′
·
Z
2
(
k
m
)
·
Γ
(
x
T
2
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
2
,
x
A
3
′′
,
x
R
,
k
m
)
·
Γ
(
x
A
3
′′
,
x
R
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
2
,
x
A
3
′′
,
x
R
)
)
}
}
︸
A
3
′′
∈
(
Π
-
2
(
t
+
Δ
t
)
-
Π
-
2
(
t
)
)
∪
(
Π
-
2
(
t
)
-
Π
-
2
(
t
+
Δ
t
)
)
×
exp
⁡
(
j
2
π
f
m
r
(
x
T
2
,
x
q
,
x
R
)
c
)
︸
A
3
′′
∈
(
Π
-
2
(
t
+
Δ
t
)
-
Π
-
2
(
t
)
)
∪
(
Π
-
2
(
t
)
-
Π
-
2
(
t
+
Δ
t
)
)
d
x
R
,whereI
target
(
·
) denotes the true target image. G
T
1
(
·
) and G
T
2
(
·
) are the ghosts corresponding to transmitter T
1 and transmitter T
2. It is noted that target obscuring is more severe in the first case than in the second case, so G
T
1
(
·
) and G
T
2
(
·
) are mainly contributed by target obscuring in the first case in real application.(B) Monostatic Through-the-Wall Mode. By using change detection method, the subtract result in monostatic through-the-wall mode is(40)
Δ
S
(
t
,
x
T
R
)
=
S
(
t
+
Δ
t
,
x
T
R
)
-
S
(
t
,
x
T
R
)
=
{
FT
t
↔
k
-
1
{
∑
n
=
1
N
A
g
(
x
U
n
(
t
+
Δ
t
)
)
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
U
n
(
t
+
Δ
t
)
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
U
n
(
t
+
Δ
t
)
)
)
}
-
FT
t
↔
k
-
1
{
∑
n
=
1
N
A
g
(
x
U
n
(
t
)
)
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
U
n
(
t
)
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
U
n
(
t
)
)
)
}
}
+
FT
t
↔
k
-
1
{
A
w
′
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
exp
⁡
⁡
(
-
j
k
m
Δ
l
b
(
x
T
R
,
x
A
3
′′
)
)
}
︸
A
3
′′
∈
(
Π
-
monos
.
(
t
+
Δ
t
)
-
Π
-
monos
.
(
t
)
)
∪
(
Π
-
monos
.
(
t
)
-
Π
-
monos
.
(
t
+
Δ
t
)
)
.Then the imaging result for an extended moving target is(41)
I
(
x
q
,
y
q
)
=
I
target
(
x
q
,
y
q
)
+
G
(
x
q
,
y
q
)
=
∬
{
FT
t
↔
k
-
1
{
∑
n
=
1
N
A
g
(
x
U
n
(
t
+
Δ
t
)
)
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
U
n
(
t
+
Δ
t
)
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
U
n
(
t
+
Δ
t
)
)
)
}
-
FT
t
↔
k
-
1
{
∑
n
=
1
N
A
g
(
x
U
n
(
t
)
)
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
U
n
(
t
)
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
g
(
x
T
R
,
x
U
n
(
t
)
)
)
}
}
exp
⁡
(
j
2
π
f
m
r
(
x
T
R
,
x
q
)
c
)
d
x
T
R
+
∬
FT
t
↔
k
-
1
{
A
w
′
·
Z
(
k
m
)
·
Γ
2
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
Λ
exterior
(
x
T
R
,
x
A
3
′′
,
k
m
)
·
exp
⁡
(
-
j
k
m
Δ
l
b
(
x
T
R
,
x
A
3
′′
)
)
}
︸
A
3
′′
∈
(
Π
-
monos
.
(
t
+
Δ
t
)
-
Π
-
monos
.
(
t
)
)
∪
(
Π
-
monos
.
(
t
)
-
Π
-
monos
.
(
t
+
Δ
t
)
)
d
x
T
R
,whereG
(
·
) is the ghost in monostatic through-the-wall mode.It is mentioned that when EM wave propagates through the wall, the velocity will be slow; that is,v
=
c
/
ε, which is determined by the electromagnetic parameters of the wall. In the processing by (41), the velocity is assumed to be propagating in free space; that is, v
=
c. Therefore, in the original imaging results, the target will be imaged behind its true position and the ghost will also be imaged behind the obscured area. The shift distance equals (
ε
-
1
)
·
d in range direction. In this paper, all the following processing is applied based on the original imaging results.


## 4. Ghost Mitigation in Through-the-Wall Moving Target Imaging

According to the proposed signal models, some conclusions can be obtained. For one thing, ghosts have a spatial relationship with the targets and transmitters in the original imaging results. In MIMO mode, as shown in Figure3, when ghost occurs in the first case, transmitters, target, and ghost are approximately on a line in the imaging result. When it occurs in the second case, mirror transmitter, ghost, and target are approximately on a line. In monostatic mode, as shown in Figure 4, the line of the target and ghost is approximately perpendicular to the antennas array. For another, ghosts are in essence the imaging residue of the back wall. So the values in range direction of their positions are constants.Spatial relationship among the antenna (or mirror antennas), target, and ghost in MIMO mode. (a) Ghost occurs in the first case. (b) Ghost occurs in the second case.
(a)
(b)Figure 4
Geometric relationship among the antenna, target, and ghost in monostatic mode.Here an approach is proposed to remove the ghosts. Firstly, adopt CFAR algorithm to obtain the detected imaging result. Then extract spatial positions of all the targets via clustering method. Finally judge the ghosts and remove them according to the spatial geometric relationship.(
1) CFAR Detection. To detect targets and suppress the clutter, a Gauss distribution based CFAR method is adopted firstly. When CFAR detection is finished, a new image, B, is obtained, with B
=
{
b
i
,
j
∣
i
=
1
,
…
,
M
B
,
j
=
1
,
…
,
N
B
}. Here, M
B and N
B are the number of pixels in range axis and cross range axis, respectively.(
2) Extract Spatial Positions via Clustering Method. By using the clustering method, spatial positions of all the targets can be extracted. During clustering processing, clustering centers are updated by iterative operation. They are computed by
(42)
x
l
′
=
1
k
l
′
{
x
l
×
k
l
+
i
×
[
1
+
b
q
max
⁡
⁡
(
b
q
)
]
}
,
y
l
′
=
1
k
l
′
{
y
l
×
k
l
+
j
×
[
1
+
b
q
max
⁡
⁡
(
b
q
)
]
}
,
k
l
′
=
k
l
+
1
+
b
i
,
j
max
⁡
(
b
i
,
j
)
,where(
x
l
,
y
l
) and (
x
l
′
,
y
l
′
) are the position of the lth clustering center before and after the updating operation, respectively. b
q denotes grey value of the pixel located at (
x
q
,
y
q
). max
⁡
(
b
i
j
) is the maximum grey value of all the pixels in the image. k
l and k
l
′ denote the number of the lth cluster before and after the updating operation, respectively.(
3) Ghosts Judgment and Mitigation. Steps for removing the ghost are as follows (taking the ith imaging result, for example).Step 1.
Detect the back wall in theith imaging result. Then calculate the symmetric positions of transmitters from the back wall. Transmitters and corresponding mirror transmitters are denoted by (
x
T
r
,
y
T
r
) and (
x
T
′
r
,
y
T
′
r
), respectively, where r is the index, with r
=
1,2
,
…
,
N
T. Here N
T is the number of the transmitters.
Using the change detection, CFAR, and clustering methods, extract spatial positions of all the targets and define them asP
=
(
x
p
i
,
y
p
i
), with p
=
1,2
,
…
,
N
P, where N
P denotes the number of the targets in the image.Step 2.
Read one of the targets in the image and define it asP
1
=
(
x
p
1
i
,
y
p
1
i
), with P
1
∈
P. Then read another target behind P
1, denoted by P
2
=
(
x
p
2
i
,
y
p
2
i
), with P
2
∈
P.(a)
If it is in the MIMO through-the-wall mode, then judge whether one of transmitters (or mirror transmitters) exists conforming to(43)
|
y
p
2
i
-
y
p
1
i
x
p
2
i
-
x
p
1
i
-
y
p
1
i
-
y
T
r
x
p
1
i
-
y
T
r
|
<
δ
or
|
y
p
2
i
-
y
p
1
i
x
p
2
i
-
x
p
1
i
-
y
p
1
i
-
y
T
′
r
x
p
1
i
-
y
T
′
r
|
<
δ
,
h
h
h
h
h
h
h
h
h
h
h
h
h
h
h
r
=
1,2
…
,
N
T
,
where the variable δ is a threshold for judging whether the target, ghost, and transmitters (or mirror transmitter) are approximately on a line. If the transmitter or mirror transmitter can be found, then go to Step 3, or update variable P
2.
(b)
If it is in the monostatic through-the-wall mode, judge whether the line ofP
1 and P
2 is perpendicular to the antennas array
(44)
|
y
p
2
i
-
y
p
1
i
x
p
2
i
-
x
p
1
i
-
x
T
r
-
x
T
r
-
1
y
T
r
-
y
T
r
-
1
|
<
μ
,
r
=
2,3
,
…
,
N
T
,
where the variable μ is the judgment threshold. If the line of P
1 and P
2 is perpendicular to the antennas array, then go to Step 3, or update variable P
2.Step 3.
Judge whether the position value in range direction ofP
2 is unchanged according to
(45)
|
y
p
2
i
-
y
p
2
i
-
1
|
≤
α
,
i
=
2,3
,
…
,
K
,
where the variable α is a threshold. y
p
2
i
-
1 denotes the position value in range direction of P
2 in the (
i
-
1
)th imaging result. K is the number of imaging results.Step 4.
IfP
2 conforms to (43)–(45), then it can be determined that it is a ghost. Once it is determined, all the pixels in the cluster of P
2 are reset to zero, otherwise they are reset to one.Step 5.
UpdateP
1 and P
2 until all the targets have been checked according to Steps 2 to 4.
At last, a new binary image is obtained, which is denoted byI
′
(
x
,
y
). In this image, pixels in the target cluster are set to one and pixels in ghost cluster are set to zero. Since I
′
(
x
,
y
) is a binary image, the final imaging result after ghost removed is computed as
(46)
I
~
(
x
,
y
)
=
I
(
x
,
y
)
⊗
I
′
(
x
,
y
)
,
where I
~
(
x
,
y
) denoted the final image after ghost removed. I
(
x
,
y
) is the original imaging result of moving target. The operation ⊗ denotes multiplication for each pixel.


## 5. Electromagnetic Simulations and Experiments

In this section, electromagnetic simulations and experiments are applied to validate the constructed signal model and the proposed method.(A) Electromagnetic Simulation about Obscuring in EM Propagation. In order to analyze the effect of target on EM propagation, we employ the FDTD algorithm. FDTD is a full-wave computational EM algorithm based on discretizing Maxwell’s time-domain equations, using finite differences in spatial and temporal dimensions. Because this algorithm operates in the time domain, results can be provided over a wide range of frequencies in one-time-marching run.The code we used in this paper is called XFDTD 7.0. Simulation model and field snapshot of EM propagation are shown in Figure5. In this simulation, a human model is placed in free space. Relative permittivity and conductivity of the model are 50 and 1 s/m, respectively. An incident plane wave propagates at 0° elevation with two polarization modes, vertical and horizontal, respectively. The incident plane wave is excitated by a Gauss impulse signal, which is shown at the left-top of Figure 5(a). 41 receivers are placed 2 m behind the human to receive the E-field. The total length and space length of receivers are 4 m and 0.1 m, respectively. Figure 5(b) is a field snapshot of EM propagation in the simulation. Incident plane wave, incident direction, human scattering, and human obscuring are marked in the figure. In order to analyze human obscuring quantitatively, received signals are discussed in Figure 6.Simulation model and field snapshot in FDTD calculation. (a) Left-top: Gauss impulse signal is used for exciting the wave; left-down: description about human model; right: simulation model in calculation. (b) Human obscuring effect in the field snapshot.
(a)
(b)There are two subplots in each figure. Left: received electronic field strength for all the receivers with frequency up to 3.3 GHz; right: electronic field strength for all the receivers at three separated frequencies: 1 GHz, 2 GHz, and 3 GHz. (a) V-V polarization. (b) H-H polarization.
(a)
(b)Figures6(a) and 6(b) show the electronic field strength for all the receiver sensors in V-V polarization and H-H polarization, respectively. The electronic field strength in full frequency range (up to 3.3 GHz) is shown in the left of each figure. Additionally, in the right of each figure, there are three subplots of electronic field strength at three separated frequencies: 1 GHz, 2 GHz, and 3 GHz. It is noted that, on the one hand, human scatting is not isotropic. The electronic field strength varies significantly at different receiver positions, especially where the receiver is totally obscured by the human. On the other hand, as an extended target, the human obscures EM propagation. As a result, electronic field strengths corresponding from the 17th to the 23rd receivers are relatively small compared to other receivers. Besides, with frequency enlargement (larger than 1 GHz in vertical polarization and 1.2 GHz in horizontal polarization, resp.), the diffraction effect can be negligible. In this case, human obscuring is significantly severe.(B) Experiments in MIMO Mode. Experiments about moving target imaging in MIMO mode have been done in National University of Defense Technology (NUDT). NUDT has developed a UWB MIMO through-the-wall radar system [31]. The radar system consists of separate transmit and receive antenna arrays. The receive array consists of fifteen elements with interelement spacing of 0.25 m. The transmit array consists of two antennas, located at either end of the receive array. The waveform used in this system is stepped frequency signal, with frequency range from 500 MHz to 2.5 GHz. The principle diagram of radar system is given in Figure 7.Figure 7
The principle diagram of through-the-wall radar system.Imaging data were collected by the radar system. The experimental imaging geometry is shown in Figure8. Radar is placed at the left side of a cinderblock building, and the antenna array is parallel to the side wall, at a distance of 8 m. The thickness of the cinderblock wall is 30 cm. The width of the side wall is 9.3 m. On the left side of the imaging scene there are some large trees. The experimental scene and the configuration inside the building are shown in Figure 8.Experimental imaging geometry and scene. (a) Photograph of the building. (b) Geometry of antenna array with the cinderblock building. (c) Configuration inside the building.
(a)
(b)
(c)There are two experiments about moving target imaging. In the first experiment, one person walked in the bedroom 2 m behind the wall during data collection. Length of his track is 3 m (shown in Figure9). The original imaging result after change detection is obtained via back projection (BP) imaging method. In the imaging, the frequency is from 1 GHz to 2 GHz, with 2 MHz increment. Taking the 27th imaging result, for example, there are three targets shown in the figure, including one true target and two ghosts (shown in Figure 10). In the figure, P1 is the true target, while P2 and P3 are the ghosts corresponding the left transmitter and right transmitter, respectively.One person walked in the bedroom. (a) Track of the movement. (b) Photo of the scene when one person was walking 2 m behind the wall in bedroom.
(a)
(b)Figure 10
The original moving target imaging result.The result after CFAR detection is shown in Figure11. The true target (P1) as well as two ghosts (P2, P3) has been detected in the figure. The energy of ghost P3 is even larger than true target P1. The ghosts severely affect the accuracy of detection.Figure 11
After CFAR detection of Figure10, ghosts are detected in the figure.Using the proposed method in this paper, the corresponding result is obtained in Figure12. Compared with Figure 10, it is shown that ghosts have been removed efficiently. Besides, the clutter has also been suppressed via the proposed method.Figure 12
Imaging result after processing by the proposed method in this paper.In the second experiment, two humans walked along the diagonal line in a bedroom during data collection (shown in Figure13). In this experiment, imaging parameters and imaging method are the same as the first one. Taking the 30th imaging result, for example, there are six targets in the figure, including two true targets and four ghosts (shown in Figure 14). In Figure 14, P1 and P2 are true targets. P3 and P4 are ghosts, which are caused by P1 corresponding to the right transmitter and the left transmitter, respectively. P5 and P6 are ghosts, which are caused by P2 corresponding to the right transmitter and the left transmitter, respectively. Besides the ghost problem, strong clutter also affects imaging quality severely.Two persons walked in the bedroom. (a) Track of the movement. (b) Photo of scene when two persons were walking in the bedroom.
(a)
(b)Figure 14
The original moving target imaging result. There are six targets in the image, including two true targets and four ghosts.Figure15 is the imaging result after CFAR detection. In Figure 15, two true targets (P1 and P2) are detected. However, two ghosts (P3, P4) are also present. Adopting the proposed method in this paper, we obtain the processing result in Figure 16. It is noted that all the ghosts are efficiently removed and the imaging quality is also significantly improved.Figure 15
Two true targets as well as three ghosts have been detected in the image after CFAR processing.Figure 16
Imaging result processed by the proposed method in this paper.(C) Electromagnetic Simulation in Monostatic Mode. To validate the performance of the proposed method in monostatic through-the-wall mode, we construct an electromagnetic simulation via XFDTD code. As shown in Figure 17, one human moves in a four-wall room. The length and width of the building are 6 and 5 m, respectively. The thickness of the wall is 0.3 m, and its relative permittivity is 4.2. The dimensions of the human model are 575 × 325 × 1875 mm. Its relative permittivity and conductivity are 50 and 1 s/m, respectively. There are 63 antennas in front of the wall to transmit and receive signals in monostatic mode. A 3 GHz band impulse signal is adopted in this simulation.Simulation scene and 3D model. (a) A human moves in a four-wall room in monostatic through-the-wall mode. (b) Top: 3 GHz band impulse signal in the time domain; left-bottom: the human model used in the simulation; right-bottom: 3D model of the scene.
(a)
(b)Resorting to the Fourier transformation and interpolation method, the frequency used in the imaging ranges from 0.5 to 1 GHz with 2 MHz steps. Taking the 15th imaging results, for example, Figure18 is the original imaging result. In this figure, P1 is the true target and P2 is the ghost. By using the CFAR detection, the true target P1 is detected but the ghost P2 is also present (shown in Figure 19). Adopting the proposed method in this paper, the processing result is given in Figure 20. Compared with Figures 18 and 19, it is shown that the ghosts are well mitigated and the target has a correspondingly higher image quality.Figure 18
The original moving target imaging result in the monostatic mode.Figure 19
The true target and the ghost have been detected in the image after CFAR processing.Figure 20
The ghost is mitigated by the proposed method in this paper.


## 6. Conclusion

In high-resolution TWI application, considering the relationship between target size and resolution, most of targets such as human and furniture are extended targets rather than point targets. Therefore, the previous TWI signal model which is based on point target assumption cannot accurately describe the real characteristics of EM propagation. In this paper, by using the extended target theory, two new TWI signal models for extended target, in MIMO and the monostatic modes, are constructed. The expressions of wall reflection and target scattering are derived.Compared with the conventional point target signal models, a new kind of ghost problem is discovered. It is caused by target obscuring in EM propagation. This problem brings severe effect on the imaging, especially on moving target imaging. Using the proposed models, we exploited the mechanism of the problem and analyzed the characteristics of the ghost.To improve the image quality in moving target imaging, an efficient approach which adopts CFAR, clustering method, and spatial geometry relationship is proposed to remove the ghosts. By using the synthetic and experimental data, the proposed models are shown in accordance with the processing results. Besides, the efficiency of the proposed method is also validated, which shows that the ghosts are efficiently mitigated and the image quality is significantly improved.It is mentioned that, in this paper, all the processing is applied based on the original imaging results. Because of the refraction phenomenon, the obscured area is not strictly on a line with the antenna (or mirror antenna) and the target in MIMO mode. However, in the original image, the imaged positions of the targets and ghosts are behind their true positions, that is, shiftingd
(
ε
-
1
) behind their true positions. Therefore, the refraction brings the errors, while the conventional BP image makes the imaged position errors. As a result, the influence of the two errors may cancel out, thus achieving a relatively satisfactory ghost-mitigated image by the proposed method.



---

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