1. Introduction
With increasing human activities in outer space, research of space targets has attracted more and more attention [
1,
2]. Usually, a target, or any structure on a target, undergoes micro-motion dynamics in addition to bulk motion. Precession is one of the typical micro-motions for space targets, such as cone-shaped targets [
3,
4,
5,
6]. The precession dynamic induces modulations on echoes, from which micro-Range (m-R) can be obtained [
7,
8]. By estimating and analyzing scattering centers’ m-R, high-resolution radar imaging, together with the motion and structure parameters, can be precisely acquired [
5].
A scattering model of the precession cone-shaped target is proposed in [
9,
10,
11]. There are three scattering centers on the model, one is located on the smooth apex of the cone and the other two are located at the intersection of the circular edge and the bottom plane, respectively. There are some weak scattering centers on the cone-shaped target whose scattering power is lower or close to noise power. Therefore, the sidelobes of dominant scattering centers and high noise interference make it difficult for m-R estimation of weak scattering centers in a low signal-to-noise ratio (SNR) condition.
Recently, many methods have been proposed to estimate m-R from wideband radar echoes. These methods can be divided into two categories: image domain methods and signal domain methods. In image domain methods, m-R curves can be extracted from the range–time image [
12,
13,
14]. In [
13], the generalized Radon transform with the CLEAN technique (GRT-CLEAN) is proposed to obtain each m-R curve, respectively, by searching for the peak of the parameter domain spectrogram. In the parameter domain spectrogram, the sinusoidal curve in range–time can be mapped to a point in the parameter domain by the generalized Radon transform (GRT) [
12,
13], and the CLEAN technique [
15] is applied to reduce the effect of unavoidable sidelobes in the range–time image. The radon transform can accomplish the connection of broken segments and the selection of overlapped segments. However, this method bears a heavy computational burden for constructing a multidimensional parameter domain spectrogram. A improved genetic algorithm with the CLEAN technique (GA-CLEAN) is proposed in [
14], which converts the complex imaging problem to a parametric optimization problem. Compared with the GRT-CLEAN method, this algorithm is not required to construct the spectrogram, and has a good convergence property. However, the limited range resolution of the range–time image limits the estimation accuracy of image domain methods.
Signal domain methods can obtain a super-resolution range estimation to improve estimation accuracy. There are two types of signal domain method. The first type of method obtains a super-resolution m-R estimation result by analyzing signal subspace [
16,
17,
18]. In [
18], estimation of signal parameters via rotational invariance techniques (ESPRIT), which exploits the rotational invariance of the underlying signal subspace, is applied to estimate m-R. However, this method is ineffective since it obtains the wrong signal subspace under high noise interference. In the second type of signal domain, trajectories of scattering centers are extracted by the sparse signal representation algorithm, and then the trajectory association algorithm is applied to associate the trajectories. In [
19], the Kalman filter is applied to associate trajectories. However, the Kalman filter suffers from the intersecting trajectory and spurious trajectory problems. In [
20], a modified Kalman filter (MKF) is proposed to associate the intersecting trajectories for instantaneous frequency estimation. In the MKF, the Kalman filter is used to associate the trajectories and the trajectory correction strategy is applied to correct wrong association. However, since trajectories of the weak scatter center cannot be extracted under low SNR circumstances, several associated segments, not the complete m-R tracks, will be obtained by the MKF. In addition, the MKF, where the process covariance initialized without any prior knowledge, is usually inaccurate or mismatched, provides poor performance in robustness, and results in wrong associations, especially under low SNR circumstances. Substantial association corrections will increase the cost time.
In this paper, a trajectory association method, which combines the adaptive Kalman filter (AKF) [
21] and the random sample consensus (RANSAC) algorithm [
22], is proposed to obtain complete m-R tracks and deal with the unknown noise covariance problem in a low SNR circumstance. In the proposed method, the AKF, which avoids the influence of the unknown covariance matrix of process noise, is applied to associate m-R trajectories. The AKF estimates the state vector by prior error covariance, which is directly reconstructed through online excavating of the posterior prior error covariance sequence, instead of needing the covariance matrix of process noise. Due to the missing trajectories, several associated segments which are parts of the m-R tracks can be obtained. Then, the RANSAC algorithm is utilized to associate the segments by estimating the parameters of each m-R track model. The RANSAC algorithm finds the parametric model for each m-R track, which contains maximum inlier trajectories, from many hypothetic parametric models, and the complete m-R tracks can finally be obtained.
There are many variants of the RANSAC algorithm. In [
23], a good sample consensus (GOODSAC), which replaces random sampling with an assessment driven selection of good samples, can foster precision and proper utilization of computational resources. However, a GOODSAC relies on the priori mathematical or geometric model to produce a list of good samples. In this paper, since there is no priori information of the mathematical or geometric model, a GOODSAC is not applicable for the connection step. In [
24], a locally optimized RANSAC improves the performance by choosing better samples. In this paper, all trajectories in one segment are used to form a sample set, and we perform a complete search to test all sample sets formed by all segments to ensure precise estimation results. Therefore, samples do not need to be chosen, and the locally optimized RANSAC is not applicable for the connection step. Maximum Likelihood SAC (MLESAC) utilizes probability distribution of error by inlier and outlier to evaluate a hypothesis [
25]. It adopts the same sampling strategy as RANSAC to generate putative solutions, but chooses the solution that maximizes the likelihood of inliers and outlier rather than the number of inliers. In this paper, since almost all the outliers have been filtered out in the associated step, MLESAC should provide almost the same results, with more computation burden compared to the RANSAC algorithm. We will offer a more detailed analysis with the comparative experiment in
Section 4. In [
26], the Randomized RANSAC (R-RANSAC) performs a preliminary test to reduce time. Hypothesis evaluation is only performed when the generated hypothesis passes the preliminary test. However, a valid hypothesis may be mistakenly rejected by the preliminary test [
27]. In summary, to obtain precise and robust results, the RANSAC algorithm is chosen to associate the segments.
Compared with the MKF, the proposed method can obtain complete m-R tracks instead of segments, and has a better performance on the estimation accuracy and time cost. Firstly, by estimating the parameters of the m-R track model in the proposed method, the complete m-R tracks rather than segments can be obtained, which overcome the shortcoming of the MKF. Secondly, when noise covariance is inaccurate or mismatched, the proposed method can obtain higher accuracy estimation results and decrease wrong association probability under low SNR circumstances. Since the covariance matrix of process noise is not required for the calculation of prior error covariance, the performance of AKF is not affected by the unknown covariance matrix of process noise. Finally, since it decreases wrong associations and association corrections, the AKF can significantly reduce time cost. Compared with ESPRIT, the proposed method has better noise robustness. Compared with image domain methods, the proposed method can obtain a super-resolution range estimation to improve estimation accuracy. Experimental results on electromagnetic simulation data demonstrate that the proposed method is more accurate and more time-saving compared with the methods introduced above.
The organization of this paper is as follows. The signal model and trajectory extraction method are provided in
Section 2. The proposed association method is developed in
Section 3. The detailed experimental results, based on electromagnetic simulation data, are given in
Section 4.
Section 5 concludes this paper.
3. Trajectory Association by Adaptive Kalman Filter
In this section, we propose an algorithm combining the AKF and the RANSAC algorithm for associating the trajectories. The associated method can associate the trajectory of the same scattering center, and the associated method can also filter out spurious trajectories extracted from the noise component. After that, the track of the scattering center, which is regarded as the m-R curve of the scattering center, can be obtained. In [
20], a trajectory association method based on the MKF is proposed to estimate the instantaneous frequency. In MKF, the Kalman filter is applied to associate trajectories, and the trajectory correction strategy is applied to correct wrong associations. For the association of m-R trajectories, the missing trajectories of weak scattering centers will make the track broken. In addition, the process noise matrix, which is initialed without any prior knowledge, is usually inaccurate, or even mismatched. To solve these problems, the AKF, which does not rely on prior and correct knowledge about the dynamical model, is applied to associate the trajectories, and the broken segments are processed to obtain the complete m-R track. The details of the proposed method are given in the following.
Due to precession dynamics, the m-R tracks are modulated as sinusoidal. Therefore, to describe the trajectory association problem, the discrete system with the linear Gaussian dynamic and measurement model is proposed, which satisfies the following equation:
where
is the state vector of
th scattering center at time step
,
is the process noise, which is a zero-mean Gaussian random process with the covariance
, and
is the state transition matrix. Additionally, the measurement
originates from
th scattering center at time step
, where
indicates whether the
th scattering center is measured at time step
.
indicates the scattering center is missing and the measurement is regarded as the random noise trajectory.
denotes the measurement noise, which is a zero-mean Gaussian random process with the covariance
.
is the identity matrix and
is the measurement matrix. The process noise
and the measurement noise
are statistically independent.
In the range–time domain, the instantaneous m-R track can be described as a smooth scattering center track, and the constant velocity (CV) model is adapted to describe its shape. In the CV model, the state vector is , where denotes the range and denotes the velocity, and the velocity is approximately constant in a short time. denotes the time between two contiguous time steps. The state transition matrix and the measurement matrix are assumed to be invariant with time. The proposed method is summarized in Algorithm 1.
Algorithm 1. The associated method combining the AKF and the RANSAC algorithm. |
Input: extracted trajectories , inlier threshold |
1: for each time step do |
2: for each track do |
3: get the prediction of the state vector . |
4: compute the inlier set . |
5: if |
6: let . |
7: else |
8: let . |
9: end if |
10: calculate , and update the parameters in AKF |
11: use the trajectory correction strategy if association mistakes occur. |
12: end for |
13: end for |
14: use the associated segments to estimate the parameters of each m-R model by the RANSAC algorithm |
For each track, the AKF is applied to associate the trajectory at each time step. Firstly, the one-step Kalman predictor is used to obtain the prediction of the state vector
at
time step (Line 3):
Then, for filtering out the spurious trajectories and finding out the associated trajectories, the trajectories in the associated inlier set
are regarded as scattering center trajectories. The associated inlier set can be obtained by
where
denotes the
th extracted trajectory at the time step
.
and
denote the inlier criterion of range and velocity, respectively. For the limitation of radar range resolution, the extracted trajectories in the prediction range cell are regarded as the trajectories of the track. Therefore, let
equal to the range resolution of radar and let
. The criterion of velocity can avoid the association error at the intersection of two tracks, where the states of two tracks are close, but the velocities of two tracks have a considerable difference. Due to the influence of the spurious trajectories, there may be more than one trajectory in the inlier set. If
is not empty, the mean value of the inlier set is regarded as the range of the corresponding scattering center. If
is empty, the scattering center is regarded to be missing at this time step where
and the predicted range
is regarded as the range of the missing trajectory (Line 4–9).
Next, update the state vector using the extracted trajectory by the AKF at the current time step. In the AKF, the corresponding prior covariance is replaced by the prior error covariance through the online excavating posterior sequence.
, which denotes the estimation of the prior error covariance, can be obtained by
where
denotes the Kalman gain matrix for
the track at
time step and
denotes the posterior estimate for state
.
is constructed by
, which denotes the estimation of prior error covariance at the previous time step. It is obvious that the calculation of
does not rely on process noise covariance
. Therefore, the AKF can eliminate the influence of the unknown covariance matrix of process noise. The parameters of the AKF at the time step
are updated as follows (Line 10):
After the AKF associates the trajectories, we use the trajectory correction strategy if association mistakes occur (Line 11), which is detailed in [
20]. Repeat the associated step by AKF (Line 3–10) and the trajectory correction strategy (Line 11) in the dwell time and, then, several m-R segments can be obtained due to the missing trajectories. The RANSAC algorithm is utilized to associate the segments by estimating the parameters of each m-R track model (Line 14). The RANSAC algorithm is a technique that estimates the parameters of a single from a batch of cluttered measurements. Here, we provide a simple introduction to this algorithm. The RANSAC algorithm consists of two repeated steps: a hypothesis generation step and a hypothesis validation step. During the hypothesis generation step, a minimum subset is used to form a hypothesis. Define a minimum subset
as
, where
is the set of cluttered measurements. The measurements in
can generate a hypothesis model
, where
denotes the parameters of this model. During the hypothesis validation step, the inlier number of the hypothesis model is quantified:
where
denotes the element in the set
and
denotes the inlier threshold. Repeat the two steps several times and find the best model
which contains the most inliers from the generated hypothesis models.
We set the inlier threshold and set the hypothesis model as a sinusoidal function, since the m-R track is modelled as sinusoidal. Since the trajectories in the segment must be the inliers of the hypothesis model of this segment due to the association by the AKF, each hypothesis model can be directly established by each segment. Here, we use the cftool MATLAB toolbox to estimate sinusoidal parameters for each segment. The quantity of segments is far less than that of trajectories, so that the proposed method generates a few hypothesis models and takes a short amount of time. To ensure precise estimation results, we calculate the inliers of every hypothesis model and use all trajectories of one segment as a set of samples.
4. Experimental Results
In this section, we use the electromagnetic computation data of the cone-shaped target in the comparison experiment to analyze the robustness of the proposed algorithm, and show the estimation and computational advantages compared with conventional algorithms. The electromagnetic computation data were obtained by FEKO using the physical optics (PO) method. The cone-shaped target in the experiments is shown in
Figure 1. With reference to other papers [
9,
10,
19,
20], we set simulation parameters of the target model and the signal model. The basic simulation parameters are listed in
Table 1. The time cost in Experiment 1 and Experiment 2 of the original paper was the executed time from the beginning of processing the echo data to obtaining the m-R estimation result. Our experiments were executed on an Intel Core CPU i7-7700 @ 3.60 GHz, and the RAM was 8 GB. The software platforms were Windows 10 and MATLAB 2017. In aiming to show good performance in solving the missing trajectory and trajectory intersection problems, we chose two sets of data under different circumstances, where
was 120 and 110 degrees, respectively, in Experiment 1 and Experiment 2. When
, the power scattered by the weak scattering center
was far less than that scattered by the other scattering centers. Therefore, a large number of missing trajectories and spurious trajectories existed, which made the track broken. When
, the trajectory of scattering center
intersected with that of scattering center
, and the intersection made it difficult to associate. We compared the performance of the proposed method with the existing methods, such as the MKF [
20], the RANSAC algorithm [
22], the Kalman filter [
19], the ESPRIT [
18], the GRT-CLEAN method [
13], and the GA-CLEAN method [
14]. For the RANSAC algorithm, we used the RANSAC algorithm to process the m-R trajectories directly instead of the associated segments. In Experiment 3, we discuss the influence of different process noise covariance matrixes on the estimation accuracy of the AKF and the MKF. In Experiment 4, we compared the performance of the RANSAC and other variants.
Experiment 1: In this experiment, we set the target line-of-sight angle . The power scattered by the weak scattering center was far less than that scattered by the other scattering centers.
Figure 3a shows the range–time image when
and the dwell time is 0.5 seconds. In this experiment, the scattering coefficients could be roughly estimated by the OMP algorithm, and the normalized estimated scattering coefficients were
,
and
. Obviously, the echo energy scattered by
was about 100 times than that scattered by the weak scattering center
.
Figure 3b shows the trajectory extraction result by the OMP algorithm. We could find that there were many missing trajectories when the dwell time was 0.25 s, making the m-R track of
disconnected.
Figure 3c,d shows the associated segments obtained by the AKF and the MKF, respectively. In these figures, curves in different colors represent different segments. The missing trajectory and the inaccurate estimation make the inlier radio decrease, and the new track is created by the RANSAC algorithm. Therefore, the track of
obtained by the MKF is broken into two segments.
Figure 4a–f shows the estimated range results and true ranges for different methods, respectively. The proposed method and the RANSAC algorithm can obtain accurate m-R estimation. ESPRIT cannot obtain an accurate m-R estimation of the weak scattering center
under a low SNR circumstance. The wrong associations by the Kalman filter occur when the missing trajectories exist. The GRT-CLEAN and GA-CLEAN algorithms can obtain complete and accurate m-R curves, but they are restricted by the range resolution of the m-R image.
To evaluate the estimation accuracy and time cost of these methods, the mean variation of estimated root-mean-square errors (MRMSEs) with different SNRs and the time cost for different methods are given in
Table 2. The MRMSE can be obtained by
where
denotes the true m-R of scatterer
at
time step and
denotes the estimated m-R of scatterer
at
time step. Since the MRMSE of the proposed method and the MKF are approximate, we use the relative MRMSE and relative time cost to analyze the performance. The relative MRMSE can be calculated by
, where
denotes the MRMSE, and relative time cost can be calculated by
, where
denotes time cost. When the relative MRMSE and the relative time cost are positive, the proposed method has a better performance in estimation accuracy and time cost. Since the MKF cannot obtain a complete track, we use the RANSAC algorithm to process the segments and obtain the complete m-R tracks for the calculation of the relative MRMSE and time cost. In
Table 2, results underlined indicate that the estimation results are the most accurate, or time costs are the lowest, and results in bold indicate that the estimation results are the second-most accurate, or the time costs are the lowest. ESPRIT cannot estimate the m-R of weak scattering center, so that we do not discuss this method in the table.
It is evident that the Kalman filter takes the least time, but its noise robustness is poor. Under the circumstance when
, it cannot associate trajectories correctly due to the influence of missing trajectories and spurious trajectories. GRT-CLEAN and GA-CLEAN algorithms obtain less accurate estimation results for relying on the resolution of the range–time image. The RANSAC can obtain the most accurate results, but the time cost of the RANSAC is far more than that of the proposed method. The proposed method can effectively deal with the missing trajectories and spurious trajectories problem, especially in the case of a low signal-to-noise ratio. In
Table 2, the proposed method has a lower MRMSE than the MKF, especially under low SNR circumstances. When the
, the proposed method cost less time than the MKF. The proposed method can obtain a higher accuracy estimation and decline the wrong association probability, and the number of association corrections by the RANSAC algorithm is decreased, which significantly reduces time cost. When
, both the MKF and AKF can avoid the existence of the wrong association, and the calculation of prior error covariance in AKF brings a little more computational burden compared with the error covariance matrix estimation step in MKF. Therefore, the proposed method costs a little more time than the MKF. Overall, the proposed method has a better performance in both accuracy and time cost than other methods.
Experiment 2: In this experiment, we set the target line-of-sight angle . The trajectory of scattering center intersected with that of scattering center .
Under the condition where the initial target line-of-sight angle is
, the range between the scattering center
and the radar along the rLOS is close to that between the scattering center
, and the radar which makes the m-R tracks of two scattering centers intersected.
Figure 5a shows the range–time image when
and the dwell time was 0.5 s.
Figure 5b shows the trajectory extraction result by the OMP algorithm. The intersection makes it difficult to associate the trajectories.
Figure 5c,d shows the associated segments obtained by the AKF and the MKF, respectively. Obviously, the intersection increased the wrong association probability and made the track broken.
Figure 6a–f shows the estimated m-R results and the true m-R curves for different methods. It is obvious that there was a wrong association obtained by the Kalman filter near the intersection, which made the trajectory of the scattering center
associate with the trajectory of the scattering center
. In the meantime, the trajectory of the scattering center
was shaken near the intersection because of the influence of spurious trajectories. The proposed method and the RANSAC algorithm could obtain accurate m-R estimation. The ESPRIT algorithm could not obtain the correct estimation results for the influence of the weak signal power of scattering center
and
. The Kalman filter obtained the wrong association near the intersection. The GRT-CLEAN and GA-CLEAN algorithms could obtain complete and accurate m-R curves, but they were restricted by the range resolution of the m-R image.
In
Table 3, results underlined indicate that the estimation results are the most accurate or the time costs are the lowest, and results in bold indicate that the estimation results are the second-most accurate or the time costs are the lowest. As is shown in this table, the Kalman filter cannot correctly associate the trajectories near the intersections in all cases of SNR, even if it costs the least time. The GRT-CLEAN and the GA-CLEAN can extract m-R curves, but the low range resolution of the range–time image will result in poor accuracy of parameter searching, and m-R estimation accuracy will be worse. The OMP algorithm in the proposed method can obtain a super-resolution m-R estimation result. Therefore, the estimation accuracy of GRT-CLEAN and the GA-CLEAN is not as good as that of our algorithm. The time cost of the RANSAC is far more than that of the proposed method. The AKF has a better performance on estimation accuracy than the MKF, and the AKF cost less time in a low SNR circumstance. When the signal-to-noise ratio is reduced, the proposed method guarantees accurate estimation results while costing only a small amount of time. In summary, its performance for estimation and time cost is significantly better than other traditional methods.
Experiment 3: The influence of different process noise covariance matrixes on the estimation accuracy of the AKF and the MKF.
In this experiment, the SNR of echoes were set to be 10 dB and the initial target line-of-sight angle was
. The process noise covariance matrixes were
, where
denotes the diagonal element of the
and
is set to range from 0.1 to 1. We used two sets of data, where the precession rates were set to be
and
.
Figure 7 shows the range–time image when the precession rates were
and
, respectively. We found that the m-R tracks on
were smoother than those on
. Therefore, the theoretical diagonal element of the process noise covariance matrix should be larger when
.
Figure 8 shows the MRMSE of the AKF and MKF with different process noise covariance matrixes, where the blue markers show the MRMSE of the MKF and the red markers show the MRMSE of the AKF. In
Figure 8a, the MRMSE is smallest when the noise variance of MKF is
, indicating that
is closest to the theoretical noise variance when
. Since the calculation of the AKF does not need the noise variance, the MRMSEs with different process noise covariance are the same, and the MRMSE of the AKF is lower than that of the MKF. In
Figure 8b, the MRMSE is the smallest when the noise variance of MKF is
, indicating that
is closest to the theoretical noise variance when
, which is a little larger than when
. The experimental results demonstrate that the AKF can obtain higher estimation accuracy than the MKF with different process noise covariance matrixes.
Experiment 4: In this experiment, we compared the performance of the RANSAC and other variants.
In this experiment, the initial target line-of-sight angle was
. Other parameters are listed in
Table 1 of the original paper. We used the AKF to obtain the associated segments and used the RANSAC, the MLESAC, and the R-RANSAC to connect the segments, respectively. As the GOODSAC and the locally optimized RANSAC were not are not applicable, we did not discuss and compare these two methods in this experiment. The MRMSE and time cost are provided in
Table 4:
As is shown in
Table 4, the MMRSE of the MLESAC and the RANSAC is the same because almost all the outliers were filtered out in the associated step by the AKF. However, the time cost of MLESAC is more than that of RANSAC for the calculation of probability distribution of error by inlier and outlier. The R-RANSAC costs less time because of the preliminary test. However, the test may reject the inlier segment, which makes the method less robust. In the proposed method, the cost time of the associated step of the AKF is far more than that of the connected step by RANSAC. To obtain precise and robust results, the RANSAC algorithm is chosen to associate the segments.