A Novel Track Initiation Method Based on Prior Motion Information and Hough Transform

. This paper deals with the problem of track initiation in dense clutters. A novel Hough transform method based on prior motion information is proposed. Firstly, with time sequence of the measurements and targets’ kinematic information considered, the single scan data accumulation effect is avoided, and the initiated tracks are more reliable. Then, measurements in observation space only vote to part of the lines that pass through them, which can effectively reduce the computational cost, and improve the initiation speed. Simulation results indicate that the proposed method is superior to the traditional algorithms in terms of detection probability, computational cost and robustness.


Introduction
Track initiation is the primary problem of multi-target tracking, especially for establishing the records for new targets and terminating the records of inexistent targets [1,2].Because of the long distance between the radar and targets, it is hard for the targets to be accurately measured by sensors, and there is no statistical rule to judge whether the measure is from target or clutter, track initiation itself is a very difficult issue.It is hoped that target tracks are initiated as early as possible, because the operators can get an early situation picture for the surveillance area.On the other hand, it is hoped that false track initiation rate is as low as possible.False tracks cause waste of radar resources, because radar beams for track maintenance are assigned to the area where there is no target [3].In addition, to detect modern dim targets, radar needs to lower the threshold.Therefore, the false alarm rate will be very high, and how to deal with the track initiation of targets moving in such circumstance has been a challenging task now.

Review of current techniques
During the past several decades, a lot of research has been contributed to track initiation, and there are two main series, one is sequential data processing technique, such as the heuristic rule method [4], logic-based method [1]; another is batch data processing technique, such as multiple hypothesis tracking method [5], integrated probabilistic data association algorithm [6], the Hough transform technique [7], and the modified Hough transform technique [8].The former method has low computational cost, and is easy to apply.However, it is just suitable for initiation in sparse clutter.Though the latter method has a better performance than the former method, and is effective in reducing false alarm probability, it has much more computational cost than the former method, in addition, it is difficult to apply.
As is known, Hough Transform has low sensitivity to local fault and good ability in suppressing noise and clutters.Because of this advantage, it has been widely applied to track initiation in dense clutter environment.As a batch data processing method, Hough transform has the disadvantage of huge computation and memory.A lot of research has been done to reduce computation and memory by researchers, and many advanced methods have been put forward, such as modified Hough transform, random Hough transform, etc.A novel method for fast track initiation in dense clutters by divide and Hough transform method is proposed in this paper and this method is very effective in decreasing computation and memory.The measurement's time sequence information and targets' kinematics information also have been effectively utilized.

Our approach
For Hough Transform, any point in data space is corresponding to the only one curve in parameter space, and points of the same line in data space are corresponding to different curves in parameter space, and those curves intersect at one point in parameters.If , , denotes as an echo in the th , ( ) , ( ) scan respectively, then and are connected to form a line C 1 and C 2 will be very close, conversely, they will be far apart.According to this principle, a track can be detected through the intensive degree of point distribution in parameter space.
Fig. 1.Three collinear points and the corresponding curves In this new method, any two observed data from different cycles is connected, which avoids loss of echo data.If only observed data from two adjacent cycles can be connected, when the target echo appears flashing, for example, a target echo is lost every other cycle, then echo data from two adjacent cycles can not be connected, which will cause loss of target information.However, not all observed data from two different cycles are from the same target, only when distance between the two points satisfies some condition, they can be from the same target.If     , , , z x y z x y is the measurement in time 12 , tt respectively, then distance between the two points is where d is the distance.Then the average speed is where v denotes as the average speed.If 1 z and 2 z can be from the same target.In formula (3), min v and max v denote as the minimum and maximum speed of the target respectively.Through this process, the number of combinations is greatly reduced.If 1 z and 2 z satisfy the above condition, then the corresponding curve through Hough Transform in parameter space is Bring formula (6) back into equations (4),  is obtained.x y x y d During the initiation period, targets are far away from sensors, so they are often thought to move in a straight line.According to the above discussion, any two measurements from a target must be mapped to the same   ,  in the parameter space.However, because of target's motion and measurement error, it is impossible that all   ,  is the same, but they will be concentrated in a small area.These small areas can be detected through the density of point distribution in the parameter space, and then the corresponding tracks are detected.
Suppose there are N intersections in parameter space, to detect target trajectory from observed data, the   parameter space is equidistantly quantized into NN   cells with distance interval   and angle interval   .Define parameter space array according to   and   , if an intersection votes to accumulated cell   Search for the peak values in accumulated cells, and ) output the corresponding observed data as the initial tracks of targets.

Simulation and results
In this section, suppose the surveillance area is 100000m long and 100000m wide.And assume 5 targets move in a straight line with constant velocity, and these targets are tracked with 2D radar.The initial position of 5 targets are (55000m, 55000m), (45000m, 45000m), (35000m, 35000m), (45000m, 25000m), (55000m, 15000m), and 5 targets move in the same speed, m s, m s.And assume radar's PRI is s, azimuth observation square root error is 0.
where N  is the number of Monte Carlo simulation, i f is number of false tracks in th i Monte Carlo simulation, and i n is number of all tracks in th i Monte Carlo simulation.In the same way, the detection probability of target t is given as where it l indicates whether track of target t is correctly initiated in th i Monte Carlo simulation, and As is seen from Table 1, with the increase of clutters, when 100

 
, there are many false tracks in the former 4 algorithms.Especially in HT method, a great number of false tracks make it impossible to detect the exact tracks for targets.All these false tracks will cause a lot of trouble for subsequent target tracking and waste of radar resources.However, the method proposed in this paper is very robust to clutter interference, which can correctly initiate targets' tracks and suppress false tracks generated by clutters.
For all these methods, the proposed method always keeps the highest probability of detection in both simulation scenarios.Even when 100

 
, it still keeps 63.2%, while the highest detection probability of the remaining methods is 45.3%.These simulation results indicate that the new method has a superior performance for track initiation, and is very robust to clutter interference.To further validate the effectiveness of the proposed algorithm, describes the computational cost for methods based on Hough transform.From the table we know that the proposed method is less computational than the HT method and the modified HT method, and is more suitable for fast track initiation in dense clutters.

Conclusions
In this paper, a novel method for Fast track initiation in dense clutters by divide and Hough transform method is proposed, which connects any two points from different cycles and calculates the exact intersection of the two points in parameter space.With the measurement's time sequence information and targets' kinematics information considered, this new algorithm avoids a single scan data accumulation effect, and can improve the reliability of formed tracks.In this new algorithm, measurements in observation space only vote to part lines instead of all in standard Hough transform that pass through them, so it can effectively reduce computation and memory, and improve the speed of track initiation.In addition, the simulation results also validate the superior performance of the new algorithm and indicate that it can effectively eliminate the false tracks, and is more robust to clutter inference.

1
ii AA  , which is corresponding to a point C 1 in parameter space with Hough Transform.In the same way, and are connected to form a line 12 ii AA  , which is corresponding to a point C 2 . in parameter space with Hough Transform.If , , are collinear, C 1 and C 2 must coincide as shown in Fig. 1.If , , are approximately collinear, then a fair comparison, we make 100 N  independent Monte Carlo runs, and the number of false measurements in the surveillance area satisfies the passion distribution with density  .Define the false track initiation probability as

Table 1
Performance comparison of 5 track initiation methods

Table 2
Comparison of computational cost for HT based methods