FR 2979709 A1 20130308 FR 1102670 A 20110902 PROCEDE POUR INITIALISER DES PISTES CARTESIENNES A PARTIR DE MESURES BISTATIQUES REALISEES PAR UN OU PLUSIEURS RECEPTEURS D'UN SYSTEME RADAR MULTISTATIQUE. L'invention concerne le domaine général du pistage dans le cadre des radars multistatiques. Elle concerne plus particulièrement le pistage multi-cibles réalisé à partir de systèmes multistatiques passifs. L'exploitation d'un système radar passif comportant une pluralité d'émetteurs et/ou de récepteurs comporte un certain nombre de difficultés. En effet chaque couple émetteur-récepteur permet de constituer une base bistatique à partir de laquelle il est possible de réaliser, de différentes manières simples et connues, la détection et le pistage bistatiques d'objets évoluant dans l'espace couvert par cette base. Cependant, lorsque l'on souhaite connaître la situation aérienne globale de la zone couverte par le système radar, il est nécessaire d'exploiter de manière conjointe les informations délivrées par les différentes bases bistatiques, généralement exploitées sous forme de pistes élémentaires, de façon à fusionner ces informations pour constituer des pistes globales ou pistes fusionnées. Une exploitation conjointe des informations produites par les différentes bases bistatiques consiste en pratique à former et à entretenir des pistes cartésiennes globales en associant les informations de détection, les plots bistatiques, produites par les différentes bases bistatiques composant le système de détection considéré, pour former et entretenir des pistes cartésiennes. Le document de la demande de brevet français déposé par la demanderesse le 01/09/2009, sous le numéro 0904144, décrit un procédé permettant une telle exploitation conjointe. Ce procédé est illustré par la figure 1. Le procédé décrit permet de former des pistes cartésiennes à partir de plots bistatiques produits par des bases bistatiques distinctes, dans le cadre d'un système de détection comportant plusieurs récepteurs, deux bases bistatiques pouvant alors être formées par le même récepteur ou par des récepteurs différents. Comme l'illustre la figure 1 et comme le mentionne le texte de la demande de brevet mentionnée précédemment, la formation et l'entretien de ces pistes cartésiennes passe naturellement par une ou plusieurs étapes d'initialisation. Diverses méthodes connues peuvent être utilisées pour réaliser l'initialisation de pistes cartésiennes à partir de plots bistatiques produit par les différentes bases bistatiques, formées à partir de récepteurs différents. Ces méthodes sont généralement intimement liées en ce qui concerne leur principe à la structure du système considéré, et à la façon dont les plots 10 bistatiques produits par les différentes bases bistatiques sont associés. Cependant, quelle que soit la méthode utilisée l'initialisation de pistes cartésiennes se heurte à la difficulté que constitue la prise en compte d'un nombre très important de plots bistatiques disponibles à chaque instant pour tenter d'initialiser des pistes. 15 Un but de l'invention est de proposer une solution permettant l'initialisation de pistes cartésiennes, solution adaptée au fonctionnement d'un système de détection multistatique passif comportant plusieurs récepteurs. Autrement dit un but de l'invention est de permettre une 20 exploitation optimale des plots bistatiques produits par les différentes bases bistatiques constituant le système, que ces plots soient produits par des bases bistatiques formées par un même récepteur ou par des bases bistatiques formées des récepteurs différents, pour réaliser l'initialisation de pistes cartésiennes. Par exploitation optimale on entend une exploitation 25 permettant d'éviter d'initialiser une piste à partir de chaque plot bistatique formé. A cet effet l'invention a pour objet un procédé pour réaliser l'association de plots bistatiques correspondant à une même cible détectée par un système radar multistatique comportant une pluralité d'émetteurs et 30 de récepteurs distants les uns des autres, un couple (émetteur, récepteur) formant une base bistatique du système, chaque plot bistatique Pi correspondant à l'écho provenant du signal émis par un émetteur Tx1 et reçu par un récepteur Ek et étant caractérisé par l'émetteur Txi ayant émis le signal dont il est l'écho. Les plots bistatiques formés étant associés aux différents émetteurs du système pour former un ensemble de couples (plot, émetteur), le procédé selon l'invention comporte: - un premier module de formation de triplets de couples (Pi, Txj), ces triplets, regroupant des couples (plot, émetteur) qui présentent une 5 compatibilité de distance et de vitesse; - un second module de formation de n-uplets de couples (Pi, Txj), ces n-uplets regroupant des triplets formés à l'étape précédente répondant à des critères de regroupement donnés qui caractérisent la vraisemblance de l'attribution des différents triplets à une même cible; 10 Les n-uplets ainsi constitués formant des paquets de plots qui définissent une même cible. Selon un mode particulier de mise en oeuvre du procédé selon l'invention, la première étape de formation de triplets de couples (plot, 15 émetteur) est une étape itérative qui considère à chaque itération trois couples (plot, émetteur), un premier couple (P1, Txi), un second couple (P2, Tx2) et un troisième couple (Préf, TXréf) pris comme un couple de référence et qui leur applique une opération d'association qui teste la possibilité d'associer le premier couple (P1, Tx1) et le second couple (P2, Tx2) au couple 20 de référence (Préf, TXréf), un couple (plot, émetteur) étant associable au couple de référence si les caractéristiques de la cible correspondant au plot considéré sont compatibles des caractéristiques de la cible correspondant au plot formant le couple de référence; trois couples (Préf, TXréf), (P1, TX1) et (P2, Tx2), pour lesquels le test d'association est réussi formant un triplet de 25 couples. Selon une variante de ce mode de mise en oeuvre, la première étape de formation de triplets de couples (plot, émetteur) comporte elle-même les opérations suivantes: 30 - un premier test d'association pour déterminer si les couples (Préf, TXréf) et (P1, Txi) sont associables; - un deuxième test d'association pour déterminer si les couples (Préf, TXréf) et (P2, Tx2) sont associables; - un troisième test de compatibilité des couples (P1, Txi) et (P2, Tx2) - un quatrième test de compatibilités des distances bistatiques des plots - un cinquième test de test de la vraisemblance de la vitesse cartésienne de la cible correspondant aux plots associés; - une sixième opération de mémorisation des triplets formés de couples (P1, Tx1) et (P2, Tx2) ayant réussi tous les tests; le test d'association mis en oeuvre lors du premier test étant identique à celui mis en oeuvre lors du second test. Selon une forme particulière de cette variante de mise en oeuvre, le test d'association mis en oeuvre lors du premier test consiste, dans le cas où les plots P1 et Préf sélectionnés ont un même récepteur Rx pour origine, à mettre en oeuvre les quatre tests indépendants suivants: - un premier test (test 1) consistant à déterminer si les vitesses 15 bistatiques Vbk(Pi) et Vbk(Préf) des plots P1 et Préf peuvent correspondre à une même vitesse de cible, le premier test étant négatif si l'on a: Vbk (Pré, ) - Vbk (P1 )I 2.Vmax Où V. représente la vitesse maximale des cibles que l'on cherche à détecter; 20 - Un deuxième test (test 2) consistant à déterminer si les valeurs des azimuts associés aux plots P1 et Préf peuvent correspondre à une même cible, compte tenu des précisions aek(Pi) et aek(Préf) des mesures d'azimuts réalisées; le deuxième test étant négatif si l'on a: lek (Préf) (P1) > 3.\/(50k (Préf ) (520k (P1) ; 25 - Un troisième test (test 3) consistant à déterminer si les ellipsoïdes correspondant aux distances bistatiques des plots P1 et Préf présentent des intersections, le troisième test étant négatif si l'on a: (Dbk (Préf )- Dbk(P1))2 - FTxréf î(Tx112 > o OÙ [X-rxi YTx1 ZTx1f et XTxréf = [XTxréf YTxréf ZTxréf 30 correspondent ici respectivement à la position cartésienne de l'émetteur Tx1 associé à P1 et à la position cartésienne de l'émetteur Txréf associé à Préf; Un Quatrième test (test 4) consistant à déterminer si les intersections des ellipsoïdes construites à partir des distances bistatiques caractérisant les plots P1 et Préf, sont compatibles des mesures d'azimuts relatives à ces plots le quatrième test étant négatif si l'on a: [cos cp - 3.6cos,p, cos 3.6cos(p J n [0,1] # 0, où cp correspond au site selon lequel la cible potentielle correspondant 5 aux plots considérés peut être observée et acos(p l'incertitude affectant la détermination de cosy); le test d'association étant positif si les résultats des quatre tests sont positifs. 10 Selon une forme particulière de cette forme particulière, le test d'association mis en oeuvre lors du second test consiste, dans le cas où les plots P2 et Préf sélectionnés ont un même récepteur Rx pour origine, à mettre en oeuvre, sur les plots P2 et Préf, les quatre tests indépendants mis en oeuvre lors du premier test. 15 Selon une autre forme particulière de la variante de mise en oeuvre décrite précédemment, le test d'association mis en oeuvre lors du premier test consiste, dans le cas où les plots P1 et Préf sélectionnés ont pour origine des récepteurs Rx1 et Rxréf distincts, à mettre en oeuvre consécutivement les 20 quatre tests suivants: - un premier test (test A) consistant à déterminer s'il est possible de définir, dans la limite de portée du système de détection considéré, une intersection des directions des azimuts des observations correspondant aux plots P1 et Préf; et, dans l'affirmative, à calculer cette intersection, en prenant 25 en compte les incertitudes des mesures des azimuts de sorte que l'on définit une aire d'intersection; - un deuxième test (test B) consistant, si le test A est positif, à déterminer si les deux sphères contenant les ellipsoïdes définis par les distances bistatiques des plots P1 et Préf présentent une intersection; 30 - un troisième test (test C) consistant, si le test B est positif, à déterminer si les deux ellipsoïdes définis par les distances bistatiques des plots P1 et Préf présentent une intersection. - un quatrième test (test D) consistant, si le test C est positif, à déterminer si les mesures des distances bistatiques associées aux plots P1 35 et Préf sont compatible de la zone d'intersection des directions des azimuts des observations correspondant aux plots P1 et Préf déterminée lors du test A; le test d'association étant positif si les résultats des quatre tests sont positifs. Selon une autre forme particulière de cette forme particulière de mise en oeuvre, le test d'association mis en oeuvre lors du second test consiste, dans le cas où les plots P2 et Préf sélectionnés ont pour origine des récepteurs Rx2 et Rxréf distincts, à mettre en oeuvre consécutivement, sur les plots P2 et Préf, les quatre tests mis en oeuvre lors du premier test (41). Selon une autre forme particulière de mise en oeuvre, le troisième test (test C) est réalisé par la méthode du multiplicateur de Lagrange appliquée aux deux ellipsoïdes considérés. Selon une autre forme particulière de la variante de mise en oeuvre décrite précédemment, dans le cas où les plots P1 et P2 sélectionnés ont un même récepteur Rx pour origine, le troisième test de compatibilité des couples (P1, Txi) et (P2, Tx2) est réalisé en mettant en oeuvre, sur les plots P1 et P2, des tests analogues à ceux mis en oeuvre lors du premier test (41) dans le cas où les plots P1 et Préf sélectionnés ont un même récepteur pour origine. Selon une autre forme particulière de la variante de mise en oeuvre décrite précédemment, dans le cas où les plots P1 et P2 sélectionnés ont pour origine des récepteurs Rx1 et Rx2 distincts, le troisième test de compatibilité des couples (P1, Txi) et (P2, Tx2) est réalisé en mettant en oeuvre, sur les plots P1 et P2, des tests analogues à ceux mis en oeuvre lors du premier test dans le cas où les plots P1 et Préf sélectionnés ont pour origine des récepteurs distincts. Selon un mode particulier de mise en oeuvre du procédé selon l'invention, le second module de formation de n-uplets de couples (plot, émetteur) comporte les étapes suivantes:, - une première étape de formation de n-uplets par analyse des triplets validés formés à partir d'un même couple de référence (Préf, TXréf), l'analyse consistant à déterminer, en comparant tous les triplets formés, les couples (P1, Tx,) pour lesquels les propriétés suivantes sont vérifiées pour tous les triplets simultanément: - le plot P, est toujours associé à un même émetteur Txj, - l'émetteur Txi est toujours associé à un même plot Pi, les couples (plot, émetteur) vérifiant ces deux propriétés étant associés dans un même ensemble appelé n-uplet - une deuxième étape de filtrage, appliquée aux n-uplets formés lors de la première étape en fonction du nombre et de la nature des couples (P1, Tx,) contenus dans ces n-uplets; les n-uplets formés comportant moins de trois couples ou les n-uplets ne comprenant pas le couple de référence (Préf, TXréf) étant éliminés; - une troisième étape ayant pour objet, pour chaque n-uplet non éliminé par la deuxième étape, de vérifier la cohérence du n-uplet, la vérification de cohérence consistant à déterminer, en utilisant les mesures de distance bistatique, de vitesse bistatique et d'azimut réalisées sur les observations radar correspondant aux plots formant le n-uplet, si les plots formant les couples du n-uplet correspondent de manière vraisemblable à une même cible; - une quatrième étape consistant à mémoriser les n-uplets ayant passé les étapes précédentes avec succès; les plots formant ces n-uplets constituant des paquets ou groupes de plots destinés à l'initialisation de 25 pistes cartésiennes. Selon une variante de ce mode de mise en oeuvre, durant la troisième étape, on calcule d'abord la position cartésienne et la vitesse de la cible correspondant aux plots formant le n-uplet, puis on réalise la transposition 30 des paramètres cartésiens de la cible ainsi déterminée en paramètre bistatiques relatifs aux différentes bases bistatiques ayant produit les plots considérés, puis encore à comparer les paramètres bistatiques ainsi déterminés à ceux des plots considérés. L'invention a également pour objet un procédé d'initialisation de pistes cartésiennes à partir de plots bistatiques correspondant à une même cible détectée par un système radar multistatique comportant une pluralité d'émetteurs et de récepteurs distants les uns des autres, un couple (émetteur, récepteur) formant une base bistatique du système, chaque plot bistatique Pi correspondant à l'écho émis par un récepteur Txi et reçu, après réflexion sur une cible, par un récepteur Ek et étant caractérisé par l'émetteur Tx1 ayant émis le signal dont il est l'écho. Le procédé selon l'invention comporte: - une étape préliminaire mettant en oeuvre le procédé d'association de plots selon l'invention, cette étape permettant de déterminer un n-uplet de plots susceptible de correspondre à une même piste; - une étape d'initialisation de pistes proprement dite, consistant à déterminer, à partir des paramètres des plots formant le n-uplet considéré, 15 les paramètres de la cible correspondante. Selon un mode de mise en oeuvre particulier de ce procédé l'étape d'initialisation de piste comporte les sous-étapes suivantes: - une première sous-étape durant laquelle on réalise une première 20 estimation grossière de la position de la cible, cette première estimation étant réalisée par mise en oeuvre de l'algorithme du gradient; - une deuxième sous-étape durant laquelle on réalise un calcul d'intersection exhaustif permettant de déterminer, à partir des mesures des distances bistatiques des plots considérés et de la précision de ces 25 distances, ainsi que des positions du récepteur et des émetteurs considérés, la position de la cible ainsi qu'une estimation de la matrice de covariance associée à cette position. - d'une troisième sous-étape durant laquelle on calcule la vitesse de la cible; 30 - d'une quatrième sous-étape durant laquelle on estime la matrice de covariance initiale C de l'état de la cible à partir des bornes de Cramer Rao déterminées grâce aux précisions des mesures en distance bistatique, vitesse bistatique et azimut. - d'une cinquième sous-étape durant laquelle la nouvelle piste est 35 construite en associant touts les attributs qui la caractérisent. Selon un mode de mise en oeuvre particulier de ce procédé, celui-ci comporte en outre une étape complémentaire, préalable à l'étape préliminaire, consistant à réaliser une première association en deux 5 dimensions des plots bistatiques formés par le système, l'analyse de la compatibilité des plots considérés étant réalisée en ne tenant pas compte de l'altitude des observations correspondant à ces plots, ni des altitudes relatives des émetteurs et des récepteurs constituant le système; les plots bistatiques associés au cours de cette étape complémentaire étant 10 directement transmis à l'étape d'initialisation de pistes, l'étape préliminaire réalisant l'association en trois dimensions des seuls plots non associés par l'étape d'association en deux dimensions. L'invention a également pour objet un procédé d'initialisation de pistes 15 radar fusionnées, associant des pistes cartésiennes externes à des piste cartésiennes formées à partir des plots bistatiques produits par un système radar multistatique comportant une pluralité d'émetteurs et de récepteurs distants les uns des autres, un couple (émetteur, récepteur) formant une base bistatique du système, chaque plot bistatique Pi correspondant à l'écho 20 émis par un récepteur Txi et reçu par un récepteur Ek et étant caractérisé par l'émetteur Txj ayant émis le signal dont il est l'écho. Le procédé selon l'invention associe le procédé d'initialisation de pistes cartésiennes selon l'invention à un traitement d'association complémentaire comportant deux étapes: 25 - une première étape qui procède à l'association des pistes cartésiennes externes pour construire des pistes externes fusionnées dont les attributs sont fonctions des attributs de pistes cartésiennes externes à l'origine de la fusion; - une deuxième étape qui a pour fonction de tenter l'association de plots 30 radar bistatiques libres, non encore associés pour initialiser des pistes cartésiennes radar, à des pistes cartésiennes fusionnées externes formées au cours de l'étape précédente. Un tel procédé, quoique de portée générale, est avantageusement 35 adapté à l'initialisation de pistes dans le cadre du procédé de pistage qui fait l'objet de la demande de brevet français 0904144 citée précédemment. Les caractéristiques et avantages de l'invention seront mieux appréciés grâce à la description qui suit, description qui s'appuie sur les 5 figures annexées qui représentent: - la figure 1, le schéma synoptique d'un exemple de procédé de pistage adapté à un système radar passif multi-récepteurs; - La figure 2, le synoptique d'un premier exemple de module 10 d'initialisation de pistes cartésiennes mettant en oeuvre le procédé d'association 3D (association en trois dimensions) de plots bistatiques selon l'invention; - la figure 3, le synoptique du procédé d'association 3D de plots bistatiques selon l'invention; 15 - la figure 4, le synoptique du module de formation de triplets de couples (Plot, Emetteur) compatibles, du procédé d'association 3D de plots bistatiques selon l'invention; - la figure 5, le synoptique de principe du module de formation de nuplets du procédé d'association 3D de plots bistatiques selon l'invention; 20 - la figure 6, le synoptique d'un second exemple de module d'initialisation mettant en oeuvre le procédé d'association 3D de plots bistatiques selon l'invention; - la figure 7, un exemple de procédé d'initialisation qui réalise une initialisation de pistes cartésiennes radar en associant des plots bistatiques 25 et des pistes cartésiennes déjà formées. La figure 1 présente le schéma synoptique d'un exemple de procédé de pistage adapté à un système radar passif comportant plusieurs récepteurs, chaque récepteur formant, avec les émetteurs d'opportunités 30 présents sur la zone surveillée, différentes bases bistatiques. Un tel procédé comporte naturellement un ou plusieurs modules d'initialisation permettant de former des pistes à partir des détections élémentaires, c'est-à-dire des plots bistatiques formés. On rappelle ici qu'un plot bistatique correspond à une observation 35 radar ayant donné lieu à la détection d'une cible éventuelle. Cette observation est associée à des mesures spécifiques telles que notamment des mesures de distance bistatique, de vitesse bistatique et d'azimut. Comme l'illustre la figure 1, le procédé de pistage bistatique pris 5 comme exemple, comporte deux niveaux de fusion d'informations conduisant à la formation et à l'entretien de pistes cartésiennes. Le premier niveau 11 a pour fonction de réaliser une fusion des informations de détection, des plots bistatiques 13, produites par les bases bistatiques constituées à partir d'un même récepteur et de former et 10 d'entretenir des pistes cartésiennes élémentaires 14 à partir de ces plots bistatiques. Le second niveau 12 à quant à lui pour fonction de réaliser une fusion des informations de détection produites par l'ensemble des bases bistatiques, tous récepteurs considérés. Le traitement de fusion réalisé 15 consiste principalement à former et à entretenir des pistes cartésiennes globales 15 à partir des pistes cartésiennes élémentaires 14 supposées correspondre à un même objet vu par un même récepteur, ainsi que des plots bistatiques restés libres. En pratique, une piste cartésienne globale 15 est formée et entretenue 20 en utilisant les plots des pistes cartésiennes élémentaires 14 associées produites au premier niveau 11, ainsi que les plots bistatiques 16 non utilisés au premier niveau pour former des pistes cartésiennes élémentaires 14 (plots non pistés) et que les plots bistatiques 17 correspondant à des pistes cartésiennes élémentaires 14 non encore fusionnées (plots pistés libres). 25 Par ailleurs, comme on peut le constater sur la figure 1, chaque niveau de fusion, de formation et d'entretien de pistes cartésiennes, comporte un module d'initialisation 18 ou 19 qui procède à la constitution de pistes nouvelles, ce module mettant en oeuvre une initialisation adaptée. Les 30 figures 2 et 6 présentent à titre d'exemple les synoptiques de deux modules d'initialisations différents 21 et 61, pouvant correspondre aux modules d'initialisation 18 et 19, mettant en oeuvre le procédé 22 d'association 3D (association en trois dimensions) de plots bistatiques selon l'invention, avant de procéder à l'initialisation 23 proprement dite. La mise en oeuvre du 35 procédé selon l'invention permet avantageusement de limiter le nombre de plots bistatiques à considérer lors de l'initialisation. De manière générale, comme l'illustre la figure 2, ce procédé utilise des plots bistatiques libres 24 ou plots bistatiques candidats, c'est-à-dire des 5 plots non encore utilisés et candidats pour former de nouvelles pistes cartésiennes 25. Selon les cas, c'est à dire selon qu'on considère un système mono-récepteur (ou l'étape 18 de pistage mono-récepteur) ou multi-récepteurs (ou l'étape 19 de pistage multi-récepteur), les pistes cartésiennes formées sont 10 dites "pistes cartésiennes élémentaires" (les pistes 14 du pistage de la figure 1 par exemple) ou "pistes cartésiennes fusionnées" (les pistes 15 du pistage de la figure 1 par exemple). Le procédé d'association 3D de plots bistatiques selon l'invention est 15 illustré par les figures 3, 4 et 5. Comme l'illustre la figure 3, le procédé 22 d'association 3D de plots bistatiques selon l'invention met en oeuvre deux modules de traitement successifs, un premier module 31 de formation de triplets de couples (plot, 20 émetteur) considérés comme compatibles et un second module 32 de formation de n-uplets de couples (plot, émetteur). Les plots constituant les nuplets ainsi constitués sont ensuite délivrés au module 23 d'initialisation de pistes proprement dit, sous forme de groupes ou paquets de plots 35. Ces deux modules de traitement font eux-mêmes appel à deux 25 modules complémentaires, un premier module complémentaire 33 qui établit la liste des plots candidats disponibles et un second module complémentaire 34 qui assure la gestion et la formation de couples (plot, émetteur). Dans la suite du texte on décrit le fonctionnement de chacun des différents modules et l'interaction entre modules. 30 Il est à noter que par triplets de couples (Plot, Emetteur) "compatibles" on entend des triplets dont les plots correspondent vraisemblablement à une même cible potentielle. Le module 31 de formation de triplets compatible, ou, plus 35 simplement, module de formation de triplets, est illustré par la figure 4. Il consiste principalement à tenter l'association de trois couples entre eux, chaque couple étant constitué par un plot candidat et un émetteur. Le plot candidat est choisi par le module 34 dans la liste des plots libres définie par le module 33. La formation des couples est réalisée par le module de gestion des couples 34 qui attribue à chaque plot candidat un émetteur parmi ceux qui constituent le système. Cette attribution est faite de manière à ce que lorsque plusieurs plots parmi les trois plots candidats choisis proviennent d'un même récepteur, ils soient associés à des émetteurs différents. Selon l'invention on tente l'association de tous les triplets de couples 10 (plot, émetteur) possibles, sachant que des triplets sont associés lorsque, bien qu'ils soient formés par des bases bistatiques différentes, ils correspondent vraisemblablement à une même cible. A cet effet le module 31 de formation des triplets comporte cinq étapes de test distinctes décrites dans la suite du texte. Ces cinq étapes sont 15 réalisées de manière itérative sur l'ensemble des couples (plot, émetteur) possibles, (Pb Txj) par référence à un couple (Préf, TXréf) pris comme référence. La première étape de test 41 consiste, à déterminer s'il est possible 20 d'associer un premier couple (P1, Txi), au couple de référence (Préf, TXréf) , le couple candidat aussi bien que le couple de référence étant déterminé par le module 34 de gestion et de formation des couples. Selon l'invention, le couple (P1, Txi), est susceptible d'être associé au couple (Préf, TXréf) si ces deux couples correspondent vraisemblablement à une même cible. 25 Selon l'invention, cette étape de test peut être mise en oeuvre de différentes façons connues. Néanmoins dans une forme de mise en oeuvre préférée elle est mise en oeuvre comme décrit ci-après. Pour cette étape, le procédé selon l'invention distingue deux cas: 30 - le cas où les plots P1 et Préf sélectionnés par le module 34 ont un même récepteur Rx pour origine. - le cas où les plots P1 et Préf ont pour origine des récepteurs Rxi et RXréf distincts. 35 Dans le premier cas, l'étape 41 de test d'association consiste, en pratique, à mettre en oeuvre les quatre tests indépendants suivants: Un premier test (test 1) consiste à déterminer si les vitesses bistatiques Vbk(Pi) et Vbk(Préf) des plots P1 et Préf peuvent correspondre à 5 une même vitesse de cible. On considère à cet effet que si on a: Vbk (Préf ) - Vbk (P1)1> 2.Vmax , alors l'association du plot P1 au plot Préf n'est pas vraisemblable. Vmax 10 représente ici la vitesse maximale des cibles que l'on cherche à détecter grâce au système multistatique considéré. Un deuxième test (test 2) consiste à déterminer si les valeurs des azimuts associés aux plots P1 et Préf (i.e. les azimuts des observations radar 15 correspondantes) peuvent correspondre à une même de cible, compte tenu des précisions aek(P1) et dek(Préf) des mesures d'azimuts réalisées. On considère à cet effet que si on a: 0k (Préf)- Ok (P1)1> (Préf ) (520k (P1) 20 alors l'association du plot P1 au plot Préf n'est pas vraisemblable. Un troisième test (test 3) consiste à déterminer si les ellipsoïdes correspondant aux distances bistatiques des plots P1 et Préf, définies 25 respectivement par rapport à l'émetteur et au récepteur formant la base bistatique relative à P1 et à l'émetteur et au récepteur formant la base bistatique relative à Préf, présentent des intersections. On considère à cet effet que si on a: 30 (Dbk (Préf) - Dbk (P1))2 -IiiTxréf )(Tx1112 > o 5 alors l'association des plots est impossible. 5(Tx1 = [XTx1 YTx1 ZTx111- correspond ici à la position cartésienne de l'émetteur Tx1 associé à P1 et XTxréf = [XTxréf YTxréf ZTxréf f correspond à la position cartésienne de l'émetteur Txréf associé à Préf - Un Quatrième test (test 4) consiste à déterminer si les intersections des ellipsoïdes construites à partir des distances bistatiques caractérisant les plots P1 et Préf, sont compatibles des mesures d'azimuts (i.e. les azimuts des observations correspondantes). A cet effet on considère le site cp selon lequel la cible potentielle 10 correspondant aux plots considérés peut être observée. (I) est défini par la relation: Préf P1 cos = , Préf .e, .cos(Ok - 001)- COS(Ok - 80 2 ) 15 Dans cette relation, P1, e1, Préf, eréf sont définis par les relations suivantes (i= 1 ou réf selon le cas): = Dbk(1:1)-L R 2.Duk'ni 20 e = ' Db, avec L. = 5-(-nd inj De même, XRX et YRX représentant les coordonnées du récepteur 25 considéré, taneol et tano,réf ont pour expressions (i= 1 ou réf selon le cas): YT - YR tan 00 = x xTxi - XRx En outre, ek a pour expression: 30 01 En revanche, si le test échoue, les étapes 41 et 42 sont répétées. Le module 31 procède à un nouvel essai d'association du couple de référence (Préf, TXréf) (ou avec un autre couple de référence) avec d'autres couples candidats (P2, Tx2) et (P1, Txi). La quatrième étape de test 44 a pour objet de déterminer si les plots associés P1, P2 et Préf correspondent de manière vraisemblable à un même objet, une même cible. Pour ce faire on utilise les mesures de distance bistatique, de vitesse bistatique et d'azimut réalisées sur les observations radar correspondant aux plots P1, P2 et Préf. On calcule d'abord la position cartésienne x1,2,réf et la vitesse v1,2,réf de l'objet correspondant aux plots associés. Puis les résultats obtenus sont utilisés pour réaliser la projection des paramètres de cette cible dans les bases bistatiques ayant produit les plots bistatiques P1, P2 et Préf considérés et comparer les paramètres bistatiques ainsi déterminés à ceux des plots P1, P2 et Préf. On s'intéresse en premier lieu aux distances bistatiques. Par suite, on compare à chaque plots P1, P2 et Préf les projections ainsi obtenues. Pour chaque plot Pi l'inégalité suivante doit être vérifiée 1Dbk (P, )- Db1, 01(1 )15_ 3..\/6Db, (Pi ) + (Pi) Où Dbp,'j(P,) est la distance bistatique obtenue par projection de réf dans la base bistatique associée au plot PI, et ci-Dbi' (P,) l'écart type associé. Si cette inégalité n'est pas vérifiée, le triplet de couple (P1, Txi) (P2, Tx2) et (Préf, Txréf) est abandonné. Comme dans le cas de l'échec de la troisième étape de test, les étapes 41, 42 et 43 sont alors réitérées avec deux autres couples deux autres couples (P1, Txi) (P2, Tx2) et le cas échéant avec un autre couple de référence (Préf, TXréf). En revanche, dans le cas contraire, le triplet de couple (P1, Txi), (P2, Tx2) et (Préf, TXréf) est conservé et l'exécution du module 31 se poursuit alors par l'étape suivante 45. Selon l'invention, les paramètres cartésiens de la cible correspondant aux plots P1, P2 et Préf, peuvent bien entendu être déterminés par toute méthode connue à partir des paramètres bistatiques des plots considérés. De même, la projection (i. e. la transposition) des paramètres cartésiens de la cible ainsi déterminée en paramètre bistatiques relatifs aux différentes bases bistatiques à l'origine de P1, P2 et Préf, peut être réalisée par toute méthode de projection connue. Néanmoins il est à noter que la position cartésienne peut être déterminée, comme décrit précédemment, à partir des intersections des ellipsoïdes correspondant aux distances bistatiques des plots P1, P2 et Préf,. Par suite, la position cartésienne de la cible étant ainsi déterminée, le 15 calcul de la vitesse cartésienne peut être réalisé de manière simple par toute méthode connue. La cinquième étape de test 45 du module 31 a pour objet de vérifier que la cible correspondant aux plots associés P1, P2 et Pré correspond bien à 20 une cible réelle. Pour ce faire on compare la valeur absolue de la vitesse cartésienne déterminée à l'étape précédente à un seuil correspondant à la vitesse maximale à laquelle une cible réelle est supposée pouvoir évoluer. Par suite si le module de la vitesse cartésienne déterminée est supérieur à cette valeur maximale, le triplet de couple (P1, Txi) (P2, Tx2) et 25 (Préf, TXréf) est abandonné. Comme dans le cas de l'échec de la troisième ou de la quatrième étape de test, les étapes 41 à 44 sont alors réitérées avec deux autres couples deux autres couples (Pl, Tx1) (P2, Tx2) et le cas échéant avec un autre couple de référence (Préf, TXréf). Dans le cas contraire le triplet de couple considéré est finalement 30 conservé. D'un point de vue fonctionnel global, le module 31 a pour objet de déterminer, en considérant un couple (plot, émetteur) donné pris comme référence (Préf, TXréf), l'ensemble des couples (Ph Txj) associables en triplets, 35 c'est-à-dire susceptibles de correspondre à une même cible. Par suite son exécution est réitérée en considérant successivement toutes les associations possibles des couples formant la liste des couples candidats avec le couple de référence (Préf, Txréf)- Finalement, l'itération du module 31 cesse lorsque toutes les associations possibles en triplets ont été testées. Par ailleurs chacune des itérations réalisées se termine, en cas de succès des tests d'association, par une étape 46 durant laquelle on procède le cas échéant à la mémorisation du triplet (P1, Txi) (P2, Tx2) et (Préf, Txréf) constitué pour cette itération. Il est à noter que la détermination à chaque instant des couples (P1, Txi) (P2, Tx2) et (Préf, Txréf) dont on doit déterminer la possibilité d'association, est réalisée par le module de gestion 34. Pour réaliser cette tâche le module 34 reçoit de la part du module 31 des informations relatives à l'échec des tests réalisé durant les étapes de test 41 à 45. La transmission de ces informations est symbolisée par les flèches entraits discontinus sur la figure 4. Ainsi, si le test réalisé par l'étape 41 ou l'étape 42 échoue, le couple (P1, Txi) ou (P2, Tx2) considéré est remplacé par un autre couple candidat et le module 31 est de nouveau exécuté. En outre, si tous les couples candidats ont été testés sans succès avec le couple de référence considéré celui-ci est remplacé, par le module de gestion 34, par un autre couple de référence. De même, si le test réalisé par l'étape 43, 44 ou 45 échoue, les couples (P1, Txl) et (P2, Tx2) considérés sont remplacés, par le module de gestion 34, par deux autres couples candidats et le module 31 est de nouveau exécuté. En outre, si tous les couples candidats ont été testés sans succès avec le couple de référence celui-ci est remplacé par un autre couple de référence. Par ailleurs, il est également à noter que, lorsque l'association de couples candidats (P1, Txi) et (P2, Tx2) et du couple de référence (Préf, Txréf) considéré est réalisée avec succès, le module 31 transmet au module 34, lors de l'étape 46, une information de fin d'exécution qui renseigne ce module sur le fait qu'une nouvelle association du couple de référence (Préf, Txréf) avec d'autres couples candidats (Pi, Txi) et (P2, Tx2) peut être tentée. En outre, si tous les couples candidats ont été testés, avec ou sans succès, avec le couple de référence, le processus de formation de triplets à partir du couple de référence (Préf, TXréf) est terminé et l'exécution du module 31 cède la place à l'exécution du module 32 de formation de n-uplets. Ainsi, le module de gestion 34 gère l'exécution itérative du module 31 5 de formation de triplets, cette exécution étant stoppée lorsque, pour le couple de référence (Préf, Txréf) considéré, l'ensemble des triplets de couples candidats ont été considérés. Les triplets validés, mémorisés durant l'étape 46 du module 31, sont 10 utilisés par le module 32 de formation de n-uplets. Le synoptique de principe du module 32, présenté sur la figure 5 montre que ce module comporte quatre étapes de traitement. La première étape 51 a pour objet d'analyser les triplets validés 15 disponibles. Ces triplets construits par implémentation du module 31 comportent un même couple de référence (Préf, TXréf). L'analyse consiste ici à déterminer, en comparant tous les triplets formés, les couples (P,, Txj) pour lesquels les propriétés suivantes sont vérifiées pour tous les triplets simultanément: 20 - le plot P; est toujours associé à un même émetteur Txj, - l'émetteur Txi est toujours associé à un même plot P1. Par suite, les couples (plot, émetteur) vérifiant ces deux propriétés sont associés dans un même ensemble appelé n-uplet, n représentant le nombre de couples (P1, Txi) formant l'ensemble considéré. 25 Ainsi, par exemple, si pour un couple de référence donné (P1, Tx1) on dispose après exécution du module 31 de quatre triplets décrits par le tableau 1 suivant: Triplet n°1 Triplet n°2 Triplet n°3 Triplet n°4 Tx1 P1 Tx1 P1 Tx1 P1 Tx1 P1 Tx2 P2 Tx2 P3 Tx4 P6 Tx4 P6 Tx3 P3 Tx4 P6 Tx5 P5 Tx7 P4 30 le n-uplet formé à partir de ces triplets conformément au procédé selon l'invention, un quadruplet dans cet exemple, sera décrit par le tableau 2 suivant: n-u plet Tx1 Pi Tx4 P6 Tx5 P5 TX7 P4 Les plots P2 et P3, sont ici écarté du quadruplet du fait que TX2 est 5 associé à deux plots différents, à savoir P2 et P3, et que P3 est associé à deux émetteurs différents, à savoir TX2 et Tx3. Ainsi, pour le couple de référence (Préf, TXréf) considéré, l'étape 51 permet de former de n-uplets de couples (Ph Txj) pour lesquels chaque plot 10 est associé bijectivement à un émetteur. La deuxième étape 52 a pour objet d'analyser de filtrer les n-uplets formés lors de l'étape 51 en fonction du nombre et de la nature des couples (Ph Txj) contenus dans ces n-uplets. 15 Ainsi, selon l'invention, si le n-uplet formé comporte moins de trois couples ou encore si ce n-uplet ne comprend pas le couple de référence (Préf, Txréf), à partir duquel les triplets ont été formés, le n-uplet formé est éliminé. Par suite, le module 32 transmet au module 34 de gestion des couples (Pb Tx,) une information conduisant le module 34 à choisir un autre 20 couple (P1, Txj) comme couple de référence et à relancer l'exécution du module 31. En revanche, si le n-uplet formé comporte trois couples ou plus, et si ce n-uplet comprend le couple de référence (Prêt, TXréf) il est conservé et utilisé lors de l'étape suivante. 25 La troisième étape 53 du module 32 a pour objet de vérifier la cohérence du n-uplet de couples formés à l'étape 51 et ayant passé l'étape 52. Cette vérification de cohérence consiste principalement à déterminer si les plots formant les couples du n-uplet correspondent de manière vraisemblable à un même objet. 30 Pour ce faire, comme dans le cas de l'étape de test 44 du module 31 de formation de triplets, on utilise les mesures de distance bistatique, de vitesse bistatique et d'azimut réalisées sur les observations radar correspondant aux plots formant le n-uplet. On calcule d'abord la position cartésienne et la vitesse de l'objet correspondant aux plots considérés. Puis les résultats obtenus sont utilisés pour réaliser la projection (i. e. la transposition) des paramètres cartésiens de la cible ainsi déterminée en paramètre bistatiques relatifs aux différentes bases bistatiques ayant produit les plots considérés, et comparer les paramètres bistatiques ainsi déterminés à ceux de ces plots. Dans le cas où la vérification est positive, le n-uplet est conservé. Dans le cas contraire, il est rejeté et le module 32 transmet au module 34 de gestion des couples (P1, Tx,) une information conduisant ce module à choisir un autre couple (P,, Txj) comme couple de référence et à relancer l'exécution du module 31 et du module 32. La quatrième étape 54 consiste quant à elle à mémoriser les n-uplets ayant passé les étapes 52 et 53 avec succès, les plots formant ces n-uplets, destinés à être utilisés pour former des pistes cartésiennes nouvelles, constituent des paquets ou groupes de plots qui sont transmis au module 23 d'initialisation de pistes proprement dit. Elle consiste également à transmettre une information au module 34 de gestion, conduisant ce module à choisir un autre couple (Ph Txj) comme couple de référence et à relancer l'exécution du module 31et du module 32. Elle consiste encore, au moment ou un n-uplet est mémorisé, à transmettre au module 33 d'actualisation de la liste des plots disponibles, la liste des plots constituant le n-uplet. Le module 33 actualise la liste des plots disponibles en conséquence et communique la liste actualisée au module de gestion 34 qui gère les couples (P1, Txj) candidats, auxquels le module 31 de formation de triplets est appliqué. D'un point de vue fonctionnel global le module 32 de formation de nuplet est exécuté pour chaque couple de référence (Préf, TXréf) à partir des triplets de couples (plot, émetteur) constitués par exécution itérative du module 31. Par suite le procédé d'association 3D de plots selon l'invention fonctionne à partir de deux boucles d'itération imbriquées, gérées à partir du module de gestion 34 et du module 33 d'actualisation. La première boucle, sur l'ensemble des couples candidats (P1, Txj), illustrée par la flèche 36 en trait discontinu sur la figure 3, concerne l'exécution itérative du module 31. La seconde boucle effectuée sur l'ensemble des couples de référence (Préf, TXréf), illustrée par la flèche 37 en trait discontinu sur la figure 3, concerne l'exécution itérative de l'ensemble constitué par les modules 31 et 32. Ainsi, la mise en oeuvre du procédé selon l'invention consiste pour un couple (Préf, TXréf) à exécuter une séquence consistant: - à déterminer, par itérations successives du module 31 sur les 10 couples (Pb Txj), l'ensemble des triplets compatibles déterminés à partir du couple (Préf, TXréf) - à former à partir de ces triplets, le n-uplet correspondant au couple (Préf, TXréf), par exécution du module 32; cette séquence étant répétée pour chaque couple (plot, émetteur) défini 15 comme couple de référence par le module 34 à partir de la liste de plots disponibles élaborée par le module 33. Les n-uplets de plots ainsi obtenus par mise en oeuvre du procédé d'association 3D de plots selon l'invention sont mis en mémoire sous forme 20 d'un ensemble de plots validés qui peuvent par exemple être utilisé, comme décrit dans la suite du document, dans un procédé global d'initialisation de pistes cartésiennes à partir de plots bistatiques dans le cadre, par exemple, du pistage réalisé par un radar multistatique comportant un ou plusieurs récepteurs formant une pluralité de bases bistatiques. 25 Comme l'illustrent les figures 2, 6 et 7, le procédé selon l'invention peut être mis en oeuvre en tant qu'étape d'un procédé d'initialisation de pistes, préliminaire à une étape d'initialisation 23 proprement dite. Les figures 2 et 6 présentent plus particulièrement des exemples de procédés 30 d'initialisation de pistes cartésiennes 25, qui ne prennent en considération que des plots bistatiques 24, procédés applicables notamment pour réaliser l'initialisation de pistes cartésiennes élémentaire 18 du procédé de pistage présenté à la figure 1. La figure 7, quant à elle, présente plus particulièrement un exemple de procédé d'initialisation qui réalise une 35 initialisation de pistes cartésiennes radar en associant des plots bistatiques et des pistes cartésiennes déjà formées, provenant par exemple d'autres moyens techniques associés au système radar multistatique considéré. Dans le procédé 21 d'initialisation illustré par la figure 2, les plots 5 bistatiques 24 produits à un instant donné, plots qui correspondent à des observations réalisées par le système radar et considérées comme relatives à des cibles potentielles, sont directement exploités par l'étape préliminaire 22 mettant en oeuvre le procédé selon l'invention. Par suite les couples (plot, Emetteur) retenus à l'issue de l'étape 22 significatifs, correspondant 10 vraisemblablement à des cibles, sont utilisés lors de l'étape 23 d'initialisation de pistes proprement dite. L'étape d'initialisation 23 peut être réalisée de différentes façons connues, dans la mesure où elle consiste à déterminer, à partir des paramètres des plots formant les n-uplets constitués à l'issue de l'étape 22, 15 les paramètres des cibles correspondantes. Cependant dans une forme préférée de mise en oeuvre, l'étape d'initialisation est réalisée en différentes sous-étapes: - une première sous-étape durant laquelle on réalise une première estimation grossière de la position de la cible. Cette première estimation est 20 réalisée par mise en oeuvre de l'algorithme du gradient ou par toute méthode appropriée connue. Ces méthodes nécessitent trois plots ou plus lorsque seule l'information de distance bistatique est utilisée. - une deuxième étape durant laquelle on réalise un calcul d'intersection exhaustif qui permet de déterminer de manière connue la 25 position de la cible, ainsi qu'une estimation de la matrice de covariance associée à cette position, à partir des mesures des distances bistatiques des plots considérés et de la précision de ces distances, ainsi que des positions du récepteur et des émetteurs considérés. Une méthode permettant de réaliser un calcul exhaustif est notamment décrite dans la demande de 30 brevet français portant le numéro FR2910131. Par suite une piste est initialisée si: - la position de la cible peut être déterminée - les écarts-types estimés des coordonnées cartésiennes(x, y, z) de la cible ne dépassent pas une valeur donnée. 35 Dans ce cas la deuxième étape est suivie: - d'une troisième sous-étape durant laquelle on calcule la vitesse de la cible; - d'une quatrième sous-étape durant laquelle on estime la matrice de covariance initiale C de l'état de la cible à partir des bornes de Cramer Rao déterminées grâce aux précisions des mesures en distance bistatique, vitesse bistatique et azimut. Pour N plots P1, P2, - - -7 PN les bornes de Cramer Rao s'écrivent de la manière suivante: (62 (P ) 0 Dbk 1 O 152 (p \ Vbk 1 / O O O O a ( ) Dbk O 0V \ bk (p N C = P Avec .p T P= raDbk(11) ax aDbk(11) av, aDbk(171) ay aDbk(Pi) av, aDbk(P1) aDbakz (P1) avz avbk(p1) ax avbk(p1) avy avbk(pi) ay avbk(pi) avy avbk(p1) avbakz (pi) avz aDbk(PN) ax aDbk(P) avx aDbk (PN ) ay aDbk(pN) av, aDbk(pN) a aDbkz(EN) avz avbk(pN ax avbk(PN) av' avbk(PN) ay avbk(PN) av, avb'(pN) avbk(pN) avz - d'une cinquième sous-étape durant laquelle la nouvelle piste est construite en associant tous les attributs qui la caractérisent. Dans le procédé d'initialisation 61 illustré par la figure 6, plus complexe que celui présenté sur la figure 2, l'étape 22 mettant en oeuvre le 20 procédé d'association 3D de plots selon l'invention est associée à une étape complémentaire 62 mettant en oeuvre un procédé d'association 2D de plots, procédé simplifié consistant à tenter d'associer les plots bistatiques 24 produits par le système radar, l'analyse de la compatibilité des plots étant dans ce cas réalisée en ne tenant pas compte de l'altitude des observations correspondant aux plots considérés, ni des altitudes relatives des émetteurs et des récepteurs constituant le système radar. Il permet généralement d'effectuer de manière simple une première association de plots qui sont directement transmis à l'étape 23 d'initialisation proprement dite. Par suite, l'étape 22 d'association 3D de plots n'est appliquée qu'aux seuls plots n'ayant pu faire l'objet d'une association lors de l'étape 62 d'association 2D. Comparée à la forme précédente, cette forme de mise en oeuvre du procédé d'initialisation, présente l'avantage de permettre de traiter dans un même temps un plus grand nombre de plots bistatiques, l'étape 62 étant d'une exécution plus rapide que l'étape 22. Le procédé de l'étape 62 d'association 2D, consiste globalement à associer a priori des plots bistatiques et à tester la validité de cette association en déterminant si une cible réelle peut correspondre à cette association, sachant que toutes les associations possibles sont étudiées. La détermination est principalement réalisée en recherchant les intersections des projections dans un même plan horizontal des ellipsoïdes définies à partir des distances bistatiques des plots associés. Différentes méthodes connues, non décrites ici, permettent de réaliser une telle association 2D. Cependant, dans une forme de mise en oeuvre préférée, l'association 2D est réalisée en mettant en oeuvre un module d'association comportant 25 deux étapes La première étape associe deux plots formés à partir d'un même récepteur; ces deux plots formant une base à laquelle sont ensuite associés successivement les plots formés à partir de tous les récepteurs. Durant cette première étape la validité de l'association des deux premiers plots formés à 30 partir d'un même récepteur est d'abord testée. Si cette validité est vérifiée, on teste la validité de l'association de ces deux plots avec les plots formés à partir de tous les récepteurs en les choisissant l'un après l'autre. La première étape est réalisée pour toutes les associations possibles. Par suite les nuplets de plots constituant des associations validées sont considérés comme 35 correspondant vraisemblablement à une même cible réelle et mémorisés. 2 9 79 709 29 La seconde étape quant à elle associe d'abord deux plots formés à partir de récepteurs différents et teste la validité de cette association, puis associe ces deux plots avec tous les autres plots issus de tous les récepteurs l'un après l'autre. 5 Par suite, comme pour la première étape, les n-uplets de plots constituant des associations validées sont considérés comme correspondant vraisemblablement à une même cible réelle et mémorisés. Cependant les triplets de plots testés sont ici ceux qui ne font partie d'aucun des n-uplets retenus lors de la première étape. 10 Dans le cas où il n'y a qu'un seul récepteur, seule la première étape est réalisée. Le procédé 71 d'initialisation de pistes radar fusionnées, illustré par la figure 7, intègre le procédé d'initialisation 61 illustré par la figure 6 dans une 15 structure fonctionnelle plus large. Le procédé présenté a pour objet de réaliser l'initialisation de nouvelles pistes cartésiennes à partir de plots bistatiques et de pistes cartésiennes déjà formées. En particulier, dans le cas où l'on dispose de pistes cartésiennes délivrées par des moyens annexes associés au système radar multistatique considéré, le procédé illustré par la 20 figure 7 a pour objet, dans un premier temps, de tenter de fusionner les pistes cartésiennes externes produite à partir de moyens divers associés au radar, puis d'associer à ces pistes cartésiennes fusionnées, dites "pistes cartésiennes fusionnées externes" des plots bistatiques radar nouvellement formés et non encore utilisés pour former des pistes cartésiennes radar. Le procédé 71 présenté ici associe au procédé d'association selon l'invention dans sa forme simple 21 (association 3D, 22) illustrée par la figure 2 ou dans sa forme plus élaborée 61 (Associations 3D, 22, et 2D, 62) illustrée par la figure 6, à un traitement d'association complémentaire comportant deux étapes 72 et 73: La première étape 72 procède à l'association des pistes cartésiennes externes 74, l'association de ces pistes cartésiennes externes entre elles pouvant être réalisée, selon l'invention, par toute méthode connue, appropriée à la fusion de pistes. La fonction de cette étape d'association est d'abord de déterminer si deux pistes cartésiennes externes sont compatibles c'est-à-dire si elles peuvent être considérées comme représentant un même objet. C'est ensuite, si tel est le cas, de construire une piste fusionnée dont les attributs sont fonctions des attributs des deux pistes cartésiennes externes reconnues comme compatibles. L'étape 72 délivre ainsi des pistes cartésiennes fusionnées dites "pistes cartésiennes fusionnées externes". La deuxième étape 73 a pour fonction de tenter l'association de plots radar bistatiques libres, non encore associés pour initialiser des pistes cartésiennes radar, à des pistes cartésiennes fusionnées externes formées au cours de l'étape précédente 72. Selon l'invention, l'association est réalisée par toute méthode d'association plot-piste connue qui permet d'entretenir une piste à l'aide de plots, un plot correspondant à la détection et à la caractérisation d'un objet donné. Elle comporte notamment la mise en oeuvre d'un test d'association plot-piste. Par suite, si une piste fusionnée externe est susceptible d'être associée à un plot bistatique, l'association est réalisée et la piste intègre ce plot. La piste ainsi formée dite "pistes cartésiennes externes fusionnées assimilées" est considérée comme une piste radar initialisée. En revanche, si une piste cartésienne fusionnée externe qui n'est susceptible d'être associée à aucun plot bistatique, est rejetée. Selon ce procédé, les plots bistatiques libres non associés sont 20 transmis au procédé d'association 3D de plots selon l'invention 21 ou 61 pour initialiser des pistes radar. Suivant cette application, comme l'illustre la figure 7, les pistes radar ainsi initialisées à l'issue par mise en oeuvre du procédé 21 ou 61 selon l'invention, ainsi que les pistes cartésiennes externes fusionnées assimilées, 25 constituent de nouvelles pistes cartésiennes radar 75, destinées à être exploitées par des moyens de pistage du radar proprement dits. REVENDICATIONS1. Procédé (22) pour réaliser l'association de plots bistatiques correspondant à une même cible détectée par un système radar multistatique comportant une pluralité d'émetteurs et de récepteurs distants les uns des autres, un couple (émetteur, récepteur) formant une base bistatique du système, chaque plot bistatique Pi correspondant à l'écho provenant du signal émis par un émetteur Tx1 et reçu par un récepteur Ek et étant caractérisé par l'émetteur Tx1 ayant émis le signal dont il est l'écho, caractérisé en ce que, les plots bistatiques formés étant associés aux différents émetteurs du système pour former un ensemble de couples (plot, émetteur), il comporte: - un premier module (31) de formation de triplets de couples (Pi, Txj), ces triplets, regroupant des couples (plot, émetteur) qui présentent une compatibilité de distance et de vitesse; - un second module (32) de formation de n-uplets de couples (Pi, Txj), ces n-uplets regroupant des triplets formés à l'étape précédente répondant à des critères de regroupement donnés qui caractérisent la vraisemblance de l'attribution des différents triplets à une même cible; Les n-uplets ainsi constitués formant des paquets de plots (35) qui définissent une même cible. 2. Procédé selon la revendication 1, caractérisé en ce que la première étape (31) de formation de triplets de couples (plot, émetteur) est une étape itérative qui considère à chaque itération trois couples (plot, émetteur), un premier couple (P1, Txi), un second couple (P2, Tx2) et un troisième couple (Préf, TXréf) pris comme un couple de référence et qui leur applique une opération d'association qui teste la possibilité d'associer le premier couple (P1, Txi) et le second couple (P2, Tx2) au couple de référence (Préf, TXréf), un couple (plot, émetteur) étant associable au couple de référence si les caractéristiques de la cible correspondant au plot considéré sont compatibles des caractéristiques de la cible correspondant au plot formant le couple de référence; trois couples (Prêt, TXréf), (P1, Tx1) et (P2, Tx2), pour lesquels le test d'association est réussi formant un triplet de couples (47). 3. .Procédé selon la revendication 2, caractérisé en ce que la première étape de formation de triplets de couples (plot, émetteur) comporte elle-même les opérations suivantes: - un premier test (41) d'association pour déterminer si les couples (Préf, TXréf) et (P1, Txi) sont associables; - un deuxième test (42) d'association pour déterminer si les couples (Préf, TXréf) et (P2, Tx2) sont associables; - un troisième test (43) de compatibilité des couples (P1, Txi) et (P2, Tx2) - un quatrième test (44) de compatibilités des distances bistatiques des plots - un cinquième test (45) de test de la vraisemblance de la vitesse cartésienne de la cible correspondant aux plots associés; - une sixième opération (46) de mémorisation des triplets (47) formés de couples (P1, Txi) et (P2, Tx2) ayant réussi tous les tests; le test d'association mis en oeuvre lors du premier test (41) étant identique à celui mis en oeuvre lors du second test (42). 4. Procédé selon la revendication 3, caractérisé en ce que, dans le cas où les plots P1 et Préf sélectionnés par le module 34 ont un même récepteur Rx pour origine, le test d'association mis en oeuvre lors du premier (41) test consiste à mettre en oeuvre les quatre tests indépendants suivants: - un premier test (test 1) consistant à déterminer si les vitesses bistatiques Vbk(Pi) et Vbk(Préf) des plots P1 et Prêt peuvent correspondre à une même vitesse de cible, le premier test étant négatif si l'on a: Vbx (Pref ) - Vbk (Pi ) > 2.Vmax Où Vmax représente la vitesse maximale des cibles que l'on cherche à détecter; - Un deuxième test (test 2) consistant à déterminer si les valeurs des azimuts associés aux plots P1 et Préf peuvent correspondre à une même cible, compte tenu des précisions aek(Pi) et aek(Préf) des mesures d'azimuts réalisées; le deuxième test étant négatif si l'on a: lek (Préf) ) - Ok (P1)1 > 3.iiGe2k (Préf) + 620, (Pl)- Un troisième test (test 3) consistant à déterminer si les ellipsoïdes correspondant aux distances bistatiques des plots P1 et Préf présentent des intersections, le troisième test étant négatif si l'on a: (Dbk (Pré, ) - Dbk (P1 ))2 - XTxréf i.(1-x1112 > où *> Txl [ _XTx1 YTx1 ZTx1 f et XTxréf = [XTxréf YTxréf ZTxréf correspondent ici respectivement à la position cartésienne de l'émetteur Tx1 associé à P1 et à la position cartésienne de l'émetteur Txréf associé à Préf; Un Quatrième test (test 4) consistant à déterminer si les intersections des ellipsoïdes construites à partir des distances bistatiques caractérisant les plots P1 et Préf, sont compatibles des mesures d'azimuts relatives à ces plots le quatrième test étant négatif si l'on a: [cos ç - 3.6's(, cos ç 3.Gcospin où ç correspond au site selon lequel la cible potentielle correspondant aux plots considérés peut être observée et acosq, l'incertitude affectant la détermination de cosy); le test d'association étant positif si les résultats des quatre tests sont positifs. 5. Procédé selon la revendication 4, caractérisé en ce que, dans le cas où les plots P2 et Préf sélectionnés par le module 34 ont un même récepteur Rx pour origine, le test d'association mis en oeuvre lors du second (42) test consiste à mettre en oeuvre, sur les plots P2 et Préf, les quatre tests indépendants mis en oeuvre lors du premier test. 6. Procédé selon la revendication 3, caractérisé en ce que, dans le cas où les plots P1 et Préf sélectionnés par le module 34 ont pour origine des récepteurs Rx1 et Rxréf distincts, le test d'association mis en oeuvre lors du premier (41) test consiste à mettre en oeuvre consécutivement les quatre tests suivants: - un premier test (test A) consistant à déterminer s'il est possible de définir, dans la limite de portée du système de détection considéré, une intersection des directions des azimuts des observations correspondant aux plots P1 et Préf; et, dans l'affirmative, à calculer cette intersection, en prenant en compte les incertitudes des mesures des azimuts de sorte que l'on définit une aire d'intersection; [0,1]# 0 ,- un deuxième test (test B) consistant, si le test A est positif, à déterminer si les deux sphères contenant les ellipsoïdes définis par les distances bistatiques des plots P1 et Préf présentent une intersection; - un troisième test (test C) consistant, si le test B est positif, à déterminer si les deux ellipsoïdes définis par les distances bistatiques des plots P1 et Préf présentent une intersection. - un quatrième test (test D) consistant, si le test C est positif, à déterminer si les mesures des distances bistatiques associées aux plots P1 et Préf sont compatible de la zone d'intersection des directions des azimuts des observations correspondant aux plots P1 et Préf déterminée lors du test A; le test d'association étant positif si les résultats des quatre tests sont positifs. 7. Procédé selon la revendication 6, caractérisé en ce que, dans le cas où les plots P2 et Préf sélectionnés par le module 34 ont pour origine des récepteurs Rx2 et Rxréf distincts, le test d'association mis en oeuvre lors du second test (42) consiste à mettre en oeuvre consécutivement, sur les plots P2 et Préf, les quatre tests mis en oeuvre lors du premier test (41). 8. Procédé selon l'une des revendications 5 ou 6, caractérisé en ce que le troisième test (test C) est réalisé par la méthode du multiplicateur de Lagrange appliquée aux deux ellipsoïdes considérées. 9. Procédé selon l'une quelconque des revendications 3 à 8, caractérisé en ce que, dans le cas où les plots P1 et P2 sélectionnés par le module 34 ont un même récepteur Rx pour origine, le troisième test (43) de compatibilité des couples (P1, Txi) et (P2, Tx2) est réalisé en mettant en oeuvre, sur les plots P1 et P2, des tests analogues à ceux mis en oeuvre lors du premier test (41) dans le cas où les plots P1 et Préf sélectionnés par le module 34 ont un même récepteur pour origine. 10. Procédé selon l'une quelconque des revendications 3 à 8, caractérisé en ce que, dans le cas où les plots P1 et P2 sélectionnés par le module 34 ont pour origine des récepteurs Rx1 et Rx2 distincts, le troisièmetest (43) de compatibilité des couples (P1, Txi) et (P2, Tx2) est réalisé en mettant en oeuvre, sur les plots P1 et P2, des tests analogues à ceux mis en oeuvre lors du premier test (41) dans le cas où les plots P1 et Préf sélectionnés par le module 34 ont pour origine des récepteurs distincts. 11. Procédé selon l'une des revendications 1 à 10, caractérisé en ce que le second module (32) de formation de n-uplets de couples (plot, émetteur) comporte les étapes suivantes:, - une première étape (51) de formation de n-uplets par analyse des triplets validés (47) formés à partir d'un même couple de référence (Préf, TXréf), l'analyse consistant à déterminer, en comparant tous les triplets formés, les couples (P' Tx1) pour lesquels les propriétés suivantes sont vérifiées pour tous les triplets simultanément: - le plot Pi est toujours associé à un même émetteur TXJ, - l'émetteur Txj est toujours associé à un même plot Pi; les couples (plot, émetteur) vérifiant ces deux propriétés étant associés dans un même ensemble appelé n-uplet - une deuxième étape (52) de filtrage, appliquée aux n-uplets formés lors de la première étape (51) en fonction du nombre et de la nature des couples (P1, Txj) contenus dans ces n-uplets; les n-uplets formés comportant moins de trois couples ou les n-uplets ne comprenant pas le couple de référence (Préf, TXréf) étant éliminés; - une troisième étape (53) ayant pour objet, pour chaque n-uplet non éliminé par la deuxième étape (52), de vérifier la cohérence du n-uplet, la vérification de cohérence consistant à déterminer, en utilisant les mesures de distance bistatique, de vitesse bistatique et d'azimut réalisées sur les observations radar correspondant aux plots formant le n-uplet, si les plots formant les couples du n-uplet correspondent de manière vraisemblable à une même cible; - une quatrième étape (54) consistant à mémoriser les n-uplets ayant passé les étapes précédentes avec succès; les plots formant ces nuplets constituant des paquets ou groupes de plots (35) destiné à l'initialisation de pistes cartésiennes. 12. Procédé selon la revendication 11, caractérisé en ce que, durant la troisième étape (53), on calcule d'abord la position cartésienne et la vitesse de la cible correspondant aux plots formant le n-uplet, puis on réalise la transposition des paramètres cartésiens de la cible ainsi déterminée en paramètre bistatiques relatifs aux différentes bases bistatiques ayant produit les plots considérés, puis encore à comparer les paramètres bistatiques ainsi déterminés à ceux des plots considérés. 13. Procédé d'initialisation de pistes cartésiennes (21) à partir de plots bistatiques correspondant à une même cible détectée par un système radar multistatique comportant une pluralité d'émetteurs et de récepteurs distants les uns des autres, un couple (émetteur, récepteur) formant une base bistatique du système, chaque plot bistatique Pi correspondant à l'écho émis par un récepteur Txj et reçu, après réflexion sur une cible, par un récepteur Ek et étant caractérisé par l'émetteur Txj ayant émis le signal dont il est l'écho; caractérisé en ce qu'il comporte: - une étape préliminaire (22) mettant en oeuvre le procédé d'association de plots selon l'une quelconque des revendications 1 à 12 cette étape permettant de déterminer un n-uplet de plots susceptible de correspondre à une même piste; - une étape (23) d'initialisation de pistes proprement dite, consistant à déterminer, à partir des paramètres des plots formant le n-uplet considéré, les paramètres de la cible correspondante. 14. procédé selon la revendication 13, caractérisé en ce que l'étape d'initialisation de piste (23) comporte les sous-étapes suivantes: - une première sous-étape durant laquelle on réalise une première estimation grossière de la position de la cible, cette première estimation étant réalisée par mise en oeuvre de l'algorithme du gradient; - une deuxième sous-étape durant laquelle on réalise un calcul d'intersection exhaustif permettant de déterminer, à partir des mesures des distances bistatiques des plots considérés et de la précision de ces distances, ainsi que des positions du récepteur et des émetteurs considérés, la position de la cible ainsi qu'une estimation de la matrice de covariance associée à cette position. 2 9 79 709 37 - d'une troisième sous-étape durant laquelle on calcule la vitesse de la cible; - d'une quatrième sous-étape durant laquelle on estime la matrice de covariance initiale C de l'état de la cible à partir des bornes de Cramer 5 Rao déterminées grâce aux précisions des mesures en distance bistatique, vitesse bistatique et azimut. - d'une cinquième sous-étape durant laquelle la nouvelle piste est construite en associant touts les attributs qui la caractérisent. 10 15. Procédé d'initialisation de pistes cartésiennes (61) selon l'une quelconque des revendications 13 ou 14, caractérisé en ce qu'il comporte en outre une étape complémentaire (62), préalable à l'étape préliminaire (22), consistant à réaliser une première association en deux dimensions des plots bistatiques formés par le système, l'analyse de la 15 compatibilité des plots considérés étant réalisée en ne tenant pas compte de l'altitude des observations correspondant à ces plots, ni des altitudes relatives des émetteurs et des récepteurs constituant le système; les plots bistatiques associés au cours de cette étape complémentaire (62) étant directement transmis à l'étape (23) 20 d'initialisation de pistes, l'étape préliminaire (22) réalisant l'association en trois dimensions des seuls plots non associés par l'étape (62) d'association en deux dimensions. 16. Procédé d'initialisation de pistes radar fusionnées (71), associant des 25 pistes cartésiennes externes à des piste cartésiennes formées à partir des plots bistatiques produits par un système radar multistatique comportant une pluralité d'émetteurs et de récepteurs distants les uns des autres, un couple (émetteur, récepteur) formant une base bistatique du système, chaque plot bistatique Pi correspondant à l'écho émis par un 30 récepteur Txi et reçu par un récepteur Ek et étant caractérisé par l'émetteur Txi ayant émis le signal dont il est l'écho; caractérisé en ce qu'il associe le procédé d'initialisation (21, 61) de pistes cartésiennes selon l'une quelconque des revendications 13 à 15 à un traitement d'association complémentaire comportant deux étapes:- une première étape (72) qui procède à l'association des pistes cartésiennes externes (74) pour construire des pistes externes fusionnées dont les attributs sont fonctions des attributs de pistes cartésiennes externes à l'origine de la fusion; - une deuxième étape (73) qui a pour fonction de tenter l'association de plots radar bistatiques libres, non encore associés pour initialiser des pistes cartésiennes radar, à des pistes cartésiennes fusionnées externes formées au cours de l'étape précédente (72).