FR 3091355 A1 20200703 FR 1874115 A 20181227 Description Titre de l'invention : Système et procédé d'identification de source radar Domaine technique de l'invention [0001] La présente invention concerne le domaine des RADAR, des sources de radiation radar et du Renseignement d'Origine Electro-Magnétique (ROEM) d'une manière générale. La présente invention concerne plus particulièrement un procédé et un système d'identification de sources radar. Technique antérieure [0002] Il existe, dans le domaine du ROEM, une difficulté persistante pour identifier les sources radar de façon rapide, efficace et fiable. [0003] Il est connu dans l'art antérieur divers systèmes ou procédés d'identification de sources radar, comme par exemple ceux décrits dans les demandes de brevet CN105631484A, CN101893704B ou CN108416290A. [0004] Les diverses solutions connues sont souvent qualifiées d'intelligence artificielle et utilisent généralement des réseaux de neurones et des méthodes telles que le « deep lemming » (« apprentissage profond » d'après la terminologie anglo-saxonne) ou d'autres méthodes de classification. [0005] Un problème important dans le domaine est lié au signal radar lui-même car il est composé de signaux radar mélangés de diverses sources qu'il n'est pas toujours aisé de discriminer, notamment à cause de divers paramètres tels que la force du signal, la localisation de la source, les perturbations environnantes, etc. [0006] Dans ce contexte, il est intéressant de proposer un système et un procédé permettant d'optimiser l'identification radar. Exposé de l'invention [0007] Un but de la présente invention est de pallier au moins certains inconvénients de l'art antérieur en proposant un système et un procédé d'identification de source radar qui soit rapide, efficace et fiable. [0008] Ce but est atteint par un procédé d'identification de source radar, à partir d'un signal généré par au moins un capteur d'impulsions radar, mis en oeuvre dans un système d'identification radar accédant à au moins une base de données de renseignement (BDR), le procédé comportant au moins un traitement des données par au moins un modèle de réseaux de neurones multi-couches, ledit capteur détectant et caractérisant les impulsions produites par les radars en présence autour du capteur en générant ledit signal sous la forme de descripteurs d'impulsions (PDW), le procédé étant caractérisé en ce qu'il comporte une étape de conditionnement préalable audit traitement, l'étape 2 de conditionnement comportant un regroupement des descripteurs d'impulsions (PDW) par comparaison d'au moins un paramètre parmi au moins la fréquence et la durée des impulsions. [0009] Ce but est également atteint par un système d'identification de source radar, comprenant des moyens informatiques de traitement de données aptes à traiter les données générées par des capteurs radar, caractérisé en ce que les moyens informatiques exécutent des instructions permettant la mise en oeuvre du procédé selon l' invention. [0010] Selon une autre particularité, le regroupement de l'étape de conditionnement est un regroupement paramétrique ou un regroupement densitométrique utilisant les données radar de base pour calculer au moins un des paramètres suivants : [0011] Dispersion fréquence Dispersion durée Dispersion niveau Dispersion de la Période de Répétition des Impulsions (PRI), Histogramme de fréquence, Histogramme de durée, Histogramme de niveau, Histogramme de la Période de Répétition des Impulsions (PRI), Données de position, Modulation la plus occurrente [0012] Selon une autre particularité, le regroupement de l'étape de conditionnement est un regroupement par tuilage utilisant les données radar de base pour calculer au moins un des paramètres suivants : [0013] La fréquence, La durée d'impulsion, Le niveau La différence (DTOA) entre les temps d'arrivée successifs des descripteurs d'impulsion (PDW). [0014] Selon une autre particularité, les modèles sont entraînés sur la base des résultats validés fournis par la base de données de renseignement (BDR). [0015] Selon une autre particularité, ledit modèle de réseaux de neurones est un modèle de type perceptron multi-couche (MLP). [0016] Selon une autre particularité, ledit modèle de réseaux de neurones est un modèle de type réseau de neurones convolutif (CNN). [0017] Selon une autre particularité, ledit modèle de réseaux de neurones est un modèle de type réseau de neurones récurrent (RNN). [0018] Selon une autre particularité, le traitement comporte les deux phases successives 3 suivantes : [0019] Un premier traitement en mode supervisé, c'est-à-dire en recherchant des émetteurs déjà connus et identifiés en bases de données. Un second traitement en mode non-supervisé, c'est-à-dire en recherchant des émetteurs non connus qui ne sont pas référencés en bases de données (BTI ou BDR) [0020] Selon une autre particularité, le traitement supervisé utilise en algorithme d'apprentissage profond (deep learning). [0021] Selon une autre particularité, le traitement non-supervisé utilise un algorithme de classification par groupe (clustering) Brève description des figures [0022] D'autres particularités et avantages de la présente invention apparaîtront plus clairement à la lecture de la description ci-après, faite en référence aux dessins annexés, dans lesquels : [0023] [fig.1] représente une vue d'ensemble d'un procédé selon certains modes de réa- lisation, [0024] [fig.21 illustre les résultats obtenus par un procédé selon certains modes de réa- lisation, [0025] [fig.34 représente une vue d'ensemble de la phase d'apprentissage d'un procédé selon certains modes de réalisation, [0026] [fig.311] représente une vue d'ensemble de la phase d'exploitation d'un procédé selon certains modes de réalisation, [0027] [fig.4] illustre les résultats obtenus par un procédé selon certains modes de réa- lisation, [0028] [fig.51 illustre des étapes d'un procédé selon certains modes de réalisation, [0029] [fig.6] illustre une étape de décision d'un procédé selon certains modes de réa- lisation, [0030] [fig.7] illustre les apprentissages des traitements supervisés et non-supervisés d'un procédé selon certains modes de réalisation, [0031] [fig.84 illustre une étape de conditionnement d'un procédé selon une première variante, [0032] [fig.813] illustre une étape de conditionnement d'un procédé selon une seconde variante, [0033] [fig.8c] illustre des étapes de conditionnement d'un procédé selon une troisième variante, [0034] [fig.8d] illustre des étapes de conditionnement d'un procédé selon une troisième variante, 100351 Ifie.91 illustre un procédé selon certains modes de réalisation. Description des modes de réalisation [0036] La présente invention concerne un système et un procédé d'identification de source radar. D'une manière générale, le système comporte des moyens informatiques exécutant des instructions permettant la mise en oeuvre du procédé. Aucun détail ne sera fourni sur les moyens informatiques bien connus de l'homme de métier, car la présente invention peut être mise en oeuvre sur divers dispositifs ou infrastructure informatique, en réseau ou non. [0037] La présente demande utilise divers acronymes parmi la liste ci-dessous, dans laquelle la définition est fournie pour chacun d'eux, afin de faciliter la compréhension du lecteur. [0038] [Tableaux 1] Terme Définition ROEM Renseignement d'Origine Electro-Magnétique BDR Base de Résultats BDTN Base de Données Techniques Navales BTI Base Technique d'Identification CERES Capacité de Renseignement Electroma2nétique Spatiale PIC Produit d'Interception CERES ELINT ELectronic INTelligence CTE Chaîne dc Traitement ELINT CTE-A Chaîne de Traitement ELINT Automatique CTE-M Chaîne de Traitement ELINT Manuel CTE-OC Chaîne de Traitement ELINT d'Outils Complémentaires CTE-RM Chaîne de Traitement ELINT de Reprise Manuelle IA Intelligence Artificielle IHM Interface Homme Machine IHS Interface Homme Système MODULE Ensemble unitaire et cohérent exécutant une suite d'instruction avec des entrées, sorties et interfaces OUTIL Moyen d'édition , de mesure et de modification des entrées/ sorties des modules DI Durée d'Impulsion DOA Direction Of Arrivai (direction d'arrivée) DTOA Delta Time of Arrivai (différence d'instants d'arrivée) ELCI B Elément Ciblé FN False Negative FP False Positive OC Graphes de Contrôle PRI Période de Répétition des Impulsions PRM Période de Répétition du Motif 15 Impulsion synthétique 6 RNN Recurent Neural Network CNN Convolutional Neural Network MLP Multi Layer Perceptron TIA Type d'Informations Attendues TOA Time of Arrivai (instant d'arrivée) TN Truc Neeative TP Truc Positive SI Système d'Information UT Groupe de Travail 100391 Le procédé de traitement détaillé dans ce document est un procédé de classification des signaux radars. 11 permet à partir de signaux enregistrés par un capteur quelconque de déterminer quels ont été les radars à l'origine de ces signaux. La classification se traduit au minimum par une identification d'une signature radar relativement à un modèle de radar. Si les données le permettent, la classification peut conduire à une identification d'une instance physique de radar. 100401 Les modèles employés pour cette classification supervisée apprennent. de résultats issus d'une base de données (BDR). Cette BDR regroupe un vaste ensemble de signaux radars pour lesquels une classification a déjà été réalisée et validée par des algorithmes ou directement par des opérateurs. D'où la notion d'apprentissage par fusion des expériences : les modèles apprennent des résultats obtenus quelle que soit leur provenance et tirent donc bénéfices de toutes les techniques déjà utilisées dans le domaine. 10041] Les fonctionnalités proposées requièrent la consultation des bases de données (BT1. BDR) d'un système de renseignement radar / ROEM. Elles bénéficient des évolutions de ces bases lors des incréments d'apprentissage. Les résultats obtenus sont réinjectés dans des bases spécifiques différente des BTI/BDR pour garantir une indépendance totale vis-à-vis du système de ROEM. Cependant, l'injection dans les bases de données (BTI. BDR) sur action explicite de l'utilisateur de la chainc manuelle (CTE M) est réalisable. [0042] Les données radar en entrée du traitement sont les données produites par un capteur radar capable de détecter et caractériser les impulsions produites par les radars en présence autour du capteur. Cette détection et caractérisation des impulsions se synthétisent sous la forme de « mots descripteurs d'impulsions » ou PDW (pour l'anglais « Pulse Description Word ») résumant le contenu des impulsions. [0043] Les PDW contiennent à minima les informations suivantes : 7 [0044] - La fréquence de l'impulsion, [0045] - Le niveau de l'impulsion, [0046] - La durée de l'impulsion, [0047] - L'instant d'arrivée de l'impulsion [0048] En fonction des capacités du capteur, les PDW peuvent contenir en supplément les informations suivantes : [0049] - Modulation de l'impulsion, [0050] - Position ou géolocalisation de l'impulsion. [0051] Les données radars se résument donc à une liste de PDW qui grossit au fur-et-à- mesure de l'arrivée au niveau du capteur des impulsions des radars environnants. La précision des informations contenue dans les PDW dépend des performances du capteur et peut varier en fonction des conditions d'interception (un radar lointain émettant à faible puissance sera généralement moins bien perçu par le capteur et donc ses impulsions seront moins bien caractérisées). [0052] Modèles multi-couches [0053] Les données relatives aux signaux radar recueillies par les systèmes ROEM sont gé- néralement traitées par des modèles de réseaux de neurones, par exemple de type MLP, CNN ou RNN. La figure 5 illustre un exemple d'un tel modèle. Le nombre de couches pour chaque modèle n'est pas fixé et à adaptable en fonction des données disponibles pour l'apprentissage. [0054] Par ailleurs, le choix du modèle peut également varier en fonction du paramètre étudié. Par exemple, dans le cas d'un conditionnement par tuilage, les images réalisées sur la fréquence peuvent alimenter des MLP alors que les images réalisées sur le niveau alimenteront des CNN. [0055] Les résultats obtenus par chaque modèle utilisé sont convoyés vers une couche de fusion puis une couche réalisant une pondération pour chaque classe de radar existante dans la base de données BDR. [0056] Les fonctions d'activation des neurones sont définies par les méthodes existantes classiquement fournies dans les librairies de modèles « Deep Learning » et « Machine Leaming » : sigmoïde, ReLu, etc... [0057] Les modèles sont entraînés sur la base des résultats validés fournis par la BDR en utilisant par descente de gradient d'une fonction de coût de type entropie croisée. D'autres types de fonction de coût pourront être utilisées en fonction des données mises à disposition. [0058] Décision [0059] Dans certains modes de réalisation, le procédé comporte une étape de décision générant des données consolidées par un seuil de décision relatif à leur cohérence. [0060] Un exemple, illustré en figure 6, d'une telle étape de décision repose sur un seuillage 8 pour sélectionner les source identifiées qui ont le meilleur taux de correspondances avec les classes radars existantes. [0061] A l'issu du traitement des données capteurs, les modèles employés fournissent en sortie un taux de correspondance à chacune des classes radars existantes dans la hase de données. Si la base contient les enregistrements de N radars différents, alors on obtient en sortie un vecteur de taille N contenant un matching score pour chaque classe radar. [0062] Le matching score peut suffire comme résultat d'interprétation des données par le modèle, cependant pour un opérateur radar il est rapidement fastidieux de visualiser ce résultat lorsque le nombre de classes radars à trouver devient grand. Le matching score est donc post-traité par un fonction de scuillage incluse dans les modèles et entraînée sur les données d'apprentissage. On a donc une couche propre à la décision dans le modèle. [0063] Les métriques possibles pour l'entraînement de cette couche de décision sont : [0064] - Le ratio de résultats exacts, Le rappel (TP / (TP + FN)) et la précision (TP / (TP + FP)) Le FI score (moyenne harmonique de la précision et du rappel). [0065] Avec TP: vrai positif, FN : faux négatif, FP : faux positif [0066] Traitements supervisé et non-supervisé [0067] Dans certains modes de réalisation, le système et le procédé utilisent un traitement en deux phases successives : un traitement supervisé et un traitement non supervisé, comme indiqué sur la figure 1. [0068] Un premier algorithme utilise l'intelligence artificielle en mode supervisé c'est-à-dire en recherchant des émetteurs déjà connus et identifiés en bases de données (BTI ou BDR) . [0069] Le second algorithme, complémentaire du premier, utilise l'intelligence artificielle en mode NON supervisé, c'est-à-dire en recherchant des émetteurs non connus qui ne sont pas référencés en bases de données (BTI ou BDR). [0070] Le traitement supervisé repose de préférence sur au moins un algorithme de type Deep Leaming. Le modèle utilisé réalise de préférence un apprentissage sur la BDR et classifie les impulsions selon les émetteurs référencés. [0071] Une analyse de complétude (vérification que toutes les impulsions ont été attribuées à une source radar) permet de récupérer les éventuelles impulsions non classifiées, afin de les injecter dans le traitement non-supervisé. [0072] Le traitement non supervisé repose de préférence sur au moins un algorithme de type Clustering. Le modèle réalise de préférence un apprentissage sur la BDR et permet de trouver de nouvelle signature dans les impulsions non classifiées. [0073] La figure 2 montre le principe du traitement en deux phases successives, dans 9 l'exemple d'un signal recueilli contenant l'imbrication des impulsions de quatre sources radars, en supposant que deux de ces radars soient référencés en base de données (BTI/BDR) et qu'il y a deux nouveaux radars, inconnus en base de données (BTI/BDR). [0074] La première phase est un traitement supervisé dont le but est de retrouver dans les signaux d'entrée des radars référencés en base de données (BTI/BDR). Cette phase permet de retrouver, dans le signal complexe (en noir, ligne du haut sur la Fig.2) les deux radars connus (un en grisé et l'autre en blanc, dans la ligne du milieu sur la Fig.2). [0075] Ensuite, dans une seconde phase, après avoir éliminé du signal les impulsions rattachées à ces deux radars connus, les impulsions résiduelles sont injectées dans un traitement non supervisé. Cc traitement non supervisé doit répartir de façon pertinente les impulsions en plusieurs radars (deux dans cet exemple), mais ceci est fait sans connaitre à priori ni le nombre, ni le type de radars. [0076] Les impulsions ainsi regroupées servent à caractériser deux nouveaux radars en- registrés en base de données (BDR) (grisé et blanc, ligne du bas sur la Fig.2), par exemple approbation d'un opérateur. [0077] Lors des prochaines interceptions de la même configuration de signal, les quatre radars pourront être détectés de façon optimale en mode supervisé car ils seront tous référencés en base de données (BTI ou BDR) [0078] Apprentissage : [0079] Le schéma présenté en figure 3b montre le traitement dans sa phase d'exploitation après apprentissage, tandis que la figure 3a montre le traitement dans sa phase d'apprentissage. [0080] Dans sa phase d'apprentissage, le traitement est exactement le même hormis le fait qu'une fonction de coût récupère en sortie du modèle les classes de sortie et calcule les variations de paramètres à appliquer au modèle pour minimiser le coût. [0081] Un exemple d'apprentissage est représenté sur la figure 7, pour les traitements supervisé et non-supervisé [0082] Traitement des signaux radar par IA : apprentissage [0083] Apprentissage : Apprentissage à partir des données fournies par la BDR (enregistrements radar dont les regroupements à faire sont étiquetés) Préparation amont : Conditionnement des enregistrements pour passer des impulsions au preclusters (les algos apprennent à classer/regrouper les preclusters) [0084] La figure 9 illustre une comparaison entre, d'une part, l'approche classique (en haut) dans laquelle les opérateurs doivent reconnaître eux-mêmes les signaux, par regroupement, caractérisation et croisement avec une BDR pour fournir une liste d'émetteurs validés et, d'autre part, une approche d'intelligence artificielle comme dans la présente invention. [0085] Dans certains modes de réalisation, il est prévu un algorithme alternatif supplé- mentaire pour la reprise manuelle, permettant à l'opérateur de combiner ces approches. La présente invention tire avantage de cette combinaison possible en apprenant à partir des entrées brutes et des sorties finales d'émetteurs validés, pour intégrer et fusionner le contenu des trois étapes de traitement : regroupement, caractérisation et croisement. [0086] Conditionnement [0087] Dans certains modes de réalisation, le système et le procédé utilisent un pré- traitement des données recueillies, appelé ici « conditionnement », comme indiqué en figure 1. [0088] Le conditionnement utilise au moins une information des PDW et de préférence au moins la fréquence et la durée de l'impulsion. Le conditionnement est une étape qui consiste à mettre en forme les données radars de manière à réduire et à homogénéiser les données qui alimenteront les modèles. Cette mise en forme se traduit par l'élaboration de features qui seront analysées par les modèles pour opérer la classification. [0089] Les impulsions ne sont pas traitées directement de façon brute, mais au contraire, des groupes (ou pré-clusters) sont construits en amont sur analyse des signaux. [0090] Les données conditionnées et leurs caractéristiques (« features ») sont construites à partir des informations de base des PDW : fréquence, niveau, durée, instant d'arrivée, modulation et position. L'intérêt du conditionnement est de soumettre les modèles à des données uniformes, quelles que soient les contextes de capture des données radars. D'autres avantages sont aussi de simplifier le problème du regroupement à quelques entités, de tirer des parti des connaissances métiers radar et de diminuer la complexité des algorithmes utilisés dans les modèles d'identification de source radar. [0091] La figure 4 montre un exemple de résultats d'un traitement utilisant le condi- tionnement selon divers modes de réalisation. [0092] Grâce au conditionnement des données (« pre-clustering »), on obtient une pluralité de potentiels groupes (pre-clusters) qui accélèrent les traitements supervisé et non-supervisé. La ligne du haut montre le signal brut à traiter, la ligne du milieu montrent les divers groupes en diverses nuances de gris et la ligne du bas montrent les quatre sources radar identifiées dans le signal (un noir, un gris, un blanc et un hachuré). [0093] 11 existe plusieurs manières de réaliser le conditionnement dans la classification su- pervisée. On en présente ici trois versions applicables dont : [0094] - Regroupement paramétrique, Regroupement densitométrique, Tuilage. 11 [0095] Ainsi, dans certains modes de réalisation, le conditionnement des données est réalisé par regroupement paramétrique, dont un exemple est représenté sur la figure 8a. [0096] Le regroupement paramétrique consiste à regrouper les PDW de départ en groupes de PDW dont les paramètres sont proches. Chaque groupe contient des PDW dont les paramètres fréquence, durée, modulation et position sont identiques à des tolérances près (les paramètres niveau et instant d'arrivée sont généralement ignorés). Le choix des tolérances est réalisé empiriquement à partir de l'analyse d'un grand nombre de données radar. [0097] Les groupes issus du regroupement paramétrique constituent des émissions radars. Une émission radar est construite à partir d'un nombre variable de PDW en fonction de la dispersion des données radars de base. Chaque émission est caractérisée par de nouvelles métriques calculées sur les PDW qu'elle contient. [0098] Une émission peut ainsi être caractérisée par les informations suivantes : [0099] - Dispersion fréquence : valeur min, valeur rnax, moyenne, écart type, Dispersion durée : valeur min, valeur max, moyenne, écart type, Dispersion niveau : valeur min, valeur max, moyenne, écart type, Dispersion PRI : valeur min, valeur max, moyenne, écart type, Histogramme de fréquence, Histogramme de durée, Histogramme de niveau, Histogramme de PRI, Données de position : position moyenne, ellipse d'incertitude, Modulation la plus occurrente. [0100] En fonction du choix des tolérances lors de la création des émissions, on peut également avoir des informations d'agilité dont : [0101] Moments de fréquence, Moments de durée, Moments de PRI. [0102] Après le conditionnement paramétrique, on passe d'une liste de PDW à une liste d'émissions dont les informations de caractérisation sont utilisées comme features par les modèles. [0103] Le conditionnement par regroupement paramétrique peut être réalisé sur l'ensemble des PDW générés par le capteur ou par tranches temporelles de manière à conserver un cadencement temporel des features. Ce cadencement temporel permet d'utiliser des modèles de type convolutifs. [0104] Dans certains modes de réalisation, le conditionnement des données est réalisé par re- groupement densitométrique dont un exemple est représenté sur la Figure 8b. [0105] Le conditionnement par regroupement densitométrique consiste à créer des groupes 12 de PDW en analysant la densité des PDW dans un espace multidimensionnel où chaque dimension représente un paramètre des PDW. L'analyse de la densité exploitent les paramètres de fréquence, de durée, de modulation et de position (les paramètres niveau et instant d'arrivée sont généralement ignorés). [0106] Contrairement au conditionnement par regroupement paramétrique, il n'y a pas de to- lérances spécifiques à chaque paramètre. La densitométrie permet d'adapter la taille des groupes réalisés en fonction de la densité des PDW dans l'espace multidimensionnel. L'intérêt de cette méthode est de pouvoir créer des groupes s'adaptant aux diverses agilités intrinsèques aux radars et d'avoir de la robustesse sur les imprécisions du capteur. [0107] La condition à la création d'un groupe est que celui-ci contient un nombre suffisant de PDW dans un espace de recherche fini. L'espace de recherche peut être adapté dynamiquement en fonction du nombre de PDW mis en groupes et permettent ainsi d'avoir des groupes de morphologies diverses. La technique s'apparente ici à une fonction de clustering par densité (type DBSCAN ou OPTICS). [0108] Les groupes ici du regroupement densitométrique constituent des clusters radars. Un cluster radar est construit à partir d'un nombre variable de PDW en fonction de la densité des données radars de base dans l'espace multidimensionnel. Chaque cluster est caractérisé par de nouvelles métriques. [0109] Un cluster peut ainsi être caractérisé par les informations suivantes : [0110] Dispersion fréquence : valeur min valeur max, moyenne, écart type, Dispersion durée : valeur min, valeur max, moyenne, écart type, Dispersion niveau : valeur min, valeur max, moyenne, écart type, Dispersion PRI : valeur min valeur max, moyenne, écart type, Histogramme de fréquence, Histogramme de durée, Histogramme de niveau, Histogramme de PRI, Données de position : position moyenne, ellipse d'incertitude, Modulation la plus occurrente. [0111] Dans certains modes de réalisation, le conditionnement des données est réalisé par tuilage, dont un exemple de fonctionnement est illustré sur les figures 8c et 8d. [0112] Le conditionnement par tuilage consiste à obtenir une représentation graphique des PDW de départ sur certains paramètres en réalisant un échantillonnage temporel. Les paramètres exploités par le tuilage sont la fréquence, la durée d'impulsion, le niveau et le DTOA. Le DTOA est la différence entre les temps d'arrivée successifs des PDW. La figure 8c illustre le calcul des grilles de tri pour chaque paramètre utilisé. [0113] Pour réaliser le tunage, il est nécessaire de construire au préalable des grilles de tri 13 pour chaque paramètre d'intérêt. Ces grilles de tri sont calculées empiriquement après analyse d'une grande quantité de données radar et permettent de définir pour chaque paramètre une liste de frontières qui serviront pour trier les PDW. Les frontières ne sont pas régulièrement espacées et sont plus concentrées pour les zones sur-représentées dans les données radars. [0114] La génération de ces frontières pour chaque paramètre est conduite de façon à minimiser l'entropie sur l'ensemble des données radars étudiés et donc d'avoir dans toutes les zones délimitées par les frontières une représentation homogène (nombre de PDW quasi identique). [0115] Une fois les grilles calculées, les données radars (PDW) sont ensuite triées selon ces grilles par tranche temporelle successive. Le découpage en tranche temporelle est fixé par un paramètre de résolution temporelle. On peut alors calculer pour chaque tranche de temps et pour chaque case d'une grille le nombre de PDW correspondants dans les données. La figure 8c illustre le tri pour chaque tranche temporelle. [0116] Une fois le tri réalisé, on obtient pour chaque paramètre une représentation graphique en deux dimensions (paramètre x temps). La valeur du pixel (niveau de gris) représente la population de PDW présents dans le pixel. [0117] De manière à obtenir le même format d'image, les données radars de base sont traitées par bloc, chaque bloc correspondant à une durée fixe, par exemple 10 secondes. Pour chacun des paramètres, la taille des images obtenues ne dépend plus que du nombre de cases de la grille et de la résolution temporelle choisie. Cette taille est fixe quelles que soient les données traitées. [0118] En résumé, on obtient tous les durées blocs (toutes les 10 secondes par exemple), quatre images représentatives des PDW contenus dans la durée bloc : [0119] Une image sur la fréquence, Une image sur la durée, Une image sur le DTOA, Une image sur le niveau. [0120] De manière optionnelle, en faisant varier la résolution temporelle, on peut obtenir une collection d'images plus grande pour chaque bloc. [0121] La présente demande décrit diverses caractéristiques techniques et avantages en référence aux figures et/ou à divers modes de réalisation. L'homme de métier comprendra que les caractéristiques techniques d'un mode de réalisation donné peuvent en fait être combinées avec des caractéristiques d'un autre mode de réalisation à moins que l'inverse ne soit explicitement mentionné ou qu'il ne soit évident que ces caractéristiques sont incompatibles ou que la combinaison ne fournisse pas une solution à au moins un des problèmes techniques mentionnés dans la présente demande. De plus, les caractéristiques techniques décrites dans un mode de réalisation 14 donné peuvent être isolées des autres caractéristiques de cc mode à moins que l'inverse ne soit explicitement mentionné. [0122] Tl doit être évident pour les personnes versées dans l'art que la présente invention permet des modes de réalisation sous de nombreuses autres formes spécifiques sans l'éloigner du domaine d'application de l'invention comme revendiqué. Par conséquent, les présents modes de réalisation doivent être considérés à titre d'illustration, mais peuvent être modifiés dans le domaine défini par la portée des revendications jointes, et l'invention ne doit pas être limitée aux détails donnés ci-dessus. [Revendication 1] [Revendication 2] [Revendication 3] [Revendication 4]