FR 3076029 A1 20190628 FR 1701382 A 20171227 PROCEDE DE TRAITEMENT POUR LA DETECTION, LA CARACTERISATION ET LE SUIVI MULTI-CIBLES ET DISPOSITIF ASSOCIE DOMAINE TECHNIQUE DE L'INVENTION Le domaine technique de l'invention est celui de la détection, de la caractérisation et du suivi multi-cibles. La présente invention concerne en particulier un procédé de traitement de mesures. La présente invention concerne également un dispositif pour la mise en oeuvre d'un tel procédé. ARRIERE-PLAN TECHNOLOGIQUE DE L'INVENTION Les algorithmes classiques de détection, de caractérisation et de suivi multi-cibles nécessitent des données fournies par l'utilisateur pour être fonctionnelles. Par exemple, un algorithme de suivi multi-cibles classique doit connaître le nombre de cibles à suivre. Il existe donc un besoin de détecter, caractériser et suivre plusieurs cibles en même temps sans qu'un utilisateur n'ait besoin de fournir d'informations a priori sur les cibles à suivre. RESUME DE L'INVENTION L'invention offre une solution aux problèmes évoqués précédemment, en permettant de détecter, caractériser et suivre plusieurs cibles dans un environnement à partir de mesures de cet environnement sans qu'un utilisateur ne 25 fournisse de données d'initialisation ou données de mise à jour. Un premier aspect de l'invention concerne un procédé de traitement de mesures obtenues par un système antennaire à un indice de temps t donné, comportant les différentes étapes suivantes : 30 - Au sein d'une première liste de clusters, un cluster comportant des 2 caractéristiques temporelles, des caractéristiques spatiales et des caractéristiques de signal, prédire, à l'indice de temps t, les caractéristiques spatiales de chaque cluster contenu dans la première liste de clusters en fonction des caractéristiques spatiales du cluster à sa dernière mise à jour 5 pour former une deuxième liste de clusters ; - Pour chaque mesure : o Former une troisième liste de clusters comprise dans la deuxième 10 liste, chaque cluster comportant les caractéristiques spatiales prédites précédemment ; o Calculer un score de vraisemblance entre la mesure et chaque cluster de la troisième liste de clusters, le score de vraisemblance 15 étant calculé à partir de deux facteurs, un premier facteur étant une distance entre la mesure et au moins une caractéristique spatiale du cluster et le deuxième facteur étant le résultat d'un traitement statistique des caractéristiques du cluster ; 20 o Sélectionner le score de vraisemblance maximal ; o Comparer le score de vraisemblance maximal avec un seuil prédéterminé : - Si le score de vraisemblance maximal est au-dessus du seuil, mettre à jour les caractéristiques de chaque cluster de la troisième liste de clusters pour former une quatrième liste de clusters intégrant les informations de la mesure ; - Sinon, créer un nouveau cluster à partir de la mesure et placer le nouveau cluster dans la deuxième liste de clusters. 3 Grâce à l'invention, les clusters permettent d'assurer le regroupement des mesures prises au cours du temps et ainsi de réaliser le suivi des cibles ayant généré ces mesures. En effet, un cluster est un ensemble de mesures, créé à partir d'une première mesure et qui évolue au fur et à mesure que de nouvelles mesures sont 5 produites. Le procédé permet donc d'assurer le suivi d'une pluralité de cibles, chaque cible correspondant à un cluster à un instant donné, sans nécessiter d'informations au préalable. On entend par mesures obtenues par un système antennaire à un indice de temps 10 t donné, des mesures obtenues à une date Dt (désignée indifféremment par instant t dans la suite de la description), les mesures obtenues précédemment (i.e. à l'indice de temps t-1) étant obtenues à une date Dm. On notera que les écarts séparant deux dates consécutives ADt= Dt, - Dt-i ne sont pas nécessairement toujours les mêmes. 15 Outre les caractéristiques qui viennent d'être évoquées dans le paragraphe précédent, le procédé selon un aspect de l'invention peut présenter une ou plusieurs caractéristiques complémentaires parmi les suivantes, considérées individuellement ou selon toutes les combinaisons techniquement possibles. 20 Avantageusement, chaque mesure présente des caractéristiques temporelles, des caractéristiques spatiales et de caractéristiques de signal qui peuvent être choisies parmi au moins le groupe formé des caractéristiques suivantes : Pour les caractéristiques temporelles : une date ; 25 Pour les caractéristiques spatiales : une position, une incertitude de position, et éventuellement une direction d'arrivée et un angle d'incertitude associé à la direction d'arrivée ; Pour les caractéristiques de signal : une fréquence centrale, une largeur de bande, une puissance, un rapport signal à bruit, un type de modulation, une 30 polarisation, une durée d'émission et trois probabilités de mode de transmission correspondant à trois modes de transmission, bursté, continu et à évasion de fréquence. 4 Avantageusement, les caractéristiques d'un cluster peuvent être choisies parmi au moins le groupe formé des caractéristiques suivantes : Pour les caractéristiques temporelles : une date de première détection, une 5 date de dernière mise à jour et un nombre de mesures ; Pour les caractéristiques spatiales : une position, une incertitude de position et éventuellement une direction d'arrivée ; Pour les caractéristiques de signal : une fréquence centrale, une largeur de bande, une puissance, un rapport signal à bruit, une durée d'émission, un 10 type de modulation, une polarisation et trois probabilités de mode de transmission correspondant à trois modes de transmission, bursté, continu et à évasion de fréquence. Avantageusement, le procédé selon un premier aspect de l'invention comporte une 15 étape de formation de la première liste de clusters avant l'étape de prédiction pour ne garder que chaque cluster d'une liste de clusters dont la position est dans une zone géographique proche définie à partir des positions des mesures ou de la connaissance de la position courant des antennes du système antennaire. Ainsi, seuls les clusters dont la position est proche de la position des mesures sont 20 conservés ce qui permet de diminuer le temps de calcul en ne traitant que les clusters les plus proches spatialement des mesures. Avantageusement, l'étape de prédiction des caractéristiques spatiales de chaque cluster de la première liste de clusters à l'indice de temps t est réalisée 25 conformément à un schéma de prédiction de Kalman à partir des informations temporelles et spatiales du cluster. Ainsi, la prédiction des caractéristiques spatiales à un instant t dépend de la prédiction des caractéristiques spatiales à l'instant t-1 ce qui permet de faire évoluer le cluster en fonction de son état précédent. 30 Avantageusement, l'étape de formation de la troisième liste de clusters comporte une première étape de filtrage pour ne garder que chaque cluster de la deuxième 5 liste de clusters répondant à des conditions de filtrage de signal basées sur les caractéristiques de signal de la mesure et les caractéristiques de signal du cluster. Ainsi, seuls les clusters dont les caractéristiques de signal sont proches des caractéristiques de signal de la mesure sont conservés ce qui permet de diminuer 5 le temps de calcul en ne traitant que les clusters les plus similaires à la mesure du point de vue du signal. Avantageusement, l'étape de formation de la troisième liste de clusters comporte une deuxième étape de filtrage pour ne garder que chaque cluster de la deuxième 10 liste de clusters répondant à une condition de filtrage spatiale basée sur les caractéristiques spatiales de la mesure et les caractéristiques spatiales du cluster. Ainsi, seuls les clusters dont les caractéristiques spatiales sont proches des caractéristiques spatiales de la mesure sont conservés ce qui permet de diminuer le temps de calcul en ne traitant que les clusters les plus similaires à la mesure du 15 point de vue spatial. Avantageusement, pendant l'étape de mise à jour, les caractéristiques de signal et les caractéristiques spatiales de chaque cluster de la troisième liste de clusters sont mises à jour à partir du score de vraisemblance normalisé du cluster associé à la 20 mesure et au cluster. Ainsi, la mesure fait évoluer chaque cluster en fonction du score de vraisemblance associé au cluster, c'est-à-dire que plus le score de vraisemblance est important, plus la mesure a de l'influence sur l'évolution des caractéristiques du cluster. 25 Un autre aspect de l'invention concerne un dispositif permettant la mise en oeuvre du procédé selon un premier aspect de l'invention comportant une unité de calcul configurée pour réaliser les étapes du procédé. L'invention et ses différentes applications seront mieux comprises à la lecture de la 30 description qui suit et à l'examen des figures qui l'accompagnent. 6 BREVE DESCRIPTION DES FIGURES Les figures sont présentées à titre indicatif et nullement limitatif de l'invention. La figure 1 montre une représentation schématique du procédé selon un premier aspect de l'invention selon un premier mode de réalisation. 5 La figure 2 montre une représentation schématique d'un cluster. La figure 3 montre une représentation schématique d'une mesure. La figure 4 montre l'étape de mise à jour des caractéristiques spatiales des clusters. 10 DESCRIPTION DETAILLEE D'AU MOINS UN MODE DE REALISATION DE L'INVENTION Sauf précision contraire, un même élément apparaissant sur des figures différentes présente une référence unique. 15 Un premier aspect de l'invention concerne un procédé 100 de traitement de mesures 1021 obtenues à partir d'acquisitions réalisées par un système antennaire dans un environnement dans lequel on veut détecter, caractériser et suivre des cibles. Dans la suite de la présente demande, on désignera par cible, un émetteur dont on veut connaître un certain nombre de caractéristiques, comme par exemple 20 la position, la vitesse ou encore la fréquence d'émission. Préférentiellement, les cibles sont mobiles, c'est-à-dire que les cibles ont une vitesse non nulle mais elles peuvent rester immobiles pendant un certain intervalle de temps. On entend par système antennaire un dispositif comportant un élément antennaire capable de capter des mesures de son environnement et d'inférer des positions. Le système 25 antennaire est préférentiellement aéroporté et mobile. Le procédé 100 selon un premier aspect de l'invention comporte plusieurs étapes dont l'enchaînement est représenté à la figure 1. Ces étapes sont mises en oeuvre par le dispositif selon un deuxième aspect de l'invention comportant une unité de 30 calcul configurée pour réaliser les étapes 101 à 107 du procédé 100. 7 A un instant t, un ensemble de mesures est produit. On désigne par ensemble de mesures 1021 une liste datée par une caractéristique temporelle 1022 de mesures comportant des caractéristiques spatiales 1023 et des caractéristiques signal 1024. La date de l'ensemble de mesures est vue comme une caractéristique 5 supplémentaire dite temporelle 1022. Par exemple, l'ensemble de mesures 1021 de la figure 2 comporte une caractéristique temporelle 1022 valant pour deux mesures (ou plus généralement une pluralité de mesures) comportant des caractéristiques spatiales 1023 et signal 1024. Préférentiellement, les caractéristiques temporelles 1022 sont des dates de mesure. Préférentiellement, les caractéristiques 10 spatiales 1023 sont par exemple une position, une incertitude de position, une direction d'arrivée et un angle d'incertitude associé à la direction d'arrivée. Préférentiellement, les caractéristiques de signal 1024 sont par exemple une fréquence centrale, une largeur de bande, une puissance, un rapport signal à bruit, une durée d'émission, un type de modulation, une polarisation et trois probabilités 15 de mode de transmission correspondant à trois modes de transmission : bursté, continu et à évasion de fréquence. On entend par mode de transmission bursté, la transmission d'un signal de puissance nulle à part à sur certains intervalles de temps courts essaimés dans le temps. On entend par mode de transmission continu, la transmission d'un signal continu. On entend par mode de transmission à 20 évasion de fréquence, la transmission d'un signal qui utilise alternativement plusieurs sous-canaux répartis dans une bande de fréquences selon une séquence donnée, par exemple pseudo-aléatoire. On dispose d'une liste de clusters générée à l'instant t-1. La liste de clusters peut 25 être vide ou comporter un ou plusieurs clusters 1011. On désigne par cluster 1011 un ensemble de caractéristiques 1015 pouvant être de trois types : des caractéristiques temporelles 1012, des caractéristiques spatiales 1013 et des caractéristiques de signal 1014 comme illustrées à la figure 3. Par exemple, le cluster 1011 de la figure 3 comporte des caractéristiques temporelles 1012, des 30 caractéristiques spatiales 1013 et des caractéristiques de signal 1014 et regroupe l'historique d'une pluralité de mesures 1021 (ici deux mesures 1021) à partir desquelles sont mises à jour ses caractéristiques 1015. Préférentiellement, les 8 caractéristiques temporelles 1022 sont par exemple une date de création, une date de dernière mise à jour et un nombre de mesures. Préférentiellement, les caractéristiques spatiales 1023 sont par exemple une position, une incertitude de position et éventuellement une direction d'arrivée. Préférentiellement, les 5 caractéristiques de signal 1024 sont par exemple une fréquence centrale, une largeur de bande, une puissance, un rapport signal à bruit, une durée d'émission, un type de modulation, une polarisation et trois probabilités de mode de transmission correspondant à trois modes de transmission : bursté, continu et à évasion de fréquence. Comme évoqué plus haut, un cluster 1011 est également un 10 regroupement de mesures 1021 comme illustré à la figure 4 dans laquelle chaque cluster 1011 englobe une pluralité de mesures 1021 à partir desquelles sont mises à jour ses caractéristiques 1015. Sur la figure 4 est illustrée la mise à jour de la position 1015 de trois clusters 1011 par l'ajout d'une nouvelle mesure 1026. Selon un premier mode de réalisation, la liste de clusters générée à l'instant t-1 est 15 une première liste de clusters 1010. Selon un deuxième mode de réalisation, la première liste de clusters 1010 est une sous-liste de la liste de clusters générée à l'instant t-1. Un filtrage est réalisé pour supprimer chaque cluster 1011 dont la position 1015 n'est pas dans une zone géographique proche des positions des mesures 1021. Il est possible par exemple 20 de calculer le barycentre des positions des mesures 1021. La zone géographique est alors par exemple un cercle ayant pour centre le barycentre et ayant un diamètre par exemple compris entre 200 et 300 km. La première liste de clusters 1010 est alors formée des clusters 1011 de la liste de clusters générée à l'instant t-1 n'ayant pas été supprimés par le filtrage. 25 Il est également possible de coupler le filtrage précédent avec un filtrage temporel qui supprime chaque cluster 1011 de la liste de clusters générée à l'instant t-1 dont la date de dernière mise à jour est trop éloignée de la date courante, ce qui signifie que la cible qui a généré les mesures 1021 regroupées dans le cluster 1011 était issue d'une erreur ou s'est éteinte depuis. 30 Une fois la première liste de clusters 1010 obtenue, l'étape 101 de prédiction des caractéristiques spatiales 1013 de chaque cluster 1011 de la première liste de 9 clusters 1010 est réalisée. En effet, la première liste de clusters 1010 comporte des clusters 1011 dont les caractéristiques spatiales 1013 ont été évaluées à l'instant t1, c'est-à-dire que, par exemple, la position 1015 de chaque cluster 1011 est celle obtenue à l'instant t-1. Comme les cibles à suivre sont mobiles, les caractéristiques 5 spatiales 1013 des clusters 1011 varient dans le temps. Cette étape 101 de prédiction permet de prédire les caractéristiques spatiales 1013 d'un cluster 1011 de la première liste de clusters 1010, comme par exemple sa position 1015, à l'instant t à partir des caractéristiques spatiales 1013 du même cluster 1011 à l'instant t-1. Préférentiellement, cette étape 101 de prédiction est réalisée selon une 10 prédiction de Kalman intégrant la position estimée en t-1, les incertitudes attachées à cette estimation et la différence entre les dates t et t-1. Les clusters 1011 de la première liste de clusters 1010 mis à jour forment une deuxième liste de clusters 1020. 15 Les étapes 102 à 107 suivantes sont réalisées pour chaque mesure 1021 obtenue à l'instant t. Selon un premier mode de réalisation, une troisième liste de clusters 1025 est formée à partir des clusters 1011 de la deuxième liste de clusters 1020 dont les 20 caractéristiques spatiales 1013 ont été prédites à l'étape 101 précédente. Selon un deuxième mode de réalisation, l'étape 102 de formation de la troisième liste de clusters 1025 comporte une première étape de filtrage pendant laquelle les clusters 1011 de la deuxième liste de clusters 1020 répondant à des conditions de 25 filtrage de signal sont ajoutés à la troisième liste de clusters 1025 initialement vide. Les conditions de filtrage de signal sont basées sur les caractéristiques de signal 1014 de chaque cluster 1011 de la deuxième liste de clusters 1020 et sur les caractéristiques de signal 1024 de la mesure 1021 en cours de traitement. Par exemple, les conditions de filtrage de signal pour un cluster 1011 sont choisies 30 parmi au moins le groupe formé des conditions suivantes : - Un mode de transmission du cluster 1011 est identique à un mode de transmission de la mesure 1021, le mode de transmission du cluster 1011 10 étant le mode transmission dont la probabilité est la plus élevée parmi les trois probabilités de mode de transmission du cluster 1011 et le mode de transmission de la mesure 1021 étant le mode transmission dont la probabilité est la plus élevée parmi les trois probabilités de mode de 5 transmission de la mesure 1021. Par exemple, si la probabilité la plus élevée parmi les trois probabilités de mode de transmission du cluster 1011 est la probabilité que le mode de transmission soit bursté, le mode de transmission est considéré bursté pour le cluster 1011 et si la probabilité la plus élevée parmi les trois probabilités de mode de transmission de la mesure 1021 est 10 la probabilité que le mode de transmission soit à évasion de fréquence, le mode de transmission est considéré à évasion de fréquence, le mode de transmission du cluster 1011 n'est alors pas identique au mode de transmission de la mesure 1021 ; La différence entre la durée d'émission de la mesure 1021 et la durée 15 d'émission du cluster 1011 est inférieure à un multiple d'une résolution temporelle associée à l'acquisition de la mesure 1021; Si le mode de transmission de la mesure 1021 est à évasion de fréquence, la fréquence de la mesure 1021 et la fréquence du cluster 1011 sont comprises dans une même bande de fréquence. 20 La différence entre la fréquence de la mesure 1021 et la fréquence du cluster 1011 est inférieure à un multiple d'une résolution fréquentielle associée à l'acquisition de la mesure 1021. A l'issue de l'étape 102 de formation de la troisième liste de clusters 1025, seuls les clusters 1011 de la deuxième liste de clusters 1020 répondant aux conditions de 25 filtrage de signal choisies sont placés dans la troisième liste de clusters 1025. Selon un troisième mode de réalisation, les clusters 1011 sélectionnés à l'issue de la première étape de filtrage ne sont pas directement placés dans la troisième liste de clusters 1025 mais dans une liste intermédiaire de clusters initialement vide. 30 L'étape 102 de formation de la troisième liste de clusters 1025 comporte alors une deuxième étape de filtrage pendant laquelle les clusters 1011 de la liste intermédiaire de clusters répondant à une condition de filtrage spatiale sont ajoutés 11 à la troisième liste de clusters 1025 initialement vide. La condition de filtrage spatiale est basée sur les caractéristiques spatiales 1013 de chaque cluster 1011 de la liste intermédiaire de clusters et sur les caractéristiques spatiales 1023 de la mesure 1021 en cours de traitement. Préférentiellement, la condition de filtrage spatiale est 5 que l'angle entre la direction d'arrivée de la mesure 1021 et la direction d'arrivée du , 3 cluster 1011 ait une valeur inférieure a -2 fois l'angle d'incertitude associé à la direction d'arrivée de la mesure 1021 en cours de traitement. A l'issue de l'étape 102 de formation de la troisième liste de clusters 1025, seuls les clusters 1011 de la liste intermédiaire de clusters répondant à la condition de filtrage 10 spatiale choisie sont placés dans la troisième liste de clusters 1025. Selon un quatrième mode de réalisation, l'étape 102 de formation de la troisième liste de clusters 1025 comporte uniquement la deuxième étape de filtrage qui forme la troisième liste de clusters 1025 à partir des clusters 1011 de la deuxième liste de 15 clusters 1020 répondant à la condition de filtrage spatiale. Une fois la troisième liste de clusters 1025 associée à la mesure 1021 en cours de traitement formée, l'étape 103 de calcul du score de vraisemblance est réalisée. Un score de vraisemblance 1030 est calculée entre la mesure 1021 et chaque cluster 20 1011 de la troisième liste de clusters 1025. Le nombre de scores de vraisemblance 1030 calculés est donc égal au cardinal de la troisième liste de clusters 1025. On entend par cardinal d'une liste, le nombre d'éléments de la liste. Le score de vraisemblance 1030 entre la mesure 1021 et un cluster 1011 de la troisième liste de clusters 1025 est calculé à partir de deux facteurs. 25 Un premier facteur définit une distance entre la mesure 1021 et le cluster 1011. Préférentiellement, cette distance est basée sur les caractéristiques de signal 1024 de la mesure 1021 et les caractéristiques de signal 1014 du cluster 1011. Plus préférentiellement, la distance est calculée à partir de trois sous-facteurs, chaque 30 sous-facteur étant lié à une unique caractéristique de signal 1024 de la mesure 1021. La distance s'exprime alors ainsi : 12 d(m, ci) = sfi x sf2 x sf3 Avec d(m, ci) la distance entre la mesure m et le i-ème cluster 1011 de la troisième liste de clusters 1025, sfi un premier sous-facteur, sf2 un deuxième sous-facteur et 5 sf3 un troisième sous-facteur. Le premier sous-facteur dépend par exemple de la fréquence : si le cluster 1011 a un mode de transmission à évasion de fréquence, il vaut par exemple 1 et si le cluster 1011 a un mode de transmission continu, il vaut : 10 M Sfi = 1 A E Lilmax Avec Afmax une fonction des paramètres d'acquisition de la mesure 1021 et avec : 15 Af = i fm - fci i Avec fm la fréquence de la mesure m 1021 et fci la fréquence du i-ème cluster 1011 de la troisième liste de clusters 1025. 20 Le deuxième sous-facteur dépend par exemple de la largeur de bande : il vaut par exemple la largeur de bande de la mesure 1021. Le troisième sous-facteur dépend par exemple de la durée d'émission : il vaut par exemple : sf3 = 1 - At Avec : I dm - cid I At = dci 25 30 13 Avec dm la durée d'émission de la mesure m 1021 et la durée d'émission du i-ème cluster 1011 de la troisième liste de clusters 1025. Finalement, le premier facteur F1 s'exprime ainsi : 1 Fi (in, ci) = d(m, ci) Le second facteur est le résultat d'un traitement statistique des caractéristiques du cluster 1011. Préférentiellement, le second facteur F2 s'exprime ainsi : exp(-s(m, ci)) n(m, ci) F2 (M, Ci) = exp(-s(m, ck)) Ek n(m, ck) Avec la somme sur tous les clusters 1011 de la troisième liste de clusters 1025 et : 5 10 15 s(m, ci) = n(m, ci) exp(-dist(m, ci)) Avec dist(m, ci) définissant une distance entre la position de la mesure m 1021 et la position prédite à l'étape 101 de prédiction du i-ème cluster 1011 et n(m, ci) permettant la prise en compte des incertitudes de position. 20 Le score de vraisemblance S 1030 s'exprime ainsi : S(m, ci) = Fi (m, ci) X F2 (M, ci) 25 Et finalement, le score de vraisemblance SN normalisé (étape 1031) s'exprime ainsi : SN(m, ci) = Ei S(m, ci) S(m, ci) 14 Une fois que chaque score de vraisemblance 1030 associé à chaque cluster 1011 de la troisième liste de clusters 1025 a été calculé, l'étape 104 de sélection consiste à sélectionner le score de vraisemblance maximal 1040 parmi tous les scores de 5 vraisemblance 1030. Une fois le score de vraisemblance maximal 1040 sélectionné, l'étape 105 de comparaison consiste à comparer le score de vraisemblance maximal 1040 à un seuil 1041 prédéterminé. Si le score de vraisemblance maximal 1040 est au-dessus 10 ou égal au seuil 1041, l'étape 106 de mise à jour des caractéristiques 1015 de chaque cluster 1011 de la troisième liste de clusters 1025 est réalisée. Sinon l'étape 107 de création d'un nouveau cluster 1060 est réalisée. L'étape 106 de mise à jour des caractéristiques 1015 de chaque cluster 1011 de la 15 troisième liste de clusters 1025 consiste à mettre à jour les caractéristiques spatiales 1013 et les caractéristiques de signal 1014 de chaque cluster 1011 à partir de son score de vraisemblance 1030, à mettre à jour la date de dernière mise à jour du cluster 1011 avec la date actuelle et à incrémenter le nombre de mesures du cluster 1011 del. 20 Préférentiellement, la mise à jour des caractéristiques spatiales 1013 de chaque cluster 1011 de la troisième liste de clusters 1025 est réalisée selon un schéma de correction de Kalman dont le gain est multiplié par le score de vraisemblance normalisé 1031 associé au cluster 1011 qui est en train d'être mis à jour. Cette mise à jour des caractéristiques spatiales 1014 est illustrée à la figure 4. Sur la figure 4, 25 trois clusters 1011 sont représentés : chaque cluster 1011 regroupe une pluralité de mesures 1021 et possède une position 1015. L'ajout d'une nouvelle mesure 1026 entraîne la mise à jour de la position 1015 de chaque cluster 1011. Préférentiellement, toutes les caractéristiques de signal 1014 d'un cluster 1011 sont mises à jour à partir de son score de vraisemblance normalisé 1031. Par exemple, 30 la fréquence du cluster 1011 mis à jour s'exprime ainsi : f1 (t) = fci(t - 1) + SN (m, ci) x fn, 15 Avec f1(t) la fréquence mise à jour du i-ème cluster 1011 de la troisième liste de clusters 1025, fci (t - 1) la fréquence avant mise à jour du i-ème cluster 1011 de la troisième liste de clusters 1025, SN(m, ci) le score de vraisemblance normalisé 1031 5 associé à la mesure m 1021 et au i-ème cluster 1011, fm la fréquence de la mesure m 1021 et N le nombre de mesures du i-ème cluster 1011. Les clusters 1011 de la troisième liste de clusters 1025 mis à jour sont ajoutés à une quatrième liste de clusters 1050 initialement vide. 10 L'étape 107 de création d'un nouveau cluster 1060 consiste à donner comme valeurs aux caractéristiques spatiales 1013 et aux caractéristiques de signal 1014 du nouveau cluster 1060 respectivement les valeurs des caractéristiques spatiales 1023 et des caractéristiques de signal 1024 de la mesure 1021. Par exemple, la position 1015 du nouveau cluster 1060 prend la valeur de la position de la mesure 15 1021. Le nombre de mesures du cluster 1011 prend la valeur 1 et la date de dernière mise à jour est initialisée à 0. Les autres caractéristiques temporelles 1012 du nouveau cluster 1060 prennent les valeurs des caractéristiques temporelles de la mesure 1021. Le nouveau cluster 1060 est ensuite placée dans la deuxième liste de clusters 1020. 20 Après avoir réalisé toutes les étapes du procédé 100 pour toutes les mesures 1021 obtenues à l'instant t, chaque cluster 1011 donne des informations sur une cible comme par exemple sa position, sa vitesse ou encore ses caractéristiques d'émission. Par exemple, à chaque instant, la position de chaque cluster 1011 25 donne la position de chaque cible de l'environnement. REVENDICATIONS1. Procédé (100) de traitement de mesures (1021) obtenues par un système antennaire à un indice de temps t donné, comportant les différentes étapes suivantes : - Au sein d'une première liste de clusters (1010), un cluster (1011) comportant des caractéristiques temporelles (1012), des caractéristiques spatiales (1013) et des caractéristiques de signal (1014), prédire, à l'indice de temps t, les caractéristiques spatiales (1013) de chaque cluster (1011) contenu dans la première liste de clusters (1010) en fonction des caractéristiques spatiales (1013) du cluster (1011) à sa dernière mise à jour pour former une deuxième liste de clusters (1020, 101) ; - Pour chaque mesure (1021) : o Former une troisième liste de clusters (1025) comprise dans la deuxième liste, chaque cluster (1011) comportant les caractéristiques spatiales prédites précédemment (102) ; o Calculer un score de vraisemblance (1030) entre la mesure (1021) et chaque cluster (1011) de la troisième liste de clusters (1025, 103), le score de vraisemblance (1030) étant calculé à partir de deux facteurs, un premier facteur étant une distance entre la mesure (1021) et au moins une caractéristique spatiale du cluster (1011) et le deuxième facteur étant le résultat d'un traitement statistique des caractéristiques (1015) du cluster (1011). o Sélectionner le score de vraisemblance maximal (1040, 104). o Comparer le score de vraisemblance maximal (1040) avec un seuil prédéterminé (1041, 105) : 17 - Si le score de vraisemblance maximal (1040) est au-dessus du seuil (1041), mettre à jour les caractéristiques (1015) de chaque cluster (1011) de la troisième liste de clusters (1025) pour former une quatrième liste de clusters (1050, 106) 5 intégrant les informations de la mesure ; - Sinon, créer un nouveau cluster (1060) à partir de la mesure (1021) et placer le nouveau cluster (1060) dans la deuxième liste de clusters (1020, 107). 10 2. Procédé selon la revendication 1, caractérisé en ce que chaque mesure (1021) présente des caractéristiques temporelles (1022), des caractéristiques spatiales (1023) et de caractéristiques de signal (1024) qui peuvent être choisies parmi au moins le groupe formé des caractéristiques 15 suivantes : Pour les caractéristiques temporelles (1022) : une date ; Pour les caractéristiques spatiales (1023) : une position, une incertitude de position, et éventuellement une direction d'arrivée et un angle d'incertitude associé à la direction d'arrivée ; 20 Pour les caractéristiques de signal (1024) : une fréquence centrale, une largeur de bande, une puissance, un rapport signal à bruit, une durée d'émission, un type de modulation, une polarisation et trois probabilités de mode de transmission correspondant à trois modes de transmission, bursté, continu et à évasion de fréquence. 25 3. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que les caractéristiques (1015) d'un cluster (1011) peuvent être choisies parmi au moins le groupe formé des caractéristiques suivantes : 30 Pour les caractéristiques temporelles (1012) : une date de première détection, une date de dernière mise à jour et un nombre de mesures ; Pour les caractéristiques spatiales (1013) : une position, une incertitude de 18 position et éventuellement une direction d'arrivée ; Pour les caractéristiques de signal (1014) : une fréquence centrale, une largeur de bande, une puissance, un rapport signal à bruit, une durée d'émission, un type de modulation, une polarisation et trois probabilités de 5 mode de transmission correspondant à trois modes de transmission, bursté, continu et à évasion de fréquence. 4. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le procédé (100) comporte une étape de formation de 10 la première liste de clusters (1010) avant l'étape (101) de prédiction pour ne garder que chaque cluster (1011) d'une liste de clusters (1010) dont la position est dans une zone géographique proche définie à partir des positions des mesures (1021) ou de la connaissance de la position courant des antennes du système antennaire. 15 5. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'étape (101) de prédiction des caractéristiques spatiales (1013) de chaque cluster (1011) de la première liste de clusters (1010) à l'indice de temps t est réalisée conformément à un schéma de 20 prédiction de Kalman à partir des informations temporelles et spatiales du cluster. 6. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'étape (102) de formation de la troisième liste de 25 clusters (1025) comporte une première étape de filtrage pour ne garder que chaque cluster (1011) de la deuxième liste de clusters (1020) répondant à des conditions de filtrage de signal basées sur les caractéristiques de signal (1024) de la mesure (1021) et les caractéristiques de signal (1014) du cluster (1011). 30 7. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'étape (102) de formation de la troisième liste de 19 clusters (1025) comporte une deuxième étape de filtrage pour ne garder que chaque cluster (1011) de la deuxième liste de clusters (1025) répondant à une condition de filtrage spatiale basée sur les caractéristiques spatiales (1023) de la mesure (1021) et les caractéristiques spatiales (1013) du cluster 5 (1011). 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que pendant l'étape (105) de mise à jour, les caractéristiques de signal (1014) et les caractéristiques spatiales (1013) de 10 chaque cluster (1011) de la troisième liste de clusters (1025) sont mises à jour à partir du score de vraisemblance normalise (1030, 1031) du cluster (1011) associé à la mesure et au cluser (1011). 9. Dispositif permettant la mise en oeuvre du procédé (100) selon l'une 15 quelconque des revendications précédentes, caractérisé en ce que le dispositif comporte une unité de calcul configurée pour réaliser les étapes du procédé.