Next Article in Journal
AI-Based Susceptibility Analysis of Shallow Landslides Induced by Heavy Rainfall in Tianshui, China
Previous Article in Journal
A Survey of Remote Sensing and Geographic Information System Applications for Flash Floods
Previous Article in Special Issue
Unsupervised Multistep Deformable Registration of Remote Sensing Imagery Based on Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Noise Robust Micro-Range Estimation Method for Precession Cone-Shaped Targets

National Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(9), 1820; https://doi.org/10.3390/rs13091820
Submission received: 30 January 2021 / Revised: 23 March 2021 / Accepted: 4 May 2021 / Published: 7 May 2021
(This article belongs to the Special Issue Signal and Image Processing for Remote Sensing)

Abstract

:
The estimation of micro-Range (m-R) is important for micro-motion feature extraction and imaging, which provides significant supports for the classification of a precession cone-shaped target. Under low signal-to-noise ratio (SNR) circumstances, the modified Kalman filter (MKF) will obtain broken segments rather than complete m-R tracks due to missing trajectories, and the performance of the MKF is restricted by unknown noise covariance. To solve these problems, a noise-robust m-R estimation method, which combines the adaptive Kalman filter (AKF) and the random sample consensus (RANSAC) algorithm, is proposed in this paper. The AKF, where the noise covariance is not required for the estimation of the state vector, is applied to associate m-R trajectories for higher estimation accuracy and lower wrong association probability. Due to missing trajectories, several associated segments which are parts of the m-R tracks can be obtained by the AKF. Then, the RANSAC algorithm is utilized to associate the segments and the complete m-R tracks can be obtained. Compared with the MKF, the proposed method can obtain complete m-R tracks instead of several segments, and avoids the influence of unknown noise covariance under low SNR circumstances. Experimental results based on electromagnetic simulation data demonstrate that the proposed method is more precise and robust compared with traditional methods.

Graphical Abstract

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.

2. Signal Model and Trajectory Extraction

2.1. Geometric Model and Signal Model

In this section, the geometric model and the signal model of the precession cone-shaped target are introduced. Figure 1 shows the model in the coordinate system O X Y Z , where the precession center of the target is located at the origin O . The target spins around the target symmetry axis at an angular velocity w s , and rotates around the axis Ω N at an angular velocity w c , which forms the precession motion. For simplicity, Ω N and the radar line-of-sight (rLOS) are on the Y O Z plane, and the symmetry is located at the Z-axis initially. The precession angle between the coning axis and the symmetric axis is α , and the initial target line-of-sight angle between the symmetric axis and the rLOS is β . The distance between the precession center and the radar is R 0 .
In this scattering center model, the target contains three main scattering centers, where the scattering center P 2 and P 3 are located at the intersection of the circular edge, and the bottom plane and P 1 is located on the smooth apex of the cone. The initial coordinate of the scattering center P 1 is p 10 = ( 0 , 0 , H ) and the initial coordinates of the scattering center P 2 and P 3 are p 20 = ( 0 , r a , h ) and p 30 = ( 0 , r a , h ) , where H denotes the distance between the precession center and the apex of the cone, h is the distance between the precession center and the bottom plane, and r a is the radius of the bottom plane. In this paper, we ignore the influence of bulk motion. According to [5], the range between scattering centers and radar can be obtained by
{ r 1 ( t ) = R 0 ( H h ) cos γ ( t ) r 2 ( t ) = R 0 + h cos γ ( t ) + r a sin γ ( t ) r 3 ( t ) = R 0 + h cos γ ( t ) r a sin γ ( t ) cos γ ( t ) = cos α cos ( α + β ) sin α sin ( α + β ) cos ( w c t )
where γ ( t ) denotes the angle between the rLOS and the symmetric axis at t time step. The micro-Range curve of P 1 is sinusoidal and the micro-Range curves of P 2 and P 3 can be described as sinusoidal curves, approximately [4].
Here, the transmitted signal s t ( t ) is the linear frequency modulated signal
s t ( t ) = r e c t ( t T p ) exp ( j 2 π f c t + j π μ t 2 )
where r e c t ( t / T p ) denotes a rectangular pulse of the pulse width T P , f c denotes the radar center frequency, and μ denotes the linear frequency modulated coefficient with the bandwidth B = μ T p . According to [28], the high-frequency scattering responses of a target can be described as a sum of responses from individual scattering centers, approximately. The echoes of the space target in a pulse are denoted as
s r ( t ) = i = 1 I a i r e c t ( t R i T p ) exp ( j 2 π f c ( t R i ) + j π u ( t R i ) 2 )
where I is the number of scattering centers on the target, R i denotes the distance between the i th scattering center and the radar which is assumed to be constant in a pulse, a i is the corresponding scattering coefficient of the i th scattering center, and c denotes the wave propagation velocity.
To obtain the range profile, the “dechirp” processing technique [29] is applied to obtain
s i f ( t ) = s r ( t ) s r e f ( t ) = i = 1 I a i r e c t ( t 2 R i / c T p ) exp ( j 4 π u c ( t 2 R r e f c ) R Δ j 4 π c f c R Δ + j 4 π u c 2 R Δ 2 )
where the reference signal s r e f ( t ) can be written as
s r e f ( t ) = r e c t ( t 2 R r e f / c T r e f ) exp ( j 2 π f c ( t 2 R r e f c ) + j π u ( t ^ 2 R r e f c ) 2 )
In the reference signal, the relative range R Δ = R i R r e f , R r e f denotes reference range, T r e f denotes the reference pulse width, which is a little greater than T p . Usually, u c 2 R Δ 2 < < f c R Δ c . Then, Equation (5) can be written as
s i f ( t ˜ ) = i = 1 I a i r e c t ( t ˜ 2 R Δ / c T p ) exp ( j 4 π R Δ c ( u t ˜ + f c ) )
where t ˜ = t 2 R r e f c , t ( 0 , T p ] . We proceed P times down-sampling in time t ˜ where the sampling interval is t = T p P . Then, the discrete frequency echoes can be expressed by
s i f ( p ) = i = 1 I a i exp ( j 4 π R Δ c f p )
where f p = f c + p B P . Let r n = n Δ r , n = 1 , 2 , , N where N Δ r denotes the scene swath. Then, the over-complete basis is constructed as follows:
X [ n , p ] = exp ( j 4 π r n f p / c )
Taking into account the noise, the echoes can be rewritten as
S 1 × P = A 1 × N X N × P + n o i s e
where A = [ 0 , 0 , , a i , 0 ] denotes the sparse coefficient vector of the echoes. Obviously A [ n ] = a i when R i R r e f = R Δ = n Δ r . The interval Δ r determines the range extraction precision which should be sufficiently small for constructing the over-complete dictionary and the super-resolution range estimation.

2.2. Trajectory Extraction via the Sparse Signal Representation Algorithm

In this section, the orthogonal matching pursuit (OMP) algorithm, which is a kind of sparse signal representation algorithm, is applied to obtain the m-R trajectories. According to the signal model of the space coning-shaped target introduced above, A denotes the scattering coefficients from the target scattering centers, and the indexes of non-zero elements denote the ranges between scattering centers and radar. A , which is a sparse vector, can be estimated by solving the following sparse optimization problem
A ^ = arg min A A 0   s . t .   S A X 2 < δ
where · 2 denotes the l 2 norm, · 0 denotes the l 0 norm, and δ denotes the noise power, which can be regarded as the heat noise from the receiver and accurately obtained by estimating the power of the non-target region in the HRRPs [30,31].
Solving Equation (10) is an NP-hard problem. The OMP algorithm can be used for (10) to attain the approximate solution. OMP is an iterative greedy algorithm. In this algorithm, the column of X , which is most strongly correlated with the residual, is chosen, and the corresponding element in A ^ can be calculated in the meantime. Then, the residual is updated by subtracting the contribution of this signal component at each iteration until residual power is lower than noise power δ . The estimated coefficient vector A ^ can be obtained at last where the non-zero elements are the scattering coefficients and the indexes of non-zero elements correspond to the range cells of scattering centers. The range of scattering centers can be regarded as the extracted scattering center trajectories. The extracted scattering center trajectories are shown in Figure 2.
In the low SNR condition, the power of echoes of the weak scattering center is close to noise power, and the inevitable estimation error makes the estimated noise power δ larger than that of the weak scattering center. Therefore, it may be impossible to extract the weak scattering centers and a large number of spurious trajectories are extracted undesirably. Besides, near the intersection of m-R curves, the multiple signal components are not separated for the limitation of resolution in the OMP algorithm. Hence, the range estimation of scattering centers is inaccurate and the scattering centers may be missing near the intersection.

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:
{ X k + 1 i = F X k i + n k i , n k i N ( 0 , R n ) Y k i = d i a g ( χ k i ) * ( C X k i ) + ( E d i a g ( χ k i ) ) v k i , v k i N ( 0 , R v )
where X k i is the state vector of i th scattering center at time step k , n k i is the process noise, which is a zero-mean Gaussian random process with the covariance R n , and F is the state transition matrix. Additionally, the measurement Y k i originates from i th scattering center at time step k , where χ k i { 0 , 1 } indicates whether the i th scattering center is measured at time step k . χ k i = 0 indicates the scattering center is missing and the measurement is regarded as the random noise trajectory. v k i denotes the measurement noise, which is a zero-mean Gaussian random process with the covariance R v . E is the identity matrix and C is the measurement matrix. The process noise n k i and the measurement noise v k i 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 X k i = [ x k i , x ˙ k i ] T , where x k i denotes the range and x ˙ k i = x k i x k 1 i Δ t denotes the velocity, and the velocity is approximately constant in a short time. Δ t denotes the time between two contiguous time steps. The state transition matrix F = [ 1 0 1 1 ] and the measurement matrix C = [ 1 0 0 0 ] 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 Y 1 : T , inlier threshold τ R
1: for each time step k do
2:  for each track do
3:   get the prediction of the state vector X k | k 1 i = F X k 1 | k 1 i .
4:   compute the inlier set J .
5:   if J =
6:    let y k i = C X ^ k i .
7:   else
8:    let y k i = m e a n ( J ) .
9:    end if
10:   calculate P ^ k | k 1 i , X ^ k | k i 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 X k i at k time step (Line 3):
X k | k 1 i = F X k 1 | k 1 i
Then, for filtering out the spurious trajectories and finding out the associated trajectories, the trajectories in the associated inlier set J are regarded as scattering center trajectories. The associated inlier set can be obtained by
J = { y k j | | y k j x ^ k i | < τ r 2 , | y k j x ^ k i t x ˙ ^ k i | < τ v }
where y k j denotes the j th extracted trajectory at the time step k . τ r and τ v 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 τ r equal to the range resolution of radar and let τ v = τ r / Δ t . 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 J is not empty, the mean value of the inlier set is regarded as the range of the corresponding scattering center. If J is empty, the scattering center is regarded to be missing at this time step where χ k i = 0 and the predicted range C x ^ k i 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. P ^ k | k 1 i , which denotes the estimation of the prior error covariance, can be obtained by
P ^ k | k 1 i = P ^ k 1 | k 2 i + Δ P ^ k 1 i *
Δ P ^ k 1 i * = ( Δ X ^ k 1 i Δ X ^ k 1 i T K k 1 i C P ^ k 1 | k 2 i ) / ( k 1 )
Δ X ^ k 1 i = X ^ k | k i X ^ k | k 1 i
where K k 1 i denotes the Kalman gain matrix for i the track at k 1 time step and X ^ k | k i denotes the posterior estimate for state X k i . P ^ k | k 1 i is constructed by P ^ k 1 | k 2 i , which denotes the estimation of prior error covariance at the previous time step. It is obvious that the calculation of P ^ k | k 1 i does not rely on process noise covariance R n . Therefore, the AKF can eliminate the influence of the unknown covariance matrix of process noise. The parameters of the AKF at the time step k are updated as follows (Line 10):
X ^ k | k i = X ^ k | k 1 i + K k i ( y k i C X ^ k | k 1 i )
K k i = P ^ k | k 1 i C T ( C P ^ k | k 1 i C T + R v ) 1
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 Z 0 as Z 0 Z , where Z is the set of cluttered measurements. The measurements in Z 0 can generate a hypothesis model f ( Γ ) , where Γ denotes the parameters of this model. During the hypothesis validation step, the inlier number of the hypothesis model is quantified:
χ = num ( { | z f ( Γ ) | < τ R , z Z } )
where z denotes the element in the set Z and τ R denotes the inlier threshold. Repeat the two steps several times and find the best model f ( Γ max ) which contains the most inliers from the generated hypothesis models.
We set the inlier threshold τ R = τ r / 2 and set the hypothesis model f ( Γ ) 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 β = 120 , the power scattered by the weak scattering center P 1 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 β = 110 , the trajectory of scattering center P 1 intersected with that of scattering center P 2 , 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 β = 120 . The power scattered by the weak scattering center P 1 was far less than that scattered by the other scattering centers.
Figure 3a shows the range–time image when S N R = 10   dB 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 a 1 = 0.09 , a 2 = 1 and a 3 = 0.30 . Obviously, the echo energy scattered by P 3 was about 100 times than that scattered by the weak scattering center P 1 . 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 P 1 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 P 1 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 P 1 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
MRMSE = m e a n i ( i = 0 N ( r t i r ^ t i ) 2 ) / N
where r t i denotes the true m-R of scatterer P i at t time step and r ^ t i denotes the estimated m-R of scatterer P i at t 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 R MKF R AKF , where R denotes the MRMSE, and relative time cost can be calculated by T MKF T AKF , where T 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 S N R 15   dB , 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 S N R 15   dB , 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 S N R > 15   dB , 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 β = 110 . The trajectory of scattering center P 1 intersected with that of scattering center P 2 .
Under the condition where the initial target line-of-sight angle is 110 , the range between the scattering center P 1 and the radar along the rLOS is close to that between the scattering center P 2 , and the radar which makes the m-R tracks of two scattering centers intersected. Figure 5a shows the range–time image when S N R = 10   dB 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 P 1 associate with the trajectory of the scattering center P 2 . In the meantime, the trajectory of the scattering center P 2 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 P 1 and P 2 . 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 β = 120 . The process noise covariance matrixes were R n = diag ( q , q ) , where q denotes the diagonal element of the R n and q is set to range from 0.1 to 1. We used two sets of data, where the precession rates were set to be w c = 4 π rad / s and w c = 6 π rad / s . Figure 7 shows the range–time image when the precession rates were w c = 4 π rad / s and w c = 6 π rad / s , respectively. We found that the m-R tracks on w c = 4 π rad / s were smoother than those on w c = 6 π rad / s . Therefore, the theoretical diagonal element of the process noise covariance matrix should be larger when w c = 6 π rad / s . 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 R n = diag ( 0.6 , 0.6 ) , indicating that q = 0.6 is closest to the theoretical noise variance when w c = 4 π rad / s . 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 R n = diag ( 0.7 , 0.7 ) , indicating that q = 0.7 is closest to the theoretical noise variance when w c = 6 π rad / s , which is a little larger than when w c = 4 π rad / s . 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 β = 120 . 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.

5. Conclusions

In this paper, we proposed a micro-Range estimation method combining the AKF and the RANSAC algorithm for a cone-shaped space target. To solve the unknown noise covariance problem and obtain complete m-R tracks in a low SNR circumstance, the AKF, which does not rely on the covariance matrix of process noise, was applied to associate the trajectories, and the RANSAC algorithm was used to obtain the complete m-R tracks by estimating the parameters of each m-R track. Compared with the MKF, the proposed method could obtain more accurate estimation results and complete m-R tracks. Compared with the Kalman filter and ESPRIT, the proposed method had a better performance on noise robustness and estimation accuracy. Compared with the GRT-CLEAN and GA-CLEAN methods, the proposed method had a fewer computational burden and obtained more precise estimation results. The proposed method can be applied in post-processing of ISAR image formation to improve precision and robustness. Its application in the imaging of a cone-shaped target and other space targets, such as rotating space debris, is part of our future research.

Author Contributions

Conceptualization, Z.Z.; methodology, Y.Z. and Z.Z.; software, Z.Z.; validation, Z.Z. and K.R.; investigation, Y.Z. and K.R.; data curation, Z.Z.; writing—original draft preparation, Z.Z.; writing—review and editing, Y.Z., L.D. and Y.L.; supervision, L.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Foundation of China under Grant 61771362 and Grant U1833203, in part by the 111 Project under Grant B18039.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Persico, A.R.; Clemente, C.; Gaglione, D.; Ilioudis, C.V.; Cao, J.; Pallotta, L.; De Maio, A.; Proudler, I.K.; Soraghan, J.J. On Model, Algorithms, and Experiment for Micro-Doppler-Based Recognition of Ballistic Targets. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 1088–1108. [Google Scholar] [CrossRef] [Green Version]
  2. Chen, V.; Li, F.; Ho, S.-S.; Wechsler, H. Micro-doppler effect in radar: Phenomenon, model, and simulation study. IEEE Trans. Aerosp. Electron. Syst. 2006, 42, 2–21. [Google Scholar] [CrossRef]
  3. Luo, Y.; Zhang, Q.; Yuan, N.; Zhu, F.; Gu, F. Three-dimensional precession feature extraction of space targets. IEEE Trans. Aerosp. Electron. Syst. 2014, 50, 1313–1329. [Google Scholar] [CrossRef]
  4. Bai, X.; Zheng, B. Imaging of rotation-symmetric space targets based on electromagnetic modeling. IEEE Trans. Aerosp. Electron. Syst. 2014, 50, 1680–1689. [Google Scholar]
  5. Bai, X.; Bao, Z. High-Resolution 3D Imaging of Precession Cone-Shaped Targets. IEEE Trans. Antennas Propag. 2014, 62, 4209–4219. [Google Scholar] [CrossRef]
  6. Bai, X.; Zhou, F.; Bao, Z. High-Resolution Three-Dimensional Imaging of Space Targets in Micromotion. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 1–13. [Google Scholar] [CrossRef]
  7. Zhang, W.; Fu, Y.; Nie, L.; Zhao, G.; Yang, W.; Yang, J. Parameter estimation of micro-motion targets for high-range-resolution radar using high-order difference sequence. IET Signal Process. 2018, 12, 1–11. [Google Scholar] [CrossRef]
  8. Fogle, O.R.; Rigling, B.D. Micro-Range/Micro-Doppler Decomposition of Human Radar Signatures. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 3058–3072. [Google Scholar] [CrossRef]
  9. Gao, H.; Xie, L.; Wen, S.; Kuang, Y. Micro-Doppler Signature Extraction from Ballistic Target with Micro-Motions. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 1969–1982. [Google Scholar] [CrossRef]
  10. He, X.; Tong, N.; Hu, X. High-Resolution Imaging and 3-D Reconstruction of Precession Targets by Exploiting Sparse Aper-tures. IEEE Trans. Aerosp. Electron. Syst. 2017, 53, 1212–1220. [Google Scholar] [CrossRef]
  11. Sun, Y.-X.; Ma, C.-Z.; Luo, Y.; Lin, Y.-Z.; Chen, Y.-A.; Zhang, Q. An interferometric-processing based three-dimensional imaging method for space rotating targets. In Proceedings of the 2016 CIE International Conference on Radar (RADAR), Guangzhou, China, 10–13 October 2016; pp. 1–5. [Google Scholar]
  12. Qi, M.W.; Xing, G.L.; Zheng, B. High-Resolution Three-Dimensional Radar Imaging for Rapidly Spinning Targets. IEEE Trans. Geosci. Remote Sens. 2008, 46, 22–30. [Google Scholar]
  13. Zhang, Q.; Yeo, T.S.; Tan, H.S.; Luo, Y. Imaging of a Moving Target With Rotating Parts Based on the Hough Transform. IEEE Trans. Geosci. Remote Sens. 2008, 46, 291–299. [Google Scholar] [CrossRef]
  14. Yu, N.; Bai, X.R.; Feng, Z.; Lei, L. Method for inverse synthetic aperture radar imaging of space debris using improved genetic algorithm. IET Radar Sonar Navig. 2017, 11, 812–821. [Google Scholar]
  15. Tsao, J.; Steinberg, B.D. Reduction of sidelobe and speckle artifacts in microwave imaging: The CLEAN technique. IEEE Trans. Antennas Propag. 1988, 36, 543–556. [Google Scholar] [CrossRef]
  16. Dai, D.H.; Wang, X.S.; Chang, Y.L.; Yang, J.H.; Xiao, S.P. Fully-Polarized Scattering Center Extraction and Parameter Esti-mation:P-ESPRIT Algorithm. In Proceedings of the 2006 CIE International Conference on Radar, Shanghai, China, 16–19 October 2006. [Google Scholar]
  17. Burrows, M.L. Two-Dimensional ESPRIT With Tracking for Radar Imaging and Feature Extraction. IEEE Trans. Antennas Propag. 2004, 52, 524–532. [Google Scholar] [CrossRef]
  18. Roy, R.; Kailath, T. ESPRIT-estimation of signal parameters via rotational invariance techniques. IEEE Trans. Acoust. Speech Signal Process. 1989, 37, 984–995. [Google Scholar] [CrossRef] [Green Version]
  19. Bai, X.; Feng, Z.; Zheng, B. High-Resolution Radar Imaging of Space Targets Based on HRRP Series. IEEE Trans. Geosci. Remote Sens. 2014, 52, 2369–2381. [Google Scholar] [CrossRef]
  20. Ren, K.; Du, L.; Lu, X.; Zhuo, Z.; Li, L. Instantaneous Frequency Estimation Based on Modified Kalman Filter for Cone-Shaped Target. Remote. Sens. 2020, 12, 2766. [Google Scholar] [CrossRef]
  21. Wang, J.; Wang, J.; Zhang, D.; Shao, X.; Chen, G. Kalman filtering through the feedback adaption of prior error covariance. Signal Process. 2018, 152, 47–53. [Google Scholar] [CrossRef]
  22. Fischler, M.A.; Bolles, R.C. Random Sample Consensus: A Paradigm for Model Fitting with Applications to Image Analysis and Automated Cartography. Read. Comput. Vis. 1981, 24, 381–395. [Google Scholar] [CrossRef]
  23. Michaelsen, E.; von Hansen, W.; Kirchhof, M.; Meidow, J.; Stilla, U. Estimating the essential matrix: Goodsac versus ransac. Photogramm. Comput. Vis. 2006, 1–6. [Google Scholar]
  24. Chum, O.; Matas, J.; Kittler, J. Locally Optimized RANSAC. In Joint Pattern Recognition Symposium; Springer: Berlin/Heidelberg, Germany, 2003; pp. 236–243. [Google Scholar]
  25. Matas, J.; Chum, O. Randomized RANSAC with Td,d test. Image Vis. Comput. 2002, 22, 837–842. [Google Scholar] [CrossRef] [Green Version]
  26. Torr, P.; Zisserman, A. MLESAC: A New Robust Estimator with Application to Estimating Image Geometry. Comput. Vis. Image Underst. 2000, 78, 138–156. [Google Scholar] [CrossRef] [Green Version]
  27. Raguram, R.; Chum, O.; Pollefeys, M.; Matas, J.; Frahm, J.-M. USAC: A Universal Framework for Random Sample Consensus. IEEE Trans. Pattern Anal. Mach. Intell. 2013, 35, 2022–2038. [Google Scholar] [CrossRef] [PubMed]
  28. Borovikov, V.A.; Kinber, B.Y. Geometrical Theory of Diffraction; IET Digital Library: Stevenage, UK, 1994. [Google Scholar]
  29. Quegan, S. Spotlight synthetic aperture radar: Signal processing algorithms. J. Atmos. Solar-Terr. Phys. 1997, 59, 597–598. [Google Scholar] [CrossRef]
  30. Du, L.; He, H.; Zhao, L.; Wang, P. Noise Robust Radar HRRP Target Recognition Based on Scatterer Matching Algorithm. IEEE Sens. J. 2015, 16, 1743–1753. [Google Scholar] [CrossRef]
  31. Tropp, J.A.; Gilbert, A.C. Signal Recovery From Random Measurements Via Orthogonal Matching Pursuit. IEEE Trans. Inf. Theory 2007, 53, 4655–4666. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Geometric model for precession target.
Figure 1. Geometric model for precession target.
Remotesensing 13 01820 g001
Figure 2. A scattering center trajectory extraction result with S N R = 10   dB .
Figure 2. A scattering center trajectory extraction result with S N R = 10   dB .
Remotesensing 13 01820 g002
Figure 3. Experimental results when β = 120 : (a) range–time image; (b) scattering center trajectory extraction result; (c) associated segments obtained by the AKF where curves in different colors represent different segments; (d) associated segments obtained by the MKF where curves in different colors represent different segments.
Figure 3. Experimental results when β = 120 : (a) range–time image; (b) scattering center trajectory extraction result; (c) associated segments obtained by the AKF where curves in different colors represent different segments; (d) associated segments obtained by the MKF where curves in different colors represent different segments.
Remotesensing 13 01820 g003aRemotesensing 13 01820 g003b
Figure 4. Estimated range versus true range when β = 120 : (a) AKF; (b) Kalman filter; (c) RANSAC; (d) ESPRIT; (e) GRT-CLEAN; (f) GA-CLEAN.
Figure 4. Estimated range versus true range when β = 120 : (a) AKF; (b) Kalman filter; (c) RANSAC; (d) ESPRIT; (e) GRT-CLEAN; (f) GA-CLEAN.
Remotesensing 13 01820 g004aRemotesensing 13 01820 g004b
Figure 5. Experimental results when β = 110 : (a) range–time image; (b) scattering center trajectory extraction result; (c) associated segments obtained by the AKF where curves in different colors represent different segments; (d) associated segments obtained by the MKF where curves in different colors represent different segments.
Figure 5. Experimental results when β = 110 : (a) range–time image; (b) scattering center trajectory extraction result; (c) associated segments obtained by the AKF where curves in different colors represent different segments; (d) associated segments obtained by the MKF where curves in different colors represent different segments.
Remotesensing 13 01820 g005
Figure 6. Estimated range versus true range when β = 110 : (a) AKF; (b) Kalman filter; (c) RANSAC; (d) ESPRIT; (e) GRT-CLEAN; (f) GA-CLEAN.
Figure 6. Estimated range versus true range when β = 110 : (a) AKF; (b) Kalman filter; (c) RANSAC; (d) ESPRIT; (e) GRT-CLEAN; (f) GA-CLEAN.
Remotesensing 13 01820 g006aRemotesensing 13 01820 g006b
Figure 7. The range–image of the targets. (a) The precession rate is w c = 4 π rad / s ; (b) The precession rate is w c = 6 π rad / s .
Figure 7. The range–image of the targets. (a) The precession rate is w c = 4 π rad / s ; (b) The precession rate is w c = 6 π rad / s .
Remotesensing 13 01820 g007
Figure 8. The MRMSE results with different process noise covariance matrixes. (a) The precession rate is w c = 4 π rad / s ; (b) The precession rate is w c = 6 π rad / s .
Figure 8. The MRMSE results with different process noise covariance matrixes. (a) The precession rate is w c = 4 π rad / s ; (b) The precession rate is w c = 6 π rad / s .
Remotesensing 13 01820 g008
Table 1. The basic simulation parameters.
Table 1. The basic simulation parameters.
ParameterValueParameterValue
Radar carrier frequency f c 10 GHzPrecession center height h 1 m
Pulse repletion frequency p r f 500 HzPrecession angle α
Bandwidth B 4 GHzPrecession rate w c 4 π rad / s
Target height H 3 mFrequency step Δ f 20 MHz
Target bottom radius r a 1 m
Table 2. MRMSE and time cost of m-R estimation when β = 120 . 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.
Table 2. MRMSE and time cost of m-R estimation when β = 120 . 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.
9 dB10 dB11 dB12 dB15 dB20 dB
MRMSE(m)KALMAN0.23770.20230.19620.19390.03740.0051
RANSAC0.00530.00530.00530.00520.00520.0050
GRT-CLEAN0.02360.02360.02350.02350.02350.0235
GA-CLEAN0.02390.02380.02380.02380.02380.0237
AKF0.00550.00545.4654 × 10−35.4306 × 10−35.3624 × 10−35.3401 × 10−3
MKF0.01930.00565.4865 × 10−35.4426 × 10−35.3652 × 10−35.3429 × 10−3
R MKF R AKF (mm)13.8560.20890.02110.01200.00280.0028
Time cost(s)KALMAN0.03040.03280.02960.03000.02950.0329
RANSAC746.75756.49742.16736.25755.67724.02
GRT-CLEAN56.02255.98755.95155.84756.00655.919
GA-CLEAN19.06319.09219.11019.05319.00819.028
AKF42.17110.9335.48872.93782.45022.4552
MKF61.75821.76386.88173.00552.46492.4000
T MKF T AKF (s)19.586810.83091.39300.06770.0147−0.0552
Table 3. MRMSE and time cost of m-R estimation when β = 110 . 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 sec-ond-most accurate or the time costs are the lowest.
Table 3. MRMSE and time cost of m-R estimation when β = 110 . 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 sec-ond-most accurate or the time costs are the lowest.
9 dB10 dB11 dB12 dB15 dB20 dB
MRMSE(m)KALMAN0.95670.83150.43600.16000.07500.0715
RANSAC0.00260.00240.00230.00230.00230.0023
GRT-CLEAN0.02710.02710.02710.02700.02700.0270
GA-CLEAN0.02590.02590.02560.02570.02540.0250
AKF0.00340.00340.00340.00320.00313.1952 × 10−3
MKF0.04420.01510.00370.00340.00323.2375 × 10−3
R MKF R AKF (mm)40.77511.6600.34320.22410.05720.0423
Time cost(s)KALMAN0.02970.03010.02930.03020.03230.0294
RANSAC771.16769.49764.25789.23775.22762.21
GRT-CLEAN55.64255.66655.54955.62555.63255.652
GA-CLEAN19.14619.14519.07919.08519.06919.062
AKF57.81420.25111.1199.77018.39328.3796
MKF78.85529.78712.58810.3148.35708.2540
T MKF T AKF (s)21.0419.53641.46890.5443−0.0362−0.1256
Table 4. MRMSE and time cost of the RANSAC, the MLESAC, and the R-RANSAC.
Table 4. MRMSE and time cost of the RANSAC, the MLESAC, and the R-RANSAC.
9 dB10 dB11 dB12 dB15 dB20 dB
MRMSE(m)MLESAC0.00540.00540.00530.00530.00530.0053
RANSAC0.00540.00540.00530.00530.00530.0053
R-RANSAC0.00780.00610.00530.00530.00530.0053
R RANSAC R MLESAC (m)0.742 × 10−60.004 × 10−60.001 × 10−6000
Time cost(s)MLESAC3.56143.25681.74241.54431.52031.1494
RANSAC3.27003.13421.64251.45621.42691.0216
R-RANSAC3.12773.01911.59481.42591.41761.0077
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhuo, Z.; Zhou, Y.; Du, L.; Ren, K.; Li, Y. A Noise Robust Micro-Range Estimation Method for Precession Cone-Shaped Targets. Remote Sens. 2021, 13, 1820. https://doi.org/10.3390/rs13091820

AMA Style

Zhuo Z, Zhou Y, Du L, Ren K, Li Y. A Noise Robust Micro-Range Estimation Method for Precession Cone-Shaped Targets. Remote Sensing. 2021; 13(9):1820. https://doi.org/10.3390/rs13091820

Chicago/Turabian Style

Zhuo, Zhenyu, Yu Zhou, Lan Du, Ke Ren, and Yi Li. 2021. "A Noise Robust Micro-Range Estimation Method for Precession Cone-Shaped Targets" Remote Sensing 13, no. 9: 1820. https://doi.org/10.3390/rs13091820

APA Style

Zhuo, Z., Zhou, Y., Du, L., Ren, K., & Li, Y. (2021). A Noise Robust Micro-Range Estimation Method for Precession Cone-Shaped Targets. Remote Sensing, 13(9), 1820. https://doi.org/10.3390/rs13091820

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop