This article is part of the series Radar and Sonar Networks.

Open Access Research

Research on micro-feature extraction algorithm of target based on terahertz radar

Zhengwu Xu*, Jian Tu, Jin Li and Yiming Pi

Author Affiliations

School of Electronic Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China

For all author emails, please log on.

EURASIP Journal on Wireless Communications and Networking 2013, 2013:77  doi:10.1186/1687-1499-2013-77


The electronic version of this article is the complete one and can be found online at: http://jwcn.eurasipjournals.com/content/2013/1/77


Received:4 December 2012
Accepted:3 February 2013
Published:19 March 2013

© 2013 Xu et al.; licensee Springer.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Micro-Doppler motion of a target is an important characteristic in high-resolution radar observation. The target feature extraction of micro motion has already been applied to many aspects of radar research. In this article, general model is established for the echo signal of a target with micro-motion. Combination of time-frequency analysis method, a method using Radon transformation to detect the parameters of lines and sinusoidal curves to estimate motion parameters of target is proposed, and the estimation of reflection coefficients of scatterers is completed through nonlinear least squares and CLEAN algorithm. The simulation result shows that this method of Radon transformation has the advantages of high precision and strong anti-noise and can extract the parameters well. The model of echo signal and method of parameter estimation are useful for radar target detection and identification.

Keywords:
Micro-Doppler; Time-frequency analysis; Feature extraction; Radon transformation; NLS algorithm

1. Introduction

Mechanical vibration or rotation of structures in a target may induce frequency modulation on returned signals and generate sidebands about the center frequency of the target’s Doppler frequency [1]. The modulation due to vibrations and rotations is called micro-Doppler phenomenon. Micro-Doppler phenomenon is very common in nature, such as the human heartbeat, vibration, or spin of missile warheads, etc.… while in terahertz band micro-Doppler phenomenon is particularly significant. The micro-Doppler effect enables us to determine the dynamic properties of the target and it offers a new approach for the analysis of target signatures. Micro-Doppler features serve as additional target features that are complementary to those made available by existing methods. The micro-Doppler effect can be used to identify specific types of vehicles, and determine their movement and the speed of their engines.

In this article, considering the micro-motion of target in the terahertz band, the echo model is established [2], time-frequency transformation and Radon transformation are applied to extract the micro characteristic parameters [3-6], Radon transformation is generally used to detect straight, based on Radon transformation, we propose a method using Radon transformation to detect the parameters of lines and sinusoidal curves to estimate the motion parameters of target. And the estimation of reflection coefficient of scatterers is completed through nonlinear least squares (NLS) and the CLEAN algorithm [2,7]. These are very helpful for radar target detection and identification [8].

2. Target echo modeling

We assume a target of multiple scatterers containing macro- and micro-motion to do a compound movement on the radar radial direction. In general, targets always do non-uniform linear motion in radical direction, which can be deemed as high-order motion, the formula. Its movement rule is a higher order polynomial of time t or the infinite series of t. Weierstrass quantitative shows that arbitrary radial rule of motion may be a finite polynomial of t approximation

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M1">View MathML</a>

(1)

where n is the number of finite polynomial, micro objectives of the target can be approximated as vibration or rotation movement, its rule of motion is

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M2">View MathML</a>

(2)

where A is the vibration amplitude, B is the frequency of vibration, φ is the initial rotation angle, when t = t, the distance between target and the radar is?

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M3">View MathML</a>

(3)

Then the echo of target can be expressed as

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M4">View MathML</a>

(4)

where σi is the scattering coefficient of the ith scatterer, f0 is the radar carrier frequency, c denotes the speed of light. Here, we only consider the condition of j = 1 and j = 2.

When j = 2, and the body of target is doing a uniformly accelerated motion, that is <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M5">View MathML</a>, then the rule of target’s motion can be approximated as

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M6">View MathML</a>

(5)

and the phase of the baseband signal is

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M7">View MathML</a>

(6)

where λ denotes radar wavelength, according to the definition of micro-Doppler, the micro-Doppler frequency of target equal to the derivative of phase with respect to time, then the instantaneous frequency of the signal can be expressed as

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M8">View MathML</a>

(7)

By using time-frequency transformation with respect to the echo of target, one can get its micro-Doppler frequency. Figure 1 is obtained in the simulation conditions as follows: radar carrier frequency f0 = 340 GHz, sampling frequency fs = 4096 Hz, sampling points N = 512, and observation time t = 6 s. Assume that there are four scatterers, and v0 = 0.03 m/s2, a = 0.01 m/s2, σ1 = 1, σ2 = 0.9, σ3 = 1, σ4 = 1, A1 = B1 = φ1 = 0, A2 = A4 = 0.1, A4 = 0.2, B2 = 3 rad/s, B3 = 2 rad/s, B4 = 1 rad/s, φ2 = φ3 = φ4 = 0.6 rad.

thumbnailFigure 1. Micro-Doppler modulation induced by target.

3. Parameter estimation

3.1. Micro-feature extraction based on Radon transformation

Previously, we established the echo model of target with micro-motion, and time-frequency transformation of the echoed signal is applied to gain the time-varying micro-Doppler frequency features of target. This section is to estimate the motion parameters of target with micro-motion, namely to extract the curves of target’s motion from the time-varying micro-Doppler frequency image.

We assume that the set of parameters to be estimated is θ = v0, a, A, B, φ, as the straight line in Figure 1 is only related to parameters (v0, a), we can divide it into a two-dimensional matrix and three-dimensional matrix to estimate, respectively, which can greatly reduce the calculation. First we will estimate (v0, a) as shown in Figure 1.

Straight line detection can be realized by Radon transformation. A projection of a two-dimensional function f(x, y) is a linear integral in a certain direction. In general, the Radon transformation of f(x, y) is defined as the linear integral of f along a line. In two-dimensional space, Radon transformation can be defined by Equation (8)

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M9">View MathML</a>

(8)

where D is the entire image plane of xy, f(x, y) is the gradation of point (x, y) in the image, ρ is the distance between the origin of coordinates and the line, θ is the angle between the line and x-axis, δ is a unit pulse function.

Radon transformation can be seen as the projection of the image in the ρθ space, of which each point corresponds to a straight line in the image space. And Radon transformation is the integral of the image pixels on each straight line, then each straight line in the image will become a bright spot in the ρθ space, which turn the line detection into detecting highlights in Radon domain.

As can be seen from Figure 2a, there are two straight lines to be detected, the result of Radon transformation can be showed at Figure 2b. On Radon domain, there are peaks which indicate that there are lines on image domain.

thumbnailFigure 2. Line detection using Radon transformation. (a) Original image. (b) Radon domain. (c) Result image.

From Figure 2b, it is easy to determine peaks’ position, which is interpreted by ρ and θ. After detecting the value of ρ and θ, we can detect the straight lines on image domain. The result is showed in Figure 2c.

As seen in the result, Radon transformation can detect a straight line very accurately in the case of low SNR. And we can detect the parameters of the line in Figure 1 that the linear parameters estimated are m = 45.5, c = 68.2 (m is the slope, c is the intercept), respectively. Taking m = 4a/λ and c = 2v0/λ into consideration we have <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M10">View MathML</a>. Estimate the parameters of the curve after estimating the parameters of the line. As the linear parameters have been estimated, while all the scatterers of target have the same macro velocity and acceleration, then the amplitude, vibration frequency, initial phase of the curve in the figure only related to parameters A, B, φ. Compensate the phase to correct the curve in the figure to a sinusoidal curve?

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M11">View MathML</a>

(9)

Now time-frequency transformation of the echo signal is carried out in the same simulation condition as Figure 1, as shown in Figure 3.

thumbnailFigure 3. Micro-Doppler modulation induced by target.

There are a few methods for micro-Doppler feature extraction such as WV peak detection method [9] and the normalized first moment method [10] that are all based on time-frequency analysis methods. Here, Radon transformation detection method is chosen.

The basic idea of Radon transformation detecting curve is to do a linear integral along a specific path in a plane.

According to Equation (7), the target micro-Doppler frequency is <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M12">View MathML</a> after compensating the phase, and the amplitude of the sinusoidal curve in Figure 3 is <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M13">View MathML</a>, the vibration frequency is B. Now we let

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M14">View MathML</a>

(10)

Then the amplitude, frequency, initial phase of the sinusoidal curve in Figure 3 is determined by the three space coordinates x, y, z, respectively.

The steps to parameterize the sinusoidal curves are given as follows:

1. Set the range of x, y, z, and establish a discrete parameter space.

2. Establish a three-dimensional accumulator array P(X, Y, Z), and initialize all elements of P(X, Y, Z) to zero.

3. Calculate the value of p for each point in the parameter space x, y, z according to Equation (10) at t.

4. Take over all values of t to calculate the linear integral <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M15">View MathML</a> of the time-varying micro-Doppler frequency image along the curve p(t).

5. Find the local peak of the accumulator to obtain the spatial coordinate position of (X, Y, Z).

After Radon transformation, each curve in the time-varying micro-Doppler frequency image will generate a corresponding peak in the parameter domain, the parameters of each scatterer can be estimated by detecting the position of peaks in the parameter domain.

Figure 4a,b shows the parameter domain image after transforming when there is no noise, and Figures 4c,d the parameter domain image with noise.

thumbnailFigure 4. The parameter domain normalized image after Radon transformation. (a) The parameter domain normalized image without noise. (b) The parameter domain image without noise. (c) The parameter domain normalized image with noise. (d) The parameter domain image with noise.

As seen from Figure 4, in the image domain, the scatterers’ echo curve entangled with each other, but in the parameter domain, there are four distinct peaks corresponded to the curves. There may be a certain deviation between the true value and the parameter value extracted from the parameter domain after Radon transformation, in order to improve the accuracy of the estimated parameter, we let the estimated value to be the initial value, and search again with a small step in its nearby space, then the accurate estimated value of the parameter can be obtained.

Simulation results show that the estimated values of the parameters are the same whether there is noise or no noise, and the space coordinates of the parameters are shown in Figure 5.

thumbnailFigure 5. The space coordinates of the estimated parameters.

The coordinates of the parameters corresponded to the four scatterers are (0,0,0), (0.1,3,0.6), (0.2,2,0.6), (0.1,1,0.6), then A1 = x1 = 0, A2 = x2 = 0.1, A3 = x3 = 0.2, A4 = x4 = 0.1, B1 = y1 = 0, B2 = y2 = 3, B3 = y3 = 2, B4 = y4 = 0, C1 = z1 = 0, C2 = C3 = C4 = 0.6. And (0, 0, 0) represents the sinusoidal line with amplitude, frequency, and initial phase are zero, which is equivalent to a straight line.

As can be seen from the above results, the Radon transformation detection method has the advantages of high precision and strong anti-noise. Its detection accuracy is determined by the size of the space grid of the three-dimensional accumulator, if the grid is small, the accuracy will be high, but the computational complexity of the algorithm grows rapidly as well.

The simulation comparison result of WV peak detection method and the normalized first moment method are shown in Figure 6.

thumbnailFigure 6. Comparison of three methods of parameters estimation. (a) Micro-Doppler modulation induced by target. (b) Radon transformation detection method. (c)WV peak detection method. (d) The normalized first moment method.

From the figure, we can know that in the condition of single scatterer without noise, the three methods can all extract micro-Doppler information very accurately.

As can be seen from Figure 7, while adding noise, both the WV peak detection method and the normalized first moment method have been unable to accurately estimate the parameters of target micro-Doppler.

thumbnailFigure 7. Comparison of three methods of parameters estimation with SNR = −6 dB. (a) Micro-Doppler modulation induced by target with SNR = −6 dB. (b) Radon transformation detection method. (c)WV peak detection method. (d) The normalized first moment method.

As seen from Figure 8, In the condition of multiple scatterers, the normalized first moment method has been completely ineffective, and the envelope of WV peak detection method is close to the target micro-Doppler frequency, but still cannot make accurate estimates. While the Radon transformation detection method proposed in the article can maintain high accuracy even in the condition of multiple scatterers with strong noise.

thumbnailFigure 8. The simulation results of three methods with two scatterers. (a) Micro-Doppler modulation induced by target with two scatterers. (b) Radon transformation detection method. (c)WV peak detection method. (d) The normalized first moment method.

3.2. Scattering coefficient estimating based on NLS

In cases when the parameters of scatterers have been estimated, the NLS algorithm may be applied to the data with previously estimated scattering coefficient responses removed via the CLEAN algorithm. The echo of target is

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M16">View MathML</a>

(11)

then the ith iteration of NLS may be expressed as

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M17">View MathML</a>

(12)

where

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M18">View MathML</a>

(13)

Such that Si−1(t) is the residual echo with the echo of (i − 1)th scatterer estimated and removed. For i = 1, S0(t) = S(t) and Si(t ; θ) denote the radar echo of ith scatterer. To solve Equation (12), set the derivative of <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M19">View MathML</a> with respect to the parameter <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M20">View MathML</a> to zero, i.e.,

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M21">View MathML</a>

(14)

Here, Si*(t; θ) denotes the conjugate of Si (t; θ), then the estimation of the reflection coefficient of that scatterer can be given by

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M22">View MathML</a>

(15)

Following the CLEAN approach, after estimating each scatterer’s reflection coefficient, the current scatterer is subtracted from the frame via (13), and (12) is re-initialized with the location of the next strongest peak. This procedure is performed until the reflection coefficient of all the scatterers has been estimated. In this simulation, the estimated reflection coefficients of the four scatterers are <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M23">View MathML</a>. In order to measure the accuracy of this algorithm, we define maximum relative error and mean square error of estimated scattering coefficient as Equations (16) and (17).

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M24">View MathML</a>

(16)

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M25">View MathML</a>

(17)

where N represents the number of scatterers, here, N = 4, σ represents the real reflection coefficient of the scatterers, <a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M26">View MathML</a> denotes the estimated value.

From Equations (16) and (17), it can be calculated that the maximum relative error is 1.98%, and the mean quare error is 2.48 × 10–4, this shows that the algorithm can be used to extract the parameters of scattering coefficient well.

Establish the detection model of target after extracting all the parameters

<a onClick="popup('http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://jwcn.eurasipjournals.com/content/2013/1/77/mathml/M27">View MathML</a>

(18)

all the parameters in the equation are given in Table 1.

Table 1. The estimated parameters of all scatterers

If j in Equation (4) is j = 1, use the above method to estimate the parameters and establish the detection model of target.

The flow chart of this algorithm is shown in Figure 9.

thumbnailFigure 9. The flow chart of this algorithm.

4. Conclusion

This article establishes a novel echo model of target with micro-motion to analyze the characteristics of micro motion and investigates methods for motion parameter estimation and micro-Doppler signature extraction from target. Estimation of micro-motion parameters is completed through time-frequency transformation of the echoed signal and Radon transformation in terahertz band, and NLS and the CLEAN algorithm are utilized to estimate the scattering coefficients of each scatterer. This simulation result proves that Radon transformation detection method has high precision and good anti-noise performance which can accurately extract the micro-parameters. By adopting this estimate method, exact parameters are obtained for given signals. Thus, greatly precise the steps of target detection and identification.

Competing interests

The authors declare that they have no competing interests.

Acknowledgement

This study was supported by the National Natural Science Foundation of China under Projects 61271287 and the Fundamental Research Funds for the Central Universities under Projects ZYGX2012J029.

References

  1. VC Chen, FY Li, SS Ho, Micro-Doppler effect in Radar phenomenon, model and simulation study. IEEE Trans AES 42(1), 2–21 (2006)

  2. OR Fogle, Micro-range-micro-Doppler feature extraction and association. IEEE Radar Conference (Kansas City, MO, 2011), pp. 167–171 2011 OpenURL

  3. J Li, Y Pi, Research on terahertz radar target detection algorithm based on the extraction of micro motion feature. J Electron Meas Instrum 24(9), 803–807 (2010). Publisher Full Text OpenURL

  4. VC Chen, H Ling, Time-Frequency Transform for Radar Imaging and Signal Analysis (Artech House, Boston, 2002)

  5. B Boashash, Time-Frequency Signal Analysis and Processing (Elsevier Ltd., Amsterdam, 2003)

  6. MI Sari, Detection and reconstruction of geometric object based on line detection using radon transform. International Conference on Electrical Engineering and Informatics (ICEEI) 2011, 17–19 July 2011 (Bandung, 2011), pp. 1–6

  7. Q Wang, M Xing, High-resolution three-dimensional radar imaging for rapidly spinning targets, geoscience and remote sensing. IEEE Trans Geosci Remote Sens 1, 22–30 (2008)

  8. JJ Lei, C Lu, Target classification based on micro-Doppler signature. Proceedings of International Conference on Radar (Washington, USA, 2005), p. 179183 2005 OpenURL

  9. C Hongyong, L Yongxiang, Analysis of micro-Doppler and parameters estimation. J Infrared Millim Waves 25(5), 360–363 (2006)

  10. S Hui-xia, L Zheng, Micro-doppler feature extraction for ballistic missile warhead. Information and Automation, 2008 (Changsha, 2008), pp. 1333–1336