1. Introduction
In recent years, the application of ultra-wideband (UWB) radar for remote sensing of human beings at short ranges through walls and ruins has led to a range of intensive studies [
1,
2,
3]. Different from thermal imaging, magnetic induction, X-ray, and some other success-restrictive methods [
4], UWB radar determines human objects by detecting his/her chest subtle motion due to respiration, without the influence of temperature, material of non-metal wall and target clothes [
2,
5]. Furthermore, because the low-frequency pulse signals occupy an extremely broad bandwidth, the electromagnetic wave generated by UWB radar is beneficial to penetrate walls and various nonmetal materials. Moreover, compared with the continuous-wave radar system, UWB radar can locate human targets with the required accuracy. Accordingly, the UWB radar is a good solution for detection and localization purposes in rescue missions after disasters such as earthquakes, mine accidents, and collapse [
2,
6,
7,
8,
9].
Short-range human targets detection using UWB radar systems through walls and ruins has been studied, e.g., in [
10,
11,
12,
13]. Up to now, the bistatic UWB radar (with one transmitting antenna and one receiving antenna) has resolved detection of a single stationary human target, yet the problem of multi-stationary target detection has been less well addressed. Several experiments for detecting multiple stationary human targets through walls have been carried out and the detection results have shown that bistatic UWB radar is able to accurately detect targets who are closest to the radar antennas only, whereas other farther targets are not able to be detected accurately. This is mainly attributed to three factors. Firstly, because of the low reflectivity of human respiration and the strong backscattered responses of obstacles, the signal to noise and clutter ratio (SNCR) of the echo signal is low. Secondly, as the transmission distance increases, the energy of the electromagnetic wave is attenuated, hence, the energy of electromagnetic waves reaching farther targets is inevitably weaker than that reaching the closest target. Thirdly, the closest target reflects the partial energy of the electromagnetic wave and, as a result of the decrease of the radar electromagnetic illumination, a shadowing region is formed and the respiration signal of farther targets located in the shadowing region are too weak to be detected [
14].
Some methods were applied for detecting the weak signal of human targets by bistatic UWB radar. For example, mean subtraction was applied to remove the stationary clutter to enhance the SNCR of the echo signal in a real scenario, but the disadvantage of mean subtraction is that this algorithm affects the scattering information of a specific human target [
13,
15]. Singular value decomposition (SVD) was also used for removing the interference of non-stationary clutter, but this method is not efficient with a low SNCR [
16]. The most common algorithm is adaptive filter [
11,
17], which needs two channels of input signals, but this algorithm ordinarily initiates an echo signal over-fitting. These signal processing algorithms above can suppress interference in a certain sense, however, for multiple stationary human targets, especially the farther targets located in the shadowing region, they do not work well.
As in previous studies, during radar detection of the respiration of stationary human targets, the distance between the radar and targets does not change. This means that the respiration signal always appears at a constant moment in the propagation time. Considering that the physical size of the human body is a certain thickness, the respiration signal of the human target always appears in several neighboring cells; thus, the signal in the area of the human target location will present a characteristic of regularity. This is different from the characteristic of noise and clutter, which are more likely to be random and disordered. According to the theory of information, entropy is a natural approach to measure the degree of order of a signal through its entropy spectrum [
18]. It is supposed that the entropy amplitude of human respiration and noise will have a certain difference. Furthermore, as a result of the radar echo signal of human respiration being non-stationary, wavelet transform (WT) is often used to analyze the signal of respiration in time-frequency decomposition. If the algorithm of entropy and WT are in mutual cooperation, wavelet entropy (WE) is expected to quantify the time dynamics with the order/disorder of the radar echo signal. In [
19], since humans and animals have different characteristic radar echos, the wavelet entropy method was used to distinguish human targets and dogs, and experimental results showed that the entropy value of the human target is much lower than a dog’s. It also indicated that the signals of human targets are more in order than that of dogs. Fugui Qi et al. also employed wavelet entropy to detect apnea in the respiratory signal of obstructive sleep apnea (OSA) by bio-radar, with surveillance results showing that the entropy value of apnea is higher than that of normal respiration [
20]. In all of the above, wavelet entropy has been applied to distinguish human targets and dogs, the normal respiration and apnea with the assistance of the time-frequency characteristic of human respiration. However, analyzing the echo signal with WE analysis for detecting multiple targets, especially in the case of the shadowing effect, has not been reported [
21,
22]. Considering that entropy is closely related to the characteristic of the time-frequency instead of its energy, WE is proposed in this paper for detecting multiple stationary human targets.
The five sections are discussed as follows: In
Section 2, the IR-UWB radar system is briefly described.
Section 3 states three experimental scenarios of multiple target detection. The method for detecting multiple human targets in the shadowing effect is proposed in
Section 4, and some data are analyzed. Then, the experiment results are presented and analyzed in
Section 5. Finally, the conclusion and some proposals for future work are provided in
Section 6.
4. Signal Processing and Analysis
4.1. Signal Pre-Processing and Analysis
The above-mentioned received waveforms were stored in a two-dimensional matrix with M × N, where M indicates the detection distance, and it was sampled at discrete instants in fast time ; each fast time corresponded with a specific distance. indicated the time signal of each fast time, and it was sampled in slow time . In order to better estimate the stationary human targets, the signal pre-processing was applied in three steps: (i) Distance accumulation was used for compressing the sample point along the fast time dimension to shorten the length of the data, and the sample point was compressed from 2048 to 200; (ii) The clutter caused by stationary objects, such as walls, usually leads to a baseline drift in those signals along the slow time dimension. For removing the clutter of stationary objects, the average of all waveforms’ values with slow time were subtracted by a slide window, and the width of the slide window was 100; (iii) A low pass filter was applied to remove substantial random noise. As the respiration signal is a kind of low-frequency narrowband signal, a finite impulse response (FIR) filter was adopted along the slow time dimension for suppressing the high-frequency noise. Then, a new matrix W could be expressed as , where row vector contained ( was sampling time in slow time, the unit was second) sampling point.
After pre-processing, the power spectrum of the data in
Figure 3a was calculated; the result was shown in
Figure 4a. Meanwhile, experimental data of
Figure 3b were processed by the same method and the power spectrum is shown in
Figure 4b. In previous studies, the magnitude of the power spectrum in the range dimension represents the degree of the signal power at different ranges [
19]. In
Figure 4a, the maximum peak appears at a distance of 3 m, which corresponds to the location of target A. However, target B cannot be detected, since he stood in the shadowing region.
Figure 4b shows that there are two peaks that appear at distances of 3 m and 6 m, respectively; the magnitude of the peak at a distance of 6 m is smaller than that at 3 m. This is because of the attenuation of the electromagnetic wave of the radar as the distance increases. Therefore, the traditional power spectrum method is not credible for detecting multiple human targets, especially in the case with the shadowing effect.
In
Figure 3a, as a result of target A partially reflecting the energy of the electromagnetic wave, the energy reflected from target B is reduced compared with the energy from the same position of a single person (only one person standing at 6 m). Fortunately, although the power spectrum of target B cannot be observed in
Figure 4, the signals of target B are also received by radar antenna.
Figure 5a shows the frequency spectrum of target B in
Figure 3a. The signals have many frequency components; a principal frequency of 0.25 corresponds to human respiration. In order to find the features of the signals of target B, who is located in the shadowing region, the signals of the same position in
Figure 3c were processed, and the result is shown in
Figure 5b. There are wider bands with more complex frequency components than
Figure 5a.
4.2. Auto-Correlation Algorithm
As previously known, correlation is closely related to the characteristic of the period of respiration rather than that of energy. In order to improve the SNCR of the signal after pre-processing, an auto-correlation algorithm is proposed, which is fundamental for wavelet entropy. The equation is expressed as follows:
Where
is the wanted signal in time dimension after the step of pre-processing,
is noise,
is a mixed signal in the time dimension. An auto-correlation algorithm is implemented for each row vector
. Since the human life signal is not correlated with noise in the time dimension,
and
are approximately zero, and as
increases, the major energy of noise focuses on the location of
,
is equal to
, and the noise is removed.
Figure 6a shows the results of the auto-correlation algorithm for the data of
Figure 5a. As mentioned above, the respiration signal appears in several cells and they have a high correlation with each other. Thus, in
Figure 6a, the stationary clutter and noise are eliminated, and the principal respiration frequency of 0.25 is reserved.
Figure 6b shows the frequency spectrum of
Figure 5b, where there is only one human target at a position of 3 m. Compared to
Figure 5b, the frequency components have been reduced, but are richer than in
Figure 6a.
4.3. Wavelet Transform
Wavelet transform is a useful tool for analyzing a non-stationary signal, such as that of radar. Continuous wavelet transform (CWT) decomposes the continuous signal into a family of functions localized in time and frequency, respectively. The decomposition is constituted by a set of functions,
, which are constructed by dilation and translation of a function
, called the mother wavelet. In this formulation, the parameter
is a dilation factor which measures the degree of compression, and the parameter
is the translation factor which determines the time location of the wavelet.
is defined as:
Compared with CWT, the discrete wavelet transform (DWT) is more computationally efficient. The main difference is that the dilation factor
and the translation factor
of the DWT can only have power of two values, i.e.,
where
is the level of DWT analysis performed and
is the parameter,
.
The coefficients of the DWT can be divided into approximation coefficients and detail coefficients. The approximation coefficients are the high-scale, low-frequency components of signal , whereas the detail coefficients are the low-scale, high-frequency components.
The discrete wavelet transform approximation coefficients of signal
at level
are expressed as follows:
where
is the scaling function associated with the wavelet function
. Similarly, the detail coefficients of level
are expressed as follows:
Based on Equations (5) and (6), the signal
is iteratively decomposed as the number of levels increases. Thus, a hierarchical set of approximations and details can be obtained through various decomposition levels. In this paper, the example of the three-level wavelet decomposition of signal
is shown in
Figure 7. Firstly,
is split into an approximation A
1 and a detail D
1. Secondly, the approximation A
1 is also split into a second-level approximation A
2 and a detail D
2. Thirdly, the approximation A
2 is also split into a third-level approximation A
3 and a detail D
3.
4.4. Wavelet Entropy Algorithm
In order to better explain the flow of signal processing by wavelet entropy algorithm, the signal vector in the slow time dimension processed by auto-correlation is considered to be
. As shown in the following, the DWT of signal
is defined as:
The wavelet coefficient
not only provides relevant information in a simple way, but also provides an estimation of the location energies at different scales. It is supposed that the sample value
corresponds to a uniform time grid. If the discrete dyadic wavelet decomposition is implemented over all resolution levels, the wavelet expansion of
can be expressed as:
where
are the wavelet coefficients of the wavelet sequence.
If the family
is an orthonormal based on the Hilbert space
, the concept of energy at the resolution level can be defined as:
The radar signal is divided into non-overlapping temporal windows of length
to follow the temporal evolution of the quantifiers. For each interval
with
samples, the appropriate signal values are assigned to the central point of the time window. The concept of the average energy of resolution level
can be defined as:
where
represents the number of wavelet coefficients at resolution
included in the time window
. The total energy of the wavelet coefficients at this time window can then be obtained by:
The relative wavelet energy that represent the energy’s probability distribution is:
Clearly,
. According to the Shannon entropy, which gives a measure of the information on any distribution for analyzing and comparing the probability distribution, the wavelet entropy is defined as:
In the frequency spectrum analysis of
Figure 6, the respiration of the human target produces a narrow-band spectrum; it indicates that the respiration signal is regular. However, the signal of the noise produces a very wide-band spectrum, which indicates that noise is disordered. According to entropy theory, a regular signal is inclined to produce lower-value entropy. Conversely, a completely random signal, such as clutter and noise, is inclined to a higher entropy value. In order to eliminate the effect of wavelet entropy, two experimental scenarios with no human target and a single human target located at 3 m are implemented. The wavelet entropy spectrum of these two scenarios is exhibited in
Figure 8.
Figure 8a shows that the entropy value of the experiment with no human target fluctuated as the range was processed. Excluding the influence of the reflection of the wall from 0 m to 1 m, the entropy value of the remaining region is from about 1.79 to 1.85. In
Figure 8b, as the target is located at 3 m, the entropy value of the target is the lowest, owing to the trailing effect of the target, and the entropy value at 3.8 m is lower than the other noise region.
5. Experiment Results and Discussion
In this section, in order to demonstrate the capability of the wavelet entropy algorithm to detect multi-stationary human targets, the adaptive filter algorithm and power spectrum were also used as references. Two cases of experiments imitating
Figure 2 were conducted and one variable was assigned in each experimental case. For the first experimental case, the variable is the distance of target B to radar antennas. Three experiment scenarios were described in the following: For experiment I, target A stood 3 m in front of the antennas, and target B stood 5 m in the shadowing region of target A. For experiment II and experiment III, target A remained in the same position as in experiment I, but the distance of target B to the radar antennas was 6 m and 7 m, respectively. For the second experimental case, the variable is the intersection angle of target A and target B. There were five experiment scenarios. Target A stood at 2.5 m in front of radar antennas and the distance of target B to the radar antennas was 4 m, but their intersection angles were 10°, 15°, 20°, 25°, and 30°.
Real experimental data were processed by the wavelet entropy algorithm and two reference algorithms based on the adaptive filter and power spectrum. The experiment results of the first and the second experimental case are shown in
Figure 9 and
Figure 10. In
Figure 9, only one periodically-varied strip that corresponds to the respiratory-motion response is obvious at the range of 3 m in
Figure 9a,d,g. The results of the power spectrum in
Figure 9b,e,h are similar to the above pseudo-color, with only one peak located at 3 m away from the radar in each Figure, respectively, which are correlated with target A. The energy of target B (at 5 m, 6 m, and 7 m in the first experimental case), who was located in the shadowing region, is negligible compared to that of target A. Thus, the method based on an adaptive filter and energy detection could easily induce a false negative result.
According to the result of
Figure 8, the entropy value displays a significant difference in the condition with a human target or without a human target. In this section, the proposed algorithm is also applied for locating the distance of multi-stationary human targets.
Figure 9c,f,i show the results of the first case of experimental data processed by wavelet entropy. In these three Figures, the entropy value at 3 m is the lowest, which indicates that the position of target A is accurately estimated. Likewise, another lower entropy value—5 m in
Figure 9c, 6 m in
Figure 9f, and 7 m in
Figure 9i—indicates the position of target B in the three experiments. However, it is noteworthy that the entropy values within 1 m behind target A are lower than the rest of the remaining region in the two Figures; it is easy to incorrectly deduce that one human target is located in this region. According to the theory of electromagnetic wave propagation, as a result of multipath effects and tailing effects, the ground reflects the respiration signals of target A, and the signals behind target A are regular and have the same frequency as the respiration of target A. Furthermore, the respiration signals appear in several neighboring cells and the physical size of the human body is not ignorable. In
Figure 9c, the entropy pit of target A is from 2.835 to 3.24 and the entropy pit of target B is from 4.905 to 5.31. In
Figure 9f, the entropy pit of target A is from 2.61 to 3.15 and the entropy pit of target B is from 5.76 to 6.25. In
Figure 9i, they are from 2.61 to 3.10 and from 6.93 to 7.33; the six ranges of pits correspond to the thickness of the human body. For the lower entropy value that is behind target A, the range of the pits does not obviously agree with the thickness of the human target. Thus, the algorithm of the wavelet entropy can detect two stationary human targets accurately despite one person being in the shadowing region.
The results of the second experimental case are shown in
Figure 10. From
Figure 10a–i, the intersection angle of target A and target B is 10°, 15°, and 20°, and target B is located in the shadowing region of target A. The results of the experiment processed by the proposed and referenced algorithms complied with the expected results. Even though the distance of the human target to the radar antennas is decreased, due to the shadowing effect of target A, target B cannot be detected by an adaptive filter and the power spectrum. Under the condition that the intersection angle of two targets reduces to 10°, target A reflects almost all of the energy of the electromagnetic wave to target B, and target B can still be detected by the wavelet entropy algorithm. As their intersection angle increases, the energy of the echo signal of target B is stronger, thus, when the intersection angle is 25° and 30°, in
Figure 10j,m, the periodically-varied strip corresponding to the respiratory-motion response of target B appears at the range of 3.7 m and 3.8 m, respectively. There are also small peaks of target B located at 3.7 m and 3.8 m away from the radar in each power spectrum of
Figure 10k,n, respectively. Thus, when the intersection angle is less than 20°, the wavelet entropy is more effective for detecting multi-stationary human targets.
6. Conclusions
In the case of multi-stationary human target detection through a brick wall using bistatic UWB radar, since the person closest to radar antennas partially reflects the energy of the electromagnetic wave, the person farther from the radar antennas is not detected accurately, especially when in the shadowing region of the closest person. Considering that the signals of the farther person are contained in the echo signal, and its frequency spectrum is different from the characteristic of noise, wavelet entropy is proposed as a new quantitative algorithm to detect human targets in the shadowing region. This algorithm utilizes the difference between periodic respiration and random noise in the time-frequency domain and calculates the entropy based on the complexity of energy after wavelet transform. In this paper, as the respiration of the human target exhibits more rhythmic behavior compared with noise, the entropy of the human target was significantly lower than in other regions. Through analyzing the two cases of experimental results, a conclusion is drawn: the reference algorithms based on adaptive filtering and the power spectrum can only detect the closest person; they are incapable of detecting the person in the shadowing region. The wavelet entropy algorithm not only detects the closest person and locates their distance, but also detects the farther human target in the shadowing region and gives an accurate distance. This is a breakthrough in detecting multi-stationary human targets. Yet, it must be admitted that the bistatic radar with one transmitter antenna and one receiver antenna restricts the expansion of the multi-stationary human target’s detection. Therefore, in future research, we will consider applying an UWB radar array for detecting multi-stationary targets. We think that information fusion of multiple channels may essentially improve the subtle signal of the remaining shadowed human targets in this kind of situation.