Brevet FR2932570A1 — publié le 20091218

Procédé et dispositif d'estimation d'au moins une composante de vitesse d'une cible mobile. La présente invention concerne un procédé et un dispositif d'estimation d'au moins une des deux composantes de vitesse d'une cible mobile à partir d'une image radar synthétisée en distance et en azimut. Les techniques de détection de cibles mobiles usuellement appelées techniques MTI de l'anglais � Moving Target Indicator � visent à détecter une cible mobile et à estimer à partir d'un ou plusieurs signaux radar la position de la cible et les composantes radiale et/ou azimutale de sa vitesse. Les techniques MTI peuvent être classées en deux types, à savoir les techniques multi-images et les techniques mono-image. Les techniques multi-images consistent à utiliser plusieurs signaux radar de la même scène enregistrés par plusieurs plate-formes ou par plusieurs antennes alignées en azimut sur une plate-forme. Ces techniques permettent une bonne estimation de la position et de la vitesse de la cible, notamment en ce qui concerne sa composante de vitesse radiale. En contrepartie, elles nécessitent l'utilisation d'une plate-forme ayant plusieurs antennes ou de plusieurs plate-formes disposées selon des contraintes géométriques. Elles sont donc liées à des contraintes systèmes fortes et coûteuses. La plupart des techniques mono-image extrait les positions et/ou les composantes de vitesse à partir des données compressées en distance mais pas en azimut. L'invention a pour but un procédé et un dispositif alternatifs d'estimation d'au moins une des deux composantes de vitesse d'une cible mobile. A cet effet, l'invention a pour objet un procédé d'estimation d'au moins une des deux composantes de vitesse d'une cible mobile à partir d'une unique image radar synthétisée en distance et en azimut, lesdites composantes de la vitesse de la cible étant une composante azimutale et une composante perpendiculaire à la composante azimutale dite composante radiale, caractérisé en ce qu'il comporte les étapes suivantes : - recherche d'un motif de corrélation spatiale de la phase généré par le mouvement de la cible, à partir de l'image radar ; - recherche d'une parabole le long du motif de corrélation spatiale de la phase et modélisation de la parabole par un polynôme du second degré ; - détermination d'un paramètre quadratique dudit polynôme du second degré ; et estimation de la composante azimutale de la vitesse de la cible à partir dudit paramètre quadratique. Suivant des modes particuliers de mise en oeuvre, le procédé d'estimation comporte l'une ou plusieurs des caractéristiques suivantes, prises isolément ou en combinaison : - l'étape de recherche d'un motif de corrélation spatiale de la phase comporte les étapes suivantes : a) extraction d'une image de phase de l'image radar synthétisée ; b) recherche d'un motif de corrélation spatiale de la phase sur l'image de phase ; - l'étape de recherche d'une parabole comporte une étape de déroulement de la phase le long du motif de corrélation spatiale de la phase ; - l'étape de modélisation comporte une étape de détermination de l'extremum de la parabole, et d'un intervalle de traitement de part et d'autre de l'extremum ; - le procédé comporte en outre les étapes suivantes : - comparaison de la composante azimutale de la vitesse de la cible avec un nombre prédéfini, et lorsque la composante azimutale de la vitesse de la cible est inférieure au nombre prédéfini ; - détermination d'un paramètre linéaire du polynôme de second degré ; - détermination de la composante radiale de la vitesse de la cible à partir dudit paramètre linéaire ; - le procédé comporte également les étapes suivantes : g) extraction d'une image d'amplitude de l'image radar synthétisée ; h) recherche d'un motif d'amplitude sur l'image d'amplitude ; - l'étape d) comporte en outre une étape de binarisation de l'image d'amplitude à partir d'un seuil prédéfini ; - le seuil est égal à la valeur d'amplitude maximale de l'image d'amplitude divisée par un nombre supérieur ou égal à 3 ; - l'étape d) comporte en outre une étape de filtrage morphologique ; - le procédé comporte en outre les étapes suivantes : - détermination de la position en azimut du centre du motif d'amplitude ; - détermination du décalage spectral généré par le déplacement radial de la cible, à partir de ladite position du centre du motif d'amplitude et y ùyo . de la relation suivante : f o = a1 -2 a2 ( ) . a dans laquelle a1 est le paramètre linéaire du polynôme de second degré, a2 est le paramètre quadratique du polynôme de second degré, fo est le décalage spectral, (Va) est la vitesse d'un capteur imageur ayant capturé l'image radar, (y'ûyo) est la coordonnée en azimut du centre du motif d'amplitude ; - détermination de la composante radiale de la vitesse de la cible à partir du décalage spectral ; - le procédé comporte en outre les étapes suivantes : - mesure de la distance entre la position en azimut du centre du motif d'amplitude et la position en azimut de l'extremum du motif de phase ; - détermination du décalage spectral généré par le déplacement radial de la cible, à partir de ladite distance mesurée et de la relation : fo = 2a2 D dans laquelle D est la distance mesurée, a2 est le paramètre quadratique du polynôme de second degré, et fo est le décalage spectral ; et - estimation de la composante radiale de la vitesse de la cible à partir du décalage spectral ; - la composante radiale de la vitesse de la cible est obtenue à partir du décalage spectral et de la relation suivante : ilfo dans laquelle : 2 sin (9 ) - VYad est la composante radiale de la vitesse de la cible mobile, - fo est le décalage spectral issu du déplacement radial de la cible mobile, - 9 est l'angle d'incidence de prise de vue du radar, et - Â est la longueur d'onde du signal émis par le radar ; et - l'image radar est une image à très haute résolution. L'invention concerne également un dispositif d'estimation d'au moins une composante de vitesse d'une cible mobile à partir d'une unique image radar synthétisée en distance et en azimut, caractérisé en ce qu'il comporte : - une unité de recherche d'un motif de corrélation spatiale de la phase généré par le mouvement de la cible, à partir de l'image radar ; - une unité de recherche d'une parabole le long du motif de corrélation spatiale de la phase ; et - une unité de modélisation de la parabole par un polynôme du second degré ; ladite unité de modélisation étant propre à déterminer un paramètre quadratique dudit polynôme du second degré ; et à estimer la composante azimutale de la vitesse de la cible à partir dudit paramètre quadratique. L'invention sera mieux comprise à la lecture de la description qui va suivre, donnée uniquement à titre d'exemple et faite en se référant aux dessins annexés, sur lesquels : - la figure 1 est une représentation en trois dimensions de l'observation selon une direction d'incidence d'une scène par un radar à synthèse d'ouverture ; - la figure 2 est un diagramme des étapes du procédé d'estimation selon l'invention ; - la figure 3 est une image représentant l'amplitude d'une image radar comportant une cible mobile ; - la figure 4 est une image représentant la phase d'une image radar comportant une cible mobile ; - la figure 5 est un graphe représentant la phase en fonction de la direction azimutale de l'image radar, le long d'un motif de phase d'une cible mobile ; - la figure 6 est l'image illustrée sur la figure 3, binarisée et filtrée pour faire apparaître la trace d'une cible mobile ; - la figure 7 est un diagramme d'une variante d'une partie du procédé d'estimation selon l'invention ; - la figure 8 est une image de phase et une image d'amplitude représentant la même cible mobile ; - la figure 9 est un graphe représentant l'amplitude en fonction de la direction azimutale de l'image le long du motif de phase et un graphe représentant la phase en fonction de la direction azimutale de l'image le long du motif d'amplitude ; et - la figure 10 est une vue schématique d'un dispositif d'estimation selon l'invention. En référence à la figure 1, le procédé d'estimation selon l'invention est réalisé à partir d'une image radar 2 obtenue par un capteur imageur 4 monté sur une plate-forme aéroportée ou une plate-forme spatiale se déplaçant à une vitesse Va par rapport à la terre et se trouvant à une altitude H. Comme visible sur la figure 1, le capteur imageur 4 est un radar à visée latérale dont l'angle d'incidence de prise de vue est référencé O. Sur la figure 1, l'axe d'azimut est représenté par l'axe y, et l'axe de distance ou axe radial est représenté par l'axe R. L'axe R est parallèle à une droite passant par la cible mobile 6 et le capteur 4. L'axe x est perpendiculaire à l'axe y. Les axes x et y forment un plan dans lequel la cible mobile 6 se déplace. L'axe x et l'axe R définissent la géométrie de l'image radar. Comme les spécialistes radar désignent la composante selon l'axe x par le terme � composante radiale � ou � composante en distance �, nous avons également désigné la composante de la vitesse selon l'axe x par le terme � composante radiale � dans la présente demande bien qu'en réalité cette dénomination soit un abus de langage. La composante de la vitesse selon l'axe x est obtenue en divisant la composante de la vitesse selon l'axe R par sin O. Cette division correspond à une projection de l'axe R sur l'axe x. Le mouvement en distance de la cible est aussi appelé mouvement radial de la cible. C'est un mouvement selon l'axe x. Le mouvement en azimut de la cible est un mouvement selon l'axe y. Le capteur imageur 4 est un radar à synthèse d'ouverture (RSO). Il comprend une antenne qui émet des impulsions hyperfréquences vers la surface de la terre et réceptionne les signaux rétrodiffusés par la surface imagée. Il transmet les signaux rétrodiffusés vers un dispositif 5 d'enregistrement et de traitement au sol. Le dispositif d'enregistrement et de traitement 5 comprend deux unités disposées l'une à bord et l'autre au sol. Il réalise un traitement appelé � synthèse d'ouverture �, puis transmet l'image radar 2 dite image radar synthétisée complexe mono-vue (ou SLC, de l'anglais � Single Look Complex �) à un dispositif 8 d'estimation des composantes azimutale et radiale de la vitesse de la cible. Le principe de la synthèse d'ouverture consiste à exploiter la modulation de fréquence produite par les variations de distance capteur 4 - cible 6 (en azimut) et la modulation de l'onde émise (en distance) pour améliorer la résolution de l'image. D'un point de vue pratique, le signal radar rétrodiffusé par la cible 6 est convolué par un signal de référence bidimensionnel linéairement modulé en fréquence avec le même taux de modulation que le signal émis. La réponse impulsionnelle d'une cible 6 ponctuelle immobile dans une image radar est un sinus cardinal en deux dimensions. Les pixels correspondant à cette réponse contiennent une information d'amplitude et une information de phase, phase formée par la somme de la phase propre de la cible et de la phase liée au trajet de l'onde émise. La synthèse d'ouverture de cibles immobiles génère une image ayant une phase aléatoire identiquement distribuée entre 0 et 2n. La bande Doppler est la bande passante du signal en azimut. La synthèse d'ouverture d'un signal rétro diffusé par une cible 6 mobile uniquement selon l'axe azimutal génère un élargissement ou une diminution de la bande Doppler de ce signal qui fait apparaître sur l'image radar 2 un motif de corrélation spatiale de la phase qui prend la forme d'une parabole de phase résiduelle. La synthèse d'ouverture d'un signal rétro diffusé par une cible 6 mobile uniquement selon l'axe radial R génère un décalage du spectre Doppler de ce signal qui fait apparaître sur l'image radar 2 un motif de corrélation spatiale de la phase qui prend la forme d'une rampe de phase résiduelle. La synthèse d'ouverture d'une telle cible génère également une délocalisation selon l'azimut de la cible 6 sur l'image. Une étude théorique permet de mettre en évidence la relation entre d'une part les composantes radiale et azimutale de la vitesse de la cible 6 et d'autre part la parabole et la rampe de phase visible sur l'image radar synthétisée. En négligeant la compression en distance, le signal radar enregistré à l'instant t par le capteur 4 et rétrodiffusé par une cible fixe située en (xo, yo) vaut : (t, -02 cl)o + 27LKa to ùt) + CQ 2 sl (t) = A exp J avec t o Yo o V a o Dans laquelle Â est la longueur d'onde du signal émis, (xo, yo) est la position de la cible fixe, to est le début d'acquisition du signal rétrodiffusé, Ro est la distance minimale entre la cible 6 et le capteur 4, Ro =. IH2 +xô , Va est la vitesse du capteur 4 , Ka est le taux de modulation en azimut et A est une constante. Ce signal sl(t) est un signal modulé linéairement en fréquence de taux de modulation Ka et de fréquence porteuse nulle. La synthèse d'ouverture de la cible 6 lorsque celle-ci est fixe, consiste à effectuer la convolution du signal enregistré sl(t) modulé linéairement en fréquence avec le taux de modulation Ka avec un signal de référence modulé linéairement en fréquence avec le même taux de modulation Ka . Le signal obtenu en sortie du filtre de synthèse d'ouverture pour la cible 6 lorsque celle-ci est fixe, vaut : s(t) = A - .exp(j(4o + rp))sin, (7-Bd (t ù to )) a dans laquelle Bd est la bande passante du signal rétrodiffusé s(t) en azimut (bande Doppler). Le signal radar s2(t) enregistré à l'instant t par le capteur 4 et rétrodiffusé par la cible 6 lorsque celle-ci se déplace uniquement selon l'azimut avec une vitesse Vaz , vaut : (t ùt)z s2 (t) = A exp j cl)o + 27rKa v 2° + rp avec ù to 0 et donc Ka v Le signal h(t) obtenu en sortie du filtre de synthèse d'ouverture pour la cible 6 lorsque celle-ci se déplace en distance, est : h(t) = a B' 1 exp(j(~0 + Co)). sin c 9iB' t ù to + K Îo~~ exp(J ot) K, K, avec B'= Bd û fo , Bd étant la bande passante du signal rétrodiffusé s3(t). s3(t)=Aexp o +~)+2j1c fo(tùto)+Ka 2 Le signal résultant h(t) est délocalisé en azimut par rapport à la position yo de la cible en début d'acquisition de la grandeur suivante : Vrad Ro sin(e) ù fo Dy =(y û yo)= =Va dans laquelle y est la position azimutale Va Ka réelle de la cible mobile. La composante radiale Vrad de la vitesse de la cible mobile 6, induit un décalage spectral fo de la réponse de la cible qui se traduit par une rampe de phase résiduelle. Ce décalage spectral vaut : 2VYad .sin(0) 2 dans laquelle : - 9 désigne l'angle d'incidence de la prise de vue du capteur 4 ; est la composante radiale de la vitesse de la cible 6, - fo est le décalage spectral induit par le déplacement en distance de la cible, et - Â est la longueur d'onde du signal émis. Ainsi à partir de la connaissance du décalage spectral fo , il est possible de trouver la composante radiale VYad de la vitesse de la cible 6 par : À.f o 'ad 2 sin (9 Lorsque la cible 6 est animée d'un mouvement en azimut et en distance, le signal i(t) résultant de la synthèse d'ouverture vaut : ~ = Aexp2 jTC f2 (t û t0 + t')+ K' (t t + t ,z 2 pour t ù t0 La composante de la vitesse radiale induit un décalage spectral f0 égal à 2V,,ad sin 9 2 Ce résultat rejoint le résultat obtenu lorsque la cible se déplace uniquement selon l'azimut et uniquement en distance. La composante de la vitesse azimutale induit une délocalisation en azimut se traduisant par un chip résiduel de modulation K' et une délocalisation en azimut provenant du mouvement en azimut et en distance de Ay = (Va ùV" ) 2 K K a,v a La phase du signal i(t) issu de la synthèse d'ouverture est ~(t)='F0+ f0(tùt0+t')+K'(tùt0+t')=a0+ait+azt2. Elle se factorise en :W= +K' tùt +t'+ f0 ù~ f o ° 2K' ) 2K' La phase ço(t) passe par une valeur extrême pour une position temporelle t0 ùt'ù2K . En conséquence, l'amplitude du signal rétrodiffusé par la cible mobile est délocalisée par rapport à la phase du signal rétrodiffusé par cette même cible. En d'autres termes, le centre du motif d'amplitude de la cible mobile est délocalisé en azimut par rapport à la valeur minimale de la phase de la cible de 2K . En d'autres termes, la composante radiale de la vitesse de la cible a pour effet de décaler de ZK la valeur extrême de la phase par rapport à la valeur extrême de l'amplitude dans l'image radar 2 dans le cas d'un mouvement à la fois radial et azimutal. Une mesure du décalage entre le centre du motif d'amplitude et le centre du motif de phase permet de déterminer le décalage spectral f0 si le taux de modulation K' est connu. La phase ço(t) associée à une cible mobile est donc modélisée par : - le terme çp1(t)= fo(t+t'ûto) qui est centré sur la coordonnée temporelle et qui évolue linéairement en fonction du temps. Ce terme n'est fonction que de la composante radiale de la vitesse, - le terme q2(t)=K'(tûto+t')2centré sur la coordonnée temporelle (t'ûto) et qui évolue quadratiquement en fonction du temps. Ce terme est fonction de la composante azimutale de la vitesse. - le troisième terme Yo est constant au cours du temps. Il est non impacté par la vitesse de la cible. Les deux composantes de vitesse d'une cible mobile peuvent ainsi être déduites des paramètres des composantes quadratique (composante azimutale) et linéaire (composante radiale) de la phase de la cible. En référence à la figure 2, le procédé d'estimation selon l'invention débute par une étape 10 d'extraction d'une image d'amplitude 12 et d'une image de phase 14 de l'image radar 2. De façon connue en soi, le signal radar réceptionné par l'antenne subit une transformation analytique pour obtenir un signal analytique ayant une composante réelle et une composante imaginaire. Ce signal analytique est traité pour produire une image d'amplitude 12 et une image de phase 14. Selon l'invention, l'image radar 2 est soit une image radar à haute résolution (résolution inférieure à 2 mètres) ou une image radar à très haute résolution (résolution inférieure à 1 mètre). Au cours d'une étape 24, un motif de corrélation spatiale de la phase, appelé motif de phase 26 est recherché sur l'image de phase 14. Cette recherche est réalisée visuellement ou par ordinateur à l'aide d'un logiciel de recherche des sinusoïdes (franges). Ce motif de phase 26 est illustré sur la figure 4. Le motif de phase 26 est généré par une cible mobile. Au cours d'une étape 28, la phase çp(t) est déroulée le long du motif de phase 26 (c'est-à-dire selon l'azimut), comme visible sur la figure 5. Au cours d'une étape 30, une parabole 32 est recherchée sur le profil de la phase déroulée çp(t). Cette parabole 32 est centrée. L'intervalle de traitement de la parabole est réduit autour de l'extrémité E pour faciliter le traitement de celle-ci, au cours d'une étape 34. Au cours d'une étape 36, la parabole 32 est modélisée par un polynôme de 5 second degré, par exemple par une méthode des moindres carrés. Cette modélisation de la phase ~p(t) à un instant t est notée çp(t)=ao+alt+a2t2. Le paramètre quadratique a2 du polynôme de second degré est déterminé à partir de la parabole 32. La composante azimutale VaZ de la vitesse de la cible est directement 10 déduite au cours d'une étape 38 du paramètre quadratique a2 du polynôme de second degré. Comme démontré précédemment, le paramètre quadratique a2 du polynôme de second degré est égal àK'. La vitesse en azimut VaZ de la cible mobile 6 est alors déduite de K' par la relation : ~îR0 2 1 LVa2 + ' dans laquelle Va est la vitesse du capteur 4. Au cours d'une étape 40, la composante azimutale Vaz de la vitesse de la cible est comparée avec un nombre prédéfini N proche de 0. Lorsque celle-ci est inférieure au nombre prédéfini N, le paramètre linéaire 20 a, du polynôme de second degré est considéré comme étant égal au décalage spectral fo , et au cours d'une étape 42, ce paramètre linéaire a, est déterminé. Au cours d'une étape 44, la composante radiale Vrad de la vitesse de la cible est obtenue à partir de la relation y ~fo °d ù 2sin(9) dans laquelle : 25 - Vrad est la composante radiale de la vitesse de la cible mobile, - fo est le décalage spectral, - 9 est l'angle d'incidence de prise de vue du capteur 4, et - Â est la longueur d'onde du signal émis. 15 Vaz = Va ù 2Ro Lorsque la composante azimutale Vaz de la vitesse de la cible n'est pas inférieure au nombre prédéfini N, il est possible d'obtenir le décalage spectral fo en recherchant le centre d'un motif d'amplitude 22 d'une image d'amplitude. A cet effet, au cours d'une étape 45, un motif d'amplitude 22, généralement appelé trace, est recherché sur l'image d'amplitude 12. Cette recherche est effectuée visuellement à partir de la position du motif de phase 26. En effet, le décalage entre le motif d'amplitude 22 et le motif de phase 26 est faible. Ainsi, les deux motifs sont très proches l'un de l'autre. En variante, cette recherche est effectuée par traitement informatique. A cet effet, l'image d'amplitude 12 est binarisée à partir d'un seuil prédéfini V. Ce seuil V est égal par exemple à trois fois la variance 6 de l'image d'amplitude additionnée à la moyenne radiométrique ,u de cette image. En variante, le seuil V est égal à la valeur d'amplitude maximale de l'image d'amplitude 12 divisée par un nombre supérieur ou égal à 3 et de préférence égal à 4. L'image d'amplitude 12 est ensuite filtrée par un filtrage morphologique d'ouverture. Ce type de filtrage est décrit dans les livres intitulés � Image analysis and mathematical morphology �, volumes 1 et 2 Academic press, Londres. Comme visible sur la figure 3, ce motif d'amplitude 22 présente la forme d'un segment de droite. Il s'étend selon la direction d'azimut. Il présente une amplitude supérieure à celle du reste de l'image d'amplitude 12, lorsque le rapport signal sur bruit est suffisant. Au cours d'une étape 46, le centre Ca du motif d'amplitude 22 est obtenu par détermination de la position y'-yo en azimut. A cet effet, les coordonnées de colonne des extrémités du motif d'amplitude 22 sont déterminées, comme visible sur la figure 6, à partir de l'image d'amplitude 12 binarisée puis filtrée par le filtre morphologique. La variable temporelle (t'-to) résulte de la relation : (t' - to) = (y'-yo) -Q dans laquelle t 1 1 1 +- fo VQ 2VQ KQ vrad Ka Au cours d'une étape 48, le décalage spectral fo est calculé à partir de la variable temporelle (t' - to), du paramètre quadratique a2, et de la relation suivante : fo = a1 -2 a2(t' - to). Enfin, au cours d'une étape 50, la composante radiale Vrad de la vitesse de la cible est obtenue à partir de la relation : y = ~fo °d ù 2 sin (9 ) Selon une variante de réalisation du procédé selon l'invention illustrée sur les figures 7 à 9, la composante radiale Vraa de la vitesse de la cible 6 comprend les étapes suivantes en lieu et place des étapes 46 à 50. Au cours d'une étape 52, la position Ca du centre du motif d'amplitude 22 est recherchée ainsi que la position de l'extremum E de la parabole 32 dans la phase déroulée. Puis, au cours d'une étape 54, le décalage ou la distance D entre d'une part la position y'-yo du motif d'amplitude 22 et la position de l'extremum E de la parabole 32 est mesurée, comme visible sur les figures 8 et 9. La position y'-yo représente la position qu'aurait la cible 6 dans l'image d'amplitude 12 si la synthèse d'ouverture était effectuée de façon à focaliser la cible mobile. Comme cette position n'est pas délocalisée dans l'image d'amplitude et compte tenu des études théoriques décrites précédemment, la distance D vaut t- . Au cours d'une étape 56, le décalage spectral fo est calculé à partir de ladite distance D mesurée, du paramètre quadratique a2 obtenu par modélisation de la parabole 32 et de la relation : fo = 2a2 D Puis, la composante radiale Vradt de la vitesse de la cible est déterminée à partir du décalage spectral, comme décrit précédemment. Avantageusement, ce procédé est mis en oeuvre à l'aide d'une unique image radar 2. Avantageusement, ce procédé est mis en oeuvre à partir d'une image radar synthétisée (SLC, de l'anglais � Single Look Complex �) et non à partir de données brutes ou de données compressées en distance mais pas en azimut. Or, il est plus facile de se procurer sur le marché actuel des images radar synthétisées plutôt que des images non compressées ou partiellement compressées. De plus, le traitement d'images non compressées nécessite une expertise particulière. Cette méthode permet d'obtenir une valeur très précise en ce qui concerne la composante azimutale de la vitesse du mobile. Elle est peu coûteuse en nombre d'images radar nécessaires et en terme de calcul. Elle peut être associée à d'autres méthodes de détermination de vitesse. Par exemple, il est possible de déterminer en premier lieu la position et/ou les deux composantes de vitesse avec cette méthode puis d'appliquer une méthode plus précise, par exemple une méthode d'autofocus, en vue d'affiner les résultats. Dans ce cas, l'utilisation de la méthode selon l'invention permet de restreindre la durée et la complexité en termes de calcul de la méthode d'autofocus. Elle peut également être associée à des techniques multi images qui généralement obtiennent des résultats de faible qualité en ce qui concerne la composante azimutale de la vitesse du mobile. Le dispositif d'estimation 8, illustré sur la figure 10 comprend une unité 62 d'extraction d'une image d'amplitude 12 et d'une image de phase 14, une unité de binarisation 64 reliée à l'unité d'extraction 62 et un filtre 66 morphologique d'ouverture relié à l'unité de binarisation 64. Le dispositif d'estimation 8 comporte en outre une unité 70 de recherche du motif d'amplitude 22 sur l'image d'amplitude et du motif de phase 24 sur l'image de phase 14, une unité 72 de recherche d'une parabole le long du motif de phase rp(t) connectée à l'unité 70, et enfin une unité de calcul 76 reliée à l'unité 70. L'unité de calcul 76 est apte à modéliser la parabole par un polynôme de second degré comme explicité ci-dessus. REVENDICATIONS1.- Procédé d'estimation d'au moins une des deux composantes (Vaz, Vrad) de vitesse d'une cible mobile (6) à partir d'une unique image radar (2) synthétisée en distance et en azimut, lesdites composantes de la vitesse de la cible (6) étant une composante azimutale (Vaz) et une composante perpendiculaire à la composante azimutale (Vaz) dite composante radiale (Vrad), caractérisé en ce qu'il comporte les étapes suivantes : - recherche (10, 24) d'un motif (26) de corrélation spatiale de la phase (çp(t)) généré par le mouvement de la cible (6), à partir de l'image radar (2) ; - recherche (28, 30) d'une parabole (32) le long du motif (26) de corrélation spatiale de la phase (çp(t)) et modélisation (34, 36) de la parabole (32) par un polynôme du second degré ; - détermination (32, 36) d'un paramètre quadratique (a2) dudit polynôme du second degré ; et estimation (38) de la composante azimutale (Vaz) de la vitesse de la cible (6) à partir dudit paramètre quadratique (a2). 2.- Procédé d'estimation selon la revendication 1, dans lequel l'étape de recherche d'un motif de corrélation spatiale (26) de la phase (ço(t)) comporte les étapes suivantes : a) extraction (10) d'une image de phase (14) de l'image radar synthétisée (2) ; b) recherche (24) d'un motif (26) de corrélation spatiale de la phase sur l'image de phase (14). 3.- Procédé d'estimation selon l'une quelconque des revendications précédentes, dans lequel l'étape de recherche (28, 30) d'une parabole (32) comporte une étape de déroulement (28) de la phase (ço(t)) le long du motif (26) de corrélation spatiale de la phase. 4.- Procédé d'estimation selon l'une quelconque des revendications précédentes, dans lequel l'étape de modélisation (34, 36) comporte une étape (34) de détermination de l'extremum (E) de la parabole (32), et d'un intervalle de traitement de part et d'autre de l'extremum (E). 5.- Procédé d'estimation selon l'une quelconque des revendications précédentes, qui comporte en outre les étapes suivantes :- comparaison (40) de la composante azimutale (Vaz) de la vitesse de la cible (6) avec un nombre prédéfini (N), et lorsque la composante azimutale de la vitesse de la cible est inférieure au nombre prédéfini (N); - détermination (42) d'un paramètre linéaire (a1) du polynôme de second degré ; - détermination (44) de la composante radiale (Vrad) de la vitesse de la cible (6) à partir dudit paramètre linéaire (a1). 6.- Procédé d'estimation selon l'une quelconque des revendications précédentes, qui comporte en outre les étapes suivantes : c) extraction (10) d'une image d'amplitude (12) de l'image radar synthétisée (2) ; d) recherche (45) d'un motif d'amplitude (22) sur l'image d'amplitude (22). 7.- Procédé d'estimation selon la revendication 6, dans lequel l'étape d) comporte en outre une étape de binarisation de l'image d'amplitude (12) à partir d'un seuil prédéfini (V). 8.- Procédé d'estimation selon la revendication 7, dans lequel le seuil (V) est égal à la valeur d'amplitude maximale de l'image d'amplitude (12) divisée par un nombre supérieur ou égal à 3. 9.- Procédé d'estimation selon l'une quelconque des revendications 6 à 8, dans lequel l'étape d) comporte en outre une étape de filtrage morphologique. 10.- Procédé d'estimation selon l'une quelconque des revendications 6 à 9, qui comporte en outre les étapes suivantes : - détermination (46) de la position (y'-yo) en azimut du centre (Ca) du motif d'amplitude (22); - détermination (48) du décalage spectral (fo) généré par le déplacement radial de la cible (6), à partir de ladite position (y'-yo) du centre (Ca) du motif Y ùYo . d'amplitude et de la relation suivante : f o = a1 -2 a2 ( ) . a dans laquelle a1 est le paramètre linéaire du polynôme de second degré, a2 est le paramètre quadratique du polynôme de second degré, fo est le décalage spectral, (Va) est la vitesse d'un capteur imageur (4) ayant capturé l'image radar(2), (y'ûyo) est la coordonnée en azimut du centre (Ca) du motif d'amplitude (22) ; - détermination (50) de la composante radiale (Vrad) de la vitesse de la cible (6) à partir du décalage spectral (fo ). 11.- Procédé d'estimation selon les revendications 4 et 6 prises en combinaison, qui comporte en outre les étapes suivantes : - mesure (52) de la distance (D) entre la position (y'-yo) en azimut du centre (Ca) du motif d'amplitude (22) et la position en azimut de l'extremum (E) du motif de phase (26) ; - détermination (54) du décalage spectral (fo) généré par le déplacement radial de la cible (6), à partir de ladite distance mesurée (D) et de la relation : fo = 2a2 D dans laquelle D est la distance mesurée, a2 est le paramètre quadratique du polynôme de second degré, et fo est le décalage spectral ; et - estimation (56) de la composante radiale (Vrad) de la vitesse de la cible (6) à partir du décalage spectral (fo ). 12.- Procédé d'estimation selon l'une quelconque des revendications 10 et 11, dans laquelle la composante radiale (Vrad) de la vitesse de la cible (6) est obtenue à partir du décalage spectral (fo) et de la relation suivante : ilf° dans laquelle : 2 sin (9 ) - VYQd est la composante radiale de la vitesse de la cible mobile, - fo est le décalage spectral issu du déplacement radial de la cible mobile, - 9 est l'angle d'incidence de prise de vue du radar, et - Â est la longueur d'onde du signal émis par le radar. 13.- Procédé d'estimation selon l'une quelconque des revendications précédentes, dans lequel l'image radar (2) est une image à très haute résolution. 14.- Dispositif (8) d'estimation d'au moins une composante (Vaz, Vrad) de vitesse d'une cible mobile (6) à partir d'une unique image radar (2) synthétisée en distance et en azimut, caractérisé en ce qu'il comporte :- une unité de recherche (62, 64, 66, 70) d'un motif (22, 26) de corrélation spatiale de la phase (çp(t)) généré par le mouvement de la cible (6), à partir de l'image radar (2) ; - une unité de recherche (72) d'une parabole (32) le long du motif (22, 26) de corrélation spatiale de la phase (qp(t)) ; et - une unité de modélisation (76) de la parabole (32) par un polynôme du second degré ; ladite unité de modélisation (76) étant propre à déterminer un paramètre quadratique (a2) dudit polynôme du second degré ; et à estimer la composante azimutale (Vaz) de la vitesse de la cible (6) à partir dudit paramètre quadratique (a2).