Abstract
MicroDoppler motion of a target is an important characteristic in highresolution 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 micromotion. Combination of timefrequency 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 antinoise 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:
MicroDoppler; Timefrequency analysis; Feature extraction; Radon transformation; NLS algorithm1. 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 microDoppler phenomenon. MicroDoppler phenomenon is very common in nature, such as the human heartbeat, vibration, or spin of missile warheads, etc.… while in terahertz band microDoppler phenomenon is particularly significant. The microDoppler effect enables us to determine the dynamic properties of the target and it offers a new approach for the analysis of target signatures. MicroDoppler features serve as additional target features that are complementary to those made available by existing methods. The microDoppler 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 micromotion of target in the terahertz band, the echo model is established [2], timefrequency transformation and Radon transformation are applied to extract the micro characteristic parameters [36], 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 micromotion to do a compound movement on the radar radial direction. In general, targets always do nonuniform linear motion in radical direction, which can be deemed as highorder 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
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
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?
Then the echo of target can be expressed as
where σ_{i} is the scattering coefficient of the ith scatterer, f_{0} 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 , then the rule of target’s motion can be approximated as
and the phase of the baseband signal is
where λ denotes radar wavelength, according to the definition of microDoppler, the microDoppler frequency of target equal to the derivative of phase with respect to time, then the instantaneous frequency of the signal can be expressed as
By using timefrequency transformation with respect to the echo of target, one can get its microDoppler frequency. Figure 1 is obtained in the simulation conditions as follows: radar carrier frequency f_{0} = 340 GHz, sampling frequency f_{s} = 4096 Hz, sampling points N = 512, and observation time t = 6 s. Assume that there are four scatterers, and v_{0} = 0.03 m/s^{2}, a = 0.01 m/s^{2}, σ_{1} = 1, σ_{2} = 0.9, σ_{3} = 1, σ_{4} = 1, A_{1} = B_{1} = φ_{1} = 0, A_{2} = A_{4} = 0.1, A_{4} = 0.2, B_{2} = 3 rad/s, B_{3} = 2 rad/s, B_{4} = 1 rad/s, φ_{2} = φ_{3} = φ_{4} = 0.6 rad.
Figure 1. MicroDoppler modulation induced by target.
3. Parameter estimation
3.1. Microfeature extraction based on Radon transformation
Previously, we established the echo model of target with micromotion, and timefrequency transformation of the echoed signal is applied to gain the timevarying microDoppler frequency features of target. This section is to estimate the motion parameters of target with micromotion, namely to extract the curves of target’s motion from the timevarying microDoppler frequency image.
We assume that the set of parameters to be estimated is θ = v_{0}, a, A, B, φ, as the straight line in Figure 1 is only related to parameters (v_{0}, a), we can divide it into a twodimensional matrix and threedimensional matrix to estimate, respectively, which can greatly reduce the calculation. First we will estimate (v_{0}, a) as shown in Figure 1.
Straight line detection can be realized by Radon transformation. A projection of a twodimensional 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 twodimensional space, Radon transformation can be defined by Equation (8)
where D is the entire image plane of x–y, 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 xaxis, δ 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.
Figure 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 = 2v_{0}/λ into consideration we have . 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?
Now timefrequency transformation of the echo signal is carried out in the same simulation condition as Figure 1, as shown in Figure 3.
Figure 3. MicroDoppler modulation induced by target.
There are a few methods for microDoppler feature extraction such as W–V peak detection method [9] and the normalized first moment method [10] that are all based on timefrequency 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 microDoppler frequency is after compensating the phase, and the amplitude of the sinusoidal curve in Figure 3 is , the vibration frequency is B. Now we let
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 threedimensional 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 of the timevarying microDoppler 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 timevarying microDoppler 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.
Figure 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.
Figure 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 A_{1} = x_{1} = 0, A_{2} = x_{2} = 0.1, A_{3} = x_{3} = 0.2, A_{4} = x_{4} = 0.1, B_{1} = y_{1} = 0, B_{2} = y_{2} = 3, B_{3} = y_{3} = 2, B_{4} = y_{4} = 0, C_{1} = z_{1} = 0, C_{2} = C_{3} = C_{4} = 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 antinoise. Its detection accuracy is determined by the size of the space grid of the threedimensional 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 W–V peak detection method and the normalized first moment method are shown in Figure 6.
Figure 6. Comparison of three methods of parameters estimation. (a) MicroDoppler modulation induced by target. (b) Radon transformation detection method. (c)W–V 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 microDoppler information very accurately.
As can be seen from Figure 7, while adding noise, both the W–V peak detection method and the normalized first moment method have been unable to accurately estimate the parameters of target microDoppler.
Figure 7. Comparison of three methods of parameters estimation with SNR = −6 dB. (a) MicroDoppler modulation induced by target with SNR = −6 dB. (b) Radon transformation detection method. (c)W–V 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 W–V peak detection method is close to the target microDoppler 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.
Figure 8. The simulation results of three methods with two scatterers. (a) MicroDoppler modulation induced by target with two scatterers. (b) Radon transformation detection method. (c)W–V 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
then the ith iteration of NLS may be expressed as
where
Such that S_{i−1}(t) is the residual echo with the echo of (i − 1)th scatterer estimated and removed. For i = 1, S_{0}(t) = S(t) and S_{i}(t ; θ) denote the radar echo of ith scatterer. To solve Equation (12), set the derivative of with respect to the parameter to zero, i.e.,
Here, S_{i}*(t; θ) denotes the conjugate of S_{i} (t; θ), then the estimation of the reflection coefficient of that scatterer can be given by
Following the CLEAN approach, after estimating each scatterer’s reflection coefficient, the current scatterer is subtracted from the frame via (13), and (12) is reinitialized 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 . 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).
where N represents the number of scatterers, here, N = 4, σ represents the real reflection coefficient of the scatterers, 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
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.
Figure 9. The flow chart of this algorithm.
4. Conclusion
This article establishes a novel echo model of target with micromotion to analyze the characteristics of micro motion and investigates methods for motion parameter estimation and microDoppler signature extraction from target. Estimation of micromotion parameters is completed through timefrequency 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 antinoise performance which can accurately extract the microparameters. 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

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

OR Fogle, MicrorangemicroDoppler feature extraction and association. IEEE Radar Conference (Kansas City, MO, 2011), pp. 167–171 2011

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

VC Chen, H Ling, TimeFrequency Transform for Radar Imaging and Signal Analysis (Artech House, Boston, 2002)

B Boashash, TimeFrequency Signal Analysis and Processing (Elsevier Ltd., Amsterdam, 2003)

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

Q Wang, M Xing, Highresolution threedimensional radar imaging for rapidly spinning targets, geoscience and remote sensing. IEEE Trans Geosci Remote Sens 1, 22–30 (2008)

JJ Lei, C Lu, Target classification based on microDoppler signature. Proceedings of International Conference on Radar (Washington, USA, 2005), p. 179183 2005

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

S Huixia, L Zheng, Microdoppler feature extraction for ballistic missile warhead. Information and Automation, 2008 (Changsha, 2008), pp. 1333–1336