Brevet FR2747787A1 — publié le 19971024

Procédé et dispositif de localisation de sources radloélec- triques L'invention concerne un procédé et un dispositif passif de localisation angulaire simultanée de sources d'ondes radioélectriques. Elle s'applique avantageusement à la localisation en site et/ou gisement, dans un repère spatial convenablement choisi, d'émetteurs par rapport à un point d'observation, par exemple d'aéronefs coopératifs autour d'un aérodrome. On sait que les procédés de localisation de sources d'ondes radioélectriques mettent en oeuvre des dispositifs d'antennes associés à un système de réception et de traitement des signaux captés. Dans le procédé le plus classique, dit à "formation de voies", on considère que chaque onde, provenant d'une source de bande passante étroite et captée par un réseau d'antennes, est plane. Les déphasages des signaux captés, résultant des différences de trajet de propagation de l'onde vers chaque antenne, sont corrigés puis les signaux sont sommés, la localisation consistant à rechercher une valeur maximale du signal somme parmi celles calculées pour toutes les directions d'observation. Ce procédé présente l'inconvénient principal de nécessiter un nombre important d'antennes pour pouvoir distinguer deux sources angulairement proches. D'autres procédés sont basés sur une approche statistique et sur des hypothèses faites sur la nature des signaux à traiter. Dans un premier procédé dit "à haute résolution", les sources sont assimilées à des processus gaussiens, et les signaux captés par les antennes peuvent alors être caractérisés par leurs moyennes et leurs covariances aussi appelées, respectivement, moment d'ordre 1 et moment d'ordre 2. Théoriquement, N antennes de réception permettent alors de déterminer jusqu'à N-l directions de sources, cependant le nombre de sources réellement localisées est très inférieur à celui estimé après avoir pris en compte les dégradations de performances prévisibles entre des résultats théoriques et expérimentaux. Dans un deuxième procédé on prend en compte la nature des signaux réels en les caractérisant non seulement par leurs moyennes et leurs covariances, mais aussi par des moments d'ordres supérieurs à 2. Dans le cas particulier de l'ordre 4, N antennes de réception peuvent permettre de détecter jusqu'à N*(N-1) directions de sources. La recherche de la valeur approchée de ces moments d'ordres supérieurs améliore donc le pouvoir de résolution angulaire, mais ce procédé nécessite des calculs nombreux et complexes limitant la localisation en temps réel et impliquant la mise en oeuvre de calculateurs coûteux. Un autre procédé, plus récent, basé aussi sur une approche statistique et sur des hypothèses faites sur la nature des signaux à traiter, consiste à construire un histogramme bidimensionnel se présentant sous forme d'une carte de points établie à partir des échantillons numérisés des signaux captés par les antennes. On peut alors estimer les directions des sources en recherchant sur l'histogramme les maxima de densité de points. Ce procédé ne s'applique qu'au traitement de signaux non gaussiens émis par un nombre restreint de sources fortement modulées. La présente invention a pour but de perfectionner ce dernier procédé en permettant de l'appliquer à une large plage de signaux non gaussiens, émis par des sources d'ondes radioélectriques, et mettant en oeuvre un petit nombre d'antennes tout en permettant la localisation angulaire d'un grand nombre de sources, avec des moyens de calcul limités. Le dispositif proposé à cet effet comprend: - au moins deux antennes, définissant au moins deux axes d'observation, - des moyens de réception, reliés aux antennes, et fournissant, à un instant donné, et pour chaque axe d'observation, une valeur numérique de l'amplitude du signal reçu, reliée à la position angulaire des sources audit instant, et - des moyens de mémorisation de ces valeurs numériques à des instants successifs. - des moyens de traitement propres à déterminer des variations de densité de l'histogramme multidimensionnel constitué, dans un repère, par des points dont les coordonnées sont lesdites valeurs numériques mémorisées simultanément. Selon un aspect de l'invention, les moyens de traitement recherchent les orientations des contours de l'histogramme dans ledit repère. Dans une version de base, on s'intéresse à un couple de signaux captés selon deux axes. L'histogramme est alors bidimensionnel. Ses contours sont des segments de droite. La recherche effectuée par les moyens de traitement comprend le dénombrement de ces segments de droite, puis l'ajustement linéaire, sur le contour de l'histogramme, d'une forme polygonale ayant le nombre voulu de segments. Le couple de signaux captés peut provenir directement de deux antennes ou également d'antennes multiples selon la technique dite de "formation de voies". Avec trois antennes, et trois signaux captés, il est possible d'accéder à un histogramme tridimensionnel. Ses contours sont alors des faces polygonales. La recherche effectuée par les moyens de traitement comprend le dénombrement de ces faces polygonales, l'ajustement, sur le contour de l'histogramme, de plans délimitant un polyèdre dont l'orientation des arêtes permet de déterminer la direction des sources. D'autres caractéristiques et avantages de l'invention apparaîtront à l'examen de la description détaillée ci-après, ainsi que des dessins annexés, sur lesquels - la figure 1 est un diagramme tridimensionnel utile à la compréhension des bases de la présente invention - la figure 2 est l'illustration schématique des diagrammes de rayonnement de deux antennes situées au même point et d'orientations perpendiculaires, sur la base desquels seront décrits les principes de la présente invention - les figures 3A et 3B sont des diagrammes montrant des exemples de distributions de signaux captés en provenance de deux sources lors de la mise en oeuvre de l'invention avec des antennes telles que définies à la figure 2 - la figure 3C est un diagramme montrant un exemple de distribution de signaux captés en provenance de trois sources lors de la mise en oeuvre de l'invention avec des antennes quelconques - les figures 4A et 4B décrivent respectivement un mode de réalisation d'un dispositif selon la présente invention et une variante de celui-ci, au niveau des antennes - la figure 5 est un mode de réalisation détaillé du dispositif selon l'invention ; - la figure 6A illustre une autre disposition, à trois antennes, utilisable pour la mise en oeuvre de l'invention ; - la figure 6B illustre les bases de définition d'un histogramme tridimensionnel avec le système d'antennes de la figure 6A ; et - la figure 7 illustre une fonction d'allongement à partir de laquelle on détermine le nombre de sources localisées et leurs directions. Les dessins annexés sont, pour l'essentiel, de caractère certain. En conséquence, ils pourront non seulement servir à mieux faire comprendre la description, mais aussi contribuer à la définition de l'invention, le cas échéant. La localisation de sources selon la présente invention utilise des systèmes récepteurs passifs, associés à plusieurs antennes, (ou réseau d'antennes) ayant des diagrammes de réception connus. Toute source (E) de rayonnement électromagnétique est située dans la bande de travail des antennes, et à distance suffisante de celles-ci pour que l'on puisse considérer que les signaux captés sont du type "champ lointain" (ondes planes). Chaque antenne Aj, définie par son diagramme de réception du(3) et par son axe d'observation xj, correspondant au maximum de gain de réception du diagramme, capte un signal Sj(t), en provenance des sources émettrices Ei à localiser selon la direction Bi, et rayonnant des signaux ei(t) d'amplitude maximale Mi, tel que Sj(t) = Si ei(t).dj(ei). (I) S'agissant d'une réception purement passive, il n'existe a priori aucune relation de phase entre la base de temps de la source émettrice et celle du système récepteur. En conséquence, pour le récepteur passif, les phases de chaque source de rayonnement électromagnétique se présentent comme des variables aléatoires, équiréparties entre 0 et 2 PI radians. Le Demandeur a en outre remarqué que de nombreux signaux sont émis sans modulation d'amplitude, ou avec une modulation d'amplitude faible : c'est le cas pour les signaux radar, certains signaux de télécommunication, ou des signaux de brouillage, notamment. Ces signaux représentent une part importante des signaux-source à traiter. La figure 1 illustre l'évolution de la forme de l'histogramme des échantillons d'un modèle de signal émis par une source en fonction de l'indice de modulation qui lui est appliqué. Sur l'axe des X figurent les valeurs de l'amplitude normée des échantillons du signal de la source d'émission. Sur l'axe des Y figure l'indice de modulation d'amplitude de 0 à 1 de la source d'émission. Sur l'axe des Z figure la densité de probabilité des échantillons de la source d'émission. Ainsi, le Demandeur observe que - pour une source non modulée, l'histogramme des niveaux émis est quasi-uniforme sur l'intervalle où il existe, hormis les valeurs limites asymptotiques qui ne jouent en pratique qu'un rôle faible, comme l'illustre le plan avant de la figure 1 (y=o); - lorsqu'on examine les différents plans de la figure 1 pour un indice de modulation qui va de 0 à 0,8, il apparaît que l'histogramme des niveaux émis ne se déforme que très peu; - c'est seulement lorsque l'indice de modulation approche de la valeur 1 que l'histogramme prend une allure différente. Le Demandeur en déduit qu'une procédure de dénombrement et de localisation adaptée à ces différents types de signaux améliore les performances de localisation, ctest-à-dire le rapport du nombre de sources localisées par rapport au nombre d'antennes utilisées. A cet effet, l'invention prévoit l'usage de l'histogramme bidimensionnel des signaux captés par un couple d'antennes. Cet histogramme se présente comme une carte (nuage) de points dont certaines caractéristiques peuvent être reliées aux directions des sources. Dans le cas de sources fortement modulées, l'histogramme se présente sous forme de lignes de points, qui peuvent être reliées aux directions des sources. Dans le cas des signaux de sources faiblement modulées, ce sont les contours de l'histogramme qui permettent de détecter la direction des sources. Les expérimentations et simulations numériques auxquelles il a été procédé montrent que ce type de traitement est robuste vis-à-vis de l'imprécision du modèle de source choisi, surtout dans le cas de sources faiblement modulées. Il est maintenant fait référence à la figure 2, pour aller plus avant dans la description du principe de l'invention. On dispose donc de deux antennes colocalisées (ayant le même centre de phase) Al et A2, possédant des diagrammes de rayonnement identiques et orthogonaux, d1 et d2, définis par les formules dl (8) = cos (8) ; d2 ( ) = sin 8 (Il) dans lesquelles 8 est l'angle existant entre la direction considérée et la direction d'observation x1 de gain maximum de l'antenne A1. En présence d'un petit nombre de sources électromagnétiques, Po, indépendantes et modulées en amplitude avec un indice de 1, l'histogramme des signaux reçus est un nuage en forme d'étoile à 2Po branches alignées deux à deux, chaque paire ainsi définie indiquant la direction d'une source. En présence de P sources électromagnétiques indépendantes, non modulées en amplitude, l'histogramme des signaux reçus est un nuage limité par P couples de droites parallèles, dont les pentes sont les directions des sources. Ces couples de droites définissent un polygone convexe. En présence de deux sources électromagnétiques, une modulée en amplitude avec un indice voisin de 1 et une non modulée, l'histogramme des signaux reçus est un nuage en forme de polygone concave. Les figures 3A et 3B illustrent les histogrammes bidimensionnels construits à partir des signaux Sl(t) en abscisses et S2(t) en ordonnées, avec P = 2 (deux sources). La figure 3A représente l'histogramme obtenu à partir de deux sources modulées en amplitude avec chacune un indice de 1. La figure 3B représente l'histogramme obtenu à partir de deux sources non modulées en amplitude. On se réfère maintenant à la figure 3B. Les signaux dus à la source El, située dans la direction e1, se placent, dans le plan (S1, S2) à l'intérieur d'un segment de droite Mll-M12, de longueur 2*M1 (deux fois M1) et d'orientation 1 ; M1 et 1 sont respectivement l'amplitude maximale et la direction de la première source. De même, les signaux dus à la source E2 se placent à l'intérieur d'un segment de droite M21-M22, de longueur 2*M2 et d'orientation 02, où M2 et 2 sont respectivement l'amplitude maximale et la direction de la seconde source. La somme vectorielle d'un signal situé sur le segment Mll-M12 avec un signal situé sur le segment M21-M22 est un point situé à l'intérieur du parallélogramme PQRS. Le procédé détaillé proposé pour localiser ces sources comprend les opérations suivantes 1) constitution de l'histogramme bidimensionnel dû aux sources en présence desquelles on se trouve 2) localisation et dénombrement des frontières du nuage de points (dans le cas de signaux faiblement modulés), ou bien de droites passant par l'origine (dans le cas de signaux fortement modulés), et 3) transformation des informations de pente desdites frontières ou desdites droites en informations de gisement des sources. Dans l'exemple qui précède, il a été fait référence à deux antennes situées au même endroit, possédant des diagrammes de rayonnement identiques particuliers, et des axes perpendiculaires l'un à l'autre. L'invention s' applique plus généralement, dès lors que l'on dispose de deux antennes ayant leurs centres de phase situés sensiblement au même point, et possédant des diagrammes de rayonnement connus, et des axes d'observation différents. La mise en oeuvre de l'invention peut être réalisée de deux manières, comme visible sur les figures 4A et 4B. Sur la figure 4A, on dispose d'un couple d'antennes Al et A2 dont les centres de phases sont éloignés d'une distance très inférieure à la longueur d'onde des signaux captés, et possédant des diagrammes de rayonnement dl et d2 d'axes Alxl et A2x2. Ce sont ces antennes Al et A2 qui fournissent directement les signaux Sl(t) et S2(t) en fonction du temps. Sur la figure 4B, on procède par formation de voies. On dispose donc d'un réseau d'antennes Al, ..., Ai, ..., An, avec des déphaseurs W1, ..., Wi, ..., Wn. Les sorties des déphaseurs sont reliées à un sommateur. En doublant le système de déphasage, il est possible de cette manière d'obtenir deux signaux S1 et S2 qui possèdent sensiblement les mêmes propriétés que ceux de la figure 4A. On se réfère maintenant à la figure 5, le couple de signaux S1 et S2 est appliqué à un convertisseur analogique/numérique 20, dont la sortie est reliée à une mémoire 30, qui stocke des couples SDl(k), SD2(k), où l'adjonction de D signifie qu'il s'agit de signaux numériques, et k est un indice chronologique variant de 1 à K, K étant le nombre de couples de valeurs mémorisées. La mémorisation des signaux SDl(k) et SD2(k) sur un temps choisi définit un histogramme bidimensionnel. L'invention met alors en oeuvre deux traitements qui peuvent être réalisés concomitamment. Le premier traitement concerne les signaux faiblement modulés, et commence par la mise en oeuvre de moyens de traitement d'image. L'étape 41, consiste à déterminer dans un repère, orthonormé par exemple, les équations des droites formant les frontières estimées du nuage de points que l'on obtient à partir du contenu de la mémoire 30. Ce nuage se présente sous la forme illustrée sur la figure 3B, dans le cas de deux sources. L'étape 42 consiste à calculer la longueur des segments-frontières portés par les droites estimées. L'étape 43 consiste à sélectionner certains segments selon leurs longueurs et mémoriser leurs pentes respectives, que l'on note bi- Et l'étape 44 consiste à convertir ces pentes en informations angulaires ei. Le second traitement concerne les signaux fortement modulés, et commence par l'étape 51 qui consiste à rechercher dans le même nuage de points des droites passant par l'origine. Ce nuage se présente sous la forme illustrée sur la figure 3A, dans le cas de deux sources. L'étape 53 consiste à estimer les pentes 9i de ces droites. Et l'étape 54 consiste à convertir ces informations de pentes en en informations angulaires Oi. L'unité d'exploitation 60 gère les informations ainsi obtenues, en recherchant notamment s'il s'agit de sources faiblement ou fortement modulées. A cet effet, on peut utiliser des informations pouvant provenir soit du dispositif selon l'invention, grâce à l'interprétation de la forme des histogrammes obtenus, soit d'une autre origine (61). D'une façon générale, les signaux captés étant à valeur moyenne nulle, les formes des histogrammes présentent une symétrie par rapport au centre du repère dans lequel ils sont représentés (figures 3A, 3B et 3C). On se réfère maintenant à la figure 3C. Pour simplifier la compréhension, le repère choisi est orthonormé, d'axes ox-oy. L'équation générale d'une droite y est de la forme y = tgb x + yg (III) où ç représente l'angle formé entre la droite considérée et l'axe ox, c'est-à-dire la pente de la droite. Dans le cas de plusieurs signaux faiblement modulés, l'histo- gramme obtenu doit être ajusté par un polygone convexe ayant des côtés parallèles deux à deux (parallélogramme dans le cas de deux sources). De la symétrie de l'histogramme, on déduit que chaque paire de côtés parallèles sera portée par un couple de droites d'équation : y = tgç x + yg avec yg é O (Iv) Ces deux droites parallèles définissent entre elles une bande de surface dont la largeur varie avec y0 L'étape 41 consiste à calculer pour chaque valeur de , comprises entre 0 et u, la valeur correspondante de yo pour que l'intersection de ladite bande et de l'histogramme comprenne un nombre fixé de points, kg, parmi les K points de l'histogramme. La valeur de kg est voisine de K, à quelques pour cent près. Les deux surfaces de l'histogramme, externes à la bande de surface délimitée par la paire de droites parallèles, ont une forme plus ou moins allongée en fonction de la valeur de b pour laquelle elles ont été définies. L'étape 42 consiste à calculer l'allongement de ces surfaces en fonction de b. Pour ce faire, on détermine d'abord, pour chaque valeur de , les coordonnées xB et yg, dans le repère ox - oy, des deux barycentres des points appartenant aux deux surfaces externes. Les équations permettant de les définir sont pour les a points parmi les K satisfaisant à la relation SD2 (k) > tgs SD1(K) + Yo (V) ainsi qu'aux équations (VI) reportées dans l'annexe de formules en fin de description, et pour les p points parmi les K satisfaisant à la relation SD2(k) On détermine ensuite la moyenne quadratique des projections, sur la droite de pente , des vecteurs délimités par le barycentre et chacun des points appartenant à la même surface externe. C'est cette moyenne qui est utilisée dans l'exemple donné comme représentative de l'allongement et que l'on note par l'équation (X) donnée en fin de description. L'ensemble des valeurs d'allongement obtenues en fonction de s dans le cas de l'histogramme de la figure 3C est représenté figure 7. Cette courbe présente des maxima locaux pour certaines valeurs de o. L'étape 43 consiste à sélectionner, parmi ces maxima locaux, ceux ayant une valeur supérieure à un seuil donné. On obtient donc P valeurs (3 dans notre exemple) de , correspondant chacune à une source localisée. L'étape 44 consiste à convertir les P valeurs de en P valeurs de e, angle de localisation des sources dans le repère des antennes. Un tableau de correspondance des valeurs de b et e est établi lors d'une étape d'apprentissage au cours de laquelle on procède à un étalonnage du dispositif en utilisant une source d'onde radioélectrique connue que l'on déplace de façon à parcourir les diagrammes de réception des antennes avec le pas désiré pour la résolution ultérieure de e. Pour chaque i, après numérisation des signaux captés, on constitue l'histogramme correspondant dans un repère du même type que celui précédemment décrit. L'histogramme se présente comme un ensemble de points alignés passant par l'origine de ce repère. La pente de la droite supportant l'ensemble des points et d'équation y = tg b x (XI) Ainsi pour chaque valeur 0, on détermine lors de l'apprentissage, un angle ç correspondant. Réciproquement en phase de localisations de sources, à chaque angle ç calculé correspond un angle 0, dans le repère des antennes, désignant la direction d'une source radioélectrique localisée. Bien entendu, lors de la phase d'exploitation 60, il est possible d'utiliser deux couples différents d'antennes, pour obtenir deux jeux d'informations angulaires différentes. Il devient ainsi possible par triangulation de connaître non seulement la direction, mais aussi la position effective des sources à localiser. Le dispositif proposé peut s'étendre immédiatement à la localisation bidimensionnelle de sources, par exemple à l'aide d'antennes embarquées sur un aéronef, ce qui donne alors accès à deux mesures angulaires, l'une de site et l'autre de gisement. Ainsi, la figure 6A montre trois antennes Al, A2, A3 de diagrammes de rayonnement identiques ou connus, et écartées entre elles d'angles égaux (leurs centres de phase sont colocalisées). L'histogramme obtenu sera alors tridimensionnel, et délimité par des plans. Les intersections de ces plans deux à deux sont des droites dont l'orientation permet de déterminer les directions des sources. En effet (figure 6B), si l'on se place dans le plan contenant les directions axiales des deux premières sources, on observe que la somme de leurs contributions s inscrit dans un parallélogramme, comme indiqué plus haut. En rajoutant la troisième source, on obtient un volume généré par le parallélogramme en se déplaçant parallèlement à lui-même en fonction de la direction de cette troisième source. Le volume ainsi créé aura une dimension 2*M3 selon cet axe, M3 étant l'amplitude maximale de la troisième source. Comme précédemment, au lieu de trois antennes, on peut également utiliser un réseau d'antennes de distribution bidimensionnelles, et synthétiser par une formation de voies classiques, les trois voies directives, décalées angulairement, nécessaires à la mise en oeuvre de cette variante de l'invention. Les premières expérimentations auxquelles s'est livré le Demandeur ont montré que le procédé et le dispositif ici proposé peuvent donner des résultats très intéressants, suivant le nombre et la nature des sources rencontrées. Lorsqu'on est en présence d'une seule source, il est en principe toujours possible d'en déterminer la direction. Dans le cas théorique où l'on serait en présence de deux sources strictement gaussiennes, la détermination des directions de celles-ci n'est plus possible. Par ailleurs, le cas le plus fréquent, rencontré avec des signaux de télécommunication ou avec des signaux radar, est celui d'un signal source d'amplitude constante et de phase aléatoire. Dans ce cas, l'histogramme s'étale presque uniformément sur le parallélogramme de la figure 3B. Il devient alors possible de localiser un grand nombre de sources. Si l'on tient maintenant compte du rapport signal à bruit, il apparaît que les résultats qui précèdent sont, pour l'essentiel, maintenus. Il est important de noter que, dans le cas précité des sources à phase aléatoire, un rapport signal à bruit de 10 dB peut être considéré comme nécessaire mais suffisant. Il faut aussi prendre en compte le fait que le nombre d'échantillons traités conjointement est à son tour lié à l'intervalle de temps sur lequel on peut considérer les caractéristiques des sources comme stables. S'il n'est pas possible d'augmenter ce temps de mesure, il convient alors d'augmenter la fréquence d'échantillonnage. Là aussi, les premières expériences conduites ont montré qu'avec les sources à phase aléatoire, un nombre d'échantillons de l'ordre de 500 peut suffire pour permettre une localisation des deux sources. Les valeurs de ces différents paramètres peuvent être définies à l'aide de l'unité d'exploitation 60. Enfin, il s'est avéré possible, avec seulement deux antennes, de discriminer angulairement un nombre de sources supérieur à deux, selon le rapport signal sur bruit et le nombre d'échantillons en présence desquels on se trouve. Ceci est à rapprocher du fait suivant : les méthodes actuellement les plus performantes (moment statistique d'ordre 2) résolvent au mieux P(P-1) sources, soit deux sources pour deux capteurs. Ceci illustre tout particulièrement l'intérêt de la présente invention. Sur un autre plan, bien que la présente invention ait actuellement donné tous ses effets avec des antennes directionnelles, certains éléments laissent penser qu'elles pourraient également fonctionner avec des aériens omnidirectionnels. ANNEXE FORMULES REVENDICATIONS 1. Dispositif pour la localisation de sources radioélectriques, comprenant: - au moins deux antennes (A1,A2), définissant au moins deux axes d'observation (Al-xl, A2-x2), - des moyens de réception (20), reliés aux antennes, et fournissant, à un instant donné, et pour chaque axe d'observation, une valeur numérique (SD1,SD2) de l'amplitude du signal reçu, reliée à la position angulaire des sources audit instant, - des moyens de mémorisation (30) de ces valeurs numériques à des instants successifs, - des moyens de traitement (40,50,60) propres à déterminer des variations de densité de l'histogramme multidimensionnel constitué, dans un repère, par des points dont les coordonnées sont lesdites valeurs numériques mémorisées simultanément, pendant un intervalle de temps choisi, caractérisé en ce que les moyens de traitement (40) recherchent les orientations des contours de l'histogramme dans ledit repère pour déterminer les directions des sources par rapport aux axes d'observation. 2. Dispositif selon la revendication 1, caractérisé en ce que l'histogramme étant bidimensionnel, ses contours sont des segments de droites, et en ce que la recherche comprend la détermination des segments de droites (41,42), puis l'ajustement, sur le contour de l'histogramme, une forme polygonale dont les côtés sont sélectionnés parmi les segments déterminés (43), chaque côté permettant de déterminer la direction d'une source (44). 3. Dispositif selon la revendication 1, caractérisé en ce que l'histogramme étant tridimensionnel, ses contours sont des faces polygonales, et en ce que la recherche comprend le dénombrement desdites faces polygonales, puis l'ajustement, sur le contour de l'histogramme, d'au moins une forme polyédrique ayant le nombre voulu de faces polygonales dont les côtés permettent de déterminer la direction d'une source. 4. Dispositif selon l'une des revendications 2 ou 3, caractérisé en ce que les formes polygonales sont convexes et ont des côtés parallèles deux à deux. 5. Dispositif selon l'une des revendications précédentes, caractérisé en ce que les moyens de traitement comprennent des moyens de traitement d'images (41,42), des moyens de dénombrement et de conversion de pentes (43,44) et des moyens d'exploitation (60) des informations provenant des couples de signaux sur chaque intervalle de temps traité. 6. Dispositif selon la revendication 5, caractérisé en ce que les moyens de traitement d'images comprennent des moyens de calcul des pentes (41) des droites portant lesdits segments. 7. Dispositif selon la revendication 5, caractérisé en ce que les moyens de traitement d'images mesurent la longueur desdits segments (42), le dénombrement des segments (43) tenant compte d'un seuil de longueur. 8. Dispositif selon la revendication 5, caractérisé en ce que les moyens de traitement comprennent des moyens de conversion (44) des pentes desdits segments en autant de directions de sources. 9. Dispositif selon l'une des revendications 5 à 8, caractérisé en ce que les moyens d'exploitation (60) tiennent compte de la distribution des points dans l'histogramme, en particulier pour rechercher le taux de modulation des sources. 10. Dispositif selon l'une des revendications 5 à 9, caractérisé en ce que les moyens d'exploitation (60) ajustent la durée dudit intervalle de temps choisi en fonction des caractéristiques de l'histogramme. 11. Dispositif selon l'une des revendications 5 à 10, caractérisé en ce que les moyens d'exploitation (60) ajustent la durée dudit intervalle de temps choisi en fonction de caractéristiques connues des sources. 12. Dispositif selon l'une des revendications 5 à 11, caractérisé en ce que les moyens d'exploitation (60) ajustent la durée dudit intervalle de temps choisi en fonction du rapport signal à bruit. 13. Dispositif selon l'une des revendications 1 à 12, caractérisé en ce que les signaux d'entrée des moyens de réception proviennent d'antennes multiples selon la technique de formation de voies (Wi,SW). 14. Procédé de localisation de sources radioélectriques, dans lequel: a) à l'aide d'au moins deux antennes (Al, A2), définissant au moins deux axes d'observation (Al-xl, A2-x2), on élabore un couple de valeurs d'amplitude (SD1,SD2) de signaux reçus, l'amplitude étant reliée à la position angulaire des sources à l'instant de réception, b) on construit (20,30) une représentation numérique de l'histogramme multidimensionnel constitué, dans un repère, par des points dont les coordonnées sont lesdites valeurs d'amplitudes, prises sur un intervalle de temps choisi, et c) on procède (40,50,60) à un traitement pour déterminer des variations de densité de cet histogramme, caractérisé en ce que l'étape c) comprend la recherche (40) des orientations des contours de l'histogramme dans ledit repère. 15. Procédé selon la revendication 14, caractérisé en ce que l'étape c) comprend: cl) la détermination des frontières du nuage de points de l'histogramme (41), c2) la sélection de certaines pentes desdites frontières (43), et c3) la transformation des informations de pente desdites frontières sélectionnées en informations de gisement des sources (44). 16. Procédé selon l'une des revendications 14 ou 15, caractérisé en ce que l'étape cl) comprend le calcul de la taille desdites frontières (42). 17. Procédé selon l'une des revendications 14 à 16, caractérisé en ce que les étapes b) et c) sont répétées en faisant varier la durée de l'intervalle de temps choisi et/ou la cadence d'échantillonnage sur cet intervalle de temps.