Next Article in Journal
Hyperspectral Super-Resolution Reconstruction Network Based on Hybrid Convolution and Spectral Symmetry Preservation
Next Article in Special Issue
Fast Frequency-Diverse Radar Imaging Based on Adaptive Sampling Iterative Soft-Thresholding Deep Unfolding Network
Previous Article in Journal
Target Detection in Hyperspectral Remote Sensing Image: Current Status and Challenges
Previous Article in Special Issue
An Anti-Intermittent Sampling Jamming Technique Utilizing the OTSU Algorithm of Random Orthogonal Sub-Pulses
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gaussian Process Gaussian Mixture PHD Filter for 3D Multiple Extended Target Tracking

1
School of Electronics and Information Engineering, Beihang University, Beijing 100191, China
2
School of Information and Electronics, Beijing Institute of Technology, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(13), 3224; https://doi.org/10.3390/rs15133224
Submission received: 19 May 2023 / Revised: 17 June 2023 / Accepted: 20 June 2023 / Published: 21 June 2023
(This article belongs to the Special Issue Advanced Radar Signal Processing and Applications)

Abstract

:
This paper addresses the problem of tracking multiple extended targets in three-dimensional space. We propose the Gaussian process Gaussian mixture probability hypothesis density (GP-PHD) filter, which is capable of tracking multiple extended targets with complex shapes in the presence of clutter. Our approach combines the Gaussian process regression measurement model with the probability hypothesis density filter to estimate both the kinematic state and the shape of the targets. The shape of the extended target is described by a 3D radial function and is estimated recursively using the Gaussian process regression model. Furthermore, we transform the recursive Gaussian process regression problem into a state estimation problem by deriving a state space model such that the estimation of the extent can be integrated into the kinematic part. We derive the predicted likelihood function of the PHD filter and provide a closed-form Gaussian mixture implementation. To evaluate the performance of the proposed filter, we simulate a typical extended target tracking scenario and compare the GP-PHD filter with the traditional Gamma Gaussian Inverse-Wishart PHD (GGIW-PHD) filter. Our results demonstrate that the proposed algorithm outperforms the GGIW-PHD filter in terms of estimating both kinematic states and shape. We also investigate the impact of the measurement rates on both filters; it is observed that the proposed filter exhibits robustness across various measurement rates, while the GGIW-PHD filter suffers under low-measurement-rate conditions.

1. Introduction

The tracking of moving targets within surveillance areas plays a crucial role in various fields such as automotive driving and robotic systems [1]. With the increasing demand for accurate tracking, high-resolution sensors such as Radar and Lidar have gained significant importance. In traditional tracking scenarios, e.g., air traffic control, each target produces at most one detection. However, in many recent applications, such as autonomous driving, this is not the case, because the resolution of the sensor is high enough that the detected object may occupy multiple resolution cells and thus generate multiple measurements. In the literature [2], such scenarios are defined as the extended object tracking (EOT) problem, which is also referred to as the extended target tracking (ETT) problem. In the ETT problem, the shape of the target is not negligible, requiring the tracking algorithm to estimate both its shape and kinematic states. Moreover, the data association and the track management in the ETT problem are far more complicated than in the point target tracking scenario. The ETT has become a key problem in the automotive system. Given its importance in automotive systems, extensive research has been conducted on the ETT problem, utilizing the random finite set (RFS) theory and nonlinear estimation methods [3,4,5,6]. Recently, deep learning-based methods have also been proposed to address the ETT problem [7,8]. A comprehensive tutorial for the ETT problem is given in [1].
The ETT problem can be divided into two parts: the non-point target measurement model and the multiple-target tracking (MTT) problem. Some commonly used non-point target measurement models are the random matrix (RM) [9,10] model and the random hypersurface model (RHM) [11]. The RM model, first introduced in [9], assumes that the measurements of extended targets are distributed over the entire ellipse surface, following a spatial Gaussian distribution. A semi-positive definite (SPD) matrix is employed to describe the contours of the ellipse. However, the original RM model neglects the measurement errors caused by the sensor, which can accumulate and lead to an overestimation of the object’s shape. To address this limitation, Koch enhances the measurement model in [10] by taking the measurement error into account. In addition, the interacting multiple model (IMM) is integrated into the RM model, enabling the filter to track maneuvering targets. Further investigation into the IMM model’s utilization for maneuvering extended target tracking is presented in [12]. A joint detection, tracking, and classification (JDTC) approach for multiple extended objects is proposed in [13]. This also uses an RM to model the ellipse. In many situations, it is not accurate enough to model the target as an ellipse. Therefore, Lan [14] proposes an approach that uses multiple sub-ellipses to model the extended targets with complex shapes. Another common method assumes that extended targets are star-convex and use the RHM to derive the measurement model. RHM-based methods assume that the measurement sources have certain spatial distributions and are located on differently scaled inner surfaces of the target. The implicit measurement model is then updated with a non-linear filter such as an extended Kalman filter (EKF). Based on the RHM, the radial function is introduced in [15] to explicitly describe the surface of the non-point target, and the Fourier series expansion is used to approximate the radial function. Sun et al. [16] propose a filter to track number-variable maneuvering extended targets using the RHM and the IMM. The B-spline model is another approach to modeling star-convex targets and is first introduced by [17]. Compared to the RHM and the RM model, the B-spline model is capable of describing arbitrary shape extended targets using the control points on the B-spline. In [18], the B-spline model is used under the Poisson multi-Bernoulli mixture (PMBM) filter framework to track non-ellipsoidal targets. The B-spline model is also used to track elongated targets such as trains and pedestrians [19]. The Gaussian process (GP) has been widely used in machine learning and other fields because of its easy-to-compute posterior probability and excellent analytical properties [20]. It is first proposed in [21] to use the GP regression to estimate the radial function. The GP model can describe arbitrary shaped extended targets and is computationally more efficient than the B-spline model. At the same time, the estimation of the target shape is more accurate than the ellipse approximation models. A new temporal covariance kernel for GP regression is proposed in [22] that could improve the accuracy of object shape estimation. The dependency of connected extended targets is utilized in [23] and the GP model is used to estimate the shape of the extended target.
The MTT problem is challenging due to the track associations, miss detections, etc. However, the RFS theory [24] proposed by Mahler provides a rigorous treatment of the MTT. Their method is completely free of explicit data associations, unlike traditional MTT techniques such as the joint probabilistic data association (JPDA) filter and the multiple hypothesis tracker (MHT). Furthermore, with the improvement of hardware [25,26,27,28], it is possible to employ algorithms based on the RFS theory. The optimal Bayesian filter derived from the RFS theory involves high-dimensional integrals. Therefore, Mahler proposes several theoretical approximation methods, such as the PHD filter [29] and the MB filter [24]. Among these methods, the PHD filter stands out as the most computationally efficient and memory-friendly, making it a practical choice given current computational resources. The first attempt to address the ETT problem within the framework of the PHD filter is proposed by [2]. A Gaussian mixture implementation of [2] is given in [30] called the extended target Gaussian mixture PHD filter (ET-GM-PHD). Nonetheless, both methods in [2,30] only estimate the kinematic states of the target and omit the shape which contains crucial information about the target. In [31], the RM approach is first used in the PHD framework to estimate the shape of an unknown number of extended targets in the presence of clutter and miss detections. Granström [32] further considers target spawning and combination and proposes the GGIW-PHD filter, which is the state-of-the-art ETT tracker and has been successfully tested in a real marine surveillance scenario [33].
Although 3D data, generated by sensors such as depth cameras and Lidars, are becoming increasingly accessible, there have been few attempts to address the ETT problem in 3D space. The relationship between the ellipse and symmetrical SPD matrix is analyzed and the ellipse fitting approach (EFA) is proposed [34]. In [35], a generalization of the RHM to 3D space is proposed; however, only cylinder targets are suitable for this method. Current approaches for 3D ETT usually use simple geometries such as cylinders and ellipses and are not able to model complexly shaped extended targets. The GP model approach to tracking 3D targets is proposed in [36,37], but it only considers the single-target situation. Therefore, to solve the 3D ETT problem, the combination of a complex-target model and multiple-target tracking filters, such as the PHD filter, is required.
In this paper, we propose a method that uses a 3D radial function to describe the shape of the target. The GP regression model is adopted to estimate the value of the radial function. We transform the regression problem into a state estimation problem. This approach allows us to use the Bayesian filter paradigm to estimate the unknown shape of the target and can be integrated into the kinematic estimation process. The computationally efficient PHD filter is combined with the measurement model to track an unknown number of extended targets in the presence of clutter. The kinematic states and the shape are then inferred by the EKF. The proposed filter is capable of tracking multiple extended targets with complex shapes in the presence of clutter. To demonstrate the performance of the algorithms, we simulate point cloud measurements in MATLAB for two dynamic objects with different shapes in the presence of clutter. We evaluate the performance of the filter using both the optimal sub-pattern assignment (OSPA) [38] and the intersection-over-union (IoU) value. The proposed method outperforms the traditional GGIW-PHD filter in the simulation experiment. We also test the stability of the two filters at different measurement rates. It is shown that the proposed filter performs better than the GGIW-PHD filter under low-measurement-rate conditions.

2. Materials and Methods

2.1. Extended Target Tracking Problem Formulation

In the multiple extended target tracking problem, the states of extended targets at time k can be represented by an RFS:
X k = { ξ k i } i = 1 N x , k F ( X ) ,
where N x , k is the number of unknown targets, ξ k i is the kinematic and extent states of the ith extended target, and F ( X ) is the multiple targets’ state space. We use operation | | to denote the cardinality of a set, and then X k = N x , k . Measurements at time k can also be represented by an RFS Z k = { z k ( j ) } j = 1 N z , k , where N z , k is the number of measurements. The ETT problem is to estimate the posterior density π k ( X k | Z 1 : k ) of the target’s state using measurements Z 1 : k = { Z 1 , , Z k } from time 1 to k .

2.2. The PHD Filter

The multi-target Bayesian filter based on RFS theory includes the prediction and the update process. According to the Chapman–Kolmogorov equation, the posterior density at time k is
π k | k 1 ( X k | Z 1 : k 1 ) = f k | k 1 ( X k | X k 1 ) π k 1 ( X k 1 | Z 1 : k 1 ) δ X k 1 ,
where f k | k 1 ( X k | X k 1 ) is the probability transition density. Given the predicted density and measurements, the updated density at time k is
π k ( X k | Z 1 : k ) = g k ( Z k | X k ) π k 1 ( X k 1 | Z 1 : k 1 ) g k ( Z k | X k ) π k | k 1 ( X k | Z 1 : k 1 ) δ X k ,
where g k ( Z k | X k ) is the likelihood density of the measurements. There is no closed solution to Equation (3) that involves high-dimensional integral. The PHD filter is an approximation of the optimal Bayesian filter that estimates multiple target states by passing the first moment of the multi-object probability density, or PHD. The PHD can be calculated by multi-object density
D ( x ) = π ( { x } X ) δ X .
The PHD filter contains a prediction and an update process to estimate the density D k + 1 | k + 1 ( ) of the target’s RFS X k . The prediction equation of the extended target PHD filter is given in [2] as
D k + 1 | k ( ξ k + 1 ) = p s ( ξ k ) p k + 1 | k ( ξ k + 1 | ξ k ) D k | k ( ξ k ) d ξ k + D k + 1 b ( ξ k + 1 ) ,
where p s ( ) is the survival probability, p k + 1 | k ( ξ k + 1 | ξ k ) is the state transition density and D k + 1 b ( ξ k + 1 ) is the birth PHD. The update equation for the extended target PHD filter is given as follows
D k | k ( ξ k | Z k ) = L Z k ( ξ k ) D k | k 1 ( ξ k | Z k 1 ) ,
where L Z k ( ) is the pseudo-measurement likelihood function. L Z k ( ) is given as follows
L Z k ( ξ k ) ( 1 e γ ( ξ k ) ) p D ( ξ k ) + e γ ( ξ k ) p D ( ξ k ) P Z k ω P W P γ ( ξ k ) | W | d W z k W ϕ z k ( ξ k ) λ k c k ( z k ) ,
where
  • γ ( ξ k ) is the expected number of measurements from targets;
  • p D ( ξ k ) is the probability of detection;
  • c k ( z k ) is the spatial distribution of clutter over the surveillance area;
  • ϕ z k ( ξ k ) is the single-target measurement likelihood function;
  • P Z denotes that P partitions the measurement set Z k into different nonempty measurement cells W . W P denotes that the cell W is in partition P .
d W and L k ( j , W ) are non-negative coefficients. We use δ denoting the Kronecker delta function, and the coefficients are given as follows
w P = W P d W P Z k W P d W ,
d W = δ W , 1 + D k + 1 | k p d γ k | W | exp ( γ k ) z k W ϕ z k ( ) λ k c k ( z k ) .

3. Gaussian Process Regression Measurement Model

The Gaussian process has been widely used in tasks such as classification and regression. A Gaussian process f ( u ) can be defined by its mean function μ ( u ) and the correlation function C ( u , u ) , which are given as
μ ( u ) = E f ( u ) ,
C ( u , u ) = E ( f ( u ) μ ( u ) ) ( f ( u ) μ ( u ) ) T .
Therefore, the GP f ( u ) can be denoted as
f ( u ) ~ G P ( μ ( u ) , C ( u , u ) ) .
The Gaussian process could be regarded as the generalization in the infinite dimension of a joint Gaussian distribution, and so the function values of N input u 1 u N are Gaussian distributed as
f ( u 1 ) f ( u N ) ~ N ( μ , K ) ,
where
μ = μ ( u 1 ) μ ( u N ) ,   K = C ( u 1 , u 1 ) C ( u 1 , u N ) C ( u N , u 1 ) C ( u N , u N ) .
We assume that a training set is given, whose input is u = [ u 1 , , u N ] T and output is z = [ z 1 , z N ] T . The measurement model of the training set is
z k = f ( u k ) + e k ,   e k ~ N ( 0 , R ) .
The Gaussian process regression is used to estimate the output f = [ f ( u 1 f ) , f ( u N f ) ] T of some given test set inputs u f = [ u 1 f , , u N f ] T . The training set, together with the measurement model, leads to the following distribution
z f ~ N 0 , K ( u , u ) + I N R K ( u , u f ) K ( u f , u ) K ( u f , u f ) .
The conditional distribution of the joint Gaussian distribution is still a Gaussian distribution that gives
p ( f | z ) = N ( A z , P ) ,
where
A = K ( u f , u ) K y ,
P = K ( u f , u f ) K ( u f , u ) K y 1 K ( u , u f ) ,
K y = K ( u , u ) + I N R .
The measurements are collected recursively rather than in the batch form in the multiple-target tracking problem, and so the Gaussian process regression needs to be improved in order to estimate the output of the test set recursively. Ref. [39] introduces basis points which are considered to be the inputs of the training set. This article adopts the same approach as [39], where a state-space model is also introduced to estimate the output recursively. u f = [ u 1 f , , u N f ] T and x k f denote the input and the output of the test set, respectively, and so the state-space model is
x k + 1 f = x k f ,
z k = H f ( u k ) x k f + e k f ,   e k f ~ N 0 , R f ( u k ) ,
where z k is the measurement at time k . We assume x f is sufficient statistics for z k . Under this assumption
p ( z k | x f , z 1 : k 1 ) p ( z k | x f ) ,
the recursive optimal Bayesian filter gives the posterior probability density function (PDF) p ( x f | z 1 : k ) as
p ( x f | z 1 : k ) p ( z k | x f , z 1 : k 1 ) p ( x f | z 1 : k 1 ) p ( z k | x f ) p ( x f | z 1 : k 1 ) .
According to Equation (17), the joint distribution of z k and x f is
z k x f ~ N 0 , K ( u k , u k ) + R K ( u k , u f ) K ( u f , u k ) K ( u f , u f ) .
Therefore, the likelihood function and the prior PDF are
p ( z k | x f ) = N ( z k ; H k f x f , R k f ) ,
p ( x f ) = N ( 0 , P 0 f ) .
where
H f ( u k ) = K ( u k , u f ) [ K ( u f , u f ) ] 1 ,
R f ( u k ) = Κ ( u k , u k ) + R K ( u k , u f ) [ K ( u f , u f ) ] 1 K ( u f , u k ) ,
P 0 f = K ( u f , u f ) .
The advantage of the method is that the extent states can be inserted into the kinematic states and be estimated using a single-state space model.

Measurement Model

In this article, we assume that the measurement sources are all on the object’s surface. A 3D radial function r = f ( θ , ϕ ) is used to model the shape of the target, whose inputs ( θ , ϕ ) are the azimuth and the elevation. For simplicity, we use γ to denote the pair ( θ , ϕ ) . The function value r of the radial function is the distance between the basis points and the center of the extended target. The radial function is assumed to be a GP, denoted as
f ( γ ) ~ G P μ r , C ( γ , γ ) + σ r 2 .
The mean function of the GP model is assumed to be some unknown constant μ r , which is Gaussian-distributed, and we assume that μ r = 0 as in [37]. The correlation function in [37] is used, which is
C ( γ , γ ) = σ f 2 exp d 2 ( γ , γ ) 2 l 2 ,
where σ f 2 represents the prior variance, l is the length scale, and d ( γ , γ ) calculates the relative distance between two inputs, given as follows
d ( γ , γ ) = arccos cos ( ϕ ) cos ( ϕ ) cos ( θ ) cos ( θ ) + cos ( ϕ ) cos ( ϕ ) sin ( θ ) sin ( θ ) + sin ( ϕ ) sin ( ϕ ) .
The output range of the distance function is [ 0 , π ] , and any coincident input is mapped to 0 , while the opposite pair is mapped to π . The employed distance d ( γ , γ ) is specified to imply a higher correlation for closer regions compared to separated regions.
Together with the radial function and the GP model, both the kinematic model and measurement model can be derived. We denote the extended target state as x k [ x ¯ k T , f k T ] T , where x ¯ k = [ x k t T , x k r T ] T is the extended target kinematic state including the translation vector x k t and the orientation vector x k r . The values of the radial function are denoted as f k . The kinematic model of the extended target is
x k + 1 = F k x k + w k ,   w k ~ N ( 0 , Q k ) ,
F k = F ¯ k F k f ,   Q k = Q ¯ k Q k f ,
F ¯ k = F k t F k r ,   Q ¯ k = Q k t Q k r .
The constant velocity (CV) model is considered for the translation of the extended target in this article, which gives
x k t [ x k c , x ˙ k c ] T ,
F k t = 1 T 0 1 I 3 ,
Q k t = T 3 / 3 T 2 / 2 T 2 / 2 T σ x I 3 .
A 3 × 3 orthogonal matrix called a rotation matrix is usually used to describe the orientation of the target. However, a rotation matrix is singular or discontinuous at some points. Quaternion is used in this article instead of a rotation matrix to describe the orientation of the target. A quaternion q is a 4D vector which contains two parts: a three-dimensional vector and a constant, denoted as
q = q ¯ T q 4 T .
The translation between the rotation matrix and the quaternion is
R G L ( q k ) = ( q 4 2 q ¯ T q ¯ ) I 3 + 2 q ¯ T q ¯ 2 q 4 [ q ¯ × ] T ,
where R ( ) represents the rotation matrix and [ q ¯ × ] is the cross-product matrix
[ q ¯ × ] = 0 q 3 q 2 q 3 0 q 1 q 2 q 1 0 .
The multiplication of two quaternions q and q is defined as
q q = q 4 q + q 4 q q ¯ × q ¯ q 4 q 4 q ¯ T q ¯ .
A useful property of quaternion is
R ( q ) R ( q ) = R ( q q ) ,
which indicates that the product of two rotation matrices can be represented by the product of quaternions. Therefore, the continuous rotation of a target can be represented by the product of quaternions. Ref. [40] introduces an approach to estimate the orientation of the target recursively using quaternions based on the property of Equation (44). This method can be augmented in the process of estimating the target’s kinematic states. The orientation of the target at time k can be expressed as
q = δ q ( a ) q ref
where q ref is the reference orientation and δ q ( a ) is the deviation from the reference orientation. To estimate the orientation recursively, the estimated orientation at time k 1 is used as the reference orientation. a is the deviation vector, and δ q ( a ) is defined using Rodrigue parameterization as
δ q ( a ) = 1 4 + a 2 a 2 .
The core idea is to estimate the deviation vector a , and the orientation at any moment can be calculated together with the reference orientation. Ref. [41] gives the dynamic model of the deviation vector
a ˙ = I 3 + 1 4 a a T + 1 2 [ a × ] ω I 3 + 1 2 [ a × ] ω .
We denote the rotation vector as x k r a T ω T T , and the dynamic model of the rotational part of the extended target can be derived using Equation (47) as
x ˙ = I 3 + 1 2 [ a × ] ω 0 3 × 1 + 0 3 I 3 α ,
where α is the acceleration of the rotation vector, and is assumed to be the Gaussian white noise in the article. Equation (48) indicates a nonlinear system, and therefore the first-moment Taylor expansion is used to derive the approximated linear model. The first-order Taylor expansion of the estimated orientation x r = x k | k r at time k is given as
x ˙ r = f ( x r ) + B α f ( x ^ r ) + A k r ( x r x ^ r ) .
The deviation vector is set to zero after each measurement update, and so the Taylor approximation in Equation (49) is
A k r = d d x r f ( x r ) x r = x ^ k | k r = 1 2 [ ω ^ k | k × ] I 3 0 3 0 3 .
The discrete system is
x k + 1 r = F k r x k r + w k r ,   w k r ~ N ( 0 , Q k r ) ,
where
F k r = exp ( A k T ) = exp T 2 [ ω ^ k | k × ] T exp T 2 [ ω k | k × ] 0 3 I 3 .
Substituting Equation (42) into Equation (52) gives
exp T 2 [ ω ^ k | k × ] = I 3 + sin T 2 ω ^ k | k ω ^ k | k [ ω ^ k | k × ] + 1 cos T 2 ω ^ k | k ω ^ k | k 2 [ ω ^ k | k × ] 2 ,
where T is the sampling time. The process noise covariance matrix is
Q k r = G k Σ α G k T .
Because the deviation vector is set to zero, G k is
G k = 0 T exp τ 2 ω ^ k | k × d τ 0 T exp τ 2 ω ^ k | k × d τ 0 3 0 T I 3 d τ .
The following are the details of the matrices
0 T exp τ 2 ω ^ k | k × d τ = T I 3 + 2 1 cos T 2 ω ^ k | k ω ^ k | k 2 ω ^ k | k × + T 2 ω ^ k | k sin T 2 ω ^ k | k ω ^ k | k 2 ω ^ k | k × 2 ,
0 T τ exp ( τ 2 [ ω ^ k | k × ] ) d τ = T 2 2 I 3 + 1 | ω ^ k | k | 2 ( 4 | ω ^ k | k | sin ( T 2 | ω ^ k | k | ) 2 T cos ( T 2 | ω ^ k | k | ) ) × [ ω ^ k | k × ] + 1 | ω ^ k | k | 2 ( T 2 2 + 2 T | ω ^ k | k | sin ( T 2 | ω ^ k | k | ) + 4 | ω ^ k | k | 2 ( cos ( T 2 | ω ^ k | k | ) 1 ) ) [ ω ^ k | k × ] 2 ,
0 T I 3 d τ = T I 3 .
We assume that the shape of the target does not change, and that therefore the dynamic model for the extent state of the extended target is
f k + 1 = f k + w k , w k ~ N ( 0 , Q k f ) ,
Q k f = ( 1 λ 1 ) P k | k f ,
which indicates that F k f = I .
Two coordinate systems are used in the article to derive the measurement model. The first one is the global coordinate denoted as the upper letter G , and the second one is the local coordinate of the target denoted as the upper letter L . Figure 1 shows the two coordinates.
We assume n k measurements { z k , i } i = 1 n k at time k are received, and so the ith measurement z k , i can be expressed as
z k , i G = x k c + p k , i f ( γ k , i L ) + e ¯ k ,   e ¯ k ~ N ( 0 , R ¯ )
where p k , i = z k , i G x k c z k , i G x k c is the direction vector. Measurements in the local coordinate can be derived using the global coordinate measurements and the reference quaternion as
z k , i L = R G L ( q k ) ( z k , i G x k c ) .
The azimuth and elevation of the measurements are
θ k L = arctan y L x L ,
ϕ k L = arctan z L ( x L ) 2 + ( y L ) 2 ,
γ k L = ( θ k L , ϕ k L ) .
Substituting Equation (21) into Equation (61) gives
z k , i G = x k c + p k , i H f ( γ k L ) f k + e k f + e ¯ k ,
which can be further simplified to obtain the pseudo-measurement function
h ( x k , z k , i G ) + e k , i = 0 ,   e k , i ~ N ( 0 , R k , i ) ,
where
h ( x k , z k , i G ) = x k c z k , i G + p k , i H f ( γ k L ) f k ,
e k , i = p k , i e k f + e ¯ k ,
R k , i = p k , i R k , i f p k , i T + R ¯ k , i .

4. Gaussian Process Gaussian Mixture PHD Filter

4.1. Prediction Process

We use G A M ( ) to denote the Gamma distribution. The prediction process of the extended target PHD filter is given as [30]
D k + 1 | k ( ξ ) = D k + 1 | k s ( ξ ) + D k b ( ξ ) ,
D k + 1 | k s ( ξ ) = j = 1 J k 1 | k 1 w k + 1 | k ( j ) G A M ( γ k + 1 | k ; α k + 1 | k ( j ) , β k + 1 | k ( j ) ) N ( x k + 1 | k ; m k + 1 | k ( j ) , P k + 1 | k ( j ) ) ,
where D k + 1 | k s ( ξ ) is the survival PHD, and D k + 1 | k b ( ξ ) is the birth PHD representing new targets. The survival PHD is given as
w k + 1 | k ( j ) = p s w k | k ( j ) ,
α k + 1 | k ( j ) = α k | k ( j ) η k , β k + 1 | k ( j ) = β k | k ( j ) η k ,
m k + 1 | k ( j ) = F k m k | k ( j ) ,
P k + 1 | k ( j ) = F k P k + 1 | k ( j ) F k T + Q k .
The birth PHD is approximated using a Gaussian mixture model
D k + 1 | k b ( ξ ) = j = 1 J b , k w b , k ( j ) G A M ( γ b , k ; α b , k ( j ) , β b , k ( j ) ) N ( x b , k ; m b , k ( j ) , P b , k ( j ) ) .

4.2. Update Process

Without considering the target spawn, the update process of the PHD filter is
D k + 1 | k + 1 ( ξ ) = D k + 1 | k + 1 ND ( ξ ) + D k + 1 | k + 1 D ( ξ ) ,
where D k + 1 | k + 1 ND ( ξ ) is the miss detection PHD given as
D k + 1 | k + 1 ND ( ξ ) = j = 1 J k + 1 | k w k + 1 | k + 1 ND , ( j ) G A M ( γ k + 1 | k + 1 ; α k + 1 | k + 1 ND , ( j ) , β k + 1 | k + 1 ND , ( j ) ) × N ( x k + 1 | k + 1 ; m k + 1 | k + 1 ND , ( j ) , P k + 1 | k + 1 ND , ( j ) )
where
w k + 1 | k + 1 ND , ( j ) = ( 1 p D ) w k + 1 | k ( j ) ,
α k + 1 | k + 1 ND , ( j ) = α k + 1 | k ( j ) , β k + 1 | k + 1 ND , ( j ) = β k + 1 | k ( j ) ,
m k + 1 | k + 1 ND , ( j ) = m k + 1 | k ( j ) ,
P k + 1 | k + 1 ND , ( j ) = P k + 1 | k ( j ) .
The updated target PHD is D k + 1 | k + 1 D ( ξ ) , which is given as
D k + 1 | k + 1 D ( ξ ) = P Z k W P j = 1 J k + 1 | k w k + 1 | k + 1 ( j , W ) G A M ( γ k + 1 | k + 1 ; α k + 1 | k + 1 ( j , W ) , β k + 1 | k + 1 ( j , W ) ) × N ( x k + 1 | k + 1 ; m k + 1 | k + 1 ( j , W ) , P k + 1 | k + 1 ( j , W ) ) .
The partition method in [30] is adopted in this article to reduce the computational load. The update of the Gamma distribution is
α k + 1 | k + 1 ( j , W ) = α k | k + 1 ( j ) + W ,
β k + 1 | k + 1 ( j , W ) = β k | k + 1 ( j ) + 1 .
The GP regression measurement model derived in the previous section is used to update the Gaussian mixture component. The pseudo-measurement model Equation (67) is a non-linear system. Therefore, the EKF is used to update the extended target states, which can be given as follows
m k + 1 | k + 1 ( j , W ) = m k | k + 1 ( j ) + K k + 1 0 h m k | k + 1 ( j ) , z k ( W ) ,
P k + 1 | k + 1 ( j , W ) = P k | k + 1 ( j ) K k + 1 H k P k | k + 1 ( j ) .
where h ( m k | k + 1 ( j ) , z k ( W ) ) is given in (68) and the Kalman gain K k + 1 is
K k + 1 = P k | k + 1 ( j ) H k T S k 1 ,
S k = H k P k | k + 1 ( j ) H k T + R k ,
H k = d d x k h k ( m k | k + 1 ( j ) , z k ( W ) ) x k = m ^ k | k + 1 ( j ) .
The weight of the gamma Gaussian mixture component is
w k + 1 | k + 1 ( j , W ) = w p p D k ( j , W ) w k | k 1 ( j ) d W β FA , k | W | ,
L k ( j , W ) = z k W ϕ ( z k | x k ) = N W ; h m k | k + 1 ( j ) , z k ( W ) , R k ,
where ϕ z k ( ) is the single extended target measurement likelihood.

5. Simulation Results and Discussion

In this section, a simulation experiment is conducted to compare the performance of the GP-PHD filter proposed in the article with the traditional GGIW-PHD filter [33]. A three-dimensional area with size [ 30   m , 30   m , 5   m ] × [ 30   m ,   3 0   m ,   5   m ] is considered. We assume that two extended targets move in the x–y–z plane. Target 1 (T1) is a cylinder with r = 1   m ,   h = 2   m and target 2 (T2) is a cube with w = 2   m ,   l = 4   m , and h = 2   m . We denote the target states as x = [ m x , m y , m z , v x , v y , v z ] T . The initial states of T1 is x 1 = [ 0 , 0 , 0 , 1 , 0 , 0 ] T and the initial states of T2 is x 2 = [ 10 , 0 , 1 , 1 , 0 , 0 ] T . T1 moves at constant speed v = 1   m / s along the y axis. T2 moves at constant speed v = 1.5   m / s in the first 10 s, and then makes a constant turn at speed ω = π 20   rad / s for 10 s. After the constant turn, T2 moves along the y axis at a constant speed. Both targets’ moving duration is 30   s and the sampling time is 0.1   s . The experiment is a simulation of a typical traffic scenario, where T1 corresponds to the pedestrian and T2 corresponds to a practical medium size vehicle. Both linear and non-linear motion models of extended targets are considered. Although we use a linear motion model, the filter is expected to be robust enough to handle such model mismatches, which are common in most tracking applications.
We assume that all the measurements come from the surface of the target and that the number of the measurements is Poisson-distributed. The Poisson rate of the measurements is λ = 20 and the covariance matrix of the measurement noise is R ¯ = 0.1 2 Ι 3   m 2 . The Poisson rate of the number of clutter is λ c = 2 , and clutter is uniformly distributed in the surveillance area. Figure 2 shows clutter and measurements collected in the surveillance area. All experiments in the following section are carried out with different realizations of the measurement noise, the measurement source, and the clutter. The result numbers of the experiments are the averages of Monte Carlo (MC) runs.
The CV model is considered for both the GGIW-PHD filter and the GP-PHD filter. The process noise of the dynamic model is σ x = 0.1 . For the GP-PHD filter, 642 basis points equally located in [ 0 , 2 π ] are used and the hyperparameters of the GP model are set to μ = 0 ,   σ f = 1 ,   σ r = 0.2 ,   l = π / 3 . The parameters for the GGIW-PHD filter are given in Table 1.
The number of the gamma Gaussian mixture component is no more than M m a x = 100 and the number of Gaussian components per target is no more than J m a x = 100 . The tracking results of the two filters are shown in Figure 3. Figure 3 shows that both filters are able to deal with the non-linear movements of the targets, while the GP-PHD filter is more capable of handling the motion model mismatches and the estimated trajectory is closer to the truth.
We use the OSPA distance to evaluate the performance of the filter, which takes into account both errors made in estimating target states and the errors in estimating target cardinality. The order and the penalty factor are set as c = 15 ,   p = 1 . The average result is obtained through 100 MC runs, which is shown in Figure 4.
The estimation of the shape is illustrated in Figure 5. The GP-PHD filter shows superior performance compared to the GGIW-PHD filter in terms of shape estimation. The GP-PHD filter uses a radial function to describe the shape of the target, which is suitable for arbitrary shape extended targets. However, the GGIW-PHD filter is limited to ellipsoidal shape targets. Therefore, both filters show satisfactory performance in estimating the shape of T1 in Figure 5a. However, the GGIW-PHD filter shows a poor performance in estimating the shape of T2, especially when the target is making a turn. The GP-PHD filter takes into account the orientation of the target. Figure 5b shows that the GP-PHD filter is able to estimate the shape of the target, even when the target is making a turn.
The performance of the estimation of the shape is evaluated using the IoU, which is defined as
IoU ( S true , S ^ ) = volume ( S true S ^ ) volume ( S true S ^ ) ,
where S true is the true target shape and S ^ is the estimated target shape. The IoU value not only accounts for the quality of the estimation of the kinematics states but also the extent. The average results of 100 MC runs are shown in Figure 6. Overall, the GP-PHD filter estimates the shape of the extended targets better than the GGIW-PHD filter. The performance of the GGIW-PHD filter is not stable, particularly when the target is making a turn.
We further examine the measurement rates’ impact on the proposed filter and the GGIW-PHD filter. The measurement rates are set as λ = 10 , 30 . The experiments are repeated for 50 MC runs, respectively, to obtain the average results shown in Figure 7 and Figure 8. It is reported that the GP-PHD filter shows good performance at both low and high measurement rates, while the measurement rate has a great impact on the performance of the GGIW-PHD filter. The GP-PHD filter is more robust compared to the GGIW-PHD filter. The main reason for the impact of the measurement rate on the GGIW-PHD filter is that there is less information about the shape that can be exploited by the GGIW-PHD filter when the measurement rate is low. However, the covariance function of the GP model assumes that the radial function has a period of π . This assumption suggests that the target is symmetric, which is the usual case in tracking. Therefore, under low-measurement-rate conditions, the GP-PHD filter can obtain more information about the target and a more accurate estimation of the shape can be expected.
All the simulations are conducted in MATLAB 2022b on a computer with an AMD Ryzen 5 3600X 6-Core processor and 16 GB RAM. The average computation time for the GP-PHD filter and the GGIW-PHD filter are 459 s and 5.11 s for 300 updates. Both the GGIW-PHD filter and the GP-PHD filter are recursively updated by the standard EKF, and therefore the computational load is mainly determined by the dimension of the state vector. The state dimension of the GP-PHD filter is dim ( x k ) = 654 , which is much greater than that of the GGIW-PHD filter. Further improvements, such as projecting the contour of the 3D target into 2D space to reduce the computational burden, can be considered.

6. Conclusions

In this paper, we propose a new approach for tracking multiple extended targets in three-dimensional space. The key contribution of our method is the utilization of a 3D Gaussian process radial function to describe the shape of the extended targets. Furthermore, we also present a state estimation method that allows for the recursive estimation of the target’s shape. A pseudo-measurement equation is derived, and a Gaussian mixture implementation of the filter is presented in the article. The proposed filter is capable of tracking the target and estimating the shape of the target simultaneously in the presence of clutter. The method is flexible enough to express a variety of different shapes. We evaluate the performance of the proposed filter and the traditional GGIW-PHD filter in the simulation experiment using the OSPA distance and the IoU value. The results demonstrate that our proposed filter outperforms the traditional GGIW-PHD filter in both kinematic states and shape estimation. Furthermore, we compare the performance of the filters under different-measurement-rate conditions. The proposed filter is robust under different-measurement-rate conditions, while the GGIW-PHD filter suffers under low-measurement-rate conditions. The algorithm provides an effective solution to the ETT problem in 3D space, which can be further used in scenarios such as autonomous driving.

Author Contributions

Conceptualization, Z.Y. and J.S.; methodology, Z.Y. and X.L.; software, Z.Y. and X.L.; validation, J.S., X.Y. and T.S.; formal analysis, Z.Y. and X.L.; investigation, X.Y.; resources, J.S.; data curation, Z.Y. and X.L.; writing—original draft preparation, Z.Y.; writing—review and editing, J.S., X.Y. and T.S.; visualization, T.S.; supervision, J.S.; project administration, J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, Grant No. 62131001 and 62171029.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Karl, G.; Marcus, B.; Stephan, R. Extended Object Tracking: Introduction, Overview, and Applications. J. Adv. Inf. Fusion 2017, 2, 139–174. [Google Scholar]
  2. Mahler, R. PHD Filters for Nonstandard Targets, I: Extended Targets. In Proceedings of the 12th International Conference on Information Fusion, Seattle, WA, USA, 6–9 July 2009; pp. 915–921. [Google Scholar]
  3. Michaelis, M.; Berthold, P.; Luettel, T.; Wuensche, H.J. Extended Target Tracking with a Particle Filter Using State Dependent Target Measurement Models. In Proceedings of the 25th International Conference on Information Fusion, Linköping, Sweden, 4–7 July 2018; pp. 259–266. [Google Scholar]
  4. De Freitas, A.; Mihaylova, L.; Gning, A.; Schikora, M.; Ulmke, M.; Angelova, D.; Koch, W. A Box Particle Filter Method for Tracking Multiple Extended Objects. IEEE Trans. Aerosp. Electron. Syst. 2019, 55, 1640–1655. [Google Scholar] [CrossRef] [Green Version]
  5. Tuncer, B.; Orguner, U.; Ozkan, E. Multi-Ellipsoidal Extended Target Tracking with Variational Bayes Inference. IEEE Trans. Signal Process. 2022, 70, 3921–3934. [Google Scholar] [CrossRef]
  6. Li, P.; Wang, W.; Qiu, J.; You, C.; Shu, Z. Robust Generalized Labeled Multi-Bernoulli Filter and Smoother for Multiple Target Tracking Using Variational Bayesian. KSII Trans. Internet Inf. Syst. 2022, 16, 908–928. [Google Scholar]
  7. Yasir, M.; Zhan, L.; Liu, S.; Wan, J.; Hossain, M.S.; Isiacik Colak, A.T.; Liu, M.; Islam, Q.U.; Raza Mehdi, S.; Yang, Q. Instance Segmentation Ship Detection Based on Improved Yolov7 Using Complex Background SAR Images. Front. Mar. Sci. 2023, 10, 1113669. [Google Scholar] [CrossRef]
  8. Meyer, F.; Win, M.Z. Scalable data association for extended object tracking. IEEE Trans. Signal Inf. Process. Netw. 2020, 6, 491–507. [Google Scholar] [CrossRef]
  9. 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]
  10. Feldmann, M.; Fränken, D.; Koch, W. Tracking of Extended Objects and Group Targets Using Random Matrices. IEEE Trans. Signal Process. 2011, 59, 1409–1420. [Google Scholar] [CrossRef]
  11. Baum, M.; Hanebeck, U.D. Random Hypersurface Models for Extended Object Tracking. In Proceedings of the 2009 IEEE International Symposium on Signal Processing and Information Technology (ISSPIT), Ajman, United Arab Emirates, 14–17 December 2009; pp. 178–183. [Google Scholar]
  12. Sun, L.; Zhang, J.; Yu, H.; Fu, Z.; He, Z. Tracking of Maneuvering Extended Target Using Modified Variable Structure Multiple-Model Based on Adaptive Grid Best Model Augmentation. Remote Sens. 2022, 14, 1613. [Google Scholar] [CrossRef]
  13. Li, Y.; Wei, P.; You, M.; Wei, Y.; Zhang, H. Joint Detection, Tracking, and Classification of Multiple Extended Objects Based on the JDTC-PMBM-GGIW Filter. Remote Sens. 2023, 15, 887. [Google Scholar] [CrossRef]
  14. Lan, J.; Li, X.R. Tracking of Extended Object or Target Group Using Random Matrix: New Model and Approach. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 2973–2989. [Google Scholar] [CrossRef]
  15. 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]
  16. Sun, L.; Yu, H.; Lan, J.; Fu, Z.; He, Z.; Pu, J. Tracking of Multiple Maneuvering Random Hypersurface Extended Objects Using High Resolution Sensors. Remote Sens. 2021, 13, 2963. [Google Scholar] [CrossRef]
  17. Yang, J.-L.; Li, P.; Ge, H.-W. Extended Target Shape Estimation by Fitting B-Spline Curve. J. Appl. Math. 2014, 2014, 741892. [Google Scholar] [CrossRef] [Green Version]
  18. Wang, Y.; Chen, X.; Gong, C.; Rao, P. Non-Ellipsoidal Infrared Group/Extended Target Tracking Based on Poisson Multi-Bernoulli Mixture Filter and B-Spline. Remote Sens. 2023, 15, 606. [Google Scholar] [CrossRef]
  19. Zea, A.; Faion, F.; Hanebeck, U.D. Tracking Elongated Extended Objects Using Splines. In Proceedings of the 19th International Conference on Information Fusion, Heidelberg, Germany, 5–8 July 2016; pp. 612–619. [Google Scholar]
  20. Rasmussen, C.E.; Williams, C.K.I. Gaussian Processes for Machine Learning; MIT Press: Cambridge, MA, USA, 2006. [Google Scholar]
  21. Wahlström, N.; Özkan, E. Extended Target Tracking Using Gaussian Processes. IEEE Trans. Signal Process. 2015, 63, 4165–4178. [Google Scholar] [CrossRef] [Green Version]
  22. Aftab, W.; Hostettler, R.; De Freitas, A.; Arvaneh, M.; Mihaylova, L. Spatio-Temporal Gaussian Process Models for Extended and Group Object Tracking with Irregular Shapes. IEEE Trans. Veh. Technol. 2019, 68, 2137–2151. [Google Scholar] [CrossRef]
  23. Akbari, B.; Zhu, H. Tracking dependent extended targets using multi-output spatiotemporal Gaussian processes. IEEE Trans. Intell. Transp. Syst. 2022, 10, 18301–18314. [Google Scholar] [CrossRef]
  24. Mahler, R. Statistical Multisource-Multitarget Information Fusion; Artech House, Inc.: Norwood, MA, USA, 2007. [Google Scholar]
  25. Liu, Y.; Xu, K.D. Millimeter-Wave Bandpass Filters Using On-Chip Dual-Mode Resonators in 0.13-μm SiGe BiCMOS Technology. IEEE Trans. Microw. Theory Tech. 2023; Early Access. [Google Scholar] [CrossRef]
  26. Xu, K.D.; Liu, Y. Millimeter-Wave On-Chip Bandpass Filter Using Complementary-Broadside-Coupled Structure. IEEE Trans. Circuits Syst. II Exp. Briefs 2023, 1. [Google Scholar] [CrossRef]
  27. Feng, Y.; Zhang, B.; Liu, Y.; Niu, Z.; Fan, Y.; Chen, X. A D-band manifold triplexer with high isolation utilizing novel waveguide dual-mode filters. IEEE Trans. Terahertz Sci. Technol. 2022, 12, 678–681. [Google Scholar] [CrossRef]
  28. Xu, K.D.; Guo, Y.J.; Liu, Y.; Deng, X.; Chen, Q.; Ma, Z. 60-GHz compact dual-mode on-chip bandpass filter using GaAs technology. IEEE Electron. Device Lett. 2021, 42, 1120–1123. [Google Scholar] [CrossRef]
  29. Mahler, R. Multitarget Bayes Filtering via First-Order Multitarget Moments. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 1152–1178. [Google Scholar] [CrossRef]
  30. Granström, K.; Lundquist, C.; Orguner, O. Extended Target Tracking Using a Gaussian-Mixture PHD Filter. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 3268–3286. [Google Scholar] [CrossRef] [Green Version]
  31. Granström, 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] [Green Version]
  32. Granström, K.; Orguner, U. On Spawning and Combination of Extended/Group Targets Modeled with Random Matrices. IEEE Trans. Signal Process. 2013, 61, 678–692. [Google Scholar] [CrossRef] [Green Version]
  33. Granström, K.; Natale, A.; Braca, P.; Ludeno, G.; Serafino, F. Gamma Gaussian Inverse Wishart Probability Hypothesis Density for Extended Target Tracking Using X-Band Marine Radar Data. IEEE Trans. Geosci. Remote Sens. 2015, 53, 6617–6631. [Google Scholar] [CrossRef]
  34. Li, B.; Mu, C.; Bai, Y.; Bi, J.; Wang, L. Ellipse Fitting Based Approach for Extended Object Tracking. Math. Probl. Eng. 2014, 2014, 632815. [Google Scholar] [CrossRef]
  35. Faion, F.; Baum, M.; Hanebeck, U.D. Tracking 3D Shapes in Noisy Point Clouds with Random Hypersurface Models. In Proceedings of the 15th International Conference on Information Fusion, Singapore, 9–12 July 2012; pp. 2230–2235. [Google Scholar]
  36. Kumru, M.; Özkan, E. 3D Extended Object Tracking Using Recursive Gaussian Processes. In Proceedings of the 21st International Conference on Information Fusion, Cambridge, UK, 10–13 July 2018; pp. 259–266. [Google Scholar]
  37. Kumru, M.; Ozkan, E. Three-Dimensional Extended Object Tracking and Shape Learning Using Gaussian Processes. IEEE Trans. Aerosp. Electron. Syst. 2021, 57, 2795–2814. [Google Scholar] [CrossRef]
  38. 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] [Green Version]
  39. Huber, M.F. Recursive Gaussian Process Regression. In Proceedings of the 2013 IEEE International Conference on Acoustics, Speech and Signal Processing, Vancouver, BC, Canada, 26–31 May 2013; pp. 3362–3366. [Google Scholar]
  40. Maeder, U.; Morari, M. Attitude Estimation for Vehicles with Partial Inertial Measurement. IEEE Trans. Veh. Technol. 2011, 60, 1496–1504. [Google Scholar] [CrossRef]
  41. Markley, F.L. Attitude Estimation or Quaternion Estimation? J. Astronaut. Sci. 2004, 52, 221–238. [Google Scholar] [CrossRef]
Figure 1. The global coordinate and the local coordinate.
Figure 1. The global coordinate and the local coordinate.
Remotesensing 15 03224 g001
Figure 2. The measurements and clutter collected in the surveillance area.
Figure 2. The measurements and clutter collected in the surveillance area.
Remotesensing 15 03224 g002
Figure 3. Tracking results of the GGIW−PHD filter and the GP−PHD filter.
Figure 3. Tracking results of the GGIW−PHD filter and the GP−PHD filter.
Remotesensing 15 03224 g003
Figure 4. The OSPA distance error for the GGIW-PHD filter and the GP-PHD filter with 100 Monte Carlo runs.
Figure 4. The OSPA distance error for the GGIW-PHD filter and the GP-PHD filter with 100 Monte Carlo runs.
Remotesensing 15 03224 g004
Figure 5. The results of shape estimation of the GGIW-PHD filter and the GP-PHD filter. (a) shape estimation of T1; (b) shape estimation of T2.
Figure 5. The results of shape estimation of the GGIW-PHD filter and the GP-PHD filter. (a) shape estimation of T1; (b) shape estimation of T2.
Remotesensing 15 03224 g005
Figure 6. The IoU plots averaged over 100 MC runs.
Figure 6. The IoU plots averaged over 100 MC runs.
Remotesensing 15 03224 g006
Figure 7. Average OSPA distance for 50 MC runs of the GGIW-PHD filter and the GP-PHD filter in different measurement rates. (a) λ = 10 ; (b) λ = 30 .
Figure 7. Average OSPA distance for 50 MC runs of the GGIW-PHD filter and the GP-PHD filter in different measurement rates. (a) λ = 10 ; (b) λ = 30 .
Remotesensing 15 03224 g007
Figure 8. Average IoU value for 50 MC runs of the GGIW-PHD filter and the GP-PHD filter in different measurement rates. (a) λ = 10 ; (b) λ = 30 .
Figure 8. Average IoU value for 50 MC runs of the GGIW-PHD filter and the GP-PHD filter in different measurement rates. (a) λ = 10 ; (b) λ = 30 .
Remotesensing 15 03224 g008
Table 1. Parameters of the GGIW-PHD filter.
Table 1. Parameters of the GGIW-PHD filter.
ParameterValue
Forgetting factor η k 1.05
Detection probability p D 0.95
Survival probability p s 0.99
Birth weight ω k ( b ) 0.11
Birth rate a k ( b ) ,   b k ( b ) 100, 5
Birth extension v k ( b ) ,   V k ( b ) 10, 10
Pruning threshold T 10 3
Merging threshold U25
Temporal decay τ120
Scaling parameter ρ1/4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yang, Z.; Li, X.; Yao, X.; Sun, J.; Shan, T. Gaussian Process Gaussian Mixture PHD Filter for 3D Multiple Extended Target Tracking. Remote Sens. 2023, 15, 3224. https://doi.org/10.3390/rs15133224

AMA Style

Yang Z, Li X, Yao X, Sun J, Shan T. Gaussian Process Gaussian Mixture PHD Filter for 3D Multiple Extended Target Tracking. Remote Sensing. 2023; 15(13):3224. https://doi.org/10.3390/rs15133224

Chicago/Turabian Style

Yang, Zhiyuan, Xiangqian Li, Xianxun Yao, Jinping Sun, and Tao Shan. 2023. "Gaussian Process Gaussian Mixture PHD Filter for 3D Multiple Extended Target Tracking" Remote Sensing 15, no. 13: 3224. https://doi.org/10.3390/rs15133224

APA Style

Yang, Z., Li, X., Yao, X., Sun, J., & Shan, T. (2023). Gaussian Process Gaussian Mixture PHD Filter for 3D Multiple Extended Target Tracking. Remote Sensing, 15(13), 3224. https://doi.org/10.3390/rs15133224

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