Next Article in Journal
RADAR-Vegetation Structural Perpendicular Index (R-VSPI) for the Quantification of Wildfire Impact and Post-Fire Vegetation Recovery
Next Article in Special Issue
DU-CG-STAP Method Based on Sparse Recovery and Unsupervised Learning for Airborne Radar Clutter Suppression
Previous Article in Journal
Automatic Deployment of Convolutional Neural Networks on FPGA for Spaceborne Remote Sensing Application
Previous Article in Special Issue
Dynamic Antenna Selection for Colocated MIMO Radar
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalized Labeled Multi-Bernoulli Multi-Target Tracking with Doppler-Only Measurements

1
School of Computer Science, Shaanxi Normal University, Xi’an 710119, China
2
Independent Consultant, Anacortes, WA 98221, USA
3
Academy of Advanced Interdisciplinary Research, Xidian University, Xi’an 710071, China
4
National Laboratory of Radar Signal Processing, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(13), 3131; https://doi.org/10.3390/rs14133131
Submission received: 5 May 2022 / Revised: 13 June 2022 / Accepted: 24 June 2022 / Published: 29 June 2022
(This article belongs to the Special Issue Small or Moving Target Detection with Advanced Radar System)

Abstract

:
The paper addresses the problem of tracking multiple targets with Doppler-only measurements in multi-sensor systems. It is well known that the observability of the target state measured using Doppler-only measurements is very poor, which makes it difficult to initialize the tracking target and produce the target trajectory in any tracking algorithm. Within the framework of random finite sets, we propose a novel constrained admissible region (CAR) based birth model that instantiates the birth distribution using Doppler-only measurements. By combining physics-based constraints in the unobservable subspace of the state space, the CAR based birth model can effectively reduce the ambiguity of the initial state. The CAR based birth model combines physics-based constraints in the unobservable subspace of the state space to reduce the ambiguity of the initial state. We implement the CAR based birth model with the generalized labeled multi-Bernoulli tracking filter to demonstrate the effectiveness of our proposed algorithm in Doppler-only tracking. The performance of the proposed approach is tested in two simulation scenarios in terms of the optimal subpattern assignment (OSPA) error, OSPA ( 2 ) error, and computing efficiency. The simulation results demonstrate the superiority of the proposed approach. Compared to the approach taken by the state-of-the-art methods, the proposed approach can at most reduce the OSPA error by 58.77%, reduce the OSPA ( 2 ) error by 43.51%, and increase the computing efficiency by 9.56 times in the first scenario. In the second scenario, the OSPA error is reduced by 62.80%, the OSPA ( 2 ) error is reduced by 43.65%, and the computing efficiency is increased by 2.61 times at most.

Graphical Abstract

1. Introduction

Tracking of multiple moving targets with Doppler measurements in multi-sensor systems has attracted interest from researchers during the past few decades [1,2,3,4,5,6,7,8,9]. Compared with other forms of measurements, Doppler is attractive for a number of reasons. For instance, Doppler sensors are much cheaper, since no hardware array is required. Another reason is that Doppler measurements are relatively accurate. In addition, the amount of information exchanged, or the bandwidth requirement between each sensor and the fusion center and/or between sensors, is significantly lower than that of other measurement types [5]. However, the Doppler measurement cannot provide complete information about the target state. This issue is often referred to as single-sensor unobservability [6]. It is a challenging problem to accurately estimate the number of targets and extract their trajectories or states from Doppler-only measurements under clutter and missed detections.
In general, multiple sensors can be used to alleviate the unobservability problem via multi-sensor data fusion, as well as achieving multi-sensor diversity gain. Chan and Jardine [2] have shown that a unique target state can be obtained using at least four Doppler measurements from different sensors. In [10], Shames et al. algebraically derived the minimum number of Doppler measurements to provide finite solutions of the target state. A dual-stage approach based on observability decomposition for centralized tracking using Doppler-only measurements was proposed in [11]. These papers, however, only study the single target case, while the multi-target case is far more complex. One major issue that comes with the multi-target case is the ghost problem, which is a consequence of the unknown data association. While multiple targets are detected by each sensor, multiple associations of the corresponding measurements from all sensors often result in fictitious or ghost states [12].
Multi-sensor multi-target tracking (MTT) is a difficult problem because of measurement origin uncertainty, the presence of false alarms, missed detections, measurement errors, and birth and death of targets [13]. As the classical MTT approaches, the joint probabilistic data association (JPDA) [14] and multiple hypothesis tracking (MHT) [15,16] algorithms were first proposed in the 1970s. The JPDA and MHT produce tracks using data association and each track has a label or identity. Multi-scan, s-dimensional (sD) assignment, or multi-frame based assignment approaches for JPDA [17] and MHT [18,19] were proposed for improved data association accuracy. Other variants of the algorithms also exist. For example, the exact nearest neighbour JPDA filter [20] was proposed to circumvent the coalescence problem of the JPDA filter, where all data association hypotheses are drastically pruned except the one with the most significant or highest probability. Another variant of JPDA named JPDA* [21] was developed based on the selective pruning, and the best data association hypothesis is chosen to calculate the measurement-to-target probabilities. There are also some methods [22,23,24] that attempt to improve the accuracy of the marginalization by optimizing the posterior density. In recent years, the random finite set (RFS) [25,26] based MTT algorithms have attracted a lot of attention [27]. For the RFS based algorithms, the posterior multi-target density encapsulates multi-target state information based on evidence from the observed data and prior information on the motion of targets [25]. As more data are processed, the posterior density of the set of real tracks will be higher than that for ghost tracks. Earlier RFS based algorithms, including probability hypothesis density (PHD) [28], cardinalized PHD [29], and multi-Bernoulli [30] filters, were not developed to output the trajectory of each target. These algorithms are approximations of the multi-target Bayesian filter but are not real trackers in principle. The first RFS-based tracker was formulated using the labeled RFS, which led to the generalized labeled multi-Bernoulli (GLMB) filter [31,32]. Its multi-scan version [33,34], as well as approximations such as the labeled multi-Bernoulli filter [35], also produce tracks. A unique feature of the GLMB approach is that the approximate error bounds can be computed analytically [32]. Moreover, the GLMB filter can achieve linear complexity in the number of measurements using Gibbs sampling [36], enabling it to handle large-scale MTT problems [37].
A number of RFS-based solutions have been proposed for MTT with Doppler-only measurements. In [5], the PHD filter was applied to Doppler-only measurements in a passive multi-static setting. Tracking of multiple cooperative targets with the multi-Bernoulli filter using Doppler-only measurements was studied in [38]. These methods use the PHD and multi-Bernoulli filters which cannot produce the target trajectories. By assuming that targets are indistinguishable, the PHD and multi-Bernoulli filters model the set of target states as the RFS and estimate the multi-target state at each time step. However, the multi-target state history does not necessarily represent the multi-target trajectory. For more details of these filters, please refer to [25,26,28,30,39,40,41]. To address the trajectory, Papi considered an application of the GLMB filter in [42] for tracking multiple targets using Doppler measurements. In [43], we studied the sensor management problem in Doppler multi-sensor systems for MTT applications using the labeled multi-Bernoulli filter. Do and Nguyen [44] used the GLMB filter for tracking multiple marine vessels in Doppler multi-static systems, and the detection probability is kept as unknown. In [45], an efficient multi-sensor GLMB filter was developed and demonstrated on various sensor types including Doppler. These RFS-based methods all require prior knowledge of the target birth distribution (or initial-state distribution) and therefore can be restrictive in practical applications.
To solve the problem of unknown target birth distribution, some researchers developed measurement-driven birth models [46,47]. The measurement-driven birth model uses the received measurements to generate target births, and hence the prior knowledge of birth distribution is not necessary. However, the notorious poor observability of Doppler measurements makes this approach difficult to apply. In general, samples are used to approximate the birth distribution, and a large number of samples is required to cover the entire surveillance area. Within the RFS filtering framework, the accept-reject (A-R) method is widely used in the state-of-the-art methods to generate birth distributions [48,49,50]. In the A-R method, samples are used to approximate the birth distribution. A large number of samples are drawn from a diffuse prior covering the entire state space, and only those samples whose Doppler shifts are similar to the generated Doppler measurement are accepted. Since a large number of samples are first drawn and then rejected, this method is apparently inefficient. In addition, there could be a large difference between the estimated birth distribution and the true one. In some cases, the under-determined solution produces a birth distribution with sufficient ambiguity and makes the subsequent MTT process difficult.
In this paper, the problem of tracking multiple distinguishable and non-cooperative targets using Doppler measurements is studied. While this is an important practical problem, to the best of our knowledge, only a few publications exist on this topic. For real-world applications, we need a realistic measurement model that accounts for the measurement noise, clutter, and missed detections. Moreover, we need a tracker that can fuse Doppler-only measurements from multiple radars to alleviate the poor observability. The tracker must operate without the prior knowledge of the birth distribution, since it is not available in practice. To address these challenges, we propose a constrained admissible region (CAR) based adaptive birth model in the GLMB filter. Recent studies have presented various CAR based approaches (e.g., see [51,52]) for the generation of prior state distributions. However, existing approaches are all designed for the space-object tracking problem in astrodynamics. For the first time, we combine a generated Doppler measurement with a set of physics-based constraints to instantiate the birth distribution. This is achieved by first analytically determining the set of samples representing possible initial states for each single Doppler measurement based on the CAR, and then eliminating the infeasible samples using the multi-sensor GLMB update. We present the resulting GLMB tracking approach in detail in the sequential Monte Carlo (SMC) implementation and demonstrate its capability.
The organization of the paper is as follows: Section 2 and Section 3 present a description of the tracking model and a brief review of the GLMB recursion. In Section 4, the proposed GLMB tracking with the CAR based adaptive birth model is described in detail, outlining the general framework and the step-by-step implementation. Numerical simulations and results are presented in Section 5. Section 6 draws conclusions.

2. Target and Sensor Models

In the 2D surveillance area, the kinematic state of a moving target at time k is described by a vector [53,54]
x k = [ p k , x p ˙ k , x p k , y p ˙ k , y ] T ,
where the target position and velocity are given as p k = [ p k , x p k , y ] T and p ˙ k = [ p ˙ k , x p ˙ k , y ] T , respectively. The measurement time interval T is assumed constant and targets move with the nearly constant velocity (NCV) motion. The direct discrete-time kinematic model [55] is used, and the dynamic equation is as follows [55]
x k = F x k 1 + w k 1 ,
where
F = I 2 1 T 0 1 ,
with ⊗ denoting the Kronecker product [56], and w k 1 is the zero-mean white Gaussian process noise with covariance Q [55]
Q = I 2 σ w 2 T 4 4 T 3 2 T 3 2 T 2 ,
where σ w is the standard deviation (SD) of the process noise. Note that, if the discretized continuous-time kinematic model is used, the covariance Q of the process noise has a different expression. For more details of the kinematic model, please refer to [55].
The targets are observed by a Doppler multi-sensor system with N s sensors. In the multi-sensor system, target x k is illuminated by the sinusoidal waveform from sensor j located at s ( j ) = [ s x ( j ) , s y ( j ) ] T . If the signal is scattered by the target, sensor j will receive it and report a Doppler measurement as [57]
z k ( j ) = h k ( j ) ( x k ) + ε k ( j ) ,
where
h k ( j ) ( x k ) = 2 p ˙ k T p k s ( j ) p k s ( j ) f c , j c
is the true Doppler shift, f c , j is the carrier frequency for sensor j, c is the speed of light, and ε k ( j ) is the Gaussian measurement noise with zero mean and SD σ ε .
The Doppler measurement (5) is detected by sensor j with the detection probability p D ( j ) ( x k ) 1 . The false detections in the surveillance area are assumed independent of the target states [31]. Specifically, for sensor j, the false detections are modeled as the Poisson RFS, and the intensity function is given as κ ( j ) [31]. Hence, the number of false detections received at sensor j is Poisson distributed, with mean λ ( j ) . The measurement set collected by sensor j at time k is denoted as Z k ( j ) . We assume that multiple sensors can simultaneously observe the targets at time k. This set of sensors is called the set of activated sensors, which is denoted by A k { 1 , , P } . The set of measurements from all activated sensors is denoted as Z k = { Z k ( j ) , , Z k ( P ) } .

3. Generalized Labeled Multi-Bernoulli Recursion

Since the number of targets is not known, the set of targets, called the multi-target state, can be modeled as an RFS. By introducing the label issue, the labeled RFS has been proposed to incorporate target identity and perform tracking [31,32]. In this paper, we use a labeled RFS named GLMB, which is conjugate with respect to the multi-target likelihood function and is closed under the multi-target Chapman–Kolmogorov equation with respect to the multi-target transition kernel [31,32]. These properties make it possible to provide analytic solutions to the multi-target filtering problems.
Each target is augmented with a label L , where L denotes a discrete countable space. For clarity, lowercase letters are used to represent single-target states, e.g., x, x , while uppercase letters are used to represent multi-target states, e.g., X, X . To distinguish the labeled states from unlabeled ones, symbols for the former are bolded, e.g., x , X , etc. The multi-target state X is a labeled RFS given as
X = { x 1 , , x N } = { ( x 1 , 1 ) , , ( x N , N ) } F ( X × L ) ,
where X and N denote the state space and the number of single-target states, respectively.
As a labeled RFS on X L , the GLMB [31,32] is important that it is closed under the Chapman–Kolmogorov equation and is also a conjugate prior [31]. The GLMB distributed according to [31,32]
π ( X ) = ( X ) c C ω ( c ) ( L ( X ) ) [ p ( c ) ] X ,
where ( X ) δ X ( | L X | ) denote the distinct label indicator, C is a discrete index set, p ( c ) is a probability density that p ( c ) ( x , ) d x = 1 , each ω ( c ) satisfies L L c C ω ( c ) ( L ) = 1 , and L : X × L L denotes the projection L ( ( x , ) ) = . For numerical implementation of multi-target trackers, it is more convenient to use the so-called δ -GLMB form [32], which is a special case of GLMB with
C = F ( L ) × Ξ ,
ω ( c ) = ω ( I , ξ ) ( L ) = ω ( I , ξ ) δ I ( L ) ,
p ( c ) = p ( I , ξ ) = p ( ξ ) ,
where Ξ is a discrete space. Therefore, the distribution of the δ -GLMB is as follows:
π ( X ) = ( X ) ( I , ξ ) F L × Ξ ω ( I , ξ ) ( L ( X ) ) [ p ( ξ ) ] X .
If the multi-target prior is a δ -GLMB, its prediction π + ( X + ) is also a δ -GLMB with the form
π + ( X + ) = ( X + ) ( I + , ξ ) F L + × Ξ ω + ( I , ξ ) × δ I + ( L ( X + ) ) [ p + ( ξ ) ] X + ,
where
ω + ( I + , ξ ) = ω B I + B ω S ( ξ ) I + L ,
p + ( ξ ) x , = 1 L p S ( ξ ) x , + 1 1 L p B x , ,
p S ( ξ ) x , = p S · , f x | · , , p ( ξ ) · , η S ( ξ ) ,
η S ( ξ ) = p S · , f x | · , , p ( ξ ) · , d x ,
ω S ( ξ ) ( L ) = η S ( ξ ) L I L 1 I L q S ξ I L ω I , ξ ,
q S ( ξ ) ( ) = q S · , , p ξ · , .
In the above equations, p S · , is the survival probability, q S x , = 1 p S x , denotes the probability of death, p S ( ξ ) · , is the density of a surviving target, p B · , is the density for the new birth target, f · | · , is the transition density, ω B is the weight of birth labels, ω B ξ is the weight of surviving labels, f , g = f ( x ) g ( x ) d x is the standard inner product of functions f and g, and 1 S X is a inclusion function, which is equal to 1, if X S and 0, otherwise. The multi-target posterior is also a δ -GLMB with the form
π ( X | Z ) = X ( I + , ξ ) F L + × Ξ θ Θ ω + ( I , ξ ) × δ I L X p ( ξ , θ ) ( · | Z ) X ,
where Θ is the space of mappings θ : L { 0 , 1 , , | Z | } such that θ ( i ) = θ ( i ) implies i = i , and
ω ( I , ξ , θ ) ( Z ) ω ( I , ξ ) [ η Z ( ξ , θ ) ] I ,
p ( ξ , θ ) ( x , | Z ) = p ( ξ ) ( x , ) ψ Z ( x , ; θ ) η Z ( ξ , θ ) ,
η Z ( ξ , θ ) = p ( ξ ) ( · , ) , ψ Z ( · , ; θ ) ,
ψ Z ( x , ; θ ) = δ 0 ( θ ( ) ) q D ( x , ) + ( 1 δ 0 ( θ ( ) ) ) p D ( x , ) g ( z θ ( ) | x , ) κ ( z θ ( ) ) ,
where p D ( · , ) is the detection probability, q D ( · , ) = 1 p D ( · , ) denotes the probability of miss detection, κ ( · ) is the clutter intensity, g ( · | · , ) is the likelihood, and δ S ( X ) is a Kronecker delta function, which is equal to 1, if X = S and 0, otherwise.
The aforementioned recursion considers the propagation of GLMB for a single sensor. For the multi-sensor case, a feasible fusion approach is the iterated-corrector strategy [28] which approximates the multi-sensor update through iterated single sensor updates. Compared with the single sensor case, the multi-sensor case provides more observations about the tracked target, resulting in better tracking performance. Although the iterated-corrector strategy has no rigorous mathematical derivation, it is still an attractive multi-sensor fusion approach because of its low computational complexity and simple implementation.

4. Multi-Sensor GLMB Filter with CAR Based Birth Model

4.1. Problem Formulation

The implementation of the standard GLMB filter is introduced in Section 3. Obviously, the filter relies on the prior knowledge of birth distributions which limits its practical applications. In Section 1, it is analyzed that measurement driven methods provide reasonable solutions for estimating birth distributions, since measurements contain observable information about targets. However, Doppler measurements are significantly less informative, which makes the problem much more challenging.
The Doppler measurement (5) is a nonlinear function of target position and velocity. We have a scalar measurement and four unknowns in the components of the target state. Therefore, many 4-dimensional target states actually exist x = [ p x , p ˙ x , p y , p ˙ y ] T that are compatible with a single Doppler measurement, even without considering the noise [12]. To visualize the high sensitivity of the Doppler measurement (5) regarding the target position, Figure 1 shows the contour lines for the Doppler measurement in the ( p x , p y ) plane, for constant velocity p ˙ x = p ˙ y = 1 m/s of a target, and the sensor position is s x = s y = 0 m with the transmitting frequency f c = 900 MHz. Figure 1 shows that a large number of different target positions can produce similar Doppler shift values.
The aforementioned problem makes target state initialization a critical issue. Successful Doppler-only MTT has either used the A-R approach [48,49,50] or assumed an accurate prior knowledge of the birth distributions [5,38,42,44]. To refine the state estimates of birth targets while tracking multiple targets simultaneously, we go a step further and propose a novel multi-sensor GLMB tracking solution. Such an approach relies on: (i) an efficient measurement classification; (ii) a smart birth distribution estimation technique; and (iii) an iterated-corrector scheme for multi-sensor information fusion.

4.2. Measurement Classification

This paper attempts to generate birth distributions adaptively using the received measurements. In order to reduce errors of the generated birth distributions caused by measurement uncertainties, the birth distribution cannot be generated by measurements originating from persistent tracks. Hence, it is necessary to divide the received measurements into two sets, i.e., the set of persistent measurements and its complementary set. This is similar to the approach proposed in [58], which developed adaptive birth models within SMC-PHD and SMC-CPHD filters. However, since there is more information available in the GLMB filter, the classification of measurements is more straightforward.
As analyzed in Section 3, the discrete space Ξ is used to estimate the multi-target posterior δ -GLMB. Ξ is actually the space of association map histories, i.e., Θ 0 : k = Θ 0 × × Θ k , where Θ t is the association map space at time t. At each time step, the association map records the relation between measurements and targets, i.e., undetected targets are assigned 0 and a target with label associated with a measurement z θ ( ) is assigned θ ( ) . Hence, using the association map, it is easily to obtain the measurement set corresponding to the persistent tracks. At time k, given that the entire measurement set is Z k and the persistent measurement set is Z s , k , the complementary set Z ˜ s , k is obtained by eliminating Z s , k from Z as
Z ˜ s , k = Z k / Z s , k .
Note that a measurement associated with an existing target can not originate from a new birth target based on the assumption that a measurement is produced by at most one target. Therefore, the complementary set Z ˜ s , k at current time has not been associated with any of persistent tracks, and hence it may from new targets. The elements of Z ˜ s , k are used to estimate the birth distribution π B ( X ) , which is then propagated to the next time step.
The birth distribution π B ( X ) is a δ -GLMB RFS [31,32]:
π B ( X ) = ( X ) L F B ω B ( L ) × δ L ( L ( X ) ) [ p B ] X ,
where L is the birth label set with weight [31,32]
ω B ( L ) = B ( 1 r B ( ) ) L 1 B ( ) r B ( ) 1 r B ( ) .
The birth δ -GLMB can be completely characterized by
r B ( ) ( z ) , p B ( ) ( · ; z ) : = B ( j ) ( z ) z Z ˜ s , k , j J ,
where r B ( ) ( z ) is the existence probability, p B ( ) ( · ; z ) is the probability density, B ( j ) ( z ) is the label assigned for the birth target initiated by z from sensor j, and J is the set of sensor numbers. To distinguish the birth distributions generated by different sensors, the label is augmented as B ( j ) = k B ( j ) , i B ( j ) , j , where k B ( j ) is the birth time and i B ( j ) is the index that distinguishes targets born at the same time.
For a measurement z Z ˜ s , k , the new-born likelihood is given by [47]
r U ( z ) = 1 ( I , ξ ) F L × Ξ θ Θ ( I ) 1 z θ ( z ) ω ( I , ξ , θ ) ,
where the inclusion function ensures that the sum is performed only on updated hypotheses that assign the measurement to a target and the weight ω ( I , ξ , θ ) is given in (21). The probability of existence for the next time for the birth distribution initiated by measurement z Z ˜ s , k depends on the new-born likelihood obtained at current time [47]:
r B ( ) ( z ) = min r B max , r U ( z ) ζ Z ˜ s , k r U ( ζ ) · λ B ,
where r B m a x [ 0 , 1 ] is the maximum probability of existence for the new birth target, and λ B denotes the expected number of target births the next time.

4.3. CAR Based Birth Distribution Estimation

In the SMC formulation of the GLMB filter, the birth probability density p B ( ) ( · ; z ) generated by the measurement z Z ˜ s , k is approximated as
p B ( ) ( · ; z ) = i = 1 M b 1 M b δ x z ( i ) ( x ) , z Z ˜ s , k
where M b is the number of samples for the birth target, and x z ( i ) is the state of the ith sample. To obtain x z ( i ) = [ p x , p ˙ x , p y , p ˙ y ] T in (31), it is assumed that the measurement noise is irrelevant since Doppler measurements are typically accurate (small σ ε ).
Referring to (6) and dropping the time index k and sensor index j, the unit vector along the radar line-of-sight (RLOS) is defined by
u : = p s p s = p s ρ ,
where ρ is the range of the target from the sensor. We define the relative position coordinates ( x r , y r )
x r : = p x s x , y r : = p y s y .
Then, the range of the target from sensor is given by
ρ = x r 2 + y r 2 .
The radial velocity is defined as the projection of the target velocity along the RLOS [59] (see Figure 2) and is given by
v r : = p ˙ T u = v cos α , 0 α 2 π
where v = p ˙ = p ˙ x 2 + p ˙ y 2 is the speed of the target, and α is the angle between the target velocity and the RLOS. Then, the Doppler shift is expressed as
z = C v r = C v cos α ,
where C 2 f c / c . The radial velocity can also be expressed as
v r = p ˙ x x r + p ˙ y y r ρ .
Next, we define physics-based constraints in the unobservable subspace of the state space to reduce the space of admissible solutions of x z ( i ) in (31).
For a tracking scenario, the minimum and maximum values of the range ( ρ min , ρ max ) and maximum possible target speed ( v max ) are usually available from the prior information [59,60]. If we assume that the random variable (RV) ρ is uniformly distributed in the interval [ ρ min , ρ max ] , then the mean and SD are given by
ρ ¯ = ( ρ min + ρ max ) / 2
σ ρ = ( ρ max ρ min ) / 12
The uniform distribution has a finite support. The mean and SD are commonly used in a prior Gaussian prior distribution, which has an infinite support. If we use these parameters in a Gaussian distribution, in some cases, the samples drawn from the Gaussian distribution may lie outside the interval and may be negative. To avoid such problems, we use the mean ρ ¯ and assume that the interval contains the RV with a high probability, e.g., 0.999. Then, σ ρ , 0.999 is given by
σ ρ , 0.999 = ( ρ max ρ min ) / 6.581 .
Therefore, the distribution of the range ρ is denoted as N ( ρ ; ρ ¯ , σ ρ , 0.999 2 ) .
The maximum value of cos α in (31) is 1. Therefore, from (36), for a given measurement z, the minimum value of speed is
v min = | z / C | .
Given the maximum possible target speed ( v max ) and the minimum value of speed, a similar argument applies to the RV speed, and the distribution of the target speed v is denoted as N ( v ; v ¯ , σ v , 0.999 2 ) .
Let θ (see Figure 2) be the target heading (azimuth angle of target velocity measured from th X-axis in the counter-clockwise direction). Since θ has a circular distribution, we assume that θ is uniformly distributed in [ 0 , 2 π ) . We draw a sample for speed from N ( v ; v ¯ , σ v , 0.999 2 ) and a sample for θ . Then, the X and Y components of velocity are given by
p ˙ x = v cos θ , p ˙ y = v sin θ .
Once we sample the distance ρ form N ( ρ ; ρ ¯ , σ ρ , 0.999 2 ) and sample the target velocity p ˙ = [ p ˙ x , p ˙ y ] T using (42), a finite number of feasible solutions for the target position p = [ p x , p y ] T can be obtained using the following proposition.
Proposition 1.
Given the range ρ, the target velocity [ p ˙ x , p ˙ y ] T , and the target-originated Doppler measurement z at the sensor located at s = [ s x , s y ] T , there are two solutions for the target position
p x = z ρ p ˙ x ± I p ˙ x p ˙ y ρ C 2 v 2 z 2 C v 2 + s x p y = z ρ p ˙ y ± p ˙ x ρ C 2 v 2 z 2 C v 2 + s y
where I p ˙ x is an indicator of the sign of p ˙ x (positive and negative). If p ˙ x 0 , I p ˙ x = 1 and if p ˙ x < 0 , I p ˙ x = 1 . The first and second pair of solutions for p x and p y use the positive and negative signs, respectively, in (43). The proof of (43) is given in the Appendix A.
For each measurement z Z ˜ s , k , the above procedure will be repeated until M b samples representing new-born states have been obtained. Example 1 below illustrates the generation of the new-born target states.
Example 1.
We assume that the transmitting frequency of a sensor s located at [ 0 , 0 ] T is f c = 900 MHz , the maximum detection range is ρ max = 25,000 m , the minimum detection range is ρ min = 500 m , and the maximum speed of the target is v max = 35 m / s . The SD of measurement noise is equal to σ ε = 1 Hz . At the current time, a 4-dimensional target with state x = [ 3000 m , 15 m / s , 3500 m , 10 m / s ] T generates a Doppler measurement z = 13.0158 Hz at sensor s. To compare the performance of the CAR and A-R methods, both methods are used to generate M b = 5000 samples representing new-born states x z ( i ) = [ p x , p ˙ x , p y , p ˙ y ] T ( i = 1 , 2 , , 5000 ) based on the measurement z.
For the CAR method, the distribution of range ρ is N ( ρ ; 14,000 m , ( 3342.95 m ) 2 ) (obtained using (38) and (40)) and the distribution of speed v is N ( v ; 18.58 m / s , ( 4.99 m / s ) 2 ) . For the A-R method, we draw samples from a multivariate Gaussian N ( x ¯ ; μ x , C x ) , where μ x = [ 0 m , 0 m / s , 0 m , 0 m / s ] T and C x = diag [ ( 15,000 m ) 2 , ( 50 m / s ) 2 , ( 15,000 m ) 2 , ( 50 m / s ) 2 ] T . Only those samples whose likelihood larger than a threshold η = 1 × 10 10 are accepted in the A-R method. The samples generated by the A-R and CAR methods are illustrated in Figure 3 and Figure 4, respectively. Comparing Figure 3 with Figure 4, it is clear that the CAR method reduces the velocity space and position space significantly.
The Doppler measurements of the samples are shown in Figure 5. It can be observed from Figure 5b that the measurement range of samples is about [ 6.51 , 19.52 ] Hz for the A-R method. Meanwhile, for the CAR method, the measurements of all samples are equal to z = 13.02 Hz . The computing time for the A-R and CAR methods are 1.31 s and 0.05 s, respectively. The computing time of the CAR method is about 26.20 times faster than the A-R method. For the A-R method, a total number of 291,204 samples are generated but only 1.72% samples are accepted. Of course, decreasing the threshold η can improve the acceptance rate of the A-R method and reduce the computing time, but will further increase the uncertainty of the samples.

4.4. Detailed Implementation Steps

The samples generated in Section 4.3 are not only from the birth targets but also include “ghost samples”. Those “ghost samples” may occupy a large part of total samples because of the high nonlinearity of the measurement model. To reduce the false samples in one iteration, we use the iterated-corrector strategy for fusion information from multi-sensor. The iterated-corrector update may also keep “ghost samples” in one iteration, while these samples will be eliminated at the following time steps. The flow chart of the proposed CAR-GLMB approach is shown in Figure 6.
Algorithm 1 shows the pseudocode of a single run k = 2 , 3 , for the proposed approach. The following parameters are always assumed available to the MTT system:
  • sensor model parameters: position s ( j ) = [ s x , s y ] T , carrier frequency f c , j , detection probability p D ( j ) ( · ) , and clutter intensity κ ( · ) for sensor j ( j = 1 , , N s ) ;
  • birth model parameters: the maximum speed v max , the minimum ρ min and maximum ρ max detection ranges of the sensor, and the number M b of birth samples;
  • likelihood function g ( z | x , ) and transition density f ( x | · , ) ;
  • survival probability function p S ( x , ) .
Algorithm 1 Step-by-step pseudocode for a single run of the CAR-GLMB method.
INPUT: → The posterior π ( X | Z ) and birth distribution π B from previous time step
OUTPUT: → The posterior π ( X | Z ) and birth distribution π B at the current time
  1:  Predict persistent and birth distributions to obtain the prior distribution π + = π S π B
  2:  Select P sensors: A k = selec sensor ( P , π + )
  3:  Collect measurement set Z k ( A k ) from sensors A k
  4:  for sensor i A k  do
  5:        if  i = 1  then
  6:             Predict π + i = π i 1 ( X | Z )
  7:       else
  8:            Pseudo-predict π + i = π i 1 ( X | Z )
  9:        end if
 10:       Calculate parameters of the posterior to obtain π i ( X | Z )
 11:  end for
 12:  The fused posterior π ( X | Z ) = π i ( X | Z )
 13:  Collect the complementary measurement set Z ˜ s , k ( A k , ) from sensors A k
 14:  for z Z ˜ s , k ( A k , ) do
 15:         n 0
 16:        while  n < M b  do
 17:             Generate a birth sample using (38)–(43).
 18:              n n + 1 .
 19:        end while
 20:        Collect the birth distribution for z as { r B ( ) ( z ) , p B ( ) ( · , · ; z ) : = B ( z ) }
 21:  end for
 22:  The birth distribution is obtained as π B = ( X + ) L F ( B ) ω B ( L ) × δ L ( L ( X + ) ) [ p B ] X +
Lines 1 of Algorithm 1 implement (13)–(19). Since the number of Doppler sensors is generally large, a subset of sensors can be selected in line 2, to ensure that the communication and real-time constraints are met. The sensor management problem is beyond the scope of this work; please refer to [61,62,63] for more information. The selected sensors in line 3 supply their measurement set Z k ( A k ) and lines 4–12 implement the iterated-corrector scheme. Finally, the generation of birth distribution is implemented in lines 13–22.

5. Numerical Simulations and Results

Numerical results are presented for 2D tracking scenarios where targets with the time-varying and unknown number are observed by Doppler-only sensors, as shown in Figure 7. The carrier frequency f c and measurement time interval T of each sensor are 900 MHz and 10 s. The probability of detection P D and measurement noise SD of each sensor are 0.95 and 1 Hz. The measurement space is Z = [ 200 , 200 ] Hz . The false detections are uniformly distributed over the measurement space with the Poisson rate λ = 2 per scan. The minimum and maximum detection ranges of a sensor are 500 m and 25,000 m, and the maximum speed of a target is 35 m/s. For the adaptive birth model, we choose r B max = 0.05 and λ B = 0.3 [47]. The survival probability is p S = 0.99 [31,32]. The A-R method (regarded as the state-of-the art method) is widely used to generate the birth distribution using Doppler-only measurements and is used as a comparative algorithm. Parameters used in the CAR and A-R methods are the same as the ones used in Example 1.
We use the OSPA [64] and OSPA ( 2 ) [37] distances to measure the tracking error which are widely used in the field of RFS based MTT. The OSPA estimates’ errors in the cardinality and localization by measuring the distance between two sets of states. Let d p ( c ) ( ϕ , ψ ) be the OSPA between two sets ϕ , ψ F ( X ) with order p and cutoff c. For ϕ = { ϕ ( 1 ) , ϕ ( 2 ) , , ϕ ( m ) } and ψ = { ψ ( 1 ) , ψ ( 2 ) , , ψ ( n ) } with m n , the OSPA d p ( c ) ( ϕ , ψ ) is measured as follows [64]:
d p ( c ) ( ϕ , ψ ) 1 n min π Π n i = 1 m d ¯ ( c ) ϕ ( i ) , ψ ( π ( i ) ) p + c p n m 1 p ,
where d ¯ ( c ) ϕ ( i ) , ψ ( π ( i ) ) = min c , ϕ ( i ) ψ ( π ( i ) ) and π is a permutation function in the set Π n of all the permutations of ψ . If m > n , then d p ( c ) ( ϕ , ψ ) = d p ( c ) ( ψ , ϕ ) . However, the OSPA error does not consider the track labeling issue [37]. The OSPA ( 2 ) [37] is an adaptation of the OSPA which contains the interpretation of a per-track per-time error and can accommodate sets of tracks. Let X = { x ( 1 ) , x ( 2 ) , , x ( m ) } and Y = { y ( 1 ) , y ( 2 ) , , y ( n ) } denote two sets of tracks, where m n . The OSPA ( 2 ) d ˜ p ( c ) ( X , Y ) between X and Y is defined as follows [37]:
d ˜ p ( c ) ( X , Y ) 1 n min π Π n i = 1 m d ˜ ( c ) x ( i ) , y ( π ( i ) ) p + c p n m 1 p ,
where d ˜ ( c ) · , · is the base-distance, which is the time average OSPA given by
d ˜ p ( c ) ( x , y ) t D x D y d ( c ) ( { x ( t ) } , { y ( t ) } ) | D x D y | , D x D y 0 , D x D y =
where D x and D y denote the time domain. If m > n , then d ˜ p ( c ) ( X , Y ) = d ˜ p ( c ) ( Y , X ) . All experiments are tested in Matlab R2010a and implemented on a computer with a 3.40 GHz processor.

5.1. Scenario 1: Single Target Simulation

We first consider a scenario for tracking a single target, which exists for time indices k = 1 , 2 , , 30 with the initial state vector x = [ 3000 m , 15 m / s , 3500 m , 10 m / s ] T . The SD of the process noise is σ w = 0.1 m / s 2 . The purpose of this scenario is to clearly observe the change of the particles representing the posterior density over time, since particles are the focus of this paper. To verify the performance of the proposed CAR based tracking method, we compare it with the A-R based tracking method.
We use a simple sensor selection approach, in which sensors are randomly selected at each measurement time. For the GLMB recursion parameters, the number of components calculated and stored in each forward propagation is 100. When the number of sensors selected at each measurement time k is two, the positions of particles representing the posterior density at different time steps are illustrated in Figure 8. It can be observed from Figure 8 that the distribution of particles is relatively scattered at the initial measurement times. Hence, the initial estimated target position is less accurate. As time progresses, the distribution of particles becomes concentrated, and the estimated target position becomes closer to the true position.
Figure 9 plots the true and CAR based estimated x and y coordinates versus time of the true trajectory and estimates of the tacking method. These results are consistent with results in Figure 8. We observe that the tracker takes some time to detect the target. However, the target’s position estimates become more accurate over time.
Next, we analyze the average tracking performance of the CAR and A-R methods. To fully study the influence of the threshold in the A-R method, three specific different thresholds η = 0 , η = 1 × 10 10 , and η = 1 × 10 3 are used. Note that using the threshold η = 0 means that the samples are generated over the entire surveillance space and all the generated samples are accepted. To have the average performance, we have performed 60 Monte Carlo runs. When two sensors are selected at each measurement time, the averaged OSPA distance ( p = 1 and c = 10,000 m) and OSPA ( 2 ) distance (with the same c, p, and window length w = 10 ) are shown in Figure 10a,b, respectively. These figures show that the CAR method is more accurate than the A-R method with different thresholds in terms of the OSPA and OSPA ( 2 ) distances. In addition, the OSPA and OSPA ( 2 ) errors of the A-R method decrease as the threshold η increases from 0 to 1 × 10 3 . A larger threshold indicates that the measurement of the sample is closer to the received measurement, i.e., the sample is more accurate. Therefore, the A-R method with a large threshold provides better tracking accuracy. However, when the threshold is large, the acceptance rate of samples is low, and the A-R method requires a long computation time to obtain a specific number of samples.
To compare the computation time and tracking accuracy of CAR and A-R methods with different number of selected sensors, the averaged computation time for executing a complete MC simulation, and the averaged OSPA and OSPA ( 2 ) errors averaging over all measurement times, are shown in Table 1. The relative computation time is determined relative to the corresponding CAR based computation time. The results demonstrate that, selecting more sensors at each measurement time, the values of the OSPA and OSPA ( 2 ) distances decrease. As the number of sensors increases, the computation time increases accordingly. When the same number of sensors is selected at each measurement time, the CAR method always performs better than the A-R method in the OSPA and OSPA ( 2 ) errors. When the number of sensors is 1, 2, and 3, our OSPA error drops by 8.10–28.92%, 30.17–44.02%, and 32.82–58.77% respectively, compared to the A-R method. Our OSPA ( 2 ) error drops by 5.29–18.92%, 20.44–32.77%, and 21.35–43.51%, respectively. For each sensor, the birth distribution generated by our approach is more accurate than that of the A-R method. This advantage becomes more significant as more sensors are used, as seen by the reduction percentage in error. When the threshold η = 0 , the computation time of the CAR method is comparable with that of the A-R method. This means that the CAR method takes little time to generate samples. When the threshold η = 1 × 10 10 , the CAR method runs about 3.74–4.92 times faster than the A-R method. When the threshold η = 1 × 10 3 , the CAR method runs about 6.87–9.56 times faster than the A-R method. These results are consistent with the theoretical analysis, showing the superiority of the CAR method with regard to the state-of-the-art A-R method.

5.2. Scenario 2: Multi-Target Simulation

In this scenario, all targets in the surveillance area travel with the NCV motion but with different velocities. A challenging tracking scenario is considered where targets move close to each other and their initial states are [ 3000 m , 15 m / s , 3500 m , 10 m / s ] T , [ 3000 m , 25 m / s , 2000 m , 5 m / s ] T , and [ 4500 m , 30 m / s , 5500 m , 0 m / s ] T . The SD of the process noise is σ w = 0.1 m / s 2 .
For the GLMB recursion parameters, the number of components calculated and stored in each forward propagation is 500. Two sensors are selected at each measurement time. Figure 11 plots the x and y coordinates versus time for the true and CAR estimated trajectories. These plots show that the CAR method is able to detect and track targets successfully.
The averaged OSPA distance ( p = 1 and c = 10,000 m) and OSPA ( 2 ) distance (with the same c, p, and window length w = 10 ) are shown in Figure 12a,b, respectively. These figures show that the CAR method is more accurate than the A-R method with different thresholds in terms of the OSPA and OSPA ( 2 ) distance.
To compare the computational efficiencies and tracking accuracies of the CAR and A-R methods as a function of number of selected sensors, the average computation time for executing a complete MC simulation, and the average OSPA and OSPA ( 2 ) errors are shown in Table 2. Consistent with the results in Scenario 1, the values of the OSPA and OSPA ( 2 ) distances are lower when more sensors are selected at each measurement time. When the number of sensors is 1, 2, and 3, our OSPA error drops by 21.90–27.03%, 23.78–44.30%, and 29.19–62.80% respectively, compared to the A-R method. Our OSPA ( 2 ) error drops by 14.07–20.40%, 15.15–27.92%, and 17.31–43.65%, respectively. There is an overall increase of the computation time with an increase in the number of sensors. The CAR method always outperforms the A-R method in terms of the OSPA and OSPA ( 2 ) distances, when the same number of sensors is selected. When the threshold is η = 1 × 10 10 , the CAR method runs about 1.57–1.66 times faster than the A-R method. When the threshold is η = 1 × 10 3 , the CAR method runs about 2.40–2.61 times faster than the A-R method. Compared with Scenario 1, Scenario 2 uses more components in the GLMB recursion to accurately track multiple targets.

6. Conclusions

MTT using Doppler-only measurements is a challenging problem in the multi-target multi-sensor tracking system due to the poor observability of the target state. In this paper, we have developed a novel CAR based birth model that adaptively generates the birth distribution using the received Doppler-only measurements. The CAR based model uses physics-based constraints in the unobservable subspace of the state space to reduce ambiguity in the initial state. We have implemented the CAR based birth model in the multi-sensor GLMB filter to demonstrate its effectiveness. Our simulation results show that the novel birth model works significantly faster than the approach used by the state-of-the-art methods, with superior tracking accuracy. Future work will investigate other types of measurements, such as Bearing-only measurements, for tracking the distinguishable and non-cooperative targets.

Author Contributions

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

Funding

This work was supported in part by the National Natural Science Foundation of China under Grant No. 62007022, the Natural Science Foundation of Shaanxi Province under Grant No. 2021JQ-209, and the Fundamental Research Funds for the Central Universities under Grant Nos. GK202103082 and JB210210.

Data Availability Statement

In this work, we have used the free RFS MATLAB code provided by Prof. Ba-Ngu Vo and Prof. Ba-Tuong Vo at http://ba-tuong.vo-au.com/codes.html, accessed on 4 May 2022.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof of Proposition 1.
We have samples of the target velocity [ p ˙ x , p ˙ y ] T and range of the target ρ , and a Doppler measurement z by a sensor located at [ s x s y ] T . Our aim is to obtain solutions for the relative position coordinates ( x r , y r ) . Once the relative position coordinates are known, then the absolute position coordinates are given by
p x = x r + s x , p y = y r + s y .
Ignoring the measurement noise, the Doppler measurement and range are given, respectively, by
z = C ρ ( p ˙ x x r + p ˙ y y r ) ,
ρ = x r 2 + y r 2 .
From (A2) we obtain
x r = 1 p ˙ x ( z ρ C + p ˙ y y r ) .
We can write (A4) as
x r = a 1 + a 2 y r ,
where
a 1 = z ρ p ˙ x C , a 2 = p ˙ y p ˙ x .
Squaring (A3), we obtain
ρ 2 = x r 2 + y r 2 .
Substituting the expression for x r from (A5) in (A7) and simplifying, we obtain the quadratic equation for y r
b 1 y r 2 + 2 b 2 y r + b 3 = 0 ,
where
b 1 = ( 1 + a 2 2 ) , b 2 = a 1 a 2 , b 3 = a 1 2 a 2 2 .
After simplifying, we obtain the solutions of the quadratic Equation (A8) as
y r = z ρ p ˙ y ± p ˙ x ρ C 2 v 2 z 2 C v 2 .
Substituting the expression of y r from (A10) in (A5) and simplifying, we obtain the two corresponding solutions for x r as
x r = z ρ p ˙ x ± I p ˙ x p ˙ y ρ C 2 v 2 z 2 C v 2 ,
where I p ˙ x is an indicator of the sign of p ˙ x (positive and negative). Then, the absolute position coordinates are determined by (A1) to obtain
p x = z ρ p ˙ x ± I p ˙ x p ˙ y ρ C 2 v 2 z 2 C v 2 + s x ,
p y = z ρ p ˙ y ± p ˙ x ρ C 2 v 2 z 2 C v 2 + s y .
Hence, we have two solution pairs for ( p x , p y ) in (A12) and (A13), and Proposition 1 is proved. □

References

  1. Salinger, S.N.; Brandstatter, J.J. Application of recursive estimation and Kalman filtering to Doppler tracking. IEEE Trans. Aerosp. Electron. Syst. 1970, 4, 585–592. [Google Scholar] [CrossRef]
  2. Chan, Y.T.; Jardine, F.L. Target localization and tracking from Doppler-shift measurements. IEEE J. Ocean. Eng. 1990, 15, 251–257. [Google Scholar] [CrossRef]
  3. Richter, Y.; Gerasimov, J.; Balal, N.; Pinhasi, Y. Tracking of evasive objects using bistatic Doppler radar operating in the millimeter wave regime. Remote Sens. 2022, 14, 867. [Google Scholar] [CrossRef]
  4. Sun, W.; Ji, M.; Huang, W.; Ji, Y.; Dai, Y. Vessel tracking using bistatic compact HFSWR. Remote Sens. 2020, 12, 1266. [Google Scholar] [CrossRef] [Green Version]
  5. Subedi, S.; Zhang, Y.D.; Amin, M.G.; Himed, B. Group sparsity based multi-target tracking in passive multi-static radar systems using Doppler-only measurements. IEEE Trans. Signal Process. 2016, 64, 3619–3634. [Google Scholar] [CrossRef]
  6. Ristic, B.; Farina, A. Recursive Bayesian state estimation from Doppler-shift measurements. In Proceedings of the International Conference on Intelligent Sensors, Sensor Networks and Information Processing, Adelaide, SA, Australia, 6–9 December 2011; pp. 538–543. [Google Scholar]
  7. Borisov, A.; Bosov, A.; Miller, B.; Miller, G. Passive Underwater Target Tracking: Conditionally Minimax Nonlinear Filtering with Bearing-Doppler Observations. Sensors 2020, 20, 2257. [Google Scholar] [CrossRef]
  8. Miller, A.; Miller, B.; Miller, G. Navigation of Underwater Drones and Integration of Acoustic Sensing with Onboard Inertial Navigation System. Drones 2021, 5, 83. [Google Scholar] [CrossRef]
  9. Asaf, T.; Itzik, K.; Reuven, K. Inertial Navigation System/Doppler Velocity Log (INS/DVL) Fusion with Partial DVL Measurements. Sensors 2017, 17, 415. [Google Scholar]
  10. Shames, I.; Bishop, A.N.; Smith, M.; Anderson, B.D.O. Doppler shift target localization. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 266–276. [Google Scholar] [CrossRef] [Green Version]
  11. Battistelli, G.; Chisci, L.; Fantacci, C.; Farina, A.; Graziano, A. A new approach for Doppler-only target tracking. In Proceedings of the 16th International Conference on Information Fusion, Istanbul, Turkey, 9–12 July 2013; pp. 1616–1623. [Google Scholar]
  12. Bosse, J.; Krasnov, O.; Yarovoy, A. Direct target localization and deghosting in active radar network. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 3139–3150. [Google Scholar] [CrossRef]
  13. Mallick, M.; Krishnamurthy, V.; Vo, B.N. Integrated Tracking, Classification, and Sensor Management: Theory and Applications; Wiley: Hoboken, NJ, USA; IEEE: Piscataway, NJ, USA, 2012. [Google Scholar]
  14. Fortmann, T.E.; Bar-Shalom, Y.; Scheffe, M. Sonar tracking of multiple targets using joint probabilistic data association. IEEE J. Ocean. Eng. 1983, 8, 173–184. [Google Scholar] [CrossRef] [Green Version]
  15. Reid, D. An algorithm for tracking multiple targets. IEEE Trans. Autom. Control 1979, 24, 843–854. [Google Scholar] [CrossRef]
  16. Kurien, T. Issues in the design of practical multitarget tracking algorithms. In Multitarget-Multisensor Tracking: Advanced Applications; Bar-Shalom, Y., Ed.; Artech House: Norwood, MA, USA, 1990; Chapter 3; pp. 43–83. [Google Scholar]
  17. Roecker, J.A. Multiple scan joint probabilistic data association. IEEE Trans. Aerosp. Electron. Syst. 1995, 31, 1204–1210. [Google Scholar] [CrossRef]
  18. Pattipati, K.; Deb, S.; Bar-Shalom, Y.; Washburn, R. A new relaxation algorithm and passive sensor data association. IEEE Trans. Autom. Control 1992, 37, 198–213. [Google Scholar] [CrossRef]
  19. Poore, A.; Rijavec, N. A Lagrangian relaxation algorithm for multidimensional assignment problems arising from multitarget tracking. SIAM J. Optim. 1993, 3, 544–563. [Google Scholar] [CrossRef]
  20. Fitzgerald, R.J. Track Biases and Coalescence with Probabilistic Data Association. IEEE Trans. Aerosp. Electron. Syst. 1985, AES-21, 822–825. [Google Scholar] [CrossRef]
  21. Blom, H.; Bloem, E.A. Probabilistic data association avoiding track coalescence. IEEE Trans. Autom. Control 2000, 45, 247–259. [Google Scholar] [CrossRef]
  22. Zhu, Y.; Liang, S.; Wu, X.; Yang, H. A random finite set based joint probabilistic data association filter with non-homogeneous Markov chain. Front. Inf. Technol. Electron. Eng. 2021, 22, 1114–1126. [Google Scholar] [CrossRef]
  23. Zhu, Y.; Wang, J.; Liang, S.; Wang, J. Covariance control joint integrated probabilistic data association filter for multi-target tracking. IET Radar Sonar Navig. 2019, 13, 584–592. [Google Scholar] [CrossRef]
  24. Zhu, Y.; Wang, J.; Liang, S. Efficient joint probabilistic data association filter based on Kullback–Leibler divergence for multi-target tracking. IET Radar Sonar Navig. 2017, 11, 1540–1548. [Google Scholar] [CrossRef]
  25. Mahler, R. Statistical Multisource-Multitarget Information Fusion; Artech House: Norwood, MA, USA, 2007. [Google Scholar]
  26. Mahler, R. Advances in Statistical Multisource-Multitarget Information Fusion; Artech House: Norwood, MA, USA, 2014. [Google Scholar]
  27. Vo, B.N.; Mallick, M.; Bar-Shalom, Y.; Coraluppi, S.; Osborne, R., III; Mahler, R.; Vo, B.T. Multitarget tracking. In Wiley Encyclopedia of Electrical and Electronics Engineering; John Wiley & Sons: Hoboken, NJ, USA, 2015; pp. 1–25. [Google Scholar]
  28. Mahler, R. Multitarget Bayes filtering via first-order multitarget moments. IEEE Trans. Aerosp. Electron. Syst. 2003, 39, 1152–1178. [Google Scholar] [CrossRef]
  29. Mahler, R. PHD filters of higher order in target number. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 1523–1543. [Google Scholar] [CrossRef]
  30. Vo, B.T.; Vo, B.N.; Cantoni, A. The cardinality balanced multi-target multi-Bernoulli filter and its implementations. IEEE Trans. Signal Process. 2009, 57, 409–423. [Google Scholar]
  31. Vo, B.T.; Vo, B.N. Labeled random finite sets and multi-target conjugate priors. IEEE Trans. Signal Process. 2013, 61, 3460–3475. [Google Scholar] [CrossRef]
  32. Vo, B.N.; Vo, B.T.; Phung, D. Labeled random finite sets and the Bayes multi-target tracking filter. IEEE Trans. Signal Process. 2014, 62, 6554–6567. [Google Scholar] [CrossRef] [Green Version]
  33. Vo, B.T.; Vo, B.N. Multi-scan generalized labeled multi-Bernoulli filter. In Proceedings of the 2018 21st International Conference on Information Fusion (FUSION), Cambridge, UK, 10–13 July 2018; pp. 195–202. [Google Scholar]
  34. Vo, B.N.; Vo, B.T. A multi-scan labeled random finite set model for multi-object state estimation. IEEE Trans. Signal Process. 2019, 67, 4948–4963. [Google Scholar] [CrossRef] [Green Version]
  35. Reuter, S.; Vo, B.T.; Vo, B.N.; Dietmayer, K. The labeled multi-Bernoulli filter. IEEE Trans. Signal Process. 2014, 62, 3246–3260. [Google Scholar]
  36. Vo, B.N.; Vo, B.T.; Hoang, H.G. An efficient implementation of the generalized labeled multi-Bernoulli filter. IEEE Trans. Signal Process. 2017, 65, 1975–1987. [Google Scholar] [CrossRef] [Green Version]
  37. Beard, M.; Vo, B.T.; Vo, B.N. A solution for large-scale multi-object tracking. IEEE Trans. Signal Process. 2020, 68, 2754–2769. [Google Scholar] [CrossRef] [Green Version]
  38. Ma, L.; Kim, D.Y.; Xue, K. Multi-Bernoulli filter for target tracking with multi-static Doppler only measurement. Signal Process. 2014, 108, 102–110. [Google Scholar]
  39. Mahler, R. Statistics 101 for multisensor, multitarget data fusion. IEEE Aerosp. Electron. Syst. Mag. 2004, 19, 53–64. [Google Scholar] [CrossRef]
  40. Mahler, R. “Statistics 102” for Multisource-Multitarget Detection and Tracking. IEEE J. Sel. Top. Signal Process. 2013, 7, 376–389. [Google Scholar] [CrossRef]
  41. Ristic, B. Particle Filters for Random Set Models; Springer: Berlin, Germany, 2013. [Google Scholar]
  42. Papi, F. Multi-sensor d-GLMB filter for multi-target tracking using Doppler only measurements. In Proceedings of the 2015 European Intelligence and Security Informatics Conference, Manchester, UK, 7–9 September 2015; pp. 83–89. [Google Scholar]
  43. Zhu, Y.; Zhao, L.; Zhang, Y.; Wu, X. Receiver selection for multi-target tracking in multi-static Doppler radar systems. EURASIP J. Adv. Signal Process. 2021, 2021, 118. [Google Scholar] [CrossRef]
  44. Do, C.T.; Nguyen, H.V. Tracking multiple targets from multistatic Doppler radar with unknown probability of detection. Sensors 2019, 19, 1672. [Google Scholar] [CrossRef] [Green Version]
  45. Vo, B.N.; Vo, B.T.; Beard, M. Multi-sensor multi-object tracking with the generalized labeled multi-Bernoulli filter. IEEE Trans. Signal Process. 2019, 67, 5952–5967. [Google Scholar] [CrossRef] [Green Version]
  46. Reuter, S.; Meissner, D.; Wilking, B.; Dietmayer, K. Cardinality balanced multi-target multi-Bernoulli filtering using adaptive birth distributions. In Proceedings of the 16th International Conference on Information Fusion Istanbul, Istanbul, Turkey, 9–12 July 2013; pp. 1608–1615. [Google Scholar]
  47. Lin, S.F.; Vo, B.T.; Nordholm, S.E. Measurement driven birth model for the generalized labeled multi-Bernoulli filter. In Proceedings of the International Conference on Control, Automation and Information Sciences, Ansan, Korea, 27–29 October 2016; pp. 94–99. [Google Scholar]
  48. Guldogan, M.B. Consensus Bernoulli filter for distributed detection and tracking using multi-static Doppler shifts. IEEE Signal Process. Lett. 2014, 21, 672–676. [Google Scholar] [CrossRef]
  49. Ristic, B.; Farina, A. Target tracking via multi-static Doppler shifts. IET Radar Sonar Navig. 2011, 7, 508–516. [Google Scholar] [CrossRef]
  50. Guldogan, M.B.; Lindgren, D.; Gustafsson, F.; Habberstad, H. Multi-target tracking with PHD filter using Doppler-only measurements. Digit. Signal Process. 2014, 27, 1–11. [Google Scholar] [CrossRef] [Green Version]
  51. Jones, B.A. PHD filter birth modeling using the probabilistic admissible region. IEEE Trans. Aerosp. Electron. Syst. 2018, 54, 1456–1469. [Google Scholar] [CrossRef]
  52. Fujimoto, K.; Scheeres, D.J.; Herzog, J.; Schildknecht, T. Association of optical tracklets from a geosynchronous belt survey via the direct Bayesian admissible region approach. Adv. Space Res. 2014, 53, 295–308. [Google Scholar] [CrossRef] [Green Version]
  53. Bar-Shalom, Y. Multitarget-Multisensor Tracking: Advanced Applications; Artech House: Norwood, MA, USA, 1990. [Google Scholar]
  54. Blackman, S.S. Multiple-Target Tracking with Radar Applications; Artech House: Norwood, MA, USA, 1986. [Google Scholar]
  55. Bar-Shalom, Y.; Li, X.R.; Kirubarajan, T. Estimation with Applications to Tracking and Navigation: Theory, Algorithms and Software; John Wiley and Sons, Inc.: New York, NY, USA, 2001; pp. 268–274. [Google Scholar]
  56. Horn, R.A.; Johnson, C.R. Matrix Analysis; Cambridge University Press: Cambridge, UK, 1985. [Google Scholar]
  57. Richards, M.A. Fundamentals of Radar Signal Processing; McGraw-Hill: New York, NY, USA, 2005. [Google Scholar]
  58. Ristic, B.; Clark, D.; Vo, B.N.; Vo, B.T. Adaptive target birth intensity for PHD and CPHD filters. IEEE Trans. Aerosp. Electron. Syst. 2012, 48, 1656–1668. [Google Scholar] [CrossRef]
  59. Mallick, M.; Bar-Shalom, Y.; Kirubarajan, T.; Moreland, M. An improved single-point track initiation using GMTI measurements. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 2697–2714. [Google Scholar] [CrossRef]
  60. Ristic, B.; Arulampalam, S.; Gordon, N. Beyond the Kalman Filter-Particle Filters for Tracking Applications; Artech House: Norwood, MA, USA, 2004. [Google Scholar]
  61. Beard, M.; Vo, B.T.; Vo, B.N.; Arulampalam, S. Void probabilities and Cauchy-Schwarz divergence for generalized labeled multi-Bernoulli models. IEEE Trans. Signal Process. 2017, 65, 5047–5061. [Google Scholar] [CrossRef] [Green Version]
  62. Panicker, S.; Gostar, A.K.; Bab-Hadiashar, A.; Hoseinnezhad, R. Tracking of targets of interest using labeled multi-Bernoulli filter with multi-sensor control. Signal Process. 2020, 171, 107451. [Google Scholar] [CrossRef]
  63. Zhu, Y.; Liang, S.; Gong, M.; Yan, J. Decomposed POMDP optimization-based sensor management for multi-target tracking in passive multi-sensor systems. IEEE Sens. J. 2022, 22, 3565–3578. [Google Scholar] [CrossRef]
  64. 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]
Figure 1. Contour plot of the Doppler measurement with a constant velocity p ˙ x = p ˙ y = 1 m/s and the sensor located at position s x = s y = 0 m.
Figure 1. Contour plot of the Doppler measurement with a constant velocity p ˙ x = p ˙ y = 1 m/s and the sensor located at position s x = s y = 0 m.
Remotesensing 14 03131 g001
Figure 2. Definition of radial velocity.
Figure 2. Definition of radial velocity.
Remotesensing 14 03131 g002
Figure 3. Results of the A-R method: (a,b) sample positions; (c,d) sample velocities.
Figure 3. Results of the A-R method: (a,b) sample positions; (c,d) sample velocities.
Remotesensing 14 03131 g003
Figure 4. Results of the CAR method: (a,b) sample positions; (c,d) sample velocities.
Figure 4. Results of the CAR method: (a,b) sample positions; (c,d) sample velocities.
Remotesensing 14 03131 g004
Figure 5. Doppler measurements of the samples: (a) the CAR method; (b) the A-R method.
Figure 5. Doppler measurements of the samples: (a) the CAR method; (b) the A-R method.
Remotesensing 14 03131 g005
Figure 6. Flow chart of the CAR-GLMB approach.
Figure 6. Flow chart of the CAR-GLMB approach.
Remotesensing 14 03131 g006
Figure 7. The locations of Doppler sensors.
Figure 7. The locations of Doppler sensors.
Remotesensing 14 03131 g007
Figure 8. The positions of particles representing the posterior density at different measurement times, and the real and estimated positions of the target are represented by red stars and blue squares, respectively: (a) k = 2 ; (b) k = 3 ; (c) k = 6 ; (d) k = 10 ; (e) k = 20 ; (f) k = 30 .
Figure 8. The positions of particles representing the posterior density at different measurement times, and the real and estimated positions of the target are represented by red stars and blue squares, respectively: (a) k = 2 ; (b) k = 3 ; (c) k = 6 ; (d) k = 10 ; (e) k = 20 ; (f) k = 30 .
Remotesensing 14 03131 g008
Figure 9. True and estimated x and y coordinates versus time for a single target.
Figure 9. True and estimated x and y coordinates versus time for a single target.
Remotesensing 14 03131 g009
Figure 10. Performance comparison in Scenario 1: (a) OSPA distance; (b) OSPA ( 2 ) distance.
Figure 10. Performance comparison in Scenario 1: (a) OSPA distance; (b) OSPA ( 2 ) distance.
Remotesensing 14 03131 g010
Figure 11. Scenario 2: True and CAR estimated x and y coordinates versus time.
Figure 11. Scenario 2: True and CAR estimated x and y coordinates versus time.
Remotesensing 14 03131 g011
Figure 12. Performance comparison in multi-target Scenario 2: (a) OSPA distance; (b) OSPA ( 2 ) distance.
Figure 12. Performance comparison in multi-target Scenario 2: (a) OSPA distance; (b) OSPA ( 2 ) distance.
Remotesensing 14 03131 g012
Table 1. Average performance comparison in Scenario 1.
Table 1. Average performance comparison in Scenario 1.
MethodNumber of SensorsComputing Time (s)Relative Computing TimeOSPA Error (m)OSPA ( 2 ) Error (m)
A-R
(threshold = 0)
15.341.025818.916883.98
220.231.073832.675375.01
336.541.073488.994758.76
A-R
(threshold = 1 × 10 10 )
120.053.844991.906257.66
270.843.743496.934934.66
3167.754.922360.853668.70
A-R
(threshold = 1 1 × 10 3 )
138.277.334500.245893.79
2130.316.873072.374542.32
3325.879.562141.503417.82
CAR15.221.004135.935581.78
218.961.002145.423613.66
334.081.001438.662688.03
Table 2. Average performance comparison in Scenario 2.
Table 2. Average performance comparison in Scenario 2.
MethodNumber of SensorsComputing Time (s)Relative Computing TimeOSPA Error (m)OSPA ( 2 ) Error (m)
A-R
(threshold = 0)
181.851.065949.456697.29
2206.511.044729.476533.19
3301.631.024184.666116.06
A-R
(threshold = 1 × 10 10 )
1127.711.665824.856600.18
2311.411.574188.666078.71
3484.161.632975.354950.81
A-R
(threshold = 1 × 10 3 )
1185.022.415558.526204.03
2475.902.403456.055550.20
3773.322.612198.644167.35
CAR176.921.004341.225331.12
2198.321.002634.284709.37
3296.511.001556.833446.18
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhu, Y.; Mallick, M.; Liang, S.; Yan, J. Generalized Labeled Multi-Bernoulli Multi-Target Tracking with Doppler-Only Measurements. Remote Sens. 2022, 14, 3131. https://doi.org/10.3390/rs14133131

AMA Style

Zhu Y, Mallick M, Liang S, Yan J. Generalized Labeled Multi-Bernoulli Multi-Target Tracking with Doppler-Only Measurements. Remote Sensing. 2022; 14(13):3131. https://doi.org/10.3390/rs14133131

Chicago/Turabian Style

Zhu, Yun, Mahendra Mallick, Shuang Liang, and Junkun Yan. 2022. "Generalized Labeled Multi-Bernoulli Multi-Target Tracking with Doppler-Only Measurements" Remote Sensing 14, no. 13: 3131. https://doi.org/10.3390/rs14133131

APA Style

Zhu, Y., Mallick, M., Liang, S., & Yan, J. (2022). Generalized Labeled Multi-Bernoulli Multi-Target Tracking with Doppler-Only Measurements. Remote Sensing, 14(13), 3131. https://doi.org/10.3390/rs14133131

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