Brevet FR2737579A1 — publié le 19970207

I Procédé de caractérisation d'un terrain survolé à partir d'un signal de sonde FM/CW. Le domaine de l'invention est celui des moyens d'aide à la navigation utilisables par des aéronefs. Plus précisément, l'invention propose un procédé d'utilisation nouvelle des signaux fournis par une sonde embarquée de type FM/CW, permettant d'effectuer une analyse des terrains survolés par un aéronef, dans le but par exemple de les identifier, ou encore de les reconnaître. Les sondes FM/CW ont des applications traditionnelles dans le domaine de l'altimétrie. Il s'agit d'émetteurs/récepteurs dont l'onde porteuse fonctionne par exemple à une fréquence de quelques GHz à quelques dizaines de GHz, et qui est modulée cycliquement par un signal modulé en fréquence. La mesure du battement entre les ondes émises et reçues fournit directement une information représenta- tive de la distance entre l'émetteur-récepteur et l'obstacle ayant réfléchi le signal. Traditionnellement, le traitement du signal de battement des sondes FM/CW a été réalisé en circuiterie analogique. Or, la densité de puissance du signal de battement produit par un tel système présente un spectre fréquentiel 11 plus ou moins étalé, comme schématisé en figure 1. Mais en fait, pour la mesure de hauteur effectuée dans l'utilisation altimétrique classique, seul une faible partie du spectre est utile, et donc effective- ment utilisée (par exemple, dans le cas simple d'une assiette de l'avion parallèle au sol, et hors effet Doppler, c'est la fréquence de battement fo correspondant au point 10 de valeur maximale Am de densité de puissance A qui est directement proportionnelle à la hauteur de l'avion par rapport au sol). En conséquence, en altimétrie classique à traitement analogique, le signal de battement est habituellement très fortement filtré de façon à éliminer les portions du spectre non exploitées, notamment la partie de décroissance 13. La présente invention se donne pour objectif nouveau d'exploiter l'ensemble de l'information contenue dans le spectre de densité de puissance du signal de battement, en considérant que le spectre est également représentatif de la nature du terrain survolé. Le spectre dépend en effet des caractéristiques de réflectivité du sol renvoyant le signal. Plus précisément, un premier objectif de l'invention est d'utiliser une sonde de type FM/CW pour réaliser l'identification du type de terrain survolé, en référence avec quelques types standards préalablement répertoriés. Un autre objectif de l'invention est de fournir un tel procédé permettant la reconnaissance d'éléments remarquables du terrain, notamment de lignes remarquables, tels que les plans d'eau, les routes, les voies ferrées, les bâtiments, les forêts, les lignes à haute tension, .. Ceci rend possible par exemple des opération de recalage de navigation par corrélation avec une carte mémorisée,..DTD: même en l'absence de relief. Un objectif complémentaire de l'invention est de fournir un tel procédé donnant une information de l'état du terrain survolé, par exemple la détermination de l'état de la mer (qui influe sur sa réflectivité) afin d'en déduire l'altitude optimale de vol d'un aéronef. Plus généralement, le procédé de l'invention a ainsi pour objectif de fournir des informations de guidage d'aéronefs, en fournissant des éléments de décision pour la navigation (par exemple un système de type à filtre de Kalman). Ces objectifs ainsi que d'autres qui apparaîtront par la suite sont atteints à l'aide d'un procédé de caractérisation d'un terrain survolé, notamment pour reconnaître un type de terrain prédéterminé, ou des transitions ou séquences de transitions remarquables sur ledit terrain, caractérisé en ce qu'il consiste à émettre en direction du terrain un signal à l'aide d'une sonde de type FM/CW, et à analyser, dans le domaine fréquentiel, la forme du spectre du signal de battement généré par la sonde pour caractériser ledit terrain. Selon l'invention, le procédé consiste à extraire de la forme du spectre dudit signal de battement des informa- tions relatives à la réflectivité dudit terrain. Préférentiellement, l'invention consiste à élaborer une quantité représenta- tive de la réflectivité du terrain survolé à partir d'au moins l'une des deux mesures suivantes: - la mesure de l'amplitude du spectre du signal de battement, comme étant représentative de la réflectivité à la perpendiculaire du sol; - la mesure de l'étalement en fréquence (inverse de la traînée) du spectre du signal de battement, comme étant représentative du diagramme de rétrodiffu- sion du terrain survolé. Avantageusement, on élabore préalablement des valeurs de référence, de façon à classer des spectres de signal de battement selon au moins trois catégories de terrains: - les terrains spéculaires, - les terrains moyennement diffusants, - les terrains fortement diffusants. Dans un mode de réalisation avantageux de l'invention, le procédé consiste à élaborer une signature du terrain survolé, à partir d'une part de ladite informa- tion de réflectivité, et d'autre part de l'information de hauteur par rapport au sol fournie par la sonde. Le procédé prévoit en outre de pouvoir utiliser comme critère de reconnaissance du spectre du signal de battement une modélisation de la décroissance dudit spectre. Dans ce cas, ladite modélisation de la pente du spectre du signal de battement est avantageusement réalisée à l'aide d'une fonction du type: y = 1/u2'l, avec p: indice de lobe d'antenne q indice de diagramme de rétrodiffusion du sol N=p+q+2 u = f/fo (fréquence de battement f réduite par la fréquence de référence fo correspondant à la réflexion à la normale au sol, hors effet Doppler; fo = 2 (h/c). (A FT)) Le coefficient N issu de la modélisation donne donc une valeur représenta- tive de q que l'on peut multiplier avec les mesures d'amplitude et d'étalement du spectre pour former un coefficient composite d'analyse de signature. La reconnaissance de la forme du spectre du signal de battement peut également comprendre la prise en compte d'au moins l'un des critères suivants: - amplitude maximale du pic principal du spectre, - mesure du contraste (amplitude maximale / amplitude moyenne), - raideur du front avant du spectre, - nombre de pics secondaires. Pour l'analyse du spectre du signal de battement, il est également avantageux de pondérer les critères de reconnaissance du spectre du signal de battement avec au moins l'un des critères complémentaires suivants: assiette de l'avion, - effet Doppler dû à la vitesse de déplacement par rapport au terrain, - diagramme de l'antenne de la sonde. Selon l'invention, la caractérisation du terrain peut également être effectuée en prenant en compte les fluctuations temporelles du spectre du signal de battement. Dans un mode de réalisation complémentaire de l'invention, la reconnais- sance dudit terrain est effectuée en combinant les informations extraites de deux spectres de signal de battement provenant de deux sondes travaillant dans des bandes de fréquence distinctes. D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture suivante de la description du principe détaillé et de divers modes de réalisation de l'invention, donnés à titres illustratifs et non limitatifs, et des dessins annexés dans lesquels: - la figure 1 schématise la forme typique du spectre du signal de battement tel qu'obtenu par une sonde de type FM/CW; - la figure 2 représente le signal en dents de scie modulant linéairement l'onde hyperfréquence d'une sonde FM/CW, de façon à illustrer le mode de production du signal de battement; - la figure 3 schématise le principe d'exploration du sol par une sonde FM/CW embarquée sur un avion, à partir duquel est mis en oeuvre le procédé de l'invention; - les figures 4A, 4B représentent respectivement les spectres successifs (fig 4A), et leur retranscription avec une amplitude binarisée (fig 4B), relevée par une sonde FM/CW embarquée sur un avion survolant un échelon de falaise dans le sens terre-mer; - les figures SA, 5B représentent respectivement le profil de hauteur et l'évolution de l'indice de réflectivité obtenus à l'aide d'une sonde FM/CW centimétrique à 4,3 GHz, lors du survol d'une digue et de sa retenue d'eau; - les figures 6A, 6B représentent respectivement le profil de hauteur et l'évolution de l'indice de réflectivité obtenus à l'aide d'une sonde FM/CW millimétrique à 60 GHz, lors du survol d'une digue et de sa retenue d'eau; - la figure 7 est un schéma permettant d'illustrer différents traitements de reconnaissance de signature réalisables sur le spectre d'un signal de battement; - les figures 8A à 8D représentent la signature théorique obtenue pour différentes configurations typiques de terrain. - la figure 9 schématise les deux étapes essentielles du procédé de l'invention: analyse spectrale, puis mise en oeuvre d'algorithmes de reconnaissance. Les sondes altimétriques FM/CW fonctionnent suivant le principe illustré en figure 2. Une onde hyperfréquence, modulée linéairement en fréquence par un signal en dents de scie 21, est émise vers le sol. Le récepteur, de type homodyne reçoit l'onde 22 réfléchie par le sol et la mélange avec l'onde émise. La fréquence de battement qui résulte d'une réflexion d'un point 31 du sol situé à une distance h de l'aéronef (voir fig.3) vaut (pour un effet Doppler nul): fb= Z. Fr/T avec: a F: fréquence de battement (fb) T: période du signal en dents de scie C: vitesse de la lumière = 2h/c Le signal de battement obtenu sur l'ensemble des réflexions issues du sol et reçues dans le lobe d'antenne est composé d'un ensemble de composantes sinusoïdales d'amplitudes différentes. L'enveloppe 11 du spectre d'un signal de battement est représentée par exemple en figure 1. La connaissance de la hauteur minimale ho entre l'aéronef 32 et le sol 30 est alors déduite de l'analyse spectrale numérique du signal de battement. Les spectres du signal de battement sont par exemple obtenus à un rythme moyen de 250 par seconde et sont composés d'environ 250 raies élémentaires de largeur égale à 500 Hz. Les critères d'extraction de la hauteur sont principalement basés sur la détection du "front avant" 12 (figure 1) de la répartition spectrale, lequel correspond sur terrain plat et en l'absence de roulis de tangage, à la hauteur ho à la verticale de l'aéronef. Dans ce type de fonctionnalité, la signature spectrale 11 n'est pas complètement utilisée. Selon l'invention, la reconnaissance de la nature du terrain survolé et plus particulièrement de sa réflectivité, suppose d'utiliser davantage la signature spectrale obtenue. Plus précisément, et comme schématisé en figure 9, la mise en oeuvre de la reconnaissance s'effectue avantageusement en deux temps: - analyse spectrale 91 du signal de battement Sb(t) - application d'algorithmes 92 de reconnaissance de la signature spectrale Sb(f)O. L'analyse spectrale peut être mise en oeuvre: - soit par technique conventionnelle comme par exemple le "corrélogram- me" (transformation fréquentielle du signal après moyennage de façon à minimiser l'effet du bruit). C'est l'application de cette technique conduisant à l'obtention d'une signature spectrale qui est détaillée ci-après. - soit par technique paramétrique, basée par exemple sur la prédiction linéaire; dans ce cas, on élabore un jeu de p coefficients de prédiction linéaire (p: ordre de mobilisation) à l'aide desquels il est possible de caractériser le terrain survolé sans passer par le spectre. Les algorithmes de reconnaissance peuvent être: - soit des techniques de caractérisation reposant sur la forme du spectre (comme détaillé ci-après en relation notamment avec les figures 7 et 8, - soit par technique de reconnaissance sans caractérisation. Une telle technique suppose par exemple d'utiliser un réseau neuronal supervisé: on réalise un apprentissage du système avec injection de signaux de battement strictement identifiés (de mer, de routes, de chemins, de bois. ) et le système peut ensuite détecter tout seul toute nouvelle excitation. Un sol (homogène) est caractérisé par sa réflectivité a(e), fonction de l'angle d'éclairement e, égale à oa0'(e) o o0 est la réflectivité à la perpendiculaire du sol et O'(e) correspond au diagramme de rétrodiffusion. On considère classiquement a'(e) = cos2q(e) o q caractérise le diagramme de rétrodiffusion..DTD: (modèle de Mulhemann). Les valeurs des paramètres oa et q sont étroitement reliées à la spécularité de la cible: si oa est grand alors q est important; de même si oao est faible alors q est faible. A titre d'exemple, on trouve aux extrêmes la forêt avec q = 2 (terrain très diffusant) et la mer calme avec q = 32. Un critère de reconnaissance peut donc être constitué en multipliant les valeurs de q et a,. Le problème consiste à caractériser, le type de terrain survolé, et notamment les valeurs des paramètres ao et q de ce terrain à travers la signature spectrale reçue par la sonde. Ceci peut être illustré par les figures 4A et 4B reproduisant graphiquement les mesures obtenues correspondant au passage d'un avion muni d'une sonde FM/CW à la verticale d'un échelon de falaise dans le sens terre (terrain diffusant) vers mer (spéculaire). La figure 4A représente les spectres successifs obtenus toutes les 22 ms sur lesquels on discerne en abscisse (axe des fréquences, donc des distances) la signature 41 peu énergétique mais très étendue de la falaise centrée sur 64 kHz et la réponse 42 très énergétique et étroite de la mer centrée sur 96 kHz. La différence entre les fréquences 64 kHz et 96 kHz est représentative de la différence de niveau de la mer et de la falaise, et donc de la différence de hauteur de l'avion par rapport à ces deux "terrains" survolés. La figure 4B représente le même scénario, avec le temps en abscisse, la fréquence en ordonnée et l'amplitude binarisée après seuillage des signaux de la Fig 4A, procurant ainsi une image en noir et blanc. Chaque spectre correspond donc à une colonne verticale vue du dessus. On remarque la différence d'étendue spectrale et de fluctuations temporelles qui sont caractéristiques des deux types de terrain, et permettent d'identifier la transition entre un type de terrain et l'autre. On remarquera en outre que la prolongation 43 de la trace 44 du signal de terre, alors que le signal de mer 45 indique que l'avion survole déjà l'étendue maritime, est caractéristique du phénomène d"'oblicité" résultant de l'effet réfléchissant du front de falaise surplombant la mer. On notera que dans le cas d'un terrain spéculaire, la puissance reçue ne dépend pratiquement plus du terrain puisque l'on considère que toute l'énergie émise est supposée provenir d'un point image de l'émetteur par rapport au sol. On définit alors une réflectivité équivalente o. = Gr/4 o G représente le gain d'antenne et r le coefficient de réflexion du sol. Ceci explique pourquoi avec une antenne à pinceau fin donc à grand gain, la dynamique d'amplitude des échos reçus est importante. Selon une caractéristique importante de l'invention, le procédé de caractérisation d'un terrain (homogène) survolé permet d'extraire les coefficients (o0,q) du terrain, à partir de la densité spectrale du signal de battement obtenu en présence d'effet Doppler. On peut démontrer que le centre de la densité spectrale de puissance du signal de battement fournit effectivement une mesure ou une tendance des coefficients o et q. Plus précisément, on peut montrer que l'ensemble des points du sol contribuant à donner l'énergie en une fréquence du spectre (ensemble appelé isofréquence) approche un cercle d'équation h2 = a+b cos avec a = r2 + ho2 + h2ax2cos2e et b= 2rho a cos 6 dont le centre et le rayon sont paramétrés par: - a (fréquence Doppler normalisée): a = 2VR/W fo avec VR: vitesse relative entre l'avion et le point réflecteur fo: fréquence de battement due au trajet à la verticale en l'absence d'effet Doppler à la verticale A: longueur de l'onde utilisée. - 6: inclinaison du vecteur vitesse par rapport à l'horizontale. Les différentes isofréquences constituent donc une famille de cercles concentriques. La densité spectrale de puissance S(u) à la fréquence réduite u(=f/fo) est obtenue en intégrant la puissance réfléchie par chaque élément de surface de l'isofréquence u, compte-tenu du diagramme d'antenne et de la variation de la réflectivité du sol avec l'incidence. On peut considérer deux cas: - a est petit, ce qui correspond au cas o la fréquence de battement f0 est importante et/ou la longueur d'onde est grande (ou encore la vitesse de l'aéronef est faible). Dans ce cas, on peut montrer que S(u) = aO2.(1/u21); la densité spectrale possède un maximum pour umi, (valeur du minimum en fréquence de battement normalisée). La décroissance de S(u) lorsque u augmente dépend presque uniquement du coefficient q de réflectivité du terrain pour une antenne donnée. Plus le terrain est diffusant (q faible) et plus le spectre est étalé. D'autre part, la densité spectrale est directement proportionnelle à ao0. - a est proche de 1, ce qui correspond au cas o la fréquence de battement f0 est faible et/ou la longueur d'onde est petite. Dans ce cas, le calcul est plus complexe et la signature est d'autant plus riche que le lobe d'antenne est important, que le terrain est diffusant. On peut en effet montrer que la tendance de la densité S(u) est donnée par: S(u) = oJ2(N-1).:k.0,.) C2k+ N1 (2a) (u2N)/(u +4ar2) avec N = p + q + 2. Malgré l'interférence d'autres paramètres, la densité spectrale reste proportionnelle à a. et sa décroissance est liée au coefficient q du terrain. La forme du spectre dépend, de plus, de la vitesse horizontale VR (dans a). Le principe de l'invention consiste à construire un indice de la réflectivité du terrain représentatif des coefficients caractéristiques ao et q de ce terrain. La complexité du problème est la conséquence des difficultés d'établissement de la forme théorique du spectre à obtenir, de l'interaction de paramètres étrangers au terrain comme le lobe d'antenne, de l'influence des attitudes avion et du vecteur vitesse, de l'influence du relief du terrain. L'identification du type de terrain repose plus particulièrement sur la mesure ponctuelle d'une quantité représentative du coefficient de rétrodiffusion moyen du terrain (indice de terrain). Ce coefficient peut être par suite moyenné ou plus généralement filtré selon une fonction de transfert du type passe-bas. L'indice est déduit de chacun des spectres instantanés du signal de battement par une technique de caractérisation du spectre, s'apparentant plus largement à la reconnaissance de forme. Les essais réalisés ont permis de montrer que la séparation obtenue entre sols à caractéristiques diffusantes et sols à tendance spéculaire est nette et sans recouvrement. Une optimisation basée sur l'influence de chacun des paramètres mesurés et débouchant sur la pondération de chacun d'eux dans la formule utilisée, permet d'améliorer encore ces résultats. Par extension, la détection de points remarquables ou encore de "coupures et de textures particulières" fait appel à l'examen temporel des mesures ponctuelles d'indice de terrain. C'est à la vue de l'évolution des indices successifs que l'on peut tenter d'identifier lesdites "coupures" de terrain (passage d'une zone diffusante à une zone spéculaire ou l'inverse par exemple) et plus encore de reconnaître les structures remarquables telles que franchissement de cours d'eau, de routes, d'autoroutes, de bois, de forêts, de lignes à haute tension... en s'aidant le plus possible de données complémentaires comme les variations simultanées de hauteur détectée, les longueurs de zones déduites,... Les critères de caractérisation de terrain par rapport à la forme du spectre du signal de battement peuvent être notamment les suivants: essais qualitatifs d'identification de terrain. L'indice utilisé est construit simplement comme le produit de l'amplitude maximum du spectre par l'inverse de l'étendue de la signature spectrale au-dessus d'un seuil (étendue appelée "traînée"). Cet indice est donc proportionnel à oao et relié à q. Il est particulièrement utile pour la détection de franchissement de zones à réflectivités différentes. - une mesure directe de oa peut être réalisée, après calibrage, en il pondérant l'amplitude du spectre avec la valeur de gain affectée à la chaîne de transmission du signal FM/CW, particulièrement dans la chaîne de réception (le gain à l'émission étant généralement standardisé). - prise en compte de la hauteur détectée comme information com- plémentaire de détection de zones remarquables, comme par exemple dans le cas de transitions clairière-bois, de digues, de lignes à haute tension. - étude de variables classifiantes complémentaires à partir de la forme théorique du spectre à obtenir même dans le cas de a élevé. On peut ainsi se référer aux figures 8A, 8B, 8C, 8D représentant respectivement la signature théorique obtenue dans les cas: - d'un sol plan horizontal spéculaire (fig 8A) - d'un sol plan incliné spéculaire (fig 8B) - d'un sol plan incliné diffusant (fig 8C) - d'un passage d'échelon (fig 8D). Comme illustré en figure 7 représentant la densité spectrale de puissance (y) en fonction de la fréquence, les critères de reconnaissance du spectre peuvent passer par des seuillages 71, des détections de maxima 72, 73, des mesures de contraste 74 contre l'amplitude maximale A,. et la valeur moyenne Vm, un filtrage de type barycentre, la mesure de la raideur r du front avant 75, l'étendue de la trainée 76, etc... Un sol est diffusant si contraste et raideur sont faibles et la trainée importante. Un sol est spéculaire si contraste et raideur sont importants et la trainée faible. - prise en compte des dégradations apportées par les différentes assiettes de l'aéronef. - prise en compte des fluctuations temporelles des spectres successifs pour caractérisation complémentaire des terrains survolés (les terrains spéculaires ont généralement une signature plus stable). Ainsi, par exemple, par la mesure de l'écart-type de fluctuation de ce maximum, on peut donc améliorer encore la fiabilité de l'identification effeectuée. - dérivée maximale du spectre. - modélisation de la décroissance du spectre (notamment par une technique de moindres carrés). A titre d'exemple, il est possible de modéliser le spectre, après seuillage, sous la forme y l/u2'1, avec p: indice de lobe d'antenne q: indice de diagramme de rétrodiffusion du sol N =p + q + 2 u = f/f0 (fréquence de battement f réduite par la fréquence de référence fo = 2 (h/c). (,F/T) Dans tous les cas, on a intérêt à travailler avec une fréquence d'échantillo- nnage élevée pour affiner la reconnaissance du spectre. On présente en figures SA et 5B les résultats obtenus à l'aide d'une sonde altimétrique classique à 4,3 GHz (À = 6 cm) et à antenne large (60 dans les deux axes). La longueur d'onde utilisée situe a proche de zéro (a=o,1 pour Vmax = 200 nm/s et fo = 80 KHz). En conséquence la forme du spectre peut s'exprimer en 1/un. En figure SA, on représente le profil de hauteur obtenu au passage au- dessus d'une digue et de la retenue d'eau qui s'y attache. La hauteur moyenne de vol est de 250 pieds. L'espacement entre chacun des spectres est de 5ms et la durée totale analysée est de 1,25 seconde, correspondant à un parcours d'une centaine de mètres. En figure 5B on trouve l'évolution de l'indice de réflectivité défini précédemment au cours des analyses successives. On présente en figures 6A et 6B les résultats obtenus avec une sonde à 60 GHz (r4 = 5 mm) et à antenne étroite (5 ). La longueur d'onde utilisée permet à a d'approcher 1. Dans l'expérimentation présente, la vitesse est d'environ 90 m/s et a vaut donc 0,5. En conséquence, la forme du spectre relève plutôt du cas n 2 traité plus haut (ca proche de 1), avec une forte influence du paramètre p (lobe d'antenne étroit). En figure 6A, on représente le profil de hauteur obtenu sur le même terrain (expérimentation simultanée) à la même hauteur. En figure 6B, on trouve l'indice de réfiectivité et son évolution temporelle. On constate que les mesures des deux types de sonde sont complémentai- res - la sonde FM/CW "centimétrique" à antenne large (Fig 6A,6B) permet une mesure correcte de la hauteur quelle que soit l'attitude de l'avion et permet bien l'identification de zones diffusantes étendues. - la sonde FM/CW "millimétrique" à pinceau étroit (Fig 6C,46) éclairant au sol une tache de faible étendue est plus sensible à l'attitude de l'avion, notamment aux grands angles pour la détermination précise de la hauteur. La dynamique d'amplitude des spectres est très importante et permet une détection aisée des zones à réflexion spéculaire de petite étendue et par ce fait des coupures de terrain (routes, chemins, cours d'eau, voies ferrées, lignes à haute tension...). REVENDICATIONS 1. Procédé de caractérisation d'un terrain survolé, notamment pour reconnaître un type de terrain prédéterminé, ou des transitions ou séquences de transitions remarquables sur ledit terrain, caractérisé en ce qu'il consiste à émettre en direction du terrain un signal à l'aide d'une sonde de type FM/CW, et à analyser, dans le domaine fréquentiel, la forme du spectre du signal de battement généré par la sonde pour caractériser ledit terrain. 2. Procédé selon la revendication 1 caractérisé en ce qu'il consiste à extraire de la forme du spectre dudit signal de battement des informations relatives à la réflectivité dudit terrain. 3. Procédé selon l'une quelconque des revendications 1 et 2 caractérisé en ce qu'il consiste à élaborer une quantité représentative de la réflectivité du terrain survolé à partir d'au moins l'une des deux mesures suivantes: - la mesure de l'amplitude du spectre du signal de battement, comme étant représentative de la réflectivité à la perpendiculaire du sol; - la mesure de l'étalement en fréquence du spectre du signal de battement, comme étant représentative du diagramme de rétrodiffusion du terrain survolé. 4. Procédé selon l'une quelconque des revendications 2 et 3 caractérisé en ce qu'il consiste à élaborer des valeurs de référence, de façon à classer des spectres de signal de battement selon au moins trois catégories de terrains: - les terrains spéculaires, - les terrains moyennement diffusants, - les terrains fortement diffusants. 5. Procédé selon l'une quelconque des revendications 2 à 4 caractérisé en ce qu'il consiste à élaborer une signature du terrain survolé, à partir d'une part de ladite information de réflectivité, et d'autre part de l'information de hauteur par rapport au sol fournie par la sonde. 6. Procédé selon l'une quelconque des revendications 1 à 5 caractérisé en ce qu'il consiste à utiliser comme critère de reconnaissance du spectre du signal de battement une modélisation de la décroissance (13) dudit spectre. 7. Procédé selon la revendication 6 caractérisé en ce que ladite modélisa- tion de la pente du spectre du signal de battement est réalisée à l'aide d'une fonction du type: y = 1/u2N1, S avec p: indice de lobe d'antenne q: indice de diagramme de rétrodiffusion du sol N=p+q+2 u = f/f (fréquence de battement f réduite par la fréquence de référence fo correspondant à la réflexion à la normale au sol, hors effet Doppler). 8. Procédé selon l'une quelconque des revendications 1 à 7 caractérisé en ce que la reconnaissance de la forme du spectre du signal de battement comprend le prise en compte d'au moins l'un des critères suivants: amplitude maximale du pic principal du spectre, - mesure du contraste (amplitude maximale / amplitude moyenne), - raideur du front avant du spectre, - nombre de pics secondaires. 9. Procédé selon l'une quelconque des revendications 1 à 8 caractérisé en ce que les critères de reconnaissance du spectre du signal de battement sont pondérés par au moins l'un des critères complémentaires suivants: attitude de l'avion, - effet Doppler dû à la vitesse de déplacement par rapport au terrain, - diagramme de l'antenne de la sonde. 10. Procédé selon l'une quelconque des revendications 1 à 9 caractérisé en ce que la caractérisation du terrain est effectuée en prenant en compte les fluctuations temporelles du spectre du signal de battement. 11. Procédé selon l'une quelconque des revendications 1 à 10 caractérisé en ce que la reconnaissance dudit terrain est effectuée en combinant les informations extraites de deux spectres de signal de battement provenant de deux sondes travaillant dans des bandes de fréquence distinctes.