# Patent

## Title
Target tracking

## Abstract

An apparatus and a method are disclosed for tracking a plurality of targets (e.g. land-based vehicles) The method can include: for a first time-step, estimating a state of each target; at a second time-step, measuring values for a state of a target; for the second time-step, estimating a state of each target using the estimated target states for the first time-step; updating the estimated target states for the second time-step using the measured values; and performing an identity management process to estimate a probability that a particular state measurement corresponds to a particular target by providing a mixing matrix, wherein an element of the mixing matrix is based on an overlap between the updated estimated target states of targets in an underlying state space of the estimated target states.

## Background
FIELD OF THE INVENTION
The present invention relates to methods and apparatus for tracking a plurality of targets.
BACKGROUND
Target tracking typically comprises inferring an unobserved state vector of a target (for example a vector including a target's position and velocity) from a given set of sensor observations and prior assumptions about target behaviour and sensor properties.
Conventional multi-target tracking algorithms tend to encounter problems when implemented in certain situations. In particular, problems tend to be encountered in situations in which the motion of targets is constrained and the paths of targets overlap or are proximate to one another (e.g. road vehicle targets may be constrained to travel on roads, and a portion of a track followed by a target may be shared by multiple targets).
SUMMARY OF THE INVENTION
In a first aspect, the present invention provides a method of tracking a plurality of targets, the method comprising, for a first time-step, estimating a state of each respective target, at a second time-step, measuring values for a state of a target, for the second time-step, estimating a state of each respective target using the estimated target states for the first time-step, and updating the estimated target states for the second time-step using the measured values; and performing an identity management process to estimate the probability that a particular state measurement corresponds to a particular target; wherein the identity management process comprises providing a mixing matrix; and an element of the mixing matrix is based on an overlap between the updated estimated target states of targets in an underlying state space of the estimated target states.
The step of providing the mixing matrix may comprise determining values for the elements {circumflex over (m)}qr, where:

{circumflex over (m)}qr=∫√{square root over (pq(x)pr(x))}{square root over (pq(x)pr(x))}dx; 

pi(x)=p(xi=x|Y) for i=q, r; 

    
    
        xi is a state of an ith target at the second time-step; and
        Y is a set of all state values and identity parameter values measured at time steps up to and including the second time-step;

and forming the elements into matrix {circumflex over (M)}=[{circumflex over (m)}qr].

    
    


Determining the mixing matrix may further comprise normalising the matrix {circumflex over (M)}.
The step of updating the estimated target states for the second time-step using the measured values may be performed using one or more of: a multi-hypothesis tracking process, a maximum likelihood data association process, and a joint probabilistic data association process.
The step of updating the estimated target states for the second time-step using the measured values may be implemented using a particle filter and/or a Kalman filter.
Motion of the targets may be constrained, to some extent.
The targets may be land-based vehicles.
A path of one of the targets may be proximate to a path of at least one different target.
The path of one of the targets may overlap to some extent the path of at least one different target.
Estimating a state of each respective of the targets may be performed using a particle filter.
A state of a target may comprise a value for a position of the target and/or a value for a velocity of the target.
In a further aspect, the present invention provides apparatus for tracking a plurality of targets, the apparatus comprising one or more processors arranged, for a first time-step, to estimate a state of each respective target; and one or more sensors arranged, at a second time-step, to measure values for a state of a target; wherein the one or more processors are arranged, for the second time-step, to estimate a state of each respective target using the estimated target states for the first time-step, update the estimated target states for the second time-step using the measured values, and perform an identity management process to estimate the probability that a particular state measurement corresponds to a particular target, wherein the identity management process comprises providing a mixing matrix, and an element of the mixing matrix is based on an overlap between the updated estimated target states of targets in an underlying state space of the estimated target states.
The elements {circumflex over (m)}qr of the mixing matrix may have the following values:

{circumflex over (m)}qr=∫√{square root over (pq(x)pr(x))}{square root over (pq(x)pr(x))}dx; 

where:

pi(x)=p(xi=x|Y) for i=q, r; 

    
    
        xi is a state of an ith target at the second time-step; and
        Y is a set of all state values and identity parameter values measured at time steps up to and including the second time-step;

and the elements are formed into matrix {circumflex over (M)}=[{circumflex over (m)}qr].

    
    


In a further aspect, the present invention provides a program or plurality of programs arranged such that when executed by a computer system or one or more processors it/they cause the computer system or the one or more processors to operate in accordance with the method of any of the above aspects.
In a further aspect, the present invention provides a machine readable storage medium storing a program or at least one of the plurality of programs according to the above aspect.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 is a schematic illustration of an example of a multi-target tracking device;
FIG. 2 is a schematic illustration of an example multi-target tracking scenario;
FIG. 3 is a process flowchart showing certain steps of an embodiment of a tracking and identification process; and
FIG. 4 is a process flowchart showing certain steps of step s6 of the process flowchart of FIG. 3.



DETAILED DESCRIPTION
FIG. 1 is a schematic illustration (not to scale) of an example multi-target tracking device, hereinafter referred to as the “target tracker 2”, that may be used to implement an embodiment of a process for tracking multiple targets and managing the identities of these targets. This process is hereinafter referred to as the “tracking and identification process”.
In this example, the target tracker 2 comprises two sensors, namely a sensor capable of measuring the position and velocity of a target (hereinafter referred to as the “position sensor 4”), and a sensor capable of sensing the identity of a target (hereinafter referred to as the “identity sensor 6”).
In other examples, the target tracker 2 comprises a different number of sensors or modules such that the same functionality is provided. For example, a single sensor may be used to produces position/velocity based and identity based information for a target. Also, in other examples identity information may be provided by processing position/velocity based information, e.g. the detection measurements may be processed by an automatic target recognition algorithm to generate target identity information (e.g. from target shape, colour, acoustic signature, etc).
The target tracker 2 further comprises a processor 8 coupled to the position sensor 4 and the identity sensor 6.
In this example, the position sensor 4 produces bearing and/or range information from the target tracker 2 to a target. The position sensor may, for example, be one or more acoustic arrays and/or electro-optical (EO) devices. In other examples, this sensor may measure other parameters related to a state of target, for example acceleration. In this embodiment, the state of a target comprises values for position and velocity parameters, i.e. variable parameters (parameters that have variable values for a given target).
The identity sensor 6 measures the acoustic signature of a target. These measurements enable the targets to be partially distinguished from each other. However, the measurements may be prone to errors. In this embodiment, the identity sensor measures a parameter (in this case acoustic signature), the values of which (ignoring errors or measurement variations) remain substantially constant for a given target (e.g. a target's colour, or shape, or acoustic signature, and so on). Such a parameter, a set of measured values of which remain substantially constant for a given target, is hereinafter referred to as an “identity parameter”. The substantially constant identity parameter values tend to allows targets to be partially distinguished. In other embodiments, the identity sensor 6 measures values of a different appropriate identity parameter for the targets.
The processor 8 implements a conventional particle filter (not shown) that utilises the measurements taken by the sensors to track a target. This is done using a Joint Probabilistic Data Association (process) as described in more detail later below with reference to FIGS. 3 and 4. Thus, the particle filter of the processor 8 is, in effect, a Joint Probabilistic Data Association Filter (JPDAF), see for example “Monte Carlo filtering for multi-target tracking and data association”, J. Vermaak, S. J. Godsill, and P. Prez, IEEE Transactions on Aerospace and Electronic Systems, 41:309-332, 2005, which is incorporated herein by reference.
FIG. 2 is a schematic illustration (not to scale) of an example multi-target tracking scenario 1 in which the target tracker 2 may be implemented.
In this example, a plurality of target trackers 2 are used to track a first target 10 and a second target 12 as the targets 10, 12 move. In this example there are four target trackers 2 at different locations.
In the scenario 1, the first target 10 travels along a first track 14. Also, the second target 12 travels along a second track 16.
In the scenario 1, the first track 14 and the second track 16 overlap to some extent. The section of the tracks 14, 16 that overlaps is hereinafter referred to as the “overlap”, and is indicated in FIG. 1 by the reference numeral 18.
The following information is useful in understanding the description of a target tracking process used by a target tracker 2 to track a target 10, 12. This process is described later below. The target tracking process may be used to track the targets 10, 12 and/or infer unobserved states (such as current and future positions and/or velocities) of the targets 10, 12.
Each target 10, 12 may be described by a state variable xj that includes the two-dimensional position and velocity of a target. The state vector is as follows:

xj=[rj,x,rj,y,vj,x,vj,y]T 

where:

    
    
        j=1 for the first target 10;
        j=2 for the second target 12;
        rj,x is the position in the x-direction of the jth target;
        rj,y is the position in the y-direction of the jth target;
        vj,x is the velocity in the x-direction of the jth target; and
        vj,y is the velocity in the x-direction of the jth target.
    
    


A set of all single target states is X, such that xjεX.
A probabilistic dynamical or motion model of the targets 10, 12 specifies how likely they are to transition from one state xjk at a time tk to another state xjk+1 at a later time tk+1=tk+Δt. This is completely specified by the following conditional probability density:

p(xjk+1|xjk)

A positional measurement of a target 10, 12 that is generated by a position sensor 4 of the ith target tracker 2 is denoted by zi in this description. An identity measurement of a target 10, 12 that generated by an identity sensor 6 of the ith target tracker 2, is denoted by ζi in this description.
For clarity and ease of understanding, it is assumed in this description that a position sensor 4 may produce false negatives but not false positives, i.e. clutter measurements are not considered. Thus, for the position sensor 14 of the ith target tracker 2, a probabilistic detection model that reflects the probability over a binary detection variable Dikε{0,1} for a given target state xj is:

p(Dik|xj)

Observations are made independently for each of the targets.

This formulation allows information about the targets' states to be updated from both detection and non-detection events. Non-detection events imply that a target 10, 12 is not likely to be in a state where it can be easily detected.
For a non-detection event (i.e. Dik=0) the observation variable is empty, i.e. zi=Ø. For a detection event (i.e. Dik=1) the observation variable is a member of the set Z, which denotes a set of all possible bearings and/or ranges measurements.
The measurement received in a detection observation is modeled by the following conditional distribution:

p(zik=z|xj,Dik=1)

where zεZ.

Thus, a full probabilistic position sensor model that represents the distribution over the measurements zikεØ∩Z is:

p(zik|xj)

This sensor model can be expanded using the modeled detection probability:

p(zik=Ø|xj)=p(Dik=0|xj)

p(zik=z|xj)=p(Dik=1|xj)p(zik=z|xj,Dik=1)

This formulation advantageously allows detection and non-detection events to be treated in a unified manner.
A detection event may generate identity information regarding the identity of a target 10, 12, i.e. an identity sensor 6 may detect a target 10, 12. An identity observation by the identity sensor 6 of the ith target tracker is represented by the variable ζiε{Ø, 1, 2}. Here ζi=1 if the first target 10 is detected, ζi=2 if the second target 12 is detected, and ζi=Ø if the identity sensor 6 does not have the ability to discriminate targets, or an automatic target recognition algorithm processing the measurements could not produce an accurate identification.
For simplicity and ease of understanding, in this description if a detection (Dk=1) is generated, the identity observation is independent of the state of the target. Thus, a model of identity observations generated by the identity sensor 6 of the ith target tracker 2 is given by a distribution over ζi given the underlying identity of the target j that is detected, i.e.:

p(ζi|xj)=p(ζi|j)

where




  
    
      p
      ⁡
      
        (
        
          
            ζ
            i
          
          =
          
            φ
            |
            j
          
        
        )
      
    
    =
    
      1
      -
      
        β
        i
      
    
  




  
    
      p
      ⁡
      
        (
        
          
            ζ
            i
          
          =
          
            j
            |
            j
          
        
        )
      
    
    =
    
      
        α
        i
      
      ⁢
      
        β
        i
      
    
  




  
    
      p
      ⁡
      
        (
        
          
            
              ζ
              i
            
            ≠
            j
          
          |
          j
        
        )
      
    
    =
    
      
        
          (
          
            1
            -
            
              α
              i
            
          
          )
        
        
          M
          -
          1
        
      
      ⁢
      
        β
        i
      
    
  



where:

    
    
        βiε[0,1] represents how probable it is that a signature observation is generated by detection event;
        αiε[0,1] represents how confident the signature observations are; and
        M is the number of targets (i.e. two in the scenario 1).
    
    


It is possible to get positional information and no signature information, but it is not possible to signature information without positional information.
A joint state of the targets 10, 12 is defined as:

X=[x1T,x2T]T 

The positional and signature observation variables are expressed together as:

y=[z,ζ]

At a kth time step that the ith target tracker 2 takes a measurement, the sensors 4, 6 of the ith target tracker 2 generate a set of observations:

yik=[yi,1k,yi,2k, . . . , yi,Likk]

where: Lik is the number of detections generated by the sensors 4, 6 of the ith target tracker at time k; and

yi,lk=[zi,lk,ζi,lk]

In the scenario 1, the tracking of the targets 10, 12 requires estimating the following probability distribution:

p(Xk|Y1:k)

where:

    
    
        Y1:k=[Y1, . . . , Yk];
        Yk=[y1k, . . . , yNk]; and
        N is the number of observations at time k.
    
    


To perform this estimation, the present invention utilises a joint probabilistic data association filter (JPDAF).
FIG. 3 is a process flowchart showing certain steps of an embodiment of the tracking and identification process to estimate target states at a kth time-step. The process comprises a JPDA process carried out by the particle filters (i.e. the JPDAFs) of the processors 8 of the target trackers 2 and an identity management process.
In this embodiment, the tracking and identification process maintains a set of marginal estimates representing the state of each target 10, 12.
The process shown in FIG. 3 is for estimating target states at a kth time-step. In practice, this process may be performed for each of a series of time-steps to estimate target states at each of these time-steps using observations/measurement as they are acquired.
At step s2, estimates of the states of the targets at time k−1 are determined. In the recursive implementation of the process of FIG. 3, these estimates are the output of the previous iteration.
At time k−1 the state of each target jε{1, 2} is estimated separately by the following distribution:

p(xjk−1|Y1:k−1)

In this embodiment, the tracking and identification process is implemented using particle filters. The above defined prior estimate of the state of each target jε{1, 2} is given by a set of S weighted samples, such that:




  
    p
    ⁡
    
      (
      
        
          x
          j
          
            k
            -
            1
          
        
        |
        
          Y
          
            1
            :
            
              k
              -
              1
            
          
        
      
      )
    
  
  ≈
  
    
      ∑
      
        s
        =
        1
      
      S
    
    ⁢
    
      
    
    ⁢
    
      
        w
        
          j
          ,
          
            (
            s
            )
          
        
        
          k
          -
          1
        
      
      ⁢
      
        δ
        ⁡
        
          (
          
            
              x
              j
              
                k
                -
                1
              
            
            -
            
              x
              
                j
                ,
                
                  (
                  s
                  )
                
              
              
                k
                -
                1
              
            
          
          )
        
      
    
  




where xj,(s)k−1 is the sth sample point (which is assigned a weight wj,(s)k−1). In this embodiment, the weights are normalised.

At step s4, observations of the targets 10, 12 are taken by the target trackers at time k (or put another way, the target states are measured at time k).
At step s5, the estimates of the states of the targets at time k are predicted (or put another way, estimates of the target states for time k are determined). In this embodiment, the estimates of the states of the targets at time k are predicted using the estimates of the states of the targets at time k−1. These predicted states are given by the following:

p(xjk|Y1:k−1)=∫p(xjk|xjk−1)p(xjk−1|Y1:k−1)dxjk−1 

In this embodiment, a set of weighted particles, that represents the predicted estimate (i.e. the left hand side of above equation), are found using the sampling importance resampling method, also known as the bootstrap filter or condensation algorithm. In this method new sample points are generated using the prediction or motion model as the proposal distribution, i.e. for each sample point xj,(s)k−1 a new sample point is produced by sampling from the prediction model:

xj,(s)k ̃p(xjk|xj,(s)k−1)

These new sample points, and the original weights, approximate the predicted distribution of the jth target:

{xj,(s)k,wj,(s)k−1}s=1S ̃p(xjk|Y1:k−1)

At step s6, a process of updating the predicted target states is performed as described below with reference to FIG. 4. The description of step s7 will follow after the description of FIG. 4.
FIG. 4 is a process flowchart showing certain steps of the process of step s6.
At step s8, all feasible data associations at time k are determined. The set of all feasible data associations are denoted by Ak. An element of this set AkεAk represents a mapping from observations to targets and is given by:

Ak=[ai,1k,ai,1k, . . . , ai,Likk]

where: ai,lkε{1, 2} specifies the target that generated the lth observation yi,lk.

For the distributed sensors (of the distributed target trackers 2 of FIG. 2) the set of feasible data associations is determined by analysing each sensor's field of view and their overlap with the predicted target estimates.
At step s10, for a given data association Ak, an estimate of the jth target's state is updated using Bayes' rule by calculating a posterior probability distribution:

p(xjk|Y1:k,Ak)=η1p(xjk|Y1:k−1)p(Yk|xjk,Ak)

where η1 is a normalisation constant.

In this embodiment, utilising a particle filter to estimate the states of the targets, this step updates the weights of each of the sample points. These new posterior weights are found for each of the feasible data associations Ak and are computed using:

wj,(s)k,Ak=η2wj,(s)k−1p(Yk|xj,(s)k,Ak)

where η2 is a normalisation constant such that the computed weight are normalised.

The superscript Ak is used to explicitly show that these weights are produced by the given data association.
This process of updating the estimates from the observed information (step s6A) can be done by just considering the state or positional information, z, however it is advantageous to explicitly include the identity information ζ. This step can be accomplished by separating the joint observation likelihood into a product over the positional information and identity information captured by the sensors. In this embodiment the resulting included identity information is carried forward into steps s10, s12 and s14. This separation gives the following:




  
    p
    ⁡
    
      (
      
        
          
            Y
            k
          
          |
          
            x
            j
            k
          
        
        ,
        
          A
          k
        
      
      )
    
  
  =
  
    
      ∏
      
        ∀
        
          i
          :
          
            j
            ∈
            
              A
              i
              k
            
          
        
      
    
    ⁢
    
      
    
    ⁢
    
      
        p
        ⁡
        
          (
          
            
              z
              
                i
                ,
                
                  l
                  ij
                
              
              k
            
            |
            
              x
              j
              k
            
          
          )
        
      
      ⁢
      
        p
        ⁡
        
          (
          
            
              ζ
              
                i
                ,
                
                  l
                  ij
                
              
              k
            
            |
            j
          
          )
        
      
      ×
      
        
          ∏
          
            ∀
            
              i
              :
              
                j
                ∉
                
                  A
                  i
                  k
                
              
            
          
        
        ⁢
        
          
        
        ⁢
        
          p
          ⁡
          
            (
            
              
                z
                
                  i
                  ,
                  
                    l
                    ij
                  
                
                k
              
              =
              
                φ
                |
                
                  x
                  j
                  k
                
              
            
            )
          
        
      
    
  




where lij is the index of the observation made by the sensor of the ith target tracker 2 that is associated with the jth target, i.e. ai,lijk=j.

The first term of the above equation finds the product of the individual detection events that are associated with the jth target, explicitly incorporating the information regarding the identity of the target, while the second term incorporates the effects of a non-uniform probability of detection over the target state space.
At step s12, the probability of a given data association p(Ak|Y1:k−1) being true is approximated using the following:

p(Ak|Y1:k)=η3p(Yk|Ak,Y1:k−1)

where η3 is a normalisation constant.

This is equivalent to assuming the prior probability p(Ak|Z1:k) is uniform. In this embodiment, the association likelihood is determined by:




  
    
      
        
          p
          ⁡
          
            (
            
              
                
                  Y
                  k
                
                |
                
                  A
                  k
                
              
              ,
              
                Y
                
                  1
                  :
                  
                    k
                    -
                    1
                  
                
              
            
            )
          
        
        =
        
        ⁢
        
          ∫
          
            
              p
              ⁡
              
                (
                
                  
                    
                      Y
                      k
                    
                    |
                    
                      X
                      k
                    
                  
                  ,
                  
                    A
                    k
                  
                
                )
              
            
            ⁢
            
              p
              ⁡
              
                (
                
                  
                    X
                    k
                  
                  |
                  
                    Y
                    
                      1
                      :
                      
                        k
                        -
                        1
                      
                    
                  
                
                )
              
            
            ⁢
            
              d
              
                X
                k
              
            
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            ∏
            
              j
              =
              1
            
            M
          
          ⁢
          
            
          
          ⁢
          
            
              p
              ⁡
              
                (
                
                  
                    
                      Y
                      k
                    
                    |
                    
                      x
                      j
                      k
                    
                  
                  ,
                  
                    A
                    k
                  
                
                )
              
            
            ⁢
            
              p
              ⁡
              
                (
                
                  
                    x
                    j
                    k
                  
                  |
                  
                    Y
                    
                      1
                      :
                      
                        k
                        -
                        1
                      
                    
                  
                
                )
              
            
            ⁢
            
              d
              
                x
                j
                k
              
            
          
        
      
    
  




In this embodiment, the probability of a given data association is given by:




  
    p
    ⁡
    
      (
      
        
          A
          k
        
        |
        
          Y
          
            1
            :
            k
          
        
      
      )
    
  
  =
  
    
      η
      4
    
    ⁢
    
      
        ∏
        
          j
          =
          1
        
        M
      
      ⁢
      
        
      
      ⁢
      
        [
        
          
            ∑
            
              s
              =
              1
            
            S
          
          ⁢
          
            
          
          ⁢
          
            w
            
              j
              ,
              
                (
                s
                )
              
            
            
              k
              ,
              
                A
                k
              
            
          
        
        ]
      
    
  




where η4 is a normalisation constant.

At step s14, for each target j, the updated estimates of the target states are determined. In this embodiment, the separate distributions over the state of each target are recovered by marginalising the joint distribution. This produces the following:




  
    
      
        
          p
          ⁡
          
            (
            
              
                x
                j
                k
              
              |
              
                Y
                
                  1
                  :
                  k
                
              
            
            )
          
        
        =
        
        ⁢
        
          ∫
          
            
              p
              ⁡
              
                (
                
                  
                    X
                    k
                  
                  |
                  
                    Y
                    
                      1
                      :
                      k
                    
                  
                
                )
              
            
            ⁢
            
              d
              
                x
                1
                k
              
            
            ⁢
            
              d
              
                x
                2
                k
              
            
            ⁢
            
              
            
            ⁢
            ...
            ⁢
            
              
            
            ⁢
            
              d
              
                x
                M
                k
              
            
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          
            ∑
            
              
                A
                k
              
              ∈
              
                A
                k
              
            
          
          ⁢
          
            
          
          ⁢
          
            
              p
              ⁡
              
                (
                
                  
                    A
                    k
                  
                  |
                  
                    Y
                    
                      1
                      :
                      k
                    
                  
                
                )
              
            
            ⁢
            
              p
              ⁡
              
                (
                
                  
                    
                      x
                      j
                      k
                    
                    |
                    
                      Y
                      
                        1
                        :
                        k
                      
                    
                  
                  ,
                  
                    A
                    k
                  
                
                )
              
            
          
        
      
    
  




This marginalisation operation is the source of the approximation within the JPDAF process that removes correlation information between the states of the targets.
In this embodiment, this marginalisation process is performed to determine posterior weight values. This is done by summing over the different association hypotheses:




  
    w
    
      j
      ,
      
        (
        s
        )
      
    
    k
  
  =
  
    
      ∑
      
        
          A
          k
        
        ∈
        
          A
          k
        
      
    
    ⁢
    
      
    
    ⁢
    
      
        p
        ⁡
        
          (
          
            
              A
              k
            
            |
            
              Y
              
                1
                :
                k
              
            
          
          )
        
      
      ⁢
      
        w
        
          j
          ,
          
            (
            s
            )
          
        
        
          k
          ,
          
            A
            k
          
        
      
    
  




After the posterior weights have been determined, an optional additional process of resampling the sample points may be performed so that unlikely points are discarded and likely points are replicated.
This completes the description of FIG. 4 (i.e. the description of step s6). Referring back to FIG. 3, the step following step s6, i.e. step s7, will now be described.
At step s7 a process of identity management is performed as follows.
This identity management process advantageously alleviates problems of mixing up the correspondence between tracks 14, 16 and targets 10, 12, e.g. when the paths of targets cross (for example at the overlap 18), or simply become close to one another. This tends to be an important problem to address when tracking a single target, or reconstructing paths of targets.
In this embodiment, for ease of understanding, the process model of the targets 10, 12 is uninformative compared to the accuracy of the sensors 4, 6. Thus, it tends to be particularly computationally efficient to maintain a flattened probabilistic estimate of where the position of the targets 10, 12 and estimate the identities of the targets separately.
In other words, in this embodiment, the states of the targets 10, 12 are estimated using two particle filters, with the problem of data association solved using the JPDA algorithm (with or without the reformulation which explicitly incorporates the available identity information to improve data association hypothesis probabilities) and the mapping between estimated tracks and actual targets is additionally estimated. This additional estimation process will be referred to as the identity management process.
In this embodiment, the identity management process, i.e. the process of maintaining a distribution over target-to-track assignments, is as described in “The identity management Kalman filter”, B. Schumitsch, S. Thrun, L. Guibas, and K. Olukotun, Proc. of Robotics: Science and Systems, 2006, and “The information-form data association filter”, B. Schumitsch, S. Thrun, G. Bradski, and K. Olukotun, pages 1193-1200, 2006, which are incorporated by reference herein. However, in other embodiments a different identity management technique may be used, for example a technique which maintains a double stochastic belief matrix that represents the marginal probabilities of the individual track-to-target assignments or others as described in “The identity management problem, a short survey”, L. J. Guibas In Proc. Int. Conf. on Information Fusion, 2008.
In this embodiment, the identity management process explicitly parameterises the joint distribution of the targets over the factorially many assignments.
In this embodiment, an extra filter is used to maintain identity information. This filter is represented as a matrix Ω, hereinafter referred to as the “identity matrix”. Each element in this matrix represents the amount of accumulated evidence that a given filtered track represents a particular target and is stored as the log probability that a given track is assigned to a given target. By representing the track-to-target association by a permutation matrix P, the information matrix defines a probability distribution over P:

p(P)=exp(trace(ΩPT))

The most likely association between tracks and targets is determined by solving the following:




  
    
      
        
          P
          *
        
        =
        
        ⁢
        
          arg
          ⁢
          
            
          
          ⁢
          
            
              max
              P
            
            ⁢
            
              p
              ⁡
              
                (
                P
                )
              
            
          
        
      
    
  
  
    
      
        =
        
        ⁢
        
          arg
          ⁢
          
            
          
          ⁢
          
            
              
                max
                ⁢
                
                  
                
              
              P
            
            ⁢
            
              trace
              ⁡
              
                (
                
                  Ω
                  ⁢
                  
                    
                  
                  ⁢
                  
                    P
                    T
                  
                
                )
              
            
          
        
      
    
  




The identity matrix maintains the correlation between the joint state of the tracks 8, 10. This is achieved by enforcing each track estimate can only represent a single target 10, 12.
In this embodiment, to maintain the identity estimates of the tracks 8, 10, the identity matrix Ω is updated when a mixing event occurs that may lead to a track swap, and when new identity information about the identity of a target 10, 12 is observed by an identity sensor 16 of a target tracker 2. The terminology “mixing event” is used herein to refer to an event, a result of which is that it is possible that the correspondence between tracks and targets may be mixed up.
A mixing event may, for example, be caused by the two targets 10, 12 coming within proximity of one another (e.g. at the overlap 18) such that they cannot be accurately identified.
In this embodiment, a mixing event is defined by a doubly stochastic matrix Mk, herein after referred to as “the mixing matrix”, that transforms the identity matrix Ω as follows:

Ωk=log(Mk exp(Ωk−1)

where log(·) and exp(·) are element wise operators.

The ijth element of the mixing matrix defines how much of the marginal of the identity of the ith target contributes to the subsequent marginal of the jth target.
In this embodiment the mixing matrix is determined by examining the overlap in the underlying state space of the track estimates, such that the elements of the mixing matrix are derived from/based on the overlap in the underlying state space of the track estimates. In other words, an element of the mixing matrix is based upon the overlap (in the underlying state space of target state estimates) of the target(s) corresponding to that element. This advantageously tends to alleviate identity mixing between targets/tracks that have little overlap in state space but a large overlap in observation space. For instance, using the position sensors 4, the targets 10, 12 may have similar positional estimates, but large differences in velocities tend not to lead to confusion in their identities. Thus, only tracks with similar state estimates are mixed.
In this embodiment, the mixing matrix is constructed by considering all feasible pairwise swaps that can occur between the targets, as follows.
The mixing matrix M is constructed from the matrix {circumflex over (M)} through appropriate normalisation using the Sinkhorn algorithm that iteratively normalises the row and column sums to create a doubly stochastic matrix.
The matrix {circumflex over (M)} is defined as follows:

{circumflex over (M)}=[{circumflex over (m)}qr]

where each element {circumflex over (m)}qr is the overlap between the qth and rth target estimates.

In this embodiment, this overlap is calculated using the Bhattacharyya coefficient:

{circumflex over (m)}qr=∫√{square root over (pq(x)pr(x))}{square root over (pq(x)pr(x))}dx 

where pt(x)=p(xtk=x|Y1:k) for t=q, r. In other embodiments, the overlap calculation may be implemented in any appropriate fashion, for example minor variations may be made to the above equation, or different equations may be used. These include the Hellinger distance defined as:





  
    
      1
      2
    
    ⁢
    
      ∫
      
        
          
            (
            
              
                
                  
                    p
                    q
                  
                  ⁡
                  
                    (
                    x
                    )
                  
                
              
              -
              
                
                  
                    p
                    r
                  
                  ⁡
                  
                    (
                    x
                    )
                  
                
              
            
            )
          
          2
        
        ⁢
        
          d
          x
        
      
    
  
  ;




and the generalised Renyi divergence:





  
    1
    
      α
      -
      1
    
  
  ⁢
  log
  ⁢
  
    
  
  ⁢
  
    ∫
    
      
        
          
            [
            
              
                p
                q
              
              ⁡
              
                (
                x
                )
              
            
            ]
          
          α
        
        
          
            [
            
              
                p
                r
              
              ⁡
              
                (
                x
                )
              
            
            ]
          
          
            α
            -
            1
          
        
      
      ⁢
      
        
          d
          x
        
        .
      
    
  




In this embodiment, when estimates completely overlap, the element corresponding to the relevant targets equals one. Also, when there is no overlap, the relevant element equal zero.
The use of the above described mixing matrix advantageously tends to allow for an accurate description of how targets may mix (or the track labels may switch).
Thus, a process of multi-target tracking utilising a JPDA process implemented using particle filters and augmented with an identity management process is provided.
The present invention tends to be advantageously capable of allowing for the tracking and identification of individual targets in situations where the respective tracks of those targets overlap or are proximate to each other. For example, the present invention is particularly useful for tracking and identifying road vehicles moving along a road.
Moreover, the present invention advantageously utilises a particle filter to overcome a problems of the JPDA filter prediction and update equations tending to be difficult or impossible to solve analytically, without imposing certain restrictions (e.g. Gaussian uncertainty and linear prediction and motion models).
A further advantage of the present invention is that it tends to enable the identifying of targets whose tracks overlap or are proximate to one another resulting in uncertainty in the true identity of the target. This is advantageous whether identity information is available or not. Moreover, these identities may be managed affectively.
Apparatus, including the one or more of the processors 18 of the target trackers 2, for implementing the above arrangement, and performing the method steps described above, may be provided by configuring or adapting any suitable apparatus, for example one or more computers or other processing apparatus or processors, and/or providing additional modules. The apparatus may comprise a computer, a network of computers, or one or more processors, for implementing instructions and using data, including instructions and data in the form of a computer program or plurality of computer programs stored in or on a machine readable storage medium such as computer memory, a computer disk, ROM, PROM etc., or any combination of these or other storage media.
Furthermore, the above embodiments may be implemented using one or more centralised processors that are provided with measurements/observations taken by sensors that are arranged to measure the targets at different points along their tracks. For example, the above embodiments may be implemented using a centralised processor module remote from position and identity sensors, the sensors being positioned at the locations of the target trackers 2 of FIG. 2.
It should be noted that certain of the process steps depicted in the flowcharts of FIGS. 3 and/or 4 and described above may be omitted or such process steps may be performed in differing order to that presented above and shown in the Figures. Furthermore, although all the process steps have, for convenience and ease of understanding, been depicted as discrete temporally-sequential steps, nevertheless some of the process steps may in fact be performed simultaneously or at least overlapping to some extent temporally.
In the above embodiments, the mixing matrix is advantageously used to describe the mixing events for performing effective identity management in a scenario in which particle filters are used to implement a JPDA process. However, in other embodiments, one or more different types of filters (e.g. a Kalman filter) may be used instead of or in addition to a particle filter. Also, in other embodiments the filters may be used to implement a different process (e.g. Multi-hypothesis tracking, or Maximum Likelihood Data Association) instead of or in addition to the JPDA process described above. Furthermore, in other embodiments, the mixing matrix is used for purposes other than identity management using an information matrix. For example, in other embodiments the mixing matrix may implemented in a process that utilises a belief matrix, such as that described in “Multiple-target tracking and identity management” I. Hwang, H. Balakrishnan, K. Roy, J. Shin, L. Guibas, C. Tomlin. Proc. IEEE Sensors, 1:36-41, October 2003 (incorporated herein by reference).
In the above embodiments, the present invention is used to tracks two targets. However, in other embodiments, the present invention may be used track any number of targets.
In the above embodiments, four target trackers are used to track the targets. However, in other embodiments any number of target trackers may be used.
In the above embodiments, targets are described by a state variable that includes the two-dimensional position and velocity of a target. However, in other embodiments the state vector of a target may comprise any appropriate parameters that described a state of a target, for example, acceleration or three-dimensional position and/or velocity.
In the above embodiments, information of failed detection (“mis-detection information”) is explicitly included in the observation likelihood. However, in other embodiments this mis-detection information is treated differently, for example it can be ignored and the information regarding failed detections not used by the tracker.
In the above embodiments, the target trackers each comprise position and identity sensors. However, in other embodiments one or more of the target trackers comprises any number of any appropriate type of sensors that are used to determine a state of a target and/or an identity of a target.
In the above embodiments, for ease of understanding, the different targets are observed independently of each other and there are no interference effects when multiple targets are observed at the same time. In other embodiments, such this is not case and problems arising from such effects are alleviated, for example by utilising a multi-target sensor model that determines the probability over observations given the states of multiple targets.
In the above embodiments, for ease of understanding, a single track is associated with a single target. However, in other embodiments the above described formulation of the JPDA process is extended to advantageously account for situations in which a single filtered track represents different targets at different times. For example, this may be accomplished by replacing the target index j with a track index t in all the filtering equations. However, in situations in which the identity of a target is not known (i.e. only estimated), the process tends to require extending. This is done by, for example, computing a set of marginal distributions over the track to target assignments. This results in a doubly stochastic belief matrix. This matrix belief matrix can be constructed from the information matrix using efficient Monte Carlo Markov chain sampling methods.

## Claims

1. A method of tracking a plurality of targets, the method comprising:
for a first time-step, estimating a state of each target;
at a second time-step:
measuring values for a state of a target; and
measuring and updating a value of an identity parameter for the target, wherein the identity parameter at least partially distinguishes the target from one or more of the other targets;

for the second time-step, estimating a state of each target using the estimated target states for the first time-step;
updating the estimated target states for the second time-step using the measured values for the states of the targets; and
performing an identity management process to estimate a probability that a particular state measurement corresponds to a particular target, wherein the identity management process comprises: providing a mixing matrix, an element of the mixing matrix being based on an overlap between the updated estimated target states of targets in an underlying state space of the estimated target states; and
wherein the providing of the mixing matrix comprises determining a value for an element of the missing matrix as a function of the state of the particular target at the second time-step, and a set of all state values and identity parameter values measured at time steps up to and including the second time-step.

2. A method according claim 1, wherein the providing of the mixing matrix comprises:
determining values for elements {circumflex over (m)}qr, where:


{circumflex over (m)}qr=∫√{square root over (pq(x)pr(x))}{square root over (pq(x)pr(x))}dx; 



pi(x)=p(xi=x|Y) for i=q, r; 

xi is a state of an ith target at the second time-step; and
Y is a set of all state values and identity parameter values measured at time steps up to and including the second time-step; and

forming the elements into matrix {circumflex over (M)}=[{circumflex over (m)}qr].

3. A method according to claim 2, wherein determining the mixing matrix comprises:
normalising the matrix {circumflex over (M)}.

4. A method according to claim 1, wherein the updating of the estimated target states for the second time-step using the measured values is performed using one or more of: a multi-hypothesis tracking process, a maximum likelihood data association process, and a joint probabilistic data association process.

5. A method according to claim 1, wherein the updating of the estimated target states for the second time-step using the measured values is implemented using one of a particle filter and a Kalman filter.

6. A method according to claim 1, wherein motion of the targets to be tracked is constrained.

7. A method according to claim 6, wherein the targets are land-based vehicles.

8. A method according to claim 1, wherein a path of one of the targets is proximate to a path of at least one different target.

9. A method according to claim 8, wherein the path of one of the targets overlaps a path of at least one different target.

10. A method according to claim 1 wherein estimating a state of each target is performed with a particle filter.

11. A method according to claim 1, wherein a state of a target comprises a value for a position of the target and/or a value for a velocity of the target.

12. Apparatus for tracking a plurality of targets, the apparatus comprising:
one or more processors arranged for a first time-step, to estimate a state of each target; and
one or more sensors arranged, at a second time-step, to:
measure values for a state of a target; and
measure and update a value of an identity parameter for the target, wherein the identity parameter at least partially distinguishes the target from one or more of the other targets;

wherein the one or more processors are configured:
for the second time-step, to estimate a state of each target using the estimated target states for the first time-step;
to update the estimated target states for the second time-step using the measured values for the states of the targets; and
to perform an identity management process to estimate a probability that a particular state measurement corresponds to a particular target; wherein the identity management process includes providing a mixing matrix, an element of the mixing matrix being based on an overlap between the updated estimated target states of targets in an underlying state space of the estimated target states; and
wherein the providing of the mixing matrix comprises determining a value for an element of the mixing matrix as a function of the state of the particular target at the second time-step, and a set of all state values and identity parameter values measured at time steps up to and including the second time-step.

13. Apparatus according to claim 12, wherein elements {circumflex over (m)}qr of the mixing matrix have the following values:


{circumflex over (m)}qr=∫√{square root over (pq(x)pr(x))}{square root over (pq(x)pr(x))}dx; 

where:


pi(x)=p(xi=x|Y) for i=q, r; 

xi is a state of an ith target at the second time-step; and
Y is a set of all state values and identity parameter values measured at time steps up to and including the second time-step; and wherein the elements are formed into matrix {circumflex over (M)}=[{circumflex over (m)}qr].

14. A non-transitory computer program product encoded with instructions that, when executed by one or more processors, causes a process to be carried out, the process comprising:
for a first time-step, estimating a state of each target;
at a second time-step:
measuring values for a state of a target; and
measuring and updating a value of an identity parameter for the target, wherein the identity parameter at least partially distinguishes the target from one or more of the other targets;

for the second time-step, estimating a state of each target using the estimated target states for the first time-step;
updating the estimated target states for the second time-step using the measured values for the states of the targets; and
performing an identity management process to estimate a probability that a particular state measurement corresponds to a particular target, wherein the identity management process comprises: providing a mixing matrix, an element of the mixing matrix being based on an overlap between the updated estimated target states of targets in an underlying state space of the estimated target states; and
wherein the providing of the mixing matrix comprises determining a value for an element of the mixing matrix as a function of the state of the particular target at the second time-step, and a set of all state values and identity parameter values measured at time steps up to and including the second time-step.

15. The computer non-transitory program product according to claim 14, wherein the updating of the estimated target states for the second time-step using the measured values is performed using one or more of: a multi-hypothesis tracking process, a maximum likelihood data association process, and a joint probabilistic data association process.

16. The computer non-transitory program product according to claim 14, wherein the updating of the estimated target states for the second time-step using the measured values is implemented using one of a particle filter and a Kalman filter.

17. The computer non-transitory program product according to claim 14, wherein motion of the targets to be tracked is constrained.

18. The computer non-transitory program product according to claim 14, wherein a path of one of the targets is proximate to a path of at least one different target.

19. The computer non-transitory program product according to claim 18, wherein the path of one of the targets overlaps a path of at least one different target.

20. The computer non-transitory program product according to claim 14, wherein estimating a state of each target is performed with a particle filter.