Label GM-PHD Filter Based on Threshold Separation Clustering

Gaussian mixture probability hypothesis density (GM-PHD) filtering based on random finite set (RFS) is an effective method to deal with multi-target tracking (MTT). However, the traditional GM-PHD filter cannot form a continuous track in the tracking process, and it is easy to produce a large number of redundant invalid likelihood functions in a dense clutter environment, which reduces the computational efficiency and affects the update result of target probability hypothesis density, resulting in excessive tracking error. Therefore, based on the GM-PHD filter framework, the target state space is extended to a higher dimension. By adding a label set, each Gaussian component is assigned a label, and the label is merged in the pruning and merging step to increase the merging threshold to reduce the Gaussian component generated by dense clutter update, which reduces the computation in the next prediction and update. After pruning and merging, the Gaussian components are further clustered and optimized by threshold separation clustering, thus as to improve the tracking performance of the filter and finally realizing the accurate formation of multi-target tracks in a dense clutter environment. Simulation results show that the proposed algorithm can form a continuous and reliable track in dense clutter environment and has good tracking performance and computational efficiency.


Introduction
Multi-target tracking (MTT) is a process of assigning the measured values to the targets, filtering them at the same time, and managing the tracks of multiple targets according to the time step [1,2]. RFS provides a unified and clear framework for MTT. Probability Hypothesis Density (PHD) [3] filtering has a good potential in solving target tracking under insufficient prior knowledge and an unknown number of targets. At present, it has been widely used in radar target tracking [4][5][6], computer vision [7], real-time positioning and map building [8], and group target tracking [9].
Traditional tracking algorithms, such as probabilistic data association (PDA) and multiple hypothesis tracking (MHT) [10][11][12][13], all transform the multi-target problem into a parallel single-target tracking problem by assigning measurements. The core of its processing method is data correlation, but when there are many targets and a large number of false alarm clutter, the correlation will bring a combination explosion and make the calculation amount increase exponentially. Correlation error and state estimation error are coupled and influence each other, resulting in large estimation error. PHD filtering technology avoids complex data association problems in the process of state estimation and can concentrate resources to deal with tracking problems. Thus, it has good application performance. Among them, Gaussian mixture (GM) and sequential Monte Carlo (SMC) are two important methods of PHD operation, which are called GM-PHD [14] and SMC-PHD [15], respectively. For the SMC single target algorithm, reference [16] proposed an adaptive model discrimination method, and does not need to be adjusted for specific problems, and has a better effect on the model uncertainty of discrete data. Reference [17] , where x n k and z m k represent the n-th target state and the m-th measurement at time k, respectively. N k and M k are the number of targets and the number of measurements at time k, respectively, and they change with time and the environment. Assuming that the prior probability of multi-target approximately obeys Poisson distribution, according to the RFS statistics theory, the PHD recurrence equation is [3] D k|k−1 (x) = p s,k|k−1 f k|k−1 (x|ζ) + β k|k−1 (x|ζ)D k−1|k−1 (ζ)dζ λc(z k )+ P D,k g k (z|ζ)D k|k−1 (ζ)dζ (2) γ k (x) and β k|k−1 (x|ζ) represent the target strength of newborn and spawned RFS respectively [31], p s,k|k−1 is the survival probability of target at time k on k − 1, P D,k represents the detection probability of target at time k, f k|k−1 (x|ζ) is the probability density function of state transition, g k (z|x) is the likelihood function of a single target, and λ is the average number of clutter, and c(z k ) is probability density functions of Poisson clutter distribution.

Gaussian Mixture Implementation
GM-PHD filter represents the prior PHD (1) and the posterior PHD (2) of multitarget as Gaussian mixture formation, and its iterative recursion can be expressed as a prediction update structure similar to Kalman filter (KF), and the algorithm is processed in a multi-target space.
Assume that the multi-target posterior PHD at time k − 1 can be expressed as a Gaussian mixture, and the equation is as follows: where J k−1 is the number of Gaussian components at time k − 1. In general, the number of J k is much greater than the number of N k . w i k−1 is the weight of the i-th Gaussian component, m i k−1 and p i k−1 represents the mean and covariance of the i-th Gaussian component. Then time k multi-target prediction PHD and time k posterior multi-target PHD can be expressed as Gaussian mixture form: where D β k|k−1 (x) and D sv k|k−1 (x) are the Gaussian mixture intensity function of spawned target and surviving target respectively, J k|k−1 is the number of Gaussian components of predicting PHD, and w i k|k−1 , m i k|k−1 and p i k|k−1 are respectively the weight, mean, and covariance of the i-th Gaussian component in the prediction intensity function. w j k , m j k , and p j k are, respectively, the weight, mean, and covariance of the j-th component in the update intensity function, F is the state transition matrix, H k is the measurement transition matrix, Q k is the process noise covariance, R k is the measurement noise covariance matrix, and I is the identity matrix.
The expected number of targetsN k|k−1 andN k associated with D k|k−1 and D k are obtained as follows:N Among them, J β,k and J γ,k , respectively, represent the number of spawned Gaussian components and newborn Gaussian components at time k.
Therefore, the terms required to represent D k at step k is However, the number of Gaussian components and the computational complexity increases very fast.

Label GM-PHD Filter
In Reference [14], a tracking label was introduced into the GM-PHD filter, which extracts the target track by assigning a label to each Gaussian component. Its algorithm includes the following steps: Step 1: initialize the Gaussian mixture posterior PHD at t = 0: Assign a unique label to each Gaussian component: τ i 0 is label of the i-th Gaussian component, and the mean is m i 0 and covariance is p i 0 .
Step 2: Predict according to Equation (4), the label set Γ k|k−1 construction method is as follows: The surviving Gaussian component still retains the original label, and the newborn or spawned Gaussian component assigns a newborn label, then the predicted label set is defined as: Among them, τ β is a spawned Gaussian component label and τ γ is a newborn Gaussian component label.
Step 3: Updating is performed according to Equation (7), and all Gaussian terms obtained after updating are assigned the same labels as their previously associated prediction terms.
Step 4: Pruning and merging is the PHD Gaussian component set updated at time k; Set a truncation threshold T r , a merging threshold T m , a maximum number of Gaussian terms J max , and two sets A and B.
Gaussian components weighted greater than the truncation threshold are put into set A: Then find out the index value of Gaussian term with the largest weight in set A: j = argmax i∈A ω i k , and make set B as If the two Gaussian components satisfy Equation (18), they are combined into one term. The combined component parameters are obtained by weighted summation: Step 5: Multi-target state extraction Besides extracting the target state, the set composed of Gaussian component labels is also extractedΓ Step 6: Track extraction At each moment, if there is a Gaussian component or its spawned Gaussian term ω i k > 0.5, it is marked as "confirmed" and other branches at ω i k ≤ 0.5 are marked as "tentative".

Robust Label GM-PHD Filter
In a dense clutter environment, due to the existence of a large number of clutters, a large number of Gaussian components with unknown sources will appear in the pruning and merging step. If not properly processed, the computational complexity of the label GM-PHD filter will be greatly increased in the subsequent prediction and update steps. The existence of dense clutter will also affect the final state extraction of the filter and then directly affect the formation of correct multi-target tracks. It is necessary to consider improving the pruning and merging method, and on this basis, optimize the state extraction process of the algorithm, thus that the output track is correct and stable, and at the same time improve the efficiency of the algorithm.

Improve Pruning and Merging Methods
As mentioned earlier, GM-PHD filter uses the Gaussian mixture distribution to model multi-target to update the posterior strength of multi-target. However, such as other PHD filters, this filter does not consider the time sequence correspondence of Gaussian components in the tracking process. Thus, it is impossible to distinguish targets. Label GM-PHD filtering method solves this problem, which uses pruning and merging process to exclude unreal targets and ensure the formation of tracks. However, the pruning and merging method does not consider the label set of Gaussian components. When the target was occluded, or the distance was close, the Gaussian components corresponding to the real target were pruned or merged, which significantly degrades the tracking performance of the algorithm. The main reason why targets cannot be tracked when they are close to each other is that the label of the Gaussian component is ignored. If the merging threshold T m is set relatively large, the number of merged Gaussian components will increase, which will lead to an increase in the probability of mixing different labels, thus that the 2 targets cannot be effectively distinguished. If the merge threshold setting T m is small, it is not conducive to the aggregation of real targets. A large number of invalid Gaussian components are still retained, which greatly increases the subsequent calculation burden. The above problems can be solved by the following improvements, that is, only Gaussian terms with the same label were allowed to be merged. This processing was not only beneficial to the problem of multi-target discrimination in dense clutter environments but can also consider setting a larger merging threshold T m to improve the target tracking performance. The expression is as follows: The combined component parameters were obtained by Equations (19)- (21). In order to reduce the operation time of each step, reference [31] proposes a tree structure, in which the initial Gaussian component was set as the origin, and if there were Gaussian components or their spawned Gaussian terms ω i k > 0.5, the tree was marked as "confirmed", and other branches unsatisfied the ω i k > 0.5 were marked as "tentative". In the process of target tracking, due to the influence of clutter and noise, the target may disappear for several frames before appearing. The confirmed Gaussian component n missed was assigned a parameter representing the number of consecutive missed detections on the branch. The method allows previously confirmed components n missed ≤ 3 marked as tentative in the current step.
The proposed algorithm flow is as follows: 1. For , the truncation threshold T r , the merging threshold T m and the maximum number of Gaussian terms J max 2.
Set the set Find the index value of the maximum weight value in set A Weighted summation to obtain Gaussian component parameters

Threshold Separation Clustering
Based on the improve pruning and merging methods, GM-PHD extracts Gaussian components with weight greater than 0.5 after pruning and merging as output and retains the remaining Gaussian components to continue iterative filtering. Due to the persistent interference of strong clutter in a dense clutter environment, it is easy to produce wrong Gaussian terms with large weights. If they are extracted and regarded as newborn targets, it will lead to high tracking errors.
It can be seen that finding the real target in a dense clutter environment and effectively separating the real target from the wrong Gaussian term are the key to solving the problem. Reference [32] proposed a threshold separation clustering method. On this basis, this paper proposes a method of clustering Gaussian components of the same label using velocity and distance information and extracting the target state, which can effectively aggregate the real target state and prevent redundant error. The Gaussian term interferes with the extraction of the target state and can reduce the amount of calculation on the basis of [32].
Firstly, the index value of the Gaussian component with the largest weight in set B was found, and the weight value was placed in the set S ω .
where N represents the number in set B after pruning and merging. The velocity difference and the corresponding Euclidean distance between the Gaussian component with the largest weight and the Gaussian component with the same label in set B were calculated.
wherein R 0 is used for the velocity information for separating, R 1 and R 2 is used for the position information for separating the target. In a linear system, R 0 is column matrices [0 1 0 1] T , and R 1 and R 2 are column matrices [1 0 0 0] T and [0 0 1 0] T , respectively. Finally, all Gaussian components satisfying the given velocity threshold M speed and Euclidean distance M d are clustered.
wherein the set S k contains all Gaussian components satisfying the conditions of Equations (27) and (28), and the clustering result is used as the final target state estimation.

Tracking Correlation to Form a Track
In the Gaussian component set S k obtained by threshold separation clustering, the Gaussian component with the largest weight in the corresponding label was found, which can effectively estimate the target track and introduce a track set k .
wherein ζ i k is the i-th track corresponding to time k, and r i k represents the corresponding track labeled τ i k . If there is a track with a label τ i k in the set k−1 and its parameter n missed = 0, m i k would be added to the labeled τ i k track and added to the track set. This Gaussian component is confirmed to match the tracking label τ i k . If there is a track in the set k−1 and has a label τ i k , but it is not added to the track set k , the number of missed detections is judged: If Gaussian component n i missed ≤ 3, m i k would be added to the labeled track τ i k . Finally, multi-target state estimation can be obtained:

Overall Steps to Improve the Algorithm
Aiming at the problem of poor tracking performance and low computational efficiency of traditional label GM-PHD in dense clutter environment, a robust label GM-PHD algorithm was proposed. The overall steps were as follows: Step 1: Initialize the parameters When t = 0, a Gaussian mixture posterior PHD and an initial label are set, such as Equations (35) and (36): Each Gaussian component is recursive by the Bayesian criterion.
Step 2: Predict step The Gaussian component is predicted according to Equation (4), and then newborn labels are assigned to the newborn and spawned Gaussian components, and the newborn label set and the spawned label set are added to the predicted label set: Step 3: Update the steps According to Equation (7), the prediction state is updated using the measurements at time k, and then all updated Gaussian terms are assigned the same label as their previously associated prediction terms.
Step 4: Pruning and merging A state-updated Gaussian mixture PHD ω i k , m i k , p i k J k i=1 is obtained; set a truncation threshold T r , a merging threshold T m , a maximum number J max of Gaussian terms, and two sets A and B.
The Gaussian component satisfying the truncation threshold was found and put into the set A, and then the Gaussian term with the largest weight j = argmax i∈A ω (i) k in the set A was found, and the Gaussian term was used for screening and merging: The combined component parameters were obtained by weighted summation: The Gaussian component of Set A that occurs in both Set A and Set B was deleted until Set A was an empty set.
Step 5: Threshold separation clustering A velocity threshold and an Euclidean distance threshold were set, and Gaussian components satisfying Equations (43) and (44) and having the same label were clustered, and then Gaussian components satisfying conditions were put into a set S k , and the clustering criteria were as follows: Step 6: Tracking association to form track If there is a track with a label τ i k in the set k−1 and its parameters n missed = 0, the corresponding Gaussian term m i k is added to the labeled track τ i k , and it is added to the track set k , at which time the Gaussian component is confirmed to match the tracking label τ i k . If there is a track in the set k−1 with a label τ i k , but it is not added to the track set k , the number of missed detections should be judged: If Gaussian component n i missed ≤ 3, m (i) k would be added to the labeled track τ i k .

Experimental Parameter Setting
In order to verify the effectiveness of the proposed algorithm in MTT scenarios, the improved label GM-PHD filter was compared with GM-PHD filter, label GM-PHD filter, and Gaussian mixture cardinalized probability hypothesis density (GM-CPHD) [28] filter. The experimental parameters were set as follows.
The tracking scene was set to multiple targets in four possible locations or spawned from other targets, and the observation area was [−1000, 1000] × [−1000, 1000] m 2 . There were 12 batches of targets in the scene. For the sake of simplicity, it was assumed that each target moves in a straight line at a uniform speed.
The state vector of the target consists of position and velocity components: x k = p x,k p y,k v x,k v y,k , and its state equation is: The sampling interval T is 1 s, the total tracking time was 100 s, and the noisy process is w k ∼ N(0, 5). The intensity of newborn target was as follows (47) γ,k = 0.03. The process noise of newborn targets obeys Gaussian distribution, the mean value was zero, and the covariance was: Q i sp,k = diag([100, 100, 100, 100]). The measurement vector is position information: z k = p z x ,k p z y ,k , and the measurement equation is given by the following equation: where measured noise was v k ∼ N(0, 5). Clutter is evenly distributed in the moving area of the target, with average scanning λ. In the Gaussian component pruning section, the truncation threshold of the Gaussian component is 10 −5 . The state extraction threshold was set to 0.5. The merging threshold was set to 10, the maximum number of Gaussian components was 100, the velocity threshold M speed = 10, and the Euclidean distance M d = 40. The number of Monte Carlo simulations was 100. Evaluating tracking quality by OSPA distance: x k is the target state vector, the two parameters of OSPA distance were set to p = 1 and c = 200, respectively. The smaller the OSPA distance, the higher the accuracy of target state estimation.

Simulation Scenario 1
In this MTT scene, clutter rate λ = 90, survival probability p s,k = 0.99, and detection probability p D,k = 0.95. Figure 1 is the real track and actual tracking measurement map of multiple targets in dense clutter and high detection probability environment. As can be seen from Figure 1b, in the MTT scene with dense clutter, the proposed improved GM-PHD algorithm can accurately obtain the target track. In Figure 1b, we can see that there are false tracks, which are caused by the clutter around the measured values of newborn targets. They will affect the stability of the filter before the track is formed. However, after clustering by speed and distance, the influence would gradually disappear with time, and the wrong target track would not appear. map of multiple targets in dense clutter and high detection probability environment. As can be seen from Figure 1b, in the MTT scene with dense clutter, the proposed improved GM-PHD algorithm can accurately obtain the target track. In Figure 1b, we can see that there are false tracks, which are caused by the clutter around the measured values of newborn targets. They will affect the stability of the filter before the track is formed. However, after clustering by speed and distance, the influence would gradually disappear with time, and the wrong target track would not appear.  Figure 2 compares the performance of various filters in dense clutter and high detection probability environment, including GM-PHD filter, label GM-PHD filter, robust label GM-PHD filter, and GM-CPHD filter. Obviously, the tracking performance of the labeled GM-PHD filter proposed in this paper was obviously better than GM-PHD, traditional labeled GM-PHD, and GM-CPHD in dense clutter environments. As can be seen from Figure 2a,b, the OSPA distance between GM-PHD and label GM-PHD filters is large, which shows that the dense clutter environment seriously affects the filter performance. As can be seen from Figure 2c,d, the improved labels GM-PHD and GM-CPHD have better performance in estimating the number of targets because the improved label GM-PHD filter combines Gaussian components with the same label and adds a threshold clustering separation step after merging and pruning, while GM-CPHD filter calculates the second moment of a posterior probability density function at the expense of operation time. Thus it is more accurate to estimate the number of multiple targets.  Figure 2 compares the performance of various filters in dense clutter and high detection probability environment, including GM-PHD filter, label GM-PHD filter, robust label GM-PHD filter, and GM-CPHD filter. Obviously, the tracking performance of the labeled GM-PHD filter proposed in this paper was obviously better than GM-PHD, traditional labeled GM-PHD, and GM-CPHD in dense clutter environments. As can be seen from Figure 2a,b, the OSPA distance between GM-PHD and label GM-PHD filters is large, which shows that the dense clutter environment seriously affects the filter performance. As can be seen from Figure 2c,d, the improved labels GM-PHD and GM-CPHD have better performance in estimating the number of targets because the improved label GM-PHD filter combines Gaussian components with the same label and adds a threshold clustering separation step after merging and pruning, while GM-CPHD filter calculates the second moment of a posterior probability density function at the expense of operation time. Thus it is more accurate to estimate the number of multiple targets.

Computational Complexity Analysis
In the process of pruning and merging the traditional GM-PHD filter and the label GM-PHD filter, all Gaussian components were pruned and merged. The algorithm proposed in this paper only processes the Gaussian components with the same label. It will undoubtedly greatly reduce the amount of calculation of the Gaussian component and avoid the problem of exponential increase in the amount of calculation. Moreover, the algorithm proposed in this paper also includes a threshold separation clustering method, which also processes Gaussian components with the same label. It can be seen in the following simulation that processing the Gaussian components with the same label will reduce the computational complexity very well. Figure 3 compares the running time of different filters in this scenario. It can be seen that the improved label GM-PHD filtering algorithm improves computational efficiency. Compared with the GM-PHD filtering algorithm, it has slightly higher computational complexity but compared with the other two filtering algorithms, it has higher computational efficiency because it uses label sets to combine Gaussian components with the same label, thus reducing the computational complexity. Figure 3 also shows that the improvement of GM-CPHD tracking performance is at the expense of computing time to some extent.
Compared with the GM-PHD filter, the operation time of the improved label GM-PHD filter is still 13.14% higher than that of the GM-PHD filter, but the computational efficiency is improved by 11.23% compared with the traditional label GM-PHD filter.

Computational Complexity Analysis
In the process of pruning and merging the traditional GM-PHD filter and the label GM-PHD filter, all Gaussian components were pruned and merged. The algorithm proposed in this paper only processes the Gaussian components with the same label. It will undoubtedly greatly reduce the amount of calculation of the Gaussian component and avoid the problem of exponential increase in the amount of calculation. Moreover, the algorithm proposed in this paper also includes a threshold separation clustering method, which also processes Gaussian components with the same label. It can be seen in the following simulation that processing the Gaussian components with the same label will reduce the computational complexity very well. Figure 3 compares the running time of different filters in this scenario. It can be seen that the improved label GM-PHD filtering algorithm improves computational efficiency. Compared with the GM-PHD filtering algorithm, it has slightly higher computational complexity but compared with the other two filtering algorithms, it has higher computational efficiency because it uses label sets to combine Gaussian components with the same label, thus reducing the computational complexity. Figure 3 also shows that the improvement of GM-CPHD tracking performance is at the expense of computing time to some extent. Compared with the GM-PHD filter, the operation time of the improved label GM-PHD filter is still 13.14% higher than that of the GM-PHD filter, but the computational efficiency is improved by 11.23% compared with the traditional label GM-PHD filter.

Simulation Scenario 2
In this MTT scene, clutter rate λ = 90, survival probability p s,k = 0.99, and detection probability p D,k = 0.8. Figure 4 shows the tracking performance comparison of each filter in dense clutter and low detection probability environment. As can be seen from Figure 4, except for the improved label GM-PHD filtering algorithm, the tracking performance of other filtering algorithms is obviously degraded under the condition of low detection probability, which indirectly proves the robustness of the proposed algorithm. After pruning and merging, the Gaussian components are separated and clustered by using velocity and distance thresholds, which effectively improves the tracking performance of multi-targets under the environment of low detection probability.

Simulation Scenario 3
In this MTT scene, the survival probability

Simulation Scenario 3
In this MTT scene, the survival probability p s,k = 0.99, detection probability p D,k = 0.9, the proposed algorithm under different clutter rate tracking performance superiority was verified by simulation. Figure 5 shows the performance comparison of algorithms with different clutter rates. As can be seen from Figure 5a-c, the tracking errors of the GM-PHD filter and the label GM-PHD filter increases linearly with the increasing clutter rate. The reason is that with the increasing clutter, the Gaussian component pruning and combining step becomes difficult to handle, resulting in error coupling of various errors. The improved label GM-PHD filtering algorithm uses the Gaussian component of the same label to process the subsequent steps, which effectively avoids these problems. It can be seen from Figure 5d that with the increase of clutter rate, the calculation amount of GM-CPHD increases linearly, while the calculation amount of the improved label GM-PHD filtering algorithm changes little, which also reflects the stability of the algorithm.

Simulation Scenario 4
In this MTT scene, the survival probability λ were simulated to prove the performance advantages of the proposed algorithm in different detection probability environments. Figure 6 shows the comparison of algorithm performance under different detection probabilities. It can be seen from Figure 6 that the improved label GM-PHD filtering algorithm had higher stability than other algorithms and had certain advantages in tracking performance under different detection probabilities.

Simulation Scenario 4
In this MTT scene, the survival probability p s,k = 0.99 and clutter rate λ= 60 were simulated to prove the performance advantages of the proposed algorithm in different detection probability environments. Figure 6 shows the comparison of algorithm performance under different detection probabilities. It can be seen from Figure 6 that the improved label GM-PHD filtering algorithm had higher stability than other algorithms and had certain advantages in tracking performance under different detection probabilities.  Simulations of the above four scenarios show that the improved label GM-PHD filtering algorithm has good tracking performance in dense clutter and low detection probability environment. It also has good stability in the environment where clutter rate and detection probability change.
This algorithm is based on GM-PHD and focuses on multi-target tracking in a dense clutter environment, and the motion of the target is also linear. In the next work, we can study non-linear high-maneuvering multi-targets and Multi-sensor observation scenes, using SMC algorithm for non-linear moving targets. By considering the tracking problems in the case of model uncertainty, the algorithms [16][17][18] can be more suitable for actual tracking scenarios.

Conclusions
Aiming at the problem of low accuracy and low computational efficiency of MTT in dense clutter environments, this paper proposes an efficient and accurate robust label GM-PHD filtering algorithm. Based on the traditional label GM-PHD, this algorithm improves the pruning and merging step to merge Gaussian components with the same label, which greatly reduces the computational burden and error merging in a dense clutter environment. According to velocity threshold and Euclidean distance, the newborn state Simulations of the above four scenarios show that the improved label GM-PHD filtering algorithm has good tracking performance in dense clutter and low detection probability environment. It also has good stability in the environment where clutter rate and detection probability change.
This algorithm is based on GM-PHD and focuses on multi-target tracking in a dense clutter environment, and the motion of the target is also linear. In the next work, we can study non-linear high-maneuvering multi-targets and Multi-sensor observation scenes, using SMC algorithm for non-linear moving targets. By considering the tracking problems in the case of model uncertainty, the algorithms [16][17][18] can be more suitable for actual tracking scenarios.

Conclusions
Aiming at the problem of low accuracy and low computational efficiency of MTT in dense clutter environments, this paper proposes an efficient and accurate robust label GM-PHD filtering algorithm. Based on the traditional label GM-PHD, this algorithm improves the pruning and merging step to merge Gaussian components with the same label, which greatly reduces the computational burden and error merging in a dense clutter environment. According to velocity threshold and Euclidean distance, the newborn state extraction method can effectively reduce the occurrence of false tracks and improve the performance of MTT in dense clutter environments. The proposed algorithm can be further improved to adapt to MTT in an unknown parameter environment.
Author Contributions: K.W. performed the simulations and wrote the paper; Q.Z. and X.H. revisedthis paper and offered some useful suggestions in methodologies and simulations. All authors have read and agreed to the published version of the manuscript.