Brevet FR3082666A1 — publié le 20191220

Procédé de mesure de la hauteur de vagues à l'aide d'un radar aéroporté. [1] L'invention se rapporte au domaine général de la surveillance maritime et des aéronefs chargés de telles opérations. [2] Elle se rapporte notamment aux hydravions de surveillance maritime 5 et aux problématiques liées aux conditions d'amerrissage et de mesure de hauteurs des vagues. CONTE.»,TE [3] Les radars de surveillances maritimes sont installés sur divers porteurs y compris sur certaines classes d'hydravions. 10 [04] Pour amerrir en toute sécurité, ces hydravions doivent contrôler la hauteur des vagues afin de vérifier qu'elle est inférieure à un seuil de sécurité. [5] Aujourd'hui, cette hauteur est en général acquise via des infrastructures de mesures positionnées au niveau de la mer ou embarquées 15 sur des navires et transmise à l'hydravion considéré. [6] Afin de gagner en autonomie, les opérationnels souhaitent disposer, à bord de l'aéronef, de moyens d'estimer de manière autonome la hauteur des vagues, notamment dans la zone d'amerrissage. Une telle autonomie permettrait plus de liberté dans la définition des missions et une plus grande 20 flexibilité du plan de vol. [7] Aujourd'hui les installations radar permettant de mesurer la hauteur des vagues sont essentiellement des installations fixes comportant un radar de surface dédié de manière plus ou moins importante à cette tâche de mesure. Cependant, ces installations opèrent avec des contraintes très 25 différentes de celles pouvant être imposées dans le contexte de l'invention. En particulier, le temps d'acquisition du signal permettant de calculer cette hauteur peut être de plusieurs minutes, durée incompatible du fonctionnement d'un radar embarqué. [8] Par suite, réaliser une mesure de la hauteur des vagues en utilisant 2 le radar embarqué à bord nécessite de définir et concevoir une chaîne de traitement intégrée au radar aéroporté présentant les caractéristiques suivantes: - Le temps d'analyse du signal doit être relativement court afin de ne 5 pas contraindre la trajectoire de l'hydravion et de ne pas mobiliser le radar de manière prohibitive. Le temps d'analyse doit à ce titre être du même ordre que celui requis pour d'autres traitements tels que les traitements d'imagerie de cibles maritimes type ISAR ("Inverse Synthetic Aperture Radar" selon la dénomination anglo-saxonne), c'est-à-dire inférieur, typiquement, à la dizaine 10 de secondes. -L'altitude de l'aéronef à laquelle la hauteur des vagues est estimée doit être suffisamment haute, supérieure à 100 pieds typiquement, pour permettre une annulation de la manoeuvre d'amerrissage en toute sécurité. - Le traitement mis en oeuvre doit fournir une mesure de la hauteur 15 des vagues avec la précision attendue (erreur maximale sur l'estimation de la hauteur de vague de l'ordre de 50 cm) et de façon indépendante dans tous les contextes d'amerrissage. Sur les systèmes statiques, la précision de l'estimation de la hauteur de vague à un endroit donné repose en effet sur une calibration de la mesure à l'aide de capteurs externes (bouée) ce qui 20 n'est pas envisageable dans le contexte d'un système embarqué que l'on souhaite autonome. [9] A l'heure actuelle il existe deux types connus de méthodes de traitements du signal radar permettant d'estimer la hauteur des vagues par analyse du signal rétrodiffusé par la surface de la mer avec un radar en 25 bande X: - Les méthodes fondées sur l'inversion de la Fonction de Transfert de Modulation, FTM, (ou " Modulation Transfer Function" (MTF) selon la dénomination anglo-saxonne); -Les méthodes fondées sur l'analyse des fréquences Doppler du 30 signal rétrodiffusé. [10] Le principe général des méthodes fondées sur l'inversion de la FTM du signal rétrodiffusé est de déterminer la hauteur des vagues à partir de la puissance des échantillons de fouillis ("clutter") de mer. [11] Pour ce faire, sur un temps de mesure donné, un radar capte un 3 signal de fouillis ("clutter") de mer à partir duquel il forme une image correspondant à la répartition (distribution) de la puissance reçue dans un repère distance/temps (ou récurrences). Sur une telle image, présentée à titre d'illustration sur la figure 1, l'intensité de chaque pixel représente la puissance de signal reçu pour l'échantillon de signal correspondant une case distance donnée pour une récurrence donnée. [12] En présence de fouillis de mer, l'image ainsi construite est structurée selon des « zébrures » verticales, orientées selon l'axe des temps (axe récurrences) traduisant la présence de vagues à la surface de l'eau. [13] Conformément à cette méthode, on déduit la hauteur des vagues à partir cette structure, la hauteur des vagues étant accessible à travers la puissance du fouillis de mer. [14] Dans le contexte considéré, l'emploi d'une telle méthode présente certains inconvénients: [15] En premier lieu il s'agit d'une méthode, coûteuse d'un point de vue calculatoire, qui nécessite une étape intermédiaire de traitement. En effet, la « prise de vue par le radar » entraine une distorsion de l'image créée, distorsion qui peut être modélisée par une fonction appelée « Modulation Transfer Function » (MTF). [16] Par suite la détermination précise de la hauteur des vagues nécessite l'estimation de cette MTF puis son inversion. La MTF inversée est ensuite appliquée à la transformée de Fourier 2D de l'image distance/temps formée à partir du signal radar reçu. Puis la transformée de Fourier inverse est appliqué à l'image ainsi obtenu par application de la MTF inversée de façon à obtenir une image distance/temps corrigée des distorsions à partir de laquelle les hauteurs des vagues sont estimées. [17] En second lieu, ce type de méthode nécessite une calibration. En effet le calcul de MTF dépend du contexte (lieu, altitude radar, etc.), de sorte que pour calibrer cette MTF, il faut des disposer de moyens sol (i.e. 30 sur la mer) de ce type de méthode peu (voire pas du tout) compatible d'une utilisation mobile et autonome. [18] En troisième lieu, ce type de méthode permet de déterminer la hauteur des vagues seulement avec une précision limitée car il est basé sur 4 l'hypothèse sous-jacente que la hauteur des vagues est proportionnelle à la puissance du fouillis de mer. Par ailleurs, la modulation de la puissance du signal due au "tilt" des vagues, est très sensible à l'altitude du capteur, et n'est donc pas particulièrement adaptée à la réalisation de mesure à partir 5 d'un aéronef [19] Enfin, en dernier lieu ce type de méthode, qui permet une estimation de la hauteur de vague pour chaque vague, apparait sur-spécifié vis-à-vis du besoin consistant à déterminer la hauteur moyenne des vagues dans une zone donnée. 10 [20] Les méthodes fondées sur l'analyse des fréquences Doppler du signal rétrodiffusé sont quant à elles basées sur le fait que la hauteur significative des vagues peut être, de manière connue, corrélée avec l'écart-type des fréquences Doppler mesurée en un point donné d'une zone considérée. 15 [21] Selon ce type de méthode, les fréquences Doppler du signal reçu en ce point sont calculées sur une durée donnée, de l'ordre de 15 minutes par exemple. Par ailleurs, de manière connue, la hauteur significative des vagues hs est reliée à l'écart-type a, des vitesses (i.e. des fréquences doppler). 20 [22] On rappelle ici que la hauteur significative des vagues h, représente la valeur moyenne de la hauteur du tiers des vague de plus grande hauteur sur la zone considérée et pendant la duré de l'observation effectuée. [23] Ce type de méthode offre une approche plus précise et plus simple à mettre en oeuvre que les méthodes classiques, basées sur l'estimation de la 25 MTF. Il présente cependant l'inconvénient de nécessiter un temps d'analyse incompatible d'une utilisation dans un contexte opérationnel tel qu'une manoeuvre d'amerrissage; la durée d'analyse du fouillis de mer, qui nécessite l'asservissement de l'antenne dans une direction donnée, étant trop importante. Par suite aucun des deux types de méthodes exposés précédemment ne donnant réellement satisfaction dans le contexte opérationnel particulier considéré ici, il existe à ce jour un besoin de trouver un moyen de traitements du signal radar permettant d'estimer de manière simple et rapide la hauteur des vagues. 5 PRES EN TA [24] Un but de l'invention est notamment de proposer un tel moyen, autrement dit, de proposer une méthode permettant d'estimer la hauteur des vagues à partir du signal capté par un radar aéroporté, ladite méthode 5 prenant en compte les contraintes opérationnelles suivantes : - le temps de réponse doit être court, typiquement inférieur à 10 secondes; - le traitement mis en oeuvre doit s'appliquer dans différents contextes maritimes sans nécessiter d'étape de calibration; 10 - la précision de l'estimation de la hauteur de vague doit être fiable et précise, et surtout éviter les risques de sous-estimation de la hauteur de vague. A cet effet l'invention a pour objet un procédé pour déterminer la hauteur des vagues au moyen d'un radar porté par un aéronef, ledit procédé 15 mettant en oeuvre les étapes suivantes : - une première étape de pointage de l'antenne du radar; - une deuxième étape de détermination du plan d'acquisition du fouillis en fonction de l'altitude de l'aéronef; - une troisième étape de détermination, pour chaque zone de fouillis 20 définie par le plan d'acquisition, de deux paramètres Doppler PARAI et PARA2 caractérisant ladite zone dans son ensemble; - une quatrième étape de calcul, pour chacun des paramètres PARAI et PARA2 calculés lors de la troisième étape, de la valeur moyenne dudit paramètre sur l'ensemble des zones considérées; 25 - une cinquième étape d'estimation proprement dite de la hauteur de vagues à partir des moyennes des paramètres PARAI et PARA2 calculées lors de la quatrième étape. La hauteur de vague ainsi estimée est destinée à être transmise à l'aéronef pour déterminer les conditions d'amerrissage de ce dernier 30 Selon la forme de mise en oeuvre considérée le procédé selon l'invention peut présenter des caractéristiques complémentaires ou alternatives listées ci-après. Ces caractéristiques peuvent être 6 considérées chacune séparément ou bien en combinaison avec une ou plusieurs autres caractéristiques. Selon une première caractéristique, la troisième étape comporte elle-même trois sous-étapes: 5 - une première sous-étape d'acquisition du signal correspondant à la zone considérée; - une deuxième sous-étape de prétraitement Doppler qui détermine la matrice des fréquences instantanées du signal reçu, représentée dans l'espace distance/temps (ou distance/récurrences); 10 - une troisième sous-étape de calcul, pour chaque zone, de la valeur de deux paramètres PARAI et PARA2 à partir de la matrice des fréquences instantanées délivrée par la deuxième sous-étape. Selon l'invention, le paramètre PARAI, correspondant à l'écart-type, sur l'ensemble des cases distance, des fréquences Doppler moyennes 15 calculées pour chaque case distance, sur l'ensemble des récurrences formant la zone considérée. le paramètre PARA2, quant à lui, correspond au rapport du premier paramètre PARAI sur l'écart-type de la fréquence Doppler moyenne du signal pour chaque récurrence. Selon une autre caractéristique, la deuxième sous-étape de la 20 troisième étape de prétraitement Doppler, met en oeuvre: - une première opération de calcul de la matrice des fréquences instantanées du signal représenté dans l'espace distance/temps (ou distance/récurrences); chaque fréquence Doppler instantanée fa de ladite matrice étant déterminée, 25 pour chaque case distance et chaque récurrence définissant la zone considérée, à partir de la relation: dcpxPRF fd 241r où PRF représente la période de récurrence du radar et où chp représente la différence de phase, pour une case distance donnée, entre les échantillons 30 de fouillis reçus par le radar pour deux récurrences consécutives; - une seconde opération de calcul d'une matrice des fréquences instantanées débruitée, ladite matrice étant obtenue en appliquant un filtre 7 passe bas à la matrice des fréquences instantanées calculée par la première opération. Selon une autre caractéristique, la cinquième étape réalise: - le calcul, pour chacun des paramètres PARAI et PARA2, d'une 5 première estimation de hauteur de vague, hestPARAl et hestpARA2 respectivement, le calcul étant réalisé pour chaque paramètre, en utilisant un estimateur donné qui exprime la hauteur de vague estimée en fonction de la valeur du paramètre considéré; - l'élaboration d'une estimation consolidée, Hest, de la hauteur de 10 vague, cette estimation étant fonction des premières estimations de hauteurs de vague hestPARAl et hestPARA2 précédemment calculées. Selon une autre caractéristique, les premières estimations de hauteurs de vague hestPARAl et hestPARA2, sont déterminées par la relation: hestp = A x p + B, 15 où p est la valeur du paramètre considéré, PARAI ou PARA2. Selon une autre caractéristique, l'estimation consolidée de la hauteur de vague est définie comme étant égale à la plus grande des premières estimations de hauteurs de vague hestPARAl et hestPARA2: Hest = max(hestPARAl, heStpARA2) 20 Selon une autre caractéristique, l'estimation consolidée de la hauteur de vague est définie comme étant égale à la moyenne des premières estimations de hauteur de vague hestPARAl et hestPARA2 Selon une autre caractéristique, l'estimation consolidée de la hauteur de vague est définie comme étant égale à une somme pondérée des 25 premières estimations de hauteur de vague hestPARAl et hestPARA2. [25] Les caractéristiques et avantages de l'invention seront mieux appréciés grâce à la description qui suit, description qui s'appuie sur les figures annexées qui présentent: 8 la figure 1, une image simplifiée, monochrome représentant graphiquement un exemple de matrice des puissances de signal reçu, mesurées pour une zone de l'espace (cases distance-récurrences) de taille donnée. 5 la figure 2, un organigramme de principe présentant les différentes étapes du procédé selon l'invention; la figure 3, un organigramme de principe détaillant la troisième étape du procédé selon l'invention; la figure 4, une image simplifiée, monochrome représentant 10 graphiquement la matrice des fréquences doppler instantanées corrigée, obtenue, selon l'invention, par prétraitement de la matrice illustrée par la figure 3; la figure 5, un organigramme de principe détaillant la cinquième étape du procédé selon l'invention. 15 Il est à noter que, sur les figures annexées, un même élément fonctionnel ou structurel porte de préférence un même symbole repère. DESCRIPTION .°h,# .L [26] Le procédé selon l'invention détermine la hauteur de la mer en se basant sur la mesure de la fréquence Doppler du fouillis de mer ou "clutter". 20 Ce mode de détermination exploite le fait que les vagues génèrent une modulation de la vitesse mesurée car elles introduisent, par gravité, des variations de vitesse mesurée à la surface de la mer. La fréquence Doppler, qui est une valeur intrinsèque du fouillis de mer, dont la valeur n'est impactée ni par le radar ni par sa position, est en 25 effet une grandeur dont la variabilité de la mesure est très faible. Sa mesure ne présente pas les inconvénients cités précédemment, de sorte qu'elle semble mieux adaptée aux conditions d'une mesure effectuée depuis un aéronef. On rappelle à ce propos que la vitesse radiale du fouillis, g., est liée à 30 la fréquence Doppler par la relation suivante: fa = Vr - (2/a) [02] où X. représente la longueur d'onde de l'onde radar émise. [27] Comme l'illustre la figure 2, le procédé selon l'invention met en 9 oeuvre les étapes de traitements suivantes : - une première étape 21 de pointage de l'antenne; - une deuxième étape 22 de détermination du plan d'acquisition du fouillis en fonction de l'altitude de l'aéronef; 5 - une troisième étape 23 de détermination, pour chaque zone de fouillis définie par le plan d'acquisition, de deux paramètres Doppler PARAI et PARA2 caractérisant ladite zone. - une quatrième étape 24 de calcul, pour chaque paramètre, PARAI ou PARA2, de la valeur moyenne sur l'ensemble des zones du paramètre 10 considéré, calculé à l'étape 23; - une cinquième étape 25 d'estimation proprement dite de la hauteur de vagues à partir des moyennes des paramètres PARAI et PARA2 calculées lors de l'étape 24. [28] La première étape 21 de pointage de l'antenne consiste à garantir 15 que la direction visée par le radar correspond à la direction d'origine de la houle (avec une tolérance de l'ordre de 20°). Cette étape, qui peut être implémentée par tout moyen connu de l'homme du métier, n'est par ailleurs pas détaillée dans le contexte de la présente invention. [29] La deuxième étape 22 réalise la détermination du plan d'acquisition 20 du fouillis. Ce plan définit, en fonction de l'altitude de l'aéronef, le nombre de zones à acquérir, la durée d'acquisition pour chaque zone ainsi que la position et la longueur des zones. [30] Selon l'invention, les dimensions des zones sont définies de façon à permettre une analyse des corrélations spatio-temporelles du fouillis. 25 Plus précisément une zone est dimensionnée de façon à contenir plusieurs pseudo-périodes de houle. Chaque zone peut ainsi par exemple présenter une longueur supérieure ou égale à 1000 mètres. [31] Par ailleurs, la durée d'acquisition pour une zone donnée est liée à l'hypothèse de stationnarité du fouillis sur la durée d'acquisition et doit donc 30 être typiquement inférieure à 0,5 secondes. [32] La troisième étape 23, illustrée par la figure 3, est une étape itérative 10 dont la fonction consiste à déterminer, pour chacune des zones définies par le plan d'acquisition, la valeur de deux paramètres PARAI et PARA2, fonctions de la fréquence doppler du fouillis mesuré sur la zone considérée. Cette troisième étape comporte elle-même trois sous-étapes: 5 - une première sous-étape 231 d'acquisition du signal correspondant à la zone considérée; - une deuxième sous-étape 232 de prétraitement Doppler qui met en oeuvre: - une première opération 2321 de calcul de la matrice des 10 fréquences instantanées du signal représenté dans l'espace distance/temps (ou distance/récurrences); - une seconde opération 2322 de calcul de la matrice des fréquences instantanées débruitée; - une troisième sous-étape 233 de calcul, pour chaque zone, de la valeur 15 de deux paramètres PARAI et PARA2 à partir de la matrice des fréquences instantanées débruitée; [33] L'acquisition, au cours de la troisième sous-étape 231, du signal réfléchi par chaque zone, est réalisée de manière classique et n'est pas détaillée ici. On précise cependant, que cette acquisition est réalisée en 20 effectuant la compensation du déplacement du radar de façon à ce qu'un même élément de fouillis soit localisé à la même distance pendant toute la durée d'acquisition. [34] Selon l'invention, le signal acquis au cours de la première sous-étape 231 de l'étape 23, qui correspond à une zone de fouillis donnée 25 caractérisée par un nombre de cases distance (Ncd) qui correspond au nombre de cases distance sur l'intervalle de distance étudié, et un nombre de récurrence (Nrec) qui correspond au nombre de récurrences émises durant la durée analysée. Pour chacun des couples case distance/récurrence, le radar détermine un échantillons de fouillis sous 30 forme d'une amplitude complexe. [35] A partir de ces échantillons, on détermine la fréquence Doppler instantanée (i.e. la fréquence doppler mesurée pour chaque case distance et 11 une récurrence donnée) en appliquant la formule suivante : clexPRF [03] fd = 2*re où le facteur PRF ("Pulse Repetition Frequency" selon la dénomination anglo-saxonne) représente la période de répétition, ou de récurrence, du 5 radar une période égale à 1950 Hz par exemple et correspond au nombre de récurrence émis en une seconde; et où clyo est la différence de phase, pour une case distance donnée, entre les échantillons de fouillis (i.e. échantillons du signal reçu par le radar) mesurés pour deux récurrences consécutives. 10 [36] La constitution de la matrice des fréquences instantanées dans l'espace distance/temps (ou distance/récurrences) formée par l'ensemble des fréquences doppler instantanées calculées pour toute les case distance de la zone considérée et pour toutes les récurrences de la période d'analyse est réalisée par la première opération 2321 de la sous-étape 232. 15 [37] La matrice ainsi obtenue est cependant difficilement exploitable en l'état. En effet, la variation de la fréquence Doppler de case à case et de récurrence à récurrence apparait fortement bruitée de sorte qu'il est difficile de déterminer une hauteur de vague à partir des variations de fréquence doppler. 20 [38] Pour pallier cet inconvénient, la sous-étape 232 du procédé selon l'invention met en oeuvre une seconde opération 2322, dont l'objet est de renforcer le contraste de l'image représentée par la matrice obtenue après mise en oeuvre de l'opération 2321. [39] Cette seconde opération réalise la construction d'une matrice des 25 fréquences instantanées débruitée, c'est-à-dire une matrice expurgée du bruit provoqué, notamment, par la forte variabilité du fouillis de mer. [40] Cette matrice expurgée est obtenue en appliquant une fonction de filtrage passe-bas consistant, par exemple, à réaliser la convolution de la matrice des fréquences instantanées obtenue après exécution de l'opération 30 1321 avec une matrice multiple de la matrice unité, dont les coefficients sont tous égaux, par exemple une matrice 5 par 11 dont tous les coefficients 12 valent 1/55. On obtient ainsi, comme l'illustre la représentation de la figure 4, une matrice des fréquences instantanées corrigée qui laisse apparaitre des variations de fréquence doppler plus franches laissant apparaitre des cases 5 distances pour lesquelles la fréquence doppler est plus élevée (lignes en noir sur la figure 4) et d'autres ou elle est plus faible . [41] Selon l'invention également, les deux paramètres PARAI et PARA2, calculés au cours de la troisième sous-étape 233, permettent avantageusement d'exploiter au mieux les corrélations entre la modulation to de la fréquence Doppler instantanée du signal de fouillis réfléchi par la mer et la hauteur de vagues. [42] Le premier paramètre, PARAI, correspond à l'écart-type des fréquences Doppler moyennes calculées pour chaque case distance. Ce paramètre est déterminé en deux temps: 15 - un premier temps durant lequel on effectue le calcul, pour chaque case distance, de la fréquence Doppler moyenne sur l'ensemble des récurrences définissant la zone considérée, fdicd, définie par la relation suivante: xgreecc.i f d icd,ir e c [04] T iCd = 1 Nr ec 20 où f d icd,irec représente la fréquence Doppler filtrée (filtrage passe bas) associée à la case distance icd et à la récurrence irec. - un second temps durant lequel on effectue le calcul du paramètre PARAI proprement dit, ce premier paramètre étant défini par la relation suivante: 25 \ 2 [05] PARAI = x E,gera (Rica - fc/cd)luta où fdcd représente la valeur moyenne des valeurs fdicd. [43] Le second paramètre, PARA2, correspond au rapport du premier paramètre PARAI sur l'écart-type de la fréquence Doppler moyenne pour chaque récurrence. Ce second paramètre est également déterminé en deux 30 temps: 13 - un premier temps durant lequel on effectue le calcul, pour chaque de la fréquence Doppler moyenne pour chaque récurrence fair, définie par la relation suivante: ecccci=Ncd [06] fdirec :cd -1 f dicd,irec 5 - un second temps durant lequel on effectue le calcul du paramètre PARA2 proprement dit, ce premier paramètre étant défini par la relation suivante: PARA2 = PARAI [07] Nrec xe,;:ge'(ratrec-i-arec)2 où fd', représente la valeur moyenne des valeurs Tairec- 10 [44] La cinquième étape 25, illustrée par la figure 5, réalise l'estimation proprement dite de la hauteur de vagues à partir des paramètres PARAI et PARA2 calculés au cours de l'étape 23 du procédé selon l'invention. La hauteur de vague 26 ainsi estimée est prise en compte par l'aéronef lors de la phase d'amerrissage, notamment pour déterminer si les conditions limites 15 de hauteur de vague permettant l'amerrissage ne sont pas dépassées. [45] Cette estimation est précédée de l'étape 24 du calcul de la valeur moyenne de chacun des paramètres PARAI et PARA2 sur l'ensemble des zones considérées, la prise en compte d'une valeur moyenne permet avantageusement d'obtenir, pour chacun des paramètres, une valeur unique 20 plus fiable. [46] Pour chacun des paramètres PARAI et PARA2, l'étape 25 du procédé selon l'invention met en oeuvre des estimateurs, 253 et 254 reliant la valeur du paramètre considéré, PARAI ou PARA2 respectivement, à la hauteur de vague. 25 [47] On effectue ainsi le calcul, pour chacun des paramètres PARAI et PARA2, d'une première estimation de la hauteur de vague, hestpARAi et hestPARA2 respectivement, le calcul étant réalisé pour chaque paramètre, en utilisant un estimateur donné qui exprime la hauteur de vague estimée en fonction de la valeur du paramètre considéré 30 [48] Selon l'invention, ces estimateurs 253 et 254 sont déterminés préalablement lors d'une phase d'apprentissage 14 Ils peuvent suivre différents modèles (linéaire, polynomial...). Dans le cas d'un estimateur linéaire, la hauteur de vague estimée, hestp, est par exemple déterminée, à partir de la valeur p du paramètre considéré (i.e. PARAI ou PARA2), par la relation suivante: 5 hestp= Axp + B, [08] où A et B représente les coefficients du modèle, préalablement déterminés. [49] Par suite, la hauteur estimée, Hest, est déterminée à partir des premières estimations hestp calculées, au cours des opérations 251 et 252, à partir des estimateurs 253 et 254 choisis. 10 [50] Selon l'invention la hauteur estimée 26 peut être obtenue de différentes façons. [51] Ainsi, par exemple: - dans une forme particulière du procédé selon l'invention, la hauteur estimée 26 est obtenue en comparant les premières estimations des 15 hauteurs de vague hest et hest délivrées par chacun des -PARAI -PARA2 estimateurs et en retenant la valeur la plus forte. On a alors: Hest =max(hestpARAi,hestpApA2); [09] - dans une forme de mise en oeuvre alternative, la hauteur de vague estimée 26 est obtenue en calculant la moyenne des premières estimations 20 des hauteurs de vague hest -PARAI et hestpARA2 délivrées par chacun des estimateurs; - dans une autre forme de mise en oeuvre alternative, REVENDICATIONS1. Procédé pour déterminer la hauteur des vagues à l'aide d'un radar porté par un aéronef, caractérisé en ce qu'il met en oeuvre les étapes suivantes : - une première étape (21) de pointage de l'antenne du radar; - une deuxième étape (22) de détermination du plan d'acquisition du fouillis en fonction de l'altitude de l'aéronef; - une troisième étape (23) de détermination, pour chaque zone de fouillis définie par le plan d'acquisition, de deux paramètres Doppler PARAI et PARA2 caractérisant ladite zone dans son ensemble; - une quatrième étape (24) de calcul, pour chacun des paramètres PARAI et PARA2 calculés lors de la troisième étape (23), de la valeur moyenne dudit paramètre sur l'ensemble des zones considérées; - une cinquième étape (25) d'estimation proprement dite de la hauteur de vagues (26) à partir des moyennes des paramètres PARAI et PARA2 calculées lors de la quatrième étape (24); la hauteur de vague (26) ainsi estimée étant transmise à l'aéronef pour déterminer les conditions d'amerrissage de ce dernier 2. Procédé selon la revendication 1, caractérisé en ce que la troisième étape (23) comporte elle-même trois sous-étapes: - une première sous-étape (231) d'acquisition du signal correspondant à la zone considérée; - une deuxième sous-étape (232) de prétraitement Doppler qui détermine la matrice des fréquences instantanées du signal reçu, représentée dans l'espace distance/temps (ou distance/récurrences); - une troisième sous-étape (233) de calcul, pour chaque zone, de la valeur de deux paramètres PARAI et PARA2 à partir de la matrice des fréquences instantanées délivrée par la deuxième sous-étape (232); le paramètre PARAI, correspondant à l'écart-type, sur l'ensemble des cases distance, des fréquences Doppler moyennes calculées pour chaque 16 case distance, sur l'ensemble des récurrences formant la zone considérée; le paramètre PARA2 correspond au rapport du premier paramètre PARAI sur l'écart-type de la fréquence Doppler moyenne du signal pour chaque 5 récurrence. 3. Procédé selon la revendication 2, caractérisé en ce que la deuxième sous-étape (232) de la troisième étape (23) de prétraitement Doppler, met en oeuvre: to - une première opération (2321) de calcul de la matrice des fréquences instantanées du signal représenté dans l'espace distance/temps (ou distance/récurrences); chaque fréquence Doppler instantanée fa de ladite matrice étant déterminée, pour chaque case distance et chaque récurrence définissant 15 la zone considérée, à partir de la relation: dexPRF [03] fd 2*n où PRF représente la période de récurrence du radar et où chp représente la différence de phase, pour une case distance donnée, entre les échantillons de fouillis reçus par le radar pour deux récurrences 20 consécutives; - une seconde opération (2322) de calcul d'une matrice des fréquences instantanées débruitée, ladite matrice étant obtenue en appliquant un filtre passe bas à la matrice des fréquences instantanées calculée par la première opération (2321). 25 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que la cinquième étape (25) réalise: - le calcul, pour chacun des paramètres PARAI et PARA2, d'une première estimation de hauteur de vague, hest -PARAI et hest -PARA2 30 respectivement, le calcul étant réalisé pour chaque paramètre, en utilisant un estimateur donné qui exprime la hauteur de vague estimée en fonction de la valeur du paramètre considéré; 17 - l'élaboration d'une estimation consolidée, Hest, de la hauteur de vague, cette estimation (26) étant fonction des premières estimations de hauteurs de vague hestpARAi et hestpARA2 précédemment calculées. 5 5. Procédé selon la revendication 4, caractérisé en ce que les premières estimations de hauteurs de vague hestpARAi et hestpARA2, sont déterminées par la relation: hestp = A x p + B, où p est la valeur du paramètre considéré, PARAI ou PARA2. 10 6. Procédé selon l'une des revendications 4 ou 5, caractérisé en ce que l'estimation consolidée de la hauteur de vague (26) est définie comme étant égale à la plus grande des premières estimations de hauteurs de vague hest PARAI et hestpARA2: Hest = max(hestpARAi, hestpARAZ) 15 7. Procédé selon l'une des revendications 4 ou 5, caractérisé en ce que l'estimation consolidée de la hauteur de vague (26) est définie comme étant égale à la moyenne des premières estimations de hauteur de vague hest PARAI et hestpARA2- 20 8. Procédé selon l'une des revendications 4 ou 5, caractérisé en ce que l'estimation consolidée de la hauteur de vague (26) est définie comme étant égale à une somme pondérée des premières estimations de hauteur de vague hestpARAi et hestpARA2. 25