FR 3153424 A1 20250328 FR 2310006 A 20230921 L’invention concerne un procédé de détermination de la direction d’arrivée des trajets de signaux ionosphériques Hautes Fréquences sur un réseau (2) de N capteurs (3), comprenant: estimation, par goniométrie parcimonieuse, du vecteur vérifiant où σ 2 est le niveau du bruit et R x est la matrice de covariance vectorisée du signal x( t) considéré en sortie du réseau de capteurs, où le vecteur a (θ,Δ) de dimension N x1 est la réponse, connue, du réseau de capteurs à un signal reçu de direction d’arrivée (θ,Δ), déduction, en fonction de l’indice q m de chaque composante non nulle dudit estimé, du gisement de trajets depuis un émetteur respectif m ; pour chaque gisement déduit, étant donné que , les composantes non nulles de type desont ensuite déterminées : l’élévation d’indice q p associée à chacune de ces composantes non nulle , est sélectionnée comme égale à une élévation d’un trajet du m ième émetteur. Figure pour l’abrégé : Fig. 1 Procédé de détermination de la direction d’arrivée des trajets de signaux ionosphériques HF sur un réseau de N capteurs, programme d’ordinateur et dispositif associés L’invention se situe dans le domaine de la goniométrie, qui est la détermination de la direction d'arrivée de signaux électromagnétiques (radio-électriques ou sonores). Elle concerne plus spécifiquement l’estimation de la direction d’arrivée de signaux incidents sur un réseau de plusieurs capteurs (quand les signaux sont radio-électriques par exemple, les capteurs sont des antennes) dont les caractéristiques sont connues. Les signaux incidents peuvent provenir d’un ou plusieurs émetteurs. Pour un même émetteur, le signal émis peut être reçu selon plusieurs trajets de propagation, donnant lieu à plusieurs signaux incidents sur le réseau de capteurs. L’invention vise à réaliser cette estimation y compris quand les signaux incidents ont la caractéristique d’être temporellement corrélés (ce qui peut être le cas lorsque les signaux sont des multi-trajets du même émetteur) et que les émetteurs, par exemple en grand nombre, génèrent au total un nombre de trajets potentiellement supérieur au nombre de capteurs du réseau de capteurs. Les situations de multi-trajets donnant lieu à la présence de trajets corrélés sont rencontrées par exemple en milieu urbain dans les réseaux GSM (Global System for Mobile Communications) ou LTE (Long Term Evolution), du fait des réflexions sur des bâtiments, ou encore dans le cadre de communications ionosphériques HF (Hautes Fréquences). On s’intéresse plus spécifiquement à ces communications HF. Pour les multi-trajets corrélés HF, les signaux de chacun des émetteurs se propagent suivant plusieurs trajets se réfléchissant sur les couches ionosphériques (en un point R1 sur la couche E et en un point R2 sur la couche F pour l’émetteur), comme illustré en FIG. 2 : les différents trajets d’un émetteur de position E arrivent généralement sur un récepteur R avec le même gisement θ m et des angles d’élévations Δ mp différents. En FIG. 2 , le nombre de trajets représentés entre l’émetteur E et le récepteur R est égal à 2 : p =1,2. Des trajets arrivent au récepteur R avec des temps de propagation différents qui sont proportionnels à la distance parcourue. On aura ainsi des trajets corrélés lorsque la différence de temps de propagation entre deux trajets est inversement proportionnelle à la bande passante de l’émetteur. C’est ainsi que les trajets depuis un émetteur en HF sont souvent corrélés entre eux. D’autre part, les techniques de goniométrie classiques nécessitent de réaliser des intercepteurs avec des sous-bandes relativement faibles de manière à ce que le nombre de trajets incidents par sous-bande soit très inférieur au nombre d’antennes du réseau d’antennes. Certaines techniques de goniométrie parcimonieuse 1D (i.e. estimation du gisement uniquement, en supposant l’élévation nulle) permettent d'estimer (ou « goniométrer ») l'incidence de potentiellement plus de trajets que d’antennes du réseau. Avec ces techniques, un intercepteur pourrait ainsi être mis en œuvre avec largement moins de sous-bandes et donc avec des bandes instantanées de traitement beaucoup plus larges. Toutefois elles ne traitent pas le problème de la goniométrie 2D (avec estimation conjointe du gisement et de l'élévation) pour une géométrie de réseau quelconque et elles sont adaptées uniquement au cas de trajets décorrélés. L’algorithme MUSIC, lui, permet d’estimer la direction d'arrivée de plusieurs trajets corrélés à condition qu'ils ne soient pas cohérents (trajets dont la différence de temps d'arrivée est quasiment nulle en étant très inférieur à l’inverse de la bande passante du signal) et qu'ils ne dépassent pas le nombre d’antennes du réseau. D’autre part, l’algorithme nécessite un temps d’observation du signal important pour bien estimer les directions des trajets corrélés. Cette durée d’observation nécessaire augmente avec le taux de corrélation des trajets. Il existe donc un besoin d’une solution de détermination des incidences 2D dans les situations de canaux de propagation à multi-trajets corrélés, tels que rencontrés notamment dans les communications de type ionosphérique HF, avec un nombre de trajets incidents supérieur au nombre d’antennes du réseau d’antennes, et qui ne présentent pas les inconvénients de l’art antérieur. A cet effet, suivant un premier aspect, la présente invention décrit un procédé de détermination de la direction d’arrivée des trajets de signaux ionosphériques Hautes Fréquences sur un réseau de N capteurs, avec N ≥ 2 et selon lequel ledit procédé est mis en œuvre par un dispositif électronique de détermination de la direction d’arrivée des trajets et est caractérisé en ce qu’il comprend les étapes suivantes : i. estimation, par goniométrie parcimonieuse, du vecteur  vérifiant l’équation : où est la matrice identité de taille NxN, σ 2 est le niveau du bruit et R x est la matrice de covariance vectorisée du signal x ( t) considéré en sortie du réseau de capteurs, avec, si est le signal en sortie du n ième capteur : où où le vecteur a (θ,Δ) de dimension N x1 est la réponse, connue, du réseau de capteurs à un signal reçu de direction d’arrivée (θ,Δ), et Q est le nombre de gisements potentiels considérés pour et Q' est le nombre d’élévations potentielles par gisement, , pour dans l’espace discrétisé d’incidence du réseau de capteurs selon une grille 2D de gisements et d’élévations , pour avec et , où ii. identification des composantes non nulles, parmi les composantes , …, dudit estimé, et déduction, en fonction de l’indice q m de chacune desdites composantes non nulles, du gisement de trajets de signaux depuis un émetteur respectif m ; iii. pour chaque gisement déduit, étant donné que , les composantes non nulles de type desont ensuite déterminées : l’élévation d’indice q p associée à chacune de ces composantes non nulle , est sélectionnée comme égale à une élévation d’un trajet du m ième émetteur. L’invention permet ainsi de déterminer les incidences dans les situations de canaux de propagation ionosphériques HF à multi-trajets corrélés, et ceci avec un nombre de trajets incidents inférieur ou supérieur au nombre de capteurs (par exemples d’antennes) du réseau de capteurs, la solution nécessitant en outre un temps de calcul et une puissance de calcul réduits. La solution proposée permet d’exploiter les avantages du modèle de la matrice de covariance vectorisée en présence de trajets corrélés qui permettent de : potentiellement traiter plus de trajets que d’antennes du réseau de goniométrie travailler avec des réseaux de goniométrie ayant un nombre d'antennes plus faible ; et ne nécessite pas la détection du nombre de trajets. La solution proposée permet en outre de déterminer les directions d’arrivées de trajets provenant d’émetteurs non « maîtrisés », i.e. dont on ne connaît pas, ou au mieux de façon imparfaite, la norme du signal émis. Dans des modes de réalisation, un tel procédé comprendra en outre l’une au moins des caractéristiques suivantes : - l’estimation de la direction d’arrivée comprend l’estimation d’un angle d’élévation et d’un angle de gisement ; – le procédé comprend au moins une des dispositions suivantes : le nombre de trajets est supérieur au nombre N de capteurs ; des trajets sont corrélés temporellement entre eux. Suivant un autre aspect, l’invention décrit un programme d’ordinateur destiné à être stocké dans la mémoire d’un bloc électronique de traitement et comprenant en outre un microcalculateur, ledit programme d’ordinateur comprenant des instructions qui, lorsqu’elles sont exécutées sur le microcalculateur, mettent en œuvre les étapes i, ii, iii et iv d’un procédé de détermination de la direction d’arrivée des trajets de signaux sur un réseau de N capteurs suivant le premier aspect de l’invention. L’invention décrit également un support non transitoire lisible par ordinateur stockant un tel programme informatique. Suivant un autre aspect, l’invention décrit un dispositif de détermination de la direction d’arrivée des trajets de signaux ionosphériques Hautes Fréquences sur un réseau de N capteurs, avec N ≥ 2, ledit dispositif électronique de détermination de la direction d’arrivée des trajets étant adapté pour estimer par goniométrie parcimonieuse, le vecteur  vérifiant l’équation : où est la matrice identité de taille NxN, σ 2 est le niveau du bruit et R x est la matrice de covariance vectorisée du signal x ( t) considéré en sortie du réseau de capteurs, avec, si est le signal en sortie du n ième capteur : où où le vecteur a (θ,Δ) de dimension N x1 est la réponse, connue, du réseau de capteurs à un signal reçu de direction d’arrivée (θ,Δ), et Q est le nombre de gisements potentiels considérés pour et Q' est le nombre d’élévations potentielles par gisement, , pour dans l’espace discrétisé d’incidence du réseau de capteurs selon une grille 2D de gisements et d’élévations , pour avec et , où ; ledit dispositif étant en outre adapté pour identifier les composantes non nulles, parmi les composantes , …, dudit estimé, et pour déduire, en fonction de l’indice q m de chacune desdites composantes non nulles, le gisement de trajets de signaux depuis un émetteur respectif m ; ledit dispositif étant en outre adapté, pour chaque gisement déduit, étant donné que , pour déterminer les composantes non nulles de type deet sélectionner l’élévation d’indice q p associée à chacune de ces composantes non nulle, comme égale à une élévation d’un trajet du m ième émetteur. Dans des modes de réalisation, un tel dispositif comprendra en outre l’une au moins des caractéristiques suivantes : - l’estimation de la direction d’arrivée comprend l’estimation d’un angle d’élévation et d’un angle de gisement ; – le procédé comprend au moins une des dispositions suivantes : le nombre de trajets est supérieur au nombre N de capteurs ; des trajets sont corrélés temporellement entre eux. L’invention sera mieux comprise et d’autres caractéristiques, détails et avantages apparaîtront mieux à la lecture de la description qui suit, donnée à titre non limitatif, et grâce aux figures annexées, données à titre d’exemple. FIG. 1 La FIG. 1 représente schématiquement un système de goniométrie 1 dans un mode de réalisation de l’invention ; FIG. 2 La FIG. 2 illustre un exemple de canal de propagation multi-trajet en HF ; FIG. 3 La FIG. 3 illustre un des modèles utilisés dans un procédé de goniométrie dans un mode de réalisation de l’invention ; FIG. 4 La FIG. 4 représente des étapes d’un procédé de goniométrie dans un mode de réalisation de l’invention. Des références identiques peuvent être utilisées dans des figures différentes lorsqu’elles désignent des éléments identiques ou comparables. La FIG. 1 représente schématiquement un système de goniométrie 1 dans un mode de réalisation de l’invention, comportant un réseau de goniométrie 2 et un bloc électronique de traitement 50. Le réseau de goniométrie 2 est un réseau de capteurs de signaux électromagnétiques, typiquement des signaux radiofréquences (RF) ou sonores. Dans l’exemple considéré ici, les signaux sont radiofréquences et le réseau de goniométrie 2 est un réseau d’antennes qui comporte N antennes de réception 3 (N≥2 pour un fonctionnement en 1 dimension (1D) et N≥3 pour un fonctionnement en 2 dimensions (2D)) disposées sur un support, ici au sol. Chaque antenne 3 est adaptée pour recevoir des signaux radiofréquences incidents (les trajets) en la position de l’antenne et les convertir en signaux électriques représentatifs des signaux incidents reçus. Les dispositions relatives des antennes 3 dépendent du type de signaux RF. Chaque signal RF incident est un signal dont la source d’émission est un émetteur radiofréquence 4 (par exemple un téléphone mobile). Le bloc de traitement 50 est adapté pour recevoir les signaux électriques en provenance des différentes antennes 3 du réseau d’antennes 2, pour effectuer des traitements du signal sur ces signaux et pour déterminer les directions d’arrivée, encore appelées « incidences »), définie par le gisement et l’élévation du trajet, de chacun des trajets reçus par les antennes 3, en mettant en œuvre les étapes décrites plus bas en référence à la FIG. 4 . Le gisement correspond à l’angle de la direction d’arrivée projetée dans le plan horizontal (le plan horizontal est le plan tangent à la surface de la Terre en la localisation du réseau d’antennes) par rapport à un axe de référence du réseau de capteurs. L’élévation correspond à l’angle formé entre le plan horizontal et la direction d’arrivée. L’élévation peut également être appelée angle de site. Dans un mode de réalisation, le bloc de traitement 50 comporte une mémoire 51 et un microprocesseur 52, et la mémoire 51 compte des instructions logicielles qui lorsqu’exécutées par le microprocesseur 52, mettent en œuvre les étapes incombant au bloc de traitement 50 et décrites plus bas en référence à la FIG. 4 . Les trajets reçus proviennent d’un ensemble de M émetteurs 4. Le nombre M est a priori non connu, comme le nombre de multi-trajets reçus en provenance d’un même émetteur. On considère que l’on est en présence a priori de trajets temporellement corrélés entre eux. Modélisation du signal reçu sur le réseau de goniométrie et déduction de la structure de la matrice de covariance vectorisée : En présence de M émetteurs ayant chacun P m trajets, m =1 à M, le signal en sortie du réseau 2 des N antennes 3 du réseau de goniométrie 1 s’écrit de la façon suivante : (équation 1) . où est le signal en sortie de la n ième antenne, le signal du m ième émetteur 4 et n (t) le vecteur bruit. Les signaux de chacun des émetteurs 4 sont supposés décorrélés (i.e. les signaux d’un quelconque des émetteurs 4 sont supposés décorrélés des signaux d’un quelconque autre des émetteurs 4). Le vecteur , de dimension N x1, caractérise la relation entre le signal mesuré par le réseau d’antennes et le signal incident unitaire de direction d’arrivée . On l’appelle réponse du réseau de goniométrie 2 à un signal reçu d’incidence Θ où Θ = (θ,Δ), θ étant le gisement et Δ l’élévation. Cette réponse est connue ; elle a par exemple été déterminée dans une phase préalable, soit par calcul, soit par calibration expérimentale, pour une valeur de Θ continue quelconque ou alors pour un ensemble prédéterminé de valeurs d’incidence. Par exemple, pour un réseau où la réponse ne dépend que de la position des capteurs, le vecteur établi par calcul s’écrit : avec , où est le vecteur de direction d’onde et la position du n ième capteur. Le p ième trajet du m ième émetteur 4 est caractérisé par son incidence Θ mp = ( θ mp , Δ mp ), son retard de propagation τ mp et son amplitude ρ mp (dans le cas particulier HF, on rappelle que le gisement des différents trajets de l’émetteur est le même pour les différents trajets). Le signal radioélectrique s m (t) d’amplitude complexe, de bande B m et de fréquence centrale f m s’exprime par exemple sous la forme : Le taux de corrélation entre deux trajets (d’indices p, q ) du m ième émetteur 4 peut s’approximer par : (équation 2) . où B m est la bande passante de cet émetteur, E [ ] est l'espérance mathématique et * désigne le conjugué. Ainsi, on considère que deux trajets, d’indices p, q, du m ième émetteur 4 sont corrélés si : (équation 3) . Dans le cas particulier, auquel on s’intéresse plus particulièrement ici, des signaux ionosophériques HF, où la bande passante est souvent de l’ordre de 3kHz, les trajets sont alors corrélés lorsque la différence de temps d’arrivée est inférieure à 0.33ms. On remarque par ailleurs que est la puissance du p ième trajet tel que . Le modèle de l’équation (1) peut se réécrire de la manière suivante (équation 4) . En supposant que le bruit soit décorrélé des signaux des émetteurs 4 et qu’il soit blanc spatialement avec où matrice identité de taille NxN et H définit l’opérateur de transposé et conjugué (σ 2 est donc le niveau du bruit), la matrice de covariance des signaux reçus s’écrit alors : (équation 5) En pratique la matrice de covariance est estimée à partir de K échantillons temporels de la manière suivante : (équation 6) . Dans ce qui suit, la matrice de covariance vectorisée est utilisée : (équation 7) . D’après les équations (4)(5) on peut montrer que : (équation 8) . où désigne le produit de Kronecker. Cette dernière expression peut se réécrire de manière suivante : (équation 9 ) . où l’amplitude du vecteur e dépend du niveau de corrélation entre les trajets. Il est ici rappelé que la vectorisation d’une matrice, par exemple la matrice R, nommée vec(R,) est un vecteur dans lequel les colonnes de ladite matrice sont disposées les unes en dessous des autres. D’après l’équation 9 et pour l’exemple de donné précédemment, on en déduit que dans cette configuration, est un vecteur de dimension dont la ème composante s’exprime avec , où et . Le vecteur correspond donc à la réponse d’un réseau virtuel où les antennes seraient positionnées aux positions . On considère un échantillonnage de l’espace d’incidence du réseau d’antennes selon une grille 2D de gisements et d’élévations contenant incidences potentielles notées pour , avec et . Par exemple, on pourra prendre un pas de grille fixe de 4° selon chacun des angles ; est alors égal à 90x22. On note alors la table de calibration suivante : (équation 10) En première approximation, on suppose que pour chaque incidence de trajet reçu par le réseau d’antennes, il existe un indice q compris entre 1 et tel que ; ce sont ces incidences-là, parmi les incidences possibles, qu’il va s’agir d’identifier. On notera ces indices q mp dans les équations 10 et 11 : Dans la pratique, sera l’incidence la plus proche de dans la grille des incidences possibles. On en déduit que (équation 11) . On obtient un modèle parcimonieux de la matrice de covariance vectorisée perturbé par un bruit dépendant des niveaux de corrélation entre les trajets de chacun des émetteurs 4 ainsi que du niveau σ 2 du bruit. La matrice de covariance vectorisée pourrait s’écrire de la manière suivante : avec (équation 12) . En prenant en compte la caractéristique particulière des incidences des multi-trajets de signaux ionosphériques HF où , on définit la table de calibration suivante : (équation 13) . avec Q le nombre de gisements potentiels considérés et Q' le nombre d’élévations potentielles par gisement, dans l’espace discrétisé d’incidence du réseau d’antennes selon une grille 2D de gisements et d’élévations pour , pour avec et et on en déduit le deuxième modèle parcimonieux suivant : avec (équation 14) où chaque vecteur de taille représente, relativement au gisement de trajets qui lui est associé, la vectorisation de la matrice de covariance des potentiels trajets et où contient des composantes non nulles seulement quand l’indice q correspond à un gisement de trajets correspondant à un émetteur d’indice m dont proviennent des signaux reçus par le réseau de capteurs ; dans ce cas q est noté q m : où les matrices R m sont de dimension Q’xQ’ et où la composante R m (i,j) de R m , est non nulle seulement lorsqu’il existe soit un trajet si , soit deux trajets (le trajet et le trajet de l’émetteur ) sinon, depuis ledit émetteur d’indice m, d’élévations associées aux indices i, j ; dans ce cas i (indice d’élévation relative au trajet ) et j (indice d’élévation relative au trajet ) sont notés respectivement et , et (dans le cas d’un seul trajet k, l est pris égal à k dans la précédente formule. . En considérant le modèle de l’équation (12), le vecteur est composé au maximum de N 2 - N +1 composantes indépendantes avec un réseau de N antennes ; la matrice est composée de ( P m ) 2 colonnes car est de dimension N x P m . Le nombre de composantes non nulles contenu dans les vecteurs { γ 1 ,..., γ M } vaut alors . On en déduit alors que la condition d’identifiabilité des trajets est la suivante : (équation 15) En supposant que le nombre de trajets par émetteur vérifie , la capacité en nombre d’émetteurs devient alors : . Principe de la goniométrie parcimonieuse Le domaine de l’estimation parcimonieuse à partir d’un modèle d’observation regroupe les principes mathématiques permettant de résoudre un système d'équations ayant plus d'inconnues que d'équations sachant que la plupart des inconnues sont nulles. L’objectif de la goniométrie parcimonieuse, dans une configuration de trajets non corrélés, est de résoudre le système d’équations pour estimer le vecteur parcimonieux γ . Selon les équations (10)(11), les incidences des trajets reçus correspondant au vecteur x (t) (tel que considéré, donc, aux K instants d’échantillonnages à utilisés pour approximer la matrice de covariance) correspondent alors aux composantes non nulles du vecteur γ où , et donc on peut en déduire que l’incidence correspondant à la q mp ème composante de est une incidence d’un trajet reçu à l’instant t, soit . Cela est illustré sur la FIG. 3 qui illustre comment l’estimation du spectre parcimonieux γ permet de déterminer quelles sont les incidences des trajets reçus parmi l’ensemble des incidences potentielles. Dit autrement : si la q ième composante de γ , i.e., composante ᵧ(q), est déterminée, à l’issue du traitement de goniométrie parcimonieuse, comme non nulle, alors l’incidence , correspondant à la q ième colonne de B , b( ) , est déterminée comme incidence d’un trajet reçu. On est toutefois en présence d’un système d’équations ayant plus d’inconnues que d’équations. Pour résoudre ce système d’équations, il est alors nécessaire d’exploiter la parcimonie du vecteur γ . Pour cela, on ajoute au critère classique des moindres carrés une régularisation ou pénalité Φ(λ, γ ) qui a pour objectif de favoriser la parcimonie de la solution devant contenir beaucoup de composantes nulles. Le critère à minimiser suivant le vecteur x a la forme suivante : (équation 17) Le poids relatifs des deux termes est réglé par un paramètre de régularisation λ afin de garantir la relation pour tous vecteurs x . Selon les résultats exposés dans « Goniométrie parcimonieuse de sources radio-électriques : modèles, algorithmes et mises en œuvre robustes », Alice Delmer, thèse de doctorat, Université Paris-Saclay, Décembre 2021, ce poids peut être fixé a priori à partir de la limite de résolution du réseau de capteurs virtuels de réponse b (Θ). La pénalité Φ(λ, x ) a pour expression : (équation 18) Les pénalités CELO et la pénalité l 0 (la pénalité l 0 étant définie comme la norme l 0 multipliée par ) utilisées dans la thèse indiquée ci-dessus ont pour point commun de vérifier . Dans le cas de l 0 , la valeur de ε tend vers 0 alors que, dans le cas de CELO, la fonction est plus large au voisinage de 0. D’autres pénalités possédant les mêmes propriétés ont été proposées dans la littérature. . L’algorithme d’estimation du vecteur γ est itératif avec une initialisation avec le vecteur x  . La relation de récurrence entre deux étapes de l’algorithme est la suivante : (équation 19) . où est le proximal de la pénalité. Dans le cas de l 0 et CELO, l’opérateur proximal vérifie : cas de l 0 : cas de CELO : avec (.)+=max(0,.). Les étapes d’estimation du vecteur parcimonieux γ sont résumées ci-dessous : Etape A.0 : Initialisation par le vecteur x  et i =0 Etape A.1 : Calcul de (correspondant à une descente de gradient sur les moindres carrés) Etape A.2 : Mise à jour du vecteur avec le proximal : Etape A.3 : Si retour à l’étape A.1 avec i = i +1. Sinon, le vecteur parcimonieux est estimé égal à . La valeur de ɛ correspond à la précision d’estimation voulue. Elle peut être reliée à un niveau de bruit moyen s’il est connu. D’autres critères d’arrêts sont aussi possibles, comme par exemple un nombre maximum d’itérations ou la décroissance du critère entre deux itérations. La pénalité CELO permet d’initialiser avec x  =0 contrairement à l 0 dans le cas général à cause de la présence de minima locaux. Les étapes d’estimation du vecteur parcimonieux γ peuvent être réalisées avec d’autres pénalités de la littérature comme MCP. Il est également envisageable d’utiliser des schémas plus complexes basés sur ces étapes simples : pénalité ou proximal variant dans le temps, paramètres dépendant des itérations, moyennage du vecteur sur plusieurs itérations… Processus de goniométrie parcimonieuse en présence de trajets corrélés En considérant maintenant que des trajets provenant d’un même émetteur 4 sont corrélés entre eux, l’approche proposée par l’invention est d’utiliser le modèle provenant des équations (12)-(14) exploitant la modélisation en gisement et élévation des multi trajets des signaux ionosphériques HF. Dans un mode de réalisation de l’invention illustré en FIG. 4 , un procédé 100 utilisant cette approche est mis en œuvre par le bloc de traitement 50 ; ce procédé 100 comprend ainsi les étapes C.0 à C.3 détaillées ci-après. Dans une étape C .0 , le bloc de traitement calcule en fonction des signaux x (t) obtenus en sortie des N antennes et en déduit y à l’aide de la formule suivante : . Puis, dans une étape C.1 , estimation du vecteur par un processus de traitement parcimonieux tel qu’exposé aux étapes A.0 à A.3 (en y remplaçant la matrice B par la matrice C définie par l’équation (13) et en y remplaçant γ par ). Etape C.2 : de , identification des vecteurs avec des composantes non nulles (ou autrement dit, de norme non nulle), parmi les sous-vecteurs , …, de ainsi estimé (rappel : est un vecteur de non nul seulement quand l’indice q correspond à un gisement des trajets d’un émetteur d’indice m dont proviennent des signaux reçus par le réseau de capteurs; dans ce cas q est noté q m ) et déduction des gisements des émetteurs incidents, ces gisements étant ceux associés aux seuls vecteurs non nuls, i.e. pour q = q m (le nombre de vecteurs ayant des composantes non nulles fournit en outre la valeur de M, le nombre d’émetteurs à l’origine des signaux reçus). Etape C.3 : pour chaque gisement déduit de l’étape C.2 (et donc pour chaque émetteur m ) , m = 1 à M, identifié à l’étape C.2, étant donné que : les composantes non nulles de type desont ensuite identifiées : chacune de ces composantes non nulles est associée à une élévation d’indice q p et celle-ci est sélectionnée comme correspondant à l’élévation d’un trajet du m ième émetteur associé au gisement déduit. L’estimation obtenue par les étapes C.0 à C.3 prend en compte à la fois la corrélation entre les trajets et la particularité des trajets HF. L’estimation est beaucoup plus précise que l’estimation obtenue par les étapes A.0 à A.3 avec le dictionnaire où la corrélation des trajets augmente sensiblement le niveau du bruit. La prise en compte de la corrélation entre les trajets sans prise en compte de la spécificité des trajets HF impliquerait d’étendre le dictionnaire par un dictionnaire qui contient tous les termes croisés et dont la taille serait prohibitive. Ces directions d’arrivée {} estimées sont ensuite exploitées, par exemple pour localiser les émetteurs 4, ce qui est utilisé par le bloc de traitement 50 ou un autre dispositif électronique, dans le contrôle du spectre, la gestion des télécommunications, le suivi des porteurs d’émetteur (recherche de victimes enfouies après des catastrophes naturelles, protection d’espèces animales …). Le procédé peut être mis en œuvre par l’exécution d’instructions logicielles sur un processeur comme décrit plus haut à titre d’exemple. Alternativement, il peut être mis en œuvre par un hardware dédié, typiquement un circuit intégré numérique, soit spécifique (ASIC) ou basé sur une logique programmable (par exemple FPGA/Field Programmable Gate Array). Procédé de détermination de la direction d’arrivée des trajets de signaux ionosphériques Hautes Fréquences sur un réseau (2) de N capteurs (3), avec N ≥ 2 et selon lequel ledit procédé est mis en œuvre par un dispositif électronique (50) de détermination de la direction d’arrivée des trajets et est caractérisé en ce qu’il comprend les étapes suivantes : i. estimation, par goniométrie parcimonieuse, du vecteur  vérifiant l’équation : où est la matrice identité de taille NxN, σ 2 est le niveau du bruit et R x est la matrice de covariance vectorisée du signal x ( t) considéré en sortie du réseau de capteurs, avec, si est le signal en sortie du n ième capteur : où où le vecteur a (θ,Δ) de dimension N x1 est la réponse, connue, du réseau de capteurs à un signal reçu de direction d’arrivée (θ,Δ), et Q est le nombre de gisements potentiels considérés pour et Q' est le nombre d’élévations potentielles par gisement, , pour dans l’espace discrétisé d’incidence du réseau de capteurs selon une grille 2D de gisements et d’élévations , pour avec et , où ii. identification des composantes non nulles, parmi les composantes , …, dudit estimé, et déduction, en fonction de l’indice q m de chacune desdites composantes non nulles, du gisement de trajets de signaux depuis un émetteur respectif m ; iii. pour chaque gisement déduit, étant donné que , les composantes non nulles de type desont ensuite déterminées : l’élévation d’indice q p associée à chacune de ces composantes non nulle , est sélectionnée comme égale à une élévation d’un trajet du m ième émetteur. Procédé de détermination de la direction d’arrivée des trajets de signaux sur un réseau (2) de N capteurs (3) selon la revendication 1, selon lequel l’estimation de la direction d’arrivée comprend l’estimation d’un angle d’élévation et d’un angle de gisement. Procédé de détermination de la direction d’arrivée des trajets de signaux sur un réseau (2) de N capteurs (3) selon l’une quelconque des revendications précédentes, comprenant au moins une des dispositions suivantes : le nombre de trajets est supérieur au nombre N de capteurs ; des trajets sont corrélés temporellement entre eux. Programme d’ordinateur, destiné à être stocké dans la mémoire d’un bloc électronique de traitement (50) et comprenant en outre un microcalculateur, ledit programme d’ordinateur comprenant des instructions qui, lorsqu’elles sont exécutées sur le microcalculateur, mettent en œuvre les étapes i, ii, iii et iv d’un procédé de détermination de la direction d’arrivée des trajets de signaux sur un réseau de N capteurs selon l’une des revendications précédentes. Dispositif (50) de détermination de la direction d’arrivée des trajets de signaux ionosphériques Hautes Fréquences sur un réseau de N capteurs, avec N ≥ 2, ledit dispositif électronique de détermination de la direction d’arrivée des trajets étant adapté pour estimer par goniométrie parcimonieuse, le vecteur  vérifiant l’équation : où est la matrice identité de taille NxN, σ 2 est le niveau du bruit et R x est la matrice de covariance vectorisée du signal x ( t) considéré en sortie du réseau de capteurs, avec, si est le signal en sortie du n ième capteur : où où le vecteur a (θ,Δ) de dimension N x1 est la réponse, connue, du réseau de capteurs à un signal reçu de direction d’arrivée (θ,Δ), et Q est le nombre de gisements potentiels considérés pour et Q' est le nombre d’élévations potentielles par gisement, , pour dans l’espace discrétisé d’incidence du réseau de capteurs selon une grille 2D de gisements et d’élévations , pour avec et , où ; ledit dispositif (50) étant en outre adapté pour identifier les composantes non nulles, parmi les composantes , …, dudit estimé, et pour déduire, en fonction de l’indice q m de chacune desdites composantes non nulles, le gisement de trajets de signaux depuis un émetteur respectif m ; ledit dispositif étant en outre adapté, pour chaque gisement déduit, étant donné que , pour déterminer les composantes non nulles de type deet sélectionner l’élévation d’indice q p associée à chacune de ces composantes non nulle, comme égale à une élévation d’un trajet du m ième émetteur. Dispositif (50) de détermination de la direction d’arrivée des trajets selon la revendication 5, dans lequel l’estimation de la direction d’arrivée comprend l’estimation d’un angle d’élévation et d’un angle de gisement. Dispositif (50) de détermination de la direction d’arrivée des trajets selon la revendication 5 ou 6, où l’une au moins des deux dispositions suivantes est vérifiée : le nombre de trajets est supérieur au nombre N de capteurs ; des trajets sont corrélés temporellement entre eux.