# Tracking Split Group withδ-Generalized Labeled Multi-Bernoulli Filter

Authors: Linhai Gan; Gang Wang
Journal: Journal of Sensors
 (2019)

---

## Abstract

As target splitting is not considered in the initial development ofδ-generalized labeled multi-Bernoulli (δ-GLMB) filter, the scenarios where the new targets appearing conditioned on the preexisting one are not readily addressed by this filter. In view of this, we model the group target as gamma Gaussian inverse Wishart (GGIW) distribution and derive a δ-GLMB filter based on the group splitting model, in which the target splitting event is investigated. Two simplifications of the approach are presented to improve the computing efficiency, where with splitting detection, we need not to predict the splitting events of all the GGIW components in every iteration. With component combination applied in adaptive birth, a redundant modeling for a newborn target or preexisting target could be avoided. Moreover, a method for labeling performance evaluation of the algorithm is provided. Simulations demonstrate the effectiveness of the proposed approach.


---

## Body

## 1. Introduction

The traditional multitarget tracking algorithm is mainly based on a standard measurement (point target) model, which assumes that one target could produce at most one measurement at a given time, and each measurement originates from at most one target [1]. However, there are many scenes; for example, a group of targets may produce a single measurement, or a single target may produce several measurements in reality, which should be handled by nonstandard measurement models. Groups are structured objects and formations of entities moving in a coordinated manner, whose number varies over time. The groups can split; in combination, they can be relatively near each other or move independently of each other [2]. As for dense group target tracking, while the number of targets in a group is large and the distribution is dense, if a point target model is used, the following problems will be encountered: ① it is hard to distinguish the targets in dense distribution, and the shelter or locomotion between targets will result in a frequency appearance or disappearance of the target, which makes it difficult to establish a stable track for each target. ② The computation resource and sensor resource will be excessively occupied in the target association and tracking process due to the large number of targets. So, a specified target model is needed for group target tracking. Though essentially different, the measurement appearance of the extended object and dense group target is similar, and their identical tracking models could be adopted [2].Tracking methods for the group with dense distribution mainly contain Poisson likelihood method [3–5], random matrix (RM) method [6–8], random hypersurface method [9–12], random finite set (RFS) method [13–15], etc. The random matrix approach initiated by Koch [16] jointly estimated the centroid state and extension state of the group with the assumption that the centroid state follows Gaussian distribution and describes the extension shape of the group with a random matrix following inverse Wishart (IW) distribution. The Gaussian inverse Wishart probability hyperthesis density (GIW-PHD) approach is proposed with the combination of the RM method and RFS theory [17]. Splitting and combination are important forms of group target motion. Group splitting and combination are modeled in [18], where an approach named GGIW- PHD is proposed, with the assumption that the measurement rate of the target follows Gamma distribution. The cardinality PHD (CPHD) approach is introduced to group target tracking in [19, 20], which could estimate the cardinality of the group target, thereby improving the performance of the approach. PHD filter and CPHD filter avoid data combination in multitarget tracking and improve the computation efficiency of the Bayes multitarget filter vastly, but could not explicitly accommodate the estimation of the target trajectory. Moreover, the PHD filter is known to produce highly uncertain estimates of target number and the CPHD filter is limited by the so-called “spooky” effect [21–24].The GLMB filter [25–29], which addresses the limitations of PHD filter and CPHD filter, could estimate the state and the trajectory of the target simultaneously. With the GLMB filter introduced in group target tracking, the GGIW-GLMB filter is proposed in [1] and thereby a better performance in accuracy is gained compared to PHD filter and CPHD filter; however, the group splitting is not addressed.This paper introduces theδ-GLMB filter for group target tracking; although it is able to estimate the track of the group, the split of the group is not readily addressed. In view of this, we consider tracking the group target based on a splitting model and derive a GGIW-δ-GLMB filter for multiple split group target tracking. Besides, simplifications concerning on splitting prediction and the target adaptive birth are presented to make our approach more feasible computationally. Moreover, the track labeling OSPA error is proposed for labeling performance evaluation.The paper is organized as follows: In Section2, the background is provided on the GGIW-GLMB filter. In Section 3, the derivation of GGIW-δ-GLMB for the multiple splitting group is developed. Two simplifications to improve computation efficiency are presented in Section 4. Section 5 presents simulation results, and Section 6 concludes the work.


## 2. Background

### 2.1. The Dynamic and Measurement Model of Group Target

Assume that the centroid state of theith group target at time k is xki, the dynamic model is given by
(1)xki=Fkk−1xk−1i+vki,where Fkk−1=F~kk−1⊗Id with F~kk−1 being the dynamic matrix in one-dimensional space. ⊗ denotes the Kronecker product, Id∈Rd×d denotes the identity matrix. vki∼Nwkj;0,Q~k∣k−1⊗Xki is the independent Gaussian process noise with Q~k∣k−1 being the covariance matrix in the one-dimensional model, and Qk∣k−1=Q~k∣k−1⊗Xki. Xki is a symmetric positive definite matrix describing the extension state of the group target.Given the measurement set at timek to be Zk=zkjj=1NZ,k, the measurement model of the group target is given by
(2)zkj=Hkxki+wkj.NZ,k is the number of measurement at time k. Hk=H~k⊗Id, specially, while d=2, H~k=1,0. wkj∼Nwkj;0,Xki+Rk is the independent Gaussian white noise with Rk being the covariance.


### 2.2. GLMB Filter for Multigroup Target

#### 2.2.1. Labeled RFS

Suppose thatX and L are the state space and discrete label space of the multitarget, respectively, let ξ∈X, l∈L, then X=ξ,li,i=1,2,3,…,X is an RFS on X×L. X denotes the number of elements in the RFS. Define the distinct label indicator function as
(3)ΔX=1,LX=X,0,LX≠X,where LX denotes the set of unique labels in X. The label consists of the time of target birth and a serial number of the target, which are distinct with each other in each realization, which means ΔX=1 is always satisfied.Define the Kronecker delta function and the set inclusion function as(4)δYX=1,X=Y,0,X≠Y,1YX=1,X⊆Y,0,X⊄Y,where X,Y may be scalars, vectors, or sets.The probability density of GLMB RFS is distributed as follows:(5)πX=ΔX∑c∈CwcLXpc·X,where C is the discrete index set, h·X=∏ξ∈Xhξ, and
(6)∑L∈L∑c∈CwcL=1,∫ξ∈Xpcξ,ldξ=1.


#### 2.2.2. The Labeled Likelihood of Group Target

The measurement likelihood function of the multigroup target is given by [1]
(7)gZX=gCZ∑i=1X+1∑UZ∈PiZθ∈ΘUZψUZ·;θX,where Z is a finite measurement set, PiZ denotes the set of all the partitions dividing Z into i groups, and UZ∈PiZ is a particular partition of Z. θ: LX→0,1,⋯,UZ denotes the association map of the target for measurement, while θl=θl′>0, then l=l′. ΘUZ is the space of θ, and UθlZ denotes the element of UZ corresponding to label l.
(8)ψUZξ,l;θ=pDξ,lg~UθlZξ,lkUθlZ,θl>0,qDξ,l,θl=0.


#### 2.2.3. The Prediction and Update of GGIW-GLMB

At timek, let γk∼GAℳγk;αk,βk be the measurement rate which follows Poisson distribution, xk∼Nxk;mk,P~k⊗Xk be the centroid state, and Xk∼ℐWXk;vk,Vk be the extension state. Then, the probability density function (PDF) of the group target state is given by
(9)pξk=pγkpxkpXk=GAℳγk;αk,βkNxk;mk,P~k⊗XkℐWXk;vk,Vk=GGℐWξk;ζk,where ξk≜γk,xk,Xk denotes the state to be estimated, ζk≜αk,βk,mk,P~k,vk,Vk denotes the parameter of the GGIW component, P~ denotes the state covariance of a single dimension, and Pk=P~k⊗Xk denotes the state covariance of all dimensions. Substitute (9) into (5), the distribution of GLMB RFS of the multigroup target could be obtained subsequently.Let the subscript “+” be the index of the next time; e.g.,s+ might denote the probability density, label, label space, etc., according to the meaning of s. Then, the expression of prediction and update of GGIW-GLMB is given by [1]
(1)
Prediction(10)π+X+=ΔX+∑c∈Cw+cLX+p+c·χ+,where w+c· and p+c· denote the weight of the predicted GLMB and the predicted probability density of the group target state, respectively.
(2)
Update(11)πX∣Z=ΔX∑c∈C∑i=1X+1∑UZ∈PiZθ∈ΘUZwUZc,θLXpc,θ·∣UZX,where wUZc,θ· and pc,θ·∣UZ denote the weight of the updated GLMB component and updated probability density of the group target sate, respectively.


## 2.1. The Dynamic and Measurement Model of Group Target

Assume that the centroid state of theith group target at time k is xki, the dynamic model is given by
(1)xki=Fkk−1xk−1i+vki,where Fkk−1=F~kk−1⊗Id with F~kk−1 being the dynamic matrix in one-dimensional space. ⊗ denotes the Kronecker product, Id∈Rd×d denotes the identity matrix. vki∼Nwkj;0,Q~k∣k−1⊗Xki is the independent Gaussian process noise with Q~k∣k−1 being the covariance matrix in the one-dimensional model, and Qk∣k−1=Q~k∣k−1⊗Xki. Xki is a symmetric positive definite matrix describing the extension state of the group target.Given the measurement set at timek to be Zk=zkjj=1NZ,k, the measurement model of the group target is given by
(2)zkj=Hkxki+wkj.NZ,k is the number of measurement at time k. Hk=H~k⊗Id, specially, while d=2, H~k=1,0. wkj∼Nwkj;0,Xki+Rk is the independent Gaussian white noise with Rk being the covariance.


## 2.2. GLMB Filter for Multigroup Target

### 2.2.1. Labeled RFS

Suppose thatX and L are the state space and discrete label space of the multitarget, respectively, let ξ∈X, l∈L, then X=ξ,li,i=1,2,3,…,X is an RFS on X×L. X denotes the number of elements in the RFS. Define the distinct label indicator function as
(3)ΔX=1,LX=X,0,LX≠X,where LX denotes the set of unique labels in X. The label consists of the time of target birth and a serial number of the target, which are distinct with each other in each realization, which means ΔX=1 is always satisfied.Define the Kronecker delta function and the set inclusion function as(4)δYX=1,X=Y,0,X≠Y,1YX=1,X⊆Y,0,X⊄Y,where X,Y may be scalars, vectors, or sets.The probability density of GLMB RFS is distributed as follows:(5)πX=ΔX∑c∈CwcLXpc·X,where C is the discrete index set, h·X=∏ξ∈Xhξ, and
(6)∑L∈L∑c∈CwcL=1,∫ξ∈Xpcξ,ldξ=1.


### 2.2.2. The Labeled Likelihood of Group Target

The measurement likelihood function of the multigroup target is given by [1]
(7)gZX=gCZ∑i=1X+1∑UZ∈PiZθ∈ΘUZψUZ·;θX,where Z is a finite measurement set, PiZ denotes the set of all the partitions dividing Z into i groups, and UZ∈PiZ is a particular partition of Z. θ: LX→0,1,⋯,UZ denotes the association map of the target for measurement, while θl=θl′>0, then l=l′. ΘUZ is the space of θ, and UθlZ denotes the element of UZ corresponding to label l.
(8)ψUZξ,l;θ=pDξ,lg~UθlZξ,lkUθlZ,θl>0,qDξ,l,θl=0.


### 2.2.3. The Prediction and Update of GGIW-GLMB

At timek, let γk∼GAℳγk;αk,βk be the measurement rate which follows Poisson distribution, xk∼Nxk;mk,P~k⊗Xk be the centroid state, and Xk∼ℐWXk;vk,Vk be the extension state. Then, the probability density function (PDF) of the group target state is given by
(9)pξk=pγkpxkpXk=GAℳγk;αk,βkNxk;mk,P~k⊗XkℐWXk;vk,Vk=GGℐWξk;ζk,where ξk≜γk,xk,Xk denotes the state to be estimated, ζk≜αk,βk,mk,P~k,vk,Vk denotes the parameter of the GGIW component, P~ denotes the state covariance of a single dimension, and Pk=P~k⊗Xk denotes the state covariance of all dimensions. Substitute (9) into (5), the distribution of GLMB RFS of the multigroup target could be obtained subsequently.Let the subscript “+” be the index of the next time; e.g.,s+ might denote the probability density, label, label space, etc., according to the meaning of s. Then, the expression of prediction and update of GGIW-GLMB is given by [1]
(1)
Prediction(10)π+X+=ΔX+∑c∈Cw+cLX+p+c·χ+,where w+c· and p+c· denote the weight of the predicted GLMB and the predicted probability density of the group target state, respectively.
(2)
Update(11)πX∣Z=ΔX∑c∈C∑i=1X+1∑UZ∈PiZθ∈ΘUZwUZc,θLXpc,θ·∣UZX,where wUZc,θ· and pc,θ·∣UZ denote the weight of the updated GLMB component and updated probability density of the group target sate, respectively.


## 2.2.1. Labeled RFS

Suppose thatX and L are the state space and discrete label space of the multitarget, respectively, let ξ∈X, l∈L, then X=ξ,li,i=1,2,3,…,X is an RFS on X×L. X denotes the number of elements in the RFS. Define the distinct label indicator function as
(3)ΔX=1,LX=X,0,LX≠X,where LX denotes the set of unique labels in X. The label consists of the time of target birth and a serial number of the target, which are distinct with each other in each realization, which means ΔX=1 is always satisfied.Define the Kronecker delta function and the set inclusion function as(4)δYX=1,X=Y,0,X≠Y,1YX=1,X⊆Y,0,X⊄Y,where X,Y may be scalars, vectors, or sets.The probability density of GLMB RFS is distributed as follows:(5)πX=ΔX∑c∈CwcLXpc·X,where C is the discrete index set, h·X=∏ξ∈Xhξ, and
(6)∑L∈L∑c∈CwcL=1,∫ξ∈Xpcξ,ldξ=1.


## 2.2.2. The Labeled Likelihood of Group Target

The measurement likelihood function of the multigroup target is given by [1]
(7)gZX=gCZ∑i=1X+1∑UZ∈PiZθ∈ΘUZψUZ·;θX,where Z is a finite measurement set, PiZ denotes the set of all the partitions dividing Z into i groups, and UZ∈PiZ is a particular partition of Z. θ: LX→0,1,⋯,UZ denotes the association map of the target for measurement, while θl=θl′>0, then l=l′. ΘUZ is the space of θ, and UθlZ denotes the element of UZ corresponding to label l.
(8)ψUZξ,l;θ=pDξ,lg~UθlZξ,lkUθlZ,θl>0,qDξ,l,θl=0.


## 2.2.3. The Prediction and Update of GGIW-GLMB

At timek, let γk∼GAℳγk;αk,βk be the measurement rate which follows Poisson distribution, xk∼Nxk;mk,P~k⊗Xk be the centroid state, and Xk∼ℐWXk;vk,Vk be the extension state. Then, the probability density function (PDF) of the group target state is given by
(9)pξk=pγkpxkpXk=GAℳγk;αk,βkNxk;mk,P~k⊗XkℐWXk;vk,Vk=GGℐWξk;ζk,where ξk≜γk,xk,Xk denotes the state to be estimated, ζk≜αk,βk,mk,P~k,vk,Vk denotes the parameter of the GGIW component, P~ denotes the state covariance of a single dimension, and Pk=P~k⊗Xk denotes the state covariance of all dimensions. Substitute (9) into (5), the distribution of GLMB RFS of the multigroup target could be obtained subsequently.Let the subscript “+” be the index of the next time; e.g.,s+ might denote the probability density, label, label space, etc., according to the meaning of s. Then, the expression of prediction and update of GGIW-GLMB is given by [1]
(1)
Prediction(10)π+X+=ΔX+∑c∈Cw+cLX+p+c·χ+,where w+c· and p+c· denote the weight of the predicted GLMB and the predicted probability density of the group target state, respectively.
(2)
Update(11)πX∣Z=ΔX∑c∈C∑i=1X+1∑UZ∈PiZθ∈ΘUZwUZc,θLXpc,θ·∣UZX,where wUZc,θ· and pc,θ·∣UZ denote the weight of the updated GLMB component and updated probability density of the group target sate, respectively.


## 3. GGIW-δ-GLMB Filter with Splitting

### 3.1. GGIW-δ-GLMB Prediction and Update

δ-GLMB is a special form of GLMB filter; let
(12)C=FL×Ξ,(13)wcL=wI,ςL=wI,ςδIL,(14)pc=pI,ς=pς,where I∈FL is the track label set of time k. Ξ is a discrete space, and ς∈Ξ represents a history of association maps up to time k−1, i.e., ς=θ1,⋯,θk−1.Substitute (12)-(14) into (10) and (11); the prediction and update of GGIW-δ-GLMB can be given by the following:
(1)
Prediction(15)π+X+=ΔX+∑I+,ς∈FL+×Ξw+I+,ςLX+δI+LX+p+ς·χ+,where
(16)w+I+,ςL+=wBL+∩BwSI+,ςL+∩L,p+ςξ+,l+=1LlpSςξ+,l++1−1LlpBξ+,l+,pSςξ+,l+=∫IdξΨSςl+,ΨSςl+=∫∫Idξ′dξ+,wSI,ςJ=ΨSςJ∑I⊆L1IJqSςI−JwI,ςJ,qSςl=∫qSξ,lpςξ,ldξ,1AB=1,B⊆A,0,others,I=pSξ,lfξ+∣ξ,lpςξ,l,where B denotes the label space for targets born at the next time; L denotes the label space for the current time, for a given label set L+; wBL+∩B represents the weight of the birth labels; and wSI+,ςL+∩L represents the weight of the surviving labels. ·,l denotes the labeled state of the group target, qSξ,l denotes the death probability, pς·,l is either the density pB·,l of a newborn group or the density pSς·,l+ of a surviving group computed from the prior density pς·,l via the single-target prediction with transition density f·∣ξ,l weighted by the probability of survival probability pSξ,l.
(2)
Update(17)πX∣Z=ΔX∑I,ς∈FL×Ξ∑i=1X+1∑UZ∈PiZθ∈ΘUZwUZI,ς,θLXδILXpς,θ·∣UZX,where
(18)wUZI,ς,θL=wI,ςLΨUZς,θL∑I,ς∈FL×Ξ∑J⊆L∑i=1X+1∑UZ∈PiZθ∈ΘUZwI,ςJΨUZς,θJ,pς,θξ,l∣UZ=pςξ,lψUZξ,l;θΨUZς,θl,ΨUZς,θl=∫pς,θξ,lψUZξ,l;θdξ,where ψUZξ,l;θ is defined as (8). While pDξ,l=pDl,
(19)ΨUZς,θl=pDlΨℵςl,UθlZkUθlZ,where ℵςl denotes the parameters of the predicted GGIW density of the group target with label l within ς. Let
(20)Ψk,γj,W=Γαk∣kj,Wβk∣k−1jαk∣k−1jΓαk∣k−1jβk∣kj,Wαk∣kj,WW!,Ψk,x,Xj,W=1πWWSk∣k−1j,Wd/2Vk∣k−1jvk∣k−1j/2Vk∣kjvk∣kj,W/2Γdvk∣kj,W/2Γdvk∣k−1j/2.Then, Ψkℵςl,W=Ψk,γℵςl,WΨk,x,Xℵςl,W.


### 3.2. GGIW-δ-GLMB Prediction and Update with Splitting

Consider the splitting model in references [18, 30], we make the following assumptions:A1: The group target splits at the prediction phase.A2: Let the dimension of the extension state be d. As the split may happen in any dimension, if a parent track with label l∈L at time k splits 2 subgroups (split pair) at time k+1, then, each subgroup takes on the label lT,i,j=l,k+1,i,j, while i=1,…,d,j=1,2, and the parameter i is the same for the split pair. Let Ll,T=lT,i,ji=1,j=1d,2, Ll,T,i=lT,i,jj=12.Different from theδ-GLMB spawning model of the point target in reference [30], neither states of the split subgroups are identical with the expected value of its parent in this paper. That is to say, the split subgroups and its parent are mutually exclusive; e.g., when the formation of the UAV cluster is split into two subgroups, the group extension and measurement rate are greatly changed. While in reference [30], one of the split subgroups has an identical expected value on state to its parent e.g., after launching the missle, the extended state of the aircraft is almost unchanged. ② In reference [30], the split and the survival of the track are in parataxis; there are 3 possible states of a current target at the next time: survival, death, and split. But in this paper, survival is the prerequisite of the split, which means only a survival track could split. Thus, the two-split model applies for different tracking scenarios.For a preexistent group, it can be survival or death; if survived, it may split or not in the next time. Suppose thatLB,LS is the label set of a newborn group and the survival group, respectively. LS denotes the label space of the survival group targets, SNT,ST are the label space of the survival but nonsplit groups and the split subgroups, respectively, then, with LT=Ll,T and LT∈ST. L+=SNT∪ST∪B=LS∪B, SNT∩ST=ST∩B=SNT∩B=∅, and L∩LS=SNT, and according to the assumption A2, all labels in LS derive from L. Given the labeled state of a group target of current time to be ξ,l, if it survived, it will either split with probability pTl, and probability density 1LSLl,T,ifTξ1:2,i,+∣ξ,l≜1LSLl,T,ifTξ1,i,+,ξ2,i,+∣ξ,l, i=1,…,d, or it does not with probability qTl=1−pTl, and evolution probability density fNTξ+∣ξ,lδll+. Assume the nonlabeled state of a split pair to be ξt,+i,κ=γt,+,xt,+,X+tt=12, i=1,…,d. ξ1:2,+i,κ,Ll,T,i≜ξt,+i,κ,Ll,T,i,jt=12 denotes the labeled state set of all possible split subgroups of track l, with 0<κ<1 denoting a split parameter. ξ+,l denotes the predicted labeled state of track l at the next time while it does not split. We model the set XS of the multigroup target which survived as conditional LMB RFS distributed according to
(21)fSXS∣X=ΔXΔXS1LSLXSΦXS;·X.LetlS∈LS
(22)ΦXS;ξ,l=∑ξ,l∈Xξ+,l+∈XSδlSl+pSξ,l∑i=1d1LXSLl,T,ipTlfTξ1:2,+i,κ∣ξ,l+δll+qTlfNTξ+∣ξ,l+1−1LXSlSqSξ,l,where lS denotes the label of the survival target. Thus, while splitting, ΦXS;ξ,l is a joint distribution.(1) Prediction. Given the multigroup target posterior to be a δ-GLMB, the predicted labeled multigroup target density is given by
(23)π+X+=ΔX+∑I+,ς∈FL+×Ξw+I+,ςLX+δI+LX+p+I+,ςX+.AsL+=LS∪B, and LS∩B=∅. Let XB=X+−XS denotes the set of a newborn group, then
(24)w+I+,ςLX+=wBI+∩BwI,ςLX1JNT∪JTLXS,p+I+,ςX+=pBXB∏l∈IΦXS;·,l,pς·,l,where JNT∈SNT and JT∈ST are originated from the current component. The proof is given in Appendix A.The subgroups of different parent tracks and the newborn group target are mutually disjoint; however, both the subgroups from a parent are conditioned on the same state, the probability density of a split pair inp+I,ςX+ is correlated, and thus the labeled multigroup target density (23) is not δ-GLMB. Reference [30] approximates the label-conditioned joint densities of the labeled multitarget density by the product of its marginal to form an approximate δ-GLMB density; however, it is not easy to realize as the state of group target is modeled as GGIW. In this paper, the analytical minimization of the Kullback-Leibler divergence (KL-DIV) is used to approximate the true density of the joint probability density. Given the parameter κ, let the approximation of the joint probability density for the split pair be p^Tςξ1:2,+,Ll,T≜p^Tςξ1,+,lT,i,1,ξ2,+,lT,i,2, and p^Tςξ1:2,+,Ll,T equating to the expected values of the joint probability density
(25)p^Tςξ1:2,+,Ll,T≈∑i=1dδiGGℐWξ1,+i,κ;ζ1,+i,κGGℐWξ2,+i,κ;ζ2,+i,κ=∑i=1dδi∏j=1l∈I2pTςξj,+i,κ,lT,i,j,where
(26)ζ1,+i,κ=α+,β+κ,m1,+i,κ,P~1,+i,κ,v1,+,κV+i,κ,ζ2,+i,κ=α+,β+1−κ,m2,+i,κ,P~2,+i,κ,v2,+,1−κV+i,κ,m1,+i,κ=m++1−κeiHTvi,m2,+i,κ=m+−κeiHTvi,P~t,+i,κ=P~+,t=1,2,vt,+i,κ=v+,t=1,2,pX+1/2≈ℐWX+1/2;v+,V+i,κ,X+1/2=11+2κκ−1X+−κκ−11+2κκ−1X+12,X+12=Hx1,+−x2,+x1,+−x2,+THT,pX+1≈ℐWX+1;v+i,κ,κV+i,κ,pX+2≈ℐWX+2;v+i,κ,1−κV+i,κ.ζ1,+i,κ and ζ2,+i,κ denote the state parameters of the subgroups split from the ith dimension of a current group target, a different value in κ represents a different split hypothesis, and ei and vi are the ith eigenvalues and eigenvectors of X+, respectively. Respectively, ξ+=γ+,x+,X+ and ζ+=α+,β+,m+,P~+,v+,V+are the prediction state and the corresponding parameters of the next time while the current group does not split. The poof of (25) can be a reference in [18].Hence, approximating the joint density of the subgroups by (25), the approximation of the predicted δ-GLMB is given by
(27)π+X+=ΔX+∑I,ς∈FL×Ξw+I+,ςLX+δI+LX+p~+I+,ς·X+,p~+I+,ς·=1STLXTpTς·+1SNTLXNTpNTς·+1BLXBpB·,where pNTς denotes that the predicted probability density of the group does not split. The proof is given in Appendix B. Note that the split pair appears or not simultaneously, and they share a common weight. Only if I+=LX+, the inner summand is nonzero.(2) Update. With the GGIW-δ-GLMB prior density given as (27), the GGIW-δ-GLMB posterior density of multiple splitting groups is given as follows:
(28)πX∣Z=ΔX∑I,ς∈ℱL×Ξ∑i=1X+1∑UZ∈PiZθ∈ΘUZwUZI,ς,θLXδILXp~ς,θ·∣UZX,where
(29)p~ς,θξ,ℓ∣UZ=p~I,ςξ,ℓψUZξ,ℓ;θΨ~UZς,θℓ,(30)Ψ~UZς,θℓ=∫p~ς,θξ,ℓ∣UZψUZξ,ℓ;θdξ.


## 3.1. GGIW-δ-GLMB Prediction and Update

δ-GLMB is a special form of GLMB filter; let
(12)C=FL×Ξ,(13)wcL=wI,ςL=wI,ςδIL,(14)pc=pI,ς=pς,where I∈FL is the track label set of time k. Ξ is a discrete space, and ς∈Ξ represents a history of association maps up to time k−1, i.e., ς=θ1,⋯,θk−1.Substitute (12)-(14) into (10) and (11); the prediction and update of GGIW-δ-GLMB can be given by the following:
(1)
Prediction(15)π+X+=ΔX+∑I+,ς∈FL+×Ξw+I+,ςLX+δI+LX+p+ς·χ+,where
(16)w+I+,ςL+=wBL+∩BwSI+,ςL+∩L,p+ςξ+,l+=1LlpSςξ+,l++1−1LlpBξ+,l+,pSςξ+,l+=∫IdξΨSςl+,ΨSςl+=∫∫Idξ′dξ+,wSI,ςJ=ΨSςJ∑I⊆L1IJqSςI−JwI,ςJ,qSςl=∫qSξ,lpςξ,ldξ,1AB=1,B⊆A,0,others,I=pSξ,lfξ+∣ξ,lpςξ,l,where B denotes the label space for targets born at the next time; L denotes the label space for the current time, for a given label set L+; wBL+∩B represents the weight of the birth labels; and wSI+,ςL+∩L represents the weight of the surviving labels. ·,l denotes the labeled state of the group target, qSξ,l denotes the death probability, pς·,l is either the density pB·,l of a newborn group or the density pSς·,l+ of a surviving group computed from the prior density pς·,l via the single-target prediction with transition density f·∣ξ,l weighted by the probability of survival probability pSξ,l.
(2)
Update(17)πX∣Z=ΔX∑I,ς∈FL×Ξ∑i=1X+1∑UZ∈PiZθ∈ΘUZwUZI,ς,θLXδILXpς,θ·∣UZX,where
(18)wUZI,ς,θL=wI,ςLΨUZς,θL∑I,ς∈FL×Ξ∑J⊆L∑i=1X+1∑UZ∈PiZθ∈ΘUZwI,ςJΨUZς,θJ,pς,θξ,l∣UZ=pςξ,lψUZξ,l;θΨUZς,θl,ΨUZς,θl=∫pς,θξ,lψUZξ,l;θdξ,where ψUZξ,l;θ is defined as (8). While pDξ,l=pDl,
(19)ΨUZς,θl=pDlΨℵςl,UθlZkUθlZ,where ℵςl denotes the parameters of the predicted GGIW density of the group target with label l within ς. Let
(20)Ψk,γj,W=Γαk∣kj,Wβk∣k−1jαk∣k−1jΓαk∣k−1jβk∣kj,Wαk∣kj,WW!,Ψk,x,Xj,W=1πWWSk∣k−1j,Wd/2Vk∣k−1jvk∣k−1j/2Vk∣kjvk∣kj,W/2Γdvk∣kj,W/2Γdvk∣k−1j/2.Then, Ψkℵςl,W=Ψk,γℵςl,WΨk,x,Xℵςl,W.


## 3.2. GGIW-δ-GLMB Prediction and Update with Splitting

Consider the splitting model in references [18, 30], we make the following assumptions:A1: The group target splits at the prediction phase.A2: Let the dimension of the extension state be d. As the split may happen in any dimension, if a parent track with label l∈L at time k splits 2 subgroups (split pair) at time k+1, then, each subgroup takes on the label lT,i,j=l,k+1,i,j, while i=1,…,d,j=1,2, and the parameter i is the same for the split pair. Let Ll,T=lT,i,ji=1,j=1d,2, Ll,T,i=lT,i,jj=12.Different from theδ-GLMB spawning model of the point target in reference [30], neither states of the split subgroups are identical with the expected value of its parent in this paper. That is to say, the split subgroups and its parent are mutually exclusive; e.g., when the formation of the UAV cluster is split into two subgroups, the group extension and measurement rate are greatly changed. While in reference [30], one of the split subgroups has an identical expected value on state to its parent e.g., after launching the missle, the extended state of the aircraft is almost unchanged. ② In reference [30], the split and the survival of the track are in parataxis; there are 3 possible states of a current target at the next time: survival, death, and split. But in this paper, survival is the prerequisite of the split, which means only a survival track could split. Thus, the two-split model applies for different tracking scenarios.For a preexistent group, it can be survival or death; if survived, it may split or not in the next time. Suppose thatLB,LS is the label set of a newborn group and the survival group, respectively. LS denotes the label space of the survival group targets, SNT,ST are the label space of the survival but nonsplit groups and the split subgroups, respectively, then, with LT=Ll,T and LT∈ST. L+=SNT∪ST∪B=LS∪B, SNT∩ST=ST∩B=SNT∩B=∅, and L∩LS=SNT, and according to the assumption A2, all labels in LS derive from L. Given the labeled state of a group target of current time to be ξ,l, if it survived, it will either split with probability pTl, and probability density 1LSLl,T,ifTξ1:2,i,+∣ξ,l≜1LSLl,T,ifTξ1,i,+,ξ2,i,+∣ξ,l, i=1,…,d, or it does not with probability qTl=1−pTl, and evolution probability density fNTξ+∣ξ,lδll+. Assume the nonlabeled state of a split pair to be ξt,+i,κ=γt,+,xt,+,X+tt=12, i=1,…,d. ξ1:2,+i,κ,Ll,T,i≜ξt,+i,κ,Ll,T,i,jt=12 denotes the labeled state set of all possible split subgroups of track l, with 0<κ<1 denoting a split parameter. ξ+,l denotes the predicted labeled state of track l at the next time while it does not split. We model the set XS of the multigroup target which survived as conditional LMB RFS distributed according to
(21)fSXS∣X=ΔXΔXS1LSLXSΦXS;·X.LetlS∈LS
(22)ΦXS;ξ,l=∑ξ,l∈Xξ+,l+∈XSδlSl+pSξ,l∑i=1d1LXSLl,T,ipTlfTξ1:2,+i,κ∣ξ,l+δll+qTlfNTξ+∣ξ,l+1−1LXSlSqSξ,l,where lS denotes the label of the survival target. Thus, while splitting, ΦXS;ξ,l is a joint distribution.(1) Prediction. Given the multigroup target posterior to be a δ-GLMB, the predicted labeled multigroup target density is given by
(23)π+X+=ΔX+∑I+,ς∈FL+×Ξw+I+,ςLX+δI+LX+p+I+,ςX+.AsL+=LS∪B, and LS∩B=∅. Let XB=X+−XS denotes the set of a newborn group, then
(24)w+I+,ςLX+=wBI+∩BwI,ςLX1JNT∪JTLXS,p+I+,ςX+=pBXB∏l∈IΦXS;·,l,pς·,l,where JNT∈SNT and JT∈ST are originated from the current component. The proof is given in Appendix A.The subgroups of different parent tracks and the newborn group target are mutually disjoint; however, both the subgroups from a parent are conditioned on the same state, the probability density of a split pair inp+I,ςX+ is correlated, and thus the labeled multigroup target density (23) is not δ-GLMB. Reference [30] approximates the label-conditioned joint densities of the labeled multitarget density by the product of its marginal to form an approximate δ-GLMB density; however, it is not easy to realize as the state of group target is modeled as GGIW. In this paper, the analytical minimization of the Kullback-Leibler divergence (KL-DIV) is used to approximate the true density of the joint probability density. Given the parameter κ, let the approximation of the joint probability density for the split pair be p^Tςξ1:2,+,Ll,T≜p^Tςξ1,+,lT,i,1,ξ2,+,lT,i,2, and p^Tςξ1:2,+,Ll,T equating to the expected values of the joint probability density
(25)p^Tςξ1:2,+,Ll,T≈∑i=1dδiGGℐWξ1,+i,κ;ζ1,+i,κGGℐWξ2,+i,κ;ζ2,+i,κ=∑i=1dδi∏j=1l∈I2pTςξj,+i,κ,lT,i,j,where
(26)ζ1,+i,κ=α+,β+κ,m1,+i,κ,P~1,+i,κ,v1,+,κV+i,κ,ζ2,+i,κ=α+,β+1−κ,m2,+i,κ,P~2,+i,κ,v2,+,1−κV+i,κ,m1,+i,κ=m++1−κeiHTvi,m2,+i,κ=m+−κeiHTvi,P~t,+i,κ=P~+,t=1,2,vt,+i,κ=v+,t=1,2,pX+1/2≈ℐWX+1/2;v+,V+i,κ,X+1/2=11+2κκ−1X+−κκ−11+2κκ−1X+12,X+12=Hx1,+−x2,+x1,+−x2,+THT,pX+1≈ℐWX+1;v+i,κ,κV+i,κ,pX+2≈ℐWX+2;v+i,κ,1−κV+i,κ.ζ1,+i,κ and ζ2,+i,κ denote the state parameters of the subgroups split from the ith dimension of a current group target, a different value in κ represents a different split hypothesis, and ei and vi are the ith eigenvalues and eigenvectors of X+, respectively. Respectively, ξ+=γ+,x+,X+ and ζ+=α+,β+,m+,P~+,v+,V+are the prediction state and the corresponding parameters of the next time while the current group does not split. The poof of (25) can be a reference in [18].Hence, approximating the joint density of the subgroups by (25), the approximation of the predicted δ-GLMB is given by
(27)π+X+=ΔX+∑I,ς∈FL×Ξw+I+,ςLX+δI+LX+p~+I+,ς·X+,p~+I+,ς·=1STLXTpTς·+1SNTLXNTpNTς·+1BLXBpB·,where pNTς denotes that the predicted probability density of the group does not split. The proof is given in Appendix B. Note that the split pair appears or not simultaneously, and they share a common weight. Only if I+=LX+, the inner summand is nonzero.(2) Update. With the GGIW-δ-GLMB prior density given as (27), the GGIW-δ-GLMB posterior density of multiple splitting groups is given as follows:
(28)πX∣Z=ΔX∑I,ς∈ℱL×Ξ∑i=1X+1∑UZ∈PiZθ∈ΘUZwUZI,ς,θLXδILXp~ς,θ·∣UZX,where
(29)p~ς,θξ,ℓ∣UZ=p~I,ςξ,ℓψUZξ,ℓ;θΨ~UZς,θℓ,(30)Ψ~UZς,θℓ=∫p~ς,θξ,ℓ∣UZψUZξ,ℓ;θdξ.


## 4. Simplification of the Approach Proposed above

### 4.1. Splitting Detection

As for GGIW-δ-GLMB with splitting derived in Section 3, if there are n GGIW in a δ-GLMB component, then the number of predicted GLMB will be 2+dκn, where κ denotes the number of the split hypotheses in each dimension of the state. Hence, as the number of the group target increases, the number of δ-GLMB increases exponentially. In practice, if it does, the split will only happen in one or just a few moments in a tracking process; splitting prediction in every iteration may be computationally infeasible for practical application. One possible solution is to assume that only one group target may split in a particular moment. As usual, the probability of group target splitting is small, and the probability of splitting more than 1 group target is much smaller; so, this assumption is reasonable. Then, for a δ-GLMB with n GGIW, the number of the predicted δ-GLMB in the next time is 2n−1n2+dκ. In this paper, we propose a splitting detection approach based on the following assumption:A3: While not splitting, the formation of the group, the target number within a group, and the spacing among them do not change significantly, e.g., formation airplane cruise.According to A3, for a group, the ratio of the estimated area to the number of measurements should not vary much, if its motion is stationary. If the ratio bigger than a given threshold, it could be a sign of splitting. Given the PDF (9), the ratio is computed by [31].
(31)US,k+1∣k=4π∏l=1dek+1∣klγk+1∣k,where ek+1∣kl denotes the lth eigenvalue of Xk+1∣k. In reference [31], a threshold UT is given, if US,k+1∣k>UT, the target is thought to split, and the target is divided into multiple ellipses; however, as the extension area and the number of group target vary much in different scenarios, it is not easy to preset a threshold UT. Thus, here we use the change of the ratio between the adjacent time as the basis of split detection
(32)RU,k+1=US,k+1∣k+1US,k∣k,where US,k∣k is computed by Xk∣k and γk∣k. Usually, if RU,k+1≤1, we let pT=0; if RU,k+1>1, we could set the split probability according to the specific condition, e.g.,
(33)pT=1,1.5<RU,k+1,0.5,1<RU,k+1≤1.5,0,RU,k+1≤1.


### 4.2. Birth Components Combination

There are two main models for target birth: one is a static model, which means that new targets could only be initiated around predetermined locations, and the other is an adaptive target birth model that allows for new targets to appear anywhere in the state space. Hence, the latter is more flexible in the scenarios with no priori information about the position of the newborn target [32].The adaptive birth density in this paper resembles the adaptive birth density of the standard LMB filter [32]. The difference is that there are many possible partitions to the measurement in multigroup target tracking and different partitions corresponding to different combination of the targets; many of the measurements in different partitions are duplicate or corresponding a same target. As each measurement set may be generated by a potential newborn target [32], modeling these duplicate measurements or measurements from the same target repeatedly would be a huge waste of computation. In view of this, we consider validating the kinematic states for the possible target. First, for the states of possible newborn targets falling inside the gate, only the one with maximum weight is remained. The gate criterion is given by
(34)Gbb=mk,b1−mk,b2TP0−1mk,b1−mk,b2≤η12,where mk,b1 and mk,b2 are potential newborn targets, while (34) is satisfied, then mk,b1 and mk,b2 are combined, and the combined state is identical to the one with greater birth probability. Gbb is a set of potential newborn targets after combination. η12 is the threshold calculated by the inverse cumulative distribution conditioned on a preset confidence probability, and we set the probability to be 0.9 in this paper.Second, for a potential newborn target inGbb and a preexistent target falling inside a gate, we think it is a preexistent one and do not model it as a newborn target anymore. The gate criterion is given by
(35)Gbe=mk,b−mk∣kTPk∣k−1mk,b−mk∣k≤η22mk,b∈Gbb,where Gbe denotes the set of possible newborn targets falling inside the gate η22 for the updated component mk∣k. Pk∣k is the state covariance. Thus, the newborn target set of time k is Gb=Gbb−Gbe.


## 4.1. Splitting Detection

As for GGIW-δ-GLMB with splitting derived in Section 3, if there are n GGIW in a δ-GLMB component, then the number of predicted GLMB will be 2+dκn, where κ denotes the number of the split hypotheses in each dimension of the state. Hence, as the number of the group target increases, the number of δ-GLMB increases exponentially. In practice, if it does, the split will only happen in one or just a few moments in a tracking process; splitting prediction in every iteration may be computationally infeasible for practical application. One possible solution is to assume that only one group target may split in a particular moment. As usual, the probability of group target splitting is small, and the probability of splitting more than 1 group target is much smaller; so, this assumption is reasonable. Then, for a δ-GLMB with n GGIW, the number of the predicted δ-GLMB in the next time is 2n−1n2+dκ. In this paper, we propose a splitting detection approach based on the following assumption:A3: While not splitting, the formation of the group, the target number within a group, and the spacing among them do not change significantly, e.g., formation airplane cruise.According to A3, for a group, the ratio of the estimated area to the number of measurements should not vary much, if its motion is stationary. If the ratio bigger than a given threshold, it could be a sign of splitting. Given the PDF (9), the ratio is computed by [31].
(31)US,k+1∣k=4π∏l=1dek+1∣klγk+1∣k,where ek+1∣kl denotes the lth eigenvalue of Xk+1∣k. In reference [31], a threshold UT is given, if US,k+1∣k>UT, the target is thought to split, and the target is divided into multiple ellipses; however, as the extension area and the number of group target vary much in different scenarios, it is not easy to preset a threshold UT. Thus, here we use the change of the ratio between the adjacent time as the basis of split detection
(32)RU,k+1=US,k+1∣k+1US,k∣k,where US,k∣k is computed by Xk∣k and γk∣k. Usually, if RU,k+1≤1, we let pT=0; if RU,k+1>1, we could set the split probability according to the specific condition, e.g.,
(33)pT=1,1.5<RU,k+1,0.5,1<RU,k+1≤1.5,0,RU,k+1≤1.


## 4.2. Birth Components Combination

There are two main models for target birth: one is a static model, which means that new targets could only be initiated around predetermined locations, and the other is an adaptive target birth model that allows for new targets to appear anywhere in the state space. Hence, the latter is more flexible in the scenarios with no priori information about the position of the newborn target [32].The adaptive birth density in this paper resembles the adaptive birth density of the standard LMB filter [32]. The difference is that there are many possible partitions to the measurement in multigroup target tracking and different partitions corresponding to different combination of the targets; many of the measurements in different partitions are duplicate or corresponding a same target. As each measurement set may be generated by a potential newborn target [32], modeling these duplicate measurements or measurements from the same target repeatedly would be a huge waste of computation. In view of this, we consider validating the kinematic states for the possible target. First, for the states of possible newborn targets falling inside the gate, only the one with maximum weight is remained. The gate criterion is given by
(34)Gbb=mk,b1−mk,b2TP0−1mk,b1−mk,b2≤η12,where mk,b1 and mk,b2 are potential newborn targets, while (34) is satisfied, then mk,b1 and mk,b2 are combined, and the combined state is identical to the one with greater birth probability. Gbb is a set of potential newborn targets after combination. η12 is the threshold calculated by the inverse cumulative distribution conditioned on a preset confidence probability, and we set the probability to be 0.9 in this paper.Second, for a potential newborn target inGbb and a preexistent target falling inside a gate, we think it is a preexistent one and do not model it as a newborn target anymore. The gate criterion is given by
(35)Gbe=mk,b−mk∣kTPk∣k−1mk,b−mk∣k≤η22mk,b∈Gbb,where Gbe denotes the set of possible newborn targets falling inside the gate η22 for the updated component mk∣k. Pk∣k is the state covariance. Thus, the newborn target set of time k is Gb=Gbb−Gbe.


## 5. Simulation Results

There are two group targets in a cluttered scenario. For group 1, initialize the kinematic state to be1900m,1700m,−100m/s,0m/s,0m/s2,0m/s2T, the measurement rate to be 30, and the extension state to be a,b=20,50m, where a,b are semi-minor and semi-major axes of the extension ellipse, respectively. The target bears a constant velocity from time k=1 through 5 and split at time k=6, the two split subgroups moving in constant turning models with turning rates π/40 rad/s and −π/40 rad/s, respectively, from time k=6 through 11, then they both move in constant velocity model, and the track is terminated at time k=14 and 15, respectively. The semi-major and semi-minor axes of the subgroups satisfy ai,bi=20,25m, with i=1,2, and the measurement rate is 15, respectively. The initialized kinematic state of the group 2 is −980m,−100m,80m/s,−20m/s,0m/s2,0m/s2T, the measurement rate is 30, and extension state satisfies a,b=15,40m. From time k=3 through 9, the target moves in a constant velocity model and is split into two subgroups with measurement rates 20 and 10, respectively, at time k=10. The first split subgroup of group 2 satisfies a1,b1=15,30m and moves with a constant turning rate with the turning rate as π/20 rad/s from time k=10 through 15, then moving in a constant velocity model, and the track is terminated at time k=17. The second one satisfies a2,b2=15,10m and moves in a constant turning model with the turning rate as −π/60 rad/s, the track being terminated at time k=18. The clutter follows a Poisson distribution with the mean value of 5 and distributed in a measurement space of −2000,3000×−3000,3000. The detection probability of the sensor is PD=0.99, and the probability of target survival is PS=0.95.
(36)F~k∣k−1=1TT2200T00e−T/θ.AndR=diag11, Qk∣k−1=Σ21−e−2T/θdiag001⊗Xk, where T=1s is the sample time, and θ=1s is the maneuvering correlation time, Σ=0.1m/s2.Figure1 shows the true tracks and the noised measurements in clutter, where “∗” denotes the clutter, “.” denotes the noised measurement of the target in the group, the ellipse denotes the distribution area of the target extension, and “→”denotes the motion orientation of the group.Figure 1
The true track and the noised measurement in cluttered condition.The OSPA error of 5 approaches is given in Figure2, where M1 denotes the δ-GLMB filter with the static birth approach, and no splitting is considered; M2 denotes the approach in reference [18], where group splitting is considered based on the PHD filter; M3 denotes a δ-GLMB filter with adaptive birth, but splitting is not modeled; M4 denotes a δ-GLMB filter based on the group splitting model with a static birth; and M5 denotes the proposed approach in this paper. For all of these approaches, if adaptive birth or a δ-GLMB based on the splitting model is adopted, the corresponding strategy in Section 4 is adapted.Figure 2
OSPA errors of M1-M5. (a)-(d) are the general OSPA error, kinematic state OSPA error, extension state OSPA error, and measurement rate OSPA error, respectively.
(a)
(b)
(c)
(d)The kinematic states of the birth group for the approaches with static birth are true kinematic states, and other parameters of these approach are initialed asP~0=diag1002502502, v0=10, V0=diag5050, α0=10, and β0=1. For the approaches with adaptive birth, the birth probability is rbirth,i=0.03,i=1,2, and the initialed kinematic state V0 is computed online; the expected number of the birth group in each iteration is λB=∑i=12rbirth,i=0.06. While the splitting model is introduced, κ=0.25,0.5,0.75. If there is no splitting detection, pT=0.03, while the splitting detection in Section 4 is conducted; if RU>1.1, pT=0.5, otherwise, pT=0.Figure2(a) shows the general OSPA error of the approaches above mentioned, where the general OSPA is a weighted sum of the kinematic state, extension state, and measurement rate with each weighted by 0.8, 0.1, and 0.1, respectively. It can be seen from Figures 2(a)–2(d) that the performance of M1 in OSPA error is the worst, which is because the approach does not model the group split and cannot track the split track efficiently. M3 and M5 have larger OSPA errors than M1, M2, and M4 in the first few periods, because there is a time delay when new targets are generated by an adaptive birth method. All of M2-M5 could track the group split, and the OSPA error increased in a few periods just behind the splitting time, because the split could not be detected by the approaches in the first few periods after group splitting. The adaptive birth approach could track the split track even though group split is not modeled, as the split track is tracked as a new track; however, there will be a big error in the initiation period of the new track. Compared with M2, M4 is superior in a whole, which indicates that the δ-GLMB filter based on a splitting model is better than the PHD one. While the target birth does not happen, M5 shows a better performance than M4 does, as the adaptive birth of M5 could reinitiate a new track when missing the track reappearance, while M4 could not.The estimated number of groups is shown in Figure3, where “truth” denotes the truth number of groups. It can be seen that M1 could not estimate the number of the split groups efficiently. For the other 4 approaches, there is a time delay while split happened, and the approaches with static birth have no time delay in target birth time. M5 has the most accurate estimation of the number of groups.Figure 3
The estimated numbers of groups.Theδ-GLMB filter could obtain the track information of the target by labeling the target track, but the labeling error may exist; to measure the accuracy and stability of the approaches in labeling, we defined the labeling OSPA error. We define a label state as tlk=tst,tsp,tdeT, where tst denotes the birth time of the track, tsp denotes the split time of the track, and tde denotes the death time of the track. Let tsp=tst if splitting does not occur. Then, we could calculate the OSPA errors of the label states. Figure 4 shows the track labeling performance of M1 and M3-M5.Figure 4
Labeling performance of the algorithms.It can be seen from Figure4 that M4 and M5 show a higher accuracy than M3, because when split happened, M4 and M5 treat the split tracks as a continuity of the original track, which is closer to the evolution of the target motion, while M3 treats the split tracks as a new track and neglects the correlations between the split tracks and the original track. We can see that a study on the label OSPA error will help the tracking system to grasp the evolution of the target motion and assist in the inference of the target intention.


## 6. Conclusions

Aδ-GLMB filter based on the split model for multigroup target tracking is proposed in this paper, which derives the δ-GLMB form for group target tracking, and the introduction of the split model makes the group split addressed. There are two simplifications: split detection and component combination for adaptive birth, of which the former avoids the prediction of the target split in every iteration, and the latter decreases the number of the birth component, and both of the two simplify the computation significantly. Finally, track labeling the OSPA error is proposed, which provides a basis for track stability and accuracy measure. As the combination and splitting usually appear as a pair in group target motion, in future, the incorporation of the combination and the proposed filter may be worth having a try. Moreover, as the characteristics of group target is hard to obtain in the context of clutter and complex background conditions in prior, the implementation of machine leaning methods [33–36] in group target tracking will be helpful for target identification, classification, and measurement partition and thus improve the tracking accuracy and efficiency; more attention is needed in the combination of machine learning method and group target tracking method in the future.



---

*Source: 9278725-2019-05-19.xml*