Next Article in Journal
AF-DHNN: Fuzzy Clustering and Inference-Based Node Fault Diagnosis Method for Fire Detection
Previous Article in Journal
NLOS Correction/Exclusion for GNSS Measurement Using RAIM and City Building Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heterogeneous Multiple Sensors Joint Tracking of Maneuvering Target in Clutter

1
Department of Automation, Nanjing University of Science and Technology, No.200, Xiaolingwei Street, Xuanwu District, Nanjing 210094, China
2
Department of Science, Nanjing University of Science and Technology, No.200, Xiaolingwei Street, Xuanwu District, Nanjing 210094, China
*
Author to whom correspondence should be addressed.
Sensors 2015, 15(7), 17350-17365; https://0-doi-org.brum.beds.ac.uk/10.3390/s150717350
Submission received: 12 May 2015 / Revised: 13 July 2015 / Accepted: 13 July 2015 / Published: 17 July 2015
(This article belongs to the Section Physical Sensors)

Abstract

:
To solve the problem of tracking maneuvering airborne targets in the presence of clutter, an improved interacting multiple model probability data association algorithm (IMMPDA-MDCM) using radar/IR sensors fusion is proposed. Under the architecture of the proposed algorithm, the radar/IR centralized fusion tracking scheme of IMMPDA-MDCM is designed to guarantee the observability of the target state. The interacting multiple model (IMM) deals with the model switching. The modified debiased converted measurement (MDCM) filter accounts for non-linearity in the dynamic system models, and reduces the effect of measurement noise on the covariance effectively. The probability data association (PDA) handles data association and measurement uncertainties in clutter. The simulation results show that the proposed algorithm can improve the tracking precision for maneuvering target in clutters, and has higher tracking precision than the traditional IMMPDA based on EKF and IMMPDA based on DCM algorithm.

1. Introduction

Target tracking is an essential requirement for the fire control system of the armed reconnaissance vehicle, which is equipped with a suite of advanced sensors to detect, locate, track, classify and automatically identify targets under all climatic conditions. The sensors fusion system has the superiority over the conventional one with a single sensor in target tracking [1]. Active and passive sensors are mutually independent or complementary to target detection and tracking. The radar is an active sensor, which has narrow beam width and high precision of range measurement. However, it is easy to be interfered with by electromagnetic field. The infrared (IR) sensor is a passive system, which is quite sensitive to atmospheric conditions and has no effect on electromagnetic interference. Furthermore, it has higher precision of angular measurements than radar [2]. The radar/IR fusion system could considerably improve the tracking precision by using their complementary characteristics. However, the measurements of the radar and IR sensors are nonlinear and the target is maneuverable. Therefore, the nonlinear filter for maneuvering target tracking should be researched for radar/IR fusion system.
A promising approach to track a maneuvering target is the interacting multiple model (IMM) algorithm. The IMM is built from a finite number of dynamic models that represent different target behavioral traits [3], which makes it natural to track maneuvering target. A converted measurement IMM filter was proposed for tracking a maneuvering target using radar/IR sensors [4]. A new distributed fusion method of radar/IR tracking system based on separation and combination of the measurements was proposed in [5]. A distributed flow of information fusion for radar/IR compound seeker was established in [6], and the federated filter was used to track the target. An adaptive grid IMM based on modified iterated extended Kalman filter for tracking a maneuvering target using radar/ IR sensors was proposed in [7].
In the radar/IR sensor compound tracking system, the tracking of an airborne target in a cluttered environment might be a challenge due to the several observations for a single airborne target, some tracking measurements do not originate from the airborne target. Therefore, the present study utilizes the probabilistic data association (PDA) filter [8,9] to assign weights to the validated measurements. The PDA filter can extend the tracking capability to a highly cluttered environment. Combining IMM with appropriate data association algorithm can realize maneuvering target tracking in clutters [10,11], such as maximum likelihood probabilistic data association (ML-PDA), IMMPDA, interacting multiple model multiple hypothesis tracking (IMM-MHT) and so on. A ML-PDA algorithm has been shown to be robust in a cluttered environment for a constant velocity target, however, it cannot be applied to the situation where targets undergo maneuvers. An adaptive update rate tracking algorithm based on modified IMMPDA is proposed to avoid tracking loss of maneuvering target tracking in clutters [12]. An interacting multiple model probability data association (IMMPDA) algorithm was proposed to support the navigation and surveillance services of the air traffic management system [13].
The dynamic of target is usually modeled and tracked in the Cartesian coordinates, whereas the measurements are provided in terms of range and angle with respect to the sensor location in the polar coordinates. Therefore, the radar/IR compound tracking becomes a kind of non-liner estimation problem. One solution to this problem is the extended Kalman filter (EKF), but would results in filter divergence [14,15]. The other solution is debiased converted measurement (DCM) Kalman filter [16], which converts the polar measurements to Cartesian coordinates and then filtering in the Cartesian coordinates. In this paper, the IMMPDA algorithm is combined with the modified debiased converted measurement (MDCM) filter to create an IMMPDA-MDCM filter for an airborne maneuvering target tracking in radar/IR fusion system. The Monte Carlo simulation results show that the proposed IMMPDA-MDCM kalman filter (IMMPDA-MDCMKF) algorithm can improve the target tracking precision, credibility and outperform the conventional algorithms.
The remainder of this paper is organized as follows. In Section 2, the sensor measurement model is derived. In Section 3, the time alignment and fusion of radar and IR sensors are derived. The IMMPDA-MDCM algorithm is proposed in Section 4. In Section 5, the simulation results demonstrate the feasibility and precision of the proposed algorithm. Conclusions are drawn in Section 6.

2. The Sensor Measurement Model

Considering an arbitrary maneuvering target in 3D Cartesian coordinates, the geometry measuring relationship between target and radar/IR platform is described in Figure 1.
Figure 1. Geometry measuring relationship between target and radar/IR platform.
Figure 1. Geometry measuring relationship between target and radar/IR platform.
Sensors 15 17350 g001
The range r R k , azimuth θ R k and elevation β R k can be measured by radar, and the radar measurement is Z R k = [ r R k , θ R k , β R k ] T .Let the state vector of target is X k = [ x k , y k , z k , x ˙ k , y ˙ k , z ˙ k , x ¨ k , y ¨ k , z ¨ k ] T . The radar measurement equation is
Z R k = [ x k 2 + y k 2 + z k 2 arctan y k x k arcsin z k x k 2 + y k 2 + z k 2 ] + [ v r R , k v θ R , k v β R , k ] = h R ( X k )+ v R k
where ( x k , y k , z k ) , ( x ˙ k , y ˙ k , z ˙ k ) , and ( x ¨ k , y ¨ k , z ¨ k ) represent target position, velocity, and acceleration components in x, y, and z direction, respectively. v r R k , v θ R k and v β R k are separately independent identically distributed (i.i.d) zero-mean Gaussian white noise, with variance σ r R k 2 , σ θ R k 2 and σ β R k 2 respectively. The measurement noise of Z R k is σ R k 2 = d i a g [ σ r R k 2 , σ θ R k 2 , σ β R k 2 ] .
The measurements of IR sensor is Z I k = [ r I k , θ I k , β I k ] T , the IR measurement equation is
Z I k = [ arctan y k x k arcsin z k x k 2 + y k 2 + z k 2 ] + [ v θ I k v β I k ] = h I ( X k )+ v I k
where v θ I k and v β I k are separately i.i.d zero-mean Gaussian white noise with variance σ θ I k 2 and σ β I k 2 respectively. The measurement noise variance of Z I k is σ I k 2 = [ σ θ I k 2 σ β I k 2 ] T .

3. Data Fusion with Radar and IR Sensors

For the convenience of discussion, we assume that radar and IR sensors lie in the same platform. Compared with the detected target, the distance between them is negligible. Therefore, it can be assumed that the two sensors are located in the same position.
Nowadays, most information fusion algorithms for the compound tracking system are using centralized fusion structure. In this paper, the centralized fusion architecture is used before IMMPDA-MDCM filtering. Because the measurements from radar and IR sensors are independent of each other, each sensor transmits data to the fusion center at different sampling period. Therefore, the measured data from different sensors should be synchronized, and time alignment is needed. The diagram of radar/IR fusion tracking architecture is shown in Figure 2.
Figure 2. Centralized fusion tracking architecture with radar and IR sensors.
Figure 2. Centralized fusion tracking architecture with radar and IR sensors.
Sensors 15 17350 g002

3.1. Time Alignment of Radar and IR

Suppose that the sampling period of radar and IR sensors are TR and TI, respectively, and TR:TI = m:n. The two sensors can be synchronized once at T (the smallest common multiple of TR and TI). Therefore, it is reasonable to select T to be the sampling period of tracking system. In time T, radar and IR sensors have n and m samples, respectively. Common algorithms of time alignment include least square method, interpolation, extrapolation and curve fitting algorithm [17]. Because of the good real-time performance, the least square algorithm is chosen in this paper.
Z ^ R k = c 1 i = 1 n Z R k i + c 2 i = 1 n i Z R k i
Z ^ I k = d 1 i = 1 m Z I k i + d 2 i = 1 m i Z I k i
where c 1 = 2 / n , c 2 = 6 / n ( n + 1 ) , d 1 = 2 / m , d 2 = 6 / m ( m + 1 ) . After alignment, the measurement noise variances of radar and IR sensors are denoted as
σ ^ R k 2 = 2 σ R k 2 ( 2 n + 1 ) n ( n + 1 )
σ ^ I k 2 = 2 σ I k 2 ( 2 m + 1 ) m ( m + 1 )

3.2. The Fusion of Synchronized Data

In this paper, the weighted average algorithm is used for the fusion of azimuth and elevation measurements. In this fusion algorithm, the constrained extremum calculated through Lagrange multiplier algorithm is chosen as the weighted coefficient, and the fusion precision can approximate the optimal. The fusion measurements and variance of azimuth are denoted as
θ k = σ ¯ θ k 2 [ θ ^ R k σ ^ θ R k 2 + θ ^ I k σ ^ θ I k 2 ]
σ ¯ θ k 2 = σ ^ θ R k 2 σ ^ θ I k 2 σ ^ θ R k 2 + σ ^ θ I k 2
The fusion measurements and variance of elevation are denoted as
β k = σ ¯ β k 2 [ β ^ R k σ ^ β R k 2 + β ^ I k σ ^ β I k 2 ]
σ ¯ β k 2 = σ ^ β R k 2 σ ^ β I k 2 σ ^ β R k 2 + σ ^ β I k 2
The time aligned range r ^ R , k and the fusion angles are merged into an augmented measurement vector as
Z k = h ( X k , v k ) = [ r ^ R k θ k β k ] + v k

4. IMMPDA-MDCM Algorithm

The IMM algorithm can estimate the state of a dynamic system with several different models that switch from one to another, and finally get a mixing output. Various nonlinear filtering algorithms can run in the IMM framework [18]. The tracking error of the single IMMPDA algorithm in clutter is large, and obvious error of peak value may appear in the period of target maneuvering. In this paper, the MDCMKF is embedded in IMMPDA architecture for maneuvering target tracking. Assuming there are r models, the target dynamics are modeled in Cartesian coordinates as
X k + 1 i = F k i X k i + G i W k i , i = 1 , , r
Where X k i is the state of the target at time k for ith model, F k i is the transition matrix of ith model, the Gi is the process noise gain matrix. Wi is the mode-dependent process noise sequences with zero mean and covariance Qi. The transition probability from model i to model j is P i j .
P = [ P 11 P 12 P 1 r P 21 P 22 P 2 r P r 1 P r 2 P r r ]

4.1. MDCM Algorithm

In the system of radar/IR compound tacking, the dynamic target is usually modeled and tracked in the Cartesian coordinates, whereas the measurements are provided in terms of range and angle with respect to the radar and IR sensors in the polar coordinates. The debiased converted measurement (DCM) Kalman filter is a popular technique for target tracking. In the spherical coordinate, the true measurements of radar are azimuth angle θ m , elevation angle β m and radial distance r m , with noise variance as σ θ 2 , σ β 2 , σ r 2 , respectively. The average true deviation u k and average true covariance R k of converted measurement are described as [19]
u k = [ u k x , u k y , u k z ] T
R k = [ R k x x R k x y R k x z R k y x R k y y R k y z R k z x R k z y R k z z ]
When measurement in the spherical coordinate is converted to be in Cartesian coordinate, the measurement is modified as
Z c = Z k u k = [ r m cos β m cos θ m r m cos β m sin θ m r m sin β m ] u k
The debiased converted measurement equation in Cartesian coordinates can be described as
Z c = H X k + v k
where H = [ 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 ] .
The covariance of the DCM is a noisy stochastic process with strong correlation to the measurement, the filter update equations are actually coupled and nonlinear, which causes the DCM to lose its unbiasedness. In this paper, a modified DCM (MDCM) filter is derived to suppress this dependence. The MDCM filter can be given as follows
Step 1: Initialization
X ¯ 0 = E [ X 0 ] , P 0 = E [ ( X 0 X ¯ 0 ) ( X 0 X ¯ 0 ) T ]
Step 2: Predict the target state
X ¯ k + 1 | k = F k X ¯ k
P k + 1 | k = F k P k F k T + G k Q k G k T
Step 3: First update of the target state
K ¯ k + 1 = P k + 1 | k H k T S ¯ k 1
S ¯ k = H k P k + 1 | k H k T + R k
X ¯ k + 1 = X ¯ k + 1 | k + K k + 1 ( Z k u k H k X ¯ k + 1 | k )
P ¯ k + 1 = ( I K k + 1 H k ) P k + 1 | k
Step 4: Second update of the target state
Step 4.1: Calculate the error covariance using the first target estimation state ( x ¯ , y ¯ , z ¯ ) and covariance ( σ x ¯ 2 , σ y ¯ 2 , σ z ¯ 2 ) .
R ¯ ( k + 1 )= E [ R | r ¯ , θ ¯ , β ¯ ] = [ R ¯ k + 1 x x R ¯ k + 1 x y R ¯ k + 1 x z R ¯ k + 1 y x R ¯ k + 1 y y R ¯ k + 1 y z R ¯ k + 1 z x R ¯ k + 1 z y R ¯ k + 1 z z ]
where R ¯ k + 1 x x , R ¯ k + 1 y y , R ¯ k + 1 z z , R ¯ k + 1 x y , R ¯ k + 1 x z and R ¯ k + 1 y z are derived by using the radar and IR measurement error. R ¯ k + 1 x x = 0.25 ( r ¯ 2 + δ r ¯ 2 + δ r 2 ) ( 1 + λ θ λ θ ¯ cos 2 θ ¯ ) ( 1 + λ β λ β ¯ cos 2 β ¯ ) 0.25 λ θ 2 λ β 2 ( r ¯ 2 + δ r ¯ 2 ) ( 1 + λ θ ¯ cos 2 θ ¯ ) ( 1 + λ β ¯ cos 2 β ¯ ) , R ¯ k + 1 y y = 0.25 ( r ¯ 2 + δ r ¯ 2 + δ r 2 ) ( 1 λ θ λ θ ¯ cos 2 θ ¯ ) ( 1 + λ β λ β ¯ cos 2 β ¯ ) 0.25 λ θ 2 λ β 2 ( r ¯ 2 + δ r ¯ 2 ) ( 1 λ θ ¯ cos 2 θ ¯ ) ( 1 + λ β ¯ cos 2 β ¯ ) , R ¯ k + 1 z z = 0.5 ( r ¯ 2 + δ r ¯ 2 + δ r 2 ) ( 1 λ β λ β ¯ cos 2 β ¯ ) 0.5 λ β 2 ( r ¯ 2 + δ r ¯ 2 ) ( 1 λ β ¯ cos 2 β ¯ ) , R ¯ k + 1 x y = 0.25 λ θ λ θ ¯ ( r ¯ 2 + δ r ¯ 2 + δ r 2 ) sin 2 θ ¯ ( 1 + λ β λ β ¯ cos 2 β ¯ ) 0.25 λ θ 2 λ β 2 λ θ ¯ ( r ¯ 2 + δ r ¯ 2 ) sin 2 θ ¯ ( 1 + λ β ¯ cos 2 β ¯ ) , R ¯ k + 1 x z = 0.5 λ θ λ β λ θ ¯ λ β ¯ ( r ¯ 2 + δ r ¯ 2 + δ r 2 ) cos θ ¯ sin 2 β ¯ 0.5 λ θ λ β 2 λ θ ¯ λ β ¯ ( r ¯ 2 + δ r ¯ 2 ) cos θ ¯ sin 2 β ¯ , R ¯ k + 1 y z = 0.5 λ θ λ β λ θ ¯ λ β ¯ ( r ¯ 2 + δ r ¯ 2 + δ r 2 ) sin θ ¯ sin 2 β ¯ 0.5 λ θ λ β 2 λ θ ¯ λ β ¯ ( r ¯ 2 + δ r ¯ 2 ) sin θ ¯ sin 2 β ¯ , r ¯ = x ¯ 2 + y ¯ 2 + z ¯ 2 , θ ¯ = arctan ( y ¯ / x ¯ ) , β ¯ = arctan ( z ¯ / x ¯ 2 + y ¯ 2 ) , λ θ ¯ = e δ θ ¯ 2 / 2 , λ β ¯ = e δ β ¯ 2 / 2 , λ θ ¯ = e 2 δ θ ¯ 2 = λ θ ¯ 4 , λ β ¯ = e 2 δ β ¯ 2 = λ β ¯ 4 , δ r ¯ 2 , σ θ ¯ 2 , σ β ¯ 2 are error covariance of r ¯ , θ ¯ and β ¯ , respectively.
Step 4.2: Update the target state
K k + 1 = P k + 1 | k H k T S k 1
S k = H k P k + 1 | k H k T + R ¯ k
X k + 1 = X ¯ k + 1 | k + K k + 1 ( Z k u k H k X ¯ k + 1 | k )
P k + 1 = ( I K k + 1 H k ) P k + 1 | k

4.2. IMMPDA-MDCM Algorithm Principle

One complete cycle of the proposed IMMPDA-MDCMKF comprises four major steps: mixing probabilities calculation, IMMPDA-MDCMKF filtering, model probability update and output mixing. Detailed steps of the proposed algorithm is given as follows
Step 1: Mixing probabilities calculation
μ k i j = p i j μ k i / C k j i , j = 1 , .... , r
where μ k i is the conditional probability of ith model at k, C k j is the normalizing constant.
C k j i = 1 r p i j μ k i    j = 1 , .... , r
Step 2: MDCM filter in clutter
Step 2.1: Input interaction. Computing the input state and covariance matrices of ith model
X k o i = i = 1 r X k i μ k i j   
P k o i = i = 1 r ( P k i + ( X k i X k o i ) ( X k i X k o i ) T ) μ k i j   
Step 2.2: State and covariance prediction
X k + 1 | k i = F k i X k o i
P k + 1 | k i = F k i P k o i ( F k i ) T + G k i Q k i ( G k i ) T
Step 2.3: Validated measurement judgment. The validation region is
d k + 1 i ( ) = ( v k + 1 i ( ) ) T ( S k + 1 i ( ) ) 1 v k + 1 i ( ) ,    = 1 , 2 n
d k + 1 1 ( 1 ) γ | | d k + 1 i ( n ) γ
where v k + 1 i ( ) and S k + 1 i ( ) are innovation vector and innovation covariance at k + 1 of ith model. v k + 1 i ( ) = Z d k + 1 H k + 1 i X k + 1 | k i , Z d k + 1 is a matrix with three rows and n columns, each column represents a set of measurements, n is the number of the measurements. Equation (37) is the validation equation, γ is the threshold corresponding to the gate probability, which can be obtained from Chi-Square tables for a chosen gate probability [20]. Once the ith measurement passes the Chi-Square test in Equation (37), it can be utilized in the rest of the probability data association filter.
Step 2.4: Converted measurement error calculation. Calculating u k + 1 and R k + 1 using Equations (14) and (15).
Step 2.5: Probabilistic data association for each validated measurement.
β k + 1 i ( l ) = e k + 1 i ( l ) b k + 1 i + l = 1 m e k + 1 i ( x ) , l = 1 , 2 , ... , m
β k + 1 i ( 0 ) = b k + 1 i b k + 1 i + l = 1 m e k + 1 i ( x )
where m is the number of validated measurements and associated with the track. β k + 1 i ( l ) is the association probability of the ith target-originated measurement. β k + 1 i ( 0 ) is the association probability of all measurements are not valid. v k + 1 i ( m d c m ) ( l ) is the innovation associated with the lth validated measurement, e k + 1 i ( x ) = exp ( 0.5 ( v k + 1 i ( m d c m ) ( l ) ) T ( S k + 1 i ) 1 v k + 1 i ( m d c m ) ( l ) ) , v k + 1 i ( m d c m ) ( l ) = Z k + 1 l u k + 1 H k + 1 i X k + 1 | k i b k + 1 i = C ( 2 π ) m / 2 | S k + 1 i | 1 / 2 ( 1 P D P G ) / P D . PD and PG are the target detection probability and the gate probability, respectively.
Step 2.6: First update of the target state and covariance
Using the combined innovation to substitute the clutter-free innovation, and calculate the gain matrix, state and covariance updating matrix.
K ^ k + 1 i = P k + 1 | k i ( H k + 1 i ) T ( S ^ k + 1 i ) 1
S ^ k + 1 i = H k + 1 i P k + 1 | k i ( H k + 1 i ) T + R k + 1
X ^ k + 1 i = X k + 1 | k i + K ^ k + 1 i v ˜ k + 1 i
v ˜ k + 1 i = l = 1 m β k + 1 i ( l ) v k + 1 i ( m d c m ) ( l )
P ^ k + 1 i = P k + 1 | k i ( 1 β k + 1 i ( 0 ) ) K ^ k + 1 i S ^ k + 1 i ( K ^ k + 1 i ) T + P ˜ ^ k + 1 i
P ˜ ^ k + 1 i = K ^ k + 1 i ( l = 1 m β k + 1 i ( l ) v k + 1 i ( m d c m ) ( l ) ( v k + 1 i ( m d c m ) ( l ) ) T v ˜ k + 1 i ( v ˜ k + 1 i ) T ) ( K ^ k + 1 i ) T
Step 2.7: Second update of the target state and covariance. Calculating R ¯ k + 1 using Equation (25).
K k + 1 i = P k + 1 | k i ( H k + 1 i ) T ( S k + 1 i ) 1
S k + 1 i = H k + 1 i P k + 1 | k i ( H k + 1 i ) T + R ¯ k + 1
X k + 1 i = X k + 1 | k i + K k + 1 i v ˜ k + 1 i
v ˜ k + 1 i = l = 1 m β k + 1 i ( l ) v k + 1 i ( m d c m ) ( l )
P k + 1 i = P k + 1 | k i ( 1 β k + 1 i ( 0 ) ) K k + 1 i S k + 1 i ( K k + 1 i ) T + P ˜ k + 1 i
P ˜ k + 1 i = K k + 1 i ( l = 1 m β k + 1 i ( l ) v k + 1 i ( m d c m ) ( l ) ( v k + 1 i ( m d c m ) ( l ) ) T v ˜ k + 1 i ( v ˜ k + 1 i ) T ) ( K k + 1 i ) T
Step 3: Model probability update
μ k + 1 i = 1 c ¯ Λ k + 1 i c i , c ¯ = i = 1 r Λ k + 1 i c i
where Λ k + 1 i is the likelihood function of ith model in IMMPDA-MDCM.
Λ k + 1 i = P D V k + 1 i ( m + 1 ) m [ b k + 1 i + l = 1 m e k + 1 i ( l ) ]
V k + 1 i = c n z γ n z / 2 | S k + 1 i | 1 / 2
c n z = { π n z / 2 ( n z / 2 ) ! ,                 n z = 2 , 4 , ... 2 ( n z + 1 2 ) ! π ( n z 1 ) / 2 ( n z + 1 ) ! , n z = 1 , 3 , ...
Step 4: Output Mixing.
The final target state estimation and covariance matrix are combined from all of the models
X k + 1 = i = 1 r X k + 1 i μ k + 1 i
P k + 1 = i = 1 r ( P k + 1 i + ( X k + 1 i X k + 1 ) ( X k + 1 i X k + 1 ) T ) μ k + 1 i

5. Simulation and Results

The following example of tracking a highly maneuvering unmanned aerial vehicle is considered. The scenario of a highly maneuvering airborne target tracking is defined as follows: the sampling rate is T = 0.1 s, the target makes five accelerating maneuver with linear segments connecting it. The initial position of the target is (10,000, 6000, 4000) m, and the velocity is (−300, −300, −100) m/s. In the first period of 1–5 s, it flies linearly by constant velocity. From 6–10 s, it makes an accelerating maneuver with (20, 50, 0) m/s2. From 11–15 s, it flies with (5, 25, 0) m/s2. From 16–20 s, it flies with (5, −25, 0) m/s2. From 21–25 s, it flies with (−25, −50, 0) m/s2. From 26–30 s, it flies with (0, 25, 0) m/s2. At last, it flies linearly from 31–35 s by constant velocity. The trajectory of target is shown in Figure 3.
Figure 3. Trajectory of target.
Figure 3. Trajectory of target.
Sensors 15 17350 g003
In this paper, the target dynamics are modeled in Cartesian coordinates as Equation (12). The constant velocity (CV) model and Singer model are put into the IMM frame, The CV model is used to describe the basic motion of the target, the Singer model is used to describe target maneuver.
The state transition matrix and noise gain matrix of CV model are defined as
F k C V = [ 1 0 0 T 0 0 0 0 0 0 1 0 0 T 0 0 0 0 0 0 1 0 0 T 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ] , G k C V = [ T 2 / 2 0 0 0 T 2 / 2 0 0 0 T 2 / 2 T 0 0 0 T 0 0 0 T 0 0 0 0 0 0 0 0 0 ]
The state transition matrix and noise gain matrix of Singer model are defined as reference [21].
F k S = [ 1 0 0 T 0 0 ϕ 17 0 0 0 1 0 0 T 0 0 ϕ 28 0 0 0 1 0 0 T 0 0 ϕ 39 0 0 0 1 0 0 ϕ 47 0 0 0 0 0 0 1 0 0 ϕ 58 0 0 0 0 0 0 1 0 0 ϕ 69 0 0 0 0 0 0 e - α x T 0 0 0 0 0 0 0 0 0 e - α y T 0 0 0 0 0 0 0 0 0 e - α z T ] , G k S = [ G 1 G 2 G 3 ]
where, ϕ 17 = ( α x T 1 + e α x T ) / α x 2 , ϕ 28 = ( α y T 1 + e α y T ) / α y 2 , ϕ 39 = ( α z T 1 + e α z T ) / α z 2 , ϕ 47 = ( 1 e α x T ) / α x , ϕ 58 = ( 1 e α y T ) / α y , ϕ 69 = ( 1 e α z T ) / α z , G 1 = d i a g ( [ σ a x [ 1 α x T e α x T + ( T 2 α x 2 / 2 ) ] / α x 3 σ a y [ 1 α y T e α y T + ( T 2 α y 2 / 2 ) ] / α y 3 σ a z [ 1 α z T e α z T + ( T 2 α z 2 / 2 ) ] / α z 3 ] ) , G 2 = d i a g ( [ σ a x ( α x T + e α x T 1 ) / α x 2 σ a y ( α T + e α y T 1 ) / α y 2 σ a z ( α T + e α z T 1 ) / α z 2 ] ) , G 3 = d i a g ( [ σ a x ( 1 e α x T ) / α x σ a y ( 1 e α y T ) / α y σ a z ( 1 e α z T ) / α z ] ) . α x = α y = α z = α , α = 0.1 is the reciprocal of the manoeuver time constant. σ a x , σ a y , σ a z are standard deviation of maneuver acceleration in x, y and z direction.
σ ( a x , a y , a z ) 2 = a ( x , y , z ) max 2 ( 1 + 4 p max p 0 ) / 3
a ( x , y , z ) max ( a x max = 25 , a y max = 50 , a z max = 0 ) is maximum acceleration of target, p max = 0.5 is the maximum probability of acceleration or deceleration, p 0 = 0.5 is the probability of without acceleration.
The measuring period of radar T R = 0.01 s , the measuring period of IR T I = 0.005 s . The range measurement variance of radar is σ r R = 100 , the azimuth and elevation angle measurement variance of radar are σ θ R = σ β R = 0.02 . The azimuth and elevation angle measurement variance of IR sensor are σ θ I = σ β I = 0.002 . The initial prior probability of the two models are μ 1 = 0.5 and μ 2 = 0.5 . That is to say, CV model has the same chance to be selected in the initialization. γ = 16 , λ = 4 e 5 , P G = 0.997 , P D = 1 . The variances of the process noise of two models are Q C V = 0.05625 I 3 m/s 2 , Q S = 3 I 3 m/s 2 . I 3 is the identity matrix of three dimensions. Considering the different process noise level, the transition probability of the system model is chosen as
p = [ 0.99 0.01 0.01 0.99 ]
The azimuth and elevation comparison after fusion are shown in Figure 4 and Figure 5. The comparison standard deviation of azimuth and elevation after time alignment and fusion are shown in Table 1.
Figure 4. The comparison of azimuth.
Figure 4. The comparison of azimuth.
Sensors 15 17350 g004
Figure 5. The comparison of elevation.
Figure 5. The comparison of elevation.
Sensors 15 17350 g005
Table 1. Standard deviation comparison.
Table 1. Standard deviation comparison.
Before Time AlignmentAfter Time AlignmentFusion
RadarIRRadarIR
Azimuth (rad)0.020.0020.020.00260.0026
Elevation (rad)0.020.0020.020.00260.0026
The tracking performances of proposed IMMPDA-MDCM algorithm, IMMPDA-DCM and IMMPDA-EKF are compared via 100 Monte Carlo simulations. All the algorithms are implemented using a personal computer (Windows 7 2009, Intel Core2 Duo CPU, 2.94 GHz, 4.0 GB of RAM, and MATLAB R2012a programming environment). The results of the root mean square error (RMSE) and runtime test of the target’s position for the three algorithms are shown in Table 2. Figure 6, Figure 7 and Figure 8 show the obtained position estimation error of three algorithms in x, y, and z direction, respectively.
Figure 6. The comparison of position error in x direction.
Figure 6. The comparison of position error in x direction.
Sensors 15 17350 g006
Figure 7. The comparison of position error in y direction.
Figure 7. The comparison of position error in y direction.
Sensors 15 17350 g007
Figure 8. The comparison of position error in z direction.
Figure 8. The comparison of position error in z direction.
Sensors 15 17350 g008
Table 2. The RMSE and runtime comparison of three different algorithms.
Table 2. The RMSE and runtime comparison of three different algorithms.
X (m)Y (m)Z (m)t (s)
IMMPDA-EKF13.36824.37925.4764.735
IMMPDA-DCM12.91121.80316.9650.767
IMMPDA-MDCM8.89113.20512.2031.396
The proposed algorithm can calculate the statistic characteristics of converted measurement errors and make the covariance to be much less noisy. As can be seen in Figure 6, Figure 7 and Figure 8 , the proposed IMMPDA-MDCM algorithm has the highest tracking precision than IMMPDA-DCM and IMMPDA-EKF algorithm, which is consistent with the results in Table 2. The total position tracking error of the proposed algorithm is reduced by 34.22% and 46.81% compared to IMMPDA-DCM and IMMPDA-EKF, respectively. In the simulation, the two update of IMMPDA-MDCM will increase thecalculation time; and the computational cost of IMMPDA-MDCM is longer than IMMPDA-DCM, but shorter than IMMPDA-EKF.

6. Conclusions

In this paper, an interacting multiple model probability data association algorithm based on modified debiased converted measurement filter (IMMPDA-MDCM) is proposed, which is capable of adaptively tracking the maneuvering airborne target. The polar measurements of radar and IR measurements are time aligned, fused and converted to Cartesian coordinates before they are applied to IMMPDA-MDCM algorithm. In the IMMPDA-MDCM algorithm, the covariance of the converted measurement is recalculated using the estimated target position information. Therefore, the proposed algorithm can reduce the effect of measurement noise on the covariance effectively. By abandoning the extended Kalman filter framework and using MDCM filter in the proposed algorithm, the linearization errors of the measurement model are avoided, and the good tracking precision is achieved with decreasing the computational complexity. Monte Carlo simulation results verify that the proposed algorithm outperforms IMMPDA-DCM and IMMPDA-EKF in terms of filtering unbiasedness and precision. The proposed algorithm is an effective algorithm for maneuvering target tracking in clutter, which can increase warfare airplane’s concealment and survival capacity.

Acknowledgments

This work is partially supported by National Nature Science Foundation of China (61301217, 61473153), Research Fund for the Doctoral Program Higher Education of China (20123219120043), Natural Science Foundation of Jiangsu Province (BK20131352) and A Project Funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Author Contributions

This work is a product of inputs and work of all the authors to the research concept and experiment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rao, G.A.; Mahulikar, S.P. New criterion for aircraft susceptibility to infrared guided missiles. Aerosp. Sci. and Technol. 2005, 9, 701–712. [Google Scholar] [CrossRef]
  2. Bar-Shalom, Y.; Li, X.R. Multitarget-Multisensor Tracking: Principles and Techniques; YBS: London, UK, 1995. [Google Scholar]
  3. Mazor, E.; Averbuch, A.; Bar-shalom, Y.; Dayan, J. Interacting multiple model in target tracking: A survey. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 103–124. [Google Scholar] [CrossRef]
  4. Yin, J.H.; Cui, B.Z.; Wang, Y.F. A novel maneuvering target tracking algorithm for radar/infrared sensors. Chin. J. Electron. 2010, 19, 752–756. [Google Scholar]
  5. Wang, Q.C.; Wang, W.F. Tracking method based on separation and combination of the measurements for radar and IR fusion System. J. Syst. Eng. Electron. 2009, 20, 241–246. [Google Scholar]
  6. Xu, Y.K.; Liang, X.G. Information fusion for radar/infrared compound seeker based on federated filter. Int. J. Dig. Conten. Technol. Appl. 2011, 5, 218–229. [Google Scholar]
  7. Wu, P.L.; Li, X.X.; Zhang, L.Z.; Bo, Y.M. Tracking algorithm with radar/IR sensors using a novel adaptive grid IMM. IET Sci. Meas. Technol. 2014, 8, 270–276. [Google Scholar] [CrossRef]
  8. Bar-Shalom, Y.; Daum, F.; Huang, J. The probabilistic data association filter estimation in the presence of measurement origin uncertainty. IEEE Control Syst. 2009, 29, 82–100. [Google Scholar] [CrossRef]
  9. Habtemariam, B.; Tharmarasa, R.; Thayaparan, T.; Mallick, M.; Kirubarajan, T.A. Multiple-detection joint probabilistic data association filter. IEEE J. Sel. Top. Sign. Process. 2013, 7, 461–471. [Google Scholar] [CrossRef]
  10. Musicki, D.; Suvorova, S. Tracking in clutter using IMM-IPDA-based algorithms. IEEE Trans. Aerosp. Electron. Syst. 2008, 44, 111–126. [Google Scholar] [CrossRef]
  11. Blackman, S.S. Multiple hypothesis tracking for multiple target tracking. IEEE Aerosp. Electron. Syst. 2004, 19, 5–18. [Google Scholar] [CrossRef]
  12. Chen, T.; He, Z.S.; Li, Y.X. Adaptive update rate tracking based on modified IMM-PDA. IEEE CIE Int. Conf. Radar 2011, 1, 804–807. [Google Scholar]
  13. Jan, S.S.; Kao, Y.C. Radar tracking with an interacting multiple model and probabilistic data association filter for civil aviation applications. Sensors 2013, 13, 6636–6650. [Google Scholar] [CrossRef] [PubMed]
  14. Alipakzad, M. Kalman filter design for time delay systems. WSEAS Trans. Syst. 2012, 11, 551–560. [Google Scholar]
  15. Julier, S.J.; Uhlmann, J.K.; Durrant-Whyte, H.F. A new method for the nonlinear transformation of means and covariances in filters and estimators. IEEE Trans. Autom. Control 2000, 45, 477–482. [Google Scholar] [CrossRef]
  16. Zhao, Z.L.; Li, X.R.; Vesselin, P.J. Best linear unbiased filtering with nonlinear measurements for target tracking. IEEE Trans. Aerosp. Electron. Syst. 2004, 40, 1324–1336. [Google Scholar] [CrossRef]
  17. Sanders-Reed, J.N. Error propagation in two-sensor 3D position estimation. Opt. Eng. 2001, 40, 627–636. [Google Scholar] [CrossRef]
  18. Lerro, D.; Bar-Shalom, Y. Interacting multiple model tracking with target amplitude feature. IEEE Trans. Aerosp. Electron. Syst. 1993, 29, 494–509. [Google Scholar] [CrossRef]
  19. Wu, P.L.; Zhang, L.Z.; Zhang, X.Y. The design of DSP/FPGA based maneuvering target tracking system. WSEAS Trans. Circ. Syst. 2014, 13, 75–84. [Google Scholar]
  20. Kirubarajan, T.; Bar-shalom, Y.; Blair, W.D.; Watson, G.A. IMMPDAF for radar management and tracking benchmark with ECM. IEEE Trans. Aerosp. Electron. Syst. 1998, 34, 1115–1134. [Google Scholar] [CrossRef]
  21. Wu, P. L.; Li, X.X.; Bo, Y.M. Iterated square root unscented Kalman filter for maneuvering target tracking using TDOA measurements. Int. J. Control Autom. Syst. 2013, 11, 761–767. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Wu, P.; Li, X.; Kong, J.; Liu, J. Heterogeneous Multiple Sensors Joint Tracking of Maneuvering Target in Clutter. Sensors 2015, 15, 17350-17365. https://0-doi-org.brum.beds.ac.uk/10.3390/s150717350

AMA Style

Wu P, Li X, Kong J, Liu J. Heterogeneous Multiple Sensors Joint Tracking of Maneuvering Target in Clutter. Sensors. 2015; 15(7):17350-17365. https://0-doi-org.brum.beds.ac.uk/10.3390/s150717350

Chicago/Turabian Style

Wu, Panlong, Xingxiu Li, Jianshou Kong, and Jiale Liu. 2015. "Heterogeneous Multiple Sensors Joint Tracking of Maneuvering Target in Clutter" Sensors 15, no. 7: 17350-17365. https://0-doi-org.brum.beds.ac.uk/10.3390/s150717350

Article Metrics

Back to TopTop