#### Abstract

Taking into account the difficulty of shape estimation for the extended targets, a novel algorithm is proposed by fitting the B-spline curve. For the single extended target tracking, the multiple frame statistic technique is introduced to construct the pseudomeasurement sets and the control points are selected to form the B-spline curve. Then the shapes of the extended targets are extracted under the Bayes framework. Furthermore, the proposed shape estimation algorithm is modified suitably and combined with the probability hypothesis density (PHD) filter for multiple extended target tracking. Simulations show that the proposed algorithm has a good performance for shape estimate of any extended targets.

#### 1. Introduction

In the traditional low resolution sensor system, each target is tracked as a single point source; that is, its extension is assumed to be neglectable in comparison with sensor resolution. With the increase of the resolution of modern radars and other detection equipment, the echo signal of a target may be distributed in a different range resolution cell; thus, the measurement is no longer equivalent to a point; that is, a single target may generate multiple measurements. Such target is referred to as an extended target in [1–4]. Recently, extended target tracking (ETT) is a hot topic in the field including the short-range applications or maritime surveillance, which has drawn a considerable attention [5–9].

In the conventional extended target tracking, the measurements are modeled as a spatial distribution model in [1], and two examples—a point target with more measurement sources and an object with infinitely thin stick—are used to prove the effectiveness of the approach. Poisson process with a spatially dependent rate parameter is introduced in [2], assuming that each target produces measurements with Poisson distributed random number. It is considered that in this measurement model, the target is sufficiently far away from the sensor, and the measurements resemble a point cluster rather than a geometric structure [7]. Random matrix (RM) is proposed in [8], which has been used to track elliptical target extension [9]. Another method is random hypersurface model (RHM) [3] which is employed for modeling the target extent. However, these methods can only effectively achieve the shape estimate for the target with similar ellipse shape. They cannot effectively estimate the irregular shape of the extended target. In [4], star-convex target extension estimation method is proposed based on RHM under the condition that a measurement source is assumed to be an element of a randomly scaled version of the shape boundary. Moreover, the one-dimension probability density needs to be specified in advance in star-convex shape estimation and it is assumed to be independent of the shape.

To solve the aforementioned problem, a novel shape estimation algorithm based on the B-spline curve fitting is proposed in this paper, and then the proposed shape estimation method is integrated into the framework of extended target probability hypothesis density (ET-PHD) filter for multiple extended target tracking [5, 6]. Simulations show that the proposed algorithm has a good performance for shape estimate of any extended targets.

#### 2. Backgrounds

##### 2.1. Kalman Filter

Assume the state equation and the measurement equation of a single target in two-dimensional plane are given by where and are the state vector and the observation vector at time , respectively. and are the transition matrix and the measurement matrix, respectively. and are the process noise and the observation noise and are uncorrelated Gaussian white noise vectors with covariance matrixes and , respectively.

Suppose that and the covariance are optimal estimations at time in the fusion center; then, the recursive steps of the Kalman Filter (KF) at time are as follows [10].(1)Prediction of state and covariance: (2)Calculating gain: (3)Update of state and covariance:

##### 2.2. The B-Spline Curve

Assume that readers are familiar with the concepts of B-spline curves. A smooth subsection curve can be obtained by fitting the control point set. The B-spline curve of order can be described as [11] where in is a control point and is the B-spline basis function, which is defined over a knot vector . The basis function can be recursively defined as [11, 12] When , we can obtain that

#### 3. Single Extended Target Shape Estimation

In this section, the Bayesian filter framework is introduced for single extended target state and shape estimates. Assume the state equation and the measurement equation are the same as (1) and (2). We define as the center of the extended target and as the control matrix including the shape information. They are not related to each other. Thus we can estimate them with two parallel KFs, which can modify them by recursion update. And the detailed steps of the proposed algorithm are described as follows.

*Step 1. *At time , initialize parameters , , , and , where and denote the state covariance and the shape covariance, respectively.

*Step 2. *When , implement the KF for state estimate:(2.1)prediction of the state and covariance according to (3) and (4), respectively;(2.2)update of the state by the latest measurement set :
where and denotes the number of the measurements at time . denotes the th measurement in . The covariance can be updated by (7).

*Step 3. *Construct a new pseudomeasurement set .(3.1)Let , and then add to the pseudomeasurement set , that is, .(3.2)Set as the largest number of pseudomeasurements. If , then let the new pseudomeasurement set . Otherwise, delete the foregoing measurements from , and then let .

*Step 4. *Update the shape of the extended target by the pseudomeasurement set .

(4.1) Divide the interval into equal angles, and generate a fixed angle set . Then partition the pseudomeasurement set by the angles, and define as the measurement set in the area of the th partition. That is,
where is the width of partition. , , , and denote the position of - and -coordinate at time . denotes a line through the origin of coordinate and along the th partition angle direction. denotes the perpendicular distance between and line . The parameters and of line satisfy . denotes a line through the origin of coordinate, and it is perpendicular to line . The parameters and of line satisfy . Figure 1 shows the th angle partition and the measurements in the area enclosed by points O1, O2, O3, and O4 which constitute the measurement set .

(4.2) Calculate which is the expectation of the perpendicular distance between the elements of and line . Then the radial extension distance can be obtained by
In this paper, we set the as the control point of the extended target shape parameter and define as the control matrix of the shape, which includes the shape information.

(4.3) Shape estimation by implementing the one-dimension (1-D) KF: assume the shape control matrix and the shape covariance have been obtained. Then the recursion steps of shape estimates are as follows.

Prediction:
where and denotes the measurement noise covariance.

Update:
where the shape covariance . denotes the measurement noise.

*Step 5. *Shape estimation according to and the B-spline curve fitting technique.(5.1)Map the control points to the Cartesian coordinates by
(5.2)Produce a closed control point set by adding the element to the end of the , and describe it as . Then the closed cubic B-spline curve can be obtained by

#### 4. Multiple Extended Target Shape Estimation

##### 4.1. Multiple Extended Target PHD Filter

The standard PHD filter for single measurement target tracking has been described in [13–15]; however, it is not suitable for METT. Recently, Mahler has derived a correct equation for extended target PHD (ET-PHD) filter based on the Poisson multitarget measurement model [5]. The Gaussian mixture implement of the multiple extended target PHD filter is presented in [6] and referred to as ET-GM-PHD filter. The detailed filter process of ET-GM-PHD can refer to [6, 16, 17].

##### 4.2. Multiple Extended Target Tracking Algorithm

In this section, we combine the proposed shape estimation algorithm into the framework of ET-GM-PHD filter, which can effectively achieve the multiple extended target tracking with different shape estimation. We refer to this algorithm as Shape-ET-GM-PHD, and its steps are as follows.

*(1) Prediction.* Assume that the state vector and its shape vector are independent, and the survival and detection probabilities are independent of them; that is, and . Assume that denotes the joint posterior PHD function at time and is approximated by the Gaussian mixture distributions. Then the predicted PHD function can be described as
where
, and denotes the mean of the shape standard variance, One has , . denotes the divided number of the shape area according to the degree of angle. and denote the control matrix of shape and its predicted covariance, and denotes the shape process noise covariance. and can also be described as and , where the definition of is the same as that of (13). denotes the predicted variance of the shape.

*(2) Update*. The updated formula of the intensity function can be described as
where denotes an intensity function of the state with shape information and denotes the pseudolikelihood function. denotes the measurement set at time . When ,
otherwise,
where denotes a partition subset of the measurement set , denotes a subset of a partition , and . denotes the measurement likelihood function of one measurement originating from an extended target , is the measurement expectation, and denotes the detection probability of the sensor. Clutter has a Poisson distribution, and its density can be described as , denotes the mean number of clutter measurements, and is the space distribution of the clutter. Consider
where denotes the Kronecker delta function and . denotes the clutter distribution of the measurement space.

In the following subsections, assume that the current estimated PHD can be approximated as a Gaussian mixture distribution. The corrected PHD can be described as where denotes the PHD of the targets without detecting cases, and it can be described as where , , , and .

denotes the PHD of the detected target cases and can be described as where where denotes the penalty coefficient. The shape parameters can be obtained by where denotes the identity matrix with Num order. denotes the measurement noise of a measurement source. denotes the decomposed function which can decompose the matrix as a set with variances. denotes the pseudocontrol point matrix which is used to update the control point matrix , where can be obtained by where denotes the measurement number of th angle direction in cell .

Notice that the D-distance partition method [17] is implemented in the part of measurement partition. Set the maximum distance as the mean size of the shape, and it can be obtained by

Finally, shapes of multiple extended targets are extracted according to and the B-spline curve fitting technique described in Section 2.2.

#### 5. Simulations

Assume that there is an extended target making a uniform motion in a two-dimensional simulation scenario, and the state equation and the measurement equation are the same as (1) and (2), respectively, where denotes the target state, which contains position and velocity information. The state transition matrix and measurement matrix . The process noise covariance with , and the measurement noise covariance with . The initial covariance and the shape noise covariance is . The other parameters , , , and . The way to generate measurements is the same as that of the RM method. The real measurement is assumed as the scattering center, that is, a measurement source. The observation measurements are generated from the scattering center with measurement noises. Notice that the measurement noise is assumed to be small compared to the target extent.

*Example 1 (star-shaped extended target). *Figure 2 shows the shape estimation by the proposed algorithm and the RM method [9], respectively. It is clear that the proposed algorithm has a higher accuracy than that of the RM method. The reason is that the RM method can only estimate the elliptical shape of the target.

Figure 3 shows the measurements of the extended target from the 1st frame to the 20th frame. As can be seen, the measurement noise is big which makes the shape hard to be estimated from the single frame. However, the proposed algorithm can extract the accurate shape by multiple frame statistic technique and the B-spline curve fitting. It is shown that the proposed algorithm has a good capacity of resisting disturbance of the noise.

Figure 4 shows the average shape estimate. It is clear that the proposed algorithm can obtain the shape features effectively.

*Example 2 (Y-shaped extended target). *Figure 5 shows the shape estimation by the proposed algorithm and the RM method, respectively. Figure 6 shows the measurements of the extended target from the 1st to the 20th frames. Figure 7 shows the average shape estimate. As can be seen, the proposed algorithm has a better performance for shape estimate of irregular extended targets.

*Example 3 (multiple extended target tracking). *The scenario of multiple extended target tracking is the same as that of [17]. There are four targets and two of them cross at time , and one target is spawned at time . Assume that the measurement noise covariance and the process noise covariance are and , respectively. Shape noise covariance is , , and . The probabilities of target survival and detection are and , respectively. The birth process is modeled on a Poisson RFS with Gaussian mixture intensity:
where , , , and , , .

The intensity of the spawned targets is , where and .

Figure 8(a) shows the tracking results of the targets in the whole tracking area, and Figures 8(b), 8(c), and 8(d) show the shape estimates of different targets in special tracking area. As can be seen from Figure 8(b), the shapes of Targets 1 and 2 are assumed as ellipses; although they are not accurate, they can be updated approximately to the real shapes of the targets at several time points. The reason is that the shape parameters in the proposed shape estimate method are updated at each time when the latest measurements arrive. Moreover, the Gaussian mixture technique is employed, which can approximately fit the real shape distribution of the extended target. In Figure 8(c), we can see that the shape estimates are not accurate when the targets cross each other, but they can also be updated approximately to the real target shapes. Figure 8(d) shows the shape estimate of a spawned target; it is showed that the proposed algorithm also has a good performance for shape estimate of the spawned target.

**(a) Measurements and tracking results**

**(b) Shape estimates of Targets 1 and 2**

**(c) Shape estimates of cross targets**

**(d) Shape estimates of the spawned target**

Figure 9 shows the number estimates and Figure 10 shows the accuracy statistic by the OSPA distance [18]. Notice that OSPA distance sharply increases at 56th and 66th seconds. The reason is that the targets make a cross with each other at 56th second and a spawned target appears. Generally, we can see that the proposed algorithm has a good performance for multiple extended target tracking with shape estimates.

#### 6. Conclusions

In this paper, a novel shape estimation algorithm is proposed based on the B-spline curve fitting. The multiple frame statistic technique is introduced to construct the pseudomeasurement sets. The selected control points are used to form the B-spline curve, and then the curve yields the shapes of the extended targets. Moreover, the proposed shape estimation algorithm is modified suitably and combined with the probability hypothesis density (PHD) filter for multiple extended target tracking. Simulations show that the proposed algorithm has a good performance for shape estimates of extended targets.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This paper is supported by the National Natural Science Foundation of China (nos. 61305017 and 61304264) and the Natural Science Foundation of Jiangsu Province (no. BK20130154).