Next Article in Journal
Detection of Pre-Impact Falls from Heights Using an Inertial Measurement Unit Sensor
Next Article in Special Issue
Implementing Deep Learning Techniques in 5G IoT Networks for 3D Indoor Positioning: DELTA (DeEp Learning-Based Co-operaTive Architecture)
Previous Article in Journal
A Real-Time Physical Progress Measurement Method for Schedule Performance Control Using Vision, an AR Marker and Machine Learning in a Ship Block Assembly Process
Previous Article in Special Issue
Leveraging Uncertainties in Softmax Decision-Making Models for Low-Power IoT Devices
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Extended Target Marginal Distribution Poisson Multi-Bernoulli Mixture Filter

Automatic Target Recognition (ATR) Key Laboratory, Shenzhen University, Shenzhen 518060, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(18), 5387; https://doi.org/10.3390/s20185387
Submission received: 3 August 2020 / Revised: 13 September 2020 / Accepted: 16 September 2020 / Published: 20 September 2020

Abstract

:
The existence of clutter, unknown measurement sources, unknown number of targets, and undetected probability are problems for multi-extended target tracking, to address these problems; this paper proposes a gamma-Gaussian-inverse Wishart (GGIW) implementation of a marginal distribution Poisson multi-Bernoulli mixture (MD-PMBM) filter. Unlike existing multiple extended target tracking filters, the GGIW-MD-PMBM filter computes the marginal distribution (MD) and the existence probability of each target, which can shorten the computing time while maintaining good tracking results. The simulation results confirm the validity and reliability of the GGIW-MD-PMBM filter.

1. Introduction

Multiple target tracking (MTT) involves estimating the state of an unknown number of targets in the presence of clutter [1,2,3]. The traditional MTT algorithm assumes each target is a point as one target generates at most one measurement in the sensor at each time step [3]. However, the rapid development of the modern high-resolution sensors makes the “point assumption” impractical because with such sensors, one target generates at least one measurement per time scan. The MTT problem with such sensors becomes an extended MTT problem [4,5]. Compared with the point target, the extended target can not only provide accurate target movement information, but also the target’ shape and size information, making it useful in the fields of robot recognition and positioning, moving crowd tracking, and tracking of close cars or large ships using the high-resolution sensors or automotive radars. Due to its wide application, the research on extended target tracking has received considerable attention from many scholars, making it a growing research hotspot.
As one of the most used extended target measurement models, The Poisson Point Process (PPP) model assumes a Poisson distributed random number of measurements are spatially distributed around the target at each time step [4,5]. To efficiently and correctly represent the spatially distributed measurements, many shape models have been developed, such as the random matrix model (RMM) [6,7,8,9,10,11,12,13,14,15], random hypersurface model (RHM) [16,17,18,19,20,21] and Gaussian process model (GPM) [22]. RMM assumes that the shape of the measurements can be approximated by an ellipse, and the measurements surrounding the centroid of the target obey a Gaussian distribution. RHM uses a more general star-convex shape to approximate the target contour. The GPM automatically learns the shape of the target through a Gaussian process and can give an analytic expression of its contour for an arbitrary shape target. Among the above three measurement models, the RMM has the smallest amount of calculation and is easy to implement, and can meet the requirements of many practical situations, so we use the RMM method in this paper.
Both in the point target-tracking filters and the extended target-tracking filters, the number of targets and the number of measurements are unknown and time-varying; thus, how to represent the targets and the measurements is a challenge. The random finite sets (RFS) theory [23], which represents the targets and the measurements as a finite variable set, is a suitable choice. To solve the point MTT problem, in the early stage, many RFS-based filters have been proposed, such as the Probability Hypothesis Density (PHD) filter [24,25,26], the Cardinalized Probability Hypothesis Density (CPHD) filter [27,28,29,30] and a series of multi-Bernoulli (MB) filters [31,32,33,34]. In recent years, scholars have proposed many RFS-based filters to solve the extended MTT problem, such as PHD for extended target tracking (ETT-PHD) [9,35,36,37], ETT-CPHD [38,39,40], gamma-Gaussian-inverse Wishart-Poisson multi-Bernoulli mixture (GGIW-PMBM) [40,41], GGIW implementation of the Labelled Multi-Bernoulli (GGIW-LMB) [42,43], and so on [22,44].
In the above-listed RFS-based filters, many kinds of conjugate priors are widely used, to give a closed-form solution of the posterior density. The prior conjugate refers to a time series of random finite set that satisfies the conjugate prior property; if the prediction is a Gaussian (Multi-Bernoulli) distribution, the update is also a Gaussian (Multi-Bernoulli) distribution. Using the conjugate prior, given enough parameters, the posterior distribution can be approximated arbitrarily.
In recent years, developing conjugate prior-based MTT filters has been a significant trend; the two kinds of most used conjugate priors in various MTT filters are the PMBM conjugate prior and the δ-Generalized Labelled Multi-Bernoulli (δ-GLMB) conjugate prior. In the δ-GLMB conjugate prior, the state of each target is added with a unique label, which is helpful to find each target’s trajectory. The PMBM prior conjugate divides the target set into two disjoint subsets, targets that have been detected and targets that have not yet been detected, and then processes the two subsets separately. This strategy can greatly reduce the target missed detection rate and improve tracking accuracy. In reference [40], it has been shown that the point target tracking filter based on PMBM conjugate prior outperforms the filter based on the GLMB conjugate prior in terms of tracking accuracy and computing time; thus, in this paper we used the PMBM conjugate prior.
In this paper, to solve the extended MTT problem and combine the PPP model and PMBM conjugate prior, we propose a GGIW implementation of a marginal distribution Poisson multi-Bernoulli mixture (MD-PMBM) filter on the basis of the GGIW-PMBM filter. The main difference between the GGIW-MD-PMBM filter and the GGIW-PMBM filter is that in the prediction step and the update step, instead of recursively propagating the joint state distribution of targets (GGIW-PMBM), the GGIW-MD-PMBM filter recursively propagates the marginal distribution and the existence probability of each target; such a strategy can effectively reduce the calculation cost while maintaining good tracking results.
The rest of the paper is organized as follows. Several model assumptions, such as the Bayesian model, PPP model, MBM model, motion model, measurement model, and PMBM conjugate prior, are discussed in Section 2. Section 3 presents the prediction step and the update step of the GGIW-MD-PMBM filter. Section 4 gives the four extended target filters’ experimental simulation results, which are summarized in Section 5.

2. Background

In this section, we outline various model assumptions, such as the Bayesian model [45,46,47], PPP model, Multi-Bernoulli (MB) process model, motion and measurement model, and PMBM conjugate prior model, which were used in the prediction and update steps of the GGIW-MD-PMBM filter.

2.1. Bayesian Model

We define the state of the i-th extended target at time step k as ξ k ( i ) , and the state of the target set at time step k can be written as
X k = { ξ k ( i ) } i = 1 N ξ , k ,
where N ξ , k =   | X k | is the number of the target, which is unknown and time-varying. The measurements (the union of target generated measurements and clutter) obtained from the sensor at time k are
Z k = { z k j } j = 1 N z , k
where N z , k =   | Z k | is the number of measurements. Let Z k denote the union of measurements set from time step 1 to k, that is, Z k =   Z 1 Z 2 Z k .
In many point target-tracking filters and extended target-tracking filters, the main purpose of using the Bayesian model is to recursively propagate the multi-target density distribution. In a Bayesian model, the multi-target set density, multi-target transition density, and the multi-target measurement likelihood are represented by f k | k ( X k | Z k ) ,   f k | k 1 ( X k | X k 1 ) , and f k | k ( Z k | X k ) , respectively. Using the Chapman–Kolmogorov equation, the multi-target set density can be defined as
f k | k 1 ( X k | Z k 1 ) = f k | k 1 ( X k X k 1 ) f k 1 k 1 ( X k 1 Z k 1 ) δ X k 1 ,
and the Bayes update as
f k k ( X k Z k ) = f k | k ( Z k X k ) f k k 1 ( X k Z k 1 ) f k | k ( Z k X k ) f k k 1 ( X k Z k 1 ) δ X k ,

2.2. PPP Model and Multi-Bernoulli (MB) Process Model

The PPP model, proposed by Gilholm et al. [4] and Granstrom et al. [5], is widely used in point target tracking filters and extended target tracking filters. In a PPP model, each target generates a Poisson distributed number of measurements, and each measurement is independent and identically distributed (i.i.d.). The Poisson rate μ and the spatial distribution f ( · ) determine the PPP intensity.
We define the spatial distribution of the target state ξ k ( i ) as f ( ξ k ( i ) ) . The Poisson density distribution of X k is
f ( X k ) = e μ ξ k ( i ) X k μ f ( ξ k ( i ) ) .
The existence of a single target can be modeled by a Bernoulli RFS [14,31,32,33,34,40,41,42]. The cardinality of a Bernoulli RFS X k ( i ) can either be 1 (with probability r k ( i ) ) or empty (with probability 1 r k ( i ) ), and the Bernoulli density distribution of X k ( i ) can be defined as
f ( X k ( i ) ) = { 1 r k ( i ) , X k ( i ) = ϕ r k ( i ) f ( ξ k ( i ) ) , X k ( i ) = { ξ k ( i ) } 0 , | X k ( i ) | 2 .
An MB RFS X k is the union of a limited number of independent Bernoulli RFS X k ( i ) , X k = i = 1 N ξ , k X k ( i ) . The MB density distribution of X k is
f ( X k ) = ( i = 1 N ξ , k ( 1 r k ( i ) ) ) 1 i 1 i n N ξ , k r k ( i 1 ) f ( ξ k ( i 1 ) ) 1 r k ( i 1 ) r k ( i n ) f ( ξ k ( i n ) ) 1 r k ( i n ) .
In this paper, we divide the targets into the undetected target set and the detected target set. The PPP represents the distribution of undetected target measurements, and the MB mixture (MBM) represents the distribution of detected target measurements.

2.3. The Motion Model and Measurement Model

To simplify the calculation, at each time step, we assume each target evolves independently from the other targets, and the measurements generated by each target are independent of each other and independent of those generated by other targets.
We assume each birth target can be described by a PPP model, and its intensity is D b ( x k ) . All the existing targets either survive with a probability P D ( x k ) or die with a probability 1 P D ( x k ) from time step k to time step k + 1 . Each extended target state consists of three variables, the measurement rate γ k , kinematic state x k , and extent state E k , that is, ξ k ( i ) = { γ k ( i ) , x k ( i ) , Ε k ( i ) } . Given γ k , x k and E k , the Markov motion model can be defined as
γ k + 1 = γ k ,
x k + 1 = f ( x k ) + w k ,
E k + 1 = M ( x k ) E k M ( x k ) T ,
where w k is Gaussian process noise with zero mean and covariance Q , f ( · ) and M ( · ) are the transformation matrix of kinematic state and extent state, respectively.
The clutter measurements at time step k can be modeled as a PPP model with measurement rate λ   and spatial Poisson distribution c ( z )   . The clutter PPP density function is k ( z ) = λ   c ( z )   . The measurements generated by each target at time step k are a Poisson distribution with a measurement rate γ k . The measurement model is
z k = H k x k + υ k ,
where H k is the known observation model function, and υ k is the Gaussian observation noise with zero mean and covariance E k .
Assuming the measurements generated by a single target follow a Gaussian distribution, that is
p ( z k | x k , υ k ) = N ( z k ; H k x k , E k ) .
The distribution of each target state can be expressed as
p ( γ k , x k , E k | Z k ) = p ( γ k | Z k ) p ( x k , E k | Z k ) p ( E k | Z k ) = G ( γ k ; α k | k , β k | k ) N ( x k ; m k | k , P k | k E k ) I W d ( E k ; v k | k , V k | k )   = G G I W ( ξ k ; ζ k | k )
where G ( γ k ; α k | k , β k | k ) is the Gamma distribution, α k | k > 0   and β k | k > 0 are its shape parameter and inverse scale parameter; N ( x k ; m k | k , P k | k E k ) is the Gaussian distribution, m k | k and P k | k E k are its mean and covariance; I W d ( E k ; v k | k , V k | k ) is the inverse Wishart distribution, v k | k and V k | k are its degrees of freedom and shape matrix. Note that ζ k | k = { α k | k , β k | k , m k | k , P k | k , v k | k , V k | k } is a shorthand for the set of GGIW density parameters. The detailed implementation of GGIW prediction and update is given in [21,23,24,25,26].

2.4. PMBM Conjugate Prior

Existing studies [40,41,42] show that developing PMBM conjugate priors-based filters is a significant trend, due to their excellent tracking performance. In a PMBM model, the extended target set X k at time k is divided into two disjoint sets, undetected targets X k u and detected targets X k d , then the PPP model is used to describe the distribution of undetected targets X k u , and the MBM model to describe the distribution of targets X k d that have been detected at least once.
X k u , X k d : X k = X k u X k d , X k u X k d = .
The PMBM set density can be defined as
f k | k ( X k | Z k ) = X k u   X k d = X k   f k | k u ( X k u | Z k ) j J ω k | k d , j f k | k d , j ( X k d | Z k ) ,
f k | k u ( X k u | Z k ) = e μ k | k u ξ X k u μ k | k u f k | k u ( ξ ) ,
f k | k d , j ( X k d | Z k ) = i I j X k i = X k d i I j f k | k d , j , i ( X k i ) ,
where f k | k d , j , i ( · ) is the Bernoulli set density. The MBM has   | J |   components, and J is the index set (or MBM component) of the MB in the MBM. I j is the Bernoulli index set of the j-th MB, that is, the j-th component has   | I j | Bernoulli components. ω d , j is the probability of the j-th MB component. The conjugacy property determines the intensities of the birth PPP and the initial undetected PPP, which are expressed by both GGIW mixture density.
Given the PMBM conjugate prior assumption, the extended target tracking is equivalent to recursively propagate the PMBM density parameters in the prediction step and the update step, which is shown in Section 3.

3. The GGIW-MD-PMBM Filter

In this section, the GGIW-MD-PMBM filter is presented. Instead of recursively propagating the joint state distribution of targets (GGIW-PMBM filter), the proposed GGIW-MD-PMBM filter recursively propagates the marginal distributions and the existence probabilities of each extended target. The GGIW-MD-PMBM filter can be modeled by the following assumptions:
(1) Clutter is uniformly distributed in the surveillance area with a given Poisson rate λ and is independent of targets’ distribution;
(2) Each extended target generates at least one measurement per time scan and evolves independently of the other targets;
(3) The birth PPP intensity is a GGIW mixture with a given Poisson rate μ b , weight ω b , and various density parameters ς b of GGIW;
D k + 1 b = μ k + 1 b j = 1 N k + 1 b ω k + 1 ( b , j ) G G I W ( ξ k + 1 ; ς k + 1 ( b , j ) ) .
(4) The initial undetected PPP intensity is also a GGIW mixture with a known Poisson rate μ 0 u , weight ω 0 u , and various density parameters ς 0 u of GGIW.
D 0 u = μ 0 u j = 1 N 0 u ω 0 ( u , j ) G G I W ( ξ 0 ; ς 0 ( u , j ) ) .
(5) The initial PMBM parameter is empty.
μ 0 | 0 u = 0   and   J 0 | 0 = 0 .
We assume the birth PPP and the initial undetected targets are all GGIW mixture intensities; according to the PMBM conjugate property, each target density in the PMBM filter is also a GGIW density.

3.1. Prediction

3.1.1. Detected Targets

According to the Chapman–Kolmogorov prediction (3) and the posterior PMBM density (15–16), the predicted distribution f k | k 1 d , j ( ξ k ) , j = 1 , , J d , k | k 1 and the existence probability r k | k 1 d , j of each detected target at time step k are
f k | k 1 d , j ( ξ k ) = G G I W ( ξ k ; ζ k | k 1 d , j ) , j = 1 , , J d , k | k 1 ,
r k | k 1 d , j = p S r k 1 | k 1 d , j , j = 1 , , J d , k | k 1 ,
where p S is the survival probability, and ζ k | k 1 d , j are the predicted parameters. The weights and number of predicted MBM components are ω d , k | k 1 j = ω d , k 1 | k 1 j   and J d , k | k 1 = J d , k 1 | k 1   , respectively.

3.1.2. Undetected Targets

To avoid missed detection, we regard previously undetected targets and current birth targets as the undetected target set, that is, the undetected targets at time step k include the undetected targets at time step k 1 and the birth targets at time step k (this paper does not consider the spawning targets). The predicted PPP for each undetected target at time step k is as follows:
f k u b , j ( ξ k )   = G G I W ( ξ k ; ζ k u b , j ) , j = 1 , , N k u b ,
f k | k 1 u u , j ( ξ k )   = G G I W ( ξ k ; ζ k | k 1 u u , j ) , j = 1 , , N k | k 1 u u .
The existence probability is ω k u b , j or ω k | k 1 u u , j .
ω k u b , j , j = 1 , , N k u b ,
ω k | k 1 u u , j = ω k 1 | k 1 u u , j p S , j = 1 , , N k | k 1 u u ,
where N k u b is the number of birth targets at time step k , and N k | k 1 u u is the number of undetected targets at time step k 1 . ζ k | k 1 u , j are the parameters ( ζ k u b , j or ζ k | k 1 u u , j ) of undetected targets at time step k . The predicted PPP for the undetected targets has a Poisson rate
μ k | k 1 u = μ k u b + P S u u μ k 1 | k 1 u u ,
where
P S u u = j = 1 N k 1 | k 1 u u ω k 1 | k 1 u u , j p S .
The existence probability of undetected targets ω k u b , j and ω k | k 1 u u , j can be normalized, as ω k u b , j = μ k u b μ k | k 1 u ω k u b , j , ω k | k 1 u u , j = μ k 1 | k 1 u u μ k | k 1 u ω k | k 1 u u , j , and j = 1 N k u b ω k u b , j +   j = 1 N k u u ω k | k 1 u u , j = 1 . For simplicity, the existence probability of each undetected target at time k can be represented as ω k | k 1 u , j , which includes ω k u b , j and ω k | k 1 u u , j .
Through the above analysis, in the prediction step, all extended target distributions and existence probabilities can be described as
f k | k 1 j ( ξ k )   , j = 1 , , N k | k 1 ,
ω k | k 1 j , j = 1 , , N k | k 1 ,
where N k | k 1 = J d , k | k 1 + N k | k 1 u u + N k u b . f k | k 1 j ( ξ k ) = f k | k 1 d , j ( ξ k ) and ω k | k 1 j = r k | k 1 d , j for 1 j J d , k | k 1 , f k | k 1 j J d , k | k 1 ( ξ k ) = f k | k 1 u u , j ( ξ k )   and ω k | k 1 j J d , k | k 1 = ω k | k 1 u u , j for J d , k | k 1 < j J d , k | k 1 + N k | k 1 u u , f k | k 1 j J d , k | k 1 N k u u ( ξ k ) = f k | k 1 u b , j ( ξ k )   and ω k | k 1 j J d , k | k 1 N k u u = ω k u b , j for j > J d , k | k 1 + N k | k 1 u u .

3.2. Update

If each target predicted density is a PMBM at time k , according to the conjugate prior property, the updated density is also a PMBM. The updated MBM is given by the formulas (31) and (32), which contain the MB components predicted by each target in the previous processing step and their related data associations.
f k | k i ( ξ k ( i ) | Z k ) = Y Z Ρ Y C C I j   ω k | k i , Ρ f k | k d , j ( ξ k ( i ) | Z k ) , j = 1 , , J d , k | k 1 + N k | k 1 u ,
ω k | k i , Ρ = ω k | k 1 j L Ρ L w Y Z Ρ Y C C I j   ω k | k i , Ρ L Ρ L w , j = 1 , , J d , k | k 1 + N k | k 1 u .
Here, the predicted likelihood L Ρ of partition and the predicted likelihood L w of each MB are as follows:
L Ρ = C C I j = | C C | > 1 l C C C I j = | C C | = 1 ( κ C C + l C ) ,
L w = j L C .
In formulas (33) and (34), for each MB component in each partition cell, the data association can be divided into the following three types:
(1) The measurement set Z can be divided into the union of two disjoint sets, Z = Y ( Z \ Y ) , where Y corresponds to the clutter and the previously undetected targets’ measurements, and Z \ Y corresponds to the previously detected targets’ measurements.
(2) For set Y , form a partition of P of non-empty cells C; each cell contains clutter or target-generated measurements.
(3) For set Z \ Y , the union of a set of index subsets   {   I j } can be used to represent the measurements related to the j-th MB component of the detected targets.
In the update step, the Bernoulli parameters of each detected target and the PPP parameters of each undetected target are determined, as given below.

3.2.1. Detected Targets

In this paper, there are two types of detected targets at time k : The targets detected for the first time in the measurement set Y , and the previously detected targets are detected again in the measurement set Z \ Y . The first detected targets are processed by the PPP model, and the second detected target is estimated with the existing MB. According to Equations (9) and (10), the spatial distribution f k | k d , j ( ξ k )   and existence probability r k | k d , j of each GGIW component of the detected target at time k are updated as follows:
f k | k d , j ( ξ k ( j ) ) = i I j ξ k i = X k d i C I j f k | k d , j , i ( ξ k ( j ) ) , j = 1 , , J d , k | k 1 ,
r k | k d , j = ω k | k 1 j C C I j r j , C j J ω k | k 1 j i C I j r j , C , j = 1 , , J d , k | k 1 .
For each partition cell, its predicted likelihood L C is related to the predicted likelihood of each GGIW component in the spatial density, such as
L C = { κ C C + l C , i f C I j =   , | C C | = 1 l C , i f C I j =   , | C C | > 1 1 r k | k 1 d , j p D + r k | k 1 d , j p D ( β k | k 1 d , j , C β k | k 1 d , j , C + 1 ) α k | k 1 d , j , C , i f C I j   , C C = r k | k 1 d , j L k u , j , C , i f C I j   , C C
The predicted likelihood L C is determined by C I j and the measurement cell   C C , as follows:
(1) When C I j =   , | C C | = 1 , the predicted likelihood L C is an approximation of the predicted likelihood of the cell composed of the undetected target-generated measurements and the clutter measurements in the first detection process. Since this cell contains the clutter and undetected targets’ measurements, L C is related to κ C C and l C , where κ C C is the predicted likelihood of clutter, and l C is the predicted likelihood of undetected target-generated measurements.
(2) When C I j =   , | C C | > 1 , it means the targets detected for the first time are divided into multiple cells. If there are undetected targets, they must be contained in the above cells, thus L C is only related to l C .
(3) When C I j   , C C = , it means L C is the predicted likelihood when the j-th MB component of the previously detected target is an empty set.
(4) When C I j   , C C , the predicted likelihood L C is the approximation when the j-th MB component of the previously detected target is not empty. The calculation L k u , j , C is presented in [39], Table 2.
The predicted likelihood l C of undetected targets can be defined as
l C = μ k | k 1 u j = 1 N k | k 1 u ω k | k 1 u , j L k u , j , C .
According to the different values of C I j and   C C , the likelihood probability r j , C can be defined as
r j , C = { l C κ C C + l C , i f C I j =   , | C C | = 1 1 , i f C I j =   , | C C | > 1 r k | k 1 d , j q D d , j , C 1 r k | k 1 d , j + r k | k 1 d , j q D j , C , i f C I j   , C C = 1 , i f C I j   , C C .
The effective probability of the missed detection can be defined as
q D d , j , C = 1 p D + p D ( β k | k 1 d , j , C β k | k 1 d , j , C + 1 ) α k | k 1 d , j , C .
The spatial distribution f k | k d , j , i ( ξ k ) of each measurement cell at time step k is as follows:
f k | k d , j , C ( ξ k ) = { ω k | k 1 u , j L k u , j , C j = 1 N k | k 1 u ω k | k 1 u , j L k u , j , C G G I W ( ξ k ; ζ k | k d , j , C ) , i f C C I j = 1 p D q D d , j , C G G I W ( ξ k ; ζ k | k 1 d , j , C ) + p D ( β k | k 1 d , j , C β k | k 1 d , j , C + 1 ) α k | k 1 d , j , C q D d , j , C   G ( γ k ; α k | k 1 u , j , C , β k | k 1 u , j , C + 1 ) × N ( x k ; m k | k 1 u , j , C , P k | k 1 u , j , C ) I W d ( X k ; v k | k 1 u , j , C , V k | k 1 u , j , C ) , i f C I j , C C = G G I W ( ξ k ; ζ k | k d , j , C ) , i f C I j , C C .
The second equation in (41) uses a gamma-mixture reduction to reduce its bi-modal GGIW distribution to a uni-modal GGIW distribution [44]. Detailed calculations of parameter ζ k | k d , j , C and prediction likelihood L k u , j , C in the GGIW update are provided in [41].

3.2.2. Undetected Targets

The detection probability p D is always less than 1, and the Poisson parameter γ k may be zero; these two cases may cause missed detection of targets. Therefore, we need to consider the above two situations, and the updated spatial density is
f k | k u , j ( ξ k ) = ( 1 p D ) ω k | k 1 u , j j = 1 N k | k 1 u q D u , j ω k | k 1 u , j G G I W ( ξ k ; ζ k | k 1 u , j ) +   p D ( β k | k 1 u , j β k | k 1 u , j + 1 ) α k | k 1 u , j ω k | k 1 u , j j = 1 N k | k 1 u q D u , j ω k | k 1 u , j G ( γ k ; α k | k 1 u , j , β k | k 1 u , j + 1 ) , × N ( x k ; m k | k 1 u , j , P k | k 1 u , j ) I W d ( X k ; v k | k 1 u , j , V k | k 1 u , j ) , j = 1 , , N k | k 1 u
where q D u , j is the effective probability of the undetected target, and μ k | k u is the PPP Poisson rate of the undetected targets,
q D u , j = 1 p D + p D ( β k | k 1 u , j β k | k 1 u , j + 1 ) α k | k 1 u , j ,
μ k | k u = μ k | k 1 u j = 1 N k | k 1 u q D u , j ,
The existence probability is
ω k | k u , j = ( 1 p D ) ω k | k 1 u , j j = 1 N k | k 1 u q D u , j ω k | k 1 u , j + p D ( β k | k 1 u , j β k | k 1 u , j + 1 ) α k | k 1 u , j ω k | k 1 u , j j = 1 N k | k 1 u q D u , j ω k | k 1 u , j , j = 1 , , N k | k 1 u .
Note: Most of the mathematical details of formulas (42) and (45) are given in the Appendix A. The Gaussian and inverse Wishart parameters in formula (42) are the same in both cases, but the gamma parameters are different. Therefore, formula (42) can be reduced by gamma mixing [44] to a single-peak GGIW distribution and the number of GGIW components remains unchanged.
The main contribution of the GGIW-MD-PMBM proposed in this paper is to recursively propagate the marginal distribution and existence probability of each target through formulas (21), (22), (29), (34), (35), (42), and (45), which is different from the GGIW-PMBM in [40] that it recursively propagates the joint target state distribution and existence probability.

3.3. Complexity Reduction and Data Association

In the tracking process, as the number of targets increases, the number of unknown data associations, the number of components in the MBM, and the number of parameters in the PMBM greatly increase, which brings the problem of “combination explosion”. It is necessary to approximate the target density through reduction methods. Possible reduction methods include gating, clustering, pruning, merging, and recycling. In this paper, we focused on developing the approximate method to recursively propagate the marginal distribution of targets in the prediction and update steps. The complexity reduction and the data association process used in this paper follow the same methods used in [40], Section V; hence the complexity reduction and the data association are briefly discussed.

4. Simulation

In this simulation, we compared the proposed GGIW-MD-PMBM filter with the GGIW-PMBM, GGIW-LMB, and GGIW-PHD filters for multi-target trajectories in the four classic scenarios, which are provided in the excellent sample MATLAB code [40,41,42,43]. Each extended target can be defined as x k = [ p k , v k ] T 4 , where p k 2 is each target’s position and v k 2 is the velocity. The motion model’s parameters used in Formulas (8)–(10) are
f ( x k ) = [ I 2 T s I 2 0 I 2 ] x k , Q = G σ a 2 I 2 G T , G = [ T s 2 2 I 2 T s I 2 ] , H k = [ I 2 0 2 ] ,
where σ a is the acceleration standard deviation and T s = 1 s is the sampling time.
To evaluate the four algorithms’ performance, we used the computing time and the generalized optimal sub-pattern assignment (GOSPA) metric [47], the cutoff parameter c = 10 , and the ordered parameter p = 1 , among which, for the distance measurement, we use the Gaussian Wasserstein Distance metric [48]. We divided the GOSPA metric into three categories: localization error, false detection error, and miss detection error. For GGIW-LMB filter and GGIW-PHD filter, we extracted the target states by taking the mean vector of all Bernoullis whose existence probability is larger than 0.5. For the GGIW-PMBM filter and GGIW-MD-PMBM filter, target state extraction was performed similarly, but only from the MB component with the highest weight.
In the four scenarios, the target detection probability p D , target survival probability p S , the clutter Poisson rate λ , and the measurement Poisson rate γ used in different scenarios are given in Table 1. In the first scenario, there are 100 time steps, and 27 highly time-varying targets are randomly generated in four positions; this scenario aims to compare the four algorithms’ tracking performances with a high clutter density and high target number scenario. In the second scenario, there are 100 time steps, and two targets are born well separated, move close to each other, and then split; this scenario aims to compare the four algorithms’ tracking performance when highly close to each other. In the third scenario, there are 10 time steps, and five targets are born closely at the same time; this scenario aims to compare the four algorithms’ tracking performances of handling dense birth. In the fourth scenario, there are 300 time steps, and two targets first get close and then they maneuver closely before splitting; this scenario aims to compare the four algorithms’ tracking performance of seriously handling the data association problem. The target trajectories of different scenarios are given in Figure 1.
Figure 2 gives the average GOSPA error of the four algorithms over 100 Monte Carlo (MC) trials. Table 2 shows the estimation errors and cycling time of the four algorithms in the four scenarios. From the simulation results in Figure 2 and Table 2, we can see that the proposed GGIW-MD-PMBM obtains the general minimum average GOSPA error, which outperforms the other three algorithms; GGIW-PMBM is the second, and GGIW-PHD obtains the highest average GOSPA error. As for the time consumption of each MC run, the GGIW-PHD has the lowest computational cost, GGIW-LMB is the second, and GGIW-MD-PMBM is the third. Compared with other results of normalized location error, the number of missed targets, and the number of false detection, the proposed GGIW-MD-PMBM algorithm is generally better than the other three algorithms. In the third scenario, the GGIW-PMBM filter outperformed the GGIW-MD-PMBM filter in terms of the location error and the number of false detections; this is because the targets in this scenario are intensively generated at the same location and at the same time. Compared with the joint state distribution (GGIW-PMBM), the marginal distribution (GGIW-MD-PMBM) has a larger deviation, but as time increases, the targets gradually move away, and the average GOSPA error of the GGIW-MD-PMBM filter is still lower than the GGIW-PMBM filter. The tracking performance of each algorithm in the fourth scenario is worse than in the other three scenarios; this is because the data association in the fourth scenario is worse than in the other three scenarios.

5. Conclusions

In this paper, we propose an efficient filter to solve the extended target tracking problem. Unlike the existing GGIW-PMBM filter that recursively propagates the joint state distribution of targets, the proposed GGIW-MD-PMBM filter recursively propagates the marginal distributions and the existence probabilities of each extended target. Comparing the other three extended target filters in four classic scenarios, the experimental results show that the proposed filter in this paper improves tracking accuracy and reduces computing time.

Author Contributions

Investigation, H.D.; Formal analysis, W.X.; Methodology, H.D.; Resources, H.D.; Simulation, H.D.; Writing—original draft, H.D.; Writing—review and editing, H.D. All authors have read and agreed to the published version of the manuscript.

Funding

This project was funded by the Shenzhen Basic Research Project, grant number JCYJ20170818102503604, and National Natural Science Foundation of China, grant number NO.61271107 and No.61703280.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Let p k , D e be the effective probability of target detection [35] (the set of target detections is nonempty with probability 1 e γ ( j ) ), then the updated PPP intensity corresponding to previously existing targets that are not detected is
D k | k u ( ξ k ) = j = 1 N k | k 1 u ( 1 p k , D e ) ω k | k u , j G G I W ( ξ k ; ζ k | k 1 u , j ) = j = 1 N k | k 1 u p ( j ) ( γ k ) N ( x k ; m k | k 1 u , j , P k | k 1 u , j ) I W d ( X k ; v k | k 1 u , j , V k | k 1 u , j )
where p ( j ) ( γ k ) is
p ( j ) ( γ k ) = ( 1 p D ) ω k | k 1 u , j G ( γ k ; α k | k 1 u , j , β k | k 1 u , j ) + p D ( β k | k 1 u , j β k | k 1 u , j + 1 ) α k | k 1 u , j ω k | k 1 u , j G ( γ k ; α k | k 1 u , j , β k | k 1 u , j + 1 ) .
At this time, all existence probabilities of the undetected targets can be normalized, and the existence probability and spatial density of each target can be as follows:
ω k | k u , j = ( 1 p D ) ω k | k 1 u , j j = 1 N k | k 1 u q D u , j ω k | k 1 u , j + p D ( β k | k 1 u , j β k | k 1 u , j + 1 ) α k | k 1 u , j ω k | k 1 u , j j = 1 N k | k 1 u q D u , j ω k | k 1 u , j , j = 1 , , N k | k 1 u
f k | k u , j ( ξ k ) = f 1 , k | k u , j ( ξ k ) + f 2 , k | k u , j ( ξ k ) , j = 1 , , N k | k 1 u
where f 1 , k | k u , j ( ξ k ) represents the spatial density in the case of a missed detection caused by p D < 1 , and f 2 , k | k u , j ( ξ k ) represents the spatial density in the case of a missed detection caused by γ k = 0 .
f 1 , k | k u , j ( ξ k ) = ( 1 p D ) ω k | k 1 u , j j = 1 N k | k 1 u q D u , j ω k | k 1 u , j G G I W ( ξ k ; ζ k | k 1 u , j )  
f 2 , k | k u , j ( ξ k ) = p D ( β k | k 1 u , j β k | k 1 u , j + 1 ) α k | k 1 u , j ω k | k 1 u , j j = 1 N k | k 1 u q D u , j ω k | k 1 u , j G ( γ k ; α k | k 1 u , j , β k | k 1 u , j + 1 ) N ( x k ; m k | k 1 u , j , P k | k 1 u , j ) I W d ( X k ; v k | k 1 u , j , V k | k 1 u , j )

References

  1. Vo, B.N.; Mallick, M.; Bar-shalom, Y.; Coraluppi, S.; Osborne, R., III; Mahler, R.; Vo, B.T. Multitarget tracking. Wiley Encycl. Electr. Electron. Eng. 2015. [Google Scholar] [CrossRef]
  2. Bar-Shalom, Y.; Li, X.R.; Kirubarajan, T. Estimation with Applications to Tracking and Navigation: Theory Algorithms and Software; John Wiley & Sons: Hoboken, NJ, USA, 2004. [Google Scholar]
  3. Granstrom, K.; Baum, M. Extended Object Tracking: Introduction, Overview and Applications. J. Adv. Inf. Fusion 2017, 12, 139–174. [Google Scholar]
  4. Gilholm, K.; Godsill, S.; Maskell, S.; Salmond, D. Poisson models for extended target and group tracking. In Proceedings of the Signal and Data Processing of Small Targets, San Diego, CA, USA, 31 July–4 August 2005; pp. 230–241. [Google Scholar]
  5. Granstrom, K.; Salmond, D. Spatial distribution model for tracking extended objects. IEEE Proc. Radarsonar Navig. 2005, 152, 364–371. [Google Scholar]
  6. Koch, J.W. Bayesian approach to extended object and cluster tracking using random matrices. IEEE Trans. Aerosp. Electron. Syst. 2008, 44, 1042–1059. [Google Scholar] [CrossRef]
  7. Feldmann, M.; FraNken, D.; Koch, W. Tracking of Extended Objects and Group Targets Using Random Matrices. IEEE Trans. Signal Process. 2010, 59, 1409–1420. [Google Scholar] [CrossRef]
  8. Feldmann, M.; Fränken, D. Advances on Tracking of Extended Objects and Group Targets using Random Matrices. In Proceedings of the 2009 12th International Conference on Information Fusion, Seattle, WA, USA, 6–9 July 2009; Volume 59, pp. 1029–1036. [Google Scholar]
  9. Granstrom, K.; Orguner, U. A phd Filter for Tracking Multiple Extended Targets Using Random Matrices. IEEE Trans. Signal Process. 2012, 60, 5657–5671. [Google Scholar] [CrossRef]
  10. Orguner, U. A variational measurement update for extended target tracking with random matrices. IEEE Trans. Signal Process. 2012, 60, 3827–3834. [Google Scholar] [CrossRef]
  11. Granstrom, K.; Orguner, U. A New Prediction Update for Extended Target Tracking with Random Matrices. IEEE Trans. Aerosp. Electron. Syst. 2014, 2, 324–336. [Google Scholar]
  12. Liu, J.; Guo, G. A Random Matrix Approach for Extended Target Tracking Using Distributed Measurements. IEEE Sens. J. 2019, 99, 1–12. [Google Scholar]
  13. Lan, J.; Li, X.R. Extended-Object or Group-Target Tracking Using Random Matrix with Nonlinear Measurements. IEEE Trans. Signal Process. 2019, 99, 1–15. [Google Scholar] [CrossRef]
  14. Wang, L.; Huang, Y.; Zhan, R.; Zhang, J. Joint Tracking and Classification of Extended Targets Using Random Matrix and Bernoulli Filter for Time-Varying Scenarios. IEEE Access 2019, 7, 129584–129603. [Google Scholar] [CrossRef]
  15. Lu, Z.; Hu, W.; Liu, Y.; Kirubarajan, T. Seamless Group Target Tracking Using Random Finite Sets. Signal Process. 2020, 176, 1–13. [Google Scholar] [CrossRef]
  16. Gao, L.; Jing, Z.; Li, M.; Pan, H. Bayesian approach to multiple extended targets tracking with random hypersurface models. IET Radar Sonar Navig. 2019, 13, 601–611. [Google Scholar] [CrossRef]
  17. Wang, X.; Li, H.G.; Kong, Y.B.; Pu, L.; Fan, P.F. SMC-PHD filter for extended target tracking based on star-convex random hypersurface models. Appl. Res. Comput. 2017, 34, 2144–2147. [Google Scholar]
  18. Liu, Z.P.; Liu, Y.J. Extended Target Tracking Algorithm Based on Star-Convex Random Hypersurface Models. Electron. Opt. Control 2017, 24, 72–76, 82. [Google Scholar]
  19. Li, Y.W. Extended Target Tracking Algorithm Based on Random Hypersurface Model with Glint Noise. Chem. Eng. Trans. 2017, 59, 685–690. [Google Scholar]
  20. Baum, M.; Hanebeck, U.D. Random Hypersurface Models for extended object tracking. In Proceedings of the IEEE International Symposium on Signal Processing and Information Technology, Ajman, UAE, 14–17 December 2009; IEEE: Piscataway, NJ, USA, 2010; pp. 178–183. [Google Scholar]
  21. Baum, M.; Hanebeck, U.D. Extended Object Tracking with Random Hypersurface Models. IEEE Trans. Aerosp. Electron. Syst. 2014, 50, 149–159. [Google Scholar] [CrossRef] [Green Version]
  22. Wahlström, N.; Özkan, E. Extended Target Tracking Using Gaussian Processes. IEEE Trans. Signal Process. 2015, 63, 4165–4178. [Google Scholar] [CrossRef] [Green Version]
  23. Mahler, R. Statistical Multisource-Multitarget Information Fusion; Artech House: Norwood, MA, USA, 2007; pp. 1–48. [Google Scholar]
  24. Lian, F.; Han, C.; Liu, W.; Chen, H. Joint spatial registration and multi-target tracking using an extended probability hypothesis density filter. IET Radar Sonar Navig. 2011, 5, 441–448. [Google Scholar] [CrossRef]
  25. Li, G.; Zhang, H.Q.; Wang, Y. An Efficient Multi-Target Tracking Algorithm Using Gaussian Mixture Probability Hypothesis Density Filter. In Proceedings of the 2018 IEEE CSAA Guidance, Navigation and Control Conference (CGNCC), Xiamen, China, 10–12 August 2018; IEEE: Piscataway, NJ, USA, 2020. [Google Scholar]
  26. Liu, L.; Ji, H.; Zhang, W.; Liao, G. Multi-Sensor Multi-Target Tracking Using Probability Hypothesis Density Filter. IEEE Access 2019, 7, 67745–67760. [Google Scholar] [CrossRef]
  27. Li, C.; Wang, R.; Hu, Y.; Wang, J. Cardinalised probability hypothesis density tracking algorithm for extended objects with glint noise. IET Sci. Meas. Technol. 2016, 10, 528–536. [Google Scholar] [CrossRef]
  28. 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]
  29. Li, B.; Yi, H.W.; Li, X.H. Innovative unscented transform–based particle cardinalized probability hypothesis density filter for multi-target tracking. Meas. Control Lond. Inst. Meas. Control 2019, 52, 1567–1578. [Google Scholar] [CrossRef]
  30. Chen, J.G.; Wang, X.H.; Ma, L.L.; Zhang, X.D.; Gong, L.M. Cardinalized Probability Hypothesis Density Smoother Using Piece-wise RTS. Comput. Eng. Appl. 2019, 55, 50–55, 95. [Google Scholar]
  31. Gostar, A.K.; Hoseinnezhad, R.; Bab-Hadiashar, A. Robust Multi-Bernoulli Sensor Selection for Multi-Target Tracking in Sensor Networks. IEEE Signal Process. Lett. 2013, 20, 1167–1170. [Google Scholar] [CrossRef]
  32. Park, W.J.; Park, C.G. Multi-target Tracking Based on Gaussian Mixture Labeled Multi-Bernoulli Filter with Adaptive Gating. In Proceedings of the 2019 First International Symposium on Instrumentation, Control, Artificial Intelligence, and Robotics (ICA-SYMP), Bangkok, Thailand, 16–18 January 2019; pp. 226–229. [Google Scholar]
  33. Cament, L.; Correa, J.; Adams, M.; Pérez, C. The Histogram Poisson, Labeled Multi-Bernoulli Multi-Target Tracking Filter. Signal Process. 2020, 176, 837–851. [Google Scholar] [CrossRef]
  34. Zhu, Y.; Wang, J.; Liang, S. Multi-Objective Optimization Based Multi-Bernoulli Sensor Selection for Multi-Target Tracking. Sensors 2019, 19, 1884–2021. [Google Scholar] [CrossRef] [Green Version]
  35. Mahler, R. PHD filters for nonstandard targets, I: Extended targets. In Proceedings of the 2009 12th International Conference on Information Fusion, Seattle, WA, USA, 6–9 July 2009; pp. 915–921. [Google Scholar]
  36. Granstrom, K.; Lundquist, C.; Orguner, U. A Gaussian mixture PHD filter for extended target tracking. In Proceedings of the International Conference on Information Fusion, Edinburgh, UK, 26–29 July 2011. [Google Scholar]
  37. Han, Y.; Zhu, H.; Han, C.Z. A Gaussian-mixture PHD filter based on random hypersurface model for multiple extended targets. In Proceedings of the 16th International Conference on Information Fusion, Istanbul, Turkey, 9–12 July 2013. [Google Scholar]
  38. Lundquist, C.; Granstrom, K.; Orguner, U. An extended target CPHD filter and a gamma Gaussian inverse Wishart implementation. J. Sel. Top. Signal Process. 2013, 7, 472–483. [Google Scholar] [CrossRef] [Green Version]
  39. Orguner, U.; Lundquist, C.; Granstrom, K. Extended target tracking with a cardinalized probability hypothesis density filter. In Proceedings of the 14th International Conference on Information Fusion, Chicago, IL, USA, 5–8 July 2011. [Google Scholar]
  40. Granstrom, K.; Fatemi, M.; Svensson, L. Poisson multi-Bernoulli conjugate prior for multiple extended object estimation. IEEE Trans. Aerosp. Electron. Syst. 2018, 3, 30–45. [Google Scholar]
  41. Granstrom, K.; Fatemi, M.; Svensson, L. Gamma Gaussian inverse-Wishart Poisson multi-Bernoulli filter for extended target tracking. In Proceedings of the 19th International Conference on Information Fusion, Heidelberg, Germany, 5–8 July 2016. [Google Scholar]
  42. Beard, M.; Reuter, S.; Granström, K.; Vo, B.T.; Vo, B.N.; Scheel, A. A Generalised Labelled Multi-Bernoulli Filter for Extended Multi-target Tracking. In Proceedings of the 18th International Conference on Information Fusion, Washington, DC, USA, 6–9 July 2015. [Google Scholar]
  43. Beard, M.; Reuter, S.; Granström, K.; Vo, B.T.; Vo, B.N.; Scheel, A. Multiple extended target tracking with labelled random finite sets. IEEE Trans. Signal Process. 2016, 64, 1638–1653. [Google Scholar] [CrossRef] [Green Version]
  44. Granstrom, K.; Orguner, U. Estimation and maintenance of measurement rates for multiple extended target tracking. In Proceedings of the 15th International Conference on Information Fusion, Singapore, 9–12 July 2012; pp. 2170–2176. [Google Scholar]
  45. Liu, Z.X.; Xie, W.X. Multi-target Bayesian filter for propagating marginal distribution. Signal Process. 2014, 105, 328–337. [Google Scholar] [CrossRef]
  46. Bar-Shalom, Y.; Willett, P.K.; Tian, X. Tracking and Data Fusion: A Handbook of Algorithms; YBS Publishing: Storrs, CT, USA, 2011. [Google Scholar]
  47. Rahmathullah, A.S.; García-Fernández, Á.F.; Svensson, L. Generalized optimal sub-pattern assignment metric. In Proceedings of the 20th International Conference on Information Fusion (Fusion), Xi’an, China, 10–13 July 2017; pp. 1–8. [Google Scholar]
  48. Yang, S.S.; Baum, M.; Granstrom, K. Metrics for performance evaluation of elliptic extended object tracking methods. In Proceedings of the International Conference on Multisensor Fusion and Integration for Intelligent Systems (MFI), Baden-Baden, Germany, 19–21 September 2016; pp. 523–528. [Google Scholar]
Figure 1. True target trajectories of four scenarios. (a) Twenty-seven targets, (b) separate/close/split, (c) dense birth, (d) nonlinear maneuver.
Figure 1. True target trajectories of four scenarios. (a) Twenty-seven targets, (b) separate/close/split, (c) dense birth, (d) nonlinear maneuver.
Sensors 20 05387 g001
Figure 2. Comparison of the simulation results of the four algorithms of GOSPA. (a) Twenty-seven targets, (b) separate/close/split, (c) dense birth, (d) nonlinear maneuver.
Figure 2. Comparison of the simulation results of the four algorithms of GOSPA. (a) Twenty-seven targets, (b) separate/close/split, (c) dense birth, (d) nonlinear maneuver.
Sensors 20 05387 g002
Table 1. The parameters used in the four scenarios.
Table 1. The parameters used in the four scenarios.
Scenario p D p S λ γ
Scenario 10.900.9960{7,8,9}
Scenario 20.980.9910{10,20}
Scenario 30.900.992010
Scenario 40.980.9910{10,20}
Table 2. Simulation results: the sum of estimation errors and cycling time; GO-GOSPA; LE-normalized location error; NM-number of missed targets; NF-number of false detection; CT-cycling time.
Table 2. Simulation results: the sum of estimation errors and cycling time; GO-GOSPA; LE-normalized location error; NM-number of missed targets; NF-number of false detection; CT-cycling time.
ScenariosGGIW-PMBMGGIW-LMBGGIW-PHDGGIW-MD-PMBM
Scenario 1GO732.81246.72873.2729
LE56.576.660.756.3
NM61.3481.02311.660.7
NF141.496.2108.5137
CT62.58.04.246.2
Scenario 2GO550.71818.24699.5550.5
LE268.0133.9562.3265.1
NM5.61479.5997.55.0
NF16.9193.83083.210.8
CT18.67.00.311.4
Scenario 3GO59.2134.8280.458.3
LE9.311.523.19.5
NM11.273.8174.710.6
NF2.112.131.82.5
CT1.20.40.11.1
Scenario 4GO2835.66266.36257.82236.7
LE1011.2869.3176.81002.8
NM253.63770.65601.2139.4
NF175.01445.2470.8106.2
CT49.76.52.240.6

Share and Cite

MDPI and ACS Style

Du, H.; Xie, W. Extended Target Marginal Distribution Poisson Multi-Bernoulli Mixture Filter. Sensors 2020, 20, 5387. https://doi.org/10.3390/s20185387

AMA Style

Du H, Xie W. Extended Target Marginal Distribution Poisson Multi-Bernoulli Mixture Filter. Sensors. 2020; 20(18):5387. https://doi.org/10.3390/s20185387

Chicago/Turabian Style

Du, Haocui, and Weixin Xie. 2020. "Extended Target Marginal Distribution Poisson Multi-Bernoulli Mixture Filter" Sensors 20, no. 18: 5387. https://doi.org/10.3390/s20185387

APA Style

Du, H., & Xie, W. (2020). Extended Target Marginal Distribution Poisson Multi-Bernoulli Mixture Filter. Sensors, 20(18), 5387. https://doi.org/10.3390/s20185387

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