Brevet FR2938074A1 — publié le 20100507

Procédé de traitement de profils verticaux mesurés de la puissance des échos renvoyés suite à une émission de signaux radar. L'invention appartient au domaine du traitement des signaux radar, en particulier pour les radars météo aéroportés. Un radar aéroporté comporte un émetteur qui envoie signaux, plus précisément des ondes radio. Celles-ci sont réfléchies par la cible et détectées par un récepteur, également incorporé au sein du radar. La position des cibles se trouvant dans le champ des ondes radio émises, est estimée grâce à leur temps de retour. Dans le cadre d'une application météorologique, une cible peut être par exemple un nuage. En d'autres termes, en interprétant les ondes renvoyées par les cibles, on en déduit leur position. Les échos renvoyés par le sol dit � échos de sol � peuvent être particulièrement gênants pour l'interprétation des signaux reçus. En effet, ces derniers peuvent être interprétés comme des signaux rétrodiffusés par des nuages. D'autant plus que selon le type de terrain observé, la puissance réfléchie par le sol est généralement aussi voire plus importante que celle réfléchie par les phénomènes météorologiques, ce qui se traduit par l'apparition d'une valeur de réflectivité élevée sur l'image radar. L'apparition de cette forte valeur de réflectivité sur une image radar va traduire un phénomène dangereux, et être indiqué en tant que tel au pilote, créant ainsi une fausse alarme. II est donc particulièrement important de s'affranchir des échos de sol, de façon à ne pas perturber le pilote. Il est connu de l'art antérieur (US 6.603.425, Collins) une méthode de détection des échos de sol à l'aide de deux balayages horizontaux, selon deux valeurs d'élévation distincte. Ces derniers sont soustraits l'un à l'autre, et la différence résultante est comparée à une valeur de référence de manière à déterminer la puissance due aux échos de sol. Toutefois, cette méthode est peu précise (en raison de faible nombre de valeur d'élévation considérée), et peu robuste dès que surviennent certains phénomènes météorologiques. En outre, de nombreux problèmes de précisions peuvent intervenir en raison du temps nécessaire pour effectuer chaque balayage horizontal. L'invention vise notamment à résoudre ce problème. A cet effet, selon un premier aspect, il est proposé un procédé de traitement de profils verticaux mesurés de la puissance des échos renvoyés suite à une émission de signaux radar, chaque profil vertical mesuré étant fonction de l'angle de balayage du faisceau radar et associé à un couple angle de pointage/case distance donné du faisceau radar. Selon une caractéristique générale de l'invention, le procédé comprend une élaboration d'un profil vertical synthétique de la puissance des échos renvoyés par le sol seulement, et pour chaque profil vertical mesuré : • une superposition du profil vertical synthétique au profil vertical mesuré considéré, pour différentes valeurs de l'angle de balayage, • un calcul de l'erreur existant entre le profil vertical mesuré et le profil vertical synthétique, pour chaque valeur de l'angle de balayage considérée, et • une exclusion du profil vertical mesuré, des valeurs pour lesquelles l'erreur calculée est inférieure à un seuil donné, dites valeurs à exclure, de façon à générer un profil vertical résultant, affranchi des échos de sol. En d'autres termes, on exclut de chaque profil vertical mesuré une zone correspondant à un profil vertical synthétique des échos de sol uniquement. Le procédé selon l'invention ne repose donc pas sur la comparaison de données, mais sur l'estimation d'un modèle qu'est le profil vertical synthétique, à partir de données observées. Cela à pour avantage de rendre le procédé beaucoup plus précis et robuste, contrairement aux solutions de l'art antérieur. De préférence, l'angle de balayage est l'angle d'élévation et l'angle de pointage est l'angle de gisement. En effet, l'utilisation de données obtenues par balayages verticaux s'avère plus fiable, car le temps nécessaire pour l'obtention d'un profil vertical est beaucoup plus court que le temps nécessaire pour obtenir une différence verticale de deux nappes, que l'on retrouve notamment dans le document US 6.603.425, Collins précité. Cela permet notamment d'éviter des problèmes de déformation ou d'ajustement de données dus au déplacement de l'avion. De préférence, ladite confrontation n'est effectuée que pour certains points du profil vertical mesuré, correspondant à des valeurs de l'angle de balayage inférieures à un seuil donné. Par exemple, ladite confrontation peut être effectuée pour l'ensemble des points du profil vertical mesuré, si le nombre de points du profil vertical mesuré est inférieur à un seuil donné. Le procédé peut comporter en outre un calcul de la hauteur du sol pour le profil vertical mesuré considéré, qui comprend : une autre superposition du profil vertical synthétique au profil vertical mesuré considéré, pour l'ensemble des valeurs de l'angle de balayage, • un calcul de l'erreur existant entre le profil vertical mesuré et le profil vertical synthétique, pour l'ensemble des valeurs de 25 l'angle de balayage, • une détermination de la valeur de l'angle de balayage du barycentre du profil vertical mesuré correspondant à l'erreur minimale calculée, • un calcul d'une hauteur dite de sol pour le profil vertical 30 mesuré considéré, à l'aide de ladite valeur de l'angle de balayage déterminée, les valeurs du profil vertical inférieures à la valeur de pointage correspondant à la hauteur de sol calculée étant alors également exclues. De préférence, lesdites valeurs à exclure et les valeurs de la hauteur de sol calculées pour un profil vertical mesuré donné, sont lissées 5 avant l'exclusion. Il est particulièrement avantageux que chaque profil vertical résultant soit filtré, de façon à éliminer les valeurs isolées. L'exclusion peut également être fonction d'un modèle numérique du sol. 10 Selon un autre aspect de l'invention, il est proposé un dispositif de traitement de signaux radar, couplé à une antenne radar apte à mesurer les profil verticaux de la puissance des échos renvoyés suite à une émission de signaux radar, chaque profil vertical mesuré étant fonction de l'angle de balayage du faisceau radar et associé à un couple angle de pointage/case 15 distance donné du faisceau radar. Selon une caractéristique générale de cet autre aspect, le dispositif comprend un moyen d'élaboration d'un profil vertical synthétique de la puissance des échos renvoyés par le sol seulement, un premier moyen de traitement apte pour chaque profil vertical mesuré 20 • à superposer le profil vertical synthétique au profil vertical mesuré considéré, pour différentes valeurs de l'angle de balayage, • à calculer l'erreur existant entre le profil vertical mesuré et le profil vertical synthétique, pour chaque valeur de l'angle de balayage considérée, et 25 un moyen d'exclusion apte à exclure du profil vertical mesuré, des valeurs pour lesquelles l'erreur calculée est inférieure à un seuil donné, dites valeurs à exclure, de façon à générer un profil vertical résultant, affranchi des échos de sol. Selon un mode de réalisation, le dispositif peut comprendre en 30 outre un deuxième moyen de traitement apte à calculer la hauteur du sol pour le profil vertical mesuré considéré, ledit moyen d'exclusion étant en outre apte à exclure les valeurs du profil vertical inférieures à la valeur de balayage correspondant à la hauteur de sol calculée. Selon un autre aspect de l'invention, il est proposé un récepteur radar incorporant un dispositif tel que décrit ci-avant. D'autres avantages et caractéristiques de l'invention apparaîtront à l'examen de la description détaillée d'un mode de mise en oeuvre d'un procédé selon l'invention ainsi que d'un mode de réalisation d'un dispositif selon l'invention, nullement limitatifs, et des dessins annexés sur lesquels : -la figure 1 illustre un exemple de déplacement d'un faisceau radar, et un profil vertical mesuré lors de ce déplacement en l'absence de phénomènes météorologiques - la figure 2a illustre un autre exemple de déplacement d'un faisceau radar, - la figure 2b illustre le un profil vertical mesuré lors du 15 déplacement représenté sur la figure 2a, en présence de phénomènes météorologiques, - la figure 3 illustre les différentes étapes mises en oeuvre d'un procédé selon l'invention, - les figures 4, 5, 6 et 7 illustrent plus précisément des étapes du 20 mode de mise en oeuvre représenté sur la figure 3, -la figure 8 détaille les différents paramètres utilisés par plusieurs étapes du mode de mise en oeuvre représenté sur la figure 7, et - la figure 9 représente un mode de réalisation d'un dispositif selon l'invention. 25 On se réfère à la figure 1. La référence AV désigne un aéroporteur, dans cet exemple un avion. Ce dernier est muni d'un radar émettant un faisceau radar FS. Sur chacun des trois schémas, respectivement A, B et C, l'angle d'élévation du faisceau radar, respectivement (PA, (PB et cpc varie. On rappelle que l'angle d'élévation (ou 30 angle de site) est l'angle formé entre l'horizon et le faisceau radar FS (dans un plan horizontal). L'espace visé par le faisceau radar est divisé en cases distance. L'une de ces case distance référencée CD est illustrée sur chacun des schémas A, B et C. Une case distance CD est une cellule d'analyse élémentaire en trois dimensions. La réponse de l'environnement à l'émission du signal radar est étudié pour chaque case distance CD. Le nombre de cases distance augmente ou diminue en fonction de la résolution du radar utilisé. Par conséquent, pour un angle de gisement fixé (on rappelle que l'angle de gisement est l'angle formé entre l'axe de l'avion et le faisceau, dans le plan horizontal), le profil vertical de la réponse de la case distance CD représentée sur chacun des schémas A, B, C est réalisé. Ce profil vertical est fonction de l'angle de balayage du faisceau radar FS, ici l'angle d'élévation. Plus précisément, le profil vertical représenté illustre la puissance 15 reçue par l'antenne radar de l'avion, pour une case distance CD donné, en fonction de l'angle d'élévation. En l'absence de tout phénomène météorologique, le profil vertical mesuré possède la forme du diagramme de l'antenne radar. La puissance reçue est alors principalement due aux échos renvoyés par le sol, dits échos 20 de sol. La figure 2a illustre le cas où le faisceau radar FS de l'avion AV rencontre un nuage NU. Dans ce cas, le profil vertical mesuré PM, (ici pour la case distance CD), possède une forme telle que celle représentée sur la courbe en pointillé de la figure 2b. La portion représentée en trait plein 25 représente la puissance reçue due aux échos de sol. De préférence, chaque profil vertical est formé à l'aide d'au moins deux voire trois valeurs distinctes de l'angle d'élévation. Un but de l'invention est de détecter au sein d'un profil vertical mesuré la portion de la puissance reçue due aux échos de sol, ce qui est 30 particulièrement délicat lorsque l'avion AV se trouve en présence de phénomènes météorologique. On se réfère à présent à la figure 3 qui illustre un mode de mise en oeuvre d'un procédé selon l'invention. Une fois les données reçues par l'antenne radar, et échantillonnées, le procédé peut comporter une première étape 10 (optionnelle) de pré-traitement des données que sont les profils verticaux mesurés. On note que l'échantillonnage peut être un suréchantillonnage pour pallier une éventuelle faiblesse de résolution angulaire des profils verticaux mesurés. Celle-ci peut consister à compenser les migrations en distance lorsque cela est nécessaire. En d'autres termes, on compense pour chaque mesure, le fait que l'avion se déplace au cours de la réception des données. Cette étape de pré-traitement 10 peut également comprendre un alignement des données de façon que chaque profil vertical mesuré corresponde à la même zone géographique du sol. Cette étape de pré-traitement peut ne pas être nécessaire dans le cas où une antenne radar de type à balayage électronique vertical, est utilisée. Celle-ci permet d'obtenir une coupe verticale beaucoup plus rapidement qu'une antenne radar à balayage mécanique, et limiter ainsi les effets de migration de distance. Une deuxième étape 20 comporte une pondération temporelle et éventuellement spatiale, de façon à opérer un lissage des données mesurées. Cette étape à pour but de rendre le traitement à suivre plus robuste, grâce à l'élimination d'échos parasites éventuels. D'éventuelles fausses alarmes sont ainsi évitées. Une troisième étape 30 comprend une sélection rapide des mesures effectuées en fonction de leur rapport signal à bruit thermique couramment appelé RSB. Si ce rapport est supérieur à un seuil appelé seuil_RSB, la mesure est exploitable et donc conservée pour la suite du traitement. A l'inverse, les mesures jugées peu fiables (noyées dans le bruit thermique) sont écartées de la suite du traitement. Le seuil seuil_RSB est déterminé de manière analytique en fonction de la probabilité de fausse alarme souhaitée et du nombre d'échantillons issus de la deuxième étape. II peut par ailleurs (comme les autres seuils mentionnés dans cette demande de brevet) être adapté suivant des variables contextuelles (zone géographique, saison, altitude de l'avion, ...). Une quatrième étape 40 comprend l'application d'une première méthode dont un objectif est de détecter au sein de chaque profil vertical mesuré la portion de la puissance reçue due aux échos de sol. Plus précisément, la première méthode est appliquée pour chaque couple case distance/angle de gisement. Cette étape 40 est détaillée sur la figure 4. Elle met en oeuvre un profil vertical synthétique des échos de sol mémorisé, référencé 41 sur la figure 3. Ce dernier correspond à la puissance des échos réfléchis uniquement par le sol, suite à l'émission d'un signal radar. Ce profil vertical synthétique est élaboré lors de la conception de l'antenne radar. Il correspond à la forme du diagramme d'antenne de l'antenne radar et/ou à partir de données observées. De préférence, le profil vertical synthétique est normalisé de façon que seule sa forme soit considérée dans la suite du traitement. Une fois élaboré, le profil vertical synthétique est mémorisé en dur. Lors d'une première étape 401 de l'application de la première méthode précitée, on superpose le profil synthétique mémorisé au profil vertical mesuré traité, c'est-à-dire pour un couple case distance/angle de gisement donné. Chaque case distance de rang i, CD; est référencée par la grandeur r; correspondant à la distance entre la cellule est l'antenne radar. L'angle de gisement de rang i est référencé e;. Cette superposition 401 est illustrée plus précisément sur la figure 5. Celle-ci représente un profil mesuré PM (en trait pointillé) en cours de traitement. Les croix situées sur l'axe de l'élévation correspondent aux différents échantillons. Sur cette figure, la première méthode n'est appliquée que sur la première partie du profil mesuré, par exemple pour les N premiers échantillons (N étant un entier). En effet, les échos de sols correspondant aux élévations les plus basses, il est préférable de se limiter à la partie basse du profil. Par conséquent, les échantillons situés à droite du seuil SM (partie hachurée) ne sont pas traités. On superpose alors sur la partie basse du profil, le profil vertical synthétique mémorisé PS. Sur cet exemple, deux positions possibles du profil synthétique par rapport au profil mesuré PM sont représentées (courbe 1 et courbe 2). On note que dans le cas où un nombre insuffisant de mesures est 10 disponible, il est possible d'appliquer la première méthode sur l'ensemble du profil mesuré. On se réfère à nouveau à la figure 4. Une fois les deux profils superposés, on calcule l'erreur existant entre eux, étape 401. Cette erreur est symbolisée par des flèches verticales 15 sur la figure 5 (courbe 2). De préférence, l'erreur est calculée en décibels de façon que chaque valeur d'élévation ait le même poids. Ces étapes 401 et 402 sont répétées pour chaque position possible du profil synthétique PS par rapport ai profil mesuré PM. Une fois les erreurs calculées, on détermine quelle est l'erreur 20 minimale 402. Sur la figure 4, l'erreur est minimale lorsque le profil vertical synthétique est à la position de la courbe 2. Autrement dit, il vient : 11 v 0(, = arg min{e(çoo )} = arg min -1 Pnesure (J) û Ps}äo,edque (J' o ) v,,, (p. N j_, où: 25 -0'0 est la valeur de l'élévation pour laquelle l'erreur est minimale ; - e est l'erreur calculée pour la valeur vo de l'élévation, - N est le nombre d'échantillons considérés, -j est le rang de l'échantillon traité, - Pmesuré(j) est la puissance du profil vertical mesuré traité pour l'échantillon de rang j, et - Psynthétique(j) est la puissance du profil vertical synthétique pour l'échantillon de rang j. Lorsque le profil vertical synthétique est à cette position, cela signifie que la partie correspondante du profil vertical mesuré correspond aux échos de sol. Ainsi, de façon à obtenir un profil vertical mesuré affranchi des échos de sol, il suffit d'exclure la zone correspondante (� valeurs à exclure �). Toutefois, dans l'exemple ci-dessous, l'exclusion est réalisée plus tardivement, dans un but d'optimisation. On se contente dans une étape 403 de mémoriser une zone correspondant à la position du profil synthétique ayant une erreur minimale, dite zone d'exclusion (ou � valeurs à exclure �). De préférence, chaque zone d'exclusion est mémorisée avec une certaine marge. On se réfère à nouveau à la figure 3. Parallèlement à l'application de la méthode 1, de façon à améliorer la définition du profil vertical résultant, une deuxième méthode est appliquée au profil vertical mesuré. De même que pour la méthode 1, comme illustré sur la figure 6, on superpose le profil vertical synthétique sur le profil vertical mesuré, pour chaque valeur d'échantillon, étape 500. Toutefois, contrairement à la première méthode, l'ensemble des échantillons du profil sont considérés pour cette deuxième méthode. Pour chaque position du profil vertical synthétique par rapport au 25 profil vertical mesuré, l'erreur entre les deux profils est calculée, 501. Une fois cette erreur calculée pour l'ensemble des échantillons, l'erreur minimale est déterminée, 502. Comme représenté sur la figure 7, l'erreur entre le profil mesuré PM et le profil synthétique PS est minimale lorsque le profil synthétique PS 30 se situe au niveau du barycentre du profil mesuré PM. Ce barycentre de la puissance du profil mesuré correspond à un maximum de réflectivité dans l'environnement de l'avion. En d'autres termes, il témoigne de la présence d'un nuage, référencé par son centre. Par conséquent, la valeur de l'élévation associée à l'erreur minimale correspond au centre d'un nuage, 503. Celle-ci référencé (pn,i pour un couple case distance/angle de gisement (r;,6;), est délivrée en sortie de la deuxième méthode. La hauteur dite de sol correspondant à l'élévation du barycentre du profil est alors calculée au cours d'une étape 60. Cette hauteur correspond à la hauteur corrigée de l'effet de courbure de la terre. Les différentes grandeurs mises en oeuvre dans ce calcul sont représentées sur la figure 8. Un demi-cercle TE symbolise grossièrement la planète terre. La grandeur ha correspond à la hauteur de l'avion. Celle-ci est conservée au sein d'une mémoire interne 61 est délivrée au cours de l'étape de calcul 60. La grandeur Rr correspond au rayon terrestre, sachant que dans cet exemple on considère le rayon terrestre radioélectrique. L'angle cf) est l'élévation du faisceau radar. La grandeur r est comme indiqué ci-avant la distance entre la case 20 distance considérée et l'antenne radar. Les grandeurs (r,cp) forment le repère sphérique lié au radar. Le point P est le centre la case distance considérée associée aux coordonnées (r,(p) La hauteur hn est la hauteur du barycentre du profil. Celle-ci est 25 lié à l'élévation du barycentre du profil par la relation : h,, = r sin(çpä ) où cpn est l'élévation correspondant au barycentre du profil détecté lors de l'étape 50. La hauteur hc est la hauteur induite par la courbure de la terre. Elle 30 est définie, comme le montre la figure 8, comme la distance entre l'horizontale touchant le sol sous l'avion AV (et orthogonale au rayon de la terre en ce point) et le point d'intersection P entre la terre et la verticale passant par un point de coordonnées (r,cp) dans le repère sphérique lié au radar. Enfin, hs est la hauteur du sol à déterminer. Ces différentes grandeurs sont liées par la relation ci-dessous : hs = (ha + he ù hä )cos a avec sin a = rcosçp RT Le théorème de Pythagore nous permet d'écrire : (R, -102 +(rcosço)2 =R,2. On a donc: h, = R,. - JRT - (rcosçp) 2 . A l'issue de l'étape de calcul 60, on connait la hauteur du sol hs pour la case-distance traitée à un angle de gisement donné. A l'issue des étapes 60 et 40, les données obtenues subissent alors une post-intégration spatiale et temporelle, de façon à lisser les 15 différentes valeurs, étape 70. A l'issue de l'étape 70 peut être généré une carte de l'altitude du sol estimée 71, à l'aide notamment des différentes valeurs hs,; pour chaque couple case distance/angle de gisement (r;,8;). Cette carte peut être utilisée ultérieurement pour d'autres applications, par exemple l'établissement d'une 20 cartographie du sol (� ground mapping � en langue anglaise). Puis l'étape 40 d'une part et les étapes 50, 60 sont réitérées pour un autre couple case distance/angle de gisement (r;,8;), jusqu'à ce que l'ensemble des couples aient été traités. Une fois la totalité des données traitées, sont supprimés pour 25 chaque profil vertical, la zone d'exclusion correspondante déterminée lors de l'étape 40, ainsi que les valeurs inférieures à l'élévation associée à la hauteur dite de sol calculée lors l'étape 60. Pour cette étape 80, il est possible d'utiliser un modèle numérique de terrain mémorisé 81. D'autres critères autres que le modèle numérique de terrain peuvent être pris pour affermir l'efficacité de la méthode : le niveau de puissance des échos retournés, la corrélation spatiale, ... De même que pour l'étape précédente, les zones de sols supprimées 82 peuvent être utilisées ultérieurement pour d'autres applications, par exemple l'établissement d'une cartographie du sol (� ground mapping � en langue anglaise). A l'issue de la suppression, les profils verticaux résultants, sont alors filtrés afin de supprimer d'éventuels points parasites isolés, étape 90. Les données correspondant aux profils verticaux filtrés (sans échos de sols) peuvent alors être affichées sur un écran, étape 100. Le procédé décrit ci-dessus peut être implémenté de façon logicielle au sein d'un dispositif de traitement. Un mode de réalisation est illustré sur la figure 9. Sur cette figure est représenté à titre d'exemple, un moyen de réception REC ici de type radar, incorporant un tel dispositif de traitement DIS selon l'invention. Le récepteur REC est couplé à une antenne ANT. Les données reçues par l'antenne ANT sont ensuite filtrées par un filtre passe-bande FPB de façon à sélectionner une bande de fréquence. Dans cet exemple, une fois filtrées, les données sont alignées par un moyen d'alignement MAL. Un moyen de post-intégration MPOST effectue ensuite une post-intégration temporelle est éventuellement spatiale sur les données alignées. Un moyen de sélection MSEL couplé en sortie du moyen de post-25 intégration MPOST, est apte à sélectionner les données dont le rapport signal à bruit dit RSB est supérieur à un seuil choisi seuil_RSB. Les données sélectionnées sont alors traitées en parallèle par deux moyens de traitement référencés MMET1 et MMET2. Chaque moyen est respectivement apte à mettre en oeuvre les méthodes 1 et 2 décrites ca- 30 avant. En d'autres termes, le premier moyen de traitement MMET1 est apte pour chaque profil vertical mesuré : • à superposer le profil vertical synthétique au profil vertical mesuré considéré, pour différentes valeurs de l'angle de pointage, • à calculer l'erreur existant entre le profil vertical mesuré et le profil vertical synthétique, pour chaque valeur de l'angle de pointage considérée. Le profil vertical synthétique, tel que décrit ci-dessus est élaboré par un moyen d'élaboration MEL couplé aux moyens de traitement MMET1 10 et MMET2. Par ailleurs, le deuxième moyen de traitement MMET2 couplé à un moyen de calcul MCAL, est apte à calculer la hauteur du sol pour le profil vertical mesuré considéré. Le moyen de calcul MCAL peut être couplé à une mémoire MEMA apte à stocker la hauteur du porteur à bord duquel est 15 embarqué le récepteur radar REC. Un moyen de post-intégration MPOST effectue une post-intégration des données traitées, de façon à lisser ces dernières. Une cartographie cartol des effets du sol (� ground mapping � en langue anglaise) peut être élaborée et délivrée par le moyen de post- 20 intégration MPOST vers d'autres modules de traitement (non représenté ici à des fins de simplification). Les données intégrées sont alors délivrées à un moyen d'exclusion MEXCL apte à exclure du profil vertical mesuré : - des valeurs pour lesquelles l'erreur calculée par le premier 25 moyen de traitement MMET1 minimale, dites valeurs à exclure, de façon à générer un profil vertical résultant, affranchi des échos de sol, et - les valeurs du profil vertical inférieures à la valeur de pointage correspondant à la hauteur de sol calculée par le deuxième moyen de traitement MMET2 et le moyen de calcul MCAL. Une cartographie carto2 des zones à exclure, peut être élaborée et délivrée par le moyen d'exclusion MEXCL vers d'autres modules de traitement (non représenté ici à des fins de simplification). Le moyen d'exclusion MEXCL peut être couplé à un mémoire 5 MEMB apte à mémoriser un modèle numérique de terrain. Chaque profil vertical résultant est alors filtré par un moyen de filtrage MFIL, puis les données filtrées sont affichées sur un écran d'affichage MAFF. 10 15 20 REVENDICATIONS1- Procédé de traitement de profils verticaux mesurés de la puissance des échos renvoyés suite à une émission de signaux radar, chaque profil vertical mesuré étant fonction de l'angle de balayage du faisceau radar et associé à un couple angle de pointage/case distance donné du faisceau radar, le procédé étant caractérisé en ce qu'il comprend une élaboration d'un profil vertical synthétique (41) de la puissance des échos renvoyés par le sol seulement, et pour chaque profil vertical mesuré : • une superposition (400) du profil vertical synthétique au profil vertical mesuré considéré, pour différentes valeurs de l'angle de balayage, • un calcul de l'erreur (401) existant entre le profil vertical mesuré et le profil vertical synthétique, pour chaque valeur de l'angle de balayage considérée, et • une exclusion (80) du profil vertical mesuré, des valeurs pour lesquelles l'erreur calculée est inférieure à un seuil donné, dites valeurs à exclure, de façon à générer un profil vertical résultant, affranchi des échos de sol. 2- Procédé selon la revendication précédente, dans lequel l'angle de balayage est l'angle d'élévation et l'angle de pointage est l'angle de gisement. 3- Procédé selon l'une des revendications précédentes, dans lequel ladite confrontation n'est effectuée que pour certains points du profil vertical mesuré, correspondant à des valeurs de l'angle de balayage inférieures à un seuil donné. 4- Procédé selon la revendication 2, dans lequel ladite confrontation est effectuée pour l'ensemble des points du profil vertical mesuré, si le nombre de points du profil vertical mesuré est inférieur à un seuil donné. - Procédé selon l'une des revendications précédentes, comportant en outre un calcul de la hauteur du sol (hs1) pour le profil vertical mesuré considéré, qui comprend : • une autre superposition (500) du profil vertical 5 synthétique au profil vertical mesuré considéré, pour l'ensemble des valeurs de l'angle de balayage, • un calcul de l'erreur (501) existant entre le profil vertical mesuré et le profil vertical synthétique, pour l'ensemble des valeurs de l'angle de balayage, • une détermination de la valeur de l'angle de balayage du barycentre du profil vertical mesuré correspondant à l'erreur minimale calculée (502), • un calcul d'une hauteur dite de sol (60) pour le profil vertical mesuré considéré, à l'aide de ladite valeur de l'angle de balayage 15 déterminée, les valeurs du profil vertical inférieures à la valeur de pointage correspondant à la hauteur de sol calculée étant alors également exclues. 6-Procédé selon la revendication précédente, dans lequel lesdites valeurs à exclure et les valeurs de la hauteur de sol calculées pour un profil 20 vertical mesuré donné, sont lissées avant l'exclusion (70). 7-Procédé selon l'une des revendications précédentes, dans lequel chaque profil vertical résultant est filtré (90), de façon à éliminer les valeurs isolées. 8-Procédé selon l'une des revendications précédentes, dans lequel 25 l'exclusion est également fonction d'un modèle numérique du sol (81). 9-Dispositif (DIS) de traitement de signaux radar, couplé à une antenne radar (ANT) apte à mesurer les profil verticaux de la puissance des échos renvoyés suite à une émission de signaux radar, chaque profil vertical mesuré étant fonction de l'angle de balayage du faisceau radar et associé à 30 un couple angle de pointage/case distance donné du faisceau radar, le dispositif (DIS) étant caractérisé en ce qu'il comprend un moyen d'élaboration (MEL) d'un profil vertical synthétique de la puissance des échos renvoyés par le sol seulement, un premier moyen de traitement (MMET1) apte pour chaque profil vertical mesuré : • à superposer le profil vertical synthétique au profil vertical mesuré considéré, pour différentes valeurs de l'angle de balayage, • à calculer l'erreur existant entre le profil vertical mesuré et le profil vertical synthétique, pour chaque valeur de l'angle de balayage considérée, et un moyen d'exclusion (MEXCL) apte à exclure du profil vertical mesuré, des valeurs pour lesquelles l'erreur calculée est inférieure à un seuil donné, dites valeurs à exclure, de façon à générer un profil vertical résultant, affranchi des échos de sol. 10- Dispositif (DIS) de traitement selon la revendication précédente, comprenant en outre un deuxième moyen de traitement (MMET2) apte à calculer la hauteur du sol pour le profil vertical mesuré considéré, ledit moyen d'exclusion étant en outre apte à exclure les valeurs du profil vertical inférieures à la valeur de balayage correspondant à la hauteur de sol calculée. 11-Récepteur radar (REC) incorporant un dispositif selon l'une des 20 revendications 9 à 10. 25 30