1. Introduction
In the traditional multi-target tracking (MTT) [
1] scenario, a single target generates, at most, one measurement at a time. In this case, the target can be regarded as a point target and its shape can be ignored. However, with the development of modern high-resolution sensors, the target will occupy multiple resolution cells of the sensor, so the point target model is no longer applicable. Instead, the target should be modeled as an extended target (ET), where a single target generates multiple measurements at each time.
In the past few years, the extended target tracking problem has attracted a lot of attention; the references [
2,
3] summarize the extended target tracking problem. At present, the measurement model of the extended target is mostly based on the non-uniform Poisson process measurement model proposed by Gilholm and Salmond [
4], which assumes that the measurement sources of the extended target are randomly distributed around the center of mass of the extended target, and the measurements received by sensors about the target and clutter at the same time all follow Poisson distribution. The extended target cannot be modeled by the point target model, and its shape model should be considered. At present, the modeling methods of extended target shape mainly include random matrix (RM) [
5,
6], random hypersurface model (RHM) [
7], and Gaussian process (GP) [
8]. The RM method assumes that the shape of the target is approximately an ellipse, while the RHM and GP methods no longer limit the shape of the target to an ellipse, but a more irregular shape; however, these two modeling methods also result in more computational complexity.
In recent years, the extended target tracking method based on random finite set (RFS) has received extensive attention. Since the target tracking framework based on RFS has a strict Bayesian theory foundation, it can simultaneously estimate the number of targets and the state without considering the data association. Therefore, it is particularly suitable for multiple extended target tracking problems. The RFS framework for MTT proposed by Mahler [
9], driven by the inherent combinatorial characteristics of Bayesian multi-objective filters, has led to the utilization of probability hypothesis density (PHD) [
10], cardinalized PHD (CPHD) [
11], and multi-Bernoulli filters [
12,
13] to approximate optimal Bayesian filtering. Based on the extended target Poisson model proposed by Gilholm, Mahler derived the measurement update equation of the PHD filter. However, the measurement update equation remains intractable.
Granstrom proposed a tractable extended target Gaussian mixture PHD filter algorithm (ET-GM-PHD) [
14], and proposed a distance partition (DP) method to partition the measurement set. The Gaussian-inverse-Wishart PHD (GIW-PHD) filtering algorithm [
6], which combines RM with the ET-PHD framework, can not only estimate the motion state of the target, but also estimate the extended state of the target.
How to partition the measurement set accurately and effectively is one of the main tasks of extended target tracking. In multiple extended target tracking, a single target generates multiple measurements, and with the interference of clutter, it is difficult to establish the correspondence between the measurements and the target. The ET-GM-PHD filtering algorithm needs to consider all partitions of the target measurement set, which is computationally difficult in practical applications. In application scenarios with a large number of targets and measurements, the DP algorithm has high computational complexity. In addition, it is essentially a hierarchical single-chain clustering algorithm, which lacks robustness and is sensitive to noise. In addition to DP, K-means++ clustering [
15], spectral clustering [
16], and expectation maximization (EM) [
6] partitioning algorithms are also included. It is difficult to deal with measurement set partition in complex scenes using single partition method. Therefore, finding a combination of multiple clustering algorithms can bring more accurate division results and provide a guarantee for subsequent high-precision filtering. A comparison of the characteristics of existing extended object tracking algorithms based on different partitioning methods is presented in
Table 1. In addition, in recent years, Ref. [
17] proposed an algorithm to create data particles in two-dimensional feature space through its geometric partition (GD). This method solves the classification problem of datasets with similar centers of gravity and overlapping objects described by different labels. Therefore, the data particle geometric partition method has broad prospects for the measurement partition problem of multiple close distance targets. However, the drawback of this method is that it can only be applied to 2D datasets.
The number of measurements in the measurement set is the main factor that affects the speed and number of partitions. How to effectively reduce the number of measurements participating in clustering without affecting the accuracy of partitioning is the focus of this paper. To solve this problem, this paper proposes an extended multi-target tracking algorithm based on temporal correlation, which is based on ET-GM-PHD and GM-PHD filtering algorithms. The measurement set is hierarchically processed at the same time, which improves the filtering accuracy and reduces the computational complexity. Experimental results show that the proposed method improves the overall filtering performance compared with the DP partition algorithm.
To summarize, the main contributions of this work are listed as follows.
- -
A hierarchical measurement partitioning method based on spatio-temporal correlation is proposed;
- -
The ET-GM-PHD filter based on spatio-temporal correlation is proposed, and the corresponding implementation scheme is given;
- -
The performance of the proposed method is assessed via simulations.
The rest of this paper is organized as follows.
Section 2 introduces relevant mathematical modeling. In
Section 3, the ET-GM-PHD filter based on the spatio-temporal correlation method is proposed. The simulation results of the proposed method are analyzed in
Section 4. Conclusions are presented in
Section 5.
Table 1.
Extended target tracking algorithm based on different measurement partition methods.
Table 1.
Extended target tracking algorithm based on different measurement partition methods.
Measurement Partition Method | Main Characteristic |
---|
Traditional method | Updating requires iterating over all possible subsets of measurement partitions |
| The computational complexity is high, and it is difficult to |
Algorithm based on distance partition [6,14] | accurately divide the measurement sub-set corresponding to the target |
| in the case of a large number of targets and lacks robustness. |
| The center point needs to be found iteratively. |
K-means++ clustering [15] | When the number of targets is large and the distance is close, it is difficult to |
| distinguish the measurement sets of different targets. |
| This method is based on adjacency matrix propagation, which can divide the |
Spectral clustering [16] | measurement set more accurately. However, this method is |
| very sensitive to clutter, and it needs to eliminate the clutter in addition, |
| and it is complex to deal with dense clutter scenes. |
| The EM measurement partitioning algorithm converges to |
| the nearest local maximum of the likelihood function. |
Expectation maximization (EM) [6] | Since the EM partitioning is initialized using the predicted GIW components, |
| the partitioning performance is sensitive to |
| the maneuverability of the motion model and prediction results. |
Spatio-temporal correlation (proposed) | The measurement set partitions are hierarchically processed simultaneously. |
| The computational efficiency and accuracy are higher. |
2. Mathematical Modeling
Multiple extended target tracking is used to estimate the number and motion state of multiple extended targets based on sensor observation data in a noisy environment. The number of targets and measurements changes over time with the birth, derivation, and death of targets and missed detection of sensors. Mahler proposed a multi-extended target tracking framework based on a Poisson multi-target model. In the RFS multi-extended filtering framework, the state and measurement of the target are, respectively, expressed as and , where is the state of a single target and is the number of targets at time k. is a single measurement value and is the number of measurements at time k.
Since the prediction formula of ET-GM-PHD is the same as that of the standard GM-PHD filter, only the update formula of the ET-GM-PHD filter is considered here. The measurement update equation can be expressed as follows:
where
represents the measurement pseudo-likelihood function and is defined as follows:
where
denotes that the measurement set
is divided into non-empty measurement units
W;
represents the measurement likelihood function generated by a single target;
represents the expected number of target measurements;
represents the probability of object detection;
denotes the average number of clutter measurements; and
is the uniform clutter density in the surveillance area. For more details of the ET-GM-PHD algorithm, we can refer to reference [
14].
It can be seen from Equation (
2) that, in the multi-extended target tracking framework, the measurement set needs to be divided, and the accuracy of the division result directly affects the subsequent filtering result. The DP algorithm proposed in reference [
14] is a widely used measurement set partitioning method. This method uses the distance between the measurements to cluster, but when two or more extended targets are close to each other, the measurements belonging to different targets will fall into the same cluster unit after partitioning, resulting in the reduction in the number of measurement units. However, the number of measurement values in the existing measurement units increases. To solve this problem, Granstrom proposed a second partition method [
14], which divided the divided measurement set again, but this method had a large amount of calculation and poor real-time tracking. In 2020, a research paper [
18] used the method based on shared nearest neighbors (SNN) of directed graph to divide the measurement values. These two methods do not consider the correlation between adjacent extended target measurements. Therefore, this paper proposes a partition method of an extended target measurement set, based on spatio-temporal correlation.
3. Partition Method Based on Spatio-Temporal Correlation
Considering the correlation between the adjacent measurements and the target motion model, this paper proposes a fast partition and hierarchical processing algorithm for the measurement set. When multiple extended targets cross or are close to each other, the algorithm can also achieve accurate cardinality estimation and target state estimation, and greatly reduces the amount of calculation. At the same time, the method is not affected by the measurement density.
The partitioning method proposed in this paper mainly includes two parts: partitioning the living target measurement set and partitioning the newborn target measurement set. Firstly, the two basic clustering algorithms used in this method are K-means++ clustering algorithm and DBSCAN clustering algorithm.
3.1.1. The K-Means++ Clustering Algorithm
In multiple extended target tracking, since the measurement set at each time contains measurements from different targets and clutter, it is necessary to use a clustering algorithm to cluster the measurement set, i.e., measurement set partition. The K-means clustering algorithm is a kind of commonly used clustering algorithm, wherein K denotes the number of clusters. Since the number of extended targets and clutter varies at different times, K is also time-varying. If the K traversal interval is (where ,), the number of measurement partitions will be quite large when the number of measurements is large, so the value range needs to be narrowed. In addition, the K-means clustering algorithm is unstable because it randomly selects the initial points.
To solve the above problems, the K-means++ algorithm was proposed in the literature [
15]. The specific algorithm steps are as follows:
Randomly select a center from the dataset ;
Select the next central point with probability , ;
Repeat step 2 until there are k central points, ;
For each , set the cluster to be the point set in X that is closer to than , where ;
For each , let be the centroid of all points in , ;
Repeat steps 4 to 5 until is no changes.
3.1.2. The DBSCAN Clustering Algorithm
DBSCAN is a density-based clustering algorithm. The algorithm divides the entire sample set quickly by using the sample sets which are closely connected in each group. Simultaneously, the algorithm groups closely connected sample sets into one class to form a cluster category, while not requiring the prior input of the number of samples for clustering. The basic idea of DBSCAN clustering algorithm is that, for each measurement in the partition cell, the neighborhood of its given radius contains no less than a given number of measurements.
Therefore, DBSCAN only needs two input parameters: the size of the radius and the minimum number of measurements in the neighborhood of each measurement in the partition cell. For a sample set of any shape, there exists some point from which any point in the sample set can be reached by a “dense” path. In other words, the points in the cluster must be density-connected.
Suppose that, at time
k, it exists a set of measurements
(subscript
k is ignored for simplicity). If a measurement is a core point, then its neighborhood
(see reference [
19] for a definition of neighborhood) contains at least
M measurements.
The steps of the DBSCAN algorithm are described as follows:
Set the cluster count and the set of core points ;
Traverse the measurement values in the measurement set Z, and determine whether it is a core point. If it is, add it to the set C, i.e., ;
Traverse the set of core points in C that have not been visited. Increment the cluster count . Set the category of as C and marked it as visited.
According to the current core point and the count of clusters , traverse all points in the neighborhood of , setting the category of as and marking them as visited. If is not visited and is a core point, set and repeat step 4 until the condition is not satisfied;
Repeat steps 3 and 4 until all the core points are visited;
The points in the measurement set Z that do not belong to the core points and have not been visited are set as noise points.
After the above steps, the measurement set Z is divided into the measurement set of the target and the noise point.
3.2. Division of Survival Target Measurement Set
For an extended target that survives at time k, according to its motion model, the target position at the previous time is correlated with the target position at the current time, and the measurements generated by the extended target at these two times are also correlated. The measurement set at time k is set to , denotes the number of measurements at time k, and the prediction of the Gaussian term of the ET-GM-PHD filter is , where is the number of Gaussian components. The process of dividing the living target measurement set is described as follows:
Select the measurement set of the surviving targets. The components with weights larger than 0.5 in the Gaussian term of time prediction are selected, which is denoted by
, and their position components are extracted as the centroid estimation of the surviving target, which is denoted by
,
. Then the candidate measurement set of the
ith extended target at time
k is given by the following:
Determine the class of the cluster. If , it means that the ith extended target at time k has died, the number of surviving extended targets is decreased by one, and the number of surviving targets at time k is finally determined as .
Cluster the living target measurements. For the set of measurements obtained in step 1, its union is . This dataset is put into the K-means++ clustering algorithm, and the ith cluster center is obtained as , , which can be regarded as the centroid estimation of the ith extended target. At this point, the tracking of the extended target is converted to the tracking of the point targets, and the estimated value of the centroid is brought into the update formula of GM-PHD, and the updated Gaussian term is .
3.3. Division of Newborn Target Measurement Set
It is assumed that the location of the newborn target and the living target are far apart in the sensor monitoring area, so that the measurements of the living target and the newborn target are less correlated. Under this assumption, the measurements of the newborn target do not exist in the set , so will only include the measurements of the newborn target (if there is a newborn target at the current time) and clutter measurements.
The steps of measurement set division for newborn targets are as follows:
- 1.
The measurement set is brought into the processing step of the DBSCAN clustering algorithm for preprocessing. The purpose of this step is to remove the clutter, where is the minimum number of measurements contained in the neighborhood () of the measurement. Furthermore, determine whether there is a cluster unit in the set; if so, go to step 2; otherwise, go to step 4.
- 2.
The set of measurements processed by the DBSCAN clustering algorithm is
, the number of clusters is
, and
only consists of measurements generated by extended targets. To this end, a novel SNN method based on directed k-nearest neighbors (kNN) graph is used to cluster, for
,
. The SNN of a directed kNN graph is defined as follows:
Following the above definition, the similarity between two centers of mass in a directed kNN graph is calculated as follows:
The specific method is described as follows:
Calculate the similarity matrix
W of the measurement set
according to the similarity defined in Equation (
5), i.e.,
Then, the dimension of W is , and W is an asymmetric matrix.
Spectral clustering is used to solve the graph partition problem for the calculated similarity matrix
W. First, we compute the unnormalized Laplacian matrix
L, as follows:
where
Calculate the normalized Laplacian matrix .
Build the matrix
U. Since the matrix
U is not a binary representation, we treat it as a new data matrix, which is denoted by the following:
where
K denotes the number of clusters,
,
,
,
,
n is the number of measurements in the measurement set
, and
is the measurement rate of the extended target.
Normalize each row of the matrix
U as follows:
is obtained. Each row in the matrix Y is clustered as a new data point in the K dimensional space, and the K-means++ clustering algorithm is used to obtain the clustering results . The measurement belongs to the jth class, when the ith row of the matrix Y is clustered into the jth class. The set of measurements is then partitioned into K classes, and . Then, a clustering result of the measurement set is obtained.
Repeat the previous step until the condition is not satisfied, and finally obtain the measurement division result at time k.
- 3.
For each partition result , the Gaussian term is updated in the update step of ET-GM-PHD filtering algorithm, and the updated Gaussian term is given by .
- 4.
The union of the updated Gaussian terms is . After pruning and merging, the target state at time k is extracted.
The proposed algorithm is shown in
Figure 1. Firstly, the measurement set is divided into the living target measurement set and the newborn target measurement set, and the computational complexity is
. For the surviving target measurement set, since the number of clusters required is known in advance, the computational complexity of K-means++ is
, where
is the number of surviving targets, and the computational complexity of GM-PHD update step is
, where
is the number of Gaussian components at the previous time. For the newly born target measurement set, the computational complexity is, at most,
. Suppose that the number of remaining measurements in the measurement set
after processing by DBSCAN algorithm is
N, the computational complexity of spectral clustering is, at most,
. When the number of measurements is large, the
k of the directed kNN graph is far less than
N, and its similarity matrix is the sparse matrix, so the calculation efficiency is high. Assume that the number of edges of the directed kNN graph is
m and the computational complexity of the eigenvector is
. The computational complexity of the ET-GM-PHD update step is
, where
denotes the number of partitions and
is the number of iterations, while the time complexity of distance partition is
.
Remark 1. When there are a large number of extended targets in the scene, and the extended targets are close to each other and move in parallel with variable speed, the distance division method only may fail to divide the measurement set. In order to improve the tracking performance, on the basis of the method proposed in this paper, multiple methods can be combined to further accurately divide the target measurement set, including gating, partitioning, assignment, etc. [6,20,21,22].
Figure 1.
Partitioning method based on spatio-temporal correlation.
Figure 1.
Partitioning method based on spatio-temporal correlation.
The flowchart of the computational complexity involved in each block of the proposed tracking method is presented in
Figure 2.
4. Simulation Results and Analysis
In order to verify the effectiveness of the proposed algorithm in this paper, this paper implements the ET-GM-PHD filtering algorithm, based on the algorithm proposed in this paper, distance partition [
14] and directed graph SNN partition [
18]. The simulation platform is MATLAB, R2018b, Win10 64-bit professional operating system,
[email protected] GHz, 16.00 GB RAM.
The bird’s eye view of the multi-object tracking scenario considered in this paper is shown in
Figure 3.
The evaluation method of OSPA [
23] is adopted to calculate the OSPA error of the filtering algorithm based on different division methods. In this simulation experiment, OSPA parameters are set as
and
, respectively. The number of Monte Carlo simulation experiments is 100. For fair comparison, the subsequent simulation experiments are the same as the simulation scenarios in the literature [
18]. In the simulation, the target kinematic model adopts the CV (constant velocity) model, and the state transition matrix is given by the following:
where the sampling interval
s in this simulation.
The process noise covariance matrix is , the measurement noise covariance matrix is , where is a identity matrix, the process noise standard deviation is set as , and the measurement noise standard deviation is set as .
The size of the area monitored by the sensor is m m. To verify the performance of the algorithm, the survival probability of each target is set as , with the detection probability of . The number of extended target measurements is a Poisson distribution with an expectation value of 10. The clutter measurements are uniformly distributed in the detection area, while the clutter count follows a Poisson distribution with an expected value of 20.
To ensure the feasibility of the algorithm, the Gaussian terms need to be pruned and merged. The maximum number of Gaussian components is , with the pruning threshold , the merging threshold , and the distance threshold .
4.1. Simulation Experiment 1
To verify the effectiveness of the proposed algorithm, two extended target cross motion scenes are constructed. The survival time of both extended target 1 and extended target 2 is 1∼100 s, and there are no derived targets in the scenario. The intensity function of the newborn target is given by the following:
where
[250 m, 250 m, 0 m/s, 0 m/s
,
[−250 m, −205 m, 0 m/s, 0 m/s
,
. The real tracks of the targets in the scene are shown in
Figure 4a, which shows that the tracks of the two extended targets cross at 56 s.
Figure 4b shows the running results of the proposed algorithm under a single simulation. It can be seen from the figure that the state estimation of the target is also accurate when the extended target tracks cross.
Figure 5a shows the OSPA distance comparison of different algorithms under 100 Monte Carlo simulations. It can be seen from the figure that the proposed algorithm OSPA distance is the minimum, even at the extended target track crossing moment. The OSPA distance gap between the algorithm, based on the directed graph partitioning of SNN, and the distance partitioning method is very small. Specifically, the OSPA distance of the directed SNN graph partition is significantly smaller in several adjacent moments when the targets cross. The comparison plot of OSPA distance verifies the effectiveness of the proposed algorithm.
Figure 5b gives a comparison plot of the extended target number estimation. It can be seen from the figure that, since the two targets intersect at 56 s, the extended target at the current moment is close to the extended target at several adjacent moments. Thus, when the measurement set is partitioned by distance partition and directed graph SNN partition methods, measurements of different extended targets fall into the same set, resulting in the potential estimation error, and the potential error of graph partition method is smaller. The method proposed in this paper takes into account the correlation of the extended targets at adjacent times, so that it can achieve accurate potential estimation when the targets are close to each other or cross.
Figure 6a,b shows the average number of partitions and the average running time of the measurement set at each time, respectively. It can be seen from the figures that distance partition has the highest average number of splits and the highest running time, because the distance partition needs to traverse multiple distance thresholds to find the best set of measurement partitions. The number of partitions of the directed graph SNN ranks second, and it first uses the DBSCAN algorithm to filter out part of the clutter, so the number of set partitions is smaller. The average number of partitions proposed in this paper is the lowest, which basically requires only one partition of the measurement set. When the extended target tracks cross, the measurement distance is close, and the number of distance division will be reduced at this time. At the same time, the running time of the algorithm reflects the number of partitions.
Figure 4.
(a) Real trajectory of the targets in Scenario 1. (b) Measurement values, true target trajectories, and estimated trajectories in Scenario 1.
Figure 4.
(a) Real trajectory of the targets in Scenario 1. (b) Measurement values, true target trajectories, and estimated trajectories in Scenario 1.
Figure 5.
(a) Measurement values, true target trajectories, and estimated trajectories in Scenario 1. (b) Comparison of OSPA errors under Scenario 1.
Figure 5.
(a) Measurement values, true target trajectories, and estimated trajectories in Scenario 1. (b) Comparison of OSPA errors under Scenario 1.
Figure 6.
(a) Comparison of the number of measurement set partitions in Scenario 1. (b) Comparison of running time of each algorithm for Scenario 1.
Figure 6.
(a) Comparison of the number of measurement set partitions in Scenario 1. (b) Comparison of running time of each algorithm for Scenario 1.
4.2. Simulation Experiment 2
The purpose of this experiment is to verify the filtering performance of the proposed algorithm for targets in scenarios where the extended targets are close to each other and moving in parallel. The survival time of extended target 1 and extended target 2 is both 1–100 s. The strength function of the newborn target in this scenario is given by the following:
where
[−600 m, −500 m, 0 m/s, 0 m/s
,
[−600 m, −600 m, 0 m/s, 0 m/s
.
The real target track of the target in scenario 2 is shown in
Figure 7a, which consists of two parallel moving extended targets that are close to each other.
Figure 7b shows the simulation results of the proposed algorithm. It can be seen from the figure that the kinematic state of the extended targets can be estimated relatively accurately even when the distance between the extended targets is close.
Figure 8a shows the OSPA distance comparison diagram of the three partitioning methods for scenario 2, and the OSPA distance index is calculated by the comprehensive calculation of the extended target state estimation and the target number estimation. It can be seen from the figure that the proposed algorithm has the smallest OSPA distance and the most accurate target number estimation at the same time. Scenes of parallel moving objects that are close to each other also verify the effectiveness of the algorithm. Meanwhile, the comparison figure of the target number estimation in
Figure 8b also reflects the effectiveness of the algorithm, and the target number estimation of the proposed algorithm is significantly better than that of the other two algorithms.
Figure 9a,b shows the average number of partitions and the average running time of the measurement set partition at each time, respectively. The figures show that the partition method proposed in this paper has the least number of partitions and the least running time. Since the proposed algorithm considers the relevance of the adjacent time of the extended target, the measurement set of the surviving target is divided by using the prior information, which greatly reduces the number of undivided measurements in the measurement set and the number of partitions.
Figure 8.
(a) Comparison of OSPA errors for Scenario 2. (b) Comparison of target number estimation in Scenario 2.
Figure 8.
(a) Comparison of OSPA errors for Scenario 2. (b) Comparison of target number estimation in Scenario 2.
4.3. Simulation Experiment 3
The goal of this experiment is to verify the effectiveness of the proposed algorithm for the birth and death scenarios of multiple extended targets. The strength function of the newborn target is given by the following:
where
100 m, 450 m, 0 m/s, 0 m/s
,
[−450 m, −450 m, 0 m/s, 0 m/s
,
450 m, −209 m, 0 m/s, 0 m/s
,
[−350 m, −350 m, 0 m/s, 0 m/s
.
Figure 10a shows the real tracks of multiple extended targets in Scenario 3.
Figure 10b shows the simulation results of the proposed algorithm.
Figure 11a shows the OSPA distance comparison of the three partition methods in scenario 3. It can be seen from the figure that the proposed method has the smallest OSPA distance. However, at the 27th and 76th seconds, the OSPA distance becomes larger, which is because there are newborn targets appearing at these two moments.
Figure 11b shows the comparison of target number estimation. As can be seen from the figure, the cardinality estimation of the proposed algorithm is significantly better than the other two algorithms when the number of targets in the scene is large.
Figure 12a,b show the average number of partitions and the average running time of the three methods. The proposed method has the lowest number of partitions and running time. Since the proposed method not only makes full use of the prior information, but also converts point target tracking into extended target tracking, the number of distance partitions increases with the number of objects in the scene.
Figure 10.
(a) Real target trajectory for Scenario 3. (b) Measurement values, true target trajectories, and estimated trajectories in Scenario 3.
Figure 10.
(a) Real target trajectory for Scenario 3. (b) Measurement values, true target trajectories, and estimated trajectories in Scenario 3.
Figure 11.
(a) Comparison of OSPA errors for Scenario 3 (b) Comparison of target number estimation in Scenario 3.
Figure 11.
(a) Comparison of OSPA errors for Scenario 3 (b) Comparison of target number estimation in Scenario 3.
Figure 12.
(a) Comparison of measurement set partition number in Scenario 3. (b) Comparison of running times of each algorithms for Scenario 3.
Figure 12.
(a) Comparison of measurement set partition number in Scenario 3. (b) Comparison of running times of each algorithms for Scenario 3.
4.4. Simulation Experiment 4
The purpose of this experiment is to verify the filtering performance of the proposed algorithm in the variable velocity case in case of variable velocity. We consider adding a constant turn (CT) model, i.e.,
where
,
denotes the rotational angular velocity, and
denotes the sampling period.
The strength function of the newborn target is given by the following:
where
The motion model for Scenario 4 is shown in
Table 2.
Figure 13a shows the real tracks of multiple extended targets in Scenario 4.
Figure 13b shows the OSPA distance comparison of the three partition methods in scenario 4. It can be seen from the figure that the proposed method has the smallest OSPA distance. However, at the 1st–5th seconds, the OSPA distance becomes larger, which is because there are newborn targets.
Figure 14a,b shows the average number of partitions and the average running time of the three methods. The proposed method has the lowest number of partitions and running time. The reason is that the proposed method not only makes full use of the prior information, but also converts point target tracking into extended target tracking.
Figure 14c shows the comparison of target number estimation. As can be seen from the figure, the cardinality estimation of the proposed algorithm is significantly better than the other two algorithms.
From the simulation results, we can see that the performance of the proposed algorithm is still better than the other two comparison methods in the variable velocity scenario with a CT model.
Figure 13.
(a) Real target trajectory for Scenario 4. (b) Comparison of OSPA errors in Scenario 4.
Figure 13.
(a) Real target trajectory for Scenario 4. (b) Comparison of OSPA errors in Scenario 4.
Figure 14.
(a) Comparison of measurement set partition number in Scenario 4. (b) Comparison of running times of each algorithms for Scenario 4. (c) Comparison of target number estimation in Scenario 4.
Figure 14.
(a) Comparison of measurement set partition number in Scenario 4. (b) Comparison of running times of each algorithms for Scenario 4. (c) Comparison of target number estimation in Scenario 4.