Next Article in Journal
Multi-Domain Airflow Modeling and Ventilation Characterization Using Mobile Robots, Stationary Sensors and Machine Learning
Previous Article in Journal
A Vibro-Acoustic Hybrid Implantable Microphone for Middle Ear Hearing Aids and Cochlear Implants
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Target Birth Intensity Multi-Bernoulli Filter with Noise-Based Threshold

School of Electronic Engineering, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(5), 1120; https://doi.org/10.3390/s19051120
Submission received: 23 January 2019 / Revised: 27 February 2019 / Accepted: 27 February 2019 / Published: 5 March 2019
(This article belongs to the Section Physical Sensors)

Abstract

:
Adaptively modeling the target birth intensity while maintaining the filtering efficiency is a challenging issue in multi-target tracking (MTT). Generally, the target birth probability is predefined as a constant and only the target birth density is considered in existing adaptive birth models, resulting in deteriorated target tracking accuracy, especially in the target appearing cases. In addition, the existing adaptive birth models also give rise to a decline in operation efficiency on account of the extra birth modeling calculations. To properly adapt the real variation of the number of newborn targets and improve the multi-target tracking performance, a novel fast sequential Monte Carlo (SMC) adaptive target birth intensity cardinality balanced multi-target multi-Bernoulli (CBMeMBer) filter is proposed in this paper. Through adaptively conducting the target birth probability in a pre-processing step, which incorporates the information of current measurements to correct the pre-setting of the target birth probability, the proposed filter can truly adapt target birth cases and achieve better tracking accuracy. Moreover, the implementation efficiency can be improved significantly by employing a measurement noise-based threshold in the likelihood calculations of the multi-target updating. Simulation results verify the effectiveness of the proposed filter.

1. Introduction

The core objective in multi-target tracking (MTT) is to estimate the states of multiple moving targets, based on obtained data with spontaneous appearance and disappearance of different targets. Through generalizing the single-target recursive filter to the multi-target tracking situation, the random finite set (RFS)-based Bayesian filter [1] was proven to be appropriate for multi-target tracking. To reduce the algorithm complexity of the RFS Bayes multi-target filter, a series of solutions such as the probability hypothesis density (PHD) filter [2,3], the cardinalized PHD (CPHD) filter [4,5], the cardinality balanced multi-target multi-Bernoulli (CBMeMBer) filter [6,7], and the generalized labeled multi-Bernoulli (GLMB) filter [8,9,10] were developed. The PHD filter propagates the first moment of the multi-target density, and the CPHD filter is a generalization of the PHD filter that jointly propagates the moment and cardinality distribution to provide better tracking accuracy. The CBMeMBer filter propagates the multi-Bernoulli distribution to approximate the multi-target density, and the GLMB filter is a generalization of the CBMeMBer filter that propagates a generalized labeled multi-Bernoulli distribution to jointly provide better performance and target trajectories. Compared to the PHD and CPHD filters which are synonymous, the multi-Bernoulli representations have advantages in cases that ask for individual target existence probabilities or sequential Monte Carlo (SMC) implementations. In the SMC multi-Bernoulli filters, the clustering algorithm to extract the states of targets can be avoided, while it is required in the SMC PHD and CPHD filters.
The RFS-based filters were widely applied in various tracking scenarios. In these scenarios, environment descriptions, e.g., the newborn target [11,12,13,14,15], the target motion models [16,17], the unknown target [18,19], the unknown process noise statistics [20,21], the unknown measurement noise statistics [22,23], the target measurement models [24,25], and the detection and clutter rate uncertainty [26,27], are of great importance. It is worth noting that the description of newborn targets attracted intensive attention over the recent years, which is typically called the target birth modeling. “Birth” means that a target appears in the observation region. In the RFS-based filters, the target birth model is expressed as an intensity component [1] and applied at each filtering step. The target birth intensity consists two parts: the target birth density and the target birth probability, which represent the positions and the corresponding possibilities that newborn targets may appear, respectively. Modeling the target birth intensity is to capture the targets that enter the field of view (FOV) in real time. However, the indeterminacy of target birth makes the problem intractable.
So far, most researches on the RFS-based Bayesian filters required a priori location information of newborn targets and took the target birth intensity as fixed values. However, in practice, such as in the area of air target tracking, if the positions of airports are unavailable or the targets emerge from other positions, serious errors may occur. Meanwhile, the rough assumption fixing the number of newborn targets as a constant is improper and may cause severe degradation of the estimation accuracy. In the existing literature, adaptive target birth modeling is rarely discussed. In Reference [11], an adaptive target birth model was applied in the PHD and CPHD filters by grouping the targets into newborn and surviving sets, whose target birth density was generated through the measurements of the previous time step. Originated from Reference [11], the adaptive birth model was applied in the CBMeMBer and GLMB filters [12,13], respectively, using a simplified structure. In Reference [14], an evaluating step of multiple measurements and states including newborn targets was proposed. In Reference [15], a newborn track detection and state estimation method was proposed using the Bernoulli RFS based on the sequential change in measurements. However, with an aim to remove the requirement of the prior location information about newborn targets, the existing adaptive birth modeling methods focus on the target birth density, while still setting the target birth probability as a constant. On the other hand, compared with the traditional fixed birth model, the time consumption of the filters applying adaptive birth models increases dramatically, since these birth models are more complex and need to be integrated into the filtering recursion at every time step.
In this paper, we propose an improved multi-target Bayesian filter based on the CBMeMBer framework which can adapt unknown target birth cases including both density and probability. The target birth density is modeled by the measurements received at the previous moment, which is similar to the method in Reference [13], while the target birth probability is calculated by a pre-processing step using the current received measurements. This CBMeMBer filter with the new adaptive target birth intensity not only eliminates the requirement of prior birth information, but also avoids the coarse assumption of the number of newborn targets. Thus, the proposed filter can achieve a better estimation accuracy and can be used in more tracking scenarios. To reduce the heavy computational burden incurred by the adaptive target birth intensity, a threshold-based method is introduced, which can eliminate most of the redundant measurement-updating calculations. Different from the widely used thresholds for single-target tracking [28,29] or requiring target–measurement association [29,30,31], the proposed threshold is determined by the measurement noise and operates on the distance between the state of particles and measurements, without any requirements for track continuity. The improved adaptive target birth intensity CBMeMBer filter with the threshold method is shown to be able to achieve better tracking performance with reduced computational burden.
The remainder of this paper is organized as follows: the fundamentals of the standard CBMeMBer filter are reviewed in Section 2. The adaption of the CBMeMBer filter to accommodate the unknown target birth intensity, and its fast filtering recursion with a novel threshold are elaborated on in Section 3 and Section 4, respectively. Simulation results are presented in Section 5, and conclusions are given in Section 6.

2. Cardinality Balanced Multi-Target Multi-Bernoulli Filter

This section reviews the standard CBMeMBer filter, which is a full Bayesian filter within the multi-Bernoulli framework and was used in many fields such as sensor networks [32], audio data [33], visual tracking [34], etc.
We assume that, at time k , the RFSs [1] of the multi-target state X k and the multi-target observation Z k in the state and observation space [1] χ and Z can be written as
X k = { x k , 1 , , x k , N ( k ) } χ ,
Z k = { z k , 1 , , z k , M ( k ) } Z .
The recursions of the multi-target posterior density are based on the framework of the multi-target Bayesian filter [1], where the RFSs are used to express the set-valued random variables.
In the CBMeMBer filters, Bernoulli and multi-Bernoulli RFSs are basic elements. A Bernoulli RFS [6] X means a dyadic distribution. With probability r and density p , X can be written as
π ( X ) = { 1 r X = r p ( x ) X = { x } .
A multi-Bernoulli RFS X = i = 1 M X ( i ) represents an integration of independent Bernoulli RFSs X ( i ) , which can be written as { ( r ( i ) , p ( i ) ) } i = 1 M . The sum of the existence probability r ( i ) describes the cardinality of the multi-Bernoulli RFS. A detailed introduction can be found in References [1,6].
The CBMeMBer filter [6] is derived based on the following assumptions:
  • The evolution of each target and the generation of each observation are all independent;
  • The clutter is independent of the observations of targets and follows a Poisson distribution;
  • One target can generate at most one observation at each scan;
  • Target birth is multi-Bernoulli and is independent of target survival.
The propagation of the CBMeMBer filter is presented below.
Suppose that the posterior multi-target density at time k 1 follows a multi-Bernoulli distribution given by
π k 1 = { ( r k 1 ( i ) , p k 1 ( i ) ) } i = 1 M k 1 .
Then, the predicted multi-target density is derived by combining the multi-Bernoulli distributions representing the newborn targets and the surviving targets, which can be written as
π k | k 1 = { ( r Γ , k ( i ) , p Γ , k ( i ) ) } i = 1 M Γ , k { ( r P , k | k 1 ( i ) , p P , k | k 1 ( i ) ) } i = 1 M k 1 ,
where { ( r Γ , k ( i ) , p Γ , k ( i ) ) } i = 1 M Γ , k is the prediction of newborn targets at time k , and { ( r P , k | k 1 ( i ) , p P , k | k 1 ( i ) ) } i = 1 M k 1 is the prediction of surviving targets at time k , with
r P , k | k 1 ( i ) = r k 1 ( i ) p k 1 ( i ) , p S , k ,
p P , k | k 1 ( i ) ( x ) = f k | k 1 ( x | ) , p k 1 ( i ) p S , k p k 1 ( i ) , p S , k ,
where p S , k ( η ) and f k | k 1 ( | η ) are the survival probability and the Markov transition density. The total number of predicted Bernoulli distributions is given by M k | k 1 = M k 1 + M Γ , k . Therefore, the predicted multi-target density at time k is still a multi-Bernoulli distribution with
π k | k 1 = { ( r k | k 1 ( i ) , p k | k 1 ( i ) ) } i = 1 M k | k 1 .
After the prediction obtained, the approximated posterior multi-target density can be written as
π k { ( r L , k ( i ) , p L , k ( i ) ) } i = 1 M k | k 1 { ( r U , k ( z ) , p U , k ( ; z ) ) } z Z k ,
where the legacy tracks are expressed by
r L , k ( i ) = r k | k 1 ( i ) 1 p k | k 1 ( i ) , p D , k 1 r k | k 1 ( i ) p k | k 1 ( i ) , p D , k ,
p L , k ( i ) ( x ) = p k | k 1 ( i ) ( x ) 1 p D , k ( x ) 1 p k | k 1 ( i ) , p D , k ,
with detection probability p D , k ( x ) . The measurement update tracks are expressed by
r U , k ( z ) = i = 1 M k | k 1 r k | k 1 ( i ) ( 1 r k | k 1 ( i ) ) p k | k 1 ( i ) , ψ k , z ( 1 r k | k 1 ( i ) p k | k 1 ( i ) , p D , k ) 2 κ k ( z ) + i = 1 M k | k 1 r k | k 1 ( i ) p k | k 1 ( i ) , ψ k , z 1 r k | k 1 ( i ) p k | k 1 ( i ) , p D , k ,
p U , k ( x ; z ) = i = 1 M k | k 1 r k | k 1 ( i ) 1 r k | k 1 ( i ) p k | k 1 ( i ) ( x ) ψ k , z ( x ) i = 1 M k | k 1 r k | k 1 ( i ) 1 r k | k 1 ( i ) p k | k 1 ( i ) , ψ k , z ,
with
ψ k , z ( x ) = g k ( z | x ) p D , k ( x ) ,
where g k ( | x ) is the measurement likelihood function, and κ k ( ) is the clutter intensity.

3. Extension of the Cardinality Balanced Multi-Target Multi-Bernoulli Filter

In the standard CBMeMBer filter of Section 2, both the target birth probability r Γ , k ( i )  and the target birth density p Γ , k ( i )  in Equation (5) are a priori assumed to be constant and independent of the measurements. However, the prior information about newborn targets is difficult to be obtained since the targets always appear irregularly over time in practice; for example, in monitoring, targets may appear from anywhere in the observation region at any time step, which makes the crude birth density and probability models invalid. In recent years, some investigations were developed for modeling the birth density [11,12,13,14,15], while the birth probability is still simply set as constant. This section firstly proposes a CBMeMBer filter with a completely measurement-driven target birth intensity, and then presents its SMC implementation.

3.1. CBMeMBer Filter Using Adaptive Target Birth Intensity

This subsection presents a new target birth intensity adaptively modeling both the target birth density [13] and the target birth probability. Figure 1 shows the framework of the adaptive target birth intensity modeling. The derivation of the improved adaptive target birth intensity CBMeMBer filter is presented below.
(1) Target birth intensity modeling
Similar to Reference [13], at each filtering step, the target birth density p Γ , k ( i ) , i = 1 , , M Γ , k is generated using measurements of the previous time step; thus, the target birth intensity { ( r Γ , k ( i ) , p Γ , k ( i ) ) } i = 1 M Γ , k can be re-expressed by
π Γ , k = { ( r Γ , k 1 ( i ) , p Γ , k 1 ( i ) ( x | z i ) ) } i = 1 | Z k 1 | ,
where the target birth probability r Γ , k 1 ( i ) is assigned equally by B Γ , k / | Z k 1 | , with B Γ , k being the expected number of target births and is usually set as a constant. This coarse determination of the expected number of newborn targets and the target birth probability is not always satisfied and can be problematic in certain cases. The birth probability that can adapt the background is quite demanding when capturing a real case. Below, a measurement-driven pre-processing method is proposed to properly adapt the real variation of the number of newborn targets.
Based on the target birth density in Equation (15), which removes the requirement for the prior position information of newborn targets, the corresponding target birth probability r Γ , k 1 ( i ) is re-calculated using the current measurements instead of being uniformly assigned as B Γ , k / | Z k 1 | .
Suppose that, at time k 1 , the posterior density of the CBMeMBer filter is a multi-Bernoulli distribution with
π k 1 = { ( r k 1 ( i ) , p k 1 ( i ) ) } i = 1 M k 1 .
Proposition 1.
In the pre-processing step, an experimental expectation B ^ Γ , k on newborn target number, e.g., B ^ Γ , k = 0.2 , is assumed and assigned to each target birth density with an identical probability r ^ Γ , k 1 ( i ) = B ^ Γ , k / | Z k 1 | . Applying the standard CBMeMBer prediction step, the predicted multi-target density of surviving targets is obtained as
π P , k | k 1 = { ( r P , k | k 1 ( i ) , p P , k | k 1 ( i ) ) } i = 1 M k 1 ,
and the experimental target birth intensity,
π ^ Γ , k = { ( r ^ Γ , k 1 ( i ) , p Γ , k 1 ( i ) ( x | z i ) ) } i = 1 | Z k 1 | ,
can be seen as a surmise of the initiation of the newborn target tracks. By correcting the experimental target birth intensity π ^ Γ , k through the prediction of the surviving targets and the current measurements z Z k , a modified target birth probability r ˜ Γ , k 1 ( i ) , i = 1 , , | Z k 1 | can be derived. Each of the modified target birth probabilities is calculated separately with full consideration of the detection of the new-born target, which can be given by
r ˜ Γ , k 1 ( i ) = r ˜ L Γ , k 1 ( i ) + r ˜ U Γ , k 1 ( i ) ,
where r ˜ L Γ , k 1 ( i ) and r ˜ U Γ , k 1 ( i ) represent the modeling of the birth probability under newborn targets being missed and detected by sensors, respectively. r ˜ L Γ , k 1 ( i ) can be expressed by
r ˜ L Γ , k 1 ( i ) = r ^ Γ , k 1 ( i ) 1 p Γ , k 1 ( i ) ( | z i ) , p D , k 1 r ^ Γ , k 1 ( i ) p Γ , k 1 ( i ) ( | z i ) , p D , k ,
which is derived similar to the legacy tracks in the standard CBMeMBer filter, and r ˜ U Γ , k 1 ( i ) can be expressed by
r ˜ U Γ , k 1 ( i ) = z Z k G k ( z ) κ k ( z ) + S Γ k ( z ) + S U k ( z ) ,
which involves operations traversing all the measurements, with
G k ( z ) = r ^ Γ , k 1 ( i ) ( 1 r ^ Γ , k 1 ( i ) ) p Γ , k 1 ( i ) ( | z i ) , ψ k , z ( 1 r ^ Γ , k 1 ( i ) p Γ , k 1 ( i ) ( | z i ) , p D , k ) 2 ,
S Γ k ( z ) = i = 1 | Z k 1 | r ^ Γ , k 1 ( i ) p Γ , k 1 ( i ) ( | z i ) , ψ k , z 1 r ^ Γ , k 1 ( i ) p Γ , k 1 ( i ) ( | z i ) , p D , k ,
S U k ( z ) = i = 1 M k 1 r P , k | k 1 ( i ) p P , k | k 1 ( i ) , ψ k , z 1 r P , k | k 1 ( i ) p P , k | k 1 ( i ) , p D , k .
G k ( ) can be seen as a function that reflects the correlation between the birth component and the measurements. S Γ k ( ) and S U k ( ) can be seen as normalization factors. With Equations (19)–(24), the modified target birth probability r ˜ Γ , k 1 ( i ) can be obtained, and explained as a rational assumption of the average number of newborn targets appearing in the corresponding birth area.
However, r ˜ Γ , k 1 ( i ) may be overestimated due to the clutter or the surviving targets, which can be mitigated by a hard-limited upper bound r max . As a result, the required target birth probability r Γ , k 1 ( i ) can be obtained as
r Γ , k 1 ( i ) = min ( r ˜ Γ , k 1 ( i ) , r max ) , i = 1 , , | Z k 1 | .
Table 1 summarizes the modeling of the target birth probability via pseudo code.
Next, we integrate the improved birth probability into the normal filtering procedure.
(2) Prediction
Upon substituting the experimental target birth probability r ^ Γ , k 1 ( i ) , i = 1 , , | Z k 1 | with r Γ , k 1 ( i ) , i = 1 , , | Z k 1 | , the new target birth intensity model { ( r Γ , k 1 ( i ) , p Γ , k 1 ( i ) ( | z i ) ) } i = 1 | Z k 1 | can be obtained and applied in the subsequent filtering step. Since the predicted density of surviving targets π P , k | k 1 is irrelevant to that of newborn targets and was computed in the pre-processing step, the re-calculation of the predicted surviving target density can be avoided in the normal filtering. Thus, the predicted density can be presented by
π k | k 1 = { ( r Γ , k 1 ( i ) , p Γ , k 1 ( i ) ( x | z i ) ) } i = 1 | Z k 1 | { ( r P , k | k 1 ( i ) , p P , k | k 1 ( i ) ) } i = 1 M k 1 ,
and can be further expressed as Equation (8) with M k | k 1 = | Z k 1 | + M k 1 .
(3) Update
Following that, the normal update step proceeds similarly to the standard CBMeMBer filter, with the new adaptive target birth intensity model involved. In updating, the calculation of the normalization function in Equation (24) can be preserved and directly imported when it is needed into Equation (12). As a result, the posterior density is obtained and the full recursion of time k is completed.
Remark 1.
Note that the target birth density is determined by the measurements of previous time steps; thus, misrecognition of clutter can be avoided. Moreover, the target birth probability is calculated in the pre-processing step without prior information requirements of the number of newborn targets. As a result, compared to the CBMeMBer filter that only adopts the adaptive target birth density, the proposed filter can achieve a better tracking accuracy.

3.2. Implementation

This subsection presents the SMC implementation of the proposed CBMeMBer filter. The SMC implementation uses particle sets to approximate the densities in filters, which relieves the restriction of linearity and Gaussian assumptions.
(1) Target birth intensity modeling
At time k 1 , the posterior density of the target is given by Equation (16), where
p k 1 ( i ) ( x ) = j = 1 L k 1 ( i ) w k 1 ( i , j ) δ x k 1 ( i , j ) ( x ) .
The probability density p k 1 ( i ) ( x ) is expressed by the weighted particle set { w k 1 ( i , j ) , x k 1 ( i , j ) } j = 1 L k 1 ( i ) , where L k 1 ( i ) is the number of particles.
Target birth density: The target birth density p Γ , k 1 ( i ) ( | z i ) of the newborn target intensity { ( r Γ , k 1 ( i ) , p Γ , k 1 ( i ) ( | z i ) ) } i = 1 | Z k 1 | is expressed by a weighted particle set x Γ , k 1 ( i , j ) , j = 1 , , L Γ , k 1 ( i ) , generated from measurement z i Z k 1 . Note that the target state can be written as x = [ p T , v T ] T , where p is the component measured by the sensor, and v is the unmeasured component. Therefore, the measurement can be written as z = h ( p ) + ε , where ε is a Gaussian measurement noise with zero mean and covariance R . Particles p ( j ) can be obtained from N ( p ; h 1 ( z ) , H R H T ) , where H is the Jacobi matrix of h 1 , and particles v ( j ) from unmeasured subspace are set a priori. The weight of each particle is identical, i.e., w Γ , k 1 ( i , j ) = 1 L Γ , k 1 ( i ) . Consequently, the target birth density can be expressed as
p Γ , k 1 ( i ) ( x | z i ) = j = 1 L Γ , k 1 ( i ) w Γ , k 1 ( i , j ) δ x Γ , k 1 ( i , j ) ( x ) .
Target birth probability: The target birth probability r Γ , k 1 ( i ) of the newborn target intensity { ( r Γ , k 1 ( i ) , p Γ , k 1 ( i ) ( | z i ) ) } i = 1 | Z k 1 | is calculated by a pre-processing step. The predicted density of surviving targets, { ( r P , k | k 1 ( i ) , p P , k | k 1 ( i ) ) } i = 1 M k 1 , can be computed as follows:
r P , k | k 1 ( i ) = r k 1 ( i ) j = 1 L k 1 ( i ) w k 1 ( i , j ) p S , k ( x k 1 ( i , j ) ) ,
p P , k | k 1 ( i ) ( x ) = j = 1 L k 1 ( i ) w P , k | k 1 ( i , j ) δ x P , k | k 1 ( i , j ) ( x ) ,
where
x P , k | k 1 ( i , j ) q k ( i ) ( | x k 1 ( i , j ) , Z k ) , j = 1 , , L k 1 ( i ) ,
w P , k | k 1 ( i , j ) = w k 1 ( i , j ) f k | k 1 ( x P , k | k 1 ( i , j ) | x k 1 ( i , j ) ) p S , k ( x k 1 ( i , j ) ) q k ( i ) ( x P , k | k 1 ( i , j ) | x k 1 ( i , j ) , Z k ) / j = 1 L k 1 ( i ) w k 1 ( i , j ) p S , k ( x k 1 ( i , j ) ) ,
and q k ( i ) ( | x k 1 ( i , j ) , Z k ) represents the importance density. To simplify the SMC implementation, it is obtained that q k ( i ) ( | x k 1 ( i , j ) , Z k ) = f k | k 1 ( | x k 1 ( i , j ) ) .
The modified target birth probability r ˜ Γ , k 1 ( i ) is calculated by Equation (19), where r ˜ L Γ , k 1 ( i ) is given by
r ˜ L Γ , k 1 ( i ) = r ^ Γ , k 1 ( i ) 1 j = 1 L Γ , k 1 ( i ) w Γ , k 1 ( i , j ) p D , k ( x Γ , k 1 ( i , j ) ) 1 r ^ Γ , k 1 ( i ) j = 1 L Γ , k 1 ( i ) w Γ , k 1 ( i , j ) p D , k ( x Γ , k 1 ( j ) ) ,
and r ˜ U Γ , k 1 ( z ) can be written as Equation (21), where
G k ( z ) = r ^ Γ , k 1 ( i ) ( 1 r ^ Γ , k 1 ( i ) ) ( 1 r ^ Γ , k 1 ( i ) ) j = 1 L Γ , k 1 ( i ) w Γ , k 1 ( i , j ) ψ k , z ( x Γ , k 1 ( i , j ) ) ( 1 r ^ Γ , k 1 ( i ) j = 1 L Γ , k 1 ( i ) w Γ , k 1 ( i , j ) p D , k ( x Γ , k 1 ( i , j ) ) ) 2 ,
S Γ k ( z ) = i = 1 | Z k 1 | r ^ Γ , k 1 ( i ) j = 1 L Γ , k 1 ( i ) w Γ , k 1 ( i , j ) ψ k , z ( x Γ , k 1 ( i , j ) ) 1 r ^ Γ , k 1 ( i ) j = 1 L Γ , k 1 ( i ) w Γ , k 1 ( i , j ) p D , k ( x Γ , k 1 ( i , j ) ) ,
S U k ( z ) = i = 1 M k 1 r ^ P , k | k 1 ( i ) j = 1 L k 1 ( i ) w P , k | k 1 ( i , j ) ψ k , z ( x P , k | k 1 ( i , j ) ) 1 r ^ P , k | , k 1 ( i ) j = 1 L k 1 ( i ) w P , k | , k 1 ( i , j ) p D , k ( x P , k | , k 1 ( i , j ) ) .
Then, the target birth probability r Γ , k 1 ( i ) can be obtained using Equation (25).
(2) Prediction
In the normal filtering step, the predicted multi-target density, integrating { ( r Γ , k 1 ( i ) , p Γ , k 1 ( i ) ( | z i ) ) } i = 1 | Z k 1 | and { ( r P , k | k 1 ( i ) , p P , k | k 1 ( i ) ) } i = 1 M k 1 , can be expressed by { ( r k | k 1 ( i ) , p k | k 1 ( i ) ) } i = 1 M k | k 1 , where
p k | k 1 ( i ) ( x ) = j = 1 L k | k 1 ( i ) w k | k 1 ( i , j ) δ x k | k 1 ( i , j ) ( x ) .
(3) Update
The updated posterior multi-target density can be expressed in the same form as Equation (9) and is computed as follows:
r L , k ( i ) = r k | k 1 ( i ) 1 α L , k ( i ) 1 r k | k 1 ( i ) α L , k ( i ) ,
p L , k ( i ) ( x ) = j = 1 L k | k 1 ( i ) w L , k ( i , j ) δ x k | k 1 ( i , j ) ( x ) ,
r U , k ( z ) = i = 1 M k | k 1 r k | k 1 ( i ) ( 1 r k | k 1 ( i ) ) α U , k ( i ) ( z ) ( 1 r k | k 1 ( i ) α L , k ( i ) ) 2 κ k ( z ) + i = 1 M k | k 1 r k | k 1 ( i ) α U , k ( i ) ( z ) 1 r k | k 1 ( i ) α L , k ( i ) ,
p U , k ( x ; z ) = i = 1 M k | k 1 j = 1 L k | k 1 ( i ) w U , k ( i , j ) ( z ) δ x k | k 1 ( i , j ) ( x ) ,
where
α L , k ( i ) = j = 1 L k | k 1 ( i ) w k | k 1 ( i , j ) p D , k ( x k | k 1 ( i , j ) ) ,
w L , k ( i , j ) = w k | k 1 ( i , j ) ( 1 p D , k ( x k | k 1 ( i , j ) ) ) / j = 1 L k | k 1 ( i ) w k | k 1 ( i , j ) ( 1 p D , k ( x k | k 1 ( i , j ) ) ) ,
α U , k ( i ) ( z ) = j = 1 L k | k 1 ( i ) w k | k 1 ( i , j ) ψ k , z ( x k | k 1 ( i , j ) ) ,
w U , k ( i , j ) ( z ) = w k | k 1 ( i , j ) r k | k 1 ( i ) 1 r k | k 1 ( i ) ψ k , z ( x k | k 1 ( i , j ) ) i = 1 M k | k 1 j = 1 L k | k 1 ( i ) w k | k 1 ( i , j ) r k | k 1 ( i ) 1 r k | k 1 ( i ) ψ k , z ( x k | k 1 ( i , j ) ) .
Remark 2.
For the SMC implementation, weight degeneration is unavoidable. To mitigate the effect of degeneration, a resampling method is applied to the tracks after the update step. The number of targets is estimated by the cardinality mean, and the states of targets are estimated by computing the means of the posterior probability densities with comparatively high existence probabilities in π k .

4. Fast Sequential Monte Carlo Cardinality Balanced Multi-Target Multi-Bernoulli Filter with Adaptive Target Birth Intensity

In this section, a threshold method used in the calculation of the measurement likelihoods is proposed and then applied to the SMC CBMeMBer filter presented in Section 3. Figure 2 shows a framework of the CBMeMBer filter integrating the threshold technique.

4.1. Improved Measurement Likelihood

In general, the proposed CBMeMBer filter with adaptive target birth intensity has a heavy computational burden compared to the standard CBMeMBer filter. This is because the target birth density is modeled by the measurement, and the target birth probability is modeled by a pre-processing step, which requires more calculation in both the birth modeling step and the subsequent update step. The SMC implementation of this filter is notably time-consuming, especially for the cases with a large number of particles to describe the target states or a large number of measurements.
To reduce the computational burden, an improved measurement likelihood is proposed. In MTT, the measurement likelihood is generally computed by
l k ( x | z ) g k ( z | x ) .
Then, a pre-determination can be used to evaluate whether this measurement likelihood is worth being preserved or not.
Proposition 2.
For each target state x , there exists a mapping on observation space based on the measurement function, i.e.,
z = h ( x ) .
If the mapping z is far from the measurement z , the corresponding target state is rarely to be matched with the measurement, and the likelihood between them will be a small value, i.e., such a target state contributes little to the filters. Hence, we can directly set this likelihood as zero to avoid extra calculations. Below, a threshold is introduced to measure the distance between the measurement and the mapping of the target. The improved measurement likelihood can be expressed as
l ˜ k ( x | z ) = { g k ( z | x ) , i f   D ( z , z ) U 0 , o t h e r w i s e ,
where U represents the threshold and D ( z , z ) represents the distance between z  and z . Note that the mapping z = [ z 1 , , z n ] T = [ h ( x ) 1 , , h ( x ) n ] T and the measurement z = [ z 1 , , z n ] T are both multi-dimensional, where n is the dimension of the measurement space. Thus, the distance D ( z , z ) and the threshold U are both vectors, and the pre-determination condition in Equation (48) can be fully formulated as
D ( z , z ) U { d ( z 1 , z 1 ) u 1 d ( z n , z n ) u n ,
where d ( z i , z i ) ,   i = 1 , , n represents the distance, such as Euclidean distance, between the corresponding elements of the mapping z and the measurement z , and u i ,   i = 1 , , n represents the corresponding threshold.

4.2. Threshold Selection

In the improved measurement likelihood, the unnecessary calculation is recognized and removed through a threshold U . A suitable threshold should be a bridge which can properly reflect the correlation between the likelihood g k ( z | x ) and the distance D ( z , z ) . Here, we take the measurement noise as the basic part of the threshold selection.
Proposition 3.
Suppose that the measurement noise follows a zero-mean Gaussian distribution with covariance matrix
R = d i a g ( [ σ 1 2 , , σ n 2 ] T ) ,
where σ i ,   i = 1 , , n represents the standard deviation. The threshold U = [ u 1 , , u n ] can be set as [ β 1 × σ 1 , , β n × σ n ] , where β i ,   i = 1 , , n represents a non-negative scale parameter. The scale parameter can be seen as a coefficient of the standard deviation σ i and can be used to zoom the threshold element u 1 . For simplification, we set the scale parameters to be identical, i.e., β 1 = = β n = β . Thus, the threshold element can be replaced by
u i = β × σ i ,   i = 1 , , n .
Obviously, under known measurement function and noise, the probability of reserving the useful information can be easily worked out based on the rule of thumb, which can be expressed by
θ β = P ( d ( z i , z i ) β × σ i ) ,   i = 1 , , n ,
where P ( d ( z i , z i ) β × σ i ) represents the probability that the distance getting a value does not exceed the threshold in the i th dimension of the observation space. Through the cumulative distribution function (CDF) of the measurement noise, the probability θ β can be easily obtained. For instance, by assuming zero-mean Gaussian measurement noise, we obtain θ β = 1 68.27 % , θ β = 2 95.45 % , θ β = 3 99.73 % , and θ β = 4 99.99 % . On the other hand, if there exists an expectation of the probability θ β , e.g., θ β = 98 % , the appropriate scale parameter β can be easily calculated by applying the inverse operation of the CDF.
Proposition 4.
To correctly integrate the proposed measurement likelihood into the Bayesian based MTT systems, the clutter density κ k ( z ) used in updating should be adjusted. For the common uniform distribution clutter model, the adjusted clutter density can be expressed by
κ ˜ k ( z ) = θ β × κ k ( z ) .
It shows that only the main factor from the clutter is considered, with a proportion θ β that is consistent with the improved measurement likelihood. In the CBMeMBer filter, the update step using the new measurement of likelihood can be rewritten as
r U , k ( z ) = i = 1 M k | k 1 r k | k 1 ( i ) ( 1 r k | k 1 ( i ) ) p k | k 1 ( i ) , ψ ˜ k , z ( 1 r k | k 1 ( i ) p k | k 1 ( i ) , p D , k ) 2 κ ˜ k ( z ) + i = 1 M k | k 1 r k | k 1 ( i ) p k | k 1 ( i ) , ψ ˜ k , z 1 r k | k 1 ( i ) p k | k 1 ( i ) , p D , k ,
p U , k ( x ; z ) = i = 1 M k | k 1 r k | k 1 ( i ) 1 r k | k 1 ( i ) p k | k 1 ( i ) ( x ) ψ ˜ k , z ( x ) i = 1 M k | k 1 r k | k 1 ( i ) 1 r k | k 1 ( i ) p k | k 1 ( i ) , ψ ˜ k , z ,
where
ψ ˜ k , z ( x ) = l ˜ k ( x | z ) p D , k ( x ) .
Remark 3.
The target states updated in this step evolve from a prediction step, which is biased to the true target states. Thus, the proportion of the reserved contributions from targets cannot be exactly measured by θ β . Under general tracking scenarios, considering that the gap between the predicted and the true target states is seldom extremely large, the probability θ β can still be a reference in choosing the scale parameter β .

4.3. Fast Sequential Monte Carlo Adaptive Target Birth Intensity Cardinality Balanced Multi-Target Multi-Bernoulli Filter

In this subsection, the threshold method is applied to our proposed SMC CBMeMBer filter to relieve the huge computational burden by reducing the likelihood calculation. The new threshold SMC CBMeMBer filter is called the fast adaptive target birth intensity SMC CBMeMBer filter.
Proposition 5.
In the fast filter, the threshold is applied in both the pre-processing step and the update step. By substituting the clutter density κ k ( z ) with κ ˜ k ( z ) and the measurement likelihood l k ( x | z ) with l ˜ k ( x | z ) , Equation (21) in the pre-processing step can be rewritten as
r ˜ U Γ , k 1 ( i ) = z Z k G ˜ k ( z ) κ ˜ k ( z ) + S ˜ Γ k ( z ) + S ˜ U k ( z ) ,
where
G ˜ k ( z ) = r ^ Γ , k 1 ( i ) ( 1 r ^ Γ , k 1 ( i ) ) j = 1 L Γ , k 1 ( i ) w Γ , k 1 ( i , j ) ψ ˜ k , z ( x Γ , k 1 ( i , j ) ) ( 1 r ^ Γ , k 1 ( i ) j = 1 L Γ , k 1 ( i ) w Γ , k 1 ( i , j ) p D , k ( x Γ , k 1 ( i , j ) ) ) 2 ,
S ˜ Γ k ( z ) = i = 1 | Z k 1 | r ^ Γ , k 1 ( i ) j = 1 L Γ , k 1 ( i ) w Γ , k 1 ( i , j ) ψ ˜ k , z ( x Γ , k 1 ( i , j ) ) 1 r ^ Γ , k 1 ( i ) j = 1 L Γ , k 1 ( i ) w Γ , k 1 ( i , j ) p D , k ( x Γ , k 1 ( i , j ) ) ,
S ˜ U k ( z ) = i = 1 M k 1 r P , k | k 1 ( i ) j = 1 L k 1 ( i ) w P , k | k 1 ( i , j ) ψ ˜ k , z ( x P , k | k 1 ( i , j ) ) 1 r P , k | k 1 ( i ) j = 1 L k 1 ( i ) w P , k | k 1 ( i , j ) p D , k ( x P , k | k 1 ( i , j ) ) .
Equation (40) in the update step can be rewritten as
r U , k ( z ) = i = 1 M k | k 1 r k | k 1 ( i ) ( 1 r k | k 1 ( i ) ) α ˜ U , k ( i ) ( z ) ( 1 r k | k 1 ( i ) α L , k ( i ) ) 2 κ ˜ k ( z ) + i = 1 M k | k 1 r k | k 1 ( i ) α ˜ U , k ( i ) ( z ) 1 r k | k 1 ( i ) α L , k ( i ) ,
where
α ˜ U , k ( i ) ( z ) = j = 1 L k | k 1 ( i ) w k | k 1 ( i , j ) ψ ˜ k , z ( x k | k 1 ( i , j ) ) .
Equation (45) can be rewritten as
w U , k ( i , j ) ( z ) = w k | k 1 ( i , j ) r k | k 1 ( i ) 1 r k | k 1 ( i ) ψ ˜ k , z ( x k | k 1 ( i , j ) ) i = 1 M k | k 1 j = 1 L k | k 1 ( i ) w k | k 1 ( i , j ) r k | k 1 ( i ) 1 r k | k 1 ( i ) ψ ˜ k , z ( x k | k 1 ( i , j ) ) .
Remark 4.
Suppose that the size of the observation region in each dimension is represented by Λ i , i = 1 , , n ; since the clutter follows a uniform distribution, the measurement likelihood calculation from the clutter can be cut down by
( 1 i = 1 n 2 u i Λ i ) × 100 % .
The percentage can be very high since the size of the observation region is much larger than that of the threshold based on the measurement noise, i.e., Λ i 2 u i ,   i = 1 , , n .
Remark 5.
Note that the fast CBMeMBer filter updates the predicted density only considering the nearby measurements and, thus, avoids a large number of unnecessary or even harmful calculations from the long-distance clutters. It can greatly reduce the negative impact of the clutters to the filter and may improve the target tracking accuracy. It is worth noting that there exists a trade-off between letting the true measurement in and keeping the clutters out. In general, the threshold will be proper, i.e., large enough to involve the true measurement, but small compared to the whole observation region, under the scale parameters being set as β = 3 5 .

5. Simulations

In this section, the optimal sub-pattern assignment (OSPA) miss-distance [35] is chosen as the tracking performance metric, which captures the cardinality error, as well as the distance of individual elements, of two finite sets. The OSPA distance of two finite sets X = { x 1 , , x m } and Y = { y 1 , , y n } is given by
d ¯ p ( c ) ( X , Y ) : = { ( 1 n ( min π Π n i = 1 m d ( c ) ( x i , y π ( i ) ) p + c p ( n m ) ) ) 1 p , m n d ¯ p ( c ) ( Y , X ) , m > n ,
where Π k is the set of permutations on { 1 , 2 , , k } , k > 0 , and the distance d ( c ) ( x , y ) : = min ( c , x y ) , while p 1 is the order of the distance and c is the cut-off parameter.
To achieve reliable results, simulation results were acquired by averaging the outcomes of 200 Monte Carlo trials.

5.1. Availability of the Threshold Method

This subsection describes the design an experiment to demonstrate the effectiveness of the proposed threshold method. In the simulation, the observation region was a circle with a diameter of 4000 m, while the kinematic and the measurement models were both non-linear. Fourteen targets appeared successively in the scenario, where 12 targets were newborn in an unknown position and moment. Figure 3 shows the true trajectories of all targets in the scenario.
The kinematic model employed a constant turn model, where the position and turn rate process noises followed zero-mean Gaussian distributions with standard deviations of σ w = 20   m / s 2 and σ u = π / 180   rad / s , respectively. The measurement model exported range and bearing values, where the range and bearing measurement noises also followed zero-mean Gaussian distributions with standard deviations of σ r = 20   m and σ θ = π / 180   rad , respectively. The clutter followed a uniform distribution, and the number of the clutter followed a Poisson distribution with a mean value of 10. The survival and detection probabilities were p S , k ( x ) = 0.95 and p D , k ( x ) = 0.95 , respectively. Each Bernoulli component was expressed by 1000 particles.
The fast adaptive target birth intensity filter proposed in Section 4 was compared with the adaptive target birth intensity filter proposed in Section 3, by applying different thresholds with scale parameters chosen from β = 1 5 . The clutter density κ ˜ used in the fast filters could be determined from Equation (53) and were expressed as 0.68268 κ , 0.95449 κ , 0.99730 κ , 0.99993 κ , and 0.99999 κ , respectively.
Figure 4 shows the estimated number of targets for the proposed filter and fast filters. It can be observed that the fast filter with the threshold of β = 1 had an obvious underestimation in target number. This was because the scale parameter β = 1 was too narrow to make the threshold completely contain the true measurement, which eliminated large amounts of useful likelihood calculations. The performance of the fast filter with β = 2 outperformed the filters with other scale parameters, since the threshold of β = 2 struck the best balance between containing the true measurement and blocking the false associated measurements. The fast filter with β = 5 achieved a target number estimation most close to the filter with traditional likelihood, since the increasing threshold contained more falsely associated likelihood information. As the scale parameter increased, the fast filter asymptotically approached the filter with traditional likelihood.
Figure 5 compares the OSPA distances between the proposed filter and fast filters. Similar to the discussions in Figure 4, the OSPA distance of the fast filter with β = 1 was the largest, and best performance could be achieved when β = 2 . Table 2 shows the total OSPA distance improvement of the fast filters with β = 1 5 compared to the filter with traditional likelihood. Obviously, with the scale parameter increasing, the OSPA distance of the fast filter first improved from worst to optimal, and then gradually approached that of the filter with traditional likelihood. This is consistent with the analysis made in Section 4.3 and is consistent with the trend in the target number estimation.
Figure 6 shows the operation time of different filters when the clutter rate varied from 1 to 50. We can observe that the threshold method greatly reduced the time consumption of the proposed adaptive target birth intensity filter, especially in high clutter rates. When the clutter rate was 50, the time consumption of the proposed filter was reduced by about 400 % . In the fast filters, the computing time increased slowly with the increase of the scale parameter but still remained at the same level. This is because all of the thresholds with given scale parameters were very small compared with the whole observation region, as shown in Equation (64). If the threshold was greatly increased to an extremely large number such as in β = 100 , the computing time was similar or equal to the filter with traditional likelihood.
The simulations demonstrate that the proposed fast filter can significantly reduce the operation time, which restricts the practical application of the proposed filter with traditional likelihood and can slightly improve the tracking accuracy. In the given scenario, the threshold with β = 2 was the best, while other thresholds with larger scale parameters were also suitable for the filter and could help achieve similar tracking performance. Of course, too large thresholds were not included, which resulted in the fast filter having no difference with the filter adopting traditional measurement likelihood.

5.2. Fast SMC Adaptive Target Birth Intensity CBMeMBer Filter

In this subsection, we illustrate the advantages of the proposed fast CBMeMBer filter in adapting unknown newborn targets. Under the same scenario as in Section 5.1, the fast adaptive target birth intensity filter was compared to the filters with the known, fixed, adaptive target birth intensities and the adaptive target birth density. The scale parameter in the fast filter was chosen as β = 2. The modeling of different birth models is shown in Table 3, where detailed descriptions on the fixed birth intensity and the adaptive birth density can be found in References [6,13], respectively.
Figure 7 shows the estimated number of targets versus time for the fast filter and the filters with different target birth models. The estimated number of the fixed birth intensity filter was slightly larger than the known birth intensity filter due to the fact that the known birth intensity, requiring all information about newborn targets, was optimal. The three filters using adaptive birth models had one time-step delay when capturing newborn targets, since the birth densities in these filters were generated by previous measurements. It is obvious that our proposed adaptive birth intensity filter outperformed the adaptive birth density filter. This was because, in the proposed adaptive birth intensity, a pre-processing procedure was introduced to model the birth probability, which resulted in the expected number of newborn targets being closer to the truth, as shown in Figure 8. Both of these filters had larger target number estimation than the fixed birth intensity filter, because the measurement information contained in the birth components could be regarded as a supplement to the targets in updating. The fast adaptive birth intensity filter further improved the performance of the adaptive birth intensity filter in target number estimation, due to the reasons explained in Section 5.1. Excluding time steps where newborn targets appeared, the fast adaptive birth intensity filter achieved a nearly equivalent target number estimation as the fixed birth intensity filter.
Figure 9 compares the OSPA distances of the filters under consideration. In time steps where targets appeared, the three filters with adaptive birth models were much worse than the known and fixed birth intensity filters in OSPA distance due to the delay in capturing newborn targets. Similar to Figure 7, the proposed adaptive birth intensity filters outperformed the adaptive birth density filter, and the known birth intensity filter outperformed the fixed birth intensity filter. Excluding time steps with newborn targets appearing, the filters with adaptive birth models had lower OSPA distances than the known and fixed birth intensity filters, which was opposite to the comparison of target number estimation. This was because the adaptive birth models supplied extra measurement information which led to a gain in target number estimation but a correction in target state estimation. In these three filters, the fast adaptive birth intensity filter achieved the best performance in OSPA distance.
Figure 10 shows the operation time of the fast filter and the filters with different birth models. The fixed birth intensity filter ran slower than the known birth intensity filter, since it included more birth components which appeared in each filtering step. The adaptive birth intensity filter ran slower than the adaptive birth density filter, since it introduced a new pre-processing step to model the birth probability. Obviously, the two filters with the adaptive birth intensity and density had very huge time consumptions, especially under dense clutter backgrounds. This was because the adaptive birth models applying measurements to capture newborn targets were very intricate in updating. Table 4 shows the time consumptions of different filters with a clutter rate of 50. By employing the threshold, the fast adaptive birth intensity filter greatly reduced the operation time; thus, it achieved a good real-time property which was even slightly better than the known and fixed birth intensity filters.
The simulations demonstrated that the fast adaptive target birth intensity CBMeMBer filter can improve the tracking performance of newborn targets without prior birth information, and it can operate well with a good processing efficiency.

5.3. Fast SMC Adaptive Target Birth Intensity CBMeMBer Filter under Specific Conditions

In this subsection, simulations are presented to further demonstrate the validity and the practical application scope of the fast CBMeMBer filter. Firstly, the the fast filter and the filters with different target birth models were compared in another scenario which included intensive target maneuvering, e.g., thrusting, accelerating, and sharp turning. The high-maneuver scenario was generated by integrating the six targets of Reference [36]. Figure 11 shows the true target trajectories in this scenario.
To adapt the maneuvering of targets, multiple target motion models including constant velocity (CV), constant acceleration (CA), and constant turn (CT) models were applied. In practice, unexpected maneuvering, i.e., none of the models matching with the motion under consideration, is often encountered. To test the filters in this case, we supposed that the sharp turning, i.e., the second turn of target 4 from time 86 to 100, and the third and fourth turns of target 6 from time 116 to 125, and from time 156 to 160, was unexpected.
Figure 12 shows the tracking performance of the filters in the case of unexpected sharp target turning. It can be observed that the tracking performance of the filters declined due to the maneuvering of targets, especially when unexpected maneuvering happened. The obvious decline in tracking performance of the filters with the fixed and known birth intensities happened at times 86, 116, and 125 until the falsely tracked targets disappeared at time 188. Different from these two filters, the tracking performance of the filters with adaptive birth models improved after the disappearance of the unexpected maneuvering, and the fast filter outperformed its counterparts. This was because the measurement information in the adaptive birth components of the fast filter could be seen as a best initiation when the motion models were matched with the maneuvering again.
After this, another simulation was performed. The scenario adopted was the same as that in Figure 3 in Section 5.1 and Section 5.2, except that this experiment included continuous missing of sensor detections, which is always encountered in real-world tracking scenarios, such as targets being sheltered, sensor interference, active stealth, etc. In this simulation, we set a surviving target from time 36 to 110 with three short moments missing detection from time 41 to 43, and another surviving target from time 46 to 110 with 10 long moments missing detection from time 51 to 60.
Figure 13 shows the tracking performance of the filters in the case that continuous missing detection occurred. We can observe that the tracking performance of the filters with the fixed and known birth intensities declined after time 40 and further declined after time 50, although the measurements were obtained again at time 44 and 61, respectively. These two filters missed the targets thoroughly once the continuous missing detection occurred. Conversely, the tracking performance of the filters with the adaptive birth models only declined at time 41 to 44 and time 51 to 61, where the decline at time 44 and 61 was the one-step delay explained in Subsection 5.2. At time 45 to 50 and after time 61, the tracking performance quickly recovered when sensor detections were obtained again, where the fast filter performed best among them.
Under these specific conditions whereby unexpected target maneuvering happened and sensor detections were continuously missed, which is common in the practical applications, the proposed fast filter with the adaptive birth intensity was valid and outperformed other compared filters. Considering the real-time ability demonstrated in Subsection 5.2, the proposed fast filter is a highly adaptable filtering method.

6. Conclusions

In this paper, an RFS-based fast adaptive target birth intensity Bayes filter was applied under a multi-Bernoulli framework. By modeling the target birth density using previous measurements [13] and calculating the target birth probability in a pre-processing step, the tracking performance could be improved without prior information requirements. By applying an improved threshold-based measurement likelihood, the computational burden introduced by the adaptive target birth intensity could be greatly reduced and the tracking performance could be further improved. The proposed fast adaptive target birth intensity CBMeMBer filter can be applied in most tracking applications and can be implemented as efficiently as the standard CBMeMBer filter.

Author Contributions

Conceptualization, X.H.; methodology, X.H.; software, X.H. and L.L.; writing—original draft preparation, X.H.; writing—review and editing, H.J. and L.L.; project administration, H.J.

Funding

This work was funded by the National Natural Science Foundation of China (61871301, 61803288), and the China Postdoctoral Science Foundation (2018M633467).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mahler, R. Statistical Multisource-Multi-Target Information Fusion; Artech House: Norwood, MA, USA, 2007. [Google Scholar]
  2. Mahler, R. Multi-target Bayes filtering via first-order multi-target moments. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 1152–1178. [Google Scholar] [CrossRef]
  3. Vo, B.N.; Ma, W.K. The Gaussian mixture probability hypothesis density filter. IEEE Trans. Signal Process. 2006, 54, 4091–4104. [Google Scholar] [CrossRef]
  4. Mahler, R. PHD filters of higher order in target number. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 1523–1543. [Google Scholar] [CrossRef]
  5. Vo, B.T.; Vo, B.N.; Cantoni, A. Analytic implementations of the cardinalized probability hypothesis density filter. IEEE Trans. Signal Process. 2007, 55, 3553–3567. [Google Scholar] [CrossRef]
  6. Vo, B.T.; Vo, B.N.; Cantoni, A. The cardinality balanced multi-target multi-Bernoulli filter and its implementations. IEEE Trans. Signal Process. 2009, 57, 409–423. [Google Scholar]
  7. Vo, B.N.; Vo, B.T.; Pham, N.T.; Suter, D. Joint detection and estimation of multiple objects from image observations. IEEE Trans. Signal Process. 2010, 58, 5129–5241. [Google Scholar] [CrossRef]
  8. Vo, B.T.; Vo, B.N. Labeled Random Finite Sets and multi-object conjugate priors. IEEE Trans. Signal Process. 2013, 61, 3460–3475. [Google Scholar] [CrossRef]
  9. Reuter, S.; Vo, B.T.; Vo, B.N.; Dietmayer, K. The labelled multi-Bernoulli filter. IEEE Trans. Signal Process. 2014, 62, 3246–3260. [Google Scholar]
  10. Vo, B.N.; Vo, B.T.; Hoang, H.G. An efficient implementation of the generalized labeled multi-Bernoulli filter. IEEE Trans. Signal Process. 2017, 65, 1975–1987. [Google Scholar] [CrossRef]
  11. Ristic, B.; Clark, D.; Vo, B.N.; Vo, B.T. Adaptive target birth intensity for PHD and CPHD filters. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 1656–1668. [Google Scholar] [CrossRef]
  12. Lin, S.; Vo, B.T.; Nordholm, S.E. Measurement driven birth model for the generalized labeled multi-Bernoulli filter. In Proceedings of the 5th International Conference on Control, Automation and Information Sciences, Ansan, South Korea, 27–29 October 2016; pp. 94–99. [Google Scholar]
  13. Reuter, S.; Meissner, D.; Wilking, B.; Dietmayer, K. Cardinality balanced multi-target multi-Bernoulli filtering using adaptive birth distributions. In Proceedings of the 16th International Conference on Information Fusion, Istanbul, Turkey, 9–12 July 2013; pp. 1608–1615. [Google Scholar]
  14. Choi, M.E.; Seo, S.W. Robust multi-target tracking scheme based on Gaussian mixture probability hypothesis density filter. IEEE Trans. Vehicl. Technol. 2016, 65, 4217–4229. [Google Scholar] [CrossRef]
  15. Choi, B.; Park, S.; Kim, E. A newborn track detection and state estimation algorithm using Bernoulli random finite sets. IEEE Trans. Signal Process. 2016, 64, 2660–2674. [Google Scholar] [CrossRef]
  16. Dong, P.; Jing, Z.; Gong, D.; Tang, B. Maneuvering multi-target tracking based on variable structure multiple model GMCPHD filter. Signal Process. 2017, 141, 158–167. [Google Scholar] [CrossRef]
  17. Liu, H.; Wu, W. Interacting multiple model (IMM) fifth-degree spherical simplex-radial cubature Kalman filter for maneuvering target tracking. Sensors 2017, 17, 1374. [Google Scholar] [CrossRef] [PubMed]
  18. Williams, J.L. An efficient, variational approximation of the best fitting multi-Bernoulli filter. IEEE Trans. Signal Process. 2015, 63, 258–273. [Google Scholar] [CrossRef]
  19. Williams, J.L. Hybrid Poisson and multi-Bernoulli filters. In Proceedings of the 15th International Conference on Information Fusion, Singapore, 9–12 July 2012; pp. 1103–1110. [Google Scholar]
  20. Jiang, H.; Cai, Y. Adaptive fifth-degree cubature information filter for multi-sensor bearings-only tracking. Sensors 2018, 18, 3241. [Google Scholar] [CrossRef] [PubMed]
  21. Hao, Y.; Xu, A.; Sui, X.; Wang, Y. A modified extended Kalman filter for a two-antenna GPS/INS vehicular navigation system. Sensors 2018, 18, 3809. [Google Scholar] [CrossRef] [PubMed]
  22. Dey, A.; Das, M.; Sadhu, S.; Ghoshal, T.K. Adaptive divided difference filter for parameter and state estimation of non-linear systems. IET Signal Process. 2015, 9, 369–376. [Google Scholar] [CrossRef]
  23. Liu, Z.; Chen, S.; Wu, H.; He, R.; Hao, L. A Student’s t mixture probability hypothesis density filter for multi-target tracking with outliers. Sensors 2018, 18, 1095. [Google Scholar] [CrossRef] [PubMed]
  24. Li, S.; Yi, W.; Hoseinnezhad, R.; Wang, B.; Kong, L. Multiobject tracking for generic observation model using labeled random finite sets. IEEE Trans. Signal Process. 2018, 66, 368–383. [Google Scholar] [CrossRef]
  25. Imani, M.; Braga-Neto, U.M. Maximum-likelihood adaptive filter for partially observed boolean dynamical systems. IEEE Trans. Signal Process. 2017, 65, 359–371. [Google Scholar] [CrossRef]
  26. He, S.; Shin, H.S.; Tsourdos, A. Joint probabilistic data association filter with unknown detection probability and clutter rate. Sensors 2018, 18, 269. [Google Scholar]
  27. Papa, G.; Braca, P.; Horn, S.; Marano, S.; Matta, V.; Willett, P. Multisensor adaptive Bayesian tracking under time-varying target detection probability. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 2193–2209. [Google Scholar] [CrossRef]
  28. Vo, B.T.; Vo, B.N.; Cantoni, A. Bayesian filtering with random finite set observations. IEEE Trans. Signal Process. 2008, 56, 1313–1326. [Google Scholar] [CrossRef]
  29. Blackman, S.S.; Popoli, R. Design and Analysis of Modern Tracking Systems; Artech House: Norwood, MA, USA, 1999. [Google Scholar]
  30. Zhang, H.; Jing, Z.; Hu, S. Gaussian mixture CPHD filter with gating technique. Signal Process. 2009, 89, 1521–1530. [Google Scholar] [CrossRef]
  31. Yoon, J.H.; Kim, D.Y.; Yoon, K.J. Efficient importance sampling function design for sequential Monte Carlo PHD filter. Signal Process. 2012, 92, 2315–2321. [Google Scholar] [CrossRef]
  32. Kim, D.Y.; Liang, M.; Jeon, M. Multi-Bernoulli filter based sensor selection with limited sensing range for multi-target tracking. In Proceedings of the 4th International Conference on Control, Automation and Information Sciences, Changshu, China, 29–31 October 2016; pp. 74–79. [Google Scholar]
  33. Gunes, A.; Guldogan, M.B. Multi-target bearing tracking with a single acoustic vector sensor based on multi-Bernoulli filter. In Proceedings of the 2015 Oceans, Genoa, Italy, 18–21 May 2015; pp. 1–5. [Google Scholar]
  34. Hoak, A.; Medeiros, H.; Povinelli, R.J. Image-based multi-target tracking through multi-Bernoulli filtering with interactive likelihoods. Sensors 2017, 17, 501. [Google Scholar] [CrossRef] [PubMed]
  35. Schuhmacher, D.; Vo, B.T.; Vo, B.N. A consistent metric for performance evaluation of multi-object filters. IEEE Trans. Signal Process. 2008, 56, 3447–3457. [Google Scholar] [CrossRef]
  36. Blair, W.D.; Watson, G.A. Benchmark Problem for Radar Resource Allocation and Tracking Maneuvering Targets in the Presence of ECM; Technical Report, Document, No. NSWCDD/TR-96/10; Naval Surface Warfare Center Dahlgren Division: Dahlgren, VA, USA, 1996. [Google Scholar]
Figure 1. Modeling of the adaptive target birth intensity.
Figure 1. Modeling of the adaptive target birth intensity.
Sensors 19 01120 g001
Figure 2. Cardinality balanced multi-target multi-Bernoulli (CBMeMBer) filter with threshold technique.
Figure 2. Cardinality balanced multi-target multi-Bernoulli (CBMeMBer) filter with threshold technique.
Sensors 19 01120 g002
Figure 3. Target trajectories in polar coordinates with start and stop positions denoted by △ and ○.
Figure 3. Target trajectories in polar coordinates with start and stop positions denoted by △ and ○.
Sensors 19 01120 g003
Figure 4. Target number estimation of the sequential Monte Carlo (SMC) adaptive target birth intensity CBMeMBer filter and the fast SMC adaptive target birth intensity CBMeMBer filter.
Figure 4. Target number estimation of the sequential Monte Carlo (SMC) adaptive target birth intensity CBMeMBer filter and the fast SMC adaptive target birth intensity CBMeMBer filter.
Sensors 19 01120 g004
Figure 5. Optimal sub-pattern assignment (OSPA) distances of the SMC adaptive target birth intensity CBMeMBer filter and the fast SMC adaptive target birth intensity CBMeMBer filter.
Figure 5. Optimal sub-pattern assignment (OSPA) distances of the SMC adaptive target birth intensity CBMeMBer filter and the fast SMC adaptive target birth intensity CBMeMBer filter.
Sensors 19 01120 g005
Figure 6. Time consumptions of the SMC adaptive target birth intensity CBMeMBer filter and the fast SMC adaptive target birth intensity CBMeMBer filter with clutter rate varied from 1 to 50.
Figure 6. Time consumptions of the SMC adaptive target birth intensity CBMeMBer filter and the fast SMC adaptive target birth intensity CBMeMBer filter with clutter rate varied from 1 to 50.
Sensors 19 01120 g006
Figure 7. Target number estimation of the fast adaptive target birth intensity CBMeMBer filter and other CBMeMBer filters with different target birth models.
Figure 7. Target number estimation of the fast adaptive target birth intensity CBMeMBer filter and other CBMeMBer filters with different target birth models.
Sensors 19 01120 g007
Figure 8. Expected number of newborn targets in constant and adaptive target birth probabilities.
Figure 8. Expected number of newborn targets in constant and adaptive target birth probabilities.
Sensors 19 01120 g008
Figure 9. OSPA distances of the fast adaptive target birth intensity CBMeMBer filter and other CBMeMBer filters with different target birth models.
Figure 9. OSPA distances of the fast adaptive target birth intensity CBMeMBer filter and other CBMeMBer filters with different target birth models.
Sensors 19 01120 g009
Figure 10. Time consumptions of the fast adaptive target birth intensity CBMeMBer filter and other CBMeMBer filters with different target birth models with clutter rates ranging from 1 to 50.
Figure 10. Time consumptions of the fast adaptive target birth intensity CBMeMBer filter and other CBMeMBer filters with different target birth models with clutter rates ranging from 1 to 50.
Sensors 19 01120 g010
Figure 11. Trajectories in Cartesian coordinates with start and stop positions denoted by △ and ○.
Figure 11. Trajectories in Cartesian coordinates with start and stop positions denoted by △ and ○.
Sensors 19 01120 g011
Figure 12. Tracking performance of the fast adaptive target birth intensity CBMeMBer filter and other CBMeMBer filters with different target birth models in the case of unexpected maneuvering: (a) target number estimation; (b) OSPA distances.
Figure 12. Tracking performance of the fast adaptive target birth intensity CBMeMBer filter and other CBMeMBer filters with different target birth models in the case of unexpected maneuvering: (a) target number estimation; (b) OSPA distances.
Sensors 19 01120 g012
Figure 13. Tracking performance of the fast adaptive target birth intensity CBMeMBer filter and other CBMeMBer filters with different target birth models in the case of continuously missing detections: (a) target number estimation; (b) OSPA distances.
Figure 13. Tracking performance of the fast adaptive target birth intensity CBMeMBer filter and other CBMeMBer filters with different target birth models in the case of continuously missing detections: (a) target number estimation; (b) OSPA distances.
Sensors 19 01120 g013
Table 1. Modeling of the target birth probability.
Table 1. Modeling of the target birth probability.
• Input: Z k 1 , B ^ Γ , k , π k 1 , Z k .
• Output: r Γ , k 1 ( i ) ,   i = 1 , , | Z k 1 | .
Step 1: p Γ , k 1 ( i ) ( | z i ) ,   i = 1 , , | Z k 1 | according to Z k 1 [13]
   r ^ Γ , k 1 ( i ) = B ^ Γ , k / | Z k 1 | ;
Step 2: π P , k | k 1 according to π k 1 , Equations (6) and (7);
Step 3: r ˜ L Γ , k 1 ( i ) according to p Γ , k 1 ( i ) ( | z i ) , r ^ Γ , k 1 ( i ) with Equation (20);
   r ^ U Γ , k 1 ( i ) according to p Γ , k 1 ( i ) ( | z i ) , r ^ Γ , k 1 ( i ) , π P , k | k 1 with Equation (21)
   r ^ Γ , k 1 ( i ) = r ^ L Γ , k 1 ( i ) + r ^ U Γ , k 1 ( i ) , i.e., Equation (19);
Step 4: r Γ , k 1 ( i ) = min ( r ^ Γ , k 1 ( i ) , r min ) ,   i = 1 , , | Z k 1 | , i.e., Equation (25).
Table 2. Total optimal sub-pattern assignment (OSPA) distance improvement of different fast filters.
Table 2. Total optimal sub-pattern assignment (OSPA) distance improvement of different fast filters.
Fast filters β = 1 β = 2 β = 3 β = 4 β = 5
Improvement 24.37 % 8.33 % 6.82 % 2.73 % 1.19 %
Table 3. Modeling of known, fixed, adaptive target birth intensities and adaptive target birth density.
Table 3. Modeling of known, fixed, adaptive target birth intensities and adaptive target birth density.
Target Birth ModelsBirth PositionBirth Probability
known birth intensitytrueTrue
fixed birth intensitysix possible appearing areasconstant, r Γ , k ( i ) = 0.03 ,   i = 1 6
adaptive birth densityprevious measurement areasconstant, B Γ , k = 0.18
adaptive birth intensityprevious measurement areasadaptively modified
Table 4. Operation time of different filters with a clutter rate of 50.
Table 4. Operation time of different filters with a clutter rate of 50.
Filtering MethodsFixed BirthKnown BirthAdaptive Birth DensityAdaptive BirthAdaptive Birth and Threshold
Opreation time73.8975 s59.7363 s208.6728 s251.4441 s52.7258 s

Share and Cite

MDPI and ACS Style

Hu, X.; Ji, H.; Liu, L. Adaptive Target Birth Intensity Multi-Bernoulli Filter with Noise-Based Threshold. Sensors 2019, 19, 1120. https://doi.org/10.3390/s19051120

AMA Style

Hu X, Ji H, Liu L. Adaptive Target Birth Intensity Multi-Bernoulli Filter with Noise-Based Threshold. Sensors. 2019; 19(5):1120. https://doi.org/10.3390/s19051120

Chicago/Turabian Style

Hu, Xiaolong, Hongbing Ji, and Long Liu. 2019. "Adaptive Target Birth Intensity Multi-Bernoulli Filter with Noise-Based Threshold" Sensors 19, no. 5: 1120. https://doi.org/10.3390/s19051120

APA Style

Hu, X., Ji, H., & Liu, L. (2019). Adaptive Target Birth Intensity Multi-Bernoulli Filter with Noise-Based Threshold. Sensors, 19(5), 1120. https://doi.org/10.3390/s19051120

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