The Scientific World Journal

Volume 2013 (2013), Article ID 489121, 7 pages

http://dx.doi.org/10.1155/2013/489121

## Sequential Bearings-Only-Tracking Initiation with Particle Filtering Method

^{1}School of Computer Science and Technology, Nanjing University of Posts and Telecommunications, Nanjing 210023, China^{2}Institute of Acoustics, Chinese Academy of Sciences, Beijing 100190, China

Received 14 August 2013; Accepted 7 October 2013

Academic Editors: N. Bouguila and N. M. Kwok

Copyright © 2013 Bin Liu and Chengpeng Hao. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The tracking initiation problem is examined in the context of autonomous bearings-only-tracking (BOT) of a single appearing/disappearing target in the presence of clutter measurements. In general, this problem suffers from a combinatorial explosion in the number of potential tracks resulted from the uncertainty in the linkage between the target and the measurement (a.k.a the data association problem). In addition, the nonlinear measurements lead to a non-Gaussian posterior probability density function (pdf) in the optimal Bayesian sequential estimation framework. The consequence of this nonlinear/non-Gaussian context is the absence of a closed-form solution. This paper models the linkage uncertainty and the nonlinear/non-Gaussian estimation problem jointly with solid Bayesian formalism. A particle filtering (PF) algorithm is derived for estimating the model’s parameters in a sequential manner. Numerical results show that the proposed solution provides a significant benefit over the most commonly used methods, IPDA and IMMPDA. The posterior Cramér-Rao bounds are also involved for performance evaluation.

#### 1. Introduction

The problem of angle/bearings-only tracking (BOT) has been studied over many years due to its tremendous importance in a variety of practical communication and signal processing applications, such as localization in wireless sensor networks [1], submarine tracking using passive sonar [2], aircraft surveillance using radars [3], and sonar-based robotic navigation [4, 5]. This paper is concerned with radar and sonar tracking. Specifically, the focus of this paper lies in a fundamental problem in many sonar and radar tracking tasks, called track initiation. This problem consists in linking sets of point observations from different time steps to fit a desired model without any *a priori* track estimates. In general, track initiation approaches can be categorized into two schemes: sequential and batch schemes [6–8]. The sequential scheme processes a sequence of measurements received during consecutive radar/sonar scans sequentially one at a time. The sequential methods are widely used in radar and sonar tracking problems. For the batch technique approach, past scans of measurements are treated simultaneously to determine feasible target trajectories (see Figure 2). The batch techniques suffer from a heavy computational load and a slow process, while they find applications in image processing and tracking in heavy clutter background [6]. This paper only concerns sequential techniques.

In the context of radar and sonar applications, the target signal, if present, is often corrupted by additive noise and comes with other measurements unrelated to the target, which originate from thermal noise and various forms of clutter such as terrain and clouds [9]. So the issue of measurement origin uncertainty has to be addressed before tracking. This is called the data association problem and has been addressed in the context of multitarget tracking (MTT) in clutter [10, 11].

The tracking initiation problem considered here has its own challenges not only due to the measurement origin uncertainty, but also the uncertainty in the target presence. In fact, for tracking initiation, it is prerequisite to determine whether the target is present before dealing with data association and state filtering. Moreover, in BOT the measurement function is highly nonlinear, which results in a non-Gaussian posterior probability density function (pdf) in the optimal estimate of the target. The consequence of this nonlinear/non-Gaussian context is the absence of a closed-form solution [12–14]. As BOT initiation involves such interactive problems, that is, the measurement origin uncertainty, the target presence uncertainty, and the nonlinear/non-Gaussian posteriors in optimal estimation, it is therefore much more difficult than any individual-related tasks such as state filtering, target detection, and data association. It is noted that, besides the simulated annealing [15, 16] and the Ants algorithm [17], there are few reports on such BOT initiation problem, not to mention sequential BOT initiation.

This paper proposes an algorithm for sequential BOT initiation. The measurement origin uncertainty, target presence uncertainty, and the nonlinear non-Gaussian factors are handled jointly within a Bayesian sequential estimation framework. Based on such Bayesian formalism, a sequential Monte Carlo (SMC), a.k.a particle filtering (PF), algorithm is derived. Performance of the proposed approach is evaluated via numerical simulations and related methods are involved for performance comparison.

The remainder of this paper is organized as follows. Section 2 formulates the sequential BOT initiation problem. Section 3 presents the Bayesian sequential estimation framework for the problem under consideration. Section 4 derives the PF algorithm under the framework developed above. In Section 5, simulations are performed to evaluate the performance of the proposed algorithm. Finally, Section 6 concludes the paper.

#### 2. Problem Formulation

##### 2.1. Dynamic Model of the Target

Consider a discrete-time dynamic model
which is used to describe the target movement in a two-dimensional (2D) plane. denotes the discrete time index, the process noise, and the target state vector. The target state consists of 2D position and velocity elements, which is defined as below:
where and denote the 2D position and velocity, respectively. The superscript denotes matrix transposition. If the movement pattern of the target is *a priori *known, then will have a specific form. For example, if the target moves with a near constant velocity, then can be modeled to be
where , , and denotes the sampling period of the measurements. For more details on dynamic models for tracking problems, see [18].

##### 2.2. Markov Chain Model for Uncertainty of Target Existence

In tracking initiation problem, the target presence is not a definitive issue either before the track under consideration is confirmed or after it is terminated. To take into account the uncertainty of the target presence, we introduce binary valued variable and let represent the target being present in the surveillance region at time , and denotes the opposite. Evolution of during the surveillance time is modeled by a Markov chain, which is characterized by the following transition probabilities:
and the initial probability of target existence, that is, . and , have specific physical meanings. is actually the probability of a new target’s birth or appearing at a time step and is the probability of the target’s death or disappeance. The transitional probability matrix is given by
where and are equivalent to probabilities of target staying alive and remaining absent, respectively. , , and are assumed *a priori *known.

##### 2.3. Measurement Model

Let denote the measurements received up to and including the time step . Here, represents the measurement received at time . It is composed of only clutters (if the target is absent currently) or clutters along with the target originated measurement element (if the target is present at this time step). Specifically, where denotes the target originated measurement, measurement elements due to clutters, and the number of clutters. Note that elements in are not ordinal. For BOT, the relationship between and is specified to be where is the measurement noise with known pdf and is independent of in (1) and the sensor’s position.

The measurement elements due to clutters are assumed to be uniformly distributed over the surveillance region, so the pdf , where denotes volume of the surveillance region. The number of clutters per scan, , is Poisson distributed with density parameter ; that is, ~.

##### 2.4. The BOT Initiation Problem

The information on target presence along with the target’s dynamic state (if target presence has been confirmed) is of interest here, while the only probe to get the desired information is the received measurements.

In Figure 1, we give an example showing the measurements generated in a BOT initiation case, where the target exists from the 8th second to the 32nd second. In the upper plot, the target-originated measurements are linked by a series of line segments, while the other points denote measurements generated due to clutters. The bottom plot represents the true scenario encountered in practice where there is no indication on the source of the measurements.

The task of track initiation just consists in taking sets of point observations from different time steps and linking together those observations that fit a desired model without any previous track estimates. As shown in Figure 1, we may draw many different lines to link many different sets of measurements, while there is only one true answer underlying these measurements. In general, such track initiation problem suffers from a combinatorial explosion in the number of potential tracks that must be evaluated.

As indicated by Figure 1, it is reasonable to input the measurements to the tracker algorithm in a batch mode, as latter observations can be used to benefit making decisions at previous time steps. That is why so many batch mode techniques are developed for track initiation, for example, the Hough-transform-based approach [6, 8, 19–21], the multiple-kd tree algorithm [22], and so on.

However, the problem of interest here is different. The focus is how to do BOT initiation in a sequential manner. We formulate this problem as follows. Given , once is observed, how to calculate immediately, . Actually, covers all the information about that we are able to get based on . It is noted that, provided that sufficient statistics of can be obtained, there is no need to reload and process for calculating upon the arrival of .

#### 3. Sequential Bayesian Estimation Framework for BOT Initiation

Here, the distribution of interest is the posterior . We present a Bayesian estimation framework for computing this posterior distribution sequentially in a recursive manner. The recursion starts from an initial distribution , where is empty, since at the very beginning, there is no measurement element received. The successive posteriors are computed as follows: Now let us examine individual parts that constitute the above integrals. First, consider the state transition part: Note that is just determined by the birth/death Markov model (6). If , the target is absent, so is undefined; otherwise, the pdf of conditional on and is given by where denotes the initial pdf of the target on its appearance.

Now, we consider the calculation of the measurement likelihood. To begin with, introduce variable to denote the number of measurement elements in . If the target is absent, that is, , then just equals the number of measurements due to clutters, that is, ; otherwise, . Now, we introduce the association variable . We use , to denote the event that the th element in is target-originated, and the others are generated due to clutters. Additionally, denotes that no elements in are target-originated. Now, the measurement likelihood is specified to be where

#### 4. Sequential Monte Carlo Implementation

In this section, SMC techniques are developed to approximately implement the Bayesian estimation framework presented in Section 3. The idea is to use the Monte Carlo sampling methodology to approximate the involved distributions by a set of weighted random samples. It is noted that convergence of SMC methods in dealing with nonlinear non-Gaussian Bayesian tracking problem has been proved [23–25].

To begin with, denote as a random measure that approximates the posterior at ; namely, . Here, is the sample size; is the importance weight of the th sample . In the following, an attempt is made to derive an evolution law for the current random set of particles to get a particle approximation to .

First, evolve the existence variable one time step further for each particle based on the birth/death Markov model (6). For example, if , then the probabilities of and are and , respectively, according to (4) and (6).

Next, we show how to generate state vector value for each particle according to (10) and (11). Let us focus on particles that are associated with positive values. They are active particles at time step . These particles can be categorized into two classes:(i)*newborn particles*: this set of particles is characterized by the transition from to ;(ii)*staying alive particles*: this is a set of particles that continues to stay active with and .According to (11), we draw state vector sample values for the *newborn particles* from a given initial pdf and draw state vector sample values for the *staying alive particles* from the dynamic transition prior density, which are determined by the target dynamic model (1). For the rest of particles that are associated with , just keep their state vector values undefined, since they are totally useless in the algorithm.

Now we see, for each particle, say the th, its associated state vector value and existence variable value are obtained, while the importance weight needs to be updated when the new measurement arrives at time step . We update the importance weights according to the solid Bayesian formalism given by (9). Given and , the (unnormalized) importance weight is just equivalent to the likelihood of . It is calculated according to (12). A normalization step is processed as follows:
to make sure that the summation of the importance weights equals 1. Then, the posterior probability of target existence, that is, , is computed as below
which satisfies . The target is declared to be present if is bigger than a given threshold, 0.6 is used in our simulation test. If the target is determined to be present, its state is estimated as follows:
To avoid the problem of degeneracy of SMC algorithm, that is, avoiding the situation that all but one of the importance weights are close to zero, a resampling procedure is performed if the effective number of particles is less than a given threshold . An estimate of the effective number of particles is computed as follows [26]:
The resampling procedure can be summarized as follows: draw particles from the current particle set with probabilities proportional to their weights and then replace the current particle set with this new one. For more discussions on resampling methods in context of PF, refer to [23, 27, 28]. The above SMC algorithm is initialized at by drawing samples , in accordance with the initial target existence probability . The initialization of the active particles’ state vector is the same as for the *newborn particles* described above.

#### 5. Simulation Study

In this section, the performance of the proposed PF algorithm for sequential BOT initiation is evaluated via simulations. The estimation accuracy of the proposed algorithm is compared to the Posterior Cramér-Rao lower bounds (CRLB) [29]. Two generally accepted approaches for sequential track initiation, namely, integrated probabilistic data association (IPDA) [9] and integrated multiple model probabilistic data association (IMMPDA) [30, 31], are also involved for performance comparison.

##### 5.1. Simulation Setting

The scenario to be investigated involves BOT of a target from an observer, as shown in Figure 1. The observer travels at a fixed speed of 10 m/s and executes 2 maneuvers during the period under surveillance. The scanning process lasts 40 seconds in total. During the first 8 seconds, the target does not emit any energy outward, so there is no target-originated elements in the received measurements. During the period from the 8th to the 32nd second, the target moves and emits energies outward, so the received measurements include the target-originated elements and the tracker algorithm is expected to be able to detect the target. After the 32nd seconds, generation of target-originated elements is terminated in the received measurements. This setting is used to test whether a tracker algorithm under examination can terminate the track in time. The target motion is simulated according to (3) subjected to an amount of process noise with . The target’s 2D position and velocity are initialized to be and , respectively. The other parameters used in this simulation are listed in Table 1. An example show of the received measurements versus discrete time index is presented in Figure 1.

##### 5.2. Performance Comparison on Target Detection

Detection performance of the involved algorithm is investigated via Monte Carlo simulations with 100 independent runs of each algorithm. The outputted probabilities of target existence at each time step are averaged and then plotted in Figure 3. As is shown, the proposed SMC algorithm beats both IPDA and IMMPDA remarkably in the rate of convergence.

Quantitative comparison is also conducted via hypothesis testing. Hypotheses and are defined to represent the events “target being absent” and “target being present,” respectively. At each time step , four possible cases may happen for each algorithm:(1) is true and the algorithm chooses ,(2) is true and the algorithm chooses ,(3) is true and the algorithm chooses ,(4) is true and the algorithm chooses .We select probability of false alarm, denoted by , and probability of detection, denoted by , as the performance metrics. Let count the number of times of case appearance in the simulation test, . Then and . So a good algorithm is expected to yield bigger and smaller . and are calculated over 100 independent runs of each algorithm and the final result is listed in Table 2. It is demonstrated again that the proposed SMC algorithm is superior to IPDA and IMMPDA in detection performance.

##### 5.3. Examination on Estimation Accuracy

The posterior Cramer-Rao lower bounds (CRLB) [29] are used to indicate the best estimation accuracy an algorithm can achieve. The performance metric, namely, the root-mean square (RMS) position error is defined as follows: where denotes the target’s 2D position at time step in the simulations, the estimation of outputted in the th test of the algorithm under consideration, and the number of independent tests under consideration. Define , as the 2D position elements of the inverse information matrix for the problem at hand; the corresponding CRLB for the metric given by (18) is calculated as below: Here, we only consider the scanning time interval between the 8th and the 32nd time steps, when the target is truly present and the proposed algorithm successfully detects it. Figure 4 shows the RMS position error of the proposed algorithm in comparison with the calculated CRLB.

#### 6. Conclusions

This paper addresses the problem of sequential BOT initiation in the context of sonar and radar applications. A sequential Bayesian estimation framework is developed for this problem. This theoretical framework addresses measurement origin uncertainty and target existence uncertainty jointly via solid Bayesian formalism. An SMC approximate algorithm is derived under this framework and its performance is evaluated via simulations. It is shown that the proposed algorithm provides a remarkable performance improvement in target detection, compared with the commonly used PDA based methods. The proposed algorithm also gives accurate estimation of the target’s state, as indicated by a comparison with the posterior CRLB.

#### Acknowledgments

The authors appreciate the insightful comments and suggestions provided by the anonymous reviewers. This work was supported by the National Natural Science Foundation (NSF) of China (Grants nos. 61302158, 61172166, 61100135, 61170065, and 61373017), Provincial Science and Technology Plan (NSF) of Jiangsu province (Grant no. BK20130869), and Natural Science Research Project for Colleges and Universities in Jiangsu province (Grant no. 13KJB520019).

#### References

- G. Mao, B. Fidan, and B. D. O. Anderson, “Wireless sensor network localization techniques,”
*Computer Networks*, vol. 51, no. 10, pp. 2529–2553, 2007. View at Publisher · View at Google Scholar · View at Scopus - L. Mattos and E. Grant, “Passive sonar applications: target tracking and navigation of an autonomous robot,” in
*Proceedings of the IEEE International Conference on Robotics and Automation*, vol. 5, pp. 4265–4270, May 2004. View at Scopus - A. Farina, “Target tracking with bearings—only measurements,”
*Signal Processing*, vol. 78, no. 1, pp. 61–78, 1999. View at Publisher · View at Google Scholar · View at Scopus - N. M. Kwok, G. Dissanayake, and Q. P. Ha, “Bearing-only SLAM using a SPRT based gaussian sum filter,” in
*Proceedings of the IEEE International Conference on Robotics and Automation*, pp. 1109–1114, April 2005. View at Publisher · View at Google Scholar · View at Scopus - P. Jensfelt, D. Kragic, J. Folkesson, and M. Björkman, “A framework for vision based bearing only 3D SLAM,” in
*Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '06)*, pp. 1944–1950, May 2006. View at Publisher · View at Google Scholar · View at Scopus - H. Leung, Z. Hu, and M. Blanchette, “Evaluation of multiple target track initiation techniques in real radar tracking environments,”
*IEE Proceedings on Radar, Sonar and Navigation*, vol. 143, no. 4, pp. 246–254, 1996. View at Google Scholar - F. Kural, F. Arikan, O. Arikan, and M. Efe, “Performance evaluation of the sequential track initiation schemes with 3D position and doppler velocity measurements,”
*Progress In Electromagnetics Research B*, no. 18, pp. 121–148, 2009. View at Publisher · View at Google Scholar · View at Scopus - Z. Hu, H. Leung, and M. Blanchette, “Statistical performance analysis of track initiation techniques,”
*IEEE Transactions on Signal Processing*, vol. 45, no. 2, pp. 445–456, 1997. View at Google Scholar · View at Scopus - D. Musicki, R. Evans, and S. Stankovic, “Integrated probabilistic data association,”
*IEEE Transactions on Automatic Control*, vol. 39, no. 6, pp. 1237–1241, 1994. View at Publisher · View at Google Scholar · View at Scopus - B. Liu, C. Ji, Y. Zhang, C. Hao, and K.-K. Wong, “Multi-target tracking in clutter with sequential Monte Carlo methods,”
*IET Radar, Sonar and Navigation*, vol. 4, no. 5, pp. 662–672, 2010. View at Publisher · View at Google Scholar · View at Scopus - J. Vermaak, S. J. Godsill, and P. Pérez, “Monte Carlo filtering for multi-target tracking and data association,”
*IEEE Transactions on Aerospace and Electronic Systems*, vol. 41, no. 1, pp. 309–332, 2005. View at Publisher · View at Google Scholar · View at Scopus - B. Liu, X.-C. Ma, C.-P. Hao, C.-H. Hou, and M. Li, “A bearings-only-tracking framework based on the EKF and UKF combined algorithm,” in
*Proceedings of the International Symposium on Intelligent Signal Processing and Communications Systems (ISPACS '07)*, pp. 184–187, December 2007. View at Publisher · View at Google Scholar · View at Scopus - N. J. Gordon, D. J. Salmond, and A. F. M. Smith, “Novel approach to nonlinear/non-gaussian Bayesian state estimation,”
*IEE Proceedings F*, vol. 140, no. 2, pp. 107–113, 1993. View at Google Scholar · View at Scopus - J. H. Kotecha and P. M. Djurić, “Gaussian particle filtering,”
*IEEE Transactions on Signal Processing*, vol. 51, no. 10, pp. 2592–2601, 2003. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - P.-Y. Ting and R. A. Iltis, “Multitarget motion analysis in a DSN,”
*IEEE Transactions on Systems, Man and Cybernetics*, vol. 21, no. 5, pp. 1125–1139, 1991. View at Publisher · View at Google Scholar · View at Scopus - B. Xu, Q. Chen, and Z. Wang, “Ants for track initiation of bearings-only tracking,”
*Simulation Modelling Practice and Theory*, vol. 16, no. 6, pp. 626–638, 2008. View at Publisher · View at Google Scholar · View at Scopus - X. R. Li and V. P. Jilkov, “Survey of Maneuvering Target Tracking. Part I: dynamic Models,”
*IEEE Transactions on Aerospace and Electronic Systems*, vol. 39, no. 4, pp. 1333–1364, 2003. View at Publisher · View at Google Scholar · View at Scopus - K. Alexiev, “Implementation of hough transform as track detector,” in
*Proceedings of the 3rd International Conference on Multisource Multisensor Information Fusion (FUSION '00)*, vol. 2, pp. THC4–TH11, IEEE, 2000. - L. Zhou, W. Xie, and L. Li, “Mean shift track initiation algorithm based on hough transform,” in
*Proceedings of the IEEE 10th International Conference on Signal Processing (ICSP '10)*, pp. 1263–1266, October 2010. View at Publisher · View at Google Scholar · View at Scopus - S. W. Yankowich and M. Farooq, “Hough transform based multisensor, multitarget, track initiation technique,”
*Optical Engineering*, vol. 37, no. 7, pp. 2064–2077, 1998. View at Google Scholar · View at Scopus - J. Kubica, A. Moore, A. Connolly, and R. Jedicke, “Fast and robust track initiation using multiple trees,” Tech. Rep. CMU-RI-TR-04-62, Robotics Institute, 2004. View at Google Scholar
- M. S. Arulampalam, S. Maskell, N. Gordon, and T. Clapp, “A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking,”
*IEEE Transactions on Signal Processing*, vol. 50, no. 2, pp. 174–188, 2002. View at Publisher · View at Google Scholar · View at Scopus - D. Crisan and A. Doucet, “A survey of convergence results on particle filtering methods for practitioners,”
*IEEE Transactions on Signal Processing*, vol. 50, no. 3, pp. 736–746, 2002. View at Publisher · View at Google Scholar · View at Scopus - X.-L. Hu, T. B. Schön, and L. Ljung, “A basic convergence result for particle filtering,”
*IEEE Transactions on Signal Processing*, vol. 56, no. 4, pp. 1337–1348, 2008. View at Publisher · View at Google Scholar · View at Scopus - J. S. Liu and R. Chen, “Sequential Monte Carlo methods for dynamic systems,”
*Journal of the American Statistical Association*, vol. 93, no. 443, pp. 1032–1044, 1998. View at Google Scholar · View at Scopus - G. Kitagawa, “Monte Carlo filter and smoother for non-Gaussian nonlinear state space models,”
*Journal of Computational and Graphical Statistics*, vol. 5, no. 1, pp. 1–25, 1996. View at Google Scholar · View at Scopus - A. Doucet, S. Godsill, and C. Andrieu, “On sequential Monte Carlo sampling methods for Bayesian filtering,”
*Statistics and Computing*, vol. 10, no. 3, pp. 197–208, 2000. View at Google Scholar · View at Scopus - P. Tichavsky, C. H. Muravchik, and A. Nehorai, “Posterior cramér-rao bounds for discrete-time nonlinear filtering,”
*IEEE Transactions on Signal Processing*, vol. 46, no. 5, pp. 1386–1396, 1998. View at Publisher · View at Google Scholar · View at Scopus - Y. Bar-Shalom, K. C. Chang, and H. A. P. Blom, “Automatic track formation in clutter with a recursive algorithm,” in
*Proceedings of the 28th IEEE Conference on Decision and Control*, pp. 1402–1408, IEEE, December 1989. View at Scopus - A. Houles and Y. Bar-Shalom, “Multisensor tracking of a maneuvering target in clutter,”
*IEEE Transactions on Aerospace and Electronic Systems*, vol. v, pp. 176–189, 1992. View at Google Scholar · View at Scopus