Micromotion feature extraction of radar target using tracking pulses with adaptive pulse repetition frequency adjustment

Yijun Chen; Qun Zhang; Changzheng Ma; Ying Luo; Tat Soon Yeo

doi:10.1117/1.JRS.8.083569

12 August 2014 Micromotion feature extraction of radar target using tracking pulses with adaptive pulse repetition frequency adjustment

Yijun Chen, Qun Zhang, Changzheng Ma, Ying Luo, Tat Soon Yeo

Author Affiliations +

Journal of Applied Remote Sensing, Vol. 8, Issue 1, 083569 (August 2014). https://doi.org/10.1117/1.JRS.8.083569

Abstract

In multifunction phased array radar systems, different activities (e.g., tracking, searching, imaging, feature extraction, recognition, etc.) would need to be performed simultaneously. To relieve the conflict of the radar resource distribution, a micromotion feature extraction method using tracking pulses with adaptive pulse repetition frequencies (PRFs) is proposed in this paper. In this method, the idea of a varying PRF is utilized to solve the frequency-domain aliasing problem of the micro-Doppler signal. With appropriate atom set construction, the micromotion feature can be extracted and the image of the target can be obtained based on the Orthogonal Matching Pursuit algorithm. In our algorithm, the micromotion feature of a radar target is extracted from the tracking pulses and the quality of the constructed image is fed back into the radar system to adaptively adjust the PRF of the tracking pulses. Finally, simulation results illustrate the effectiveness of the proposed method.

1. Introduction

In recent years, space targets (e.g., rocket bodies, fragments of satellites, ballistic warheads, etc.), especially space debris, are fast becoming a threat to aerospace activities.¹^,² Generally, a space target has a complex micromotion, such as spinning, precessing, and rolling, in addition to the body translation.³^,⁴ These micromotions will produce additional frequency modulation on the radar echo signal, which is called the micro-Doppler effect.⁵^–⁷ The micro-Doppler signature describes the characteristics of the micromotion target and reflects the transient behavior of the Doppler frequency shift, which provides important information for space target recognition.⁸ Therefore, since the concept of micro-Doppler was proposed, micromotion feature extraction has attracted much attention in research.⁹^–¹² Micro-Doppler feature extraction of the micromotion is usually implemented based on a time-frequency analysis technique. Based on this, an estimation method of the motion parameters was proposed in Ref. 13, and an adaptive chirplet representation for feature extraction from targets with rotating parts was analyzed in Ref. 14. A feature extraction based on empirical-mode decomposition (EMD) was proposed in Refs. 15 and 16. Also, image processing algorithms, such as the Radon transform¹⁷ and Hough transform,¹⁷^–¹⁹ have been introduced for the micromotion feature extraction. However, these methods need a long observation time and high pulse repetition frequency (PRF), because the time-frequency distribution (TFD)-Hough requires at least one-half or one full rotating period to successfully extract the micromotion feature.¹⁸ For example, assuming that the rotating radius is 5 m, the rotating frequency is 3 Hz, and the carrier frequency of the radar is 1 GHz, then the PRF must be higher than 1260 Hz to avoid the frequency-domain aliasing phenomenon of the micro-Doppler signal.

However, a phased array radar has to simultaneously serve different activities and the long-observation time and high PRF needed would place a great demand on the radar resources. On the other hand, the micromotion feature extraction is independent of the signal bandwidth. Hence, if the micromotion feature extraction of the radar target can be implemented with a narrow-band tracking pulse, then there is no need to allocate additional fixed radar resources for this task. Therefore, it will lighten the demand on the radar resources and the efficiency of the radar can be significantly improved. This is the motivation of this paper.

The radar target tracking period is usually relatively long. As such, the bandwidth of the micro-Doppler signal is relatively wide, which will create frequency-domain aliasing and make micromotion feature extraction difficult. To solve the problem, a varying PRF method²⁰ can be utilized to eliminate the aliasing. Broadly speaking, although the PRF for the tracking pulse is relatively low, micromotion feature extraction can be implemented through changing the PRF of the tracking pulse. However, in the case of a varying PRF, existing micromotion feature extraction algorithms, such as the time-frequency techniques, the Hough transform, the EMD algorithm, etc., will not function properly. In this paper, an orthogonal matching pursuit (OMP) algorithm is used to complete the micromotion feature extraction of targets via constructing a corresponding atom set of the micro-Doppler signals. Furthermore, since the rotating radius and initial phase of a micromotion scatterers can be estimated by the OMP algorithm, we can obtain the two-dimensional (2-D) distributed image of the target scatterers.

Because the rotating radius and rotating frequency are not known a priori, it is not possible to develop an algorithm to determine the optimal PRF. Thus, the simplest way to find the most suitable PRF is by first estimating the dominant micro-Doppler frequencies of all micromotions and using them to formulate a scheme for a varying PRF. However, this action itself is also a waste of the radar resources. Therefore, an adaptively varying PRF method is proposed in this paper to further improve radar efficiency. In this method, through constructing the atoms set at the prevailing tracking data rate, the target image quality is assessed and fed back to the radar system to guide the adjustment of the PRF of the tracking pulses and the reconstruction of the dictionary atoms.

This paper is organized as follows: The analysis of the radar echo of micromotion target and the micromotion feature extraction with varying PRF sampling based on OMP are illustrated in Sec. 2. The micromotion feature extraction method using tracking pulses with an adaptive PRF is proposed in Sec. 3. Simulations are presented in Sec. 4 to validate the effectiveness of the proposed method, and some conclusions are made in the last section.

2. Micromotion Feature Extraction of Radar Target with Varying PRF Sampling

2.1.

Analysis of the Radar Echo of Micromotion Target

Taking a spinning target as an example, the analysis of the narrow-band radar echo is as follows: Assume that the translational motion of a spinning target has been compensated and the geometry of the radar and the spinning target is shown in Fig. 1, where the LOS represents the line-of-sight direction of the radar, the target is rotating around the $Z$ -axis at an angular velocity $\vec{ω}$ , and $α$ denotes the intersection angle between the LOS and $\vec{ω}$ . In fact, the rotation velocity $\vec{ω}$ can be divided into ${\vec{ω}}_{e}$ and ${\vec{ω}}_{R}$ , where ${\vec{ω}}_{e}$ is perpendicular to the LOS and ${\vec{ω}}_{R}$ is parallel to the LOS. Obviously, the rotation produced by ${\vec{ω}}_{R}$ will not produce any radial motion and will not produce micro-Doppler modulation on the echo signal. On the other hand, ${\vec{ω}}_{e}$ will produce radial motion and also micro-Doppler modulation. Therefore, ${\vec{ω}}_{e}$ is called the effective rotation vector.²¹ The effective imaging plane is perpendicular to ${\vec{ω}}_{e}$ , and $P^{'}$ is the projection of the target scatterer $P$ in the imaging plane.

Fig. 1

Geometry of the spinning target.

Assume that the signal transmitted by the radar is $p (t) = \exp (j 2 π f_{c} t)$ , then the baseband echo signal from the scatterer $p$ on the target is

Eq. (1)

s_{p} (τ) = σ_{p} \exp (- j \frac{4 π f_{c} R_{p} (τ)}{c}) τ \in [0, T_{c}],

where

τ

is the slow time,

f_{c}

denotes the carrier frequency, and

σ_{p}

is the reflection coefficient of the

p

th scatterer.

R_{p} (t)

is the instantaneous distance between the

p

th scatterer and the radar at slow time

τ

. Under the far-field condition and based on the plane wave approximation, the instantaneous distance

R_{p} (t)

can be expressed as

Eq. (2)

R_{p} (τ) = r_{p} \sin ((ω_{p} + f_{t r}) τ + θ_{p}) \sin (α) = r_{p}^{'} \sin ((ω_{p} + f_{t r}) τ + θ_{p}),

where the

(r_{p}, θ_{p})

is the polar coordinate of the

p

th scatterer in the rotating plane,

ω_{p} = | {\vec{ω}}_{p} |

denotes the length of the target rotation angular vector, namely, the rotation frequency around the

Z

-axis,

r_{p}^{'}

is the effective rotation radius in the imaging plane, and

f_{tr}

is the relative rotating speed caused by the target translational movement. Generally, the angular variation caused by target spinning is far greater than the relative rotation angle caused by

f_{tr}

. For instance, assume the distance from the target to the radar is 1000 km, the horizontal flight speed is 5000 m/s, and the relative rotating speed caused by the translational movement is approximately only 0.005 rad/s. Under these conditions,

f_{tr}

can be ignored. Assume that the target consists of

P

scatterers, then the baseband echo signals from the target can be represented as

Eq. (3)

s (τ) = \sum_{p = 1}^{P} σ_{p} \exp (- j \frac{4 π f_{c}}{c} \cdot r_{p}^{'} \sin (ω_{p} τ + θ_{p})) .

For convenience, Eq. (3) can be written in the discrete form as follows:

Eq. (4)

s (m) = \sum_{p = 1}^{P} σ_{p} \exp (- j \frac{4 π f_{c}}{c} \cdot r_{p}^{'} \sin (ω_{p} \cdot m \cdot PRI + θ_{p})) m = 1,2, \dots, N,

where

N = PRF \cdot T_{c}

,

PRI = 1 / PRF

, indicating the pulse repetition interval.

If the radar transmits the pulses with varying PRF, namely the radar transmits the pulses to the target with ${PRF}_{n}$ within the interval $[t_{n}, t_{n + 1}] (n = 1,2, \dots, N, t_{1} = 0, t_{N + 1} = T_{c})$ , then the baseband echo signals of the target can be written as

Eq. (5)

s_{s} (m) {\begin{cases} = \sum_{p = 1}^{P} σ_{p} \exp [- j \frac{4 π f_{c}}{c} \cdot r_{p}^{'} \sin (ω_{p} \cdot m \cdot {PRI}_{1} + θ_{p})] & m = 1, \dots, M_{1} \\ = \sum_{p = 1}^{P} σ_{p} \exp [- j \frac{4 π f_{c}}{c} \cdot r_{p}^{'} \sin (ω_{p} \cdot M_{1} \cdot {PRI}_{1} + ω_{p} \cdot (m - M_{1}) \cdot {PRI}_{2} + θ_{p})] & m = M_{1} + 1, \dots, M_{1} + M_{2} \\ ⋮ \\ = \sum_{p = 1}^{P} σ_{p} \exp [- j \frac{4 π f_{c}}{c} \cdot r_{p}^{'} \sin (ω_{p} \sum_{n = 1}^{N - 1} M_{i} \cdot {PRI}_{n} + ω_{p} \cdot (m - M_{all} + M_{N}) \cdot {PRI}_{N} + θ_{p})] & m = M_{all} - M_{N} + 1, \dots, M_{all} \end{cases},

where

M_{n} = (t_{n + 1} - t_{n}) \cdot {PRF}_{n}

and

M_{all} = M_{1} + M_{2} + \dots + M_{N}

.

2.2.

Micromotion Feature Extraction Method Based on OMP

In recent years, the advantages of high efficiency and easy implementation have seen the wide use of the OMP algorithm for sparse signal decomposition and reconstruction. From the principle of the OMP algorithm, we know that the decomposition and reconstruction of the signals can be implemented by constructing the atoms in the dictionary according to the intrinsic characteristics of the signal itself.²²^,²³

It should be pointed out that when extracting the micromotion feature of a target, we need to satisfy $PRF > 8 π f R_{\max} / λ$ to avoid the frequency-domain aliasing of the micro-Doppler signal,¹⁶ where $f$ represents the target rotation frequency, $R_{\max}$ denotes the maximum rotation radius of the target, and $λ$ is the wavelength. Obviously, the low PRF of tracking pulses usually cannot satisfy this requirement, which leads to under-sampling of the micro-Doppler signal and produces frequency-domain aliasing. Thus, the effective results of micromotion feature extraction cannot be obtained based on the OMP for a low-PRF tracking radar. However, varying the PRF sampling can enlarge the nonambiguous interval of the micro-Doppler signal.

When observing the target with the varying PRF as described in Eq. (5), the peak of the frequency spectrum of the signal with frequency $F$ is $f_{n} (n = 1,2, \dots, N)$ with ${PRF}_{n} (n = 1,2, \dots, N)$ sampling and the nonambiguous interval is ${PRF}_{n} (n = 1,2, \dots, N)$ . Only when $F$ satisfies $f_{n} = F \mod ({PRF}_{n}) (n = 1,2, \dots, N)$ , can it be mapped to the corresponding atoms in the dictionary. At this time, the nonambiguous interval of the signal is the least common multiple of ${PRF}_{n} (n = 1, \dots N)$ , which significantly enlarges the nonambiguous interval in the micro-Doppler frequency domain. As a result, the micromotion feature extraction of the radar target can be realized based on the OMP. It should be pointed out that ${PRF}_{n} (n = 1, \dots N)$ should be co-prime to make the nonfuzzy interval of the signal as large as possible.

From Eq. (5), it can be seen that the echo of each micromotion point is determined by the parameter $(r_{p}^{'}, ω_{p}, θ_{p})$ . Therefore, the atom can be constructed by setting different $(r^{'}, ω, θ)$ values

Eq. (6)

\vec{d} (r^{'}, ω, θ) = \sum_{n = 1}^{N} \exp [- j \frac{4 π f_{c}}{c} \cdot r^{'} \frac{R_{\max}^{'} - R_{\min}^{'}}{N_{r^{'}}} \cdot \sin (ω \frac{Ω_{\max} - Ω_{\min}}{N_{ω}} \cdot (\sum_{i = 0}^{n - 1} M_{i} \cdot {PRI}_{i} + (U (m - M_{0} - \dots - M_{n - 1} - 1) - U (m - M_{0} - \dots - M_{n} - 1)) \cdot (m - M_{0} - \dots - M_{n - 1}) \cdot {PRI}_{n}) + θ \frac{Θ_{\max} - Θ_{\min}}{N_{θ}})],

where

M_{0} = 0

,

{PRI}_{0} = 0

, and

U (x) = {\begin{matrix} 1 x \geq 0 \\ 0 x < 0 \end{matrix}

;

[R_{\min}^{'}, R_{\max}^{'}], [Ω_{\min}, Ω_{\max}], [Θ_{\min}, Θ_{\max}]

are the value ranges of

r^{'}

,

ω

, and

θ

, respectively, and

(R_{\max}^{'} - R_{\min}^{'}) / N_{r^{'}}

,

(Ω_{\max} - Ω_{\min}) / N_{ω}

, and

(Θ_{\max} - Θ_{\min}) / N_{θ}

are the increment steps of the searches for

r^{'}

,

ω

, and

θ

, respectively.

Applying normalized processing to each atom: $\vec{d} (r^{'}, ω, θ) = \vec{d} (r^{'}, ω, θ) / {‖ \vec{d} (r^{'}, ω, θ) ‖}_{2}$ , then the atom set can be represented as

Eq. (7)

D = [{\vec{d} (1,1, 1), \dots, \vec{d} (1,1, N_{θ}), \vec{d} (1,2, 1), \dots, \vec{d} (1,2, N_{θ}), \dots, \vec{d} (1, N_{ω}, N_{θ}), \vec{d} (2,1, 1), \dots, \vec{d} (2, N_{ω}, N_{θ}), \dots, \vec{d} (N_{r^{'}}, N_{ω}, N_{θ})]}_{M \times N_{r^{'}} N_{ω} N_{θ}} .

For the sake of simplicity, $D$ is recorded as $D = {[{\vec{d}}_{1}, {\vec{d}}_{2}, \dots, {\vec{d}}_{N_{r^{'}} N_{ω} N_{θ}}]}_{M \times N_{r^{'}} N_{ω} N_{θ}}$ . The signal $s_{s}$ with a varying PRF sampling can be denoted as

Eq. (8)

s_{s} = D β,

where

β

is the projection coefficient of

s_{s}

under the atom set

D

.

The detailed steps of the micromotion feature extraction algorithm with varying PRF sampling based on the OMP can be described as follows:

Step 1: Initialize the parameters: the residual $s_{r_{0}} = s_{s}$ , the index set ${pos}_{0} = \emptyset$ , the matched atoms recorder matrix $Π_{0} = \emptyset$ , the power threshold $δ > 0$ , the iteration counter $h = 1$ , the maximum iteration number $H$ , and $β$ is the zero vector with $N_{r^{'}} N_{ω} N_{θ} \times 1$ dimensions.
Step 2: Calculate the inner product between $s_{r_{h - 1}}$ and all atoms in the atom set $D$ : ${< s_{r_{h - 1}}, d_{l} > | l = 1,2, \dots, N_{r^{'}} N_{ω} N_{θ}}$ .
Step 3: Find the number of the atom in the atom set corresponding to the maximum inner product: ${pos}_{h} = \arg \max_{l = 1,2, \dots N_{r^{'}} N_{ω} N_{θ}} | < s_{r_{h - 1}}, d_{l} > |$ ; then set ${pos}_{h} = {pos}_{h - 1} \cup {{pos}_{h}}$ .
Step 4: Record the atom corresponding to the maximum inner product in $Π$ , $Π_{h} = [Π_{h - 1}, d_{{pos}_{h}}]$ , and remove it from the atom set $D$ .
Step 5: Solve the optimization problem $\hat{β} = \arg \min_{β} ‖ s_{s} - Π_{h} β ‖_{2}$ with the least squares method: $\hat{β} = {(Π_{h}^{H} Π_{h})}^{- 1} Π_{h}^{H} s_{s}$ .
Step 6: Update the residual: $s_{r_{h}} = s_{s} - Π_{h} \hat{β}$ .
Step 7: Set $h = h + 1$ , if $h < H$ and ${‖ s_{r_{h}} ‖}_{2} > δ$ , repeat Steps 2 to 6; if $h = H$ or ${‖ s_{r_{h}} ‖}_{2} \leq δ$ , iteration ends, proceeds to Step 8.
Step 8: Extract the micromotion feature of the target: the micromotion parameter $(r^{'}, ω, θ)$ of the micromotion point can be extracted according to the atom number in the index set pos. The detailed steps can be described as follows:
First, transform the atom number in ${pos}_{H}$ into the number of the parameter $r^{'}, ω, θ$ :
Eq. (9)
$index_r^{'} (i) = ⌈ p o s (i) / N_{ω} N_{θ} ⌉, i = 1, \dots, H$
Eq. (10)
$index_ω (i) = ⌈ (pos (i) - (index_r^{'} (i) - 1) \cdot N_{ω} N_{θ}) / N_{θ} ⌉, i = 1, \dots, H$
Eq. (11)
$index_θ (i) = pos (i) - (index_r^{'} (i) - 1) \cdot N_{ω} N_{θ} - (index_ω (i) - 1) \cdot N_{θ}, i = 1, \dots, H .$
Then, the micromotion feature of the micromotion point can be obtained
Eq. (12)
$r^{'} (i) = \frac{R_{\max}^{'} - R_{\min}^{'}}{N_{r^{'}}} \cdot index_r^{'} (i)$
Eq. (13)
$ω (i) = \frac{Ω_{\max} - Ω_{\min}}{N_{ω}} \cdot index_ω (i)$
Eq. (14)
$θ (i) = \frac{Θ_{\max} - Θ_{\min}}{N_{θ}} \cdot index_θ (i) .$
Step 9: Because $(r^{'}, θ)$ is the polar coordinate of a scatterer in the imaging plane, the 2-D distributed information of the target points can be determined according to the micromotion feature parameter $r^{'} (i)$ and $θ (i)$ extracted in Step 8, and $\hat{β} (i)$ is the reflection coefficient of the point, thus we get the image of the target.

3. Micromotion Feature Extraction of Radar Target from Tracking Pulse with Adaptive PRF Adjustment

Since the tracking pulses are usually narrow-band, the micromotion feature extraction of the radar target can be implemented from tracking pulses with varying PRF samplings based on the OMP as mentioned in Sec. 2. However, without knowing the rotating radius and rotating frequency a priori, it is difficult to judge if the ${PRF}_{n}$ ( $n = 1, 2, \dots, N$ ) used can effectively eliminate the frequency-domain aliasing of the micro-Doppler signal. Therefore, we need to feed the feature extracted from the micromotion and the target image quality back to the radar system and adaptively guide the adjustment of the PRF and the dictionary atoms.

Set the tracking data rate set ${PRF}_{1} < {PRF}_{2} < \dots < {PRF}_{N}$ and ensure that ${PRF}_{n} (n = 1, \dots N)$ are co-prime. Set the initial target tracking data rate $f_{t, 1} = {PRF}_{1}$ , extract the micromotion feature of the target based on the OMP algorithm at the periodic interval $T$ , and select the correlation coefficient of the target images $F_{n - 1} (i, j)$ and $F_{n} (i, j)$ as an image quality assessment metric, which is defined as

Eq. (15)

α = \frac{\sum_{i} \sum_{j} F_{n - 1} (i, j) ⊙ F_{n} (i, j)}{\sqrt{‖ F_{n - 1} (i, j) ‖_{2} ‖ F_{n} (i, j) ‖_{2}}},

where

⊙

represents the Hadamard product, and

{‖ \cdot ‖}_{2}

denotes the

l_{2}

-norm for the matrix.

Obviously, if the correlation coefficient is relatively small, this means the similarity of target images from adjacent intervals based on the OMP algorithm is relatively low and the image is not identical to the distribution of the scattering centers of the target, thus the PRF needs to be further increased. On the contrary, if the correlation coefficient is relatively large, this indicates that the images are almost identical and the micromotion feature extraction of radar target has been effectively achieved.

Therefore, we should select an appropriate threshold $T_{α}$ . If $α$ is less than the threshold, then it means that the required precision of the micromotion feature extraction has not been attained and we should increase the tracking data rate in the next interval $T$ : $f_{t, n + 1} = {PRF}_{n + 1}$ ; Otherwise, we have the required precision for the micromotion feature extraction and we shall maintain the target tracking data rate $f_{t, n + 1} = {PRF}_{n}$ .

It needs to be pointed out that the dictionary atoms of the OMP algorithm should be updated with a change of the target tracking data rate.

The atom in the dictionary during the $n$ ’th micromotion feature extraction of the target is recorded as ${\vec{d}}_{n} (r^{'}, ω, θ)$ . If the correlation coefficient of the target image obtained from this time interval and the previous time interval does not exceed the threshold $T_{α}$ , then we need to increase the tracking data rate in the next interval: $f_{t, n + 1} = {PRF}_{n + 1}$ . Then, the updated equation of the atom in the dictionary can be written as

Eq. (16)

{\vec{d}}_{n + 1}^{'} (r^{'}, ω, θ) = {\vec{d}}_{n} (r^{'}, ω, θ) + \exp [- j \frac{4 π f_{c}}{c} \cdot r^{'} \frac{R_{\max}^{'} - R_{\min}^{'}}{N_{r^{'}}} \cdot \sin (ω \frac{Ω_{\max} - Ω_{\min}}{N_{ω}} \cdot (\sum_{i = 0}^{n} M_{n} \cdot {PRI}_{n} + (U (m - M_{0} - \dots - M_{n} - 1) - U (m - M_{0} - \dots - M_{n + 1} - 1)) \cdot (m - M_{0} - \dots - M_{n}) \cdot {PRI}_{n + 1}) + θ \frac{Θ_{\max} - Θ_{\min}}{N_{θ}})],

where

M_{n + 1} = T \cdot {PRF}_{n + 1}

.

In conclusion, the flowchart of the micromotion feature extraction method from tracking pulse with the adaptive PRF proposed in this paper is shown in Fig. 2.

Fig. 2

Micromotion feature extraction of radar target from tracking pulse with adaptive pulse repetition frequency (PRF) adjustment.

What should be pointed out is that with the condition of high-frequency scattering, the shielding effect can be easily induced by the spin motion because of the geometry of the target. In this case, a nonlinear shielding function $F_{p}$ needs to be introduced to describe the target scattering state, then the reflection coefficient of the $p$ th scatterer is $σ_{p} \cdot F_{p}$ . Without the loss of generality, assuming $F_{p}$ obeys the Bernoulli distribution

Eq. (17)

P {F_{p} (τ) = 1} = 1 - ξ P {F_{p} (τ) = 0} = ξ,

where

ξ

is the shielding rate and denotes the possibility of scattering shielding. However, the shielding effect does not change the intrinsic characteristics of the radar echo signal, so the atom set

D

is still reasonable and the proposed micromotion feature extraction method from the tracking pulse can also successfully extract the micromotion features of the target.

4. Experimental Results

To validate the effectiveness of the micromotion feature extraction method with varying PRF sampling based on the OMP, simulations are given in this section. The pulse duration is $T_{p} = 1 μ s$ and the carrier frequency is $f_{c} = 1 GHz$ . The radar is located at the origin $O$ of the coordinate system, and the coordinate of the target reference point is (0, 1000, 1700) km. The target consists of three rotational scatterers whose coordinates are $(10.4, - 6,0)$ , $(- 10.4, - 6,0)$ , and $(0,10, 0)$ , respectively, in meters relative to the reference point in the rotating plane. The flight speed of the target is $(0,0,1000) m / s$ , the target is rotating around the $Z$ -axis, and its rotation velocity is $(0,0, π) rad / s$ . The distribution of the target scatterer in the imaging plane is shown in Fig. 3.

Fig. 3

Geometry of the target.

The observation duration is $T_{c} = 2 s$ . The signals are sampled with a $PRF = 13 Hz$ at the first 0.6 s, sampled with a PRF of 15 Hz between 0.6 to 1.2 s, and sampled with a PRF of 19 Hz for the final 0.8 s. Obviously, each PRF does not satisfy $PRF > 8 π f R_{\max} / λ = 252$ , and the time-frequency analysis of the echo signal while varying the PRF sampling is as shown in Fig. 4.

Fig. 4

Time-frequency distribution with varying PRF sampling.

From Fig. 4, it can be seen that the micromotion feature of the target cannot be extracted due to aliasing. However, the micromotion feature extraction of the target can be implemented with the method proposed in Sec. 2 of this paper. Setting the maximum iteration count at $H = 4$ , the threshold of the residual signal energy is $δ = 0.1$ . The result of the micromotion feature extraction of the target is shown in Table 1. The target’s image can be obtained as shown in Fig. 5(a), and the time-frequency distribution of the reconstructed signal from the decomposition result is shown in Fig. 5(b).

Table 1

Micromotion features extracted with varying PRF sampling.

	$r^{'} (m)$	$ω (rad / s)$	$θ (rad)$	Coefficient $x$
Feature parameter 1	6.0	$π$	$11 π / 6$	23.58
Feature parameter 2	5.0	$π$	$π / 2$	21.37
Feature parameter 3	5.9	$π$	$7 π / 6$	21.01
Feature parameter 4	6.0	$π$	$7 π / 6$	3.42

Fig. 5

Image of the target reconstructed using Orthogonal Matching Pursuit (OMP) method and the time-frequency distribution of the reconstructed signal for noise-free case. (a) The target image by OMP. (b) The time-frequency distribution of the reconstructed signal from the decomposition result.

From Table 1, we can see that the algorithm can successfully extract the micromotion features of the three scatterers with $r^{'}$ be 5.0, 5.9, and 6.0 m, $ω$ be $π rad / s$ , and $θ$ be $π / 2 rad$ , $7 π / 6 rad$ , and $11 π / 6 rad$ , respectively, which are close to the original real parameters. It should be pointed out that four groups of feature parameters are extracted since the maximum iteration number $H$ in the OMP algorithm has been set as 4. Table 1 indicates that as the iteration increases, the energy of the extracted signal decreases correspondingly, and the fourth group of feature parameters corresponds to the residual signal with the lowest energy. The target image shown in Fig. 5(a) is almost identical to the scatterers shown in Fig. 3. The time-frequency distribution of the reconstructed signal from the decomposition result also reinforces the correctness of the results.

To validate the robustness of the algorithm, Gaussian noise with $SNR = 0 dB$ is added into the echo signal. The echo signal of the target is sampled with varying PRF, and then the micromotion feature extraction of the target is implemented with the proposed method. The result is shown in Table 2. The image reconstructed from the extracted parameter $(r^{'}, θ)$ and coefficient $\hat{β}$ is shown in Fig. 6(a), and the time-frequency distribution of the reconstructed signals from the decomposition result is shown in Fig. 6(b).

Table 2

Micromotion features extracted with varying PRF sampling (SNR=0 dB).

	$r^{'} (m)$	$ω (rad / s)$	$θ (rad)$	Coefficient $x$
Feature parameter 1	6.0	$π$	$11 π / 6$	24.16
Feature parameter 2	6.0	$π$	$7 π / 6$	21.09
Feature parameter 3	5.0	$π$	$π / 2$	19.87
Feature parameter 4	5.9	$π$	$11 π / 6$	3.13

Fig. 6

Image of the target reconstructed using OMP method and the time-frequency distribution of the reconstructed signal for $SNR = 0 dB$ . (a) The target image by OMP. (b) The time-frequency distribution of the reconstructed signal from the decomposition result.

From Table 2, we can see that the algorithm can successfully extract the micromotion features of the three scatterers with a low $SNR = 0 dB$ . The image shown in Fig. 6(a) is almost identical to the distribution of the target scatterers. However, when SNR is reduced to $- 5 dB$ , the micromotion features cannot be correctly extracted with the algorithm.

The previous simulations have been carried out with varying PRFs (0.6 s at 13 Hz, 0.6 s at 15 Hz, and 0.8 s at 19 Hz). Next, we carried out the micromotion feature extraction method with an adaptive PRF. We stepped the correlation of the images through the set of tracking data rate: ${PRF}_{1} = 5 Hz$ , ${PRF}_{2} = 7 Hz$ , ${PRF}_{3} = 9 Hz$ , ${PRF}_{4} = 11 Hz$ , ${PRF}_{5} = 13 Hz$ , ${PRF}_{6} = 17 Hz$ , and ${PRF}_{7} = 19 Hz$ . The interval of the micromotion feature extraction of the target based on the OMP algorithm is $T = 0.4 s$ , and the threshold of the correlation coefficient of the target images from the adjacent times is $T_{α} = 0.8$ .

The plot of the correlation coefficient of the images from each adjacent two time-intervals is shown in Fig. 7. We can see that the correlation coefficient reaches the threshold after the fifth micromotion feature extraction, and the target tracking data rate is $f_{t, 5} = 13 Hz$ , much smaller than the 19 Hz used in the previous simulations; here, the extracted micromotion feature is shown in Table 3 and the image reconstructed from the extracted parameter is shown in Fig. 8.

Fig. 7

Correlation coefficient of the target’s images.

Table 3

Micromotion features extracted from the tracking pulse with adaptive PRF.

	$r^{'} (m)$	$ω (rad / s)$	$θ (rad)$	Coefficient $x$
Feature parameter 1	6.0	$π$	$11 π / 6$	23.11
Feature parameter 2	5.9	$π$	$7 π / 6$	20.84
Feature parameter 3	5.0	$π$	$π / 2$	20.80
Feature parameter 4	5.1	$π$	$π / 2$	2.93

Fig. 8

The fifth target image by OMP.

To highlight the advantage of the proposed algorithm, the traditional micromotion feature extraction method based on TFD-Hough transform is used for comparison. Assume $PRF = 1000 Hz$ that satisfies $PRF > 8 π f R_{\max} / λ$ to avoid the spectrum-aliasing problem. The micromotion features extracted by the TFD-Hough transform are shown in Table 4. From Table 4, we can see that the precision of the TFD-Hough transform is not satisfactory because the method is limited by the time-frequency resolution. Also, the TFD-Hough method needs 2-s continuous time resource for target observation with $PRF = 1000 Hz$ . Comparatively, the proposed micromotion feature extraction method using tracking pulses can obtain more precise extraction results without continuous time-resource allocation. As a result, the radar efficiency can be improved.

Table 4

Micromotion features extracted based on TFD-Hough.

	$r^{'} (m)$	$ω (rad / s)$	$θ (rad)$
Feature parameter 1	4.4	$π$	0
Feature parameter 2	5.3	$π$	$3 π / 4$
Feature parameter 3	5.0	$π$	$π / 8$
Feature parameter 4	5.7	$π$	0
Feature parameter 5	4.9	$π$	$π / 12$
Feature parameter 6	5.5	$π$	$π / 4$
Feature parameter 7	7.2	$π$	$19 π / 12$
Feature parameter 8	6.2	$π$	$π / 2$

Finally, we discussed the shielding effect of scatterers. Assume $ξ = 0.5$ and $PRF = 1000 Hz$ , the time-frequency distribution of the radar echo signal is shown in Fig. 9. We can see that the shielding effect produces an adverse influence on the time-frequency distribution. In the final simulation, we implement the micromotion feature extraction method with an adaptive PRF under the shielding effect. The set of tracking data rate is the same as previous simulations, and the plot of the correlation coefficient of the images from each adjacent two time-intervals is shown in Fig. 10. We can see that the correlation coefficient reaches the threshold after the sixth micromotion feature extraction, and the target tracking data rate is $f_{t, 6} = 17 Hz$ . Here, the extracted micromotion feature is shown in Table 5 and the time-frequency distribution of the reconstructed signals from the decomposition result is shown in Fig. 11. The simulation shows the good performance of our proposed method under the shielding effect.

Fig. 9

Time-frequency distribution with shielding effect.

Fig. 10

Correlation coefficient of the target’s images.

Table 5

Micromotion features extracted from the tracking pulse with adaptive PRF under shielding effect.

	$r^{'} (m)$	$ω (rad / s)$	$θ (rad)$	Coefficient $β$
Feature parameter 1	6.0	$π$	$11 π / 6$	21.13
Feature parameter 2	5.9	$π$	$7 π / 6$	17.94
Feature parameter 3	5.0	$π$	$π / 2$	17.26
Feature parameter 4	5.1	$π$	$π / 2$	2.16

Fig. 11

The time-frequency distribution of the reconstructed signal from the decomposition result.

5. Conclusion

To lighten the demand on radar resources, a micromotion feature extraction method using tracking pulses with an adaptive PRF is proposed in this paper. With this method, the atom set of the micro-Doppler signals is constructed according to the prevailing tracking data rate, and the correlation of the images (constructed at the current and previous PRFs) is fed back to the radar system and used to adjust the PRF of the transmitted pulse and the dictionary atoms. The micromotion feature extraction of the target is then implemented using tracking pulses based on the OMP method. The steps of the method are described in detail, and the simulations are given to illustrate its effectiveness. Although the spinning target is taken as an example in this paper, the proposed method is also suitable for other micromotion targets such as vibration and precession targets.

In fact, recognition of a target can be achieved according to the result of the micromotion target feature extraction of the target. On this basis, the tracking data rate can be adjusted again according to the different types of the targets. Meanwhile, the tracking property can be fed back to the radar system to adjust the tracking period.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China under Grant 61201369 and 61172169 and in part by Natural Science Foundation Research Program of Shaanxi Province under Grant 2013JQ8008.

References

1.

X. Baiet al., “High-resolution three-dimensional imaging of spinning space debris,” IEEE Trans. Geosci. Remote Sens., 47 (7), 2352 –2362 (2009). http://dx.doi.org/10.1109/TGRS.2008.2010854 IGRSD2 0196-2892 Google Scholar

2.

Y. Luoet al., “Micro-motion feature extraction of target in inverse synthetic aperture radar imaging with sparse aperture,” J. Electromagn. Waves Appl., 27 (14), 1841 –1849 (2013). http://dx.doi.org/10.1080/09205071.2013.825890 JEWAE5 0920-5071 Google Scholar

3.

J. Heet al., “Micro-Doppler effect analysis and feature extraction in inverse synthetic aperture imaging LADAR imaging,” J. Appl. Remote Sens., 5 (1), 051502 (2011). http://dx.doi.org/10.1117/1.3652706 1931-3195 Google Scholar

4.

K. Y. Guoet al., “Influence of migratory scattering phenomenon on micro-motion characteristics contained in radar signals,” IET Radar Sonar Navig., 7 (5), 579 –589 (2013). http://dx.doi.org/10.1049/iet-rsn.2012.0058 1751-8784 Google Scholar

5.

Y. Luoet al., “Micro-Doppler effect analysis and feature extraction in ISAR imaging with stepped-frequency chirp signals,” IEEE Trans. Geosci. Remote Sens., 48 (4), 2087 –2098 (2010). http://dx.doi.org/10.1109/TGRS.2009.2034367 IGRSD2 0196-2892 Google Scholar

6.

V. C. Chenet al., “Micro-Doppler effect in radar: phenomenon, model and simulation study,” IEEE Trans. AES, 42 (1), 2 –21 (2006). http://dx.doi.org/10.1109/TAES.2006.1603402 IEARAX 0018-9251 Google Scholar

7.

C. Wanget al., “Micro-Doppler feature extraction using single-frequency radar for high-speed targets,” in Radar Conf. 2013, IET Int., 1 –6 (2013). Google Scholar

8.

Y. Renet al., “A space target recognition method based on compressive sensing,” in 2011 Sixth Int. Conf. on Image and Graphics (ICIG), 582 –586 (2011). Google Scholar

9.

T. ThayaparanL. StankovićI. Djurović, “Micro-Doppler-based target detection and feature extraction in indoor and outdoor environments,” J. Franklin Inst., 345 (6), 700 –722 (2008). http://dx.doi.org/10.1016/j.jfranklin.2008.01.003 JFINAB 0016-0032 Google Scholar

10.

H. SunZ. Liu, “Micro-Doppler feature extraction for ballistic missile warhead,” in Int. Conf. on Information and Automation, 1333 –1336 (2008). Google Scholar

11.

H. Gaoet al., “Micro-Doppler signature extraction from ballistic target with micro-motions,” IEEE Trans. AES, 46 (4), 1969 –1982 (2010). http://dx.doi.org/10.1109/TAES.2010.5595607 IEARAX 0018-9251 Google Scholar

12.

D. Zhanget al., “The study of three-dimensional micro-motion feature extraction of ballistic target based on netted radar,” Appl. Mech. Mater., 530–531 534 –539 (2014). http://dx.doi.org/10.4028/www.scientific.net/AMM.530-531 AMMPDO 1662-7482 Google Scholar

13.

T. Thayaparan, “Analysis of radar micro-Doppler signatures from experimental helicopter and human data,” IET Radar, Sonar & Navig., 1 (4), 289 –299 (2007). http://dx.doi.org/10.1049/iet-rsn:20060103 1751-8784 Google Scholar

14.

J. LiH. Ling, “Application of adaptive chirplet representation for ISAR feature extraction from targets with rotating parts,” IEE Proc. Radar Sonar Navig., 150 (4), 284 –291 (2003). http://dx.doi.org/10.1049/ip-rsn:20030729 IRSNE2 1350-2395 Google Scholar

15.

Y. LiL. DuH. Liu, “Hierarchical classification of moving vehicles based on empirical mode decomposition of micro-Doppler signatures,” IEEE Trans. Geosci. Remote Sens., 51 (5), 3001 –3013 (2013). http://dx.doi.org/10.1109/TGRS.2012.2216885 IGRSD2 0196-2892 Google Scholar

16.

L. Duet al., “Robust classification scheme for airplane targets with low resolution radar based on EMD-CLEAN feature extraction method,” Sens. J., 13 (12), 4648 –4662 (2013). http://dx.doi.org/10.1109/JSEN.2013.2272119 ISJEAZ 1530-437X Google Scholar

17.

X. Chenet al., “Maneuvering target detection via radon-fractional Fourier transform-based long-time coherent integration,” IEEE Trans. Signal Process., 62 (4), 939 –953 (2014). http://dx.doi.org/10.1109/TSP.2013.2297682 ITPRED 1053-587X Google Scholar

18.

Q. Zhanget al., “Imaging of a moving target with rotating parts based on the Hough transform,” IEEE Trans. Geosci. Remote Sens., 46 (1), 291 –299 (2008). http://dx.doi.org/10.1109/TGRS.2007.907105 IGRSD2 0196-2892 Google Scholar

19.

Y. Lianget al., “Analysis and extraction of micro-Doppler features in FMCW-ISAR using a modified extended Hough transform,” in 2011 IEEE CIE Int. Conf. on Radar, 1652 –1655 (2011). Google Scholar

20.

Y. L. WangY. N. PengZ. Bao, “Space-time adaptive processing for airborne radar with various array orientations,” IET Radar Sonar Navig., 144 (6), 330 –340 (1997). http://dx.doi.org/10.1049/ip-rsn:19971606 1751-8784 Google Scholar

21.

M. RueggE. MeierD. Nuesch, “Vibration and rotation in millimeter-wave SAR,” IEEE Trans. Geosci. Remote Sens., 45 (2), 293 –304 (2007). http://dx.doi.org/10.1109/TGRS.2006.887025 IGRSD2 0196-2892 Google Scholar

22.

D. NeedellR. Vershynin, “Signal recovery from incomplete and inaccurate measurements via regularized orthogonal matching pursuit,” IEEE J. Sel. Top. Signal Process., 4 (2), 310 –316 (2010). http://dx.doi.org/10.1109/JSTSP.2010.2042412 1932-4553 Google Scholar

23.

Y. Luoet al., “Micro-Doppler feature extraction for wideband imaging radar based on complex image orthogonal matching pursuit decomposition,” IET Radar Sonar Navig., 7 (8), 914 –924 (2013). http://dx.doi.org/10.1049/iet-rsn.2012.0327 1751-8784 Google Scholar

Biography

Yijun Chen received her MS degree in electrical engineering from the Institute of Information and Navigation, Air Force Engineering University (AFEU), Xi’an, China, in 2013, where she is currently working toward the PhD degree in electrical engineering. She is currently with the Radar and Signal Processing Laboratory, Institute of Information and Navigation, AFEU. Her research interests include signal processing and cognitive radar.

Qun Zhang received his MS degree in mathematics from Shaanxi Normal University, Xi’an, China, in 1988, and the PhD degree in electrical engineering from Xidian University, Xi’an, in 2001. He is a senior member of the Chinese Institute of Electronics (CIE), a committee member of the Radiolocation Techniques Branch, CIE. He was also a vice co-chair of the Technical Program Committee of Workshop for Space, Aeronautical and Navigational Electronics (WSANE), 2009.

Changzheng Ma received his BS degree from Lanzhou University, Lanzhou, China, in 1989, the MS degree from the Chinese Academy of Sciences, Beijing, China, in 1992, and the PhD degree from Xidian University, Xi’an, China, in 1999. He is currently a research fellow with the Department of Electrical and Computer Engineering, National University of Singapore, Singapore. His research interests include signal processing, radar imaging, and wireless communication.

Ying Luo received his MS degree in electrical engineering from the Institute of Information and Navigation, Air Force Engineering University (AFEU), Xi’an, China, in 2008, and the PhD degree in electrical engineering from AFEU, Xi’an, in 2013. He is currently with the Radar and Signal Processing Laboratory, Institute of Information and Navigation, AFEU. His research interests include signal processing and ATR in SAR and ISAR.

Tat Soon Yeo received his B.Eng. (Honors I) degree from the University of Singapore, Singapore, in 1979, the M.Eng. degree from the National University of Singapore (NUS), Singapore, in 1981, and the PhD degree from the University of Canterbury, New Zealand, in 1985, under a Colombo Plan Scholarship. He received the MINDEF-NUS Joint R&D Award in 1997 and the IEEE Millennium Medal in 2000.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.

Citation Download Citation

Yijun Chen, Qun Zhang, Changzheng Ma, Ying Luo, and Tat Soon Yeo "Micromotion feature extraction of radar target using tracking pulses with adaptive pulse repetition frequency adjustment," Journal of Applied Remote Sensing 8(1), 083569 (12 August 2014). https://doi.org/10.1117/1.JRS.8.083569

Published: 12 August 2014

Access the abstract

JOURNAL ARTICLE
14 PAGES

DOWNLOAD PAPER SAVE TO MY LIBRARY

GET CITATION

CITATIONS

Cited by 3 scholarly publications and 1 patent.

Explore citations on Lens.org

KEYWORDS

Feature extraction

Radar

Chemical species

Detection and tracking algorithms

Time-frequency analysis

Associative arrays

Image quality

1.

Introduction

2.

Micromotion Feature Extraction of Radar Target with Varying PRF Sampling

2.1.

Analysis of the Radar Echo of Micromotion Target

Fig. 1

Eq. (1)

Eq. (2)

Eq. (3)

Eq. (4)

Eq. (5)

2.2.

Micromotion Feature Extraction Method Based on OMP

Eq. (6)

Eq. (7)

Eq. (8)

Eq. (9)

Eq. (10)

Eq. (11)

Eq. (12)

Eq. (13)

Eq. (14)

3.

Micromotion Feature Extraction of Radar Target from Tracking Pulse with Adaptive PRF Adjustment

Eq. (15)

Eq. (16)

Fig. 2

Eq. (17)

4.

Experimental Results

Fig. 3

Fig. 4

Table 1

Fig. 5

Table 2

Fig. 6

Fig. 7

Table 3

Fig. 8

Table 4

Fig. 9

Fig. 10

Table 5

Fig. 11

5.

Conclusion

Acknowledgments

References

Biography

Show All Keywords

Keywords/Phrases

Search In:

Publication Years