Next Article in Journal
Cloud Based IoT Network Virtualization for Supporting Dynamic Connectivity among Connected Devices
Next Article in Special Issue
Elimination of Motion-Induced Phase Based on Double-Time Switching Scheme for SAA FMCW Radar
Previous Article in Journal
Reduction of Induction Motor Energy Consumption via Variable Velocity and Flux References
Previous Article in Special Issue
915-MHz Continuous-Wave Doppler Radar Sensor for Detection of Vital Signs
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multisensor RFS Filters for Unknown and Changing Detection Probability

1
Electronics & Information Engineering, Beihang University, Beijing 100191, China
2
Department of Engineering, University of Cambridge, Cambridge CB12PZ, UK
*
Author to whom correspondence should be addressed.
Electronics 2019, 8(7), 741; https://doi.org/10.3390/electronics8070741
Submission received: 16 May 2019 / Revised: 19 June 2019 / Accepted: 28 June 2019 / Published: 30 June 2019
(This article belongs to the Special Issue Radar Sensor for Motion Sensing and Automobile)

Abstract

:
The detection probability is an important parameter in multisensor multitarget tracking. The existing multisensor multi-Bernoulli (MS-MeMBer) filter and multisensor cardinalized probability hypothesis density (MS-CPHD) filter require that detection probability is a priori. However, in reality, the value of the detection probability is constantly changing due to the influence of sensors, targets, and other environmental characteristics. Therefore, to alleviate the performance deterioration caused by the mismatch of the detection probability, this paper applies the inverse gamma Gaussian mixture (IGGM) distribution to both the MS-MeMBer filter and the MS-CPHD filter. Specifically, the feature used for detection is assumed to obey the inverse gamma distribution and is statistically independent of the target’s spatial position. The feature is then integrated into the target state to iteratively estimate the target detection probability as well as the motion state. The experimental results demonstrate that the proposed methods can achieve a better filtering performance in scenarios with unknown and changing detection probability. It is also shown that the distribution of the sensors has a vital influence on the filtering accuracy, and the filters perform better when sensors are dispersed in the monitoring area.

1. Introduction

The objective of multitarget tracking is to jointly estimate the number and state of multiple targets from sensor observations. Recently, the tracking algorithm based on random finite set (RFS) has gained increasing attention in the field of multitarget tracking [1,2,3]. Unlike traditional multitarget tracking algorithms, such as multiple hypotheses tracking (MHT) [4] and joint probabilistic data association (JPDA) [5], this method directly avoids the data association step. Specifically, by treating target states and sensor measurements as RFSs, the target tracking problem is transformed into a set-valued estimation problem under the Bayesian framework. On the basis of the finite set statistics (FISST) [6], a probability hypothesis density (PHD) filter has been proposed in [7]. The PHD filter uses the moment approximation and propagates the first-order moment of the multitarget probability density, which is assumed to be Poisson-distributed. However, the cardinality variance estimated by the PHD filter grows with the number of the targets. As an improvement to the PHD filter, the cardinalized PHD (CPHD) filter [8] jointly propagates the first-order moment and the cardinality distribution. Unlike the above RFS filters based on the moment approximation, the multitarget multi-Bernoulli (MeMBer) filter proposed in [6] directly approximates the multitarget posterior with multi-Bernoulli components, and thus more density information is preserved in the Bayesian recursion. However, the first-order Tayler linear approximation in the updating process introduces a cardinality estimation bias [9], and thus a cardinalized MeMBer (CBMeMBer) filter is proposed in [10] to eliminate this bias. The limit of the multi-Bernoulli filter is that it can only accommodate target states. Therefore, the δ-generalized labeled multi-Bernoulli (δ-GLMB) filter which can output target tracks has been developed [11], and two approximations of the δ-GLMB filter, labeled the multi-Bernoulli (LMB) filter and the marginalized δ-GLMB filter, have been designed to reduce the computational cost [12,13]. Nevertheless, when the track estimation is not needed, the multi-Bernoulli filter and the CPHD/PHD filters have been shown to be more efficient compared with the labeled RFS filters; thus, they have been widely used [14,15,16].
A number of RFS tracking methods have also been studied in the multisensor field [17,18,19,20,21,22]. Mahler [17] derived a generalized PHD (G-PHD) filter for the case of two sensors, and Delande et al. [18] extended it to scenarios with an arbitrary number of sensors. To avoid the combinatorial computational complexity of the G-PHD filter, the iterated corrector PHD (IC-PHD) filter was proposed in [19]. The parallel-combination approximate multisensory (PCAM) PHD filter [20] is a theoretical approximation of multisensor PHD filters and is independent of the sensor order. In [21], the multisensor CPHD (MS-CPHD) filter was introduced, and the multisensor MeMBer (MS-MeMBer) filter was proposed in [22]. Based on the labeled RFS, the centralized multisensor δ-GLMB (MS-δ-GLMB) with extended association maps was introduced in [23].
The detection probability is usually considered as a known parameter in general tracking algorithms. Whereas in real scenarios, the detection probability is affected by sensors, targets, the environment, and the features used for detection [24]. For instance, Mahler et al. [25] proposed CPHD/PHD filters that can accommodate parameter mismatch in the clutter rate and detection profile in which the detection probability is obtained by the Beta Gaussian mixture (BGM) model. Subsequently, a multitarget multi-Bernoulli filter [26,27] was proposed to adaptively learn the nonhomogeneous clutter intensity and detection probability. Meanwhile, a multitarget tracker based on the labeled RFS was proposed in [28], which can estimate the clutter rate and detection profile. All these filters assume that the detection probability is Beta-distributed and calculate the detection probability through the accumulation of observation effects in the BGM model. However, the filtering performance deteriorates when the detection probability is low, and poor model initialization parameters can adversely affect the filtering performance. For this reason, the inverse gamma Gaussian mixture (IGGM) model is used in CPHD/PHD filters [29]. Specifically, the feature used for detection is an inverse gamma distribution, and by combining the feature into the target motion state, the performance loss caused by the mismatch of the detection probability parameters is reduced.
In the field of multisensor and multitarget tracking, the MS-CPHD filter and MS-MeMBer filter are two important methods, and the ways they are generalized to handle scenarios with an unknown and time varying detection of probability are similar. Therefore, in this paper, we propose IGGM-MS-CPHD/MeMBer filters for the unknown and varying detection probability case. Initially, we incorporated the detection feature into the target motion state. On the basis of the augmented state, the IGGM model was applied to both the MS-CPHD filter and the MS-MeMBer filter to iteratively estimate the detection probability and motion parameters. Then the estimated detection probability was updated to the next step for further filtering. The simulation results show that compared to GM-based MS-CPHD/MeMBer filters, the proposed filters achieve a better performance when the detection probability is unknown and dynamically changing and when sensors are dispersed in the observation area.
This paper is organized as follows: Section 2 gives a review of the MS-CPHD/MeMBer filters and the IGGM model, the proposed IGGM-MS-CPHD/MeMBer filters are derived in Section 3, the simulation results are presented in Section 4, and the conclusions are given in Section 5.

2. Background

In this section, we present the basic materials for the MS-CPHD/MeMBer filters and the IGGM model.

2.1. Multisensor Measurement Partitioning

Measurement partitioning means that the measurement set is divided into a finite number of mutually disjoint subsets. The MS-CPHD/MeMBer filters have high computational complexity since all measurement partitioning forms are required for calculating the accurate solution. To improve the efficiency, Nannuru et al. [21] proposed a greedy measurement partitioning algorithm that obtains a finite number of partitioning hypotheses to derive an approximate solution of the MS-CPHD filter. Subsequently, Sancan et al. applied this partitioning algorithm to the MS-MeMBer filter in [22]. The objective of the greedy measurement partitioning algorithm is to recombine the multisensor observations and obtain several measurement partitioning hypotheses. Assume that Z k 1 : s = { Z k 1 , Z k 2 ,…, Z k s } is the observation set generated by all the s sensors at time k, where Z k j is the measurement set generated by the jth sensor. Then by using the greedy measurement partitioning algorithm, Z k 1 : s can be reconstructed as Z k 1 : s = P = { W 1 : s 1 , W 1 : s 2 , , W 1 : s | P | 1 , V } , where P is the measurement partitioning hypothesis. W 1 : s j denotes the measurements generated by target j from all sensors at time k, V is the clutter measurements from all sensors, and the divided parts follow Equations (1)–(3):
W 1 : s i W 1 : s j = , 1 i j | P | 1 ,
W 1 : s i V = , 1 i | P | 1 ,
W 1 : s 1 W 1 : s 2 W 1 : s | P | 1 V = Z k 1 : s .
Define Q to be the set composed of the measurement partitioning hypothesis P, T W 1 : s j = { ( i , l ) | z l i W 1 : s j } is a mapping function, and z l i is the lth measurement in the observation set generated by sensor i.

2.2. The MS-CPHD Filter

The prediction process of the MS-CPHD filter is the same as that of the single-sensor CPHD filter in [8]. Define x as the state vector of the target. Assume that the first-order moment of the multitarget posterior at time k−1 is D k 1 | k 1 ( x ) and the cardinality distribution is p k 1 | k 1 ( n ) , then the predicted PHD function is given by
D k | k 1 ( x ) = ρ s v ( x ) f k | k 1 ( x | · ) , D k 1 | k 1 + b k ( x ) ,
and the predicted cardinality distribution follows
p k | k 1 ( n ) = j = 0 n p b , k ( n j ) l = j C j l p k 1 | k 1 ( l ) D k 1 | k 1 , ρ s v j D k 1 | k 1 , 1 ρ s v l j D k 1 | k 1 , 1 l ,
where b k ( x ) and p b , k ( n j ) are the first-order moment and cardinality distribution of the newborn targets at time k, respectively. f k | k 1 ( x | · ) is the transfer function. The combination coefficient follows C j l = l ! / ( j ! ( l j ) ! ) , and p s v ( x ) is the survival probability of the target with state x at time k.
The updated PHD equation of the MS-CPHD filter is
D k | k ( x ) r k | k 1 ( x ) = α 0 j = 1 s q d j ( x ) + P Q α P ( j = 1 | P | 1 ρ W 1 : s j ( x ) ) ,
where
α 0 = P Q ( K P M ( | P | ) ( γ ) j = 1 | P | 1 d W 1 : s j ) P Q ( K P M ( | P | 1 ) ( γ ) ) j = 1 | P | 1 d W 1 : s j ,
α P = K P M ( | P | 1 ) ( γ ) j = 1 | P | 1 d W 1 : s j U Q ( K U M ( | U | 1 ) ( γ ) j = 1 | U | 1 d W 1 : s j ) ,
ρ W 1 : s j ( x ) = ( ( i , l ) T W 1 : s j p d i ( x ) g i ( z l i | x ) ) ( i , ) T W 1 : s j q d i ( x ) r k | k 1 ( x ) ( ( i , l ) T W 1 : s j p d i ( x ) g i ( z l i | x ) ) ( i , ) T W 1 : s j q d i ( x ) d x .
The normalized function for the predicted PHD D k | k 1 ( x ) is r k | k 1 ( x ) . Assume α P as the weight of the measurement partitioning hypothesis P in the set Q, and the all-clutter partitioning hypothesis has an associated weight α 0 . d W 1 : s j denotes the weight of W 1 : s j in hypothesis P, and U is also a measurement partitioning hypothesis. p d i ( x ) is the detection probability of the target with state x by sensor i, and q d i ( x ) = 1 p d i ( x ) . g i ( z l i | x ) is the measurement likelihood of sensor i. The probability that the target with state x is not detected by any sensor is given by γ = r k | k 1 ( x ) j = 1 s q d j ( x ) d x . K P = i = 1 s C i ( m i | P | i ) ( 0 ) . C i ( n ) ( ) and M i ( n ) ( ) are the nth-order derivatives of the probability generating functions (PGFs) of the clutter cardinality distribution and the predicted cardinality distribution, respectively.
After the updating process, the posterior cardinality distribution is
p k | k ( n ) = p k | k 1 ( n ) P Q | P | n + 1 ( K P n ! ( n | P | + 1 ) ! γ n | P | + 1 j = 1 | P | 1 d W 1 : s j ) P Q K P M ( | P | 1 ) ( γ ) j = 1 | P | 1 d W 1 : s j .

2.3. The MS-MeMBer Filter

Since no observation information is required in the prediction process, the predicted density of the MS-MeMBer filter is the same as that of the MeMBer/CBMeMBer filters [9,10]. Suppose that at time k–1, the multitarget density is a multi-Bernoulli of the form
π k 1 | k 1 = { ( r k 1 | k 1 ( i ) , p k 1 | k 1 ( i ) ) } i = 1 M k 1 | k 1 ,
where M k 1 | k 1 is the component number. Then the predicted multitarget density is also a multi-Bernoulli and is formed by the union of surviving sets { ( r S , k | k 1 ( i ) , p S , k | k 1 ( i ) ) } i = 1 M k 1 | k 1 and newborn sets { ( r b , k ( i ) , p b , k ( i ) ) } i = 1 M b , k , i.e.,
{ ( r k | k 1 ( i ) , p k | k 1 ( i ) ) } i = 1 M k | k 1 = { ( r S , k | k 1 ( i ) , p S , k | k 1 ( i ) ) } i = 1 M k 1 | k 1 { ( r b , k ( i ) , p b , k ( i ) ) } i = 1 M b , k ,
where r ( i ) in Equation (12) is the existence probability of target i, and p ( i ) ( · ) is the corresponding density function. The surviving multi-Bernoulli components are given by Equations (13) and (14):
r S , k | k 1 ( i ) = r k 1 | k 1 ( i ) p k 1 | k 1 ( i ) , ρ s v ,
p S , k | k 1 ( i ) ( x ) = f k | k 1 ( x | · ) , p k 1 | k 1 ( i ) ρ s v p k 1 | k 1 ( i ) , ρ s v .
As the updated multitarget posterior density π k | k no longer has the original form, an approximate multi-Bernoulli distribution π ^ k | k = P Q j = 1 M k | k 1 { ( r k | k ( P , j ) , p k | k ( P , j ) ) } with an equal PHD function is derived in [16], where
r k | k ( P , j ) = { α P r k | k 1 ( j ) p k | k 1 ( j ) , γ 1 r k | k 1 ( j ) + r k | k 1 ( j ) p k | k 1 ( j ) , γ W 1 : s j =   α P W 1 : s j ,
p k | k ( P , j ) ( x ) = { p k | k 1 ( j ) ( x ) γ ( x ) p k | k 1 ( j ) , γ W 1 : s j = p k | k 1 ( j ) ( x ) f ( W 1 : s j | x ) p k | k 1 ( j ) ( x ) f ( W 1 : s j | x ) d x W 1 : s j .
In Equation (16), the probability that the target with state x is not detected by any sensor is given by γ ( x ) = i = 1 s ( 1 p d i ( x ) ) . Define f ( W 1 : s j | x ) to be the multisensor likelihood for the single target with state x:
f ( W 1 : s j | x ) ( i , l ) T W 1 : s j p d i ( x ) g i ( z l i | x ) c ( z l i ) ( i , ) T W 1 : s j ( 1 p d i ( x ) ) .
c ( · ) is the density function of the clutter. The weight of the measurement partitioning hypothesis P in set Q is
α P K P j = 1 M k | k 1 φ W 1 : s j j [ 1 ] U Q K U j = 1 M k | k 1 φ W 1 : s j j [ 1 ] ,
where
φ W 1 : s j j [ u ] { 1 r k | k 1 ( j ) + r k | k 1 ( j ) p k | k 1 ( j ) , u γ W 1 : s j = r k | k 1 ( j ) u ( x ) p k | k 1 ( j ) ( x ) f ( W 1 : s j | x ) d x W 1 : s j .

2.4. The IGGM Model

The CPHD/PHD filters for scenarios with an unknown detection probability were proposed in [18]. This method assumes that the detection probability is Beta-distributed, and the detection probability is calculated through the accumulation of observation effects. However, the filtering performance decreases significantly when the detection probability is low. In practice, the target detection probability depends on the features used for detection. For example, the gray-scale characteristics and structures of the infrared target [30] need to be considered, and, in radar detection, the signal-to-noise ratio (SNR) and the amplitude are crucial for determining the target detection probability [31,32,33]. Generally, these features are usually non-negative and non-Gaussian. In addition, the stronger the SNR and echo amplitude the higher the target detection probability. Therefore, the IGGM model and its analytical form of CPHD/PHD filters were proposed in [29]. In the IGGM model, it is assumed that the single-target state x contains the kinematic state x ˜ and also the feature a used for detection, i.e., x = [ x ˜ , a ] T . Then the measurement state can be denoted as z = [ z ˜ , h ] T , where z ˜ is the position measurement and h is the feature measurement. In a short time, rather than the target’s position, the SNR and echo amplitude are mainly dependent on the radar reflection cross section (RCS) and the target’s shape. Thus it can be assumed that the feature a is statistically independent of the kinematic state x ˜ , and the detection probability is determined by the feature a, i.e., p d , k ( x ˜ , a ) = p d , k ( a ) .
The probability density of the inverse gamma distribution I G ( x ; u , v ) is
I G ( x ; u , v ) = v u Γ ( u ) x u 1 exp ( v x ) ,
where variable x > 0, shape parameter u > 0, and scale parameter v > 0. Γ ( · ) denotes the gamma function. The mode of the inverse gamma distribution is v / u + 1 , and the mean value is v / u 1 .
Accordingly, the probability density of the gamma distribution follows
G ( x ; u , v ) = v u Γ ( u ) x u 1 exp ( v x ) ,
where x > 0, u > 0, and v > 0. In Equation (21), the mode and mean values are ( u 1 ) / v and u / v , respectively.
Assume that at time k–1, feature a obeys the inverse gamma distribution with
f k 1 | k 1 ( a ) = I G ( a ; u k 1 , v k 1 ) ,
the Markov transfer density L k | k 1 ( a | a ) with feature a is determined as follows:
f k | k 1 ( a ) = L k | k 1 ( a | a ) f k 1 | k 1 ( a ) d a ,
f k | k 1 ( a ) = I G ( a ; u k | k 1 , v k | k 1 ) ,
u k | k 1 = k u u k 1 ( 0 < k u < 1 ) ,
v k | k 1 = v k 1 u k 1 1 ( u k | k 1 1 ) ,
and the measurement likelihood g k ( h | a ) obeys the gamma distribution, i.e.,
g k ( h | a ) = G ( h ; ξ , ξ a ) .
Thus, feature a is also an inverse gamma distribution after the prediction and update step:
f k | k ( a ) A z ( h ) I G ( a ; u k | k , v k | k ) ,
where
A z ( h ) = v k | k 1 u k | k 1 Γ ( u k | k 1 ) ξ ξ Γ ( ξ ) Γ ( u k | k ) v k | k u k | k h ξ 1 ,
u k | k = ξ + u k | k 1 ,
v k | k = v k | k 1 + ξ h .
The larger the parameter ξ is the greater the influence the measurement h has on the estimation.
Assume that the single-target kinematic states x ˜ , the Markov transfer density L k | k 1 ( x ˜ | x ˜ ) with state x ˜ , and the measurement likelihood function g ˜ k ( z ˜ | x ˜ ) at time k–1 are all Gaussian,
f k 1 | k 1 ( x ˜ ) = N ( x ˜ ; m k 1 , P k 1 ) ,
L k | k 1 ( x ˜ | x ˜ ) = N ( x ˜ ; F k 1 x ˜ , Q k 1 ) ,
g ˜ k ( z ˜ | x ˜ ) = N ( z ˜ ; H k x ˜ , R k ) ,
where N ( · ; m , P ) denotes a Gaussian density with mean m and covariance P. F k 1 is the state transition matrix, Q k 1 is the process noise covariance, H k is the observation matrix, and R k denotes the observation noise covariance. Then the form of the kinematic state x ˜ remains unchanged after the prediction and update step:
f k 1 | k ( x ˜ ) = N ( x ˜ ; m k | k 1 , P k | k 1 ) ,
f k | k ( x ˜ ) q z ( z ˜ ) N ( x ˜ ; m k , P k ) ,
where
m k | k 1 = F k 1 m k 1 ,
P k | k 1 = F k 1 P k 1 F k 1 T + Q k 1 ,
q z ( z ˜ ) = N ( z ˜ ; m k | k 1 , R k + H k P k | k 1 H k T ) ,
m k ( z ˜ ) = m k | k 1 + K k ( z ˜ H k m k | k 1 ) ,
P k = P k | k 1 K k H k P k | k 1 ,
K k = P k | k 1 H k [ R k + H k P k | k 1 H k T ] 1 .

3. Multisensor Filters Based on IGGM Model

The prior density and the likelihood function of the feature a obey the inverse gamma distribution and gamma distribution, respectively, which are a pair of conjugate distributions. Thus, we can derive the analytical forms of the MS-CPHD/MeMBer filters which estimate the detection probability based on the IGGM model.

3.1. The IGGM-MS-CPHD Filter

3.1.1. Prediction

Assume that at time k–1 the posterior PHD has the IGGM form
D k 1 | k 1 ( x ) = i = 1 J k 1 | k 1 w k 1 | k 1 ( i ) N ( x ˜ ; m k 1 | k 1 ( i ) , P k 1 | k 1 ( i ) ) I G ( a ; u k 1 | k 1 ( i ) , v k 1 | k 1 ( i ) ) .
The PHD of the birth RFS at time k is
D b , k ( x ) = i = 1 J b , k w b , k ( i ) N ( x ˜ ; m b , k ( i ) , P b , k ( i ) ) I G ( a ; u b , k ( i ) , v b , k ( i ) ) ,
and the survival probability of the target with state x is constant, i.e., ρ s v ( x ) = ρ s v . w k 1 | k 1 ( i ) and w b , k ( i ) are the weights of the target state. J k 1 | k 1 and J b , k are component numbers. Then the predicted PHD also comprises IGGM with
D k | k 1 ( x ) = i = 1 J k 1 | k 1 p s v w k 1 | k 1 ( i ) N ( x ˜ ; m k | k 1 ( i ) , P k | k 1 ( i ) ) I G ( a ; u k | k 1 ( i ) , v k | k 1 ( i ) ) + i = 1 J b , k w b , k ( i ) N ( x ˜ ; m b , k ( i ) , P b , k ( i ) ) I G ( a ; u b , k ( i ) , v b , k ( i ) )
where u k | k 1 ( i ) , v k | k 1 ( i ) , m k | k 1 ( i ) , and P k | k 1 ( i ) are given in Equations (25), (26), (37), and (38).

3.1.2. Update

Suppose that at time k the predicted PHD of the multitarget is
D k | k 1 ( x ) = i = 1 J k | k 1 w k | k 1 ( i ) N ( x ˜ ; m k | k 1 ( i ) , P k | k 1 ( i ) ) I G ( a ; u k | k 1 ( i ) , v k | k 1 ( i ) ) ,
where J k | k 1 = J k 1 | k 1 + J b , k . Then the form of the multitarget PHD remains unchanged after the update step:
D k | k ( x ) = i = 1 J k | k 1 α 0 j = 1 s q d j ( a ) w ˜ k | k 1 ( i ) N ( x ˜ ; m k | k 1 ( i ) , P k | k 1 ( i ) ) I G ( a ; u k | k 1 ( i ) , v k | k 1 ( i ) ) + P Q j = 1 | P | 1 i = 1 J k | k 1 w k | k ( W 1 : s j , i ) N ( x ˜ ; m k | k ( W 1 : s j , i ) , P k | k ( W 1 : s j , i ) ) I G ( a ; u k | k ( W 1 : s j , i ) , v k | k ( W 1 : s j , i ) )
where the updated weight is
w k | k ( W 1 : s j , i ) = w ˜ k | k 1 ( i ) ( u , l ) T W 1 : s j p d u ( a ) q z ( z ˜ l u ) A z ( h l u ) ( v , ) T W 1 : s j q d v ( a ) .
w ˜ k | k 1 ( i ) = w k | k 1 ( i ) / j = 1 J k | k 1 w k | k 1 ( j ) gives the normalized weight. The updated Gaussian function and inverse gamma function are determined by Equations (49) and (50):
N ( x ˜ ; m k | k ( W 1 : s j , i ) , P k | k ( W 1 : s j , i ) ) N ( x ˜ ; m k | k 1 ( i ) , P k | k 1 ( i ) ) ( u , l ) T W 1 : s j g ˜ u ( z ˜ l u | x ˜ ) ,
I G ( a ; u k | k ( W 1 : s j , i ) , v k | k ( W 1 : s j , i ) ) I G ( a ; u k | k 1 ( i ) , v k | k 1 ( i ) ) ( u , l ) T W 1 : s j g u ( h l u | a ) .
Moreover, the updated cardinality is given in Equation (10) with
γ = n = 1 J k | k 1 w ˜ k | k 1 ( n ) j = 1 s ( 1 p d j ( a k | k 1 ( n ) ) ) ,
d W 1 : s j = ( n = 1 J k | k 1 w k | k 1 ( n ) ( i , l ) T W 1 : s j p d i ( a k | k 1 ( n ) ) q z ( z ˜ l i ) × A z ( h l i ) ( i , ) T W 1 : s j q d i ( a k | k 1 ( n ) ) ) ( i , l ) T W 1 : s j c ˜ ( z ˜ l i ) A z ( h l i ) ,
where, in Equations (51) and (52), the predicted feature a is estimated by the mean value, i.e., a k | k 1 ( n ) = v k | k 1 ( n ) / ( u k | k 1 ( n ) 1 ) . c ˜ ( · ) is the spatial density function of the clutter.

3.2. The IGGM-MS-MeMBer Filter

3.2.1. Prediction

Assume that at time k–1 the multitarget posterior has an IGGM form with J k 1 | k 1 ( i ) components:
p k 1 | k 1 ( i ) ( x ) = n = 1 J k 1 | k 1 ( i ) w n , k 1 | k 1 ( i ) N ( x ˜ ; m n , k 1 | k 1 ( i ) , P n , k 1 | k 1 ( i ) ) I G ( a ; u n , k 1 | k 1 ( i ) , v n , k 1 | k 1 ( i ) ) .
The survival probability of the target ρ s v , k ( x ) = ρ s v , k is constant, and the density function of the newborn target is
p b , k ( i ) ( x ) = n = 1 J b , k ( i ) w b , n , k ( i ) N ( x ˜ ; m b , n , k ( i ) , P b , n , k ( i ) ) I G ( a ; u b , n , k ( i ) , v b , n , k ( i ) ) ,
where w n , k 1 | k 1 ( i ) and w b , n , k ( i ) are the weights of the target state. Then the predicted density of each Bernoulli component also has the IGGM form, which can be obtained from Equations (13), (14), (53), and (54):
r S , k | k 1 ( i ) = r k 1 | k 1 ( i ) ρ s v ,
p S , k | k 1 ( i ) ( x ) = n = 1 J k 1 | k 1 ( i ) w n , k 1 | k 1 ( i ) N ( x ˜ ; m n , k | k 1 ( i ) , P n , k | k 1 ( i ) ) I G ( a ; u n , k | k 1 ( i ) , v n , k | k 1 ( i ) ) ,
where m n , k | k 1 ( i ) , P n , k | k 1 ( i ) , u n , k | k 1 ( i ) , and v n , k | k 1 ( i ) are respectively shown in Equations (25), (26), (37), and (38).

3.2.2. Update

Suppose that at time k the predicted density of each Bernoulli component is
p k | k 1 ( i ) ( x ) = n = 1 J k | k 1 ( i ) w n , k | k 1 ( i ) N ( x ˜ ; m n , k | k 1 ( i ) , P n , k | k 1 ( i ) ) I G ( a ; u n , k | k 1 ( i ) , v n , k | k 1 ( i ) ) ,
where J k | k 1 ( i ) = J k 1 | k 1 ( i ) + J b , k ( i ) . Then this IGGM form remains unchanged after the update step, with
p k | k ( P , i ) ( x ) = n = 1 J k | k 1 ( i ) w n , k | k ( P , i ) N ( x ˜ ; m n , k | k ( P , i ) , P n , k | k ( P , i ) ) I G ( a ; u n , k | k ( P , i ) , v n , k | k ( P , i ) ) ,
r k | k ( P , i ) = { α P r k | k 1 ( i ) n = 1 J k | k 1 ( i ) γ ( a n , k | k 1 ( i ) ) w n , k | k 1 ( i ) 1 r k | k 1 ( i ) + r k | k 1 ( i ) n = 1 J k | k 1 ( i ) γ ( a n , k | k 1 ( i ) ) w n , k | k 1 ( i ) W 1 : s i =   α P W 1 : s i ,
where the Gaussian component, the inverse gamma component, and the updated weight after the update process are given as follows:
N ( x ˜ ; m n , k ( P , i ) , P n , k ( P , i ) ) { N ( x ˜ ; m n , k | k 1 ( i ) , P n , k | k 1 ( i ) ) W 1 : s i = N ( x ˜ ; m n , k | k 1 ( i ) , P n , k | k 1 ( i ) ) ( j , l ) T W 1 : s i g ˜ j , k ( z ˜ l j | x ˜ ) W 1 : s i ,
I G ( a ; u n , k ( P , i ) , v n , k ( P , i ) ) { I G ( a ; u n , k | k 1 ( i ) , v n , k | k 1 ( i ) ) W 1 : s i = I G ( a ; u n , k | k 1 ( i ) , v n , k | k 1 ( i ) ) ( j , l ) T W 1 : s i g j , k ( h l j | a ) W 1 : s i ,
w n , k | k ( P , i ) { γ ( a n , k | k 1 ( i ) ) w n , k | k 1 ( i ) m = 1 J k | k 1 ( i ) γ ( a n , k | k 1 ( i ) ) w m , k | k 1 ( i ) W 1 : s i = ( w n , k | k 1 ( i ) ( j , l ) T W 1 : s i p d j ( a n , k | k 1 ( i ) ) q z ( z ˜ l j ) A z ( h l j ) c ( z l j ) · ( j , ) T W 1 : s i ( 1 p d j ( a n , k | k 1 ( i ) ) ) ) ( m = 1 J k | k 1 ( i ) w m , k | k 1 ( i ) ( j , l ) T W 1 : s i p d j ( a m , k | k 1 ( i ) ) q z ( z ˜ l j ) A z ( h l j ) c ( z l j ) · ( j , ) T W 1 : s i ( 1 p d j ( a m , k | k 1 ( i ) ) ) ) W 1 : s i .
a n , k | k 1 ( i ) = u n , k | k 1 ( i ) / ( v n , k | k 1 ( i ) 1 ) is the mean value of the feature a, and the clutter density function is c ( z ) = c ˜ ( z ˜ ) A z ( h ) . The detection probability of sensor j is determined by
p d j ( a ) = { ε 1 ( exp ( ( a F t h ) / δ 1 ) ε 2 ) a < F t h ε 1 ( 2 exp ( ( a F t h ) / δ 2 ) ε 2 ) a F t h ,
where F t h is the gate, ε 1   =   ( 2 exp ( F t h / δ 1 ) ) 1 , and ε 2 = exp ( F t h / δ 1 ) .

4. Simulation Results and Analysis

In this section, we explain the two simulation experiments we carried out to verify the proposed IGGM-MS-CPHD/MeMBer filters. First, the proposed filters were compared with the GM-MS-CPHD/MeMBer filters. Then, we analyzed the performance of the proposed algorithms under different sensor distributions. The optimal subpattern assignment (OSPA) distance [34] was used as the error metric:
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 order p = 1 and truncation threshold c = 100. X = { x 1 , , x m } and Y = { y 1 , , y n } . d ¯ ( c ) ( x , y ) = min ( d ( x , y ) , c ) denote the distance between vector x and y cutting off at threshold c. Пn is the set of permutations on {1, 2, …, n}.

4.1. System Model

Figure 1 shows a two-dimensional scene with six targets in a uniform linear motion. Figure 2 shows the detection probability of targets and their detection feature a. In this paper, the amplitude of the target echo signal was taken as the detection feature. Each feature obeyed the sinusoidal distribution with a period of 40 s. The observation process lasted for 100 frames and the sampling period was Δt = 1 s.
The target state x contains the kinematics state x ˜ and the detection feature a, where x ˜ k = [ p x , k , p ˙ x , k , p y , k , p ˙ y , k ] T denotes the position and velocity. The measurement z also consists of the position measurement z ˜ and the feature measurement h, i.e., z = [ z ˜ , h ] T , and z ˜ = [ r k , θ k ] T is the range and bearing information. In Equation (33), the state transition matrix and the process noise covariance matrix are respectively
F k 1 = [ 1 Δ t 0 0 0 1 0 0 0 0 1 Δ t 0 0 0 1 ] Q k 1 = [ Δ t 3 3 Δ t 2 2 0 0 Δ t 2 2 Δ t 0 0 0 0 Δ t 3 3 Δ t 2 2 0 0 Δ t 2 2 Δ t ] σ v 2 ,
where the standard deviation of the process noise is σ v = 8 m/s2.
For the nonlinear observation model, the extended Kalman filter (EKF) method [35] was used. The observation matrix and the measurement noise covariance matrix of sensor i are given as follows:
H k i = [ p ˜ x , k | k 1 p ˜ x , k | k 1 2 + p ˜ y , k | k 1 2 0 p ˜ y , k | k 1 p ˜ x , k | k 1 2 + p ˜ y , k | k 1 2 0 p ˜ y , k | k 1 p ˜ x , k | k 1 2 + p ˜ y , k | k 1 2 0 p ˜ x , k | k 1 p ˜ x , k | k 1 2 + p ˜ y , k | k 1 2 0 ] , R k = [ σ r 2 0 0 σ θ 2 ] .
The measurement noise in range and bearing are σ r = 30 m and σ θ = 1°, respectively. p ˜ x , k | k 1 = p x , k | k 1 + τ x i and p ˜ y , k | k 1 = p y , k | k 1 + τ y i , where [ τ x i , τ y i ] T is the coordinate offset of the sensor i with respect to the fusion center.
In the simulation, the target survival probability ρ s v , k was set to 0.99. As shown in Table 1, the birth model was set according to the target initial states. For the detection feature a in the prediction process, we set the transfer parameter ku to 0.7, and the parameter in measurement likelihood was set to ξ = 3. In Equation (63), we set F t h to 5.5, δ 1 to 4, and δ 2 to 2. We assumed that the clutter measurements generated by the sensors were Poisson with the uniform spatial distribution c ˜ ( z ˜ ) and mean clutter rate λ c , i.e., κ k ( z ) = λ c c ˜ ( z ˜ ) I G ( a ; u k κ , v k κ ) g k ( h | a ) d a .
For efficiency purposes, the pruning and merging methods in [10,29] were adopted in this paper. For the greedy measurement partitioning algorithm, the maximum number of the measurement subsets Wmax was set to 4, and the maximum number of partitioning hypotheses Pmax was set to 4. In the GM-MS-CPHD/MeMBer filters, the detection probability p d was fixed at 0.9.

4.2. Simulation 1

This section verifies the performance of the proposed IGGM-MS-CPHD/MeMBer filters in scenarios where the detection probability is unknown and dynamically changing. As shown in Figure 1, the sensors were located at [0 km, 0 km], [10 km, 10 km], [−10 km, 10 km], [−10 km, −10 km], and [10 km, −10 km], respectively. We took Sensor 1 as the fusion center, and the scopes of the sensors are shown in Table 2.
The average clutter intensity λ c was set to 10, with shape parameter u k κ set to 31 and scale parameter v k κ set to 280. Figure 3 shows the observation results of the sensors, in which the black dots represent the measurements.
Specifically, we compared the traditional GM-MS-CPHD/MeMBer filters with the IGGM-MS-CPHD/MeMBer filters derived in this paper. The average OSPA and the cardinality estimation of the four algorithms with 100 Monte Carlo runs are shown in Figure 4 and Figure 5. It can be seen that at the beginning of the filtering process the four filters achieve similar filtering performance. The detection probability of Target 1 was lower than 0.4 at about 30 s, while the parameter was fixed at pd = 0.9. As the detection probability decreases the filtering performance of the GM-MS-CPHD/MeMBer filters begins to deteriorate with a cardinality estimation bias. In contrast, the proposed IGGM-MS-CPHD/MeMBer filters have lower OSPA distances and obtain a smaller deviation of the estimated cardinality.
Table 3 shows the average single-step running time of the four algorithms with 100 Monte Carlo simulations. It can be seen that the proposed IGGM-MS-CPHD/MeMBer filters have a higher computational cost, and the IGGM-MS-CPHD filter demands the longest time.
The filtering performances under different clutter intensities are shown in Figure 6 and Figure 7. We observed that a higher clutter intensity leads to a larger OSPA distance and more computational cost. By comparison, the IGGM-MS-MeMBer filter is the most efficient method.

4.3. Simulation 2

This section analyzes the performance of IGGM-MS-CPHD/MeMBer filters under different sensor distributions. Suppose that there are four sensors in the two-dimensional space for observation. All sensors have the same scope with bearing θ (0–360°) and range r (0–25 km), and the other parameters are the same as those in Simulation 1. Three distributions of sensors are given as follows:
  • Case 1: assume that the four sensors are located at [−10 km, 10 km], [−10 km, −10 km], [10 km, −10 km], and [10 km, 10 km], respectively;
  • Case 2: assume that the four sensors are located at [−10 km, 0 km], [0 km, −10 km], [10 km, 0 km], and [0 km, 10 km], respectively;
  • Case 3: assume that the four sensors are located at [−9 km, −10 km], [−3 km, −10 km], [3 km, −10 km], and [9 km, −10 km], respectively.
The performance of IGGM-MS-CPHD/MeMBer filters under different sensor distributions is shown in Figure 8. It can be seen that the distribution of sensors has an important effect on the filtering errors, while the running time is barely affected by the sensor distributions. We observed that the proposed IGGM-MS-CPHD/MeMBer filters and the GM-MS-CPHD/MeMBer filters can achieve a better filtering accuracy when the sensors are dispersed in the monitoring area (e.g., Cases 1 and 2).

5. Conclusions

In this paper, a multisensor CPHD filter and a multisensor MeMBer filter were derived for the unknown and varying detection probability case. By using the IGGM model, the proposed filters can estimate the unknown and changing detection probability. The simulation results show that the proposed methods have improved filtering performances as well as higher computational costs with unknown and dynamically changing detection probability. We also conclude that the distribution of sensors has an important influence on the filtering accuracy, and the filters perform better when the sensors are dispersed in the monitoring area. Future work will explore the multisensor multitarget filter, which can estimate the clutter intensity and the detection probability adaptively, and refine the proposed algorithms in this paper to output target tracks.

Author Contributions

Z.Z. and J.S. conceived the idea, Z.Z. and Q.L. conducted the experiments and wrote the initial draft of the paper, and J.S. and Q.L. revised the manuscript. All authors read and approved the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (61471019).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mahler, R.P.S. Advances in Statistical Multisource-Multitarget Information Fusion; Artech House: Norwood, MA, USA, 2014. [Google Scholar]
  2. Quan, H.W.; Li, J.H.; Peng, D.L. Multi-target Tracking Based on Random Set Theory. Fire Control Command Control 2015, 40, 49–56. [Google Scholar]
  3. Shen, T.H.; Xue, A.K.; Zhou, Z.L. Multi-sensor Gaussian Mixture PHD Fusion for Multi-target Tracking. Acta Autom. Sin. 2017, 43, 1028–1037. [Google Scholar]
  4. Reid, D. An algorithm for tracking multiple targets. IEEE Trans. Autom. Control 1979, 24, 843–854. [Google Scholar] [CrossRef]
  5. Kirubarajan, T.; Bar-Shalom, Y. Probabilistic data association techniques for target tracking in clutter. Proc. IEEE 2004, 92, 536–557. [Google Scholar] [CrossRef]
  6. Mahler, R.P.S. Statistical Multisource-Multitarget Information Fusion; Artech House: Norwood, MA, USA, 2007. [Google Scholar]
  7. Vo, B.N.; Ma, W.K. The Gaussian Mixture Probability Hypothesis Density Filter. IEEE Trans. Signal Process. 2006, 54, 4091–4104. [Google Scholar] [CrossRef]
  8. Vo, B.T.; Vo, B.N.; Cantoni, A. Analytic Implementations of the Cardinalized Probability Hypothesis Density Filter. IEEE Trans. Signal Process. 2007, 55, 3553–3567. [Google Scholar] [CrossRef]
  9. Vo, B.T.; Vo, B.N.; Cantoni, A. On multi-Bernoulli approximations to the Bayes multi-target filter. In Proceedings of the International Colloquium on Information Fusion, Xi’an, China, 22–25 August 2007. [Google Scholar]
  10. 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]
  11. Vo, B.T.; Vo, B.N. Labeled random finite sets and multi-object conjugate priors. IEEE Trans. Signal Process. 2013, 61, 3460–3475. [Google Scholar] [CrossRef]
  12. 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]
  13. Fantacci, C.; Vo, B.T.; Papi, F.; Vo, B.N. The Marginalized delta-GLMB Filter. arXiv 2015, arXiv:1501.00926. [Google Scholar]
  14. Mertens, M.; Ulmke, M. Ground moving target tracking using signal strength measurements with the GM-CPHD filter. In Proceedings of the 2012 Workshop on Sensor Data Fusion: Trends, Solutions, Applications (SDF), Bonn, Germany, 4–6 September 2012. [Google Scholar]
  15. Maggio, E.; Taj, M.; Cavallaro, A. Efficient Multitarget Visual Tracking Using Random Finite Sets. IEEE Trans. Circuits Syst. Video Technol. 2008, 18, 1016–1027. [Google Scholar] [CrossRef] [Green Version]
  16. Kim, D.Y.; Jeon, M. Robust multi-Bernoulli filtering for visual tracking. In Proceedings of the International Conference on Control, Automation and Information Sciences (ICCAIS 2014), Gwangju, Korea, 2–5 December 2014; pp. 47–51. [Google Scholar]
  17. Mahler, R.P.S. The multisensor PHD filter: I. General solution via multitarget calculus. Proc. Spie Int. Soc. Opt. Eng. 2009, 7336. [Google Scholar]
  18. Delande, E.; Duflos, E.; Heurguier, D.; Vanheeghe, P. Multi-Target PHD Filtering: Proposition of Extensions to the Multi-Sensor Case; Research Report RR-7337; HAL-INRIA: Le Chesnay, French, 2010; p. 64. [Google Scholar]
  19. Mahler, R.P.S. The multisensor PHD filter: II. Erroneous solution via Poisson magic. In Proceedings of the SPIE International Conference on Signal Processing, Sensor Fusion, Target Recognition, Orlando, FL, USA, 13–17 April 2009. [Google Scholar]
  20. Mahler, R.P.S. Approximate multisensor CPHD and PHD filters. In Proceedings of the International Conference on Information Fusion, Edinburgh, UK, 26–29 July 2010. [Google Scholar]
  21. Nannuru, S.; Blouin, S.; Coates, M.; Rabbat, M. Multisensor CPHD filter. IEEE Trans. Aerosp. Electron. Syst. 2016, 52, 1834–1854. [Google Scholar] [CrossRef]
  22. Saucan, A.A.; Coates, M.J.; Rabbat, M. A Multisensor Multi-Bernoulli Filter. IEEE Trans. Signal Process. 2017, 65, 5495–5509. [Google Scholar] [CrossRef]
  23. Wei, B.; Nener, B.; Liu, W.; Ma, L. Centralized multi-sensor multi-target tracking with labeled random finite sets. In Proceedings of the International Conference on Control, Automation and Information Sciences (ICCAIS), Ansan, Korea, 27–29 October 2016. [Google Scholar]
  24. Ren, X.Y. A Novel Multiple Target Tracking Algorithm and Its Evaluation; Graduate University of Chinese Academy of Sciences: Beijing, China, May 2012. [Google Scholar]
  25. Mahler, R.P.S.; Vo, B.T.; Vo, B.N. CPHD filtering with unknown clutter rate and detection profile. IEEE Trans. Signal Process. 2011, 59, 3497–3513. [Google Scholar] [CrossRef]
  26. Vo, B.T.; Vo, B.N.; Hoseinnezhad, R.; Mahler, R.P.S. Multi-Bernoulli filtering with unknown clutter intensity and sensor field-of-view. In Proceedings of the 45th Annual Conference on Information Sciences and Systems (CISS), Baltimore, MD, USA, 23–25 March 2011. [Google Scholar]
  27. Vo, B.T.; Vo, B.N.; Hoseinnezhad, R.; Mahler, R.P.S. Robust Multi-Bernoulli filtering. IEEE J. Sel. Top. Signal Process. 2013, 7, 399–409. [Google Scholar] [CrossRef]
  28. Punchihewa, Y.; Vo, B.T.; Vo, B.N.; Kim, D.Y. Multiple Object Tracking in Unknown Backgrounds with Labeled Random Finite Sets. IEEE Trans. Signal Process. 2017, 66, 3040–3055. [Google Scholar] [CrossRef]
  29. Li, C.; Wang, W.; Kirubarajan, T.; Sun, J.; Lei, P. PHD and CPHD Filtering with Unknown Detection Probability. IEEE Trans. Signal Process. 2018, 66, 3784–3798. [Google Scholar] [CrossRef]
  30. Qian, K.; Zhou, H.; Qin, H.; Rong, S.; Zhao, D.; Du, J. Guided filter and convolutional network based tracking for infrared dim moving target. Infrared Phys. Technol. 2017, 85, 431–442. [Google Scholar] [CrossRef]
  31. Lerro, D.; Bar-Shalom, Y. Automated Tracking with Target Amplitude Information. In Proceedings of the American Control Conference, San Diego, CA, USA, 23–25 May 1990; pp. 2875–2880. [Google Scholar]
  32. Clark, D.; Ristic, B.; Vo, B.N.; Vo, B.T. Bayesian Multi-Object Filtering With Amplitude Feature Likelihood for Unknown Object SNR. IEEE Trans. Signal Process. 2010, 58, 26–37. [Google Scholar] [CrossRef]
  33. Skolnik, M. Radar Handbook, 3rd ed.; McGraw-Hill: New York, NY, USA, 2008. [Google Scholar]
  34. 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]
  35. Grewal, M.S.; Andrews, A.P. Kalman Filtering: Theory and Practice Using MATLAB, 3rd ed.; Wiley: Hoboken, NJ, USA, 2008. [Google Scholar]
Figure 1. Truth trajectories of the target.
Figure 1. Truth trajectories of the target.
Electronics 08 00741 g001
Figure 2. Target detection probability and the feature a.
Figure 2. Target detection probability and the feature a.
Electronics 08 00741 g002
Figure 3. Measurement of the sensors.
Figure 3. Measurement of the sensors.
Electronics 08 00741 g003
Figure 4. Optimal subpattern assignment (OSPA) distance.
Figure 4. Optimal subpattern assignment (OSPA) distance.
Electronics 08 00741 g004
Figure 5. Estimation of the cardinality.
Figure 5. Estimation of the cardinality.
Electronics 08 00741 g005
Figure 6. OSPA distances for different clutter intensities.
Figure 6. OSPA distances for different clutter intensities.
Electronics 08 00741 g006
Figure 7. Average running times for different clutter intensities.
Figure 7. Average running times for different clutter intensities.
Electronics 08 00741 g007
Figure 8. Performance under different sensor distributions.
Figure 8. Performance under different sensor distributions.
Electronics 08 00741 g008
Table 1. The target initial states.
Table 1. The target initial states.
Target State   x ˜ uvSurvival Time (s)
1[6000, −80, 6000, −80]T51400[1 70]
2[6000, −100, 6000, −15]T51400[10 70]
3[−6000, −15, 6000, −100]T41360[15 80]
4[−6000, 110, 6000, −75]T41360[20 100]
5[0, −50, −6000, 75]T51500[30 100]
6[0, 70, −6000, 60]T51500[40 100]
Table 2. The scopes of the sensors.
Table 2. The scopes of the sensors.
SensorRange rBearing θ
10–10 km0–360°
20–25 km180–270°
30–25 km270–360°
40–25 km0–90°
50–25 km90–180°
Table 3. Average single-step running time.
Table 3. Average single-step running time.
IGGM-MS-CPHDIGGM-MS-MeMBerGM-MS-CPHDGM-MS-MeMBer
180 ms75 ms127 ms58 ms

Share and Cite

MDPI and ACS Style

Zhang, Z.; Li, Q.; Sun, J. Multisensor RFS Filters for Unknown and Changing Detection Probability. Electronics 2019, 8, 741. https://doi.org/10.3390/electronics8070741

AMA Style

Zhang Z, Li Q, Sun J. Multisensor RFS Filters for Unknown and Changing Detection Probability. Electronics. 2019; 8(7):741. https://doi.org/10.3390/electronics8070741

Chicago/Turabian Style

Zhang, Zhiguo, Qing Li, and Jinping Sun. 2019. "Multisensor RFS Filters for Unknown and Changing Detection Probability" Electronics 8, no. 7: 741. https://doi.org/10.3390/electronics8070741

APA Style

Zhang, Z., Li, Q., & Sun, J. (2019). Multisensor RFS Filters for Unknown and Changing Detection Probability. Electronics, 8(7), 741. https://doi.org/10.3390/electronics8070741

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