FR 2910131 A1 20080620 FR 0610961 A 20061215 1 LOCALISATION PASSIVE ENSEMBLISTE EN MODES TOA/TDOA DOMAINE DE L'INVENTION La présente invention concerne la résolution des problèmes liés à localisation passive de cibles en modes TOA (Time Of Arrivai) ou TDOA 5 (Time Difference Of Arrivai). Pour le TOA, l'invention permet de localiser une cible quelconque en exploitant les temps d'arrivée, sur un même récepteur, des ondes émises par un ou plusieurs émetteurs et réfléchies par la cible. 10 Pour le TDOA, l'invention permet de localiser une cible émettrice par mesure des différences de temps d'arrivée de l'onde émise sur plusieurs récepteurs synchronisés et dispersés dans l'espace. La présente invention concerne plus particulièrement (mais non 15 exclusivement) le domaine des radars passifs. La présente invention, de par son caractère générique, permet d'adresser également toutes les combinaisons de mode (Multi-TOA, Multi-TDOA, ou encore Mixte TOA/TDOA). CONTEXTE DE L'INVENTION - ART ANTERIEUR Les méthodes de localisation passive, que ce soit en mode TOA ou TDOA, ont pour principe de base de déterminer le positionnement de cibles 25 en utilisant les informations fournies simultanément par différentes sources d'informations. En mode TOA ou mode "Time Of Arrivai" selon la dénomination anglo-saxonne, on s'intéresse aux signaux émis par un ou plusieurs émetteurs et à ces mêmes signaux réfléchis par une cible. La mesure, pour un émetteur 30 donné, du retard entre le trajet direct (émetteur -+ récepteur) et le trajet réfléchi (émetteur -* cible --* récepteur) permet définir une courbe (ou une surface si le problème est traité en 3D) de localisation qui prend la forme d'une ellipse ou d'un ellipsoïde. 20 2910131 2 Ainsi, si l'on exploite plusieurs émetteurs simultanément et pourvu que l'on soit en mesure de recevoir les trajets direct et réfléchi, il est possible de déterminer la position de la cible recherchée en déterminant les points/zones d'intersections des différentes courbes/surface de localisation entre elles. En mode TDOA ou mode "Time Difference Of Arrivai" selon la dénomination anglo-saxonne, on s'intéresse à la localisation de cibles émettrices au moyen de plusieurs récepteurs synchronisés (au moins deux), un récepteur principal et un (ou plusieurs) récepteur(s) secondaire(s). Comme dans le mode TOA la position de la cible émettrice est alors déterminée de manière géométrique à partir des courbes de localisation (hyperboles ou hyperboloïdes) établies à partir des mesures des retards entre les différents signaux provenant de la cible et arrivant à des instants différents sur chacun des récepteurs. IJn des problèmes posés par la localisation passive provient de la prise en compte de l'incertitude de la mesure délivrée par les récepteurs. En effet pour réaliser une localisation précise de la cible il est nécessaire de connaître le plus précisément possible la zone de localisation compatible avec les incertitudes affectant les mesures. Par suite les courbes de localisation sont en réalité des zones de localisation, chaque zone étant située entre deux courbes extrêmes dont l'espacement est fonction de la précision des récepteurs. De la sorte, les points d'intersections des différentes courbes de localisations obtenues idéalement en utilisant plusieurs sources émettrices (mode TOA) ou plusieurs récepteurs (mode TDOA) sont remplacées dans la pratique par des zones d'intersection à l'intérieur desquelles se situent ces points. Dès lors, la recherche des zones de localisation compatibles avec les rnesures passe d'abord par la caractérisation mathématique des dites zones puis par leur recherche dans un espace d'intérêt (i.e. l'espace dans lequel on recherche la présence d'une cible). Les méthodes usuelles utilisées pour réaliser cette recherche sont généralement des méthodes à grilles qui consistent à mailler finement tout l'espace dans lequel on recherche les cibles, c'est à dire un espace 2910131 3 suffisamment vaste pour contenir la zone d'incertitude, et à analyser systématiquement chaque maille pour vérifier si elle appartient à la zone de localisation. Ce sont également des méthodes algébriques du type moindre carré, ou encore des méthodes probabilistes dont la complexité (et donc la 5 complexité de mise en oeuvre) croît avec le nombre de source d'information. Hormis les méthodes à grille, toutes ces méthodes ne permettent pas de restituer finement les zones d'incertitudes liées aux erreurs de mesure (elles ne donnent que des ellipses ou des ellipsoïdes d'erreur). En outre les méthodes à grilles nécessitent une capacité de calcul importante pour traiter 10 chaque maille de manière suffisamment rapide et déterminer suffisamment rapidement la zone de localisation. PRESENTATION DE L'INVENTION 15 Un but de l'invention est de bénéficier des avantages des traitements par grille notamment en terme de résolution, sans pour autant en subir les inconvénients. A cet effet l'invention a pour objet un procédé de localisation passif d'une 20 cible en mode TOA ou TDOA mettant en oeuvre un découpage successif en pavés d'un espace initial (dans lequel on cherche à localiser une cible). L'ensemble des pavés est analysé de manière itérative. A chaque itération chaque pavé d'intérêt est découpé en sous-pavés identiques de taille inférieure. Un pavé d'intérêt est, selon l'invention, un pavé dont au moins un 25 point appartient à la zone de localisation recherchée. Le processus itératif s'arrête lorsque la taille (résolution) des sous-pavés obtenus à l'itération en cours correspond à la résolution souhaitée. L'ensemble des pavés résultant de ce processus 'fournit une approximation de la forme de la zone de localisation recherchée (si cette dernière existe 30 dans l'espace initial). Plus précisément, elle a pour objet un procédé de localisation passif d'une cible en mode TOA ou TDOA mettant en oeuvre un découpage successif en pavés d'un espace initial (dans lequel on cherche à localiser 35 une cible) et une recherche au sein de chaque pavé de la présence de points 2910131 4 appartenant à la zone de localisation recherchée. Cette invention est en outre caractérisée par le fait que le découpage et la recherche sont réalisés sous forme d'étapes itératives sur une sélection de pavés candidats, modifiée itération après itération, de sorte qu'à chaque itération on effectue, 5 parmi les pavés de la sélection obtenue à l'issue de l'itération précédente, la recherche des pavés contenant au moins un point appartenant à la zone de localisation recherchée. Les pavés ne contenant aucun point sont par suite exclus de la sélection, tandis que les pavés contenant au moins un point sont re-découpés en sous-pavés et remplacés dans la sélection par les sous-pavés ainsi formés. La sélection obtenue à chaque itération définit la zone de localisation avec une résolution croissant à chaque itération. Selon l'invention, les itérations sont interrompues lorsque les sous-pavés formant la sélection définissent la zone de localisation avec la 15 résolution souhaitée. Le procédé selon l'invention comporte, en autre, une étape d'initialisation dans laquelle on définit un premier pavé [X ] correspondant à un espace de recherche "a priori", ce pavé constituant la sélection initiale. Selon un mode de mise en oeuvre préféré, le procédé selon l'invention comporte principalement: -une étape d'initialisation dans laquelle on définit un premier pavé 25 IX ] correspondant à un espace de recherche a priori et qui conduit à la constitution d'une liste initiale L formée du pavé [X ], - une étape itérative consistant à: - décomposer chaque pavé [Xn] de la liste courante L '"' en quatre sous-pavés jointifs [Xn1], [X,,2], [X~3] et [Xn4] 30 -rechercher et sélectionner parmi les sous-pavés ainsi constitués ceux contenant au moins un point faisant partie de la zone de localisation, -actualiser la liste courante pour obtenir une liste L ' dans laquelle: 20 2910131 5 -a) les pavés [Xr,] pour lesquels aucun sous-pavé n'a été sélectionné sont supprimés, -b) les autres pavés sont remplacés par les sous-pavés [Xn'] sélectionnés, l'étape itérative comporte en outre une opération de comparaison de la taille ÉL'o(1) des pavés constituant la liste Loi, le procédé étant interrompu lorsque la taille AL'0(1) est supérieure à l'objectif de résolution A bl désiré. 10 Selon ce mode de mise en oeuvre, un pavé [XI est sélectionné si au moins un de ses points vérifie le critère défini comme suit: [x1 [yi]- - 0 E J([XiD avec J([XJ]) _ + en mode "TOA 2D", [a]2 [b]2 - 0 E J([XiD avec J([XiD = [x[yir en mode "TDOA 2D", ù [c]2 [d]2 - 0 E J([XiD avec J([XiD = [xj12. + [y'f + [z1J -1 en mode "TOA [a] [b]2 [b]2 3D", r 12 [yiJ [ziI2 - 0 E J([XiD avec 4)01 4 x _ ù --1 en mode [c]12 [d]2 [d]2 "TDOA 3D" L'invention repose sur l'approche ensembliste utilisée pour localiser la cible. Les mesures issues des différents senseurs sont modélisées par des intervalles (i.e., elles intègrent une erreur bornée). Sachant ces mesures, 25 l'invention permet, par un processus itératif basé sur un critère ensembliste ad hoc (i.e. dépendant du problème, TOA/TDOA mono/multi senseurs) de trouver et d'approcher avec la résolution voulue, toutes les zones de l'espace susceptibles, aux vues des mesures, de contenir une cible. La localisation 5 15 20 2910131 6 est ensembliste dans le sens où l'invention fournit un ensemble solution au problème de la localisation (i.e. un ensemble de position de cible contenant de manière garantie la vraie position de la cible). Contrairement à une méthode à grille, il n'est pas nécessaire de mailler tout 5 l'espace, le processus itératif utilisé dans l'invention permet de se concentrer directement sur les régions d'intérêt. DESCRIPTION DES FIGURES 10 Les caractéristiques et avantages de l'invention seront mieux appréciés grâce à la description qui suit, description qui expose l'invention au travers d'un mode de réalisation particulier pris comme exemple non limitatif et qui s'appuie sur les figures annexées, figures qui représentent: - la figure 1, une configuration géométrique mono-émetteur typique 15 permettant la mise en oeuvre du mode TOA , - la figure 2, une représentation du jeu de courbes de localisation idéalement obtenu en mode TOA avec trois sources émettrices, - la figure 3, une configuration géométrique bi-récepteur typique permettant la mise en oeuvre du mode TDOA, 20 - la figure 4, une représentation du jeu de courbes de localisation idéalement obtenu en mode TDOA avec trois récepteurs secondaires, - la figure 5, la représentation d'une zone de localisation réellement obtenue en mode TOA avec une source émettrice, - la figure 6, la représentation d'une zone de localisation réellement 25 obtenue en mode TDOA avec un récepteur, - la figure 7, un organigramme de principe du procédé selon l'invention, - la figure 8, une illustration de l'application du procédé selon l'invention au mode TOA 2D avec un système mono-émetteur, 30 - la figure 9, une illustration de l'application du procédé selon l'invention au mode TDOA 2D avec un système mono-récepteur secondaire, - la figure 10, une illustration de l'application du procédé selon l'invention au mode TOA 2D avec un système bi-émetteurs, - la figure 11, une illustration de l'application du procédé selon 35 l'invention au mode TDOA 2D avec un système bi-récepteurs secondaires, 2910131 DESCRIPTION DETAILLEE On s'intéresse dans un premier temps à la figure 1 qui présente de 5 manière schématique le principe de fonctionnement idéal du mode de localisation TOA. Pour des raisons de clarté de l'exposé le modèle illustré ici est un modèle à deux dimensions correspondant à l'exploitation des signaux issus d'un seul émetteur (modèle mono-émetteur 2D). Comme l'illustre la figure, ce mode de fonctionnement "mono- émetteur" met en jeu une source émettrice 12 (émetteur), un récepteur 11, et une cible 13 que l'on cherche à localiser. Pour localiser la cible, le récepteur 11 mesure le retard temporel existant entre l'onde reçue en provenance directe de l'émetteur 12 (onde directe), et l'onde reçue provenant de la réflexion de l'onde émise par la source 12, sur la cible 13 que l'on cherche à 15 localiser (onde réfléchie). Dans une telle configuration, si on appelle L la distance entre la source émettrice 12 et le récepteur 11, RT la distance entre la source 12 et la cible 13 et RR la distance entre la cible et le récepteur, la localisation de la cible 13 20 consiste à effectuer la mesure du retard entre les ondes directe et réfléchie et à déterminer à partir de ce retard, par toute méthode appropriée connue, la distance Rb parcourue par l'onde réfléchie, définie par: Rb = RT + RR La cible 13 est alors localisée comme l'illustre la courbe 21 de la figure 2 sur une ellipse d'ambiguïté ayant pour foyers 22 et 23 la position de la source émettrice 12 et celle du récepteur 11 et dont le grand axe (axe semi-majeur) a pour longueur Rb/2. 30 Par suite, si on définit un repère cartésien xOy centré sur le milieu du segment de droite [Rx, Tx] 14 reliant la sourcel2 et le récepteur 11 et dont le vecteur Ox est colinéaire au vecteur RXTx , l'équation de cette ellipse de localisation 21 dans le repère xOy a pour expression la relation suivante: 7 25 35 2910131 x2 y2 a2 +b2 avec: 5 a = Rb/2 et b = /a2 ûL2/4 Pour affiner la localisation de la cible 13, il est bien entendu nécessaire de 10 disposer de plusieurs source d'émission. De la sorte pour une même cible 13 on peut associer à chaque source une ellipse de localisation, les intersections de ces ellipses définissant les positions possibles de la source à localiser. Les courbes 21 et 24 de la figure 2 illustrent les résultats obtenus avec une configuration "bi-émetteurs 2D" à partir de laquelle on obtient deux 15 ellipses de localisation dont les foyers 22, 23 et 25 sont respectivement le récepteur et la première source pour l'ellipse 21, et le récepteur et la deuxième source pour l'ellipse 24. L'ensemble des lieux possibles de localisation de la cible 13 est alors constitué des deux points d'intersection 26 et 27 des deux courbes. 20 Par suite, pour déterminer la position de la cible 13 sans ambiguïté il est nécessaire de disposer d'au moins une source d'émission supplémentaire (configuration "multi-émetteurs 2D") l'intersection des trois ellipses 21, 24, et 28 de localisation définissant un point commun unique, le point 26 par exemple, sur lequel est située la cible 13. 25 On s'intéresse ensuite à la figure 3 qui présente de manière schématique le principe de fonctionnement idéal du mode de localisation TDOA. Pour des raisons de clarté de l'exposé le modèle illustré ici est, comme pour la présentation du mode TOA, un modèle à deux dimensions correspondant à 30 l'exploitation des signaux issus d'une cible émettrice 31 et reçus par deux récepteurs 32 et 33 distants l'un de l'autre, un récepteur de référence dit récepteur principal et un récepteur, dit récepteur secondaire, synchronisé sur le récepteur de référence (modèle mono-récepteur secondaire 2D). 8 [1] 2910131 9 Comme l'illustre la figure 3, ce mode de fonctionnement "mono-récepteur secondaire" met en jeu une cible émettrice 31 (émetteur) que l'on cherche à localiser, un récepteur 32 de référence, et un récepteur secondaire 33. Pour localiser la cible, on exploite alors le retard temporel existant entre 5 l'instant de réception par le récepteur de référence de l'onde émise par la cible 31 et l'instant de réception de cette même onde par le récepteur secondaire 33. Dans une telle configuration, si on appelle L la distance entre les deux 10 récepteurs 32 et 33, RR la distance entre la cible émettrice 31 et le récepteur de référence 32 et RR1 la distance entre la cible émettrice 31 et le récepteur secondaire 33, la localisation de la cible 31 consiste à effectuer la mesure du retard entre les ondes reçues par les deux récepteurs et à déterminer à partir de ce retard, par toute méthode appropriée connue, l'écart de distance 15 parcourue Re définie par: Rd1 = RR1 û RR La cible 31 est alors localisée comme l'illustre la courbe 41 de la figure 4 sur 20 une hyperbole ayant pour foyers 42 et 43 la position du récepteur de référence et celle du récepteur secondaire et dont la distance inter-sommets a pour valeur Rd1/2. Par suite, si on définit un repère cartésien xOy centré sur le milieu du segment de droite [RX, Rx1] 34 reliant le récepteur de référence 32 et le 25 récepteur secondaire 33 et dont le vecteur Ox est colinéaire au vecteur RXRx1 , l'équation de cette hyperbole de localisation 41 dans le repère xOy a pour expression: x2 y2 c2 avec: c = Rd-1/2 [2] 30 2910131 et d = 1jL2/4 û c2 Pour affiner la localisation de la cible 31, il est bien entendu nécessaire de disposer de plusieurs récepteurs secondaires. De la sorte pour une même 5 cible 31 on peut associer à chaque couple récepteur de référence/récepteur secondaire, une hyperbole de localisation. Les intersections de ces hyperboles définissent alors les positions possibles de la cible. Les courbes 41 et 44 de la figure 4 illustrent les résultats obtenus avec une configuration "bi-récepteurs secondaires 2D" à partir de laquelle on obtient deux 10 hyperboles de localisation dont les foyers 42, 43 et 45 sont respectivement le récepteur Rx de référence et le récepteur secondaire Rd pour l'hyperbole 41 et le récepteur de référence Rx et un second récepteur secondaire Rx2 pour l'hyperbole 44. L'ensemble des lieux possibles de localisation de la cible 31 est alors constitué des points d'intersection 46 et 47 des deux courbes. 15 Par suite, pour déterminer la position de la cible 31 sans ambiguïté il est nécessaire de disposer d'au moins un récepteur secondaire supplémentaire (configuration "multi-récepteurs secondaires 2D") l'intersection des trois hyperboles 41, 44, et 48 de localisation définissant un point commun unique, le point 46 dans l'exemple, sur lequel est située la cible 31. 20 Le principe théorique de détermination exposé dans les paragraphes précédents au travers d'une localisation en deux dimensions (localisation "2D") peut naturellement être généralisé à une localisation dans l'espace (ie en localisation 3D). 25 I En mode TOA l'ellipse de localisation en mode "mono-émetteur" est remplacée par un ellipsoïde dans l'espace, ellipsoïde qui peut être représenté par l'équation suivante: x2 y2 z2 a2 +b2 +b2 =1 avec: a=Rb/2 10 [3] 30 2910131 et b =.Ja2 ùL2/4 De même, en mode TDOA, l'hyperbole de localisation en mode "mono-récepteur secondaire" est remplacée par un hyperboloïde dans 5 l'espace, hyperboloïde qui peut être représenté par: x2 y2 z2 =1 c2 d2 d2 [4] avec: c = Rdl/2 et d=,jL2/4ùc2. 15 On s'intéresse à présent aux figures 5 et 6 qui illustrent au travers des cas simples de localisation en mode TOA 2D (figure 5) et en mode TDOA 2D (figure 6) le problème posé par la précision des mesures réelles obtenues avec les récepteurs. 20 Comme tout dispositif de mesure, les récepteurs utilisés par les systèmes de localisation passive fournissent des mesures entachées d'une certaine imprécision que l'on peut supposer bornée. Cette imprécision conduit à représenter les mesures de distance effectuée en mode TOA ou TDOA, non plus par des valeurs exactes, mais par des intervalles dont la 25 taille correspond à l'erreur maximale de mesure. Ces intervalles [x] sont définis par la relation suivante: [x] = [x-, x+ ] = {x e Î/x- x x+ } [5] 30 où f représente l'ensemble des nombres réels. Si Rb représente l'erreur maximale sur la mesure de Rb en mode TOA, une mesure de Rb à erreur bornée se défini par l'intervalle: 11 10 2910131 12 [Rb-, Rb ] := [ Rb -QRb , Rb +GRb ] [6] Par suite, la prise en compte d'une erreur bornée sur la mesure de Rb en mode "TOA 2D" conduit à la localisation de la cible considérée non pas 5 sur une courbe mais sur une surface de localisation délimitée par deux ellipses confocales 51 et 52 de paramètres respectifs Rb- et Rb+, telle que celle illustrée sur la figure 5. Cette surface 53 peut être décrite par la forme paramétrique suivante: {(x y)e92 x2 y2 1E [a]2 +[b]2 [7] 10 15 où: [a] _ [Rb ]/2 et [b] _ \[a]2 - L2/4 sont évalués en appliquant les règles du calcul par intervalle. 20 De manière analogue, si l'on effectue une localisation dans l'espace ("TOA 3D"), la prise en compte des erreurs bornées sur la mesure de Rb conduit à la localisation de la cible considérée dans un volume défini par deux ellipsoïdes confocaux, volume qui peut être décrit par la forme paramétrique suivante: 25 (x, y, z) e 9i 3 x2 y2 z2 1E [a]2 + [b]2 + [b]2 [8] Comme pour le mode TOA, la prise en compte d'une erreur bornée sur la mesure de Rdi en mode "TDOA 2D" conduit à la localisation de la cible 30 considérée non pas sur une courbe mais sur une surface de localisation délimitée par deux hyperboles confocales 61 et 62 de paramètres respectifs Rd1- et Rd1+ définis par l'intervalle [RC1 , Rd1+] := [ Rdl -QRdi , Rdl +QRd, ], telle 2910131 13 que celle illustrée à la figure 6. Cette surface 63 peut être décrite par la forme paramétrique suivante: 1E 5 avec: (x,y,z)E912 x2 ù y2 L ~c]2 [d]2 J [9] [c] _ [R dl 1/2 10 et [d] _./L2/4 ù[c]2 Dans le cas d'une localisation dans l'espace ("TDOA 3D"), la prise en compte des erreurs bornées sur la mesure de Rdl conduit à la localisation de la cible considérée dans un volume défini par deux hyperboloïdes confocaux, volume 15 qui peut être décrit, de manière analogue, par la forme paramétrique suivante: 1E x2 y2 z2 [c]2 [d]2 [d]2 (x, y, z) E i 3 [10] 20 On s'intéresse à présent à la figure 7 qui présente sous forme schématique le principe du procédé selon l'invention. Afin de rendre l'exposé du principe de fonctionnement du procédé selon l'invention plus clair, ce principe est décrit ici de manière détaillée pour le cas particulier de la recherche d'une zone d'une localisation au moyen de systèmes de 25 localisation de type "mono-émetteur" (localisation en mode TOA) ou "mono-récepteur secondaire" (localisation en mode TDOA). Le principe de base du procédé consiste en un affinage progressif de la zone de localisation. Il consiste à découper l'espace que l'on désire 30 analyser (i.e le pavé [xo]x[yo] initial de recherche) en sous pavés jointifs dans lesquels la présence d'une cible est évaluée au moyen d'un critère ad-hoc. Les sous pavés dans lesquels la présence d'une cible est avérée sont à leur 2910131 14 tour découpés, les autres sont rejetés. Le processus itératif ainsi créé se répète tant que la présence d'une cible est avérée dans les pavés en cours d'analyse et qu'un critère d'arrêt (correspondant à un objectif de largeur des pavés) n'est pas atteint. 5 Le procédé selon l'invention comporte à cet effet plusieurs étapes: - une étape 71 d'initialisation - une étape 72 de calcul itératif 10 L'étape 71 d'initialisation consiste à définir un pavé initial, [X0 ], correspondant à un espace de recherche a priori. En mode 2D on a ainsi: [Xo]=[xo]x{yo], 15 et en mode 3D : [Xo]_ [xo]x [y ]x [zo]. 20 L'étape 71 consiste également à initialiser une liste Lo, de taille variable au cours de la mise en oeuvre du procédé, comportant la liste des pavés à étudier. Le contenu de Lo est initialisé avec [Xo ] L'étape 71 consiste encore à définir un critère d'arrêt du procédé. Ce critère 25 d'arrêt est donné ici par la résolution "objectif'Aobi avec laquelle on souhaite finalement caractériser la zone de localisation. Cette résolution est naturellement limitée par la précision des mesures fournies par les récepteurs, mais peut être arbitrairement fixée dans cette limite. En mode de localisation 2D elle est peut être définie dans un repère xOy par 30 4 obj ( obJ LobJ ). x, y L'étape 71 est suivie d'une étape 72 réalisant un traitement itératif qui comporte deux boucles de traitement imbriquées, une boucle principale 73 et une boucle secondaire 74. 35 2910131 15 La boucle principale 73 consiste à actualiser la liste Lo établie lors de l'initialisation avec les résultats du traitement réalisé par la deuxième boucle. Ainsi, à chaque itération i de la boucle principale on dispose d'une liste Lo réactualisée, notée L'o. 5 A chaque itération, chaque pavé [Xôk]= L'o(k) constituant la liste L'o est découpé en N sous-pavés jointifs (N=4 en mode 2D, N=8 en mode 3D) et regroupés en une listes Li. La boucle secondaire 74 consiste quant à elle, à supprimer de LI , les pavés 10 incompatibles avec la zone de localisation recherchée (c'est à dire les pavés ne validant pas le critère). Par suite, l'élément [Xô k j de Lia est remplacé, dans L'o , par la liste Li. 15 Selon l'invention, la méthode choisie pour déterminer si un sous-pavé [Xi] inclut une partie de la zone de localisation consiste à déterminer si un ou plusieurs points du sous-pavé appartiennent à cette zone. Pour ce faire dans les systèmes "mono-émetteur" ou "mono-récepteur secondaire" on considère simplement la grandeur J([Xi]), définie selon le 20 mode de localisation considéré par: - en mode "TOA 2D": JQXiD = [xi +Lys ù [a] [b]2 - en mode "TOA 3D": JQXJD = [xir +[y'ï [zi2 -1 [a] [b]2 [b]2 [12] - en mode "TDOA 2D": JQXiD= [13] [y' -1 [c]2 û [d]2 25 - en mode "TDOA 3D": JQXiD = [x'1 _ [y'1 ù [zu]2 -1 [14] [c] [d] [d] 2910131 16 avec [Xi]= [x1]x [y-] pour les modes "2D", et [Xi]= [xi]x [yi]x [zi] pour les modes "3D" 5 Par suite pour vérifier si au moins un point d'un pavé [Xi] appartient à la zone de localisation, il suffit, selon l'invention, de vérifier si: OEJ([XiD 10 La boucle principale est par nature une boucle sans fin pour laquelle il est nécessaire de déterminer une condition d'arrêt et une opération de test de l'apparition de cette condition d'arrêt. Selon l'invention cette condition d'arrêt est initialisée lors de l'étape 71 et porte sur la résolution A bJ avec 15 (laquelle on souhaite définir la zone de localisation. Aussi la boucle principale '73 comporte une opération de test 75 exécutée en fin de traitement à chaque itération i. Cette condition est définie par la relation suivante: 20 où A1-o représente la résolution du premier pavé constituant la liste L'o , chaque pavé de la liste ayant une résolution identique à ce stade. Ainsi, tant que la condition de la relation [11] se trouve vérifiée, l'étape 72 de calcul itératif se répète. L'itération ne cesse que lorsque la résolution 25 souhaitée est atteinte. Ainsi, appliqué à un traitement de localisation TOA ou TDOA à deux dimensions et pour un découpage itératif de chaque pavé en P = 4 sous-pavés le procédé selon l'invention peut être décrit par l'enchaînement 30 d'actions suivant: 1. Définition de la taille initiale d'une cellule d'analyse (pavé): AL' o( 1) REVENDICATIONS 1. Procédé de localisation passif d'une cible en mode TOA ou TDOA mettant en oeuvre un maillage (découpage) en pavés de l'espace dans lequel se situe la zone de localisation et une recherche au sein de chaque pavé de la présence de points appartenant à cette zone de localisation considérée, caractérisé en ce que le découpage et la recherche sont réalisés sous forme d'étapes itératives sur une sélection de pavés candidats, modifiée itération après itération, de sorte qu'à chaque itération on effectue, parmi les pavés de la sélection obtenue à l'issue de l'itération précédente, la recherche des pavés contenant au moins un point appartenant à la zone de localisation, les pavés ne contenant aucun point étant par suite exclus de la sélection, tandis que les pavés contenant au moins un point sont redécoupés en sous-pavés et remplacés dans la sélection par les sous-pavés ainsi formés, la sélection obtenue à chaque itération définissant la zone de localisation avec une précision croissant à chaque itération. 2. Procédé selon la revendication 1, comportant en outre une étape d'initialisation dans laquelle on définit un premier pavé [X ] correspondant à un espace de recherche a priori, ce pavé constituant la sélection initiale. 3. Procédé selon l'une des revendications 1 ou 2, dans lequel les itérations sont interrompues lorsque les sous-pavés formant la sélection définissent la zone de localisation avec la résolution souhaitée. 4. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'il comporte: - une étape d'initialisation (71) dans laquelle on définit un premier pavé [X ] correspondant à un espace de recherche a priori et qui conduit à la constitution d'une liste initiale L0 formée du pavé [X ], - une étape itérative (72) consistant à: 2910131 23 - décomposer (73) chaque pavé [Xn] de la liste courante Loi-1 en quatre sous-pavés jointifs [Xn1], [Xn2], [xn'i et [x,4], - rechercher et sélectionner (73) parmi les sous-pavés ainsi constitués ceux contenant au moins un point faisant partie de la zone 5 de localisation, - actualiser (73) la liste courante pour obtenir une liste Loi dans laquelle: - les pavés [Xn] pour lesquels aucun sous-pavé n'a été sélectionné sont supprimés, 10 - les autres pavés sont remplacés par les sous-pavés [Xn'] sélectionnés, l'étape itérative comportant un test (75) de comparaison de la résolution do(l) des pavés constituant la liste Loi, le procédé étant 15 interrompu lorsque la résolution do(l) est supérieure à la résolution Dobj désirée. 5. Procédé selon la revendication 4, pour lequel, en mode de localisation "TOA 2D mono émetteur", un sous pavé [Xn'] est 20 sélectionné s'il vérifie le critère défini par: - 0 E J([XJDavec J([XjD= Lx1J+ Ly, -1 [a]2 [b]2 6. Procédé selon la revendication 4, pour lequel, en mode de 25 localisation "TDOA 2D mono-récepteur secondaire", un sous pavé [Xn'] est sélectionné s'il vérifie le critère défini comme suit: - 0 E J([X'Iavec JQXj [Xi 2 ù [yj]2 -1 [c]2 [d]2 30 7. Procédé selon la revendication 4, pour lequel, en mode de localisation "TOA 3D mono émetteur", un sous pavé [Xn'] est sélectionné s'il vérifie le critère défini comme suit: 2910131 24 - 0 E J([Xi1) avec JExO = LX1f2 + [yi 2 + ù [a] [b]2 [b]2 8. Procédé selon la revendication 4, pour lequel, en mode de 5 localisation "TDOA 3D mono-récepteur secondaire", un sous pavé [Xn'] est sélectionné si au moins un de ses points vérifie le critère défini comme suit: [Xj~ û [y;J [Z;]2 - 0 E J([Xi-) avec JQXJIl) _ -1 [c]2 [d]2 [d]2 10 9. Procédé selon la revendication 4, pour lequel, en mode de localisation "TOA multi-émetteurs" ou en mode de localisation "TDOA multi-récepteurs secondaire", un sous pavé [Xn'] est sélectionné s'il vérifie le critère A défini par la relation : 15 -A=ra. =N dans laquelle N représente le nombre d'émetteurs ou de récepteurs secondaires constituant le système de localisation TOA ou TDOA 20 considéré et dans laquelle les termes ai sont définis comme suit: 1si0EJ1([X]) a' Il 0 sinon f 1 si 0 J2 ([X]) a2 l0 sinon {1 si 0 e JN([X]) aN= 0 sinon 25 2910131 25 Où J;([X]) représente pour chaque système unitaire "mono-émetteur" (Rx, Txi) ou "mono-récepteur secondaire" (Rx, Rxi), la grandeur définie par: _ XJ 2 f y1]~ + L 2 -1 en mode "TOA 2D", [a] [b]2 [ 2 ' 1 2 + y 2 + [Z ]2 -1 en mode "TOA 3D", [a] [b] [b] [xi ]2 yi ù -1 en mode "TDOA 2D", [c] [d]2 1- 2 i2 [Zf 12 [X~ `y-ù -1 en mode "TDOA 3D". [c] [d]2 [d] 10 10. Procédé selon la revendication 4, pour lequel, en mode de localisation "TOA multi-émetteurs" ou en mode de localisation "TDOA multi-récepteurs secondaire", pour un nombre N d'émetteurs ou de récepteurs secondaires constituant le système de localisation TOA ou 15 TDOA considéré, un sous pavé [Xni] est sélectionné s'il vérifie le critère A défini par la relation : N - A=. >M J=1 20 dans laquelle M représente un nombre déterminé, inférieur à N, d'émetteurs ou de récepteurs secondaires constituant le système de localisation TOA ou TDOA considéré et dans laquelle les termes aj sont définis comme suit: 1 si 0 E J, ([X]) a' 0 sinon 5 - J([XiD = - J([XiD = 25 2910131 26 5 1 si 0 E J2([X]) a2 0 sinon {1 si 0 E JN([X]) aN 0 sinon Où Ji([X]) représente pour chaque système unitaire "mono-émetteur" (RX, TXi) ou "mono-récepteur secondaire" (RX, RXi), la grandeur définie par: [x112+[y,ù1 [a] [b]2 en mode "TOA 2D", -J([Xii)= 10 [x'- ]2 +[yi2 [z'12 ]2 -1 [a]2 [b] [b]en mode "TOA 3D", -J([XiJ) [xj- ]2 [y'J [c]2 [d]2 en mode "TDOA 2D", 1 - J([Xi]) _ [x'- ]2 [yif-[Z']2-1 [c]2 [d]2 [d]2en mode "TDOA 3D". 15 11. Application du procédé selon la revendication 9 ou 10 à la localisation passive d'une cible par un système global constitué de systèmes unitaires "mono-émetteur" et "mono-récepteur secondaire" et fonctionnant dans un mode mixte TOAITDOA. 20