Multi-target data processing for multi-static and multi-channel passive radars

• US 8,654,007 B2
• Filed: 08/28/2012
• Issued: 02/18/2014
• Est. Priority Date: 07/09/2009
• Status: Active Grant

• Current Assignee: Thales SA
• Original Assignee: Thales SA
• Inventors: Allam, Sbastien, De Gramont, Emmanuel, Deltour, Jean-Claude, Klein, Mathieu, Greverie, Wilfried
• Primary Examiner(s): Sotomayor, John B
• Assistant Examiner(s): Windrich, Marcus
• Application Number:
  US13/596,327
• Publication Number:
  US 20130076558A1
• Time in Patent Office: 539 Days
• Field of Search: 342/59, 342 89- 90, 342107-108, 342/463
• US Class Current: 342/195
• CPC Class Codes: G01S 13/003 Bistatic radar systems; Multistatic radar systems
• CPC Class Codes: G01S 13/726 Multiple target tracking
• CPC Class Codes: G01S 13/878 Combination of several spaced transmitters or receivers of known location for determining the position of a transponder or a reflector G01S13/874 takes precedence

Abstract

The present invention relates to the general field of radar tracking applied to multi-static radar systems. It consists of a method for coherently merging the individual tracks generated from the various bistatic bases that make up the system so as to generate and maintain global tracks so that an object detected by a single bistatic base is represented by the corresponding track and an object detected by a number of bistatic bases and giving rise to the creation of a number of individual tracks is represented only by a single global track.

Claims

1. A tracking method for multi-channel radar systems, the system comprising a receiver exploiting the transmissions of a number of utility transmitters, each transmitter forming with the receiver a bistatic base, the transmission produced by each of the transmitters exhibiting a spectrum consisting of different frequency whose treatment leads to each measurement instant in the creation of blips, each blip being characterized by a vector of three measurement constituted by a bistatic distance measurement, a bistatic speed measurement and an azimuth measurement, as well as a 3 ×

3 diagonal matrix of the errors associated to these measurements, comprising;

a first upstream merging step;

applied separately at each bistatic base, operating directly in the space of the components of the vectors characterizing the blips, and implementing an operation to form individual tracks, carrying the identification of the blips of different frequency channels relative to the same target and the instantaneous merging of these blips, the operation to form individual tracks comprising for each measurement instant;

a first blip merging substep for associating the individual blips of the various reception channels of each bistatic base which exhibit, for a given measurement instant, measurement vectors y_j that are close to one another and forming merged blips, each merged blip j being characterized by a measurement vector yj =(rb, vb, θ

) where rb, vb and θ

respectively represent the bistatic distance measurement, the bistatic speed measurement and the azimuth measurement of the detected object by the considered bistatic base;

a second substep for constructing individual tracks from the merged blips, an individual track being formed by the succession over time of the merged blips relating to one and the same target, each detected target being represented by a bistatic individual track which characteristics change over time depending on blips associated to it;

the blip merging substep using a blip aggregation module which associates the individual blips of various reception channels of one and the same bistatic base which exhibit, for a given measurement instant, measurement vectors y_j that are close to one another, taking the measurements accuracy into account, and forms merged blips, the aggregation of two individual blips p₁ and p₂ is achieved if the measurement vectors y1 and y2 characterizing these blips satisfy the proximity criterion defined by the following three relations;

2. The method of claim 1, wherein the blip merging substep further comprises the following modules:

a module for merging the aggregated individual blips which generates, for each aggregate , a single resultant measurement vector y(r_b, V_b, θ

) forming a merged blip, this measurement vector being made up from the measurement vectors yj of the aggregated individual blips, the merged blip being used to form a given individual track; and

a module which checks the state of the individual tracks formed to determine the likelihood of an individual track considering the number of individual blips forming the aggregate from which the merged blip used to form the track originate.

3. The method of claim 1, wherein the blip merging substep further comprises the following modules:

a module for merging the aggregated individual blips which generates, for each aggregate, a single resultant measurement vector y(r_b, V_b,θ

), from measurement vectors y_{j (rb}, V_b, θ

) of the blips forming the aggregate, y being defined by the following equalities;

5. The method of claim 1, wherein the module for initializing new global tracks of the downstream merging step successively comprises the following five steps:

a first step which tests the azimuth consistency of the candidate individual tracks;

a second step which tests the distance consistency of the candidate individual tracks;

a third step which tests the speed consistency of the candidate individual tracks;

a fourth step which proceeds with the initialization of a global track;

a fifth step which consists in aggregating, with the initialized global track, individual tracks coming from the bistatic bases not involved in the initialization of the global track.

7. The method of claim 1, wherein the second downstream merging step consists in nearing the zone of intersection of the location ellipsoids by a juxtaposition of blocks, each block corresponding to a parallelepiped characterized by the position of its center and the width of each of its sides, the union of these blocks encompassing the area of intersection, and extracting from the set of blocks which include the area, an average position value and an associated covariance matrix, each block being associated with a probability characterizing its coherence with bistatic distances to the origin of the zone of intersection.

9. The method of claim 1, wherein the module to maintain the formed global tracks of the second downstream merging step implements the calculation of the statistical distance γ

separating the projection of a global track in a bistatic base b and an individual track of this same bistatic base b, candidate to the merging, the statistical distance γ

being function of the bistatic components, distance, speed and azimuth characterizing the association-candidate individual track and the components of the projection of the state x of the global track in the bistatic base b.

11. The method of claim 1, wherein the non-linear filtering module of the second downstream merging step proceeds, for each measurement instant k, with the prediction {tilde over (X)}_k of the components of the state vector x corresponding to each global track based on the estimation {circumflex over (X)}_k−

1 of the components of the state vector x at the instant k−

1, this estimation itself being performed on the basis of the observation vector Y<sub>k−

1</sub> generated by means of the measurements relating, for the instant k−

1, to the individual tracks associated with the global track concerned and the prediction {tilde over (X)}<sub>k−

1</sub> of the components of the state vector x at the instant k−

1, the components of the observation vector consisting of the result of the merging of the bistatic distances of each individual track, generated in the locating module for the global tracks of the downstream merging step and for the bistatic speeds of the same individual tracks obtained on completion of the upstream merging step.

13. The method of claim 1, wherein the non-linear filtering module of the second downstream merging step proceeds, for each measurement instant k, with the prediction {tilde over (X)}_k of the components of the state vector x corresponding to each global track based on the estimation {circumflex over (X)}_k−

1 the components of the state vector x at the instant k−

1, this estimation itself being performed on the basis of the observation vector Y<sub>k−

1</sub> generated by means of the measurements relating, for the instant k−

1, to the individual tracks associated with the global track concerned and the prediction {tilde over (X)}<sub>k−

1</sub> of the components of the state vector x at the instant k−

1, the components of the observation vector consisting of the result of the merging of the bistatic distances of each individual track, generated in the locating module for the global tracks of the downstream merging step and for the bistatic speeds of the same individual tracks obtained on completion of the upstream merging step, the estimation {circumflex over (X)}_k and the prediction {tilde over (X)}<sub>k−

1</sub> of the observation vector x being determined by means of an extended Kalman filter.

4. A tracking method for multi-channel radar systems, the system comprising a receiver exploiting the transmissions of a number of utility transmitters, each transmitter forming with the receiver a bistatic base, the transmission produced by each of the transmitters exhibiting a spectrum consisting of different frequency whose treatment leads to each measurement instant in the creation of blips, each blip being characterized by a vector of three measurement constituted by a bistatic distance measurement, a bistatic speed measurement and an azimuth measurement, as well as a 3 ×

3 diagonal matrix of the errors associated to these measurements, comprising;

a first upstream merging step, applied separately at each bistatic base, operating directly in the space of the components of the vectors characterizing the blips, and implementing an operation to form individual tracks, carrying the identification of the blips of different frequency channels relative to the same target and the instantaneous merging of these blips, the operation to form individual tracks comprising for each measurement instant;

a first blip merging substep for associating the individual blips of the various reception channels of each bistatic base which exhibit, for a given measurement instant, measurement vectors y_j that are close to one another and forming merged blips, each merged blip j being characterized by a measurement vector yj =(rb, vb, θ

) where rb, vb and θ

respectively represent the bistatic distance measurement, the bistatic speed measurement and the azimuth measurement of the detected object by the considered bistatic base;

a second substep for constructing individual tracks from the merged blips, an individual track being formed by the succession over time of the merged blips relating to one and the same target, each detected target being represented by a bistatic individual track which characteristics change over time depending on blips associated to it;

the blip merging substep using;

a first blip aggregation module which associates the individual blips of various reception channels of one and the same bistatic base which exhibit, for a given measurement instant, measurement vectors y_j that are close to one another, taking the measurements accuracy into account, and forms merged blips,

a second module for merging the aggregated individual blips which generates, for each aggregate, a single resultant measurement vector y(r_b, V_b, θ

), from measurement vectors y_j(r_b, V_b, θ

) of the blips forming the aggregate, y being defined by the following equalities;

6. The method of claim 4, wherein the module for initializing new global tracks of the downstream merging step successively comprises the following five steps:

a first step which tests the azimuth consistency of the candidate individual tracks;

a second step which tests the distance consistency of the candidate individual tracks;

a third step which tests the speed consistency of the candidate individual tracks;

a fourth step which proceeds with the initialization of a global track;

a fifth step which consists in aggregating, with the initialized global track, individual tracks coming from the bistatic bases not involved in the initialization of the global track.

8. The method of claim 4, wherein the second downstream merging step consists in nearing the zone of intersection of the location ellipsoids by a juxtaposition of blocks, each block corresponding to a parallelepiped characterized by the position of its center and the width of each of its sides, the union of these blocks encompassing the area of intersection, and extracting from the set of blocks which include the area, an average position value and an associated covariance matrix, each block being associated with a probability characterizing its coherence with bistatic distances to the origin of the zone of intersection.

10. The method of claim 4, wherein the module to maintain the formed global tracks of the second downstream merging step implements the calculation of the statistical distance γ

separating the projection of a global track in a bistatic base b and an individual track of this same bistatic base b, candidate to the merging, the statistical distance γ

being function of the bistatic components, distance, speed and azimuth characterizing the association-candidate individual track and the components of the projection of the state x of the global track in the bistatic base b.

12. The method of claim 4, wherein the non-linear filtering module of the second downstream merging step proceeds, for each measurement instant k, with the prediction {tilde over (X)}_k of the components of the state vector x corresponding to each global track based on the estimation {circumflex over (X)}_k−

1 of the components of the state vector x at the instant k−

1, this estimation itself being performed on the basis of the observation vector Y<sub>k−

1</sub> generated by means of the measurements relating, for the instant k−

1, to the individual tracks associated with the global track concerned and the prediction {tilde over (X)}<sub>k−

1</sub> of the components of the state vector x at the instant k−

1, the components of the observation vector consisting of the result of the merging of the bistatic distances of each individual track, generated in the locating module for the global tracks of the downstream merging step and for the bistatic speeds of the same individual tracks obtained on completion of the upstream merging step.

14. The method of claim 4, wherein the non-linear filtering module of the second downstream merging step proceeds, for each measurement instant k, with the prediction {tilde over (X)}_k of the components of the state vector x corresponding to each global track based on the estimation {circumflex over (X)}_k−

1 of the components of the state vector x at the instant k−

1, this estimation itself being performed on the basis of the observation vector Y<sub>k−

1</sub> generated by means of the measurements relating, for the instant k−

1, to the individual tracks associated with the global track concerned and the prediction {tilde over (X)}<sub>k−

1</sub> of the components of the state vector x at the instant k−

1, the components of the observation vector consisting of the result of the merging of the bistatic distances of each individual track, generated in the locating module for the global tracks of the downstream merging step and for the bistatic speeds of the same individual tracks obtained on completion of the upstream merging step, the estimation {circumflex over (X)}_k and the prediction {tilde over (X)}<sub>k−

1</sub> of the observation vector x being determinated by means of an extended Kalman filter.