FR 3029300 A1 20160603 FR 1402710 A 20141128 PROCEDE DE LOCALISATION PASSIVE D'UN EMETTEUR NON MOBILE La présente invention concerne le domaine de l'écoute électromagnétique. La présente invention concerne plus particulièrement un 5 procédé de localisation passive d'un émetteur non mobile. Dans un contexte aéroporté ou spatial, la localisation passive d'émetteur peut être effectuée par différentes méthodes comme par exemple, par des mesures d'angle d'arrivée (ou AOA pour "Angle Of 10 Arrivai"), d'évolution de phase sur une grande base (ou LBI pour "Long Base Interferometer" ou LBPDE pour "Long Base Phase Difference Evolution"), de différences de temps d'arrivée (ou TDOA pour " Time Difference Of Arrivai") entre au moins deux porteurs, ou encore par Différence de Temps de Passage de Lobe (ou DTPL). 15 Les mesures d'angle d'arrivée sont couramment utilisées pour la localisation passive d'émetteurs. Elles nécessitent que les récepteurs soient équipés soit d'un goniomètre d'amplitude, soit d'un interféromètre. Le défaut principal de la localisation par mesure d'angle d'arrivée est son manque de 20 précision quand le temps alloué à la localisation est court. Si l'on dispose d'un interféromètre sur grande base, on peut utiliser des mesures de LBPDE cependant, ces mesures sont sujettes à des ambigüités. De plus, les techniques actuelles de localisation par évolution de 25 phase sur une grande base ou par différence de temps d'arrivée, ne prennent pas en compte le caractère périodique des formes d'onde radar les plus fréquemment rencontrées. Il en résulte, pour ces formes d'ondes, des performances sous-optimales aussi bien en termes de précision de localisation qu'en termes de charge de calcul et d'utilisation de la bande 30 passante entre les porteurs. Il est connu de l'art antérieur, notament par "An Algebraic Solution of the GPS Equations" de Stephen Bancroft (IEEE Transactions AES VOL. AES-21, NO. 7 JANUARY 1985), une méthode de localisation dite de 35 Bancroft. Dans cette publication, une estimation explicite (non itérative) en moindres carrés de la position est exposée dans un contexte similaire du GPS. La méthode ne fait pas d'hypothèse particulière sur la nature du signal. Cependant, si on tente d'étendre directement la méthode présentée au cas d'un train d'impulsion périodique, on trouve un système de trois équations du second degré qui n'est pas résoluble explicitement en général. De plus, de nombreuses solutions non valables doivent être éliminées si le système est résolu numériquement. Un problème avec la localisation passive d'émetteurs à l'aide des temps d'arrivée des impulsions (ou TOA pour " Time Of Arrivai") vient du nombre de mesures à traiter. Lors de l'interception d'un signal, une mesure de TOA est associée à chaque impulsion reçue. Le nombre d'impulsions interceptées, donc de mesures de temps d'arrivée d'impulsions, devient rapidement très grand car ces impulsions sont émises à une cadence très élevée. De plus, avec les solutions existantes, calculer une localisation à partir des TOA implique d'estimer le temps d'arrivée de chaque impulsion. La complexité du problème d'estimation va donc croitre rapidement avec le nombre d'impulsions. Par exemple, pour un algorithme explicite basé sur la méthode des moindres carrés linéaires, la complexité asymptotique est 20 donnée par 0(12M) avec I le nombre de paramètres inconnus à estimer et M le nombre de mesures. Dans le cas d'estimation avec les temps d'arrivée on a M=P (nombre de récepteurs) x N (nombre d'impulsions) et 1=3 (dimensions spatiales) +N (nombre de dates d'émission inconnues, soit nombre d'impulsions) soit une complexité asymptotique en 0(N3P) qui est prohibitive 25 au vu du nombre d'impulsions susceptibles d'être reçues. Concernant la localisation basée sur les différences de temps d'arrivée des impulsions, comme on calcule une différence entre les mesures, on n'a plus à estimer les dates d'émission de chaque impulsion. 30 Cependant, le calcul de cette différence augmente la variance des mesures de TDOA par rapport à celles de TOA et fait apparaitre des corrélations entre les mesures. La matrice de covariance des mesures de différences de temps d'arrivée des impulsions est : _ ,2 2 1 1 ETDOA 'TOA - 1 2 11 1 1 2] 1 Au lieu de : TOA = tr7'0A [1 ... 0 1 0 0 0 -.. 0 0 ... L'augmentation des termes diagonaux et l'introduction de termes hors-diagonale dans la matrice de covariance des mesures conduit à une légère 5 dégradation des performances de localisation. De plus, l'élimination du paramètre "date d'émission " correspondant à la date de début d'émission des impulsions, enlève la possibilité d'introduire une connaissance a priori sur ce paramètre. Or un a priori permet d'introduire de l'information supplémentaire dans l'estimateur ce qui permet d'améliorer 10 les performances de la localisation finale. Un autre problème de la localisation passive par TOA ou TDOA est lié aux communications entre les plateformes. L'une ou l'autre des deux méthodes précédentes implique une communication entre les plateformes 15 pour calculer les résultats de la localisation. Les informations transmises d'une plateforme vers l'autre sont les temps d'arrivée des impulsions reçues. Comme le nombre d'impulsions reçues est important et que les temps d'arrivée doivent être transmis précisément (codés sur un nombre de bits importants), on a donc une forte contrainte sur le débit de la liaison de 20 données. De plus ces transmissions pénalisent la discrétion des plateformes. Un but de l'invention est notamment de corriger tout ou partie des inconvénients de l'art antérieur en proposant une solution permettant d'améliorer la précision de localisation d'un émetteur non mobile sans 25 pénaliser la discrétion des détecteurs. A cet effet, l'invention a pour objet un procédé de localisation passive d'un émetteur non mobile au sol mis en oeuvre par un groupe d'au moins une station réceptrice, chacune des stations réceptrices comprenant un radar passif et une référence de temps, l'ensemble des références de temps étant synchronisées entre elles et ledit émetteur émettant un ensemble d'impulsions périodiques, ledit procédé étant caractérisé en ce qu'une première estimation de la position dudit émetteur est réalisée par la méthode de Bancroft à partir des temps moyens d'arrivée des impulsions émise par l'émetteur au niveau de chaque station du groupe d'au moins une station réceptrice, le résultat obtenu étant ensuite utilisé comme point d'initialisation d'une méthode de maximum de vraisemblance afin de converger vers la position dudit émetteur. Selon une variante de mise en oeuvre, le procédé est mis en oeuvre par au moins deux stations réceptrices et ledit procédé comprend : Une étape Etpl de mesure des temps d'arrivés des impulsions de l'émetteur par chaque station réceptrice, Une étape Etp2 de calcul du temps moyen d'arrivée des impulsions par chaque station réceptrice, Une étape Etp3 de transmission des temps moyens d'arrivée des impulsions calculés à une station réceptrice prédéterminée, Une étape Etp4 de calcul, par ladite station réceptrice prédéterminée, de la localisation de l'émetteur à partir de tous les temps moyens d'arrivée des impulsions à l'aide de la méthode de Bancroft puis de la méthode du maximum de vraisemblance. Selon une variante de mise en oeuvre, la ou les stations réceptrices sont sont des stations réceptrices aéroportées. Selon une variante de mise en oeuvre, la ou les stations réceptrices sont sont des stations réceptrices spatiales. Les avantages fournis, outre l'optimalité de l'estimation dans les cas des émetteurs périodiques, sont relatifs à la quantité restreinte de données échangées entre les stations réceptrices ainsi qu'à la précision de localisation atteignable en un temps aussi court qu'un train d'impulsions radar. Par ailleurs un avantage important apporté par l'invention est une initialisation pertinente des algorithmes itératifs d'estimation de la position à partir des mesures de temps d'arrivée. D'autres particularités et avantages de la présente invention 5 apparaîtront plus clairement à la lecture de la description ci-après, donnée à titre illustratif et non limitatif, et faite en référence aux dessins annexés, dans lesquels : La figure 1 illustre un exemple de mise en oeuvre du procédé de 10 localisation selon l'invention ; La figure 2 représente les temps d'arrivées d'impulsions au niveau d'une station réceptrice à différents instants ; La figure 3 illustre des étapes possibles du procédé de localisation selon l'invention ; 15 La figure 4 représente la configuration utilisée pour les résultats présentés ; La figure 5 représente des exemples d'ellipses à 90% de confiance obtenues à l'aide du procédé selon l'invention. 20 Le principe de la localisation passive selon l'invention repose sur une estimation de la position d'un émetteur non mobile au sol en deux temps à l'aide de deux méthodes de calcul différentes. Classiquement, le calcul de la position d'un émetteur est effectué à 25 l'aide d'une estimation du maximum de vraisemblance. Le problème est que cette méthode est difficile à appliquer directement car elle implique des équations non linéaires compliquées à résoudre. L'estimateur maximise une fonction qui comporte plusieurs maxima locaux. De plus, cet algorithme est très sensible aux conditions initiales et a besoin d'avoir une bonne 30 initialisationpour converger vers le maximum global de vraisemblance. L'invention propose d'effectuer en deux temps l'optimisation du critère de maximum de vraisemblance. Une première estimation sommaire est faite avec la méthode de Bancroft connue de l'homme du métier notamment par "An Algebraic Solution of the GPS Equations" de Stephen Bancroft (IEEE Transactions AES VOL. AES-21, NO. 7 JANUARY 1985). Cette première estimation est effectuée en supposant que les mesures ne sont pas bruitées. La méthode de Bancroft a pour avantage d'être explicite et donc non itérative. Elle permet de réduire les équations à un problème des moindres carrés. L'estimation obtenue est ensuite utilisée comme point d'initialisation pour l'algorithme de maximum de vraisemblance qui lui est itératif. L'utilisation d'une première estimation sommaire a pour avantage de fournir un point d'initialisation, pour l'algorithme de maximum de vraisemblance, qui est suffisamment proche de la solution pour converger vers le maximum global et non vers un maximum local. De plus, cela permet de réduire le nombre d'itérations de l'algorithme de maximum de vraisemblance. En référence aux figures 1 et 2, on considère P stations réceptrices 10, comme par exemple des stations réceptrices aéroportées ou spatiales, avec P un entier strictement positif quelconque. Chaque station réceptrice 10 comprend un radar passif et une référence de temps configurée de sorte que l'ensemble des références de temps sont synchronisées entre elles. On suppose que l'émetteur 15 émet un train de N impulsions périodiques (N représentant un entier strictement positif) avec une période de récurrence notée T. Dans ces conditions, le temps d'émission de la nième impulsion peut s'écrire : to + nT où to représente le temps du milieu de la rafale d'impulsions et où 1-N N-1 n = - ...1-2. Afin de simplifier les expressions, on considèrera le 2 nombre d'impulsions N impair. Bien entendu, le raisonnement peut être généralisé à un nombre entier quelconque N strictement positif. Le temps d'arrivée tr,,p de l'impulsion d'indice n sur la pième station réceptrice (avec p=1,...,P) est le temps d'émission de l'impulsion ajouté au 30 temps de propagation vers la station, avec un bruit de mesure. Ce temps d'arrivée peut s'écrire : tn,p = rp t\n,p tO nT + wn,p (1) Où : c représente la vitesse de propagation, rp(t,p) représente la distance entre l'émetteur et la pième station réceptrice au temps d'arrivée tn,p to représente le temps d'émission de l'impulsion, T représente la période de répétition des impulsions wn,p représente le bruit de mesure supposé gaussien, centré et indépendant en n (d'une mesure à l'autre) et en p (d'un récepteur à l'autre) de variance notée 010A - On suppose que la durée d'un train d'impulsions est suffisamment courte pour que chacune des P stations réceptrices ait toujours une trajectoire 12 assimilable à un mouvement rectiligne uniforme. Sous ces hypothèses on peut vérifier que : tn,p ap nT + wn3, , (2) 0 tù : a - + 0 , P c 1-0,p représente la distance entre l'émetteur et la pième station réceptrice au milieu de sa trajectoire, T représente la période de récurrence ou période de répétition des impulsions ; wn,p représente le bruit lié à la nième impulsion sur la pième station réceptrice. Le modèle (2) va permettre d'estimer la période de récurrence des impulsions T et le temps moyen d'arrivée des impulsions sur la station 25 réceptrice d'indice p, p variant de 1 à P. La relation non linéaire entre (ap)p.i...p et les coordonnées de l'émetteur, implicites dans ro,p, pourra ensuite être inversée par un algorithme explicite ce qui fournira une initialisation pertinente des algorithmes itératifs 30 d'estimation fine de la position de la cible. On considère le cas particulier de la localisation passive d'un émetteur 15 par une seule station réceptrice 10. Si P=1 on peut mettre les équations (2) sous la forme vectorielle suivante : Avec : t_ 1 N-] [v_ N-1 il T1 = 111 = 2 01 = [T] et = 2 W N-1 t N-1 1 1 2, 2 Les vecteurs ri, 0 et w1 sont de taille N et la matrice H1 est de taille Nx2. Les composantes du bruit w1 étant Gaussiennes, centrées et indépendantes, l'estimation du maximum de vraisemblance pour les paramètres 01 est donnée par la résolution du problème des moindres carrés dont la solution est classiquement donnée par : = (1/11/0-11/Ir1. La matrice Hi*Hi est diagonale : 0 1-111-11 = [0 (N-1)N(N+1) 12 D'où les estimations : et N-1 2 al = N-1 n=- 2 N-1 T= 2 12 - l)N(N + 1) ntn'i -1 n=- Où â1 représente l'estimation du temps moyen d'arrivée des impulsions sur la station réceptrice et l'émetteur et D représente l'estimation 20 de la période de récurence. On considère maintenant le cas d'un nombre de stations réceptrices 10 au moins égal à 2. Si P>1, on peut mettre les équations (2) sous une forme vectorielle similaire à l'équation (3) précédente de la manière suivante : Ti) = HpOp Wp . Où, par blocs de taille N associés à chacune des P stations réceptrices, on a: Les vecteurs Tp, Op et wp sont de taille NP et la matrice Hp est de taille 10 NPx(P+1). De même que dans le cas précédent où P=1, les composantes du bruit wp sont Gaussiennes, centrées et indépendantes. Si leur variance est la même pour toutes les stations, l'estimation du maximum de vraisemblance pour les paramètres Op est encore donnée par la résolution du problème des moindres carrés dont la solution est : (gp = (J;Hp)-1H;Tp 15 La matrice Hp*Hp est diagonale : Tp = , Hp = Ft_ N1 -1 2-, t 1 2 Ft_13 N;1, t 2 Op = [a'pl et wp = Hplp = 0 NO" --- 0 1'.I 0 0 P(N-1).N(N+1) 12 D'où les estimations : 12 P 1 N-1 p=1 Où etp représente l'estimation du temps moyen d'arrivée des impulsions sur la station réceptrice d'indice p et l'émetteur et D représente 5 l'estimation de la période de récurrence des impulsions. Les expressions obtenues montrent que le procédé d'estimation peut être mis en oeuvre en calculant le temps moyen d'arrivée sur chacune des stations puis en mettant en commun les résultats obtenus pour en déduire 10 une estimation de la position de la cible. L'estimation de la période de récurrence est donnée quant à elle par la moyenne des P estimations faites individuellement par chaque station. Si la variance du bruit n'est pas la même d'une station à l'autre, l'estimation de ep est modifiée en la résolution du problème des moindres 15 carrés pondérés, ce qui ne change que l'expression finale de D où la moyenne devient une moyenne pondérée, mais pas celles de 'ai, , p =1...P. Comme vu précédemment, l'estimation de la position de l'émetteur est réalisée en deux temps. Une première estimation grossière est calculée par 20 la méthode de Bancroft. Le résultat obtenu est ensuite utilisée comme point d'initialisation pour une méthode du maximum de vraisemblance (également connue sous le sigle anglo saxon MLE pour " Maximum Likelihood Estimator"). Les expressions de ai, obtenues précédemment correspondent à des '.2 "TOA réalisations de variables aléatoires de moyenne 1:21-' + to et de variance La méthode de Bancroft consiste à ajuster les paramètres inconnus, à savoir la position (x,y,z) de l'émetteur et le temps d'émission to, aux valeurs de tip en négligeant le bruit, par l'intermédiaire de la relation liant les distances rop aux positions de l'émetteur et des P stations réceptrices. La relation entre rop, la position (x,y,z) de l'émetteur, et la position (xp,yp,zp) de la piè" station est : rop = .\/(x - xp)2 + (y - Yp )2 + (Z Zp )2 lo D'autre part, en absence de bruit, on a rop = c(cîp - to). En élevant les deux expressions de rop au carré, on obtient pour tout p=1...P : x?, +y72, +z?, - (ctip)2 +x2 +y2 + z2 - (ct0)2 = 2 (xxp yyp zzp - (ctip)(cto)). Cette dernière expression fait apparaître la fonctionnelle de Minkowski notée telle que pour tous vecteurs de dimension quatre v=(vi,v2,v3,v4) 15 et w=(wi,w2,w3,w4) on ait : (171w) = viwi + v2w2 + v3w3 - v4w4 - Si on place les informations liées à une station réceptrice p dans un vecteur sp et les inconnues dans un vecteur u: xP i F [ x I Yp Y s ;up = z = P Z C - ei, c . to En introduisant la matrice B suivante : Yi z1 -c - [xi Y2 z2 -c - ai B='2 B = 2 . Xp YP zp -c - tii; et en posant : (su s1) [(s2ls2) 1 k =1- 1 1 A = -2 (ulu) , e =[I, 2 : (sp Isp) 1 les P équations sur les inconnues s'écrivent de manière compacte : Bu = Ae + k. Le système ci-dessus est sur-dimensionné pour un nombre de stations réceptrices supérieur à 4 (P équations pour 4 inconnues). Ce système peut être résolu avec la pseudo-inverse B+ de la matrice B en supposant A connu dans un premier temps : 13+ = (B*B)-113* . u = B+ (k + Ae). En utilisant cette dernière expression de u pour développer l'expression de A= , on obtient un trinôme du second degré en A: A2(B+e,B+e)+ A(2(B+k,B+e)- 1) + (Bk, lek) = O. Ce dernier fournit deux solutions candidates pour A, notées A1 et A2, conduisant à deux solutions possibles ui et u2 au problème : = B+(k + Aie) ; u2 = B+(k + A2e) . lo Une de ces deux solutions sera cohérente par rapport à la situation. Par exemple, l'altitude estimée sera proche de la surface de la terre. Comme vu précédemment, la position trouvée par la méthode de 15 Bancroft est ensuite utilisée comme point d'initialisation pour une estimation du maximum de vraisemblance. Dans cette deuxième estimation, On prend maintenant en compte le bruit wn dans les mesures de temps d'arrivée des impulsions. L'expression des devient : 20 OÙ : c représente la vitesse de propagation des impulsions, (x, y, z) représente les coordonnées de la position de l'émetteur, (xp, yp, zp) représente les coordonnées de la position de la pième station réceptrice, to représente la date de début d'émission des impulsions, fp représente une fonction non linéaire qui relie la position de l'émetteur et la date de début d'émission. Cette fonction fp représente le temps d'arrivée de l'impulsion sans bruit au niveau de la pleine station, wn représente le bruit de mesure considéré comme un bruit 2 gaussien de variance cr rT ,?, = cN0A avec o-72.0A la variance instrumentale de la mesure des TOA et N le nombre d'impulsions périodiques considérés. Suivant un mode de mise en oeuvre de l'invention la méthode du maximum de vraisemblance peut contenir un a priori sur le paramètre to correspondant à la date d'émission de l'impulsion au centre du train d'impulsions émis par l'émetteur. On peut considérer, sans perte de généralité, que les impulsions sont datées en prenant comme origine la date de réception de l'impulsion centrale par la première station réceptrice (10) : CF17=0. Ainsi si l'on considère que les distances émetteur-stations réceptrices sont très supérieures aux distances entre les stations réceptrices (cas des bases courtes : par exemple la distance entre l'émetteur et la stations peut être de l'ordre de 200 km alors que la distance entre les stations peut être de l'ordre du km), on peut considérer que la distance do = - to) = -c to correspondant à une distance grossière entre les stations réceptrices et l'émetteur peut être modélisée par une variable aléatoire d'espérance notée do et de variance notée 40 relativement élevée. On choisit do proche de la limite de portée des récepteurs. Cela revient à éliminer des positions trop éloignées, qui seraient incohérentes car 30 dans ce cas le signal n'aurait pas été capté, et trop proches car on suppose que le signal aurait été capté auparavant. En pratique, même si l'émetteur est très proche des stations réceptrices, cet a priori n'introduit quasiment pas de biais car l'information qu'il apporte devient de plus en plus négligeable comparée à l'information apportée par les mesures réelles au fur et à mesure que l'émetteur se rapproche. Pour les P stations réceptrices on a le vecteur de mesures suivant : I di -1V, (xi - ,X)2 + (Yi - y)2+(zi - z) + to c f i(x, y, z, to) = 1 , + w = L p(x, y, z, to) + w -cV (xp - x)2 + (Y p - Y)2+(z p - z)2 + to --Croic to 43 La matrice de covariance liée à ces mesures (supposées indépendantes) est: 0 a 0 I E = ,r2 L' do / %.* Si on fait une approximation linéaire au point (xi, yi, zi, toi) pour les fonctions f p on obtient : X-Xi aP = fp(xbybzi t ) (afP afP afp lafp\ X - , 43X aY aZ C Ôt0) Yi I z zi c(to - avec (-"P),X1 xi) xi,yi,zi,toi C\i(Xp - xi)2 (yp - yi)2+(Zp - Zi)2 f9 f\ Yi Yp ay xbyi,zi,toi (a f\ zi - zp az ) xi,yi,zi,toi c,\1(xp - xi)2 (Yp - Yi)2+ (Zp - Zi)2 ( _1 f = ato ) xi,yi,zi,toi On considère une variable intermédiaire E définie par : 0 0 e = afi afi ()fi 1- xbyi,zi,to, x 15 - _ x - xi 1 c y - h - z - zi . c(to - toi) 1 _ c 1 - C- ax a:fp ay az 0:fp 0:1"p Ox Oy Oz 0 0 0 I di1 Ffi(xt,Yi,zt,toi)1\ - ap fp(x1,yi,zi,t03 \-do/c to 1 i = MA-R a toi ax '3Y az c(to - toi) T 0 0 0 1 Jj , A= z - Y Yi x - xi et N = f p yi, Zi, to I) Cli \-o/c xi,yi,zbtoi Cette variable matricielle E représente l'erreur d'approximation à l'itération courante. La première colonne représente l'erreur d'approximation en x, la deuxième en y, la troisième en z et la dernière en to. On minimise le critère c défini par : C = eTE-1e = (MA - R)TE-1(MA - R) Cela revient à minimiser la norme du vecteur d'erreur pondéré par la matrice E qui contient les variances des estimateurs ep. On obtient donc une nouvelle estimation de la position avec : A = (mTE-im)-1mTE-iR On réitère la procédure jusqu'à ce que A, le vecteur de déplacement vers la nouvelle position, soit inférieur à un seuil prédéterminé. Suivant un 15 mode de mis en oeuvre, on peut choisir comme contrainte, par exemple, 11All N = 30 Paramètres D = 200 km D = 60 km P =4 E-TOA : 4,86 km E-TDOA : 5,26 km E-TOA : 404 m E-TDOA : 412 m P = 8 E-TOA : 3,04 km E-TDOA : 3,35 km E-TOA : 236 m E-TDOA : 261 m Tableau 1 N = 3 Paramètres D = 200 km D = 60 km P = 4 E-TOA : 12,3 km E-TDOA : 16,7 km E-TOA : 1,28 km E-TDOA : 1,33 km P = 8 E-TOA : 9,21 km E-TDOA : 11,6 km E-TOA : 799 m E-TDOA : 873 km Tableau 2 La figure 5 illustre des exemples d'ellipses à 90% de confiance obtenues à l'aide du procédé selon l'invention notée E-TOA, à l'aide d'un estimateur utilisant la méthode du maximum de vraisemblance classique notée E-TDOA et avec la méthode de Bancroft notée E-BAN. Dans cet exemple, P=8 stations réceptrices, N=10 impulsions, et D=200 km. Le centre de la base contenant les stations réceptrices est situé en (0, 0). Les tableaux 1 et 2 ainsi que la figure 5 mettent en évidence la meilleure précision de localisation du procédé selon l'invention par rapport aux méthodes classiques. Comme vu précédemment, le fait de transmettre les temps moyen d'arrivée des impulsions plutôt que toutes les impulsions reçues permet de réduire la quantité de données échangées entre les stations réceptrices et permet donc une plus grande discrétion des stations. Le tableau ci-dessous présente le nombre de transfert de données 10 entre deux stations réceptrices pour N=3 et N=100 impulsions dans le cas du procédé selon l'invention (E-TOA) et celui d'un estimateur utilisant la méthode du maximum de vraisemblance classique (E-TDOA). Estimateur N = 3 N = 100 E-TOA 1 1 E-TDOA 3 99 Tableau 3 15 La présente invention a été décrite à travers une application aéroportée, cette description n'est nullement limitative et peut être généralisé à d'autres domaines comme par exemple le domaine spatial. 20 REVENDICATIONS1. Procédé de localisation passive d'un émetteur (15) non mobile au sol mis en oeuvre par un groupe d'au moins une station réceptrice (10), chacune des stations réceptrices (10) comprenant un radar passif et une référence de temps, l'ensemble des références de temps étant synchronisées entre elles, ledit émetteur (15) émettant un ensemble d'impulsions périodiques, ledit procédé étant caractérisé en ce qu'une première estimation de la position dudit émetteur (15) est réalisée par la méthode de Bancroft à partir des temps moyens d'arrivée des impulsions émise par l'émetteur (15) au niveau de chaque station (10) du groupe d'au moins une station réceptrice (10), le résultat obtenu étant ensuite utilisé comme point d'initialisation d'une méthode de maximum de vraisemblance afin de converger vers la position dudit émetteur. 2. Procédé selon la revendication précédente mis en oeuvre par au moins deux stations réceptrices (10) dans lequel ledit procédé comprend : Une étape Etp1 de mesure des temps d'arrivés des impulsions de l'émetteur (15) par chaque station réceptrice (10), - Une étape Etp2 de calcul du temps moyen d'arrivée des impulsions par chaque station réceptrice (10), - Une étape Etp3 de transmission des temps moyens d'arrivée des impulsions calculés à une station réceptrice (10) prédéterminée, Une étape Etp4 de calcul, par ladite station réceptrice (10) prédéterminée, de la localisation de l'émetteur (15) à partir de tous les temps moyens d'arrivée des impulsions à l'aide de la méthode de Bancroft puis de la méthode du maximum de vraisemblance. 3. Procédé selon une des revendications précédentes dans lequel la ou les stations réceptrices (10) sont des stations réceptrices aéroportées. 4. Procédé selon une des revendications 1 ou 2 dans lequel la ou les stations réceptrices (10) sont des stations réceptrice spatiales.