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
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 . 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 , 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 .
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
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, f0 is the radar carrier frequency, c denotes the speed of light. Here, we only consider the condition of j = 1 and j = 2.
and the phase of the baseband signal is
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
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.
Figure 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)
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 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.
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 = 2v0/λ 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?
Figure 3. Micro-Doppler modulation induced by target.
There are a few methods for micro-Doppler feature extraction such as W–V peak detection method  and the normalized first moment method  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 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 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.
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 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 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 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) Micro-Doppler 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 micro-Doppler 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 micro-Doppler.
Figure 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)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 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.
Figure 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)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
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 with respect to the parameter to zero, i.e.,
Here, Si*(t; θ) denotes the conjugate of Si (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 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 . 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).
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.
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.
The authors declare that they have no competing interests.
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.
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
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