Next Article in Journal / Special Issue
The Positioning Accuracy of BAUV Using Fusion of Data from USBL System and Movement Parameters Measurements
Previous Article in Journal
Development of Gated Pinned Avalanche Photodiode Pixels for High-Speed Low-Light Imaging
Previous Article in Special Issue
Regular Deployment of Wireless Sensors to Achieve Connectivity and Information Coverage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sensor Fusion of Gaussian Mixtures for Ballistic Target Tracking in the Re-Entry Phase

School of Automation Science and Electrical Engineering, Beihang University, Beijing 100191, China
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(8), 1289; https://doi.org/10.3390/s16081289
Submission received: 12 July 2016 / Revised: 5 August 2016 / Accepted: 9 August 2016 / Published: 15 August 2016
(This article belongs to the Special Issue Advances in Multi-Sensor Information Fusion: Theory and Applications)

Abstract

:
A sensor fusion methodology for the Gaussian mixtures model is proposed for ballistic target tracking with unknown ballistic coefficients. To improve the estimation accuracy, a track-to-track fusion architecture is proposed to fuse tracks provided by the local interacting multiple model filters. During the fusion process, the duplicate information is removed by considering the first order redundant information between the local tracks. With extensive simulations, we show that the proposed algorithm improves the tracking accuracy in ballistic target tracking in the re-entry phase applications.

1. Introduction

The problem of ballistic target tracking in the re-entry phase has attracted much attention due to both of its theoretical and practical significance. Technically speaking we need to develop a stochastic nonlinear filter for state estimation with respect to the ballistic target dynamics. Practical applications are in the fields of surveillance for safety against the re-entry of space debris produced by old satellites or spacecraft at the end of their lifetime [1,2,3]. Tracking a ballistic target has been identified as a stressful filtering problem due to the strong nonlinearity exhibited by the forces acting on the target [2,4,5]. Besides, the nonlinearity also results from the bearing and range measurements which are given with nonlinear measurement functions in relation to the target [6]. In the literature several nonlinear filtering methods have been applied to this tracking problem under the assumption that the ballistic coefficient is already known [7], including extended Kalman filter [8], unscented Kalman filter [9], ensemble Kalman filter [10], cubature Kalman filter [11] and particle filter [3]. The tracking problem becomes more difficult when the ballistic coefficient is known crudely or totally unknown. In such a case, an approach which has been investigated to be effective is the multiple model approach [12,13,14,15]. The structure of multiple model filter consists of a bank of local filters, with each filter modelled with a different value of ballistic coefficient. The output of the multiple model filter is in terms of Gaussian mixtures obtained by weighted mixing of the estimate from each local filter.
The technology of track-to-track fusion (T2TF) makes it possible to achieve further improvement of tracking accuracy than in the single sensor case. A basic T2TF system [16] consists of two sensors and a fusion center. Each sensor generates a local track by observing the target of interest and estimating with filtering. The fusion center fuses local tracks to obtain a fused track with improved accuracy. There has been a great deal of work in developing T2TF algorithms [17,18,19,20], however, most of these algorithms have been designed for systems with Gaussian density based filters (e.g., Kalman filter, extended Kalman filter and unscented Kalman filter) and may not be applicable to systems with multiple model filters. Demand for tracking ballistic targets with unknown ballistic coefficients poses a new requirement on using multiple model filters, which makes it necessary to fuse local tracks in terms of the Gaussian mixtures model [21] for tracking accuracy improvement. In the case of Gaussian mixtures, an empirical study into the use of Chernoff information has been investigated in [22], and an approximate approach for Chernoff fusion of Gaussian mixtures has been proposed in [23]. One drawback of this approach is that an exhaustive search must be performed to find the optimal weighting parameter [24], which makes it computationally intensive and difficult to implement in practical systems.
In this study, we consider a distributed ballistic target tracking scenario where the estimates of the local sensors are to be fused. We assume that the local sensors run IMM filters for handling unknown ballistic coefficient and the output of the local estimators are Gaussian mixtures. As a result, the track-to-track fusion problem we consider involves the fusion of Gaussian mixtures densities. The main contribution of this paper is to obtain the common information between the local estimates by considering the first order redundant information and derive the required fusion expressions to be employed in the fusion center. The proposed fusion rule is computationally efficient since it yields an analytical fused estimate for Gaussian mixtures. Simulation results show that the track-to-track fusion with Gaussian mixtures provides estimation accuracy improvement over the single Gaussian case. The remainder of the paper is organized as follows: the problem is formulated in Section 2. The Gaussian mixtures fusion rule is presented in Section 3. The simulation results are given in Section 4. The conclusions are presented in Section 5.

2. Problem Formulation

2.1. Target Dynamics

The kinematics of the ballistic target in the re-entry phase is derived under the following assumptions [7]: the forces acting on the target include Earth’s gravity and the aerodynamic drag; the effects of centrifugal acceleration, Coriolis acceleration, wind, lift force and spinning motion are ignored due to their minor effect on the target trajectory. Thus the system dynamics of a ballistic target [3,6] is given by:
x k = Ψ ( x k 1 ) + G [ 0 g ] + w k
where x k 1 = [ x p , k 1 , x v , k 1 , y p , k 1 , y v , k 1 ] T is the target state vector at time step k 1 . x p , k 1 and y p , k 1 represent positions, x v , k 1 and y v , k 1 represent velocities in Cartesian coordinates ( x , y ) . w k is a zero-mean white Gaussian process noise vector with covariance Q given by:
Q = q [ Θ 0 0 Θ ]  with  Θ = [ T 3 / 3 T 2 / 2 T 2 / 2 T ]
where q is a parameter related to the process noise intensity. The nonlinear function Ψ ( · ) in Equation (1) is given by:
Ψ ( x k 1 ) = F x k 1 + G f ( x k 1 )
where matrices F and G are given by:
F = [ 1 T 0 0 0 1 0 0 0 0 1 T 0 0 0 1 ] ,    G = [ T 2 / 2 0 T 0 0 T 2 / 2 0 T ]
and f ( ) is the aerodynamic drag given by
f ( x k 1 ) = g ρ ( y p , k 1 ) x v , k 1 2 + y v , k 1 2 [ x v , k 1 y v , k 1 ] 2 β
where β is the ballistic coefficient, g is the gravity acceleration, ρ ( ) is the air density [25] which decays exponentially with the altitude y p , k 1 as
ρ ( y p , k 1 ) = c 1 e c 2 y p , k 1

2.2. Sensor Configurations and Measurement Model

The sensor configuration involves two ground-radar sensors S 1 and S 2 , which are located at ( x R S l , y R S l ) for l = 1 , 2 respectively. The measurement equation is given by:
z k S l = h ( x k ) + v k S l  for  l = 1 , 2
where z k S l = [ r k S l , θ k S l ] T is the measurement vector, v k = [ v r , k S l , v θ , k S l ] T is the measurement noise for sensor S l . The sensor measures the range r k S l and the bearing θ k in relation to the target. Thus the measurement model can be written as
r k S l = ( x p , k S l x R S l ) 2 + ( y p , k S l y R S l ) 2 + v r , k S l θ k S l = tan 1 ( y p , k S l y R S l x p , k S l x R S l ) + v θ , k S l  for  l = 1 , 2
where
v r , k S l N ( 0 , R r S l ) v θ , k S l N ( 0 , R θ S l )

2.3. Local Filter

Since the actual ballistic coefficient is unknown to the sensors, an interactive multiple model (IMM) filter [26] with r = 2 models is applied at the single sensor level. The models are based on the extended Kalman filter (EKF), but having different value of ballistic coefficient β i S l , i = 1 , , r . In general, the IMM filter merges state estimates computed under each possible model using local filters, with each filter using a different combination of the previous model-conditioned estimates [27]. The recursion equations for sensor S l consist of four major steps in each cycle as follows.
Step 1: Calculation of the mixing probabilities for i , j = 1 , , r
μ i | j , k 1 | k 1 S l = 1 c ¯ j S l p i j S l μ i , k 1 S l
where c ¯ j S l = i = 1 r p i j S l μ i , k 1 S l .
Step 2: Mixing for j = 1 , , r :
x ^ 0 j , k 1 | k 1 S l = i = 1 r x ^ i , k 1 | k 1 S l u i | j , k 1 | k 1 S l P 0 j , k 1 | k 1 S l = i = 1 r μ i | j , k 1 | k 1 S l [ P i , k 1 | k 1 S l + ( x ^ i , k 1 | k 1 S l x ^ 0 j , k 1 | k 1 S l ) ( x ^ i , k 1 | k 1 S l x ^ 0 j , k 1 | k 1 S l ) T ]
Step 3: Mode-matched filtering for j = 1 , , r
x ^ j , k | k 1 S l = Ψ ( x ^ 0 j , k 1 | k 1 S l ) + G [ 0 g ] P j , k | k 1 S l = ( F + G f J , j , k S l ) P 0 j , k 1 | k 1 S l ( F + G f J , j , k S l ) T + Q
x ^ j , k | k S l = x ^ j , k | k 1 S l + K k S l ( z k S l h ( x ^ j , k | k 1 S l ) ) P j , k | k S l = ( I K k S l h J , j , k S l ) P j , k | k 1 S l
where K k S l = P j , k | k 1 S l h J , j , k T ( h J , j , k P j , k | k 1 S l h J , j , k T + R S l ) 1 and R S l = [ R r S l 0 0 R θ , k S l ] . f J , j , k S l is the Jacobian of f S l ( ) calculated at the estimated state x ^ 0 j , k 1 | k 1 S l , h J , j , k S l is the Jacobian of h S l ( ) calculated at the estimated state x ^ j , k | k 1 S l . The detailed derivation of f J , j , k S l and h J , j , k S l for models given in Equations (5) and (8) are shown in Appendix A and Appendix B respectively.
Step 4: Mode probability update for j = 1 , , r
Λ j , k S l = N ( z k S l ; h ( x ^ j , k | k 1 S l ) , h J , j , k P j , k | k 1 S l h J , j , k T + R S l ) μ j , k S l = 1 c S l Λ j , k S l c ¯ j S l c S l = j = 1 r Λ j , k S l c ¯ j S l
A block diagram for one cycle of the IMM filter with two models ( r = 2 ) for sensor S l is given in Figure 1.

2.4. Fusion Architecture

Among the various choices for multiple model filters, the IMM filter has been shown to be one of the most effective schemes for hybrid systems [28]. In this study we assume that the local sensors run IMM filters for handling unknown ballistic coefficient, therefore the output of each local sensor is a Gaussian mixtures model. In this regard, the track-to-track fusion (T2TF) problem we consider here involves the fusion of Gaussian mixtures densities. As shown in Figure 2, the fusion architecture presented includes two ground-radar sensors S1 and S2 with local measurements obtained periodically and local tracks updated synchronously. In addition, the fusion center may transmit the latest fused track back to one local sensor. When the local sensor receives the fused track, it will be used to replace the local one. In practice, this condition usually happens when the fusion center is collocated with one local sensor [29]. In this study, we assume that after each fusion, the local track of sensor S1 is replaced with the fused track from the fusion center, while sensor S2 operates independently without information feedback from the fusion center.

3. Gaussian Mixtures Fusion

3.1. Basic Fusion Process and Redundant Information

The basic fusion process [19] is formulated as:
p f ( x ) = 1 c p S 1 ( x ) p S 2 ( x ) p c ( x )
where p f ( x ) is the fused estimate, p S 1 ( x ) and p S 2 ( x ) are the estimates from local sensors S 1 and S 2 respectively, c is a normalization constant. The common information between the local estimates is given in the denominator by p c ( x ) , which is to be subtracted out. While the removal of duplicate information is straightforward in the theoretical formulation [30], identification of duplicate information for distributed estimation system can be difficult in practical implementation. In this paper, the first order redundant information [16,31] is considered. In this way, the fusion center only needs to keep track of the previous data received from sensor S2 at the previous communication time step and remove it when fusing the current estimates from sensor S1 and S2.
Specifically, when both p S 1 ( x ) and p S 2 ( x ) are Gaussian mixtures, namely:
p S l ( x ) = i = 1 r μ i , k S l N ( x ; x ^ i , k | k S l , P i , k | k S l )  for  l = 1 , 2
where r = 2 since the IMM filter is with two models for each sensor. Then the common information at time step k can be obtained as:
p c ( x ) = i = 1 2 μ i , k c N ( x ; x ^ i , k | k c , P i , k | k c )
where:
μ i , k c = j = 1 r μ j , k 1 S 2 μ j | i , k 1 | k 1 S 2 x ^ i , k | k c = Ψ ( x ^ 0 i , k 1 | k 1 S 2 ) + G [ 0 g ] P i , k | k c = ( F + G f J , i , k S 2 ) P 0 i , k 1 | k 1 S 2 ( F + G f J , i , k S 2 ) T + Q  for  i = 1 , 2
and f J , i , k S 2 ( ) is the Jacobian calculated at the estimated state x ^ 0 i , k 1 | k 1 S 2 .

3.2. Fusion of Gaussian Mixtures

In a distributed fusion problem, assuming two Gaussian mixtures, p S 1 ( x ) = i = 1 2 μ i , k S 1 N ( x ; x ^ i , k | k S 1 , P i , k | k S 1 ) and p S 2 ( x ) = i = 1 2 μ i , k S 2 N ( x ; x ^ i , k | k S 2 , P i , k | k S 2 ) are to be fused with a common prior distribution p c ( x ) . With the standard Bayesian fusion formula, i.e., Equation (15), the fused probability density function (PDF) can be obtained as:
p f ( x ) = 1 c i = 1 2 j = 1 2 μ i , k S 1 μ j , k S 2 N ( x ; x ^ i , k | k S 1 , P i , k | k S 1 ) N ( x ; x ^ j , k | k S 2 , P j , k | k S 2 ) p c ( x )
To obtain an analytical fused result and avoid the potential complexity, one idea is to approximate the denominator p c ( x ) with a single Gaussian PDF using moment matching [27], namely:
p c ( x ) = i = 1 2 μ i , k c N ( x ; x ^ i , k | k c , P i , k | k c ) N ( x ; x ^ k | k c , P k | k c )
where:
x ^ k | k c = i = 1 2 μ i , k c x ^ i , k | k c P k | k c = i = 1 2 μ i , k c [ P i , k | k c + ( x ^ i , k | k c x ^ k | k c ) ( x ^ i , k | k c x ^ k | k c ) T ]
With this approximation, the fusion expressions to be applied in IMM filter which is the major contribution of the current work are derived as follows:
p f ( x ) 1 c i = 1 2 j = 1 2 μ i , k S 1 μ j , k S 2 N ( x ; x ^ i , k | k S 1 , P i , k | k S 1 ) N ( x ; x ^ j , k | k S 2 , P j , k | k S 2 ) N ( x ; x ^ k | k c , P k | k c ) = 1 c i = 1 2 j = 1 2 μ i , k S 1 μ j , k S 2 c i j , k | k N ( x ; x ^ i j , k | k , P i j , k | k )
where
P i j , k | k 1 = ( P i , k | k S 1 ) 1 + ( P j , k | k S 2 ) 1 ( P k | k c ) 1 P i j , k | k 1 x ^ i j , k | k = ( P i , k | k S 1 ) 1 x ^ i , k | k S 1 + ( P j , k | k S 2 ) 1 x ^ j , k | k S 2 ( P k | k c ) 1 x ^ k | k c
and
c i j , k | k = exp ( ζ i , k | k S 1 + ζ j , k | k S 2 ζ k | k c ζ i j , k | k )
with
ζ i , k | k S 1 = 1 2 ( d log 2 π log | ( P i , k | k S 1 ) 1 | + ( x ^ i , k | k S 1 ) T ( P i , k | k S 1 ) T x ^ i , k | k S 1 ) ζ j , k | k S 2 = 1 2 ( d log 2 π log | ( P j , k | k S 2 ) 1 | + ( x ^ j , k | k S 2 ) T ( P j , k | k S 2 ) T x ^ j , k | k S 2 ) ζ k | k c = 1 2 ( d log 2 π log | ( P k | k c ) 1 | + ( x ^ k | k c ) T ( P k | k c ) T x ^ k | k c ) ζ i j , k | k = 1 2 ( d log 2 π log | P i j , k | k 1 | + x ^ i j , k | k T P i j , k | k T x ^ i j , k | k )
The derivation of Equations (23)–(25) is presented in Appendix C.

3.3. Gaussian Mixtures Reduction

As one can see from Equation (22), the fused Gaussian mixtures PDF has an exponentially growing number of components as more Gaussian mixtures are multiplied in the long run. Thus it is necessary to manage the components growth with a method of Gaussian mixtures reduction. Suppose that we are given a Gaussian mixtures model with n components, and we wish to approximate it with a mixture of r components ( r < n ). In general, the Gaussian mixtures reduction algorithm can be operated in the following manner.
While more than r components remain, choose two components that in a sense to be least dissimilar and replace them by their moment matching merge as:
w m = w i + w j
x ^ m = w i w m x ^ i + w j w m x ^ j
P m = i = 1 2 w i w m [ P i + ( x ^ i x ^ m ) ( x ^ i x ^ m ) T ]
where ( w i , x ^ i , P i ) and ( w j , x ^ j , P j ) are two weighted Gaussian components to be merged, and ( w m , x ^ m , P m ) is their moment matching approximation.
For the dissimilarity measure between two components of a Gaussian mixtures model, we adopt a metric that is proposed in Reference [32] as an upper bound on the discrimination of the Gaussian mixtures after the merge from the Gaussian mixtures before the merge. The dissimilarity measure is given as follows:
D ( ( w i , x ^ i , P i ) , ( w j , x ^ j , P j ) ) ( w i + w j ) ln ( det ( P m ) ) w i ln ( det ( P i ) ) w j ln ( det ( P j ) ) 2
Thus in each iteration, a Gaussian mixtures model with n components is constructed after the Gaussian mixtures fusion, then we operate the procedure of Gaussian mixtures reduction in an iterative manner until n = r , so as to meet the requirement that the IMM filter is set to be with r models. In our case, the fused Gaussian mixtures consist of n = 4 components after fusion in Equation (22), namely:
p f ( x ) = i = 1 2 j = 1 2 w i j N ( x ; x ^ i j , P i j )
where:
w i j = μ i S 1 μ j S 2 c i j c
Note that the subscript k | k is omitted for the sake of brevity. We set ( w 11 , x ^ 11 , P 11 ) and ( w 22 , x ^ 22 , P 22 ) as major components, since they are consistent with model 1 and 2 in a straightforward manner. For the rest two cross-components, namely ( w 12 , x ^ 12 , P 12 ) and ( w 21 , x ^ 21 , P 21 ) , they are supposed to be merged with one of the major components according to the dissimilarity measure given in Equation (29). More specifically, for ( w 12 , x ^ 12 , P 12 ) , if D ( ( w 12 , x ^ 12 , P 12 ) , ( w 11 , x ^ 11 , P 11 ) ) < D ( ( w 12 , x ^ 12 , P 12 ) , ( w 22 , x ^ 22 , P 22 ) ) , then we merge ( w 12 , x ^ 12 , P 12 ) with ( w 11 , x ^ 11 , P 11 ) , otherwise we merge ( w 12 , x ^ 12 , P 12 ) with ( w 22 , x ^ 22 , P 22 ) ; for ( w 21 , x ^ 21 , P 21 ) , if D ( ( w 21 , x ^ 21 , P 21 ) , ( w 11 , x ^ 11 , P 11 ) ) < D ( ( w 21 , x ^ 21 , P 21 ) , ( w 22 , x ^ 22 , P 22 ) ) , then we merge ( w 21 , x ^ 21 , P 21 ) with ( w 11 , x ^ 11 , P 11 ) , otherwise we merge ( w 21 , x ^ 21 , P 21 ) with ( w 22 , x ^ 22 , P 22 ) . In the end, a flow chart of the proposed fusion methodology is given in Figure 3 as follows:

4. Simulation Results

The proposed algorithm is applied in a re-entry ballistic target tracking scenario to verify the performance through a series of simulation runs. To the best author’s knowledge, the Gaussian mixtures fusion methodology has not been applied in a distributed ballistic target tracking scenario for track-to-track fusion before. The system dynamic model is represented as Equation (1), where we use T = 2 s, q = 1 m2·s−3 and g = 9.81 m·s2. The initial state x 0 is:
x 0 = [ 23200 m 2290 cos ( 190 ° ) m s 1 88000 m 2290 sin ( 190 ° ) m s 1 ]
with initial covariance:
P 0 = [ 1000 2 m 2 0 0 0 0 20 2 m 2 s 2 0 0 0 0 1000 2 m 2 0 0 0 0 20 2 m 2 s 2 ]
For the air density ρ ( ) , we have c1 = 1.227, c2 = 1.093 × 10−4 for yp < 9144 m and, and c1 = 1.754, c2 = 1.49 × 10−4 for yp > 9144 m. The actual target ballistic coefficient is β = 40,000 kg·m−1·s−2.
Figure 4 shows the target trajectory in the X-Y plane. Figure 5 and Figure 6 show the velocity of the ballistic target, and the aerodynamic drag acceleration against time, respectively. It can be observed that the velocity decreases with the increment of aerodynamic drag f ( ) .
For the measurement model, we assume that the two sensors are homogeneous, in the sense that R r S 1 = R r S 2 = 100 2   m 2 and R θ S 1 = R θ S 2 = 0.05 2   rad 2 . Besides, we consider that the sensors are located at ( x R S 1 , y R S 1 ) = ( 0 , 0 ) and ( x R S 2 , y R S 2 ) = ( 50,000 , 0 ) , respectively.
For the local filters, it is assumed that the actual ballistic coefficient of the target is unknown to both sensors, so that β 1 S 1 = β 1 S 2 = 60,000  kg · m 1 · s 2 and β 2 S 1 = β 2 S 2 = 10,000  kg · m 1 · s 2 are used as the IMM models for both sensor S 1 and S 2 . Besides, the Markov chain transition matrix was taken as:
[ p i j ] = [ 0.95 0.05 0.05 0.95 ]
The proposed Gaussian mixtures fusion method is verified by comparing the results between sensors S1 and S2. Note that the estimate from sensor S1 is identical to the fused result, for the reason that after each fusion the local track from sensor S1 is replaced with the fused track from the fusion center. Unlike sensor S1, the estimate from sensor S2 is solely obtained from one singe IMM filter, since sensor S2 operates by itself without information feedback from the fusion center.
Figure 7 and Figure 8 show the estimated model probabilities for sensors S1 and S2, respectively. It can be seen that during 0–50 s the probabilities for models 1 and 2 are equal to each other, for the reason that the aerodynamic drag is zero which makes it not able to judge which model matches better. The ballistic coefficient becomes observable from t = 50 s, at that time the aerodynamic drag begins to emerge. It can be observed that the model probabilities for sensor S1 is more stable than that for sensor S2, namely the proposed fusion method improves the stability of the estimated model probabilities.
We further calculate the estimated ballistic coefficient with the total probability theorem [25] as:
β ^ S l = i = 1 2 β i S l μ i S l  for  l = 1 , 2
where μ i S l is the estimated model probability for β i S l . Figure 9 and Figure 10 plot the estimated value of ballistic coefficient from sensors S1 and S2, respectively. It can be seen that after the drag force emerges, the estimated value of the ballistic coefficient for sensor S1F, namely the fused one, appears to be more accurate and stable than that for sensor S2, in the sense that the first one ranges around [3.6,4.4] × 104, and the latter one ranges around [1.5,6.0] × 104.
In the following, we compare the position estimation accuracy in terms of root mean square error (RMSE) between the proposed Gaussian mixtures fusion method, the single IMM filter and the covariance intersection (CI) method. The CI method [18] is the most well-known fusion technique which yields consistent estimates by optimizing a nonlinear cost function associated with the fused covariance. For the reason that CI cannot be applied to Gaussian mixtures fusion directly, it requires one to approximate the output of IMM filter with the single Gaussian density using moment matching before track fusion. Figure 11 and Figure 12 and Table 1 show the RMSE of the algorithms in X and Y axes, respectively, by performing 1000 runs of Monte Carlo simulation. It can be seen that both of the fusion methods, namely the proposed Gaussian mixtures fusion, and CI provide performance improvement over the single IMM filter. Furthermore, the performance of the proposed Gaussian mixtures fusion is better than CI. The first reason is that we remove the first order redundant information between the local tracks during Gaussian mixtures fusion process, however, the CI only provides a conservative estimate due to its ignorance of the cross correlations between the local tracks.
The second reason is that instead of fusing Gaussian mixtures directly as for the proposed fusion method, for CI one needs to approximate the original local Gaussian mixtures estimate with mean and covariance before track fusion, thus the resulting approximation error could degrade the estimation performance.

5. Conclusions

A Gaussian mixtures fusion algorithm for the track-to-track fusion problem is proposed in this study. The common information to be reduced during the fusion process is approximated with the first order redundant information between the local tracks. The proposed fusion algorithm is applied to tracking a ballistic target with unknown ballistic coefficient using IMM filters. A series of Monte Carlo simulations are conducted to evaluate the sensor fusion performance. The results indicate that the proposed algorithm improves the estimation accuracy in terms of the root mean square error. There are two issues to be addressed in the future. First, a more detailed analysis of the effects of the first order approximation must be carried out. Second, experiments should be conducted to explore the effect of using higher order approximation to formulate the common information between the local estimates.

Acknowledgments

This study was partially supported by the National Natural Science Foundation of China (Grant nos. 61273349 and 61573042).

Author Contributions

K.L. wrote the paper, R.Z. conceived the experiments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

f J , j , k S l ( 1 , 1 ) = 0
f J , j , k S l ( 1 , 2 ) = g ρ ( x ^ 0 j , k 1 | k 1 S l ( 3 ) ) 2 β S l 2 ( x ^ 0 j , k 1 | k 1 S l ( 2 ) ) 2 + ( x ^ 0 j , k 1 | k 1 S l ( 4 ) ) 2 ( x ^ 0 j , k 1 | k 1 S l ( 2 ) ) 2 + ( x ^ 0 j , k 1 | k 1 S l ( 4 ) ) 2
f J , j , k S l ( 1 , 3 ) = g c 2 ρ ( x ^ 0 j , k 1 | k 1 S l ( 3 ) ) 2 β S l x ^ 0 j , k 1 | k 1 S l ( 2 ) ( x ^ 0 j , k 1 | k 1 S l ( 2 ) ) 2 + ( x ^ 0 j , k 1 | k 1 S l ( 4 ) ) 2
f J , j , k S l ( 1 , 4 ) = g ρ ( x ^ 0 j , k 1 | k 1 S l ( 3 ) ) 2 β S l ( x ^ 0 j , k 1 | k 1 S l ( 2 ) ) ( x ^ 0 j , k 1 | k 1 S l ( 4 ) ) ( x ^ 0 j , k 1 | k 1 S l ( 2 ) ) 2 + ( x ^ 0 j , k 1 | k 1 S l ( 4 ) ) 2
f J , j , k S l ( 2 , 1 ) = 0
f J , j , k S l ( 2 , 2 ) = g ρ ( ( x ^ 0 j , k 1 | k 1 S l ( 3 ) ) ) 2 β S l ( x ^ 0 j , k 1 | k 1 S l ( 2 ) ) ( x ^ 0 j , k 1 | k 1 S l ( 4 ) ) ( x ^ 0 j , k 1 | k 1 S l ( 2 ) ) 2 + ( x ^ 0 j , k 1 | k 1 S l ( 4 ) ) 2
f J , j , k S l ( 2 , 3 ) = g c 2 ρ ( x ^ 0 j , k 1 | k 1 S l ( 3 ) ) 2 β S l × ( x ^ 0 j , k 1 | k 1 S l ( 4 ) ) ( x ^ 0 j , k 1 | k 1 S l ( 2 ) ) 2 + ( x ^ 0 j , k 1 | k 1 S l ( 4 ) ) 2
f J , j , k S l ( 2 , 4 ) = g   ρ ( x ^ 0 j , k 1 | k 1 S l ( 3 ) ) 2 β S l · x ^ 0 j , k 1 | k 1 S l ( 2 ) + 2 x ^ 0 j , k 1 | k 1 S l ( 4 ) ( x ^ 0 j , k 1 | k 1 S l ( 2 ) ) 2 + ( x ^ 0 j , k 1 | k 1 S l ( 4 ) ) 2

Appendix B

h J , j , k S l ( 1 , 1 ) = x ^ j , k | k 1 S l ( 1 ) x R S l ( x ^ j , k | k 1 S l ( 1 ) x R S l ) 2 + ( x ^ j , k | k 1 S l ( 3 ) y R S l ) 2
h J , j , k S l ( 1 , 2 ) = 0
h J , j , k S l ( 1 , 3 ) = x ^ j , k | k 1 S l ( 3 ) - y R S l ( x ^ j , k | k 1 S l ( 1 ) - x R S l ) 2 + ( x ^ j , k | k 1 S l ( 3 ) - y R S l ) 2
h J , j , k S l ( 1 , 4 ) = 0
h J , j , k S l ( 2 , 1 ) = x ^ j , k | k 1 S l ( 3 ) y R S l ( x ^ j , k | k 1 S l ( 1 ) x R S l ) 2 + ( x ^ j , k | k 1 S l ( 3 ) y R S l ) 2
h J , j , k S l ( 2 , 2 ) = 0
h J , j , k S l ( 2 , 3 ) = x ^ j , k | k 1 S l ( 1 ) x R S l ( x ^ j , k | k 1 S l ( 1 ) x R S l ) 2 + ( x ^ j , k | k 1 S l ( 3 ) y R S l ) 2
h J , j , k S l ( 2 , 4 ) = 0

Appendix C

The multivariate Gaussian PDF N ( x ; x ^ , P ) can be written in canonical notation [33] as:
p ( x ) = N ( x ; x ^ , P ) = exp ( ζ + y T x 1 2 x T Y x )
where:
Y = P 1 y = P 1 x ^ ζ = 1 2 ( d log 2 π log | Y | + y T Y 1 y )
d is the dimensionality of x , | Y | det ( Y ) .
Define:
N ( x ; x ^ i , k | k S 1 , P i , k | k S 1 ) = exp ( ζ i , k | k S 1 + ( y i , k | k S 1 ) T x 1 2 x T Y i , k | k S 1 x ) N ( x ; x ^ j , k | k S 2 , P j , k | k S 2 ) = exp ( ζ j , k | k S 2 + ( y j , k | k S 2 ) T x 1 2 x T Y j , k | k S 2 x ) N ( x ; x ^ k | k c , P k | k c ) = exp ( ζ k | k c + ( y k | k c ) T x 1 2 x T Y k | k c x )
so we have:
N ( x ; x ^ i , k | k S 1 , P i , k | k S 1 ) N ( x ; x ^ j , k | k S 2 , P j , k | k S 2 ) N ( x ; x ^ k | k c , P k | k c ) = exp [ ( ζ i , k | k S 1 + ζ j , k | k S 2 ζ k | k c ) + ( y i , k | k S 1 + y j , k | k S 2 y k | k c ) T x 1 2 x T ( Y i , k | k S 1 + Y j , k | k S 2 Y k | k c ) x ]
We further define:
Y i j , k | k Y i , k | k S 1 + Y j , k | k S 2 Y k | k c y i j , k | k y i , k | k S 1 + y j , k | k S 2 y k | k c ζ i j , k | k 1 2 ( d log 2 π log | Y i j , k | k | + y i j , k | k T Y i j , k | k 1 y i j , k | k )
Then we have:
N ( x ; x ^ i , k | k S 1 , P i , k | k S 1 ) N ( x ; x ^ j , k | k S 2 , P j , k | k S 2 ) N ( x ; x ^ k | k c , P k | k c ) = exp [ ( ζ i , k | k S 1 + ζ j , k | k S 2 ζ k | k c ) ζ i j , k | k + ζ i j , k | k + ( y i , k | k S 1 + y j , k | k S 2 y k | k c ) T x 1 2 x T ( Y i , k | k S 1 + Y j , k | k S 2 Y k | k c ) x ] = exp [ ( ζ i , k | k S 1 + ζ j , k | k S 2 ζ k | k c ) ζ i j , k | k + ζ i j , k | k + y i j , k | k T x 1 2 x T Y i j , k | k x ] = exp ( ζ i , k | k S 1 + ζ j , k | k S 2 ζ k | k c ζ i j , k | k ) exp ( ζ i j , k | k + y i j , k | k T x 1 2 x T Y i j , k | k x )
Therefore N ( x ; x ^ i , k | k S 1 , P i , k | k S 1 ) N ( x ; x ^ j , k | k S 2 , P j , k | k S 2 ) N ( x ; x ^ k | k c , P k | k c ) is a scaled Gaussian PDF as
N ( x ; x ^ i , k | k S 1 , P i , k | k S 1 ) N ( x ; x ^ j , k | k S 2 , P j , k | k S 2 ) N ( x ; x ^ k | k c , P k | k c ) = c i j , k | k N ( x ; x ^ i j , k | k , P i j , k | k )
where:
P i j , k | k 1 = ( P i , k | k S 1 ) 1 + ( P j , k | k S 2 ) 1 ( P k | k c ) 1 P i j , k | k 1 x ^ i j , k | k = ( P i , k | k S 1 ) 1 x ^ i , k | k S 1 + ( P j , k | k S 2 ) 1 x ^ j , k | k S 2 ( P k | k c ) 1 x ^ k | k c
and the scaling factor is:
c i j , k | k = exp ( ζ i , k | k S 1 + ζ j , k | k S 2 ζ k | k c ζ i j , k | k )
with:
ζ i , k | k S 1 = 1 2 ( d log 2 π log | ( P i , k | k S 1 ) 1 | + ( x ^ i , k | k S 1 ) T ( P i , k | k S 1 ) T x ^ i , k | k S 1 ) ζ j , k | k S 2 = 1 2 ( d log 2 π log | ( P j , k | k S 2 ) 1 | + ( x ^ j , k | k S 2 ) T ( P j , k | k S 2 ) T x ^ j , k | k S 2 ) ζ k | k c = 1 2 ( d log 2 π log | ( P k | k c ) 1 | + ( x ^ k | k c ) T ( P k | k c ) T x ^ k | k c ) ζ i j , k | k = 1 2 ( d log 2 π log | P i j , k | k 1 | + x ^ i j , k | k T P i j , k | k T x ^ i j , k | k )

References

  1. Reali, F.; Palmerini, G.B.; Farina, A.; Graziano, A.; Giompapa, S. Parametric analysis of ballistic target-tracking problem by multiple model approach. IET Radar Sonar Navig. 2013, 7, 459–469. [Google Scholar] [CrossRef]
  2. Benavoli, A.; Chisci, L.; Farina, A. Tracking of a Ballistic Missile with A-Priori Information. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 1000–1016. [Google Scholar] [CrossRef]
  3. Farina, A.; Ristic, B.; Benvenuti, D. Tracking a ballistic target: comparison of several nonlinear filters. IEEE Trans. Aerosp. Electron. Syst. 2002, 38, 854–867. [Google Scholar] [CrossRef]
  4. Saulson, B.; Chang, K.C. Nonlinear Estimation Comparison for Ballistic Missile Tracking. Opt. Eng. 2004, 43, 1424–1438. [Google Scholar] [CrossRef]
  5. Dodin, P.; Minvielle, P.; Cadre, J.P.L. Estimating the Ballistic Coefficient of a Re-Entry Vehicle. IET Radar Sonar Navig. 2007, 1, 173–183. [Google Scholar] [CrossRef]
  6. Zhao, Z.; Chen, H.; Chen, G.; Kwan, C.; Li, X.R. Comparison of Several Ballistic Target Tracking Filters. In Proceedings of the American Control Conference, Minneapolis, MN, USA, 14–16 June 2006; pp. 2197–2202.
  7. Farina, A.; Benvenuti, D.; Ristic, B. Estimation Accuracy of a Landing Point of a Ballistic Target. In Proceedings of the Fifth International Conference on Information Fusion, Annapolis, MD, USA, 8–11 July 2002; pp. 2–9.
  8. Ristic, B.; Farina, A.; Benvenuti, D.; Arulampalam, M.S. Performance Bounds and Comparison of Nonlinear Filters for Tracking a Ballistic Object on Re-Enty. IEE Proc. Radar Sonar Navig. 2003, 150, 65–70. [Google Scholar] [CrossRef]
  9. Julier, S.J.; Uhlmann, J.K. Unscented Filtering and Nonlinear Estimation. Proc. IEEE 2004, 92, 401–422. [Google Scholar] [CrossRef]
  10. Singh, N.K.; Bhaumik, S.; Bhattacharya, S. Tracking of Ballistic Target on Re-Entry Using Ensemble Kalman Filter. In Proceedings of the Annual IEEE India Conference, Kochi, India, 7–9 December 2012; pp. 508–513.
  11. Tseng, C.H.; Lin, S.F.; Jwo, D.J. Fuzzy Adaptive Cubature Kalman Filter for Integrated Navigation Systems. Sensors 2016, 16, 1167. [Google Scholar] [CrossRef] [PubMed]
  12. Zhu, W.; Wang, W.; Yuan, G. An Improved Interacting Multiple Model Filtering Algorithm Based on the Cubature Kalman Filter for Maneuvering Target Tracking. Sensors 2016, 16, 805. [Google Scholar] [CrossRef] [PubMed]
  13. Barrios, C.; Motai, Y.; Huston, D. Intelligent Forecasting Using Dead Reckoning with Dynamic Errors. IEEE Trans. Ind. Inf. 2016. [Google Scholar] [CrossRef]
  14. Lee, S.J.; Motai, Y.; Choi, H. Tracking Human Motion with Multichannel Interacting Multiple Model. IEEE Trans. Ind. Inf. 2013, 9, 1751–1763. [Google Scholar] [CrossRef]
  15. Himberg, H.; Motai, Y.; Bradley, A. A Multiple Model Approach to Track Head Orientation with Delta Quaternions. IEEE Trans. Cybern. 2013, 43, 90–101. [Google Scholar] [CrossRef] [PubMed]
  16. Chang, K.; Chong, C.Y.; Mori, S. Analytical and Computational Evaluation of Scalable Distributed Fusion Algorithms. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 2022–2034. [Google Scholar] [CrossRef]
  17. Grime, S.; Durrant-Whyte, H. Communication in Decentralized Systems. IFAC Control Eng. Pract. 1994, 2, 849–863. [Google Scholar] [CrossRef]
  18. Hurley, M.B. An Information Theoretic Justification for Covariance Intersection and Its Generalization. In Proceedings of the Fifth International Conference on Information Fusion, Annapolis, MD, USA, 8–11 July 2002; pp. 505–511.
  19. Liggins, M., II; Hall, D.; Llinas, J. Handbook of Multisensor Data Fusion: Theory and Practice; CRC Press: Boca Raton, FL, USA, 2008; pp. 424–427. [Google Scholar]
  20. Martin, T.W.; Chang, K.C. A Distributed Data Fusion Approach for Mobile Ad Hoc Networks. In Proceedings of the 8th International Conference on Information Fusion, Philadelphia, PA, USA, 25–28 July 2005; pp. 1062–1069.
  21. Weng, Y.; Xiao, W.; Xie, L. Diffusion-Based Em Algorithm for Distributed Estimation of Gaussian Mixtures in Wireless Sensor Networks. Sensors 2011, 11, 6297–6316. [Google Scholar] [CrossRef] [PubMed]
  22. Julier, S.J. An Empirical Study into the Use of Chernoff Information for Robust, Distributed Fusion of Gaussian Mixture Models. In Proceedings of the 9th International Conference on Information Fusion, Florence, Italy, 10–13 July 2006; pp. 1–8.
  23. Gunay, M.; Orguner, U.; Demirekler, M. Approximate Chernoff Fusion of Gaussian Mixtures Using Sigma-Points. In Proceedings of the 17th International Conference on Information Fusion (FUSION), Salamanca, Spain, 7–10 July 2014; pp. 1–8.
  24. Üney, M.; Clark, D.E.; Julier, S.J. Distributed Fusion of PHD Filters via Exponential Mixture Densities. IEEE J. Sel. Top. Sign. Proces. 2013, 7, 521–531. [Google Scholar] [CrossRef]
  25. Shin, S.J. Re-Entry Vehicle Tracking with a New Multiple Model Estimation Applicable to Highly Non-Linear Dynamics. IET Radar Sonar Navig. 2015, 9, 581–588. [Google Scholar] [CrossRef]
  26. Farrell, W. Interacting Multiple Model Filter for Tactical Ballistic Missile Tracking. IEEE Trans. Aerosp. Electron. Syst. 2008, 44, 418–426. [Google Scholar] [CrossRef]
  27. Bar-Shalom, Y.; Li, X.R. Estimation with Applications to Tracking and Navigation; John Wiley & Sons: New York, NY, USA, 2001; pp. 453–457. [Google Scholar]
  28. Li, W.; Jia, Y. An Information Theoretic Approach to Interacting Multiple Model Estimation. IEEE Trans. Aerosp. Electron. Syst. 2015, 51, 1811–1825. [Google Scholar] [CrossRef]
  29. Tian, X.; Bar-Shalom, Y. On Algorithms for Asynchronous Track-to-Track Fusion. In Proceedings of the 13th Conference on Information Fusion (FUSION), Edinburgh, UK, 26–29 July 2010; pp. 1–8.
  30. Hall, D.; Chong, C.Y.; Llinas, J.; Liggins, M. Distributed Data Fusion for Network-Centric Operations; CRC Press: Boca Raton, FL, USA, 2012; pp. 105–107. [Google Scholar]
  31. Chang, K.; Chong, C.Y.; Mori, S. On Scalable Distributed Sensor Fusion. In Proceedings of the 11th International Conference on Information Fusion, Cologne, Germany, 30 June–3 July 2008; pp. 1–8.
  32. Runnalls, A.R. Kullback-Leibler Approach to Gaussian Mixture Reduction. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 989–999. [Google Scholar] [CrossRef]
  33. Bromiley, P. Products and Convolutions of Gaussian Probability Density Functions. Available online: http://www.tina-vision.net/docs/memos/2003-003.pdf (accessed on 14 August 2014).
Figure 1. Structure of the IMM algorithm for sensor S1.
Figure 1. Structure of the IMM algorithm for sensor S1.
Sensors 16 01289 g001
Figure 2. Fusion architecture.
Figure 2. Fusion architecture.
Sensors 16 01289 g002
Figure 3. Flow chart of the fusion methodology.
Figure 3. Flow chart of the fusion methodology.
Sensors 16 01289 g003
Figure 4. Ballistic target trajectory.
Figure 4. Ballistic target trajectory.
Sensors 16 01289 g004
Figure 5. Velocity of ballistic target against time.
Figure 5. Velocity of ballistic target against time.
Sensors 16 01289 g005
Figure 6. Aerodynamic drag acceleration against time.
Figure 6. Aerodynamic drag acceleration against time.
Sensors 16 01289 g006
Figure 7. Model probabilities for sensor S1.
Figure 7. Model probabilities for sensor S1.
Sensors 16 01289 g007
Figure 8. Model probabilities for sensor S2.
Figure 8. Model probabilities for sensor S2.
Sensors 16 01289 g008
Figure 9. Estimated value of the ballistic coefficient for sensor S1.
Figure 9. Estimated value of the ballistic coefficient for sensor S1.
Sensors 16 01289 g009
Figure 10. Estimated value of the ballistic coefficient for sensor S2.
Figure 10. Estimated value of the ballistic coefficient for sensor S2.
Sensors 16 01289 g010
Figure 11. Comparison of RMSE in X axis between the proposed Gaussian mixtures (GM) fusion, local IMM and CI.
Figure 11. Comparison of RMSE in X axis between the proposed Gaussian mixtures (GM) fusion, local IMM and CI.
Sensors 16 01289 g011
Figure 12. Comparison of RMSE in Y axis between the proposed Gaussian mixtures (GM) fusion, local IMM and CI.
Figure 12. Comparison of RMSE in Y axis between the proposed Gaussian mixtures (GM) fusion, local IMM and CI.
Sensors 16 01289 g012
Table 1. Comparison of RMSE (m) in X and Y axis (64 s–120 s).
Table 1. Comparison of RMSE (m) in X and Y axis (64 s–120 s).
S1S2CI
XYXYXY
478.9953.8660.31163578.21102
425.3820.7625.81081537.11014
345.4644.8563.2927.4474.9868.6
233.9410488.8786.9404.2710.1
151.9235.1427.7687319531.9
109.1150.3385.2647.5258.2426.8
80.25100.3336.2574.2212.8349
56.4871.63295.6500.6185319.1
42.1667.96273.1429.3163.9294.3
35.566.76266.1410.4160.7280.4

Share and Cite

MDPI and ACS Style

Lu, K.; Zhou, R. Sensor Fusion of Gaussian Mixtures for Ballistic Target Tracking in the Re-Entry Phase. Sensors 2016, 16, 1289. https://doi.org/10.3390/s16081289

AMA Style

Lu K, Zhou R. Sensor Fusion of Gaussian Mixtures for Ballistic Target Tracking in the Re-Entry Phase. Sensors. 2016; 16(8):1289. https://doi.org/10.3390/s16081289

Chicago/Turabian Style

Lu, Kelin, and Rui Zhou. 2016. "Sensor Fusion of Gaussian Mixtures for Ballistic Target Tracking in the Re-Entry Phase" Sensors 16, no. 8: 1289. https://doi.org/10.3390/s16081289

APA Style

Lu, K., & Zhou, R. (2016). Sensor Fusion of Gaussian Mixtures for Ballistic Target Tracking in the Re-Entry Phase. Sensors, 16(8), 1289. https://doi.org/10.3390/s16081289

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