Brevet FR3050279A1 — publié le 20171020

Calcul automatique d'une dimension d'une plateforme mobile La présente invention concerne un procédé de calcul automatique d'une dimension d'une plateforme mobile, ladite dimension étant une dimension de longueur ou une dimension de largeur, et un système associé. L'invention se situe dans le domaine de la surveillance, en particulier de la surveillance maritime lorsque la plateforme considérée est un navire. Dans le contexte de la surveillance maritime, les navires fournissent des informations auto-déclaratives qui donnent l'identification, la position, la route et les dimensions (longueur et largeur) d'un navire, selon un système d'identification automatique, connu sous l'acronyme AIS (« Automated Identification System »). Il est utile d'être en mesure de vérifier ces informations fournies par un navire, afin de détecter des navires suspects, dans un contexte civil ou militaire. A cet effet, on connaît l'utilisation d'informations fournies par un système radar. De manière connue, un système radar comporte un émetteur apte à émettre des ondes électromagnétiques radio et, un récepteur d'ondes électromagnétiques radio réfléchies, et des moyens de calcul permettant de déterminer, la puissance rétrodiffusée par une cible en fonction de la distance entre cette dernière et le radar. Un tel mode de fonctionnement d'un système radar est connu sous le nom de mode « range-profile ». La définition du range profile correspond à la représentation de la puissance rétrodiffusée par la plateforme détectée en fonction de la distance. Le Range Profile d'une cible ne peut s'établir qu'à partir du moment où elle a été détectée auparavant à l'aide d'un mode de détection. Dans l'état de la technique, il est connu de fournir un signal de type « range- profile » à un opérateur, sur une interface graphique de visualisation, l'opérateur étant alors apte à déplacer des curseurs sur l'interface graphique pour effectuer une estimation manuelle d'une dimension, par exemple la longueur, du navire observé. En effet, en fonction du positionnement du système radar par rapport à l'axe du navire observé, le système radar permet de détecter la longueur ou la largeur du navire, ou plus généralement de la plateforme observée. Il est possible alors pour l'opérateur de comparer la dimension calculée à une dimension correspondante, obtenue à partir de l'identifiant AIS reçu du navire observé. Cependant, une telle opération manuelle prend du temps à l'opérateur, et de plus peut présenter des imprécisions. En plus, dans le contexte maritime, les mesures radar sont bruitées à cause de l'environnement (présence de fouillis de mer), et donc il est difficile pour un opérateur humain d'exploiter efficacement un signal bruité affiché sur une interface graphique. 3050279 2 Il est donc souhaitable d'automatiser l'estimation de dimensions d'une plateforme mobile telle qu'un navire, pour éviter les aléas précités. De manière plus générale, dans le domaine de la surveillance à distance, il est souhaitable de pouvoir estimer automatiquement des dimensions d'une plateforme en 5 mouvement, notamment dans des milieux dans lesquels les mesures radar sont bruitées. A cet effet, selon un premier aspect, l'invention propose, un procédé de calcul automatique d'une dimension d'une plateforme mobile, ladite dimension étant une dimension de longueur ou une dimension de largeur, comportant une étape d'émission, par un système radar, d'ondes radioélectriques en direction de ladite plateforme selon un 10 axe d'émission prédéterminé et d'acquisition d'au moins un signal numérique de profile de puissance représentatif d'une puissance de signal réfléchi reçu en fonction d'une distance radiale selon ledit axe d'émission par rapport à un point de référence dudit système radar. Le procédé comporte des étapes de : - application d'un opérateur de filtrage sur le signal numérique de profil de 15 puissance acquis permettant d'obtenir un signal numérique filtré, - détermination par calcul d'une première distance radiale correspondant à un premier pic de variation du signal numérique filtré et d'une deuxième distance radiale correspondant à un deuxième pic de variation du signal numérique filtré, - calcul d'une dimension radiale de la plateforme en fonction desdites première et 20 deuxième distance radiale. Avantageusement, le procédé de l'invention permet, grâce à l'application d'un opérateur de filtrage, d'obtenir un signal filtré qui permet de minimiser l'influence du bruit et du fouillis de mer, et de traiter automatiquement ce signal filtré pour obtenir une dimension de la plateforme observée. 25 Le procédé de calcul automatique d'une dimension d'une plateforme mobile selon l'invention peut présenter une ou plusieurs des caractéristiques ci-dessous, prises indépendamment ou dans toutes leurs combinaisons techniquement acceptables. Il comprend l'acquisition d'une pluralité de signaux numériques de profil de puissance, et les étapes d'application d'un opérateur de filtrage et de détermination d'une 30 première et d'une deuxième distance radiale sont appliquées pour chaque signal numérique de profil de puissance acquis. Le calcul d'une dimension radiale de la plateforme comprend le calcul d'une valeur de longueur radiale estimée pour chaque signal numérique de profil de puissance acquis en effectuant une différence entre la deuxième distance radiale et la première distance 35 radiale, et la mémorisation de la valeur de longueur radiale estimée à partir de chaque signal numérique de profil de puissance acquis. 3050279 3 Le procédé comprend un calcul d'une longueur radiale estimée finale à partir des valeurs de longueur radiale estimées mémorisées. La longueur radiale estimée finale est égale à la médiane des valeurs de longueur 5 radiale estimées mémorisées. Le procédé comprend une étape de validation statistique des valeurs de longueur radiale estimées mémorisées, et une étape d'attribution d'un score de confiance associé à la longueur radiale estimée finale. L'étape de validation statistique comprend un calcul d'histogramme des valeurs de 10 longueur radiale estimées, et un calcul d'au moins un paramètre de dispersion associé audit histogramme. La première distance radiale est la distance radiale la plus petite correspondant à un pic de variation du signal numérique filtré, et la deuxième distance radiale est la distance radiale la plus grande correspondant à un pic de variation du signal numérique 15 filtré. L'étape de détermination par calcul d'une première distance radiale correspondant à un premier pic de variation du signal numérique filtré et d'une deuxième distance radiale correspondant à un deuxième pic de variation du signal numérique filtré comprend un calcul, effectué par application d'une fenêtre glissante, d'une valeur d'intégrale discrète 20 partielle du signal numérique filtré compris dans ladite fenêtre glissante et une comparaison de la valeur d'intégrale partielle à un seuil de variation prédéterminé. Le seuil de variation prédéterminé est calculé en fonction de l'intégrale discrète totale du signal numérique filtré. Le procédé comporte en outre une étape de calcul d'une dimension de longueur 25 ou de largeur en fonction de la dimension radiale calculée et une étape d'estimation de corrélation entre la dimension de longueur ou de largeur calculée et une dimension de longueur ou de largeur prévue, et une levée d'alerte en fonction du résultat de l'étape d'estimation de corrélation. Selon un autre aspect, l'invention concerne un programme d'ordinateur 30 comportant des instructions logicielles qui, lorsqu'elles sont mises en oeuvre par un dispositif programmable, mettent en oeuvre un procédé de calcul automatique d'une dimension d'une plateforme mobile, ladite dimension étant une dimension de longueur ou une dimension de largeur, tel que brièvement décrit ci-dessus. Selon un autre aspect, l'invention concerne un support d'enregistrement 35 d'informations, comportant des instructions logicielles pour l'exécution d'un procédé de 3050279 4 calcul automatique d'une dimension d'une plateforme mobile tel que brièvement décrit ci-dessus, lorsque ces instructions sont exécutées par un dispositif programmable. Selon un autre aspect, l'invention concerne un système de calcul automatique d'une dimension d'une plateforme mobile apte à mettre en oeuvre un procédé tel que 5 brièvement décrit ci-dessus. Ce système comporte une unité de calcul et de traitement comportant un processeur et une ou plusieurs unités de mémoire aptes à stocker des données et des instructions de code de programme, mettant en oeuvre : -un module d'application d'un opérateur de filtrage sur le signal numérique de 10 profil de puissance acquis permettant d'obtenir un signal numérique filtré, -un module de détermination par calcul d'une première distance radiale correspondant à un premier pic de variation du signal numérique filtré et d'une deuxième distance radiale correspondant à un deuxième pic de variation du signal numérique filtré, et 15 -un module de calcul d'une dimension radiale de la plateforme en fonction desdites première et deuxième distance radiale. D'autres caractéristiques et avantages de l'invention ressortiront de la description qui en est donnée ci-dessous, à titre indicatif et nullement limitatif, en référence aux 20 figures annexées, parmi lesquelles : - la figure 1 est une vue schématique de dessus d'un cas d'utilisation de la présente invention, comprenant une plateforme et une antenne d'émission/réception de dispositif radar ; - la figure 2 est un exemple de signal radar acquis avant et après filtrage de 25 débruitage ; - la figure 3 est un synoptique des principales étapes d'un procédé de calcul de dimension d'une plateforme selon un mode de réalisation ; - la figure 4 est un synoptique des principales étapes mises en oeuvre pour l'estimation automatique de distances radiales correspondant à une plateforme dans un 30 mode de réalisation ; - les figures 5 à 7 illustrent des histogrammes de valeurs de longueur radiale estimées. La figure 1 illustre schématiquement en vue de dessus une plateforme 2, par 35 exemple un navire en cours de navigation en mer, dont on cherche à évaluer une des dimensions parmi sa longueur L et sa largeur I. 3050279 5 La plateforme 2 est schématisée, et représentée sous forme d'un rectangle en vue du dessus. On définit la longueur L et la largeur I d'une plateforme réelle par les dimensions correspondantes d'un rectangle circonscrit à la plateforme, en vue du dessus. 5 Lorsque la plateforme 2 est en mouvement, ce mouvement est caractérisé par un vecteur de vitesse P> . Une antenne 4 d'émission/réception d'ondes radio électromagnétiques, faisant partie d'un système radar 8, est apte à émettre des ondes radio électromagnétiques en direction de la plateforme 2 et à recevoir des ondes réfléchies par toute structure 10 réfléchissante rencontrée. Dans l'exemple de la figure 1, l'axe d'émission ou direction de pointage de l'antenne 4 est référencé Z. L'axe Z forme avec la direction du vecteur de vitesse T7 de la plateforme 2 un angle O. 15 L'antenne 4 est couplée à un dispositif radar 6, l'ensemble des éléments 4 et 6 formant le système radar 8. De manière connue, le système radar 8 est adapté à déterminer la valeur du vecteur de vitesse e de la plateforme 2 projeté sur l'axe Z. Le système radar 8 est adapté à acquérir un signal numérique dit « range profile », 20 représentant la puissance du signal réfléchi en réponse à un signal radio électromagnétique émis, en fonction de la distance par rapport à un point de référence de l'antenne 4. Le dispositif radar 6 comporte notamment des modules émetteur/récepteur 10, des modules 12 de modulation/démodulation des signaux, une unité 14 de calcul et de 25 traitement comportant un processeur et apte à effectuer des calculs lorsque le dispositif 6 est mis sous tension, une ou plusieurs unités de mémoire 16 aptes à stocker des données et des instructions de code de programme. Optionnellement, le dispositif radar 6 comporte également ou est connecté à une interface homme-machine 18, comportant un écran de visualisation et des moyens 30 d'interaction avec un opérateur. Le système radar 8 est par exemple embarqué à bord d'un porteur de surveillance (non représenté), qui peut être fixe ou en mouvement. Dans un mode de réalisation, une partie des modules du dispositif 6 est implémentée indépendamment, par exemple l'interface homme-machine, les unités de 3050279 6 mémoire 16 et l'unité de calcul 14 sont externes au dispositif 6 et font partie d'un dispositif programmable dédié, par exemple un ordinateur. Dans un mode de fonctionnement, le dispositif radar 6 fournit un signal radio électromagnétique à émettre selon la direction Z de pointage de l'antenne 4, et en 5 fonction du signal réfléchi obtenu par le module récepteur 10, génère un signal numérique représentatif de la puissance du signal reçu par l'antenne en fonction d'une distance radiale, appelé ci-après signal de profil de puissance acquis. La figure 2 montre un exemple de signal numérique de profil de puissance acquis Sl, représentatif de la puissance normalisée (axe des ordonnées) du signal reçu par 10 l'antenne 4 en fonction de la distance radiale d (axe des abscisses), entre un point de référence de l'antenne 4 et un point spatial de réflexion sur un axe parallèle à l'axe Z. Le signal numérique de profil de puissance acquis est bruité et une détection directe d'une dimension radiale de la plateforme à partir de ce signal bruité risque d'être erronée. 15 Le signal numérique de profil de puissance S1 comporte des zones dans lesquelles il a une faible amplitude, correspondant à du bruit et du fouillis de mer et des zones d'augmentation d'amplitude, correspondant à des obstacles réfléchissants sur le trajet des ondes électromagnétiques émises. En particulier dans le cas où la plateforme 2 est un navire, les augmentations de 20 puissance correspondent à des éléments de structure du navire, qui s'élèvent selon la direction Y perpendiculaire à la direction Z par rapport à la structure de base d'un navire qui est au niveau de la mer. Par exemple, des éléments de structure métalliques réfléchissant s'élèvent à la proue et à la poupe d'un navire, permettant d'indiquer les extrémités du navire selon l'axe 25 Z d'observation et d'en déduire la dimension radiale correspondant. Avantageusement, selon l'invention, le signal numérique de profil de puissance acquis est filtré par un opérateur de filtrage permettant de minimiser l'influence du bruit comme expliqué en détail ci-après. La figure 2 illustre schématiquement une version débruitée par filtrage médian Sr 30 correspondant au signal numérique de profil de puissance acquis S. Le signal Sr est représenté dans le même référentiel que le signal numérique de profil de puissance acquis S1, comportant la puissance normalisée du signal sur l'axe des ordonnées et la distance radiale par rapport à un point de référence de l'antenne 4 sur l'axe des abscisses. 3050279 7 La figure 3 est un synoptique des principales étapes d'un procédé de calcul automatique d'une dimension parmi une dimension de longueur et une dimension de largeur d'une plateforme mobile selon un premier mode de réalisation. Ce procédé est par exemple mis en oeuvre par un logiciel apte à être exécuté par 5 un processeur d'un dispositif programmable lorsqu'il est mis sous tension. Lors d'une première étape 30, un nombre N de signaux numériques de profil de puissance est acquis par le système radar 8. De préférence, dans le cas d'application de mesure automatique de dimension de navire dans le domaine maritime, le nombre N est strictement supérieur à 1, et de 10 préférence supérieur à 10. En effet, l'acquisition des signaux étant bruitée, il est avantageux de disposer d'une pluralité de signaux pour obtenir une estimation fiable. On obtient ainsi un ensemble de signaux numériques de profil de puissance S'15i5.N. Chaque signal de profil de puissance S, est composé d'échantillons S (k) , pour k 15 compris entre 1 et P, P étant la longueur du signal S,, chaque indice k correspondant à une distance radiale par rapport à un point de référence de l'antenne comme expliqué ci-dessus. L'étape 30 d'acquisition est suivie d'une étape 32 de calcul d'une dimension radiale de la plateforme observée à partir de chaque signal numérique de profil de 20 puissance acquis. L'étape 32 est elle-même composée de plusieurs étapes. D'abord, un opérateur de filtrage de réduction de bruit est appliqué au signal numérique de profil de puissance acquis à l'étape 34. Dans le cas d'application à la mesure automatique des dimensions d'un navire, le « bruit » désigne le fouillis de mer et 25 le bruit thermique du RADAR. Le filtrage de réduction de bruit est, par exemple, un filtrage par un filtre médian, dont le principe est de définir une fenêtre glissante de taille r choisie, de considérer le sous-signal intercepté par la fenêtre et de remplacer l'échantillon central de ce sous-signal par la valeur médiane du sous-signal. 30 De préférence, r est un nombre entier impair. A l'issue de l'étape de filtrage, on obtient un signal numérique filtré Sr de P échantillons 1S ,F (k),1 5_ k 5 Pl. Ainsi, les valeurs d'amplitude extrêmes, correspondant à des perturbations de mesure, sont éliminées. 3050279 8 Plus la taille r de la fenêtre de filtrage est grande, plus le signal est lissé, ce qui permet d'atténuer le bruit, mais au risque de perdre des informations sur la cible détectée. De préférence, r est inférieur ou égal à 10% de P. Par exemple, r=15. Avantageusement, le filtrage médian est un filtrage d'implémentation simple, qui 5 ne nécessite pas de ressources calculatoires importantes et pour lequel le temps de calcul est court. En alternative, on applique un autre opérateur de filtrage que le filtrage médian, par exemple un filtre moyenne. Ensuite, lors d'une étape 36, on détecte les deux extrémités de la plateforme 10 surveillée selon l'axe d'observation, ces extrémités correspondant par hypothèse à des variations importantes du signal numérique filtré issu de l'étape 34. Plus précisément, on détermine des distances radiales d'extrémité, à savoir une première distance radiale dri correspondant à un premier pic de variation du signal numérique filtré et une deuxième distance radiale dr2 correspondant à un deuxième pic de 15 variation du signal numérique filtré. La première distance dr1 est la distance radiale la plus petite correspondant à un pic de variation du signal numérique filtré, et la deuxième distance dr2 est la distance radiale la plus grande correspondant à un pic de variation du signal numérique filtré. Pour cela, un calcul par intégrale discrète et seuillage est mis en oeuvre dans un 20 mode de réalisation illustré plus en détail à la figure 4. La figure 4 décrit les principales étapes de la détermination des première et deuxième distances radiales selon un mode de réalisation de l'invention, pour un signal numérique filtré SiF de P échantillons, chaque valeur d'échantillon étant égale à une puissance reçue normalisée pour une distance radiale du point de référence de l'antenne 25 précédemment mentionné. Lors d'une première étape 60, l'intégrale discrète hett(SiF ) du signal numérique SiF est calculée : Int,(S,F )= E S iF (k) k=1 Ensuite, une taille de fenêtre de calcul d'intégrale partielle est définie à l'étape 62, 30 en fonction de la longueur P du signal numérique filtré : x De préférence, x est inférieur ou égal à 10% de P. M=L X I' 100 3050279 9 Lors d'une étape 64 de détermination de seuil de variation, un seuil de variation Intseuili est fixé. Dans un mode de réalisation, le seuil de variation Int si epi, est une valeur constante Int seuil obtenue à partir d'une mémoire ou fournie par un opérateur, et utilisée pour chaque 5 signal de profil de puissance Si. En variante, le seuil de variation Int euil est calculé à partir de la valeur d'intégrale discrète Int,(Sr ) : Intiset'.,=-L- 100xInt,(S F ) où y est un nombre réel. De préférence, y est inférieur ou égal à 10% de P. 10 A l'étape 66 d'initialisation, un indice de boucle n est initialisé à 1, et les valeurs de distances radiales dri et dr2 sont initialisées à 0. Il est vérifié si n est inférieur à P-M à l'étape suivante 68, et en cas de vérification négative, le procédé prend fin. En cas de vérification positive, l'intégrale discrète partielle, égale à l'intégrale 15 discrète du signal numérique filtré compris dans la fenêtre de M échantillons partant de l'échantillon d'indice n, est calculée à l'étape de calcul d'intégrale partielle 70 : n+M Int p(S,F ,n)=IS,F (k) 1c=n On vérifie ensuite à l'étape 72 si l'intégrale discrète partielle calculée dépasse le seuil de variation Int ,i,' préalablement calculé. 20 En cas de vérification positive à l'étape 72, on vérifie si dri=0 à l'étape 74, et en cas de vérification positive à l'étape 74, on passe à l'étape 76 de détermination de la première distance radiale à laquelle on met dri à la distance radiale correspondant à l'indice pi = n +[ indice ndice correspondant au milieu de la fenêtre de taille M de calcul de 2 l'intégrale discrete partielle. 25 En cas de vérification négative à l'étape 74, on passe à l'étape 78 de détermination de la deuxième distance radiale à laquelle met la valeur dr2 à la distance radiale correspondant à l'indice pi = n+[ 2 I - En cas de vérification négative a l'étape 72, cette étape 72 est suivie de l'étape 80 d'incrémentation de l'indice n de un, elle-même suivie de l'étape 68 précédemment 30 décrite. 3050279 10 En outre, les étapes 76 et 78 sont toutes deux suivies de l'étape 80 d'incrémentation de l'indice n. Ce procédé permet de déterminer la première distance dr1 correspondant à la première valeur de variation dépassant le seuil de variation rencontrée en traitant les 5 données dans l'ordre croissant des distances radiales, et la deuxième distance dr2 correspondant à la dernière valeur de variation dépassant le seuil de variation rencontrée en traitant les données dans l'ordre croissant des distances radiales. Bien entendu, diverses variantes d'implémentation à la portée de l'homme du métier sont envisageables pour le calcul des première et deuxième distances radiales. 10 De retour à la figure 3, on mémorise à l'issue de l'étape 36 la première distance dri correspondant à une première extrémité de la plateforme, et la deuxième distance dr2 correspondant à une deuxième extrémité de la plateforme. On en déduit, à l'étape 38 d'estimation d'une dimension de longueur radiale de la plateforme visée, une valeur LR à partir du signal numérique Si par le calcul simple 15 suivant : La longueur radiale estimée LR est mémorisée dans une mémoire du dispositif programmable mettant en oeuvre le procédé de l'invention. Les étapes 34 à 38 sont effectuées pour chaque signal numérique de profil de 20 puissance acquis Si, donc l'étape 38 est suivie de l'étape 34 avec une valeur de i incrémentée de un tant que i reste inférieur ou égal au nombre d'acquisitions N. Une fois tous les signaux numériques de profil de puissance acquis traités, l'étape 38 est suivie d'une étape 40 de calcul d'une longueur radiale estimée finale à partir des valeurs de longueur radiale estimées, et d'une étape 42 de validation statistique de la 25 longueur radiale estimée finale. Lors de l'étape 40, la longueur radiale estimée finale LR est calculée comme étant égale à la valeur médiane des valeurs de longueur radiale estimées préalablement mémorisées. En variante, on utilise un autre estimateur statistique pour déterminer la longueur 30 radiale estimée finale LR à partir des longueurs radiales estimées LR , par exemple la moyenne. Ensuite, lors de l'étape de validation statistique 42, un histogramme des valeurs de longueur radiale estimées LR est calculé, et la dispersion des valeurs de l'histogramme par rapport à la valeur médiane est évaluée. 35 Des paramètres de dispersion sont calculés. 3050279 11 Dans un mode de réalisation, la dispersion est caractérisée par le paramètre E1 défini à partir de l'écart interquartiles, dont le calcul est défini de la manière suivante : El = Q3 Où Q1 représente le premier quartile, c'est-à-dire la première valeur (dans l'ordre 5 croissant) d'une série supérieure ou égale à 25% des données, et Q3 représente le troisième quartile, c'est-à-dire la première valeur (dans l'ordre croissant) d'une série supérieure ou égale à 75% des données. En d'autres termes, l'écart interquartile indique la largeur de l'intervalle dans lequel se trouvent les 50% des valeurs centrales. Selon une variante, une autre méthode consiste à évaluer l'écart entre la moyenne 10 et la médiane des valeurs de l'histogramme. S'il y a une différence significative entre la moyenne et la médiane, c'est qu'il y a une forte dispersion et potentiellement plusieurs pôles ne permettant pas de statuer sur la longueur radiale estimée finale. Bien entendu, on peut utiliser toute variante à la portée de l'homme du métier pour l'estimation de la dispersion de l'histogramme. 15 Ainsi, la cohérence des valeurs de longueur radiale estimées est évaluée, ce qui permet de valider ou d'invalider les estimations, et d'associer un score de confiance à la longueur radiale estimée finale LR de la plateforme. Dans un mode de réalisation, plusieurs scores de confiance sont prédéfinis, un premier score correspondant à une confiance élevée, un deuxième score correspondant à 20 une confiance moyenne et un troisième score correspondant à une confiance faible. Lorsque le paramètre El d'écart interquartiles est utilisé pour estimer la dispersion, on définit deux valeurs d'écart interquartiles de seuil, E1,1 et E1,2. Par exemple, les valeurs de seuil E1,1 et E1,2 sont déterminées préalablement, soit par calcul théorique, soit à partir d'observations de test. 25 Lorsque la valeur de E1 est inférieure à E1,1, le premier score est associé à la longueur radiale estimée finale LR. Lorsque la valeur de E1 est supérieure à E1,1 mais inférieure à une deuxième valeur de seuil E1,2, le deuxième score est associé à la longueur radiale estimée finale LR. Lorsque la valeur de E1 est supérieure à E1,2 le troisième score est associé à la 30 longueur radiale estimée finale LR. Plusieurs cas de figure sont illustrés aux figures 5 à 7. Selon un premier cas de figure, l'histogramme comporte un seul pôle centré sur la valeur médiane, comme illustré schématiquement à la figure 5. Dans ce premier cas, les estimations sont validées et l'écart interquartiles E1 est 35 inférieur au premier seuil E1,1. On peut déduire que la médiane des valeurs de longueur 3050279 12 radiale estimées est une estimation représentative de la longueur radiale LR de la plateforme. Le premier score de confiance est associé à la longueur radiale estimée finale. Selon un deuxième cas de figure, l'histogramme comporte deux pôles, dont un 5 pôle majoritaire centré sur la valeur médiane, et un pôle minoritaire, comme illustré schématiquement à la figure 6. Dans ce deuxième cas, l'écart interquartiles El est supérieur au premier seuil E1,1 mais inférieur au deuxième seuil E1,2, la longueur radiale estimée finale est validée, car elle appartient au pôle majoritaire, mais avec un score de confiance moins bon que dans 10 le premier cas. On lui associe le deuxième score de confiance, inférieur au premier score de confiance. Selon un troisième cas de figure, l'histogramme comporte plusieurs pôles sans pôle majoritaire, et donc présente une dispersion importante, comme illustré schématiquement à la figure 7. 15 Dans ce troisième cas, l'écart interquartiles E1 est supérieur au deuxième seuil E1,2, la valeur médiane de la distribution des valeurs de longueur radiale estimées n'est pas une estimation représentative de la longueur radiale de la plateforme, les estimations sont bruitées. Un plus faible score de confiance que dans les autres cas sera associé à la longueur radiale estimée finale. On lui associe le troisième score de confiance, inférieur 20 au deuxième score de confiance. Ensuite, une dimension de la plateforme, parmi sa longueur ou sa largeur, est calculée à l'étape 44 en fonction de la longueur radiale estimée finale, LR. La sélection de la dimension estimée parmi la longueur L et la largeur I dépend de l'angle 0 entre l'axe Z et la direction du vecteur de vitesse i/ de la plateforme 2. 25 Si la valeur de l'angle 9 est inférieure ou égale à 45°, on estime la longueur L de la plateforme par : L= e cos(0) Si la valeur de l'angle 0 est strictement supérieure à 45°, on estime la largeur I de la plateforme par : 30 / = LR x sin(0) Pour la dimension de la plateforme, L ou I, calculée ainsi automatiquement, on conserve le score de confiance préalablement évalué. Ensuite, la dimension de la plateforme et son score de confiance sont fournis à une étape 46 de mise en oeuvre de ces valeurs pour la surveillance. 3050279 13 Dans un mode de réalisation, lorsque la plateforme considérée est un navire, l'étape 46 comprend l'obtention, lors d'une étape 48, de la valeur correspondante prévue de longueur, notée Linit ou de largeur, notée à partir des données de type AIS correspondant à un identifiant transmis par le navire 2. 5 Selon un mode de réalisation, l'étape 48 n'est réalisée que si le score de confiance attribué à la dimension de la plateforme calculée est supérieur ou égal à un seuil, par exemple supérieur ou égale au deuxième score décrit ci-dessus. Une corrélation entre la valeur de longueur L (respectivement de largeur I) calculée automatiquement et la valeur de longueur Linit (respectivement la valeur de 10 largeur linit) prévue est effectuée à l'étape 50, par exemple en calculant la valeur absolue de la différence entre ces deux valeurs, et en comparant cette différence à un seuil d'alerte prédéterminé. Si la différence est supérieure au seuil d'alerte prédéterminé, une alerte est remontée à l'étape 52. 15 Par exemple, une alerte sonore ou visuelle est remontée via une interface homme- machine du dispositif programmable mettant en oeuvre l'invention. Selon une variante, l'étape 46 comprend l'affichage sur un écran de visualisation destiné à un opérateur de la dimension estimée de la plateforme et du score de confiance associé. 20 Avantageusement, le procédé de l'invention permet de calculer automatiquement à distance une dimension d'une plateforme en mouvement parmi sa longueur et sa largeur, et ainsi de valider ou non l'auto-déclaration fournie par cette plateforme. Avantageusement, la dimension calculée automatiquement a un score de confiance associé, permettant d'améliorer la sécurité de la surveillance effectuée. 25 REVENDICATIONS1.- Procédé de calcul automatique d'une dimension d'une plateforme mobile, ladite dimension étant une dimension de longueur ou une dimension de largeur, comportant une étape d'émission, par un système radar, d'ondes radioélectriques en direction de ladite plateforme selon un axe d'émission prédéterminé et d'acquisition d'au moins un signal numérique de profile de puissance représentatif d'une puissance de signal réfléchi reçu en fonction d'une distance radiale selon ledit axe d'émission par rapport à un point de référence dudit système radar, caractérisé en ce qu'il comporte des étapes de : - application (34) d'un opérateur de filtrage sur le signal numérique de profil de puissance acquis (Si) permettant d'obtenir un signal numérique filtré, - détermination (36) par calcul d'une première distance radiale (dr1) correspondant à un premier pic de variation du signal numérique filtré et d'une deuxième distance radiale (dr2) correspondant à un deuxième pic de variation du signal numérique filtré, - calcul (38-44) d'une dimension radiale de la plateforme en fonction desdites première et deuxième distance radiale. 2.- Procédé selon la revendication 1, caractérisé en qu'il comprend l'acquisition d'une pluralité de signaux numériques de profil de puissance, et en ce que les étapes d'application d'un opérateur de filtrage et de détermination d'une première et d'une deuxième distance radiale sont appliquées pour chaque signal numérique de profil de puissance acquis. 3.- Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que le calcul d'une dimension radiale de la plateforme comprend le calcul (38) d'une valeur de longueur radiale estimée pour chaque signal numérique de profil de puissance acquis (Si) en effectuant une différence entre la deuxième distance radiale et la première distance radiale, et la mémorisation de la valeur de longueur radiale estimée à partir de chaque signal numérique de profil de puissance acquis. 4.- Procédé selon la revendication 3, caractérisé en ce qu'il comprend un calcul (40) d'une longueur radiale estimée finale à partir des valeurs de longueur radiale estimées mémorisées. 3050279 15 5.- Procédé selon la revendication 4, caractérisé en ce que la longueur radiale estimée finale est égale à la médiane des valeurs de longueur radiale estimées mémorisées. 5 6.- Procédé selon l'une des revendications 4 ou 5, caractérisé en ce qu'il comprend une étape de validation statistique (42) des valeurs de longueur radiale estimées mémorisées, et une étape d'attribution d'un score de confiance associé à la longueur radiale estimée finale. 10 7.- Procédé selon la revendication 6, caractérisé en ce que l'étape de validation statistique (42) comprend un calcul d'histogramme des valeurs de longueur radiale estimées, et un calcul d'au moins un paramètre de dispersion associé audit histogramme. 8.- Procédé selon l'une quelconque des revendications 1 à 7, caractérisé en ce 15 que la première distance radiale (dr1) est la distance radiale la plus petite correspondant à un pic de variation du signal numérique filtré, et la deuxième distance radiale (dr2) est la distance radiale la plus grande correspondant à un pic de variation du signal numérique filtré. 20 9.- Procédé selon l'une quelconque des revendications 1 à 8, caractérisé en ce que l'étape de détermination par calcul d'une première distance radiale correspondant à un premier pic de variation du signal numérique filtré et d'une deuxième distance radiale correspondant à un deuxième pic de variation du signal numérique filtré comprend un calcul (70), effectué par application d'une fenêtre glissante, d'une valeur 25 d'intégrale discrète partielle du signal numérique filtré compris dans ladite fenêtre glissante et une comparaison (72) de la valeur d'intégrale partielle à un seuil de variation prédéterminé. 10.- Procédé selon la revendication 9, caractérisé en ce que le seuil de variation 30 prédéterminé est calculé en fonction de l'intégrale discrète totale du signal numérique filtré. 11.- Procédé selon l'une quelconque des revendications 1 à 10, caractérisé en ce qu'il comporte en outre une étape de calcul (44) d'une dimension de longueur ou de largeur en fonction de la dimension radiale calculée et une étape d'estimation de corrélation (50) entre la dimension de longueur ou de largeur calculée et une dimension 3050279 16 de longueur ou de largeur prévue, et une levée d'alerte (52) en fonction du résultat de l'étape d'estimation de corrélation. 12.- Programme d'ordinateur comportant des instructions logicielles qui, 5 lorsqu'elles sont mises en oeuvre par un dispositif programmable, mettent en oeuvre un procédé de calcul automatique d'une dimension d'une plateforme mobile, ladite dimension étant une dimension de longueur ou une dimension de largeur, conforme aux revendications 1 à 11. 10 13.- Support d'enregistrement d'informations, caractérisé en ce qu'il comporte des instructions logicielles pour l'exécution d'un procédé conforme aux revendications 1 à 11, lorsque ces instructions sont exécutées par un dispositif programmable. 14.- Système de calcul automatique d'une dimension d'une plateforme mobile, 15 ladite dimension étant une dimension de longueur ou une dimension de largeur, comportant une antenne d'émission d'ondes radioélectriques en direction de ladite plateforme selon un axe d'émission prédéterminé et un module d'acquisition d'au moins un signal numérique de profile de puissance représentatif d'une puissance de signal réfléchi reçu en fonction d'une distance radiale selon ledit axe d'émission par rapport à un 20 point de référence de l'antenne, caractérisé en ce qu'il comporte une unité (14) de calcul et de traitement comportant un processeur et une ou plusieurs unités de mémoire (16) aptes à stocker des données et des instructions de code de programme, mettant en oeuvre : - un module d'application d'un opérateur de filtrage sur le signal numérique de 25 profil de puissance acquis (S1) permettant d'obtenir un signal numérique filtré, - un module de détermination par calcul d'une première distance radiale (dr1) correspondant à un premier pic de variation du signal numérique filtré et d'une deuxième distance radiale (dr2) correspondant à un deuxième pic de variation du signal numérique filtré, 30 -un module de calcul d'une dimension radiale de la plateforme en fonction desdites première et deuxième distance radiale.