FR 2759784 A1 19980821 FR 9702007 A 19970220 La présente invention concerne un procédé d'évaluation du déphasage perturbateur des signaux radar d'un radar à synthèse d'ouverture et un dispositif pour la mise en oeuvre de ce procédé, ainsi que des applications dudit procédé. On sait que l'observation de la terre à l'aide de radars à synthèse d'ouverture embarqués sur un aéronef est la source d'un nombre élevé d'applications, tant dans le domaine civil (agriculture, géologie, océanographie,...) que dans le domaine militaire (observation d'un champ de bataille, recalage de la navigation d'un missile,...). Un tel radar à synthèse d'ouverture présente de très nom- breux avantages, notamment pour réaliser une cartographie d'un terrain survolé, qui expliquent le grand intérêt et le développement important de ce type de radar, en parti- culier: - il peut être utilisé par tout temps; - il permet d'élaborer des images de très haute résolution; - il présente un encombrement réduit; et - ses traitements sont peu dépendants de la distance d'ob- servation. De façon connue, un tel radar à synthèse d'ouverture compor- te des moyens pour réaliser successivement les opérations suivantes: émettre un rayonnement électromagnétique en direction d'un territoire observé et survolé par l'aéronef porteur dudit radar; - détecter un signal correspondant audit rayonnement élec- tromagnétique renvoyé par ledit territoire survolé et observé; - effectuer un traitement de détermination de distance, à partir dudit signal détecté; - effectuer un traitement de focalisation destiné à corriger un déphasage perturbateur existant dans ledit signal détecté et dû notamment au déplacement dudit radar à synthèse d'ouverture de manière à obtenir un signal corrigé; - effectuer un traitement de détermination de fréquence à partir dudit signal corrigé; et - effectuer un traitement de transposition destiné à trans- poser, dans un repère de référence lié audit territoire observé, les résultats obtenus à partir des opérations précédentes et définis dans un repère de distance et de fréquence lié au radar, de manière à former une image radar définie dans ledit repère de référence lié au territoire observe. Toutefois, un tel radar à synthèse d'ouverture présente l'inconvénient que certaines des opérations citées précédem- ment, notamment le traitement de focalisation, nécessitent pour leur mise en oeuvre des données de navigation très précises, relatives à l'aéronef porteur dudit radar. Ces données de navigation sont généralement déterminées par des équipements auxiliaires de navigation embarqués sur l'aéronef. Aussi, pour être performante, la mise en oeuvre dudit radar à synthèse d'ouverture impose de lourdes contraintes sur la précision desdits équipements auxiliai- res, ce qui entraîne des conséquences négatives quant au coût, à la masse et au volume de ces derniers. Pour remédier à cet inconvénient, on connaît des moyens qui permettent de compenser les erreurs dues aux mouvements de l'aéronef porteur du radar, et qui sont notamment suscepti- bles de déterminer ledit déphasage perturbateur, sans utiliser lesdits équipements auxiliaires, en analysant simplement le signal radar détecté. Ces moyens connus mettent en oeuvre des solutions complexes de traitement de signaux qui permettent de déterminer un gradient Doppler représentatif dudit déphasage perturbateur. Une première solution connue, divulguée par la demande de brevet français n 95 14555 déposée le 8 décembre 1995 par la demanderesse, propose de mesurer le décalage fréquentiel moyen entre deux images radar et de déterminer directement ledit gradient Doppler à partir de ce décalage fréquentiel moyen. Cette première solution, qui permet un traitement rapide, est toutefois peu robuste. Une seconde solution connue consiste à engendrer une plura- lité d'images radar et à sélectionner l'image optimale en ce qui concerne le contraste, ce qui permet d'en déduire le gradient Doppler. Cette seconde solution, bien que simple et robuste, néces- site une charge de calcul élevée. Aussi, pour diminuer un peu cette charge de calcul, on effectue généralement le calcul de façon itérative, au détriment de la robustesse de la solution (apparition de problèmes d'initialisation et de convergence vers des maxima locaux du contraste, maîtrise réduite de la charge de calcul,...). La présente invention a pour objet de remédier à ces in- convénients. Elle concerne un procédé de mise en oeuvre rapide et robuste, permettant d'évaluer le déphasage pertur- bateur des signaux radar d'un radar à synthèse d'ouverture, sans utiliser d'informations relatives au mouvement de l'aéronef porteur dudit radar. A cet effet, selon l'invention, ledit procédé est remarqua- ble en ce que:30 - on effectue une pluralité d'opérations complexes dans chacune desquelles: A/ on mesure N signaux radar successifs, N étant un nombre entier prédéfini, lesdits signaux radar étant des signaux complexes définis dans le domaine temporel; B/ on répète i fois les opérations suivantes, i variant de 1 à N/(2K), K étant un entier prédéfini inférieur ou égal à N/2 et de préférence supérieur ou égal à 10: a) à partir des N signaux radar successifs mesurés, on calcule le produit conjugué d'un signal radar initial et d'un signal radar décalé de it, r étant une variable temporelle; B) on calcule une expression mathématique correspondant au double du carré de la transformée de Fourier dudit produit conjugué; y) on réalise une transposition de manière à faire dépendre cette expression mathématique d'un gradient Doppler; C/ on calcule la racine carrée de la somme des résultats obtenus à l'étape y) pour les N/(2K) répétitions, de manière à obtenir l'écart-type du spectre de puissance des signaux radar; D/ à partir des signaux radar mesurés à l'étape A/, on détermine la puissance moyenne desdits signaux; E/ à partir dudit écart-type du spectre de puissance o et de ladite puissance moyenne m, on calcule le contraste X à partir de la relation X=(a/m); et - à partir de la pluralité des valeurs de contraste ainsi calculées, on détermine le maximum dudit contraste, de manière à mettre en évidence un gradient Doppler optimal qui est représentatif dudit déphasage perturbateur. Ainsi, grâce à l'invention, on évalue ledit déphasage perturbateur directement à partir des signaux radar détec- tés, donc sans utiliser d'informations fournies par des équipements auxiliaires de navigation. De plus, ledit procédé qui préconise une optimisation du contraste permet une mise en oeuvre rapide et robuste. En outre, au contraire de la seconde solution connue précitée préconisant également une optimisation du contraste, le procédé conforme à l'invention s'affranchit de toute approche itérative, ce qui est bien entendu avantageux pour la vitesse de mise en oeuvre et la robustesse (pas de problème d'initialisation ou de convergence vers des maxima locaux du contraste, parfaite maîtrise de la charge de calcul). 1O De plus, de façon avantageuse: - à l'étape y), pour réaliser la transposition, on remplace dans ladite expression mathématique la fréquence f par l'expression (M.i.r), B étant ledit gradient Doppler; et/ou - à l'étape D/, la puissance moyenne m est calculée à partir de l'expression: j=N m = Si 1 2 j=1 les signaux sj étant les signaux radar mesurés à l'étape A/. En outre, avantageusement, le maximum du contraste est déterminé par moyenne de ladite pluralité des valeurs de contraste calculées. La présente invention concerne également un dispositif pour la mise en oeuvre du procédé précité. Selon l'invention, ledit dispositif est remarquable en ce qu'il comporte: - des premiers moyens pour calculer l'écart-type du spectre de puissance des signaux radar; - des deuxièmes moyens pour calculer la puissance moyenne desdits signaux radar; et - des troisièmes moyens pour calculer, à partir des résul- tats obtenus par lesdits premiers et deuxièmes moyens, le contraste d'une image radar, pour déterminer le maximum dudit contraste et pour en déduire ledit gradient Doppler optimal. Par ailleurs, la présente invention concerne deux applica- tions du procédé précité. Une première application est relative à un radar à synthèse d'ouverture embarqué à bord d'un aéronef et comportant de façon connue: - des moyens pour émettre un rayonnement électromagnétique en direction d'un territoire observé et survolé par ledit aéronef; - des moyens pour détecter un signal correspondant audit rayonnement électromagnétique renvoyé par ledit territoire observé et survolé; - des moyens pour effectuer un traitement de détermination de distance, à partir dudit signal détecté; - des moyens de focalisation destinés à corriger un dépha- sage perturbateur existant dans ledit signal détecté et dû notamment au déplacement dudit radar à synthèse d'ouver- ture de manière à obtenir un signal corrigé; - des moyens pour effectuer un traitement de détermination de fréquence, à partir dudit signal corrigé; et - des moyens de transposition destinés à transposer, dans un repère de référence lié audit territoire observé, les résultats obtenus à partir des opérations précédentes et définis dans un repère de distance et de fréquence lié au radar, de manière à former une image radar définie dans ledit repère de référence lié au territoire observe. Selon l'invention, ledit radar est remarquable en ce que: - dans un premier mode de réalisation, lesdits moyens de focalisation utilisent des informations obtenues par la mise en oeuvre du procédé précité conforme à l'invention, par exemple pour affiner la correction du déphasage perturbateur réalisée de façon usuelle, à l'aide d'équipe- ments auxiliaires de navigation; et - dans un second mode de réalisation, pour déterminer ledit déphasage perturbateur, lesdits moyens de focalisation mettent en oeuvre exclusivement le procédé conforme à l'invention. Dans ce second mode de réalisation, de préférence, les moyens de transposition utilisent également des informations obtenues par la mise en oeuvre dudit procédé, ce qui permet de réaliser un radar à synthèse d'ouverture complètement autonome vis-à-vis des équipements de navigation. Une seconde application du procédé conforme à l'invention concerne un dispositif pour déterminer des informations cinématiques d'un aéronef porteur d'un radar à synthèse d'ouverture, ledit dispositif comportant, de façon connue: - des premiers moyens de calcul pour déterminer le gradient Doppler du déphasage perturbateur des signaux radar détectés par ledit radar à synthèse d'ouverture; et - des seconds moyens de calcul pour déduire dudit gradient Doppler lesdites informations cinématiques, par exemple la vitesse et/ou l'accélération radiale de l'aéronef. Dans ce cas, selon l'invention, lesdits premiers moyens de calcul mettent en oeuvre le procédé précité. Les figures du dessin annexé feront bien comprendre comment l'invention peut être réalisée. Sur ces figures, des réfé- rences identiques désignent des éléments semblables. La figure 1 est le schéma synoptique d'un dispositif confor- me à l'invention. La figure 2 est le schéma synoptique d'une première applica- tion du dispositif conforme à l'invention. La figure 3 est le schéma synoptique d'une seconde applica- tion du dispositif conforme à l'invention. Le dispositif 1 conforme à l'invention et représenté schéma- tiquement sur la figure 1 est destiné à évaluer le déphasage perturbateur des signaux radar d'un radar à synthèse d'ou- verture, non représenté sur cette figure 1. De façon connue, le signal radar s(t) détecté par un radar à synthèse d'ouverture peut s'écrire sous la forme: s(t) = A(t).exp[jP(t)]. exp(japot2) dans laquelle: - A(t) représente une amplitude; - ó(t) est la phase; et - l'expression exp(jupot2) représente le déphasage perturba- teur dû notamment au déplacement dudit radar. Par conséquent, pour connaître ledit déphasage perturbateur, il suffit de déterminer le gradient Doppler Mo. Pour ce faire, selon l'invention, ledit dispositif 1 compor- te: - des moyens 2 pour calculer l'écart-type du spectre de puissance des signaux radar précisés ci-dessous et reçus par une branche 3A d'une liaison 3 dédoublée; - des moyens 4 pour calculer la puissance moyenne desdits signaux radar reçus par une branche 3B de la liaison 3; et - des moyens 5 reliés aux moyens 2 et 4 respectivement par l'intermédiaire de liaisons 6 et 7, pour calculer, à partir des résultats obtenus par lesdits moyens 2 et 4, le contraste d'une image radar, pour déterminer le maximum dudit contraste et pour en déduire ledit gradient Doppler po. Les moyens 5 peuvent transmettre le résultat obtenu, par l'intermédiaire d'une liaison 8, à un dispositif utilisateur non représenté sur la figure 1. Pour la mise en oeuvre du dispositif 1, ce dernier reçoit, pour chacune d'une pluralité d'opérations complexes, N signaux radar du radar à synthèse d'ouverture, lesdits signaux étant des signaux complexes définis dans le domaine temporel et précisés ci-dessous, et N étant un entier prédéfini. Ces signaux sont transmis auxdits moyens 2 et 4. o10 Selon l'invention, lesdits moyens 2 comprennent notamment des éléments de calcul 9 à 12 qui réalisent i fois les mêmes opérations à partir de signaux radar différents, comme indiqué cidessous, i variant de 1 à N/(2K) et K étant un entier prédéfini, inférieur ou égal à N/2 et de préférence15 supérieur ou égal à 10. On décrit ci-après les traitements réalisés par lesdits éléments de calcul 9 à 12 reliés ensemble respectivement par l'intermédiaire de liaisons 13 à 15, pour une opération i quelconque:20 - l'élément de calcul 9 relié à la liaison 3A calcule le produit conjugué du signal radar initial s(t) et d'un signal radar s(t-ir) décalé de ir, t étant une variable temporelle; l'élément de calcul 10 calcule la transformée de Fourier dudit produit conjugué; - l'élément de calcul 11 calcule une expression mathématique correspondant au double du carré de ladite transformée de Fourier; et l'élément de calcul 12 réalise une transposition de ladite expression mathématique qui dépend de la fréquence f de manière à la faire dépendre du gradient Doppler p véri- fiant la relation p=f/(i.r). Les résultats des calculs de l'élément de calcul 12 sont transmis par une liaison 16 à un élément de calcul 17 qui calcule la racine carrée de la somme des différents résul- tats obtenus par les i opérations successives de manière à former l'écart-type du spectre de puissance. Par conséquent, l'écart-type a calculé par lesdits moyens 2 vérifie l'expression: /) = I2i=N/(2K) * a(p) = V 2 Z IF[s(t).s (t-iz](f=ir) 2 (1) s (t-ir) étant le conjugué de s(t-ir). On précise à présent la justification mathématique de cette expression (1). De façon connue, l'écart-type o du spectre de puissance est défini par la relation: (H) = _J (Pfoc(f))2df - m2 dans laquelle: - Pfoc(f) est le spectre de puissance; - f est la fréquence; - m est la puissance moyenne; et - p est le gradient Doppler. On notera que dans les modes de calcul usuels, d'optimisa- tion du contraste, l'écart-type o qui dépend du gradient y doit être calculé pour toutes les valeurs de p existantes, ce qui n'est pas le cas dans la présente invention, comme on le précisera ci- dessous. A l'aide du théorème de Parceval, on obtient la relation: +CO+ o(p) j= V_ Pfoc(f)J2df - m2 = _ CIF l[Pfoc(f)]12dr - m2 -1 F étant la transformée de Fourier inverse et r la variable temporelle. Or, la transformée de Fourier inverse du spectre de puis- sance est égale à la fonction d'autocorrélation du signal temporel (théorème de Wiener-Kintchine) qui vérifie donc: -1 Rfoc(t) F [Pfoc(f)] = - J sfoc(t)s foc(t-z) dt, sfoc(t) étant le signal radar complexe dans le domaine temporel, après le traitement de focalisation et s foc(t-t) étant le conjugué de sfoc(t-r), ce qui permet de déduire l'expression suivante de l'écart-type o: a(p) ( = __ JRfoc(r)]2dr - m2 A partir des propriétés suivantes de la fonction d'autocor- rélation: Rfoc(0) = m puissance moyenne Rfoc(-r) = R foc(t) propriété de symétrie, R foc(r) étant le conjugué de Rfoc(r), on obtient l'expres- sion simplifiée: o(P) = V 2 >0 IRfoc(t)]12 dr On calcule à présent la valeur de Rfoc(r) à partir d'une analyse fréquentielle du signal s(t) non focalisé. On sait que le traitement de focalisation du signal s(t) consiste à effectuer l'opération suivante: sfoc(t) = s(t) exp(-jrpt2). En intégrant cette expression dans la fonction d'autocorré- lation Rfoc(r), on obtient la relation: +a, Rfoc(r) = -J sfoc(t)s foc(tr) dt = exp(jurpt2) -f s(t)s (t-r)exp(-j2nptt) dt = exp(jrrur2)F[s(t).s (t-t)]f= F étant la transformée de Fourier. On en déduit finalement l'expression suivante (2) de l'écart- type: (10 o() =/2 t>0 I F[s(t).s (t-r)]f =12 dr (2) De plus, par la transposition de cette expression (2) dans le domaine numérique, on obtient l'expression (1) précitée utilisée pour la mise en oeuvre de la présente invention. Par conséquent, au contraire des méthodes connues de calcul d'optimisation du contraste qui nécessitent N calculs pour l'écart- type à partir des N signaux radar détectés, le procédé conforme à l'invention est mis en oeuvre à partir uniquement de N/(2K) calculs de l'écart-type, ce qui en- traîne une diminution importante de la charge de calcul, par exemple d'un facteur 20 lorsque K est égal à 10 et d'un20 facteur 200 lorsque K est égal à 100. Par ailleurs, les moyens 4 calculent la puissance moyenne m à partir de l'expression: j=N m= E Isju2 (3) j=1 les signaux sj étant les signaux radar reçus par la liaison 3B. Pour justifier l'expression (3), on notera que, de façon connue, la puissance moyenne s'écrit: m = -J Pfoc(f)df = J Sfoc(f)12 df, Pfoc(f) et Sfoc(f) correspondant respectivement au spectre de puissance et au signal complexe, dans le domaine fréquen- tiel, après traitement de focalisation. En utilisant le théorème de Parceval, on obtient: r+C m = J Isfoc(t) I2 dt r+C - J Is(t)exp(-jrrpt2) I2 dt, s(t) et sfoc(t) correspondant respectivement au signal radar complexe dans le domaine temporel, avant et après traitement15 de focalisation. On peut en déduire l'expression de la puissance moyenne à partir du signal s(t) non focalisé: m= J Is(t) 2 dt (4), indépendant de M. La puissance moyenne m peut donc être calculée une fois pour toute, uniquement à partir du signal radar complexe non focalisé s(t). Finalement, en transposant cette expression (4) dans le domaine numérique, on obtient l'expression (3) utilisée par les moyens 4. Par ailleurs, les moyens 5 comportent: - un élément de calcul 18 qui calcule le contraste X à partir des valeurs a et m reçues respectivement par les liaisons 6 et 7, au moyen de la relation X=ao/m; - un élément de calcul 19 qui est relié par une liaison 20 à l'élément de calcul 18 et qui calcule un contraste moyen xm à l'aide de la relation Xm=x/N; et - un élément de calcul 21 qui est relié par une liaison 22 à l'élément de calcul 19 et qui détermine, de façon connue, le maximum dudit contraste par moyenne des contrastes moyens de la pluralité desdites opérations complexes, de manière à mettre en évidence le gradient Doppler optimal Mo recherché. Ce dernier peut être transmis à un dispositif utilisateur par la liaison 8. Ainsi, grâce à l'invention, on évalue ledit gradient Doppler Mo représentatif du déphasage perturbateur directement à partir des signaux radar détectés, donc sans utiliser d'informations fournies par des équipements auxiliaires de20 navigation. De plus, le dispositif 1 conforme à l'invention est rapide et robuste. Ledit dispositif 1 peut être mis en oeuvre dans de nombreu- ses applications. Une première application représentée schématiquement sur la figure 2 concerne un radar D1 à synthèse d'ouverture. Ce radar D1 est monté sur un aéronef non représenté et est utilisé pour observer un territoire survolé par ledit aéronef. De façon connue, un tel radar D1 à synthèse d'ouverture comporte: - des moyens 23 pour émettre un rayonnement électromagnéti- que S1 en direction d'un territoire observé et survolé par ledit aéronef; - des moyens 23 pour détecter un signal S2 correspondant audit rayonnement électromagnétique renvoyé par ledit territoire. Lesdits moyens 23 d'émission et de réception sont réalisés, de préférence, sous forme d'un système électronique 23A comportant une antenne 23B d'émission et de réception; - des moyens 24 reliés par une liaison 25 auxdits moyens 23 et susceptibles d'effectuer un traitement de détermination de distance, à partir dudit signal détecté S2 de manière à obtenir le signal radar s(t) indiqué précédemment; - des moyens 26 reliés par une liaison 27 auxdits moyens 24 et susceptibles d'effectuer un traitement de focalisation destiné à corriger le déphasage perturbateur existant dans ledit signal s(t) et dû notamment au déplacement du radar D1 à synthèse d'ouverture embarqué à bord dudit aéronef, de manière à obtenir un signal corrigé; des moyens 28 reliés par une liaison 29 auxdits moyens 26 et susceptibles d'effectuer un traitement de détermination de fréquence, à partir dudit signal corrigé; et - des moyens 30 reliés par une liaison 31 auxdits moyens 28 et susceptibles d'effectuer un traitement de transposition destiné à transposer, dans un repère de référence lié audit territoire observé, les résultats obtenus à partir des opérations précédentes et définis dans un repère de distance et de fréquence lié au radar, de manière à former une image radar définie dans ledit repère de référence lié au territoire observé. Selon l'invention, le dispositif 1 est associé audit radar D1 ou intégré dans ce dernier. Plus précisément, la liaison35 3 d'entrée du dispositif 1 est liée à la liaison 27 et la liaison 8 de sortie est notamment raccordée aux moyens 26 réalisant le traitement de focalisation. Cette application peut être mise en oeuvre selon deux modes de réalisation différents. Dans un premier mode de réalisation, lesdits moyens 26 utilisent à la fois les informations transmises par le dispositif 1 et des informations transmises par des équipe- ments auxiliaires de navigation, non représentés. Dans ce cas, le dispositif 1 peut être utilisé pour affiner la correction du déphasage perturbateur, réalisée de façon usuelle, ce qui permet de diminuer les contraintes imposées auxdits équipements auxiliaires. Dans un second mode de réalisation, lesdits moyens 26 réalisent la correction du déphasage perturbateur uniquement à partir des informations transmises par le dispositif 1 conforme à l'invention. On sait en outre que les moyens 30 nécessitent également pour leur traitement des données de navigation fournies par des équipements auxiliaires. Aussi, selon l'invention, pour que le radar D1 soit complè- tement indépendant de tels équipements auxiliaires, un perfectionnement dudit second mode de réalisation prévoit de transmettre également auxdits moyens 30 les informations déterminées par le dispositif 1. Pour le traitement mis en25 oeuvre dans ce cas, on peut se référer à la demande de brevet français précitée n 95 14555. Ainsi, on obtient un radar D1 à synthèse d'ouverture entiè- rement autonome, vis-à-vis des équipements de navigation. Dans une seconde application représentée schématiquement sur la figure 3, le dispositif 1 conforme à l'invention est associé à un dispositif D2 susceptible de déterminer des informations cinématiques d'un aéronef porteur d'un radar à synthèse d'ouverture, dont on a uniquement représenté les moyens 23 et 24 nécessaires pour la mise en oeuvre de la présente application. De façon connue, un tel dispositif D2 comporte: - des premiers moyens de calcul pour déterminer le gradient Doppler du déphasage perturbateur des signaux radar détectés par ledit radar à synthèse d'ouverture; et - des seconds moyens de calcul 34 pour déduire dudit gra- dient Doppler, de façon usuelle, lesdites informations cinématiques, en particulier la vitesse de l'aéronef et/ou son accélération radiale. Selon l'invention, lesdits premiers moyens de calcul sont réalisés sous forme du dispositif 1 conforme à l'invention. On notera que ledit dispositif D2 peut transmettre les informations ainsi déterminées, par l'intermédiaire d'une20 liaison 35, à un dispositif utilisateur non représenté. Ces informations peuvent notamment être utilisées pour le recalage de la navigation, par exemple pour recaler une centrale inertielle embarquée sur ledit aéronef. REVENDICATIONS 1. Procédé d'évaluation du déphasage perturbateur des signaux radar d'un radar à synthèse d'ouverture, caractérisé en ce que: - on effectue une pluralité d'opérations complexes dans chacune desquelles: A/ on mesure N signaux radar successifs, N étant un nombre entier prédéfini, lesdits signaux radar étant des signaux complexes définis dans le domaine temporel; B/ on répète i fois les opérations suivantes, i variant de 1 à N/(2K), K étant un entier prédéfini inférieur ou égal à N/2: a) à partir des N signaux radar successifs mesurés, on calcule le produit conjugué d'un signal radar initial et d'un signal radar décalé de it, t étant une variable temporelle; B) on calcule une expression mathématique correspondant au double du carré de la transformée de Fourier dudit produit conjugué; y) on réalise une transposition de manière à faire dépendre cette expression mathématique d'un gradient Doppler; C/ on calcule la racine carrée de la somme des résultats obtenus à l'étape y) pour les N/(2K) répétitions, de manière à obtenir l'écart-type du spectre de puissance des signaux radar; D/ à partir des signaux radar mesurés à l'étape A/, on détermine la puissance moyenne desdits signaux; E/ à partir dudit écart-type du spectre de puissance o et de ladite puissance moyenne m, on calcule le contrastex à partir de la relation X=(o/m); et - à partir de la pluralité des valeurs de contraste ainsi calculées, on détermine le maximum dudit contraste, de manière à mettre en évidence un gradient Doppler optimal qui est représentatif dudit déphasage perturbateur. 2. Procédé selon la revendication 1, caractérisé en ce que ledit entier K est supérieur ou égal à 10. 3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce qu'à l'étape y), pour réaliser la transpo- sition, on remplace dans ladite expression mathématique la fréquence f par l'expression (p.i.t), p étant ledit gradient Doppler. 4. Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce qu'à l'étape D/, la puissance moyenne m est calculée à partir de l'expression: j-N m = Isj]2 J-1 les signaux sj étant les signaux radar mesurés à l'étape A/. 5. Procédé selon l'une quelconque des revendications précé- dentes, caractérisé en ce que le maximum du contraste est déterminé par moyenne de ladite pluralité des valeurs de contraste calculées. 6. Dispositif pour la mise en oeuvre du procédé spécifié sous l'une quelconque des revendications 1 à 5, caractérisé en ce qu'il comporte: - des premiers moyens (2) pour calculer l'écart-type du spectre de puissance des signaux radar; - des deuxièmes moyens (4) pour calculer la puissance moyenne desdits signaux radar; et des troisièmes moyens (5) pour calculer, à partir des résultats obtenus par lesdits premiers et deuxièmes moyens (2, 4), le contraste d'une image radar, pour déterminer le maximum dudit contraste et pour en déduire ledit gradient Doppler optimal. 7. Radar à synthèse d'ouverture embarqué à bord d'un aéronef et comportant: - des moyens (23) pour émettre un rayonnement électromagné- tique (Sl) en direction d'un territoire observé et survolé par ledit aéronef; - des moyens (23) pour détecter un signal (S2) correspondant audit rayonnement électromagnétique renvoyé par ledit territoire observé et survolé; - des moyens (24) pour effectuer un traitement de détermina- tion de distance, à partir dudit signal détecté (S2); - des moyens de focalisation (26) destinés à corriger un déphasage perturbateur existant dans ledit signal détecté et dû notamment au déplacement dudit radar à synthèse d'ouverture de manière à obtenir un signal corrigé; - des moyens (28) pour effectuer un traitement de détermina- tion de fréquence, à partir dudit signal corrigé; et - des moyens de transposition (30) destinés à transposer, dans un repère de référence lié audit territoire observé, les résultats obtenus à partir des opérations précédentes et définis dans un repère de distance et de fréquence lié au radar, de manière à former une image radar définie dans ledit repère de référence lié au territoire observé, caractérisé en ce que lesdits moyens de focalisation (26) utilisent des informations obtenues par la mise en oeuvre du procédé spécifié sous l'une quelconque des revendications 1 à 5. 8. Radar à synthèse d'ouverture selon la revendication 7, caractérisé en ce que, pour déterminer ledit déphasage perturbateur, lesdits moyens de focalisation (26) mettent en oeuvre exclusivement le procédé spécifié sous l'une quel- conque des revendications 1 à 5. 9. Radar à synthèse d'ouverture selon la revendication 8, caractérisé en ce que lesdits moyens de transposition (30) utilisent des informations obtenues par la mise en oeuvre du procédé spécifié sous l'une quelconque des revendications 1 à 5. 10. Dispositif pour déterminer des informations cinématiques d'un aéronef porteur d'un radar à synthèse d'ouverture, ledit dispositif (D2) comportant: - des premiers moyens de calcul (1) pour déterminer le gradient Doppler du déphasage perturbateur des signaux radar détectés par ledit radar à synthèse d'ouverture; et - des seconds moyens de calcul (34) pour déduire dudit gradient Doppler lesdites informations cinématiques, caractérisé en ce que lesdits premiers moyens de calcul (1) mettent en oeuvre le procédé spécifié sous l'une quelconque des revendications 1 à 5.