FR 3046326 A1 20170630 FR 1502678 A 20151223 1 PROCEDE ET SYSTEME POUR LA DETECTION DE SIGNAUX DANS UN RESEAU MULTI SOURCES L'invention concerne un procédé et un système permettant 5 notamment de détecter des signaux issus de sources par exemple dans un réseau de communication. L'invention s'applique aussi dans le domaine de la géolocalisation. L'invention concerne notamment le traitement d'antennes qui permet d'analyser ou d'estimer les signaux des sources émettrices à partir 10 du traitement des signaux en sortie des capteurs (ou antennes) des différentes stations. Dans un contexte électromagnétique, les capteurs sont des antennes et les sources radioélectriques se propagent avec une polarisation dépendante de l'antenne d'émission. Dans un contexte acoustique, les capteurs sont des microphones et les sources (ou émetteurs) 15 sont sonores. La figure 1 montre que la i-ième station 10; est composée d'un réseau de capteurs élémentaires 11k; où les différentes sources 12q (ou émetteurs) arrivent suivant des angles Oqi qui sont différents, avec q l'indice de la source (figure 2). Ces angles se nomment aussi en abrégé anglo-saxon AoA (Angle of Arrivai). Le temps de propagation entre une source et une 20 station se nomme en abrégé anglo-saxon ToA (Time of Arrivai). Pour une source arrivant sur deux stations, la différence entre les temps de propagation ToAs de deux stations se nomme en abrégé anglo-saxon TDoA (Time Difference of Arrivai). L'un des objectifs des techniques de traitement d'antennes est d'exploiter la diversité spatiale d'un réseau multi-antennes 25 recevant les signaux de plusieurs sources. Le réseau d'une station permet alors de différencier les sources par un contraste spatial qui dépend de l'angle d'arrivée (AoA) ou encore par leur polarisation incidente. Le réseau global de géolocalisation composé des antennes de toutes les stations permet, d'un point de vue traitement d'antennes, de différencier les sources 30 par leurs temps d'arrivée (ToA), en plus des angles d'arrivée AoA sur chacune des stations. 3046326 2 L'estimation paramétrique permet de déterminer les paramètres des sources : les coordonnées de la position, les angles d'arrivée AoAs ou les différences de temps d'arrivée TDoAs des sources incidentes. Les AoAs et TDoAs sont dépendants des positions des sources et des stations. Une application particulière de l'estimation paramétrique est la goniométrie dont l'objectif est de déterminer les angles d'arrivée (ou AoAs) des sources sur chacune des stations, à partir des signaux observés en sortie des capteurs du réseau de la station. L'invention concerne aussi la géolocalisation de sources en exploitant les signaux reçus en sortie des antennes 11k; d'une station 10i multi-capteurs (ou multi-antennes). Les techniques de géolocalisation permettent d'estimer la position géographique paramétrée par les vecteurs de coordonnées pl, p2, pq des sources 121, 122, 12q, en x, y et z dans l'espace (repère cartésien) ou en x, y dans le plan horizontal 13. Les signaux transmis par les émetteurs (ou les sources) sont inconnus des systèmes de réception des stations qui sont généralement des réseaux multi-antennes. Les figures 2 et 3 définissent l'angle d'incidence 0 d'un émetteur de type mobile 12 arrivant sur les capteurs d'une station 10i. Les angles d'incidence sont paramétrés en une dimension 1D par l'azimut cp et en deux dimensions 2D par le vecteur 0 dont les composantes sont l'azimut cp et l'élévation A. La figure 3 montre par ailleurs que l'angle d'incidence d'une source 12 dépend à la fois de la position b de la station 10; et aussi de celle de la position de la source 12p qui est notée par p. Dans le cadre des techniques de traitement d'antennes, les méthodes à haute résolution de type MUSIC, par exemple la méthode décrite dans le document de RO.SCHMIDT intitulé « Multiple emitter location and Signal Parameters Estimation », IEEE Trans on Antennas Propagation, vol 34, pp276-280, 1986, permettent d'estimer des paramètres de type AoAs (angles d'arrivée) à partir des signaux en sortie de plusieurs capteurs. Ces méthodes ont l'avantage d'être non biaisées en présence de plusieurs sources. Ces techniques ont été adaptées dans le cadre de la géolocalisation en utilisant 3046326 3 les signaux de toutes les antennes de chacune des stations. Cela permet d'utiliser conjointement les angles d'arrivée (AoAs) des sources sur chacune des stations ainsi que les temps d'arrivée (ToA) qui sont porteurs d'une information de distance entre les stations et les sources. 5 L'un des problèmes techniques de la géolocalisation de sources à partir d'un ensemble de stations est celui de la transmission d'informations entre les stations qui est nécessaire pour déterminer la position des sources. Cette information est aussi bien un paramètre (AoA ou un ToAs) que les signaux temporels reçus sur les stations. L'objectif est de transmettre le 10 minimum de quantité d'information tout en garantissant des performances quasi-optimales. En effet, les paramètres estimables au niveau d'une station ne permettent pas de déterminer la position d'une source. En particulier, à partir d'une seule station, il est possible de mesurer la direction (ou l'AoA) d'une source. Toutefois cette information angulaire n'est pas suffisante pour 15 estimer la position. Pour cela, il faudrait connaître la distance entre l'émetteur et la station en plus de sa direction d'arrivée. Pour obtenir la position de la source, il est donc nécessaire d'obtenir des mesures supplémentaires à l'aide d'au moins une autre station : dans les systèmes de géolocalisation par triangulation, les angles d'arrivée sont mesurés par au moins deux stations. 20 Dans le cadre d'un système de télécommunication, le temps d'arrivée ToA est estimable à partir de la connaissance des signaux pilotes. Cela permet d'avoir accès à une information qui permettrait d'estimer une distance émetteur-récepteur. Toutefois, comme l'instant d'origine de transmission est généralement inconnu sur une station, il est nécessaire de mesurer les 25 temps d'arrivée ToAs ou les angles d'arrivée AoAs sur d'autres stations pour lever l'ambiguïté de distance et déterminer la position des sources. Deux grandes approches ont été envisagées dans le domaine de la géolocalisation avec un système composé de plusieurs stations multi ou mono-antenne(s) : 30 - Les techniques dites en deux étapes sont illustrées sur la figure 4. Cette approche est celle qui est utilisée dans la plupart des systèmes 3046326 4 de géolocalisation. Dans la première étape 41, un ensemble des paramètres de type AoAs( ê,,) ou ToAs( ) est estimé sur chacune des stations 101,1=1.. L pour des sources 12q, q=1..Q, l'estimation des ToAs est effectuée à partir de la connaissance d'un signal pilote. Les 5 paramètres de type TDoA, correspondant à la différence de temps d'arrivée entre deux stations, nécessitent la transmission d'un signal entre les deux stations. Les paramètres (AoA, ToA) estimés sur les différentes stations [{Ôq/},{i-sq/},...1 sont alors transmis 42 à une station centralisée 43 où la seconde étape de détermination de la 10 position des sources 44, 13, = ÿ, 2)T , est appliquée en utilisant les paramètres estimés lors de la première étape. Les techniques dites en une étape, par exemple comme celle décrite dans le document J.Bosse, A.Ferreol et P.Larzabal, A Space Time Array Processing for Passive Geolocalization of Radio Transmitters, 15 IEEE ICASSP 2011-Prague, pp 2596-2509, sont illustrées sur la figure 5. Cette approche nécessite tout d'abord la transmission des signaux multi-antennes 51 reçus sur chacune des stations 101 vers un module (t) de centralisation des signaux 52, qui reçoit un signal y(t)- ye,(t) où les vecteurs xi(t) sont composés des signaux reçus sur les 20 antennes de la i-ième station. Des techniques à haute résolution (HR) de type MUSIC peuvent alors être appliquées pour déterminer directement, module d'estimation des positions, 53 la position des sources,pq5,,2)T pour les Q sources. Les deux approches connues de l'art antérieur ont chacune leurs 25 propres inconvénients, en particulier : - Les techniques en deux étapes comme la triangulation (AoAs) ou encore les méthodes TDoAs dites hyperboliques sont souvent 3046326 5 ambiguës en contexte mufti-sources. Pour les méthodes AoAs (ou triangulation), l'information de corrélation ou non des signaux entre stations n'est pas utilisée. Les méthodes de types TDoAs ont l'inconvénient d'être biaisées ou non résolues en multi-sources. En 5 résumé, les paramètres en sortie de la première étape donnent généralement une information incomplète des sources à localiser. Les performances d'estimation de la position à la seconde étape sont fortement dépendantes de celles de l'estimation des paramètres AoAs et ou TDoAs de la première étape. Cela a pour conséquence de 10 limiter les performances de géolocalisation soit en termes de précision, soit en termes de résolution. - Les techniques de géolocalisation en une étape ont notamment l'inconvénient de nécessiter la transmission de tous les signaux des capteurs (ou antennes) de chacune des stations vers une station 15 centralisée. Toutefois, les méthodes TDoA qui sont en deux étapes nécessitent aussi une transmission de signaux entre deux stations de base. Cependant, la quantité de signaux à transmettre est moindre que les techniques en une étape, car la méthode traditionnelle qui consiste en une mesure des maxima de la fonction d'intercorrélation 20 entre un seul signal de chaque base (on inter corrèle les signaux de deux capteurs éloignés), requiert la transmission d'un seul signal de l'une des deux stations vers l'autre station. L'un des problèmes non résolu par l'art antérieur est, notamment, celui de la mise en oeuvre des méthodes de géolocalisation en une étape, 25 sans avoir besoin de transmettre tous les signaux en sortie de toutes les antennes de chacune des stations. De manière plus générale, ce problème se trouve aussi dans les procédés de détection de signaux issus de plusieurs sources dans un réseau de communication. Dans la suite de la description le mot « projection » désigne la 30 même étape que le mot « compression » et le mot « signal » le même objet que le mot « observation ». 3046326 6 L'un des objectifs du procédé selon l'invention est, notamment, de diminuer le nombre de signaux à transmettre sans perdre les informations et les paramètres nécessaires pour la reconstruction de signaux. Ainsi, en contexte mono-source, le nombre de signaux transmis est égal à celui la 5 technique traditionnelle de géolocalisation de type TDoA. L'invention concerne un procédé de détection de signaux issus d'une ou de Q sources dans un réseau de communication comprenant au moins deux stations, une station étant équipée de plusieurs capteurs élémentaires ou antennes de réception en liaison avec un processeur de 10 traitement de données, comportant en combinaison au moins les étapes suivantes : - Séparer à partir des observations reçues sur une station, les signaux des Q sources avec le signal du bruit, - Déterminer un espace signal ES' de dimension égal au nombre de 15 sources pour une station, - Projeter l'observation reçue sur une station I, x/ (t) de dimension NI en utilisant l'espace signal E5' afin d'obtenir une observation xl (t)=Eswx/ (t) de dimension Qxl, - Transmettre le signal projeté x, (t) = E:Hx/ (t) ainsi que le projecteur à 20 un module de reconstruction de signal afin d'obtenir un signal y/ (t) = Esiii(t) de dimension Nixl Selon une variante de réalisation, le réseau comporte plusieurs stations locales et une station centrale, l'étape de projection i/ (t)=Esufx/ (t) est exécutée au niveau d'une station locale I et l'étape de reconstruction du 25 signal y/ (t) = E515e/ (t) est exécutée au niveau d'un module de reconstruction de la station centrale. L'étape de projection du signal reçu sur l'espace signal (t)=Esnixi(t) et l'étape de transmission du projecteur sont, par exemple, exécutées au niveau d'une station locale. 3046326 7 L'étape de séparation du signal utile et du bruit peut être réalisée en exécutant au moins les étapes suivantes : - Calcul de la matrice de covariance du signal RX = Elx/ (t)x/H (t)] , - Décomposition en valeur singulière de IK pour obtenir les valeurs 5 propres associées aux vecteurs propres ell. Détermination du nombre Q de sources avec un algorithme connu utilisant les valeurs propres , - Construction de la matrice de compression [eu eQ,] " Lorsque l'on met en oeuvre une méthode TARGET connue de 10 l'homme du métier, on détermine en plus la matrice bruit ainsi que de la valeur propre de bruit 62 correspondant à la moyenne des N, -Qplus faibles valeurs propres de Rix et on transmet en plus le niveau de bruit 62 et les valeurs propres signal .11, - - - /10 pour la reconstruction du signal. Selon une variante de réalisation le procédé comporte au moins 15 les étapes suivantes : - Acquisition des signaux de la seconde station x2(t,) avec un échantillonnage tel que ti=ixT, où T, est supérieur et proche de 1/B , où B est la bande des récepteurs, - Transfert des signaux x2 (t,) de la seconde station vers la première 20 station selon les étapes du procédé pour recevoir le signal x2 (t), la matrice diagonale g de l'équation Q oftx/) = 2 eqieH=E:gEsni et l'espace signal q=1 2q/ - 6 E52 = [eu eQ2] sur la première station, - Sur-échantillonnage d'un facteur K des signaux x2 (t,) et xi (t, ) pour 25 obtenir les signaux i2(i,) et x1(ii) avec i,=Kxt,, 3046326 8 5 - Calcul des matrices 1 k12 (r) T (t)i2H (t-r) et ft21(-r) -TEi2 (t (t) à partir des t=i signaux x2 (t1)et xi (Ti), - Calcul des matrices G, (r) et G2 (l-) selon G1 (r)=Im -(Rxi)+ 1Z12 (r)S121221 (-r) G2 (r) = Im2 -Es2n2ft21 (-T)(Rxi (r)E52H - Calcul et optimisation du critère 11 in a, JTARGET (1919 r) = (6s"r)-4- 1.(e 2'r) avec *Ir' rj 0)(0)a, (0)(8) 2 02 = f (81,r) pour obtenir les couples AOA-TDOA (01(pq),1-2(pq)--r1 (pq)) de chacune des sources, 10 - Déduction de la position pq de chacune des sources à partir de (01 (Pq ), T2 (Pg Z-1 (Pq )) et en utilisant r1 (pq = HPq k(01 (p )) = Pq Ipq -bill La station centrale choisie est, par exemple, la première station du réseau global et le procédé comporte au moins les étapes suivantes : 15 - Pour 1=2 jusqu'à L, nombre de stations dans le réseau de communications, - Transfert des signaux de la /-ième station vers la première station selon les étapes du procédé de l'invention afin d'obtenir le signal y/ (t) sur la première station, 20 - Construction de l'observation y (t) multi-stations avec y (t) = 3046326 9 - Application d'une méthode en une étape pour obtenir les positions pq des sources incidentes : choisie parmi la liste suivante DPD, spatio-temporelle, TDOA. A partir du signal y1 (t) = (t) de dimension Nix1 on réalise une 5 géolocalisation de la source émetteur. L'invention concerne aussi un système pour détecter des signaux issus d'une ou de Q sources dans un réseau de communication comprenant au moins deux stations, une station étant équipée de plusieurs capteurs élémentaires ou antennes de réception en liaison avec un processeur de 10 traitement de données, comportant au moins les éléments suivants : - Une station est équipée d'un module adapté pour déterminer un espace signal E) de dimension égal au nombre de sources pour une station, et à projeter l'observation reçue sur une station I , x/(t) de dimension N, en utilisant l'espace signal ES' afin d'obtenir une observation 15 x, (t)=Esmxi(t) de dimension Qx1, - Un module de transmission des observations projetées, - Transmettre le signal projeté i/(t)=Esulyel(t) ainsi que le projecteur à un module de reconstruction de signal afin d'obtenir un signal (t) = E51i1(t) de dimension N,x1. 20 Le système peut comporter un dispositif de traitement d'antennes en une étape recevant les signaux yi (t) (t) de dimension Nix1 pour les L stations. D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description d'exemples donnés à titre 25 illustratif et nullement limitatifs annexée des figures qui représentent : - Figure 1, une illustration de contexte d'une réception de deux sources avec L stations éloignées, - Figure 2, une représentation de la propagation de plusieurs sources en LoS vers une station multi-capteurs, 3046326 10 - Figure 3, la représentation de l'incidence 0-(ço, A) d'une source sur un réseau de capteur où ço est l'azimut et A l'élévation, - Figure 4, une représentation des méthodes de géolocalisation en deux étapes, triangulation, temps d'arrivée ou ToA (time of arrivai), fréquence 5 d'arrivée ou FoA (Frequence of Arrivai) avec transmission des paramètres intermédiaires, - Figure 5, une illustration des méthodes de géolocalisation en une étape, avec transmission des signaux bruts des stations, - Figure 6, une représentation d'un exemple d'étapes de la méthode de 10 géolocalisation en une étape selon l'invention, avec transmission des signaux compressés des stations, le nombre des signaux transmis par station étant au plus égal au nombre de sources à géo-localiser, - Figure 7, une illustration d'une méthode de géolocalisation en une étape, cas du DPD avec une décomposition du signal en K sous- 15 canaux, - Figure 8, une représentation d'une méthode de géolocalisation en une étape, cas du LOST/LOST-FIND avec la construction d'une observation spatio-fréquentielle, et - Figure 9, une illustration d'une méthode de géolocalisation en une 20 étape, cas de la méthode large bande TARGET. Le procédé mis en oeuvre par l'invention consiste notamment à projeter une observation multi-antennes sur un sous-espace signal obtenu par une technique de décomposition en éléments propres telle que celle décrite dans la méthode MUSIC connue de l'homme du métier dont la 25 dimension est égale au nombre de sources ou en utilisant une technique équivalente. Pour mettre en oeuvre le procédé selon l'invention, figure 6, une station locale à un réseau de télécommunication est équipée de modules adaptés. Le contexte est un réseau global 10 comprenant plusieurs stations 30 101. Une station 101 comporte plusieurs antennes 10k1 de réception et une antenne d'émission/réception 6011 ayant pour fonction de recevoir et/ou de 3046326 11 transmettre les signaux traités par l'étape de compression ou de décompression, par exemple, un module de traitement 6021 adapté à exécuter, par exemple, les étapes du procédé selon l'invention, et un module de communication 605. Le module de communication 605 assure la 5 transmission des signaux i(t) et la communication au sein du réseau entre les différentes stations. Les signaux sont ensuite reçus sur une station qui comprend par exemple un module 67 de décompression des signaux, suivi d'un module de centralisation des signaux et d'estimation de la position, 68, correspondant par exemple à un algorithme de traitement d'antennes en une 10 étape. Dans les exemples exposés ci-après les données sont transmises à la première station qui est considérée comme la station centrale au réseau global. Sans sortir du cadre de l'invention, on pourrait envisager de transmettre les données en provenance des stations vers une station extérieure au réseau qui comporterait un module de traitement d'antennes en 15 une étape, par exemple. Dans l'exemple qui suit, un opérateur est considéré comme une station. Une station centrale, par exemple la station 101, est désignée au démarrage comme étant une station qui va recevoir tous les signaux qui ont été ou sont reçus sur les autres stations du réseau global. La station centrale 20 sera, par exemple, la station du réseau de communications qui comporte le plus grand nombre d'antennes. La station centrale étant déterminée, le procédé va exécuter au niveau de chacune des stations locales du réseau, une projection des signaux reçus (observations) sur un espace signal propre à une station (station qui a reçu ces observations), puis une étape de 25 transmission vers la station centrale des signaux projetés de chacune des stations, ainsi que le projecteur utilisé et calculé à partir des signaux de la station. La station centrale sera alors en mesure de reconstruire le signal correspondant à une observation multi-antennes. Quelques rappels sur la représentation mathématique du signal en 30 sortie d'une station multi-antennes sont donnés. Le signal sur la 1-ième station peut être représenté de la manière 3046326 12 suivante, en présence de Q sources de signal sq (t), le signal reçu en sortie de la /-ième station s'écrit : x, (t) = x,,'(t) où x, (t) est un vecteur signal de dimension /V, tel que xh (t) est le signal en sortie de la i-ième antenne de la /-ième station. Le vecteur n, West le bruit 5 additif qui est supposé spatialement blanc. Le signal sq, (t) de la q-ième source arrive sur la /-ième station avec une atténuation complexe ,o,q et un ToA rl(pq) vérifiant : 11Pq c-1)1 (P q)= où b, est le vecteur position de la /-ième station, c la célérité de la lumière et ull la norme du vecteur u . Dans l'équation (1), le vecteur a, (0) est la lo réponse du réseau de capteurs de la /-ième station à une direction A , le signe - est une simplification dans l'écriture indiquant que ce vecteur dépend de la source indicée par q et de la station indicée par I. La direction est l'AoA de la q-ième source incidente à la l-ième station vérifiant : et (Pq) (3) . k(e,(p )) = pq Dpq où k(0) est le vecteur d'onde de la figure 3, et cos Mcos (à) k(e)= sin(ço)cos(A) sin (A) 15 où ço est l'azimut et A l'élévation. En présence de L stations, les méthodes précitées de géolocalisation en une étape nécessitent la construction de l'observation multi station suivante : âq, = (0, (pq)) =EQ àq,sq, (t)+ ni (t) où S (t) = PiqS q(t - r (p q)) (2) . avec 0=(gp, A) (4) . 3046326 13 x(t) = (5) xL (t) _ _ qui requiert les signaux reçus sur ou par toutes les stations du système ou réseau global 10. En supposant que la première station 101 soit la station centrale, il est nécessaire de transmettre les vecteurs x2 (t), xL (t) reçus par les stations correspondant à N2 + NI, signaux à transmettre pour 5 appliquer les méthodes en une étape. Dans cet exemple, l'objectif est de transmettre, au maximum, un nombre de signaux valant QxL, nombre qui est inférieur au nombre N2 ± NL, en présence de Q sources. D'après l'équation (1), l'observation du signal sur la /-ième station peut se ré-écrire de la manière suivante : À, =là,/ ào] xi(t)= Àisi(t)+ni(t) où si/ (t) s (t) = s (t) 10 où À1 est une matrice définissant l'espace signal des sources sur la 1-ième station. Cette matrice À/ caractérise entièrement les caractéristiques en termes d'AOAs des sources incidente à cette station. En supposant que le bruit additif soit blanc on sait que la matrice de covariance de n, (t)vérifie = E [ni (t)niii (t)] = cr2IN, où .H désigne le transposé conjugué , E[..] 15 l'espérance mathématique et IN, la matrice identité de dimension /V/ xN, Dans ces conditions la matrice de covariance de l'observation x/(t) s'écrit : le' (t)]= ÀiRsi + cr2IN avec Rsi = E[s (t)s111 (t)]. En présence de Q sources non cohérentes la décomposition en valeurs singulières de Ri. s'écrit : (6) . (7) 3046326 14 Q( Rxl = E 621-1 où H = ee q=1 \i=Q+1 où H est le projecteur bruit de l'algorithme MUSIC et où les valeurs 9, >0-2 (bruit) sont les valeurs propres de l'espace signal associé au vecteur propre eql " D'après le principe de la décomposition en éléments propres exploité par MUSIC, la matrice Es1= [eu - - - e de l'espace signal vérifie : À1 = avec fi, = Es1/1À, où la matrice fi/ est de dimension QxQ. 5 Cela veut dire que les colonnes de la matrice Es' = [eu ' - eol engendrent le même espace que les colonnes de la matrice À. . Cet espace est l'espace signal des sources qui contient toute l'information utile. Il est alors possible d'envisager de projeter l'observation x1 (t) sur la matrice signal Es pour obtenir l'observation projetée : î, (t)=Es'xi (t)= 13/ (t) où (t)=Esilini (t) 10 où x, (t) est une observation compressée (projetée) de dimension Qx1 inférieure à la dimension de x1(t) qui est de dimension Nix1. Selon l'équation (9) il est possible de construire une observation y i(t) de dimension Nix1 proche de xl (t), à partir du signal compressé (projeté) xl (t) ainsi que de la matrice de compression (projection) Es1 . L'observation y1(t) décompressée 15 (reconstruite en utilisant le projecteur et en réalisant l'opération inverse de la projection sur l'espace signal) s'écrit alors : y1 W.=Es1xr(t) = Àis, (t)+ à, (t) où à, (t)(t) = Es'Esuin, (t) Les étapes de compression et décompression du signal sur la 1-ième station du procédé selon l'invention sont alors les suivantes : Etape-1.1 : Calcul de la matrice de covariance du signal Ri. = E[x, (t)x/11 (t)] , 20 Etape-1.2 : Décomposition en valeur singulière de Rix pour obtenir les (8) (9) (10) 3046326 15 valeurs propres À, associées aux vecteurs propres ed . Détermination du nombre Q de sources avec un algorithme connu de l'homme du métier utilisant les valeurs propres 11,d , Etape-1.3 : Construction de la matrice de compression ou de projection 5Esi = [eu eo] ainsi que de la valeur propre de bruit o2 lorsque la méthode est TARGET qui est la moyenne des N, -Q plus faible valeurs propres de RI. , Etape-1.4 : Compression du signal en effectuant l'opération X, (t) = Esmx/ (t), 602i, 10 Etape-1.5 : Transmission, 605, du signal X1(t), ainsi que de la matrice de compression Es vers la première station 101 qui est la station centralisée. Pour la méthode TARGET décrite à la figure 9, on transmet, en plus, le niveau de bruit 62 et les valeurs propres signal .1.0 - , Etape-1.6 : Décompression du signal sur la première station 101 en 15 effectuant : y1(t)=Esliel(t), 67, module de décompression. L'étape suivante consiste à exécuter un algorithme connu de l'homme du métier, permettant d'estimer les positions des sources, module de centralisation et d'estimation des positions des signaux 68. Dans cet exemple, la station centrale est équipée d'un algorithme de géolocalisation 20 de traitement d'antennes en une étape. Sans sortir du cadre de l'invention, il serait aussi possible de transmettre les différentes données à un module de traitement d'antennes en une étape situé sur une station extérieure au réseau. En supposant que la première station 101 est la station centrale, 25 les méthodes de géolocalisation en une étape, avec la transmission décompressée des signaux des stations vers une station centrale, sont alors applicables avec l'observation multi-station suivante : 3046326 16 y(t)= (12) y, (t) Les étapes permettant d'appliquer les méthodes en une étape, par exemple, la méthode de détermination directe de position DPD (Direct Position Determination) décrite à la figure 7 et la famille spatio-temporelle : LOST et LOST-FIND décrite à la figure 8 avec une compression du signal à partir de L 5 stations éloignées sont par exemple décrites ci-après : Etape-2.1 : Initialisation du processus à /=2, Etape-2.2 : Transfert des signaux de la /-ième station vers la première station 101 selon les étapes 1.1 à 1.6 du procédé selon l'invention pour obtenir le signal 371(t) sur la première station, puis incrémenter 1=1+1, 10 Etape-2.3 : Si / REVENDICATIONS1 - Procédé de détection de signaux issus d'une ou de Q sources dans un réseau de communication (10) comprenant au moins deux stations (101), une station (101) étant équipée de plusieurs capteurs élémentaires (11kI) ou antennes de réception en liaison avec un processeur de traitement de données, comportant en combinaison au moins les étapes suivantes : - Séparer à partir des observations reçues sur une station (101), les signaux des Q sources avec le signal du bruit, - Déterminer un espace signal Es1 de dimension égal au nombre de sources pour une station (101), - Projeter l'observation reçue sur une station I , x, (t) de dimension N1 en utilisant l'espace signal Es1 afin d'obtenir une observation (t)=Esrxxi (t) de dimension Qx1, (6021), - Transmettre le signal projeté = (t) (t) ainsi que le projecteur à un module de reconstruction de signal afin d'obtenir un signal y, (t) = (t) de dimension N1x1. 2 - Procédé selon la revendication 1 caractérisé en ce que le réseau comporte plusieurs stations locales (102, .., 10L) et une station centrale (101) et en ce que l'étape de projection X, (t) = ES`Hx,(t) est exécutée au niveau d'une station locale I et l'étape de reconstruction du signal y1 (t) = (t) est exécutée au niveau d'un module (67)-(6021) de la station centrale. 3 - Procédé selon la revendication 1 caractérisé en ce que l'étape de projection du signal reçu sur l'espace signal i/(t)=Esuixi(t) et l'étape de transmission du projecteur sont exécutées au niveau d'une station locale. 4 - Procédé selon l'une des revendications 1 à 3 caractérisé en ce que 3046326 21 l'étape de séparation du signal utile et du bruit est réalisée en exécutant au moins les étapes suivantes : - Calcul de la matrice de covariance du signal Rix = E[xi(t)x1H (t)1, - Décomposition en valeur singulière de lex pour obtenir les valeurs 5 propres .1.,/ associées aux vecteurs propres eu. Détermination du nombre Q de sources avec un algorithme connu utilisant les valeurs propres /1d , - Construction de la matrice de compression Es/ = [eu eQ11 " 10 5 - Procédé selon la revendication 4 caractérisé en ce que pour une méthode TARGET on détermine en plus la matrice bruit ainsi que de la valeur propre de bruit 62 qui est la moyenne des Ni -Q plus faibles valeurs propres de Rix et on transmet en plus le niveau de bruit 62 et les valeurs propres signal .1v - - - .1.0 pour la reconstruction du signal. 15 6 - Procédé selon la revendication 5 caractérisé en ce que l'on exécute au moins les étapes suivantes : Acquisition des signaux de la seconde station x2 (t,) avec un échantillonnage tel que t,=ixT, où Te est supérieur et proche de 11B , où B est la bande des 20 récepteurs, Transfert des signaux x2(t, ) de la seconde station vers la première station selon les étapes du procédé pour recevoir le signal i2(t), la matrice fiz Q 1 diagonale f22 de l'équation i') .E, _ _2 eq/eq/H = EsigEsm et l'espace q=1 ql signal Est = [e12 " eQ2 sur la première station, 25 Sur-échantillonnage d'un facteur K des signaux x2 (t,) et x1 (t,) pour obtenir 3046326 22 les signaux z2 (0 et xi (0 avec Ti =Kxti , Calcul des matrices ft12(r)= 1 -Ex (t)% (t-r) et 1 2 (t-r)x,H(t) à T t=1 t=i partir des signaux i2(0 et x1(0 , Calcul des matrices G1 (r) et G2 (r) selon 5 G, (r) =IM--( iti2 (7)ÇI2R21 (-1") G2 (r) = - Es21-22i121(-r)(Rxi it12 (r)Es2H Calcul et optimisation du critère (8 ,r)+ JT,(02,r) JTARGET '1-\ 7; 2 avec les couples AOA-TDOA (01 (pj JT in \ (8)G i(r)a,(0) pour obtenir j aiH (0)a,(0) 02 = f (81,r) (pq)) de chacune des sources, Déduction de la position pq de chacune des sources à partir de 10 (01 (p,), /-2 (pj - (p,)) et en utilisant ri (pq) OPq -b/11 k (0, )) - P q b cI), -bill 7 - Procédé selon l'une des revendications 1 à 5 caractérisé en ce que la station centrale est la première station du réseau global (10) et en ce qu'il comporte au moins les étapes suivantes : 15 - Pour 1=2 jusqu'à L, nombre de stations dans le réseau de communications, - Transfert des signaux de la /-ième station vers la première station selon les étapes du procédé selon l'une des revendications 1 à 5 afin d'obtenir le signal y, (t) sur la première station, 20 - Construction de l'observation y(t) multi-stations avec : y(t) x, (t) Y2 (t) Y L(t) 3046326 23 - Application d'une méthode en une étape pour obtenir les positions pq des sources incidentes : choisie parmi la liste suivante DPD, spatio-temporelle, TDOA. 5 8 - Procédé selon l'une des revendications 1 à 7 caractérisé en ce qu'à partir du signal yi (t) = (t) de dimension Nix1 on réalise une géolocalisation de la source émetteur. 9 - Système pour détecter des signaux issus d'une ou de Q sources dans un 10 réseau de communication (10) comprenant au moins deux stations (101), une station (101) étant équipée de plusieurs capteurs élémentaires (11k1) ou antennes de réception en liaison avec un processeur de traitement de données, comportant au moins les éléments suivants : - Une station est équipée d'un module adapté pour déterminer un espace 15 signal E5' de dimension égal au nombre de sources pour une station (101), et à projeter l'observation reçue sur une station I (t) de dimension N, en utilisant l'espace signal Es! afin d'obtenir une observation X, (t) = Esmxi (t) de dimension Qx1, - Un module de transmission (605) des observations projetées, 20 - Transmettre le signal projeté z, (t) = Esnixi (t) ainsi que le projecteur à un module (67)(6021) de reconstruction de signal afin d'obtenir un signal y1 (t) = (t) de dimension Nix1. 10 - Système selon la revendication 9 caractérisé en ce qu'il comporte un 25 dispositif de traitement d'antennes en une étape recevant les signaux (t) = Esii/(t) de dimension Nix1 pour les L stations.