Brevet FR2901613A1 — publié le 20071130

ci, les solutions suivantes peuvent être utilisées seules ou de manière complémentaire les unes des autres. Par exemple, une seconde solution consiste à utiliser des traitements adaptatifs à haute résolution sur le signal reçu. Ces traitements présentent le double inconvénient d'être difficiles à mettre en oeuvre et de donner de mauvais résultats si les hypothèses faites sur le modèle varient, même très légèrement. En outre, ils exigent une capacité de calcul importante, ce qui augmente le coût du système. Une troisième solution réside dans l'augmentation de la résolution distance et doppler. Cette augmentation de résolution ne peut cependant être effectuée que si le radar est positionné à un niveau suffisamment élevé par rapport au sol. Ceci n'est pas toujours possible compte tenu du contexte d'emploi qui contraint fortement la hauteur des antennes radar, toujours par exemple, pour un radar mobile. Une quatrième solution vise à augmenter la diversité des fréquences d'émission. Néanmoins, diversifier les fréquences nécessite de disposer d'une plage de fréquences importante afin que ceci soit efficace. Ce qui est souvent incompatible avec l'application en jeu et n'est, par ailleurs, pas toujours autorisé par les organismes d'allocation des fréquences. Un but de l'invention est notamment de pallier les inconvénients précités. A cet effet, l'invention a pour objet un procédé de détermination de l'élévation d'une cible proche du sol par un système de détection électromagnétique. Ce procédé extrait la mesure de l'élévation à partir de l'amplitude du signal d'interférence produit par un signal émis directement par la cible et par un signal émis par la cible en direction du sol puis réfléchi par le sol vers le radar. Le système émet un signal de détection vers la cible. Le signal d'interférence est produit par le signal ré-émis directement par la cible et le signal ré-émis par la cible en direction du sol puis réfléchi par le sol vers le système de détection. L'antenne de réception est constituée par exemple de réseaux d'éléments rayonnants. Le signal d'interférence reçu par cette antenne peut être échantillonné au niveau des sous-réseaux. La mesure de l'élévation est alors établie à partir de l'amplitude du signal d'interférence échantillonné. Les sous-réseaux correspondent par exemple chacun à une ligne d'éléments rayonnants et peuvent être disposés de manière horizontale, parallèle et uniforme sur l'ensemble de l'antenne. Les sous-réseaux peuvent se situer également dans un plan sensiblement vertical. La fréquence f du signal échantillonné, est par exemple la valeur qui minimise la fonction suivante : J (U ) = I ri+k + rr-k ù 21"1 cos (27z kA i,k avec • r;+k représentant le résultat de la soustraction entre la valeur du signal reçu au niveau du sous-réseau d'ordre i+k et la valeur du signal reçu 10 au niveau du sous-réseau d'ordre i+k+1, r;_k représentant le résultat de la soustraction entre la valeur du signal reçu au niveau du sous-réseau d'ordre i-k et la valeur du signal reçu au niveau du sous-réseau i-k-1 • A représentant la distance entre les centres de phases de deux sous-15 réseaux consécutifs. Le rapport entre la hauteur h2 de la cible par rapport au sol et la distance d projetée au sol entre le radar et la cible s'exprime en fonction de la fréquence f du signal d'interférence calculée selon la relation suivante : h 2/d = /f / 2 où X est la longueur d'onde du signal émis. 20 Avantageusement la hauteur d'élévation h2 peut être déterminée pour un sol non plan et non horizontal à partir : • du point de réflexion R du signal ré-émis par la cible vers le sol sur ce dernier, • de la hauteur d'élévation déterminée pour un sol plan et horizontal. 25 L'expression de l'égalité entre les angles d'incidence du rayon ré-émis par la cible vers le sol et celui réfléchi par le sol avec le plan tangent au sol au point de réflexion permet d'obtenir les coordonnées (XR, ZR) du point de réflexion R ainsi que l'angle aR entre le plan tangent au sol au point R et l'horizontale à partir de la fonction de coût suivante : 30 C(r) = h, -zR +xR xaR h2 x R d - xR où h1 représente la hauteur par rapport au sol de l'antenne, h2 représente l'estimation de la hauteur d'élévation calculée dans le cas d'un sol plan et horizontal, les valeurs (xR, ZR, OR) minimisant cette fonction sont les coordonnées du point de réflexion R. L'obtention de la valeur de l'angle îR entre le plan tangent en R au terrain et l'horizontale permet notamment de corriger la hauteur de la cible par rapport au plan tangent en une valeur H2 de la manière suivante : H2 = h2 + zR + (d ù xR) x aR, h2 étant la hauteur calculée de la cible par rapport au sol dans le cas d'un sol plan et horizontal et d étant la distance 10 projetée au sol de la cible à l'antenne. Le traitement selon l'invention peut ensuite être amélioré par un traitement CHA. L'antenne utilisée peut être également une antenne lacunaire. La cible peut comporter un émetteur actif ou bien se comporter comme un 15 réflecteur passif. La cible peut être un drône en phase d'atterrissage. La comparaison de l'amplitude du coefficient de réflexion avec un seuil fixé est par exemple préalablement effectuée. La hauteur h2 de la cible peut alors être déterminée par un traitement de type monopulse si la valeur de l'amplitude du coefficient de réflexion est inférieure à ce seuil. 20 L'invention a notamment pour principaux avantages qu'elle s'adapte à des modèles de radars existants, qu'elle peut être mise en oeuvre sur différents types d'antennes, qu'elle utilise un traitement simple à mettre en oeuvre et qu'elle nécessite des moyens de calcul peu important. 25 D'autres caractéristiques et avantages de l'invention apparaîtront à l'aide de la description qui suit, faîte en regard des dessins annexés qui représentent : • la figure 1, un aéronef en phase d'atterrissage où figurent les 30 éléments d'un dispositif utilisé par le procédé de l'invention, • la figure 2, le phénomène de réflexion de l'image de la cible consécutivement à l'émission d'une onde radar, • la figure 3, le diagramme d'antenne constitué par la cible et son image, 35 • la figure 4a, le signal d'interférence reçu au niveau de l'antenne, • la figure 4b, les amplitudes minimums et maximums du signal d'interférence reçu au niveau de l'antenne. • la figure 5, les différentes phases de traitement du procédé, • la figure 6, un exemple d'échantillonnage du signal reçu en amplitude sur 1/2 période, • la figure 7, un exemple d'échantillonnage du signal reçu en amplitude sur y 1eme de période, • la figure 8a, le phénomène de réflexion dû à l'image de la cible dans le cas d'un sol horizontal, io • la figure 8b, le phénomène de réflexion dû à l'image de la cible dans le cas où le sol n'est pas horizontal, • la figure 9, le phénomène de réflexion dans le cas où le sol n'est ni plan ni horizontal. 15 La figure 1 présente un aéronef 1 en phase d'atterrissage sur une piste 2. Cet aéronef, par exemple un drône, est guidé par un système composé d'un radar 3, d'une balise 4 embarquée à bord du drône, et d'une balise 5 au sol. Ce système de guidage est notamment décrit dans la demande de brevet français n 04 12313. Pour participer au guidage de 20 l'aéronef 1, le radar doit connaître la position de celui-ci, c'est-à- dire ses coordonnées en azimuth et en élévation. A cet effet, le radar émet un signal 6 en direction de l'aéronef 1, que l'on appellera cible par la suite. La cible peut se comporter soit comme un réflecteur passif soit comme un émetteur grâce à la balise 4 embarquée. La description de la présente 25 invention se place dans le cas de l'utilisation d'un émetteur à bord de la cible, mais l'invention peut aussi bien s'appliquer au cas d'un émetteur passif. Il a été indiqué précédemment que des signaux ré-émis par la cible pouvaient par ailleurs se réfléchir sur le sol avant d'être captés par le radar. L'état de la piste 2 peut être tel que la réflexion est quasi-spéculaire, ces 30 réflexions produisant des mesures aberrantes. La figure 2 illustre les chemins des signaux ré-émis par la cible 1 de manière directe et indirecte. Sur cette figure, seule l'antenne 21 du radar 3 est représentée. Cette antenne peut être utilisée en émission/réception ou 35 seulement en réception. Dans cet exemple, l'antenne 21 est une antenne réseau dite complète c'est-à-dire qu'elle est constituée de N sous-réseaux 22 espacés régulièrement dans le plan vertical et respectant les conditions d'interdistance nécessaires afin de ne pas créer de lobes ambigus. La présente invention peut par ailleurs s'appliquer également au cas d'une antenne lacunaire. Chacun des sous-réseaux de l'antenne 21 correspond par exemple à une ligne de plusieurs éléments rayonnants, dont la résultante correspond à un centre de phase dans le plan élévation. Sur la figure 2, les sous-réseaux sont représentés par des lignes constituées d'éléments rayonnants qui reçoivent un signal de même phase selon le plan d'élévation. Ces lignes peuvent être horizontales comme c'est le cas sur la figure 2 ou bien non horizontales. Par la suite, les lignes sont considérées comme horizontales. Afin de simplifier le schéma, seuls les rayonnements émis et/ou reçus par un sous-réseau 23 sont représentés. La ligne 23 correspondante est située à une hauteur hi par rapport au sol. Cette ligne émet et/ou reçoit un signal 24 en direction de la cible 1. La cible, quant à elle, est située à une hauteur h2 du sol et à une distance projetée au sol d du radar. Le signal résultant, ré-émis par la cible, se compose entre autre des signaux suivants captés par la ligne 23 : • d'une part un signal direct 25 qui part en direction de la ligne 23 et forme avec l'horizontale un angle 6c aussi nornmé angle d'élévation de la cible vue du radar, • d'autre part un signal 26 lui-même réfléchi par le sol 30. Ce dernier signal 26 est réfléchi de manière quasi-spéculaire par le sol. Il se réfléchit en un signal 27 en direction de la ligne 23. L'angle 0; formé entre le signal 27 réfléchi par le sol et l'horizontale est appelé angle d'élévation de l'image vue par le radar. En effet, le signal 27 réfléchi par le sol peut être assimilé à un signal 28 émis par l'image 29 de la cible par rapport au sol. Le système ainsi constitué par la cible et son image forme un dipôle cible- image. Le signal capté par la ligne 23 se compose donc du signal 25 ré-émis par la cible et du signal 27 ré-émis par l'image. La figure 3 représente un diagramme d'antenne à large faisceau du type radar de veille. Plus particulièrement, cette figure illustre le diagramme d'antenne modélisé correspondant à un dipôle composé de la cible 1 et de son image 29. Cette modélisation revient à interpréter la totalité du signal, de forme sinusoïdale, provenant du bipôle cible-image au niveau de l'antenne comme s'il s'agissait d'un signal ré-émis par un système ponctuel situé au centre du dipôle. Les grandeurs physiques caractérisant le diagramme à un instant donné sont les suivantes : • d : la distance projetée au sol entre la cible et le radar, • 0 : l'angle de visée par rapport au plan horizontal du sol, • h1 : la hauteur par rapport au sol de la ligne considérée sur l'antenne, • h2 : la hauteur de la cible par rapport au sol. Le diagramme d'interférence résultant de la modélisation du signal ré-émis décrite précédemment est de la forme, à une constante multiplicative près : 11(o) 12= 1+p2+2pcos((4iih2sin(0))/X)$) (1) 15 Où • p représente l'amplitude du coefficient de réflexion sur le sol, • 4 est un déphasage entre le signal incident et le signal ré-émis, • est la longueur d'onde du signal émis. Etant donné que l'antenne radar se trouve proche du sol, que la cible est en 20 phase d'atterrissage, et donc qu'elle se trouve également proche du sol, sa hauteur par rapport au sol est donc faible devant la distance d et on peut donc considérer que l'on a de faibles valeurs pour l'angle 0. Ceci permet donc de simplifier l'équation (1) afin d'obtenir pour l'interférence la formule suivante, en assimilant sin(0) à 0 : 25 11(0)12=1 +p2+2pcos((47th20)/X)+$) (2) Si h1 est la hauteur d'une ligne de l'antenne, toujours pour de faibles valeurs de 0, on a la relation suivante: h1=dx0 (3) 30 Ce qui donne pour l'équation (2) la forme suivante : 1I(h1)12= 1+p2+2pcos((47th2h1)/Xd)+ $4)) (4) 35 Cette relation (4) exprime l'allure de la courbe d'interférence en fonction de la hauteur h1. La figure 4a illustre le signal d'interférence 41 produit par le signal direct 25 et le signal réfléchi 27 au niveau de l'antenne 21 tel qu'exprimé par la relation (4). L'invention exploite ce signal d'interférence 41 présent au niveau de l'antenne. Plus particulièrement, le procédé selon l'invention extrait la hauteur h2 à partir de ce signal d'interférence. Ce signal d'interférence est composé d'une partie constante 1+p2 et d'une sinusoïde fonction de la hauteur de la cible h2. Le procédé selon l'invention utilise avantageusement l'existence d'un tel io signal d'interférence 41 au niveau de l'antenne, l'amplitude et la période de ce signal n'étant pas connues à priori. Pour caractériser ce signal, un échantillonnage 42 : S1, S2,, SN, est effectué au niveau des lignes rayonnantes 22, 23 de l'antenne. L'échantillonnage peut être effectué de manière régulière dans l'espace si la distance A entre deux lignes 15 rayonnantes consécutives est constante. Il peut également être adapté à un réseau de lignes non uniforme. L'échantillonnage peut donc dépendre de la structure de l'antenne. La figure 5 illustre les différentes étapes possibles pour la 20 détermination de la hauteur h2 à partir du signal d'interférence 41. Une première étape 51 permet de déterminer la valeur de l'amplitude du coefficient de réflexion p sur le sol. Les amplitudes minimum et maximum du signal d'interférence illustrées sur la figure 41) sont données par les relations suivantes : I,nax = (1+p)' et Iman =0 ùp)2 . En conséquence, le donné Imin (1ù p) 2 2 rapport des amplitudes minimum et maximum est par Imax (1+ p) d'où l'on déduit l'amplitude du coefficient de réflexion : p = (1- VImin/Imax)/(1 + -Jlmin/lmax) (6) Une fois la valeur de p ainsi calculée, dans une deuxième étape 52, cette valeur est comparée à un seuil donné qui peut être fixé par exemple à 30 0,5. Si la valeur de p est inférieure à ce seuil, on considère qu'il n'y a pas d'interférences dues à la réflexion du signal réfléchi par le sol. La localisation de la cible en élévation est alors réalisée au moyen d'un traitement classique d'écartométrie 53 de type monopulse par exemple. Dans le cas contraire où 25 une réflexion due au dipôle cible-image est bien détectée, un traitement adapté au signal résultant de la cible et de son image est effectué. Ce traitement commence par une étape 54 d'échantillonnage spatial du signal reçu sur chacune des lignes d'antenne tel qu'illustré par exemple sur la figure 4a. Cet échantillonnage nécessite avantageusement peu de valeurs, il peut notamment être réalisé pour une dizaine de valeurs. Dans la suite, le nombre des échantillons collectés sera représenté par n. Dans une étape suivante 55, la fréquence spatiale du signal d'interférence va être estimée à partir de l'échantillonnage du signal réalisé i o dans l'étape précédente. Le principe de base de la détermination de la fréquence est le suivant : Soit V un signal générique de la forme : V(x)=cos(2fx+ça) (7) Avec : 15 • (p, représentant une inconnue, • x , la position d'une ligne sur le réseau de l'antenne • f, la fréquence du signal reçu. Afin d'estimer la fréquence f du signal global reçu, on considère par exemple un ensemble de n échantillons représenté par la grandeur S = (si,..., sä) du 20 signal d'interférence V (x) échantillonné sur les positions xi = xo + iO avec i prenant les valeurs de 1 à n, xo représentant la position d'une ligne quelconque servant de référence sur l'antenne et A étant le pas d'échantillonnage correspondant à la distance entre deux lignes rayonnantes consécutives. Ainsi pour deux lignes quelconques on obtient à 25 partir de l'équation (7) la relation suivante quelque soit i et quelque soit k: Si+k + Si-k = cos (2 7Zfx i+k + q))+ cos (2 fx i-k + ç) (8) Or par échantillonnage xi+k = xi + kO et de même xi+k = xi ù kO . 30 La relation (8) devient alors quelque soit i et quelque soit k : Si+k + Si_k = cos (2 fxi + ça + 27 fkO) + cos (21Lt i + ça ù 27zfkO ) (9) Puis à l'aide de la relation suivante cos (X û Y )+ cos (X + Y) = 2 cos (X )cos (Y) (10) on obtient à partir de (9) la relation Si+k + Si-k = 2 cos (27zfci + ço )cos (27JkA) (11) Comme d'après l'équation (7), si = cos (2rcfxi + ço), on obtient donc la relation : si+k + si_k = 2si cos (21rfko) (12) ce qui mène à l'équation suivante : Si+k + si k û 2si cos (2e0) = 0 (13) La résolution de cette équation permet d'obtenir une estimation de f . Cette estimation est obtenue par minimisation de la fonction J définie par : J(W ) = ISi+k + si_k û 2si cos(2rrWk (14) i,k f est alors estimée par la valeur de Wo qui minimise la fonctionnelle J . Soit : min J(W) = J(Wo ) L'estimée f de f est alors donnée par : f = arg(min J(W )) = Wo w Dans le cadre du problème traité, un modèle du signal ré-émis provenant du bipôle cible-image peut être le suivant : 25 30 (15) (16) E(x) = A+acos(27fx+ço)+b(x) (17) où (ço,A,a) sont inconnues et b(x) représente un bruit blanc Gaussien de variance .2. Dans ce cas, on considère n échantillons de ce signal -rl=E(xo+1A) R= rä = E(xo + nA) Les échantillons du bruit blanc Gaussien sont supposés indépendants entre eux et de même variance 62 . Ils seront donc négligés par la suite. Une première étape consiste par exemple à éliminer le terme A constant représentant l'amplitude du signal d'interférence 41 présent au niveau de l'antenne pour se ramener à un modèle du type sinusoïde plus bruit. A cette fin un estimateur de la forme suivante peut être utilisé : 1f A = ù min ri + max ri 2 i=[1,11] i=[1,11] / qui conduit aux échantillons RBis = principe de base de la détermination de f explicité précédemment. De la 20 même manière, l'estimation de f est donc effectuée en minimisant la fonctionnelle suivante : J(U )= Iri+k +ri_k -2ricos(2%zUkA)I (19) i,k 25 La valeur de f recherchée est la valeur qui minimise la fonction J(U). Il est à noter que certains des n échantillons retenus doivent être proches du minimum et du maximum du signal afin d'obtenir une bonne estimation A de A. Dans le cas contraire, l'estimation A risque d'être (18) Les échantillons R sont exprimés en fonction de l'estimateur A de A . Ce sur lesquels est appliqué le räûA biaisée notamment si les n points de mesure couvrent moins d'une période du signal. Ceci est mis en évidence dans l'exemple qui suit. A titre d'exemple, on utilise pour chaque résultat n échantillons du signal 17 reçus sur l'antenne avec par exemple n = 11. La forme du signal présenté 5 est la suivante : E(x) = A+acos(27zfx+~p)+b(x) (17) Sur chacune des figures 6 et 7, on retrouve en abscisse 60 les n io échantillons, et en ordonnée 61 l'amplitude A du signal d'interférence présent au niveau de l'antenne. Chaque figure présente les résultats de l'estimation de A pour deux échantillons différents. La première figure 6 utilise deux échantillonnages 63 et 64 réalisés sur 1/2 période de signal de référence. Les résultats 65 et 66 obtenus respectivement pour l'échantillon 63 et 64 donnent 15 des estimations de A qui diffèrent de près de 10%. La seconde figure 7 utilise sur 7/101ème de période cieux échantillonnages 70 et 72 donnant respectivement les résultats 71 et 73 pour la valeur estimée de A. La différence entre les deux résultats est ici de l'ordre de 4%. Ceci montre bien que les résultats de l'estimation de A peuvent être 20 fortement biaisés si les échantillons couvrent seulement une petite partie de la période du signal. Ceci entraîne donc aussi une erreur importante sur l'estimation de f . Il est donc nécessaire de choisir correctement les échantillons utilisés, ce qui n'est pas toujours possible. C'est pourquoi, de manière préférentielle, le terme constant A peut 25 être éliminé par soustraction directe du signal reçu sur une ligne au signal reçu sur la ligne suivante. On obtient ainsi un nouvel échantillonnage R' de la forme suivante : R' ù r'1 = r1 ù r2 r'2 = r2 ù r3 . Les traitements d'estimation de r',_1 = rnù1 ù rn f sont ensuite les mêmes que ceux décrit précédernment. Une fois l'estimation obtenue pour f , l'étape suivante 56 permet 30 d'obtenir le rapport h2/d, hauteur de la cible sur la distance entre le radar et la cible projetée au sol, de la façon suivante : En considérant à un instant donné comme référence de phase et d'amplitude le chemin direct pour la première ligne, celle-ci reçoit un signal de valeur S1 : S1 =1+pxe"axe J(2z12 )2('1h2/d) (20) Avec : • p : le module du coefficient de réflexion sur le sol • cp : la phase du coefficient de réflexion sur le sol • : la longueur d'onde du signal émis par la cible Les hypothèses sont les suivantes : • p et cp sont constants pour les n signaux ré-émis, • Oc est identique pour les n signaux directs (on suppose que les signaux sont parallèles) En posant A comme la distance entre la première et la seconde ligne, le signal reçu par la seconde ligne du réseau s'exprime de la façon suivante : s2 =ei(2z/2)AxSln(~)(1+pxeJ0xe J(27c/2)2(112(h1+A)/d)) (21) Ce qui donne pour la (k-1)ième ligne: s =e1(27z12)kxAxsln(0)(1+pxe1(0xe i(2z1.2)2(h2(h1+kxA)/d)) 1+k (22) Et donc pour la kiéme ligne : s k e>(2TC/ü)(kù1)xAxsin(~)(1+ pxei0 xe J(2~z/Â)2(h2(hl + (k -1) x A) / d) ) k (23) Le module du signal reçu pour chaque ligne est donc donné par la relation suivante : 2 = Sk x Sk* =1 + p2 + 2p cos((27r l2)2(h2 (h1 + (k --1) x A) / d) + ça) (24) Sk En posant q)' = (2z/ À)2(h2 (h1 ù a) / d) + rp , on obtient la relation suivante : Sk 2 =1 + p 2 + 2p cos((2z / 2) 2(h2 x k x A) / d) + ço) (25) Cette relation fait apparaître une modulation d'amplitude sur la hauteur de 5 l'antenne dont la fréquence f dépend directement de la hauteur de la cible par rapport au sol : f = 2h 2/2d (26) 10 ce qui donne le rapport hauteur de la cible sur la distance projeté au sol entre le radar et la cible suivant : h2/d = Âf /2 (27) 15 Le rapport h2/d est donc défini pour une valeur de la longueur d'onde À donnée. La dernière étape 57 a notamment pour objectif de donner une estimation de la hauteur de la cible ainsi que des angles d'élévation. En considérant que la cible évolue à très basse altitude, les angles 20 d'élévation ont notamment une faible valeur. Ainsi., une estimation d' de d peut être obtenue, par exemple par une mesure radar. Cette estimation d', du fait de la faible valeur des angles, peut être assimilée à la distance radiale r entre le centre de phase de l'antenne radar et la cible. Cette distance r étant mesurée par le radar. L'erreur commise en faisant cette estimation tend 25 par ailleurs à diminuer au fur et à mesure que la cible se rapproche du sol. Les angles d'élévation décrits dans la figure 2 sont équations suivantes obtenues par construction : • pour l'angle d'élévation 9c de la cible vue du radar : 8 /h h = artct 2 ' gd d' 30 • pour l'angle d'élévation 0; de l'image vue du radar : B, = ùartctg / d + d' i Le rapport h2/d ayant été calculé à l'étape précédente, les angles 0c et 0; se calculent facilement. Par ailleurs, une fois 0c déterminé, h2 s'en déduit de la manière suivante : h2=rsin(0c)+hi (30) 5 r étant la distance radiale mesurée par l'antenne radar. Les équations (28), (29), (30) permettent donc d'obtenir les valeurs estimées de Oc, 0;, et h2. Les figures 8a et 8b mettent en évidence la différence entre les 10 résultats obtenus dans le cas où l'on a un sol horizontal et dans le cas où il ne l'est pas. La figure 8a représente le cas d'un sol horizontal et notamment les grandeurs caractéristiques de l'élévation : Oc, 0;, et h2. Le principe de construction de la réflexion quasi-spéculaire du rayon ré-émis par la cible 26 est que les angles d'incidence du signal 26 et de réflexion du signal 27 par 15 rapport au plan de réflexion sont supposés égaux à un angle a,. La hauteur h2, qui intervient dans les équations précédentes, correspond en réalité à la hauteur de la cible par rapport au plan de réflexion, qui ici est confondu avec le sol, et non pas à la hauteur de la cible par rapport au sol à sa verticale. Si l'on considère la figure 8b, lorsque le sol n'est pas horizontal, l'image de la 20 cible n'a pas la même position que sur la figure 8a. En effet, l'image est la projection symétrique de la cible par rapport au plan de réflexion qui, dans le cas de la figure 8b, n'est pas horizontal. La hauteur h2 de la cible par rapport au plan de réflexion sera donc différente dans le cas d'un sol non horizontal. L'estimation de l'élévation peut ainsi être biaisée et il peut être nécessaire 25 d'affiner le résultat. La figure 9 représente le phénomène de réflexion en présence d'un sol 90 non plan et non horizontal, le profil du sol entre le radar et la cible étant connu. Les éléments suivants de la figure 1 et 2 sont représentés sur la 30 figure : • la cible 1 représentée de manière ponctuelle, • la ligne 23 de l'antenne située à une hauteur hi du sol, • la distance radiale r entre le radar et la cible, • la distance d horizontale entre le radar et la cible, • le signal 24 émis par la ligne 23 et ré-émis en un signal 25 de manière directe par la cible 1, • le signal 26 ré-émis en direction du sol par la cible et sa réflexion 27 sur le sol. La hauteur H2 représente la hauteur de la cible par rapport au plan 91 tangent au sol en un point R appelé point de réflexion. La détermination du point de réflexion passe par le calcul de sa position dans le référentiel du radar. La position de R est caractérisée par ses coordonnées XR et ZR respectivement selon l'axe horizontal 93 et l'axe vertical 94 ainsi que par l'angle aR formé entre la pente 91 du sol en ce point et l'horizontale 92. Le profil du sol est supposé connu, ce qui donne accès à la connaissance de zR et aR pour toutes les positions potentielles du point R c'est à dire que zR=f(xR), la fonction f étant connue, et aR=g(xR), la fonction g étant également connue. Le principe utilisé pour calculer les coordonnés du point R est basé sur le fait que les angles de réflexion et d'incidence du signal ré-émis par la cible en direction du sol sont égaux en ce point. Ils sont représentés respectivement sur la figure par les angles a1 et a2. En considérant que les angles sont faibles, l'hypothèse suivante est prise: tan(a) = sin(a) = a, a étant l'un quelconque des angles représentés sur la figure. Par construction et en prenant l'hypothèse énoncée ci-avant on obtient les relations suivantes pour a1 et a2 : h, - zR + xR xaR (31) xR et a2 = h2 (32) d - x,Z Comme a1 = a2, on obtient alors la relation suivante : h, -zR +xR xaR h2 xR d - xR avec h1 la hauteur de l'antenne par rapport au sol et h2, la hauteur d'élévation calculée dans le cas d'un sol plan et horizontal. a, = (33) Les coordonnées du point R sont obtenues par minimisation de la fonction de coût suivante : h,-zR+x12xaR h, (34) xR r ù xR Les coordonnées du point R ainsi obtenues permettent de définir le plan tangent au sol au point de réflexion, ou plan de réflexion, qui est par la suite assimilé au sol lui-même afin deramener le cas du sol non plan et non 1 o horizontal au cas du sol non horizontal plus simple à traiter. Une fois les coordonnées du point R calculées et notamment l'angle OCR, une correction est appliquée sur la hauteur d'élévation obtenue pour un sol plan et horizontal. La hauteur h2 par rapport au sol horizontal devient alors la hauteur H2 par rapport au plan de réflexion. Cette hauteur H2 est donnée par 15 la relation suivante : H2 = h2 + zR + (d ù xR) x aR (36) où d représente la distance entre la cible et le radar projetée sur l'horizontale 20 92. On obtient ainsi la valeur corrigée pour la hauteur de la cible par rapport au plan de réflexion pour un sol non plan et non horizontal. L'estimation ainsi obtenue pour l'élévation peut être avantageusement améliorée par un calcul de corrélation en amplitude phase du signal reçu 25 avec une réplique prédéfinie aussi nommé CHA ou Correlation Height Algorithm. Un des avantages du procédé selon l'invention est d'améliorer la localisation en élévation tout en conservant une dimension d'antenne constante, et une possibilité de veille dans un domaine angulaire étendu. La 30 stabilité des résultats obtenus dans le domaine de couverture sans élargissement de lobe et sans affaiblissement est également appréciable. D'autre part, la simplicité du traitement mis en oeuvre permet son application à moindre coût. II nécessite également une capacité de calcul raisonnable étant donné que l'on peut limiter le nombre de voix de réception. En effet un petit nombre de mesures suffit pour obtenir une bonne estimation de l'élévation. Un autre avantage de l'invention est que le procédé peut être utilisé aussi bien sur une antenne réseau complète que sur une antenne lacunaire. De plus elle peut aussi bien s'appliquer pour une antenne comportant des réseaux uniformes ou non uniformes. Dans l'exemple de mise en oeuvre de l'invention, l'antenne de réception des signaux ré-émis par la cible correspond à l'antenne du radar émettant le signal de détection. Cette antenne pourrait être différente de l'antenne d'émission. REVENDICATIONS 1. Procédé de détermination de l'élévation d'une cible (1) proche du sol par un système de détection électromagnétique (3) caractérisé en ce qu'il extrait la mesure de l'élévation à partir de l'amplitude du signal d'interférence (41) produit par un signal (25) émis directement par la cible et par un signal (26) émis par la cible en direction du sol puis réfléchi par le sol vers le radar. 2. Procédé selon la revendication 1, caractérisé en ce que le système émet un signal de détection vers la cible, le signal d'interférence étant produit par le signal (25) ré-émis directement par la cible et le signal (26) ré-émis par la cible en direction du sol puis réfléchi par le sol vers le système de détection. 3. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'antenne de réception étant constituée de réseaux d'éléments rayonnants, le signal d'interférence (41) reçu par cette antenne est échantillonné au niveau des sous-réseaux (22), la mesure de l'élévation étant établie à partir de l'amplitude du signal d'interférence échantillonné. 4. Procédé selon la revendication 3, caractérisé en ce que les sous-réseaux (22) correspondant chacun à une ligne d'éléments rayonnants sont disposés de manière horizontale, parallèle et uniforme sur l'ensemble de l'antenne. 5. Procédé selon l'une quelconque des revendications 3 ou 4, caractérisé en ce que les sous-réseaux se situent dans un plan sensiblement vertical. 6. Procédé selon l'une quelconque des revendications 3 à 5, caractérisé en ce que la fréquence f du signal échantillonné, est la valeur qui minimise la 30 fonction suivante : J (U ) = ~j I Ni+k + ri-k ù tri cos (22rUkA i,k avec• r;+k représentant le résultat de la soustraction entre la valeur du signal reçu au niveau du sous-réseau d'ordre i+k et la valeur du signal reçu au niveau du sous-réseau d'ordre i+k+1, r;_k représentant le résultat de la soustraction entre la valeur du signal reçu au niveau du sous- réseau d'ordre i-k et la valeur du signal reçu au niveau du sous-réseau i-k-1 • A représentant la distance entre les centres de phases de deux sous-réseaux consécutifs. 7. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le rapport entre la hauteur h2 de la cible par rapport au sol et la distance d projetée au sol entre le radar et la cible s'exprime en fonction de la fréquence f du signal d'interférence calculée selon la relation suivante : h2/d = Àf / 2 où est la longueur d'onde du signal émis. 8. Procédé selon la revendication 7, caractérisé en ce que la hauteur d'élévation h2 est déterminée pour un sol non plan et non horizontal à partir : • du point de réflexion R du signal ré-émis par la cible vers le sol (26) sur ce dernier, • de la hauteur d'élévation déterminée pour un sol plan et horizontal, L'expression de l'égalité entre les angles d'incidence du rayon ré-émis par la cible vers le sol (26) et celui réfléchi par le sol avec le plan tangent au sol au point de réflexion (27) permettant d'obtenir les coordonnées (xR, ZR) du point de réflexion R ainsi que l'angle îR entre le plan tangent au sol (91) au point R et l'horizontale (80) à partir de la fonction de coût suivante : h, - zR + xR x aR h2 où h, représente la hauteur par rapport au sol de l'antenne, h2 représente l'estimation de la hauteur d'élévation calculée dans le cas d'un sol plan et horizontal, et les valeurs (xR, zR, îR) minimisant cette fonction sont les coordonnées du point de réflexion R. 9. Procédé selon la revendication 8, caractérisé en ce que l'obtention de la valeur de l'angle OCR entre le plan tangent en R au terrain et l'horizontale permet de corriger la hauteur de la cible par rapport au plan tangent en une valeur H2 de la manière suivante : Hz = h2 + zR -f- (d û xR) x aR avec h2 la hauteur calculée de la cible par rapport au sol dans le cas d'un sol plan et horizontal et d étant la distance projetée au sol de la cible à l'antenne. 10. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le traitement selon l'invention est amélioré par un 1 o traitement CHA. 11. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'antenne utilisée est une antenne lacunaire. 15 12. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la cible comporte un émetteur actif. 13. Procédé selon l'une quelconque des revendications 1 à 14, caractérisé en ce que la cible est un réflecteur passif. 14. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la cible est un drône en phase d'atterrissage. 15. Procédé selon les revendications 1 ou 2, caractérisé en ce que la 25 comparaison de l'amplitude du coefficient de réflexion avec un seuil fixé est préalablement effectuée, la hauteur h2 de la cible étant déterminée par un traitement de type monopulse si la valeur de l'amplitude du coefficient de réflexion est inférieure à ce seuil. 20