Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016

arXiv:1603.03450v1 [stat.ME] 9 Mar 2016

Multisensor–Multitarget Bearing–Only Sensor
Registration
Ehsan Taghavi, R. Tharmarasa
and T. Kirubarajan
McMaster University
Hamilton, Ontario, Canada
email: {taghave, tharman, kiruba}@mcmaster.ca
Mike McDonald
Defence Research and Development Canada
Ottawa, Ontario, Canada
email: mike.mcdonald@drdc-rddc.gc.ca

Abstract
Bearing–only estimation is one of the fundamental and challenging problems in target tracking. As in the case of radar tracking,
the presence of offset or position biases can exacerbate the challenges in bearing–only estimation. Modeling various sensor biases
is not a trivial task and not much has been done in the literature specifically for bearing–only tracking. This paper addresses the
modeling of offset biases in bearing–only sensors and the ensuing multitarget tracking with bias compensation. Bias estimation
is handled at the fusion node to which individual sensors report their local tracks in the form of associated measurement reports
(AMR) or angle-only tracks. The modeling is based on a multisensor approach that can effectively handle a time–varying number
of targets in the surveillance region. The proposed algorithm leads to a maximum likelihood bias estimator. The corresponding
Cramér–Rao Lower Bound to quantify the theoretical accuracy that can be achieved by the proposed method or any other algorithm
is also derived. Finally, simulation results on different distributed tracking scenarios are presented to demonstrate the capabilities
of the proposed approach. In order to show that the proposed method can work even with false alarms and missed detections,
simulation results on a centralized tracking scenario where the local sensors send all their measurements (not AMRs or local
tracks) are also presented.
Index Terms
Bias estimation, bearing–only tracking, target motion analysis, triangulation, filtering, multisensor–multitarget

I. I NTRODUCTION
Multisensor–multitarget bearing–only tracking is a challenging problem with many applications [4], [5], [13], [27]. Some
applications of bearing–only tracking are in maritime surveillance using sonobuoys, underwater target tracking using sonar
and passive ground target tracking using Electronic Support Measures (ESM) or Infra–red Search and Track (IRST) sensors.
In such applications, it is of interest to find the target position as well as any biases that may affect estimation performance.
From the early works in [1], [2], [36] to the more recent works in [29] and the references therein, the focus has been only
on tracking the targets based on measurements from a bearing–only sensor. However, due to the limitations of single sensor
bearing–only tracking, i.e., due to the need for own–ship maneuvers for the observability of state parameters [28], the issue of
biases in passive single sensor tracking has not been addressed in the literature. The main focus of this paper is multisensor
bearing–only tracking in the presence of biases. In multisensor bearing-only tracking, observability is no longer a major issue.
However, in the case of port–starboard ambiguity, the problem of observability was discussed in detail in [10]. Besides, the
presence of sensor biases that are often unaccounted for can degrade the estimation results significantly. Most of the works on
bias estimation have been about radar tracking (see [20], [42], [43] and the references therein) or using other measurements
besides bearing information [6], [7]. For example, when the elevation information is available, one can estimate the offset
biases as in [7], [14].
With the objective of providing a combined bias estimation and target tracking algorithm that is both theoretically sound
and practical, the problem of multisensor bearing-only multitarget tracking is considered in this paper. Having more than one
passive sensor in the surveillance region ensures the observability of the state parameters, i.e., position and velocity of the
target, without the need for maneuvers [33], [35], [44]. One of the issues that can complicate bearing–only tracking is the bias
in the sensors. For example, in maritime surveillance using sonobuoys, which are usually dropped from an aircraft or thrown
from a ship close to an area of interest, the exact locations of the sonobuoys are not known. This leads to position biases [25].
This is also an isuue in modern systems such as autonomous underwater vehicles (AUV) [9]. In addition, the impact with the
water surface and the waves can result in systemic offset biases [6]. In wide area surveillance using airborne IRST sensors,
uncertain platform motion can contribute to biases as well. Offset bias can be modeled as an additive constant term affecting
the measurement equation and the sensor position uncertainty can be modeled using a random walk [31].

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016
Negligible biases can be treated as residual errors. This residual error can be used in the form of additional uncertainty
in the measurements later in the filtering step. However, if offset biases are larger than the noise standard deviation of the
bearing–only measurements, a mirror of the target’s bearing is sent to fusion center instead of its actual value, which will
result in totally erroneous estimates. Further, these erroneous measurements can worsen state estimation results when fused
with measurements from similarly biased sensors. In order to benefit from the information available from multiple bearing–only
sensors with offset bias, one needs to model, estimate and then correct the biases. This is precisely the motivation for this
paper.
The focus of this paper is offset bias estimation. In order to model offset biases, one can transform the measurement
space of sensor data to Cartesian coordinates followed by covariance matrix transformation. This transformation will make it
possible to find an exact model for the biases that can be used in bias estimation and correction. However, the full position
information is not available in a single bearing–only measurement. The process of measurement transformation is done by
paring measurements from different sensors in the surveillance region. That is, this transformation is done through triangulation
[21]. With the new pseudo–measurement in Cartesian coordinates, position and velocity of the targets can be estimated over
time as new measurements are generated from paired sensors at subsequent times [3], [19], [24], [35]. Assuming that these
state estimates also carry the effects of offset biases, it is possible to find such biases, if any, and correct them. In the case
of bistatic passive sensors, one can use the methods in [37], [44]. However, previous work on multisensor bearing–only bias
estimation is still limited. The method presented here gives a comprehensive analysis of offset bias modeling in multisensor
passive bearing–only sensors.
The proposed method gives an exact model for bearing–only biases in Cartesian coordinates. In addition, the formulation
of an appropriate likelihood function enables the use of maximum likelihood estimators to find the biases. Also, the proposed
model is robust against large sensor noise standard deviation. Finally, as shown through simulation results, large bearing biases
can be estimated accurately, which leads to correspondingly accurate target state estimation results.
The goal of this paper is to present a step–by–step approach for designing a multisensor–multitarget tracking system based
on biased bearing–only measurements and give a practical solution to the problem of bias estimation. In Section II, a detailed
model of the multisensor bearing–only estimation problem with bias is given. Section III is devoted to modeling the offset biases
in Cartesian coordinates. In Section IV, a practical solution for bias estimation is proposed. Section V presents the derivation of
Cramér–Rao lower bounds. Simulation results are shown in Section VI along with discussions on different scenarios. Finally,
Section VII ends the paper with conclusions.
II. B EARING –O NLY E STIMATION P ROBLEM
Bearing–only sensors with operating ranges of hundreds of meters to a few kilometers are one of the most crucial sensors in
maritime or ground surveillance applications. These sensors can actively or passively detect the directions of arrival of signals
emitted by the targets of interest. While underwater surveillance is the common application of bearing–only tracking, it is also
used in surface and air target tracking. For example, ESM and IRST sensors also use bearing–only sensors for tracking. As
shown in Figure 1, bearing–only sensors can be on the own–ship or deployed separately in the surveillance region. Moreover,
they can operate under different environmental conditions as shown in Figure 1 [8], [26].

Sonar sensor

𝜌2
𝜌1

spherical spreading loss

cylindrical spreading

Fig. 1: Bearing–only sensors and signal propagation types [26]
The bearing–only measurement from passive sensors is written as
zs (k)

=

θs (k) + bs + ws (k),

(1)

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016
where zs (k) is the direction of arrival at sensor s, θs is the true bearing of the target, bs is the constant bias in the measurements
of sensor s and ws is an additive zero–mean white Gaussian noise with variance σθ2 . It is assumed that there are S bearing–only
sensors in the surveillance region at positions (xs , ys ) for s = 1, 2, ..., S, and, they record targets’ bearings at time instants
k = 1, ..., K. Note that there is no index to denote target ID, but wherever such clarification is needed, it will be included.
In this paper, bias estimation is handled at the fusion node to which individual sensors report their local tracks in the form
of associated measurement reports (AMR) [16], [46] or angle-only tracks. That is, a distributed tracking system is considered
as in the case of [31], [32], [42]. However, the difference is that in these earlier works local full-state tracks were sent to the
fusion node for bias estimation whereas AMRs or bearing-only tracks are sent in the present case where the local sensors do
not have full observability due to bearing-only measurements. As in [31], [32], [42], working with tracks or AMRs obviates
the need at the fusion center to address false alarms and missed detections that are inevitable at the local sensors although
ghost tracks may be present among local tracks. However, in order to show that the proposed method can work even with false
alarms and missed detections in a centralized tracking system, in one experiment in Section VI we assume that the local sensors
send all their measurements (including false alarms and missed detections instead of AMRs or local bearing-only tracks) and
evaluate the performance of the proposed method. In [37] and [44], where bias estimation at measurement level (rather than at
track level as in [31], [32], [42]) is considered, false alarms and missed detections are not addressed at all. In [30], a joint data
association and bias estimation method was proposed for linear measurement models, which is not applicable for bearing-only
systems. A general case of multistatic passive radar system with false alarms and missed detections was considered in [17]
and [15], but, the bias problem was not addressed.
The goal of bearing–only tracking is to find the bias in each sensor and estimate each target’s position as accurately as
possible based on the model given in (1), either as decoupled parameters or jointly [11], [45]. Due to the computational burdens
of joint tracking and parameter estimation methods [22], a decoupled bias and state estimation is presented in this paper.
It is assumed that each target is following the Discretized Continuous White Noise Acceleration (DCWNA) or the nearly
constant velocity (NCV) model [5]. As a result, a target’s state vector in 2D Cartesian coordinate is given by

T
x(k) vx (k) y(k) vy (k)
x (k) =
(2)
with (x(k), y(k)) being the position and (vx (k), vy (k)) being the velocity. The motion model can further be defined as
x (k + 1)

=

F (k)x (k) + ν(k)

(3)

where the state transition matrix is


1
 0
F (k) , 
 0
0

T
1
0
0

0
0
1
0


0
0 

T 
1

and the covariance matrix of ν(k) which is an additive zero–mean white Gaussian noise vector is


Qx 0
Q =
0 Qy

(4)

(5)

in which Qx and Qy are defined as

Qx

=


Qy

=

1 3
3T
1 2
2T
1 3
3T
1 2
2T

1 2
2T



T
1 2
2T

T

q̄x

(6)

q̄y

(7)



where q̄x and q̄y are noise intensities with dimension m2 /s3 [4].
In order to transform the measurements from polar to Cartesian coordinates, it is initially assumed that S is even and that
the sensors are paired into S2 one–to–one pairs. Note that this constraint is relaxed in Section IV. For a pair of sensors with a
single target in the common field of view, the best estimate for the location of the target, independent of its previous location,
can be obtained through triangulation [14], [21], [38]. The triangulated estimates of the target position at time k using sensor
pair (i, j), ignoring measurement noise, are given by
x̂ij (k)

=

yj − yi + xi tan (θi (k) + bi )
xj tan (θj (k) + bj )
−
tan (θi (k) + bi ) − tan (θj (k) + bj ) tan (θi (k) + bi ) − tan (θj (k) + bj )
(8)

ŷij (k)

=

yj tan (θi (k) + bi ) − yi tan (θj (k) + bj ) (xi − xj ) tan (θi (k) + bi ) tan (θj (k) + bj )
+
tan (θi (k) + bi ) − tan (θj (k) + bj )
tan (θi (k) + bi ) − tan (θj (k) + bj )
(9)

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016
where x̂ij (k) and ŷij (k) are the X and Y Cartesian estimates, respectively.

T
In addition, by defining the stacked covariance matrix in the bearing–only coordinate for the stacked measurement θi (k) θj (k)
as
"
#
σθ2i (k)
0
Rij (k) =
(10)
0
σθ2j (k)

T
one can calculate the covariance matrix of the transformed vector xij (k) yij (k)
as

−1
−1
xy
(k) =
J(k)T (Rij (k)) J(k)
Rij

(11)

T
xij (k) yij (k)
and

T
Ji (k) Jj (k)
J(k) =

where J(k) is the Jacobian matrix with respect to



(12)

Further, the elements of J(k) can be written as
#
"
(xij (k) − xi )
− (yij (k) − yi )
Ji (k) =
2
2 ,
2
2
(xij (k) − xi ) + (yij (k) − yi )
(xij (k) − xi ) + (yij (k) − yi )

(13)

and
"
Jj (k)

=

− (yij (k) − yj )
2

(xij (k) − xj ) + (yij (k) − yj )

2

,

#

(xij (k) − xj )
2

2

(xij (k) − xj ) + (yij (k) − yj )

(14)

If no bias estimation is needed, these pseudo–measurements in Cartesian coordinates can be used in a Kalman filter with their
associated covariance matrices to recursively estimate the target’s position [12].
III. B IAS M ODELING IN C ARTESIAN C OORDINATES
As discussed in Section II, bearing–only biases are in the form of the additive constants. Although additive constant biases
have already been dealt with in the case of radar measurements [42] and [43], it is not possible to generate similar pseudo–
measurements with bearing–only data directly. One way to formulate pseudo–measurements with biases is to model in Cartesian
coordinates. In this section, a step–by–step procedure to model the biases and to separate them from original track positions
in Cartesian coordinates is given. In Section IV, the pseudo–measurement generation is discussed in detail.
In Section II, the process of mapping from bearing–only measurements to Cartesian was given. In order to model the biases,
one can start with separating the bias terms in (8) and (9) from the original track position in Cartesian coordinates. This
separation of bias terms provides the necessary information to create a pseudo–measurement that properly addresses the biases
as in the case of radar measurements. Once the pseudo–measurements are generated, it is possible to estimate the biases and
remove them. The process of finding the bias terms that contribute to (8) and (9) starts with expanding the tan (·) function as
tan (α + β)

=

tan (α) + tan (β)
1 − tan (α) tan (β)

(15)

Applying (15) to (8) and (9), and separating the terms related to the bias from those related to the target position will give a
new set of equations to define the position of the target in Cartesian coordinates. To make the parameter separation easier to
follow, the common terms are defined and named first. For the common terms in x̂ij (k) and ŷij (k), one can define
D

=

tan (θi (k)) − tan (θj (k))

(16)

B

=

1 + tan (θi (k)) tan (θj (k))

(17)

Further, define the following for the terms in x̂ij (k):
Dx

= xi tan (θi (k)) − xj tan (θj (k)) + (yj − yi )

(18)

Bxi

= − (yj − yi ) tan (θi (k)) + xj tan (θi (k)) tan (θj (k)) + xi

(19)

Bxj
Bxij

= − (yj − yi ) tan (θj (k)) − xi tan (θi (k)) tan (θj (k)) + xj

(20)

= xj tan (θi (k)) − xi tan (θj (k)) + (yj − yi ) tan (θi (k)) tan (θj (k))

(21)

Similarly, for the terms in ŷij (k), define
Dy

= yj tan (θi (k)) − yi tan (θj (k)) + (xj − xi ) tan (θi (k)) tan (θj (k))

Byi

=

(xj − xi ) tan (θj (k)) + yi tan (θi (k)) tan (θj (k)) + yj

(23)

Byj
Byij

=

(xj − xi ) tan (θi (k)) − yj tan (θi (k)) tan (θj (k)) − yi

(24)

= yi tan (θi (k)) − yj tan (θj (k)) + (xj − xi )

(22)

(25)

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016
With these factorizations, bias terms can be separated from the target state values in Cartesian coordinates. It can be seen that
T
the vector x̂ij (k) ŷij (k)
can be written as1



 u
xij (k) + βx (θi (k), θj (k), bi , bj )
x̂ij (k)
=
u
yij
(k) + βy (θi (k), θj (k), bi , bj )
ŷij (k)
(26)
where
xuij (k)

=

u
yij
(k)

=

Dx
D
Dy
D

(27)
(28)

and the bias terms can be written in as
βx (θi (k), θj (k), bi , bj )

βy (θi (k), θj (k), bi , bj )

=

=

Bxi tan (bi ) + Bxj tan (bj ) + Bxij tan (bi ) tan (bj )
D + D tan (bi ) + B tan (bj ) + D tan (bi ) tan (bj )
Dx B tan (bi ) − Dx B tan (bj ) + Dx D tan (bi ) tan (bj )
− 2
D + BD tan (bi ) − BD tan (bj ) + D2 tan (bi ) tan (bj )
Byi tan (bi ) + Byj tan (bj ) + Byij tan (bi ) tan (bj )
D + D tan (bi ) + B tan (bj ) + D tan (bi ) tan (bj )
Dy B tan (bi ) − Dy B tan (bj ) + Dy D tan (bi ) tan (bj )
− 2
D + BD tan (bi ) − BD tan (bj ) + D2 tan (bi ) tan (bj )

(29)

(30)

Note that in the above formulations it is assumed that the biases and the true bearings are available. In practice, only the noisy
or estimated values are available. Assuming small values for the bias and noise terms, one can use the above formulation
without significant loss in accuracy. Similar assumptions has been made in the previous works on bias estimation [31], [32],
[42]. A technical discussion on the range of bias and noise values for which the above formulation is valid is given in Section
VI.
IV. B EARING –O NLY T RACKING AND R EGISTRATION
Bearing–only sensor registration is a challenging problem in target tracking that has been addressed in [34], [37], [41].
In order to find the biases and correct the measurements, one should first look into the observability of the bias variables.
Note that in (26), provided that the target is not on the line that connects the two sensors used in the triangulation or in the
vicinity of one of the sensors, the state parameters are observable [4]. In addition, if there are two pairs of sensors tracking
the same target, the biases become observable as it is shown in IV-A. To estimate the biases decoupled from the state vector, a
pseudo–measurement that can address the bias vector directly must be defined. In this section, a new formulation is proposed
to create a pseudo–measurement that can be used for bias estimation with bearing–only data. The key requirement of this
method in order to ensure observability of all parameters is to have at least two sensor pairs in the surveillance region. In the
following, two practical scenarios that can be expanded to a more general formulation to handle varying number of sensors
and targets are discussed in detail.
A. Pseudo–measurement of bearing–only measurements
To handle practical bearing–only scenarios, two different cases are analyzed here. In each case, a separate pseudo–measurement
model is proposed along with its associated covariance matrix. The main idea is to use two different position approximations
to create a pseudo–measurement as discussed below.
1) Four–sensor (or any even number of sensors) case: Assuming that there are two pairs of sensors in the surveillance
region, a vector of nonlinear pseudo–measurements can be defined by subtracting the target positions based on the pairs (i, j)
and (m, n) as


x̂ij (k) − x̂mn (k)
zb (k) =
+ w(k)
(31)
ŷij (k) − ŷmn (k)
where w(k) is the additive zero–mean white Gaussian noise associated with the pseudo–measurement and its covariance matrix
is defined as
Rw (k)
1 The

superscript “u” indicates that the parameter is unbiased.

=

xy
xy
Rij
(k) + Rij
(k)

(32)

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016
Using the fact that in the absence of bias and noise terms, measurements from any two sensors point to the same target location
regardless of the sensor locations, the pseudo–measurement zb (k) can be written as


βx (bi , bj ) − βx (bm , bn )
zb (k) =
+ w(k)
(33)
βy (bi , bj ) − βx (bm , bn )
for two uncorrelated pairs of sensors. This can be applied to any even number of sensors.
2) Three–sensor (or any odd number of sensors) case: In this case, one must create two position approximations from
triangulation to be able to create a pseudo–measurement for biases. Since there are only three sensors, the possible pairs are
(i, j) and (i, m). Then, the pseudo–measurement can be approximated by


x̂ij (k) − x̂im (k)
zb (k) =
+ w(k)
(34)
ŷij (k) − ŷim (k)
where w(k) is approximately additive zero–mean white Gaussian noise associated with the pseudo–measurement and its
covariance matrix is defined as
xy
xy
= Rij
(k) + Rim
(k)

Rw (k)

(35)

Because of the correlation between the two tracks from three sensors in Cartesian coordinates, the noise is not white anymore
and this formulation is only an approximation.
As in the case for four sensors, the pseudo–measurement zb (k) can be written as


βx (bi , bj ) − βx (bi , bm )
zb (k) =
+ w(k)
(36)
βy (bi , bj ) − βx (bi , bm )
Note that for simplicity, the arguments θi (k), θj (k), θm (k) and θn (k) have been dropped from (33) and (36). With an arbitrary
odd number of sensors, one sensor need to be paired with two other, resulting in in a similar approximation. As for the more
general case of n arbitrary sensors, methods similar to [42] can be adopted to handle the situation.
B. Batch maximum–likelihood estimator
To apply a batch estimator for bias estimation, one needs to form a likelihood function. Following (35) and (33), and
assuming that the noise is white, zero–mean and Gaussian, the likelihood function of the bias parameters given two pairs of
sensors is


1
1
T
exp
−
z
(k)
−
h(b)
p (zb (k) | b) =
b
1
−
2
2π |Rw (k)| 2

−1
(Rw (k)) [zb (k) − h(b)]
(37)
where the nonlinear function h(b) of the bias vector is given by


βx (bi , bj ) − βx (bm , bn )
h(b) =
βy (bi , bj ) − βx (bm , bn )

(38)

Assuming independence over time, one can write the likelihood function over k = 1, ..., K as
p (Zb | b)

=

K
Y

p (zb (k) | b)

(39)

k=1

where
Zb

{zb (1), zb (2), . . . , zb (K)}

=

(40)

Finally, the vector b̂ that maximizes the likelihood function can be written as
b̂

=

arg max p (Zb | b)
b

(41)

The above assumes that there is only one target, but it can be extended to the multitarget case using the stacked parameter
estimator.

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016
V. C RAM ÉR –R AO L OWER BOUND FOR B EARING ONLY BIAS ESTIMATION
This section is devoted to the calculation of the Cramér–Rao Lower Bound (CRLB) on the estimation accuracy of bias
parameters by using the pseudo–measurements introduced in Subsection IV-A. Note that based on (38), the measurement
equation for target t at time k is
 t

βx (bi , bj ) − βxt (bm , bn )
t
hb (k) =
+ wt (k)
(42)
βyt (bi , bj ) − βxt (bm , bn )
Assuming that there are N targets in the surveillance region and that the bias parameters are constant over time k = 1, ...K,
the stacked measurement vector can be written as
Z

= h(b) + w

(43)

where


zb1 (1)T

=



h1b (1)T

=



Z =
h(b)
w

w1 (1)T

···

zbN (1)T

···

T
hN
b (1)

···

···

zb1 (K)T

···

h1b (K)T

wN (1)T

Further, the covariance matrix of the noise vector w is
 1
R (1) · · ·
R = diag

RN (1)

···
···

···
···

w1 (K)T

···

R1 (K) · · ·

zbN (K)T

T

T

wN (K)T
RN (K)

(44)

T

T
hN
b (K)



(45)
(46)

(47)

The covariance matrix of an unbiased estimator b̂ is bounded from below by [4]


T 
b̂ − b b̂ − b
E
≥ J−1

(48)

where J is the Fisher Information Matrix (FIM) given by
n
o
T
J = E [∇b log p(Y | b)] [∇b log p(Y | b)]
|b=b̃
n
o
T
= E [∇b λ] [∇b λ]
|b=b̃

(49)

in which b̃ is the true value of the bias vector b, p(Z | b) is the likelihood function of b, λ = − ln p(Y | b) and ∇ is
gradient operator. Based on the model for the stacked measurement vector in (43), one can define the Jacobian matrix of htb (k)
evaluated at the true value b [39] as

T
H̃bt (k) = 5b htb (k)T
(50)
Then, defining
H̃b

=



H̃b1 (1)T

···

H̃bN (1)T

···

H̃b1 (K)T

···

H̃bN (K)T

T

(51)

one can write
J

=

n
o
−1
H̃b
E H̃T
bR

(52)

VI. S IMULATION R ESULTS
To evaluate how the proposed method performs in practical scenarios, a series of simulations is presented in this section.
The implementation details on bias estimation, filtering and fusion are also discussed. Simulation results on different scenarios
are given with discussions on the advantages and disadvantages of the proposed bias estimator.
A. Motion models and measurement generation
A tracking scenario with four bearing–only sensors and sixteen targets is considered as shown in Figure 2. It is assumed
that all sensors are synchronized and that the bias ranges are between −0.05 rad and 0.05 rad. The standard deviation of
measurement noise is σθ = 0.0261 rad = 1.5◦ for target bearing measurements, which is higher than what was previously
assumed in the literature [6]. The true motion of the targets is modeled using the DCWNA or the NCV model [4] with
qx = qy = 0.001m2 /s3 . In the local trackers, DCWNA and Continuous Wiener Process Acceleration (CWPA) are used with
qx = qy = 0.1m2 /s3 to create different scenarios for the simulation. In all scenarios, the sampling time is T = 1s. To validate
the proposed bias model and to quantify its performance at different bias values, the biases are set to both positive and negative
values in different ranges as follows:

T
0.04 rad −0.02 rad 0.03 rad −0.02 rad
btest1 =
(53)

T
−0.04 rad −0.02 rad −0.03 rad −0.02 rad
btest2 =
(54)

T
0.04 rad 0.02 rad 0.03 rad 0.02 rad
btest3 =
(55)

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016

10000
Target locations
Sensor locations

8000

y (m)

6000

4000

2000

0

−2000
0

0.5

1
x (m)

1.5

2
x 10

4

Fig. 2: Initial locations of the targets and sensors
TABLE I: Parameter settings for the Genetic Algorithm
Parameter
Lower bound
Upper bound
Number of generations
Tolerance value

TABLE II: Comparison of

TABLE III: Comparison of

Value
−0.05 rad
0.05 rad
100
1 × 10−15

p
CRLB {[b]i } and RMSE of GA output for bias estimation of btest1 in a distributed system
Bias parameter

q

CRLB [b]i (rad)

RMSE (rad) of the GA bias estimate

bi
bj
bm
bn

0.071 × 10−3
0.071 × 10−3
0.537 × 10−3
0.462 × 10−3

3.32 × 10−3
3.28 × 10−3
3.92 × 10−3
1.45 × 10−3

p

CRLB {[b]i } and RMSE of the GA output for bias estimation of btest2 in a distributed system

Bias parameter

q

CRLB [b]i (rad)

bi
bj
bm
bn

10−3

0.0904 ×
0.0941 × 10−3
0.4243 × 10−3
0.3787 × 10−3

RMSE (rad) of the GA bias estimate
1.72 × 10−3
1.81 × 10−3
2.95 × 10−3
1.73 × 10−3

B. Bias estimation and the Genetic Algorithm
In this paper, the Genetic Algorithm (GA) [18] is used to solve the optimization problem in (41). The Genetic Algorithm is an
efficient optimization algorithm for highly nonlinear objective functions [40] that is widely used in different applications [18].
Note that although the GA is a batch ML estimator, the length of the window can be varied depending on user criteria to meet
the real–time requirements. The parameters used in the simulations are shown in Table I. The algorithms were implemented
TM
on a computer with Intel R Core i7-3720Qm 2.60GHz processor and 8GB RAM.
C. Bias estimation: Four–sensor distributed problem
In this scenario, all four sensors defined earlier are used to implement the GA. Four out of sixteen targets are used for
performance evaluation. AMRs or local bearing–only tracks are collected over 100 time steps and the GA is applied to the
whole data in batch mode. The GA is run with the settings in Table I and the final results are gathered after the termination of
the GA. Then, the estimated bias vector is used over 100 Monte Carlo runs to calculate the Root Mean Square Error (RMSE)
for
p comparison. As the benchmark, the CRLB is also calculated based on the derivations in Section V. The RMSE values and
CRLB {[b]i } of the ML estimates with the three different
p sets of bias parameters are shown in Tables II, III and IV.
Although there is a difference between the RMSE and CRLB {[b]i }, the RMSE results are nearly an order of magnitude
smaller than the bias values, which indicates that any correction made based on the estimated biases will result in better
position estimates. Note that the CRLB in bearing–only tracking problems can be overly optimistic and may even approach

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016
TABLE IV: Comparison of

p
CRLB {[b]i } and RMSE of GA output for bias estimation of btest3 in a distributed system
Bias parameter

q

CRLB [b]i (rad)

RMSE (rad) of the GA bias estimate

10−3

0.0973 ×
0.0954 × 10−3
0.5347 × 10−3
0.4880 × 10−3

bi
bj
bm
bn

2.73 × 10−3
3.75 × 10−3
2.43 × 10−3
2.66 × 10−3

zero (i.e., perfect estimates) in a network of bearing–only sensors [23]. Thus, the difference between the theoretical CRLB
and the empirical RMSE is not of major concern.
To show how much the proposed bias estimation method can help in correcting the target tracks, another simulation is
conducted. In this simulation, it is assumed that the tracker has access to the final estimated bias vector (output of the GA)
and then a Kalman filter is run with the bias estimates in hand. To use the estimated bias parameters, one should, first, correct
the bearing–only measurements with the estimated values. The correction must be done both in the measurement vector and
its associated covariance matrix. Since the estimated biases do not have the covariance information, a scaled version of the the
calculated CRLB of the bias parameters is used instead. The scaling factor can be determined through experiments. Then, the
tracker can be run with these corrected measurements to find the position and velocity estimates of all targets in the surveillance
region. The position RMSE of the original tracks before correction, the RMSE of the corrected estimates and the Cramér–Rao
lower bounds are shown in Figures 3 and 4.
3

Position RMSE (m)

10

RM SEOriginal
RM SECorrected
RM
p SEBiasf ree
√ ΣOriginal
pΣCorrected
ΣBiasf ree
2

10

1

10

20

40
60
Time Step k

80

100

Fig. 3: Position RMSE with distributed tracking for corrected and original tracks of target 2 (set 1)
As shown in Figures 3 and 4, for case 1, Figures 5 and 6 for case 2 and Figures 7 and 8 for case 3, the position error
is reduced significantly in terms of the position RMSE. This demonstrates the effectiveness of the bias model proposed and
the ML estimation algorithm, i.e., the GA, that is applied to the data. Although, according to the results, the correction factor
varies based on the sensor–target orientation, the corrected track and its associated covariance do follow the bias–free values
accurately, which demonstrates the capability of the proposed algorithm in estimating the biases. To process a batch of data
with K = 100, the computational time is 44.1s in MATLAB.
D. Bias estimation: Three–sensor distributed problem
To show the accuracy of bias estimation when there are only three sensors in the surveillance region, the same Genetic
Algorithm is used to solve the ML problem in (41). The primary issue with three sensors is that one of them, here sensor
i, is used to pair with both sensors j and m. This leads to correlation between the two tracks over time. Assuming that the
correlation is negligible, the same method is applied to estimate the biases. Further, for the case of three sensors, the bias
values are set to

T
0.04 rad 0.02 rad 0.03 rad
btest =
(56)

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016

3

Position RMSE (m)

10

RM SEOriginal
RM SECorrected
RM
p SEBiasf ree
√ ΣOriginal
pΣCorrected
ΣBiasf ree
2

10

1

10

20

40
60
Time Step k

80

100

Fig. 4: Position RMSE with distributed tracking for corrected and original tracks of target 3 (set 1)
3

Position RMSE (m)

10

RM SEOriginal
RM SECorrected
RM
p SEBiasf ree
√ ΣOriginal
pΣCorrected
ΣBiasf ree
2

10

1

10

20

40
60
Time Step k

80

100

Fig. 5: Position RMSE with distributed tracking for corrected and original tracks of target 2 (set 2)

to be able to distinguish between the original and corrected tracks. The constraints on the lower and upper bounds of the
optimization algorithm are set to −0.05 rad and 0.05 rad, respectively. The tolerance value for the GA is also set to 1 × 10−15 .
Figures 9 and 10 show the result for position estimates in Cartesian coordinates.
The corrected tracks are sill better in terms of RMSE, which means that the estimated biases are accurate enough in spite
of the correlation.
E. Real–time window–based Genetic Algorithm
It is important for the proposed method to be able to work in real–time. For this purpose, the GA can be set to run, in each
iteration, for a specific window size or duration. The settings of the GA for real time scenarios are given in Table V.
The final estimates and the population matrix in one window can be used as the initial conditions for the next window.
Thus, the estimates of the biases can be used to correct the measurements at the end of processing each window of data. The
results of the simulations using this approach for the case of four sensors (test 1) are given in Figures 11 and 12.
As shown in Figures 11 and 12, the GA is still able to find the biases with a smaller window size and fewer generations.
Note that updating the biases with smaller window sizes enables the use of methods similar to [42], [43] for arbitrary number

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016

3

Position RMSE (m)

10

RM SEOriginal
RM SECorrected
RM
p SEBiasf ree
√ ΣOriginal
pΣCorrected
ΣBiasf ree
2

10

1

10

20

40
60
Time Step k

80

100

Fig. 6: Position RMSE with distributed tracking for corrected and original tracks of target 3 (set 2)

3

Position RMSE (m)

10

RM SEOriginal
RM SECorrected
RM
p SEBiasf ree
√ ΣOriginal
pΣCorrected
ΣBiasf ree
2

10

1

10

20

40
60
Time Step k

80

100

Fig. 7: Position RMSE with distributed tracking for corrected and original tracks of target 2 (set 3)

TABLE V: Parameter settings for real time genetic algorithm
Parameter name
Lower bound
Upper bound
Number of generations
Tolerance value
Window size

Value
−0.05 rad
0.05 rad
50
1 × 10−15
10

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016

3

Position RMSE (m)

10

RM SEOriginal
RM SECorrected
RM
p SEBiasf ree
√ ΣOriginal
pΣCorrected
ΣBiasf ree
2

10

1

10

20

40
60
Time Step k

80

100

Fig. 8: Position RMSE with distributed tracking for corrected and original tracks of target 3 (set 3)
3

Position RMSE (m)

10

RM SEOriginal
RM SECorrected
RM
p SEBiasf ree
√ ΣOriginal
pΣCorrected
ΣBiasf ree
2

10

1

10

20

40
60
Time Step k

80

100

Fig. 9: Position RMSE of corrected and original tracks for the three–sensor distributed tracking case (Target 2)

of sensors in the surveillance region. With this setting, biases can be updated every 9.25 s in MATLAB. This simple example
shows that even when the processing time is a crucial parameter in the design, one can still handle bias estimation in real–time
using the proposed method.
F. Bias estimation with false measurements: Four–sensor centralized problem
To demonstrate how the proposed method performs in the presence of false alarms and missed associations in a centralized
fusion framework where local sensors send all their measurements instead of AMRs or bearing–only tracks, a simulation is
presented in this section. The probability mass function of the number of false alarms or clutter points in surveillance volume
V as a function of their spatial density λ is defined as
m

(λV )
(57)
m!
where m is the number of false alarms [5]. With bearing–only measurements, the volume is V = 2π (rad). It is assumed
that the average number of false alarms per unit volume in a scan, i.e., λ, is 0.5. Also, it is assumed that PD = 0.7 for each
sensor. Note that in this centralized case, the local sensors send all measurements (rather than AMRs or local bearing–only
µFA (m)

= e−λV

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016

3

Position RMSE (m)

10

RM SEOriginal
RM SECorrected
RM
p SEBiasf ree
√ ΣOriginal
pΣCorrected
ΣBiasf ree
2

10

1

10

20

40
60
Time Step k

80

100

Fig. 10: Position RMSE of corrected and original tracks for the three–sensor distributed tracking case (Target 3)
3

Position RMSE (m)

10

RM SEOriginal
RM SECorrected
RM
p SEBiasf ree
√ ΣOriginal
pΣCorrected
ΣBiasf ree
2

10

1

10

20

40
60
Time Step k

80

100

Fig. 11: Position RMSE for corrected and original tracks for the four–sensor distributed tracking case and window size of 10
(Target 2)

tracks) to the fusion node. Because the proposed method is a batch estimator and uses measurements from both sensor pairs to
create a pseudo–measurement for bias estimation, false tracks are often removed prior to generating the pseudo–measurement
vector, which then is sent to the bias estimator. Typically, false tracks do not exist for more than a few time steps as they
are dependent on all four sensors creating false alarms at the same time steps, in the same region, and for a reasonably long
interval of time.
To show the accuracy of bias estimation
p in a centralized system, the same Genetic Algorithm is used to solve the ML
problem in (41). The RMSE values and CRLB {[b]i } of the ML estimates with the bias parameters as defined in (53) are
shown in Table VI. Note that the CRLB values are optimistic because they do not factor in the false alarms or the missed
detections and that the ML estimator does not factor in the false alarms or the missed detections explicitly. A comprehensive
centralized bias estimator is under development. The focus of this paper is the decentralized one.
VII. C ONCLUSIONS
In this paper, a new mathematical model for bearing–only bias estimation in distributed tracking systems was proposed.
This model was based on triangulation using the associated measurement reports or local bearing–only tracks from different

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016

3

Position RMSE (m)

10

RM SEOriginal
RM SECorrected
RM
p SEBiasf ree
√ ΣOriginal
pΣCorrected
ΣBiasf ree
2

10

1

10

20

40
60
Time Step k

80

100

Fig. 12: Position RMSE for corrected and original tracks for the four–sensor distributed tracking case and window size of 10
(Target 3)
TABLE VI: Comparison of CRLB and GA output for bias estimation of btest1 with measurement origin uncertainty in a
centralized system (λ = 0.5 and PD = 0.7)
Bias parameter

CRLB

GA bias estimate

bi
bj
bm
bn

0.2123 × 10−3
0.2123 × 10−3
0.7674 × 10−3
0.5779 × 10−3

7.419 × 10−3
12.85 × 10−3
6.143 × 10−3
4.352 × 10−3

sensor pairs. It was shown that the proposed bias model has the advantage of being practical in scenarios with multiple
sensors. In particular, the proposed algorithm is effective when the sensor noise level and bias values are high. In addition,
previously proposed algorithms were dependent on target–sensor maneuvers and/or limited to certain noise levels. The new bias
model can handle any type of target–sensor motion and it is effective against 5◦ of offset bias in each sensor and uncertainty
levels up to 1.5◦ of noise standard deviation, which is higher than what was assumed in the literature previously. Also, the
proposed method can handle false alarms and missed detections in a centralized architecture. That is, the proposed algorithm
is practical in scenarios with realistic sensor parameter values. Finally, a batch ML estimator was proposed to solve the bias
estimation problem along with simulation results. A comprehensive centralized bias estimation algorithm with data association
for bearing–only sensors is in progress.
R EFERENCES
[1] V. J. Aidala, “Kalman filter behavior in bearings–only tracking applications”, IEEE Transactions on Aerospace and Electronic Systems, vol. 1, pp.
29–39, 1979.
[2] V. J. Aidala and S. E. Hammel, “Utilization of modified polar coordinates for bearings–only tracking”, IEEE Transactions on Automatic Control, vol.
28, no. 3, pp. 283–294, 1983.
[3] M. S. Arulampalam, B. Ristic, N. Gordon, and T. Mansell, “Bearings–only tracking of manoeuvring targets using particle filters”, EURASIP Journal
on Applied Signal Processing, pp. 2351–2365, 2004.
[4] Y. Bar-Shalom, X. R. Li, and T. Kirubarajan, Estimation with Applications to Tracking and Navigation: Theory, Algorithms and Software, Wiley, NY,
2001.
[5] Y. Bar-Shalom, P. K. Willett, and X. Tian, Tracking and Data Fusion: A Handbook of Algorithms, YBS Publishing, Storrs, CT, 2011.
[6] D. Belfadel, R. W. Osborne, and Y. Bar-Shalom, “Bias estimation for optical sensor measurements with targets of opportunity”, 16th International
Conference on Information Fusion, pp. 1805–1812, Istanbul, Turkey, July 2013.
[7] D. Belfadel, R. W. Osborne, and Y. Bar-Shalom, “A minimalist approach to bias estimation for passive sensor measurements with targets of opportunity”,
Proc. SPIE Conference on Signal and Data Processing of Small Targets, Vol. 8857, 2013.
[8] D. A. Blank and A. E. Bock, Introduction to Naval Engineering, Naval Institute Press, Annapolis, MD, U.S.A, 2005.
[9] P. Braca, R. Goldhahn, G. Ferri, and K. LePage, “Distributed information fusion in multistatic sensor networks for underwater surveillance”, IEEE
Sensors Journal, 2015.
[10] P. Braca, P. Willett, K. Lepage, S. Marano, and V. Matta, “Bayesian tracking in underwater wireless sensor networks with port-starboard ambiguity”,
IEEE Transactions on Signal Processing, vol. 62, no. 7, pp. 1864–1878, April 2014.
[11] M. Bugallo, T. Lu, and P. Djuric, “Bearings–only tracking with biased measurements”, 2nd IEEE International Workshop on Computational Advances
in Multi-Sensor Adaptive Processing, pp. 265–268, 2007.
[12] D. E. Catlin, Estimation, Control, and the Discrete Kalman Filter, vol. 71, Springer Science & Business Media, 2012.

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016
[13] Y. T. Chan and S. W. Rudnicki, “Bearings–only and Doppler–bearing tracking using instrumental variables”, IEEE Transactions on Aerospace and
Electronic Systems, vol. 28, no. 4, pp. 1076–1083, 1992.
[14] H. Chen and F. Lian, “Bias estimation for multiple passive sensors”, International Conference on Measurement, Information and Control (MIC), vol.
2, pp. 1081–1084, May 2012.
[15] S. Choi, D. F. Crouse, P. Willett, and S. Zhou, “Approaches to Cartesian data association passive radar tracking in a DAB/DVB network”, IEEE
Transactions on Aerospace and Electronic Systems, vol. 50, no. 1, pp. 649–663, 2014.
[16] R. L. Cooperman, “Tactical ballistic missile tracking using the interacting multiple model algorithm”, Fifth International Conference on Information
Fusion, vol. 2, pp. 824–831, Annapolis, Maryland, USA, 2002.
[17] M. Daun, U. Nickel, and W. Koch, “Tracking in multistatic passive radar systems using DAB/DVB-T illumination”, Signal Processing, vol. 92, no. 6,
pp. 1365–1386, 2012.
[18] L. Davis, Handbook of Genetic Algorithms, vol. 115, Van Nostrand Reinhold, New York, 1991.
[19] A. Farina, “Target tracking with bearings–only measurements”, Signal Processing, vol. 78, no. 1, pp. 61–78, 1999.
[20] S. Fortunati, A. Farina, F. Gini, A. Graziano, M. S. Greco, and S. Giompapa, “Least squares estimation and Cramér–Rao type lower bounds for relative
sensor registration process”, IEEE Transactions on Signal Processing, vol. 59, no. 3, pp. 1075–1087, 2011.
[21] G. Hendeby, R. Karlsson, F. Gustafsson, and N. Gordon, “Recursive triangulation using bearings–only sensors”, IEE Seminar on Target Tracking:
Algorithms and Applications, pp. 3–10, Birmingham, England, March 2006.
[22] M. A. Hopkins and H. F. Vanlandingham, “Pseudo-linear identification: Optimal joint parameter and state estimation of linear stochastic MIMO systems”,
IEEE American Control Conference, pp. 1301–1306, 1988.
[23] P. R. Horridge and M. L. Hernandez, “Performance bounds for angle–only filtering with application to sensor network management”, 6th International
Conference on Information Fusion, pp. 695–703, Cairns, Queensland Australia, July 2003.
[24] R. A. Iltis and K. L. Anderson, “A consistent estimation criterion for multisensor bearings–only tracking”, IEEE Transactions on Aerospace and
Electronic Systems, vol. 32, no. 1, pp. 108–120, 1996.
[25] K. Johansson, K. Jöred, and P. Svensson, “Submarine tracking using multi–sensor fusion and reactive planning for the positioning of passive sonobuoys”,
Hydroakustik, vol. 97, 1997.
[26] N. Z. Kolev, Sonar Systems, InTech, Croatia, 2011.
[27] T. R. Kronhamn, “Bearings–only target motion analysis based on a multihypothesis Kalman filter and adaptive own–ship motion control”, IEE
Proceedings–Radar, Sonar and Navigation, vol. 145, no. 4, pp. 247–252, 1998.
[28] B. La Scala and M. Morelande, “An analysis of the single sensor bearings–only tracking problem”, 11th International Conference on Information
Fusion, pp. 1–6, Cologne, Germany, June 2008.
[29] P. H. Leong, S. Arulampalam, T. A. Lamahewa, and T. D. Abhayapala, “A Gaussian–sum based cubature Kalman filter for bearings–only tracking”,
IEEE Transactions on Aerospace and Electronic Systems, vol. 49, no. 2, pp. 1161–1176, 2013.
[30] Z. Li, S. Chen, H. Leung, and E. Bosse, “Joint data association, registration, and fusion using EM–KF”, IEEE Transactions on Aerospace and Electronic
Systems, vol. 46, no. 2, pp. 496–507, 2010.
[31] X. Lin, Y. Bar-Shalom, and T. Kirubarajan, “Exact multisensor dynamic bias estimation with local tracks”, IEEE Transactions on Aerospace and
Electronic Systems, vol. 1, no. 40, pp. 576–590, 2004.
[32] X. Lin, Y. Bar-Shalom, and T. Kirubarajan, “Multisensor multitarget bias estimation for general asynchronous sensors”, IEEE Transactions on Aerospace
and Electronic Systems, vol. 41, no. 3, pp. 899–921, 2005.
[33] M. Mallick and T. Kirubarajan, “Multi-sensor single target bearing–only tracking in clutter”, NASA STI/Recon Technical Report N, vol. 3, pp. 15779,
2001.
[34] D. W. McMichael and N. N. Okello, “Maximum likelihood registration of dissimilar sensors”, The First Australian Data Fusion Symposium, pp. 31–34,
November 1996.
[35] D. Mušicki, “Bearings–only multi-sensor maneuvering target tracking”, Systems & Control Letters, vol. 57, no. 3, pp. 216–221, 2008.
[36] S. C. Nardone, A. G. Lindgren, and K. F. Gong, “Fundamental properties and performance of conventional bearings–only target motion analysis”, IEEE
Transactions on Automatic Control, vol. 29, pp. 775–787, 1984.
[37] N. Okello and B. Ristic, “Maximum likelihood registration for multiple dissimilar sensors”, IEEE Transactions on Aerospace and Electronic Systems,
vol. 39, no. 3, pp. 1074–1083, July 2003.
[38] Y. Qi, Z. Jing, and S. Hu, “Modified maximum likelihood registration based on information fusion”, Chinese Optics Letters, vol. 5, no. 11, pp. 639–641,
2007.
[39] B. Ristic, S. Arulampalam, and N. Gordon, Beyond the Kalman Filter: Particle Filters for Tracking Applications, vol. 685, Artech House, Boston, MA,
2004.
[40] K. C. Sharman and G. D. McClurkin, “Genetic algorithms for maximum likelihood parameter estimation”, International Conference on Acoustics,
Speech, and Signal Processing, pp. 2716–2719, Glasgow, Scotland, May 1989.
[41] Q. Song and Y. He, “A real–time registration algorithm for passive sensors with TOA and angle biases”, 3rd International Congress on Image and
Signal Processing, vol. 9, pp. 4170–4173, October 2010.
[42] E. Taghavi, R. Tharmarasa, T. Kirubarajan, and Y. Bar-Shalom, “A practical bias estimation algorithm for multisensor–multitarget tracking”, IEEE
Transaction on Aerospace and Electronic Systems, Accepted for publication, Ocotber 2015.
[43] E. Taghavi, R. Tharmarasa, T. Kirubarajan, and Y. Bar-Shalom, “Bias estimation for practical distributed multiradar-multitarget tracking systems”, 15th
International Conference on Information Fusion, pp. 1304–1311, Istanbul, Turkey, July 2013.
[44] B. Xu and Z. Wang, “Biased bearings–only parameter estimation for bistatic system”, Journal of Electronics (China), vol. 24, no. 3, pp. 326–331, 2007.
[45] B. Xu, Z. Wu, and Z. Wang, “On the Cramér–Rao lower bound for biased bearings–only maneuvering target tracking”, Signal Processing, vol. 87, no.
12, pp. 3175–3189, 2007.
[46] J. Yosinski, N. Coult, and R. Paffenroth, “Network–centric angle only tracking”, SPIE Optical Engineering+ Applications, pp. 74450O–74450O,
International Society for Optics and Photonics, 2009.

Accepted for publication in IEEE Transactions on Aerospace and Electronics Systems, March 2016
Ehsan Taghavi received the M.Sc. degree in communication engineering in 2012 from Chalmers University of Technology, Gothenburg, Sweden, where he worked on particle filter smoother. He is currently pursuing the Ph.D. degree in computational science and
engineering at McMaster University, Hamilton, Canada. His research interests include estimation theory, scientific computing, signal
processing, parameter estimation, mathematical modeling and algorithm design.

Ratnasingham Tharmarasa was born in Sri Lanka in 1975. He received the B.Sc.Eng. degree in electronic and telecommunication
engineering from University of Moratuwa, Sri Lanka in 2001, and the M.A.Sc and Ph.D. degrees in electrical engineering from
McMaster University, Canada in 2003 and 2007, respectively.
From 2001 to 2002 he was an instructor in electronic and telecommunication engineering at the University of Moratuwa, Sri Lanka.
During 2002-2007 he was a graduate student/research assistant in ECE department at the McMaster University, Canada. Currently
he is working as a Research Associate in the Electrical and Computer Engineering Department at McMaster University, Canada. His
research interests include target tracking, information fusion and sensor resource management.

Thiagalingam Kirubarajan (S’95M’98SM’03) was born in Sri Lanka in 1969. He received the B.A. and M.A. degrees in electrical
and information engineering from Cambridge University, England, in 1991 and 1993, and the M.S. and Ph.D. degrees in electrical
engineering from the University of Connecticut, Storrs, in 1995 and 1998, respectively.
Currently, he is a professor in the Electrical and Computer Engineering Department at McMaster University, Hamilton, Ontario.
He is also serving as an Adjunct Assistant Professor and the Associate Director of the Estimation and Signal Processing Research
Laboratory at the University of Connecticut. His research interests are in estimation, target tracking, multisource information fusion,
sensor resource management, signal detection and fault diagnosis. His research activities at McMaster University and at the University
of Connecticut are supported by U.S. Missile Defense Agency, U.S. Office of Naval Research, NASA, Qualtech Systems, Inc.,
Raytheon Canada Ltd. and Defense Research Development Canada, Ottawa. In September 2001, Dr. Kirubarajan served in a DARPA
expert panel on unattended surveillance, homeland defense and counterterrorism. He has also served as a consultant in these areas
to a number of companies, including Motorola Corporation, Northrop-Grumman Corporation, Pacific-Sierra Research Corporation,
Lockhead Martin Corporation, Qualtech Systems, Inc., Orincon Corporation and BAE systems. He has worked on the development of a number of engineering
software programs, including BEARDAT for target localization from bearing and frequency measurements in clutter, FUSEDAT for fusion of multisensor data
for tracking. He has also worked with Qualtech Systems, Inc., to develop an advanced fault diagnosis engine.
Dr. Kirubarajan has published about 100 articles in areas of his research interests, in addition to one book on estimation, tracking and navigation and two
edited volumes. He is a recipient of Ontario Premiers Research Excellence Award (2002).

Michael McDonald received a B.Sc (Hons) degree in Applied Geophysics from Queens University in Kingston, Canada in 1986 and
a M.Sc. degree in Electrical Engineering in 1990, also from Queen’s University. He received a Ph.D in Physics from the University
of Western Ontario in London, Canada in 1997. He was employed at ComDev in Cambridge, Canada from 1989 through 1992 in
their space science and satellite communications departments and held a post-doctoral position in the Physics department of SUNY
at Stony Brooke from 1996 through 1998 before commencing his current position as Defence Scientist in the Radar Systems section
of Defence Research and Development Canada, Ottawa, Canada. His current research interests include the application of STAP
processing and nonlinear filtering to the detection of small maritime and land targets as well as the development and implementation
of passive radar systems.

