Brevet FR3005168A1 — publié le 20141031

PROCEDE D'ELIMINATION D'ECHOS PARASITES EN IMAGERIE SAR La présente invention concerne le domaine général des radars aéroportés ou spatiaux à synthèse d'ouverture. La présente invention concerne plus particulièrement un procédé d'élimination d'échos parasites sur les images SAR. Les radars à synthèse d'ouverture (ou SAR pour Synthetic Aperture Radar selon la terminologie anglo-saxonne) sont largement utilisés dans les activités d'analyse à distance des caractéristiques physiques des zones terrestres et maritimes à partir de leur réflectivité radar. Ces activités sont appelées " Radar remote sensing " en anglais. Une des activités dans ce domaine est la détection de nappes d'hydrocarbures sur la mer par différence de niveau de réflectivité radar entre les zones polluées et les zones non polluées. Des études ont été réalisées pour évaluer la perturbation apportée de la présence de pluie pour cette analyse. C'est le cas par exemple de l'étude " Rain-induced modification of SAR performance» de P. Ferrazzoli et G. Schiavon - Adv Space Res vol 7 N°11 1987. Dans cette étude, les défauts liés à la présence de pluie sont présentés mais l'analyse effectuée ne traite pas le fait que les échos de pluie peuvent avoir un spectre Doppler décalé par rapport au spectre des échos de sol de la zone observée. De plus, cette étude ne propose pas de solution pour réduire la perturbation. D'autres études traitent de l'impact de la présence d'échos mobiles dans la zone observée par un radar SAR (voir notamment " Principles and applications of imaging radars'> - Manual of Remote Sensing - Third Edition - Volume 2) et proposent des solutions pour pouvoir détecter des échos mobiles à partir de traitements SAR. Cependant, ces techniques, qui s'apparentent aux techniques de détection de cibles mouvantes (ou MTI pour Moving Target Indicator selon la terminologie anglo saxonne), ne sont effectives que pour détecter des cibles ayant une vitesse radiale suffisante pour que leur fréquence Doppler sortent du spectre Doppler des échos fixes observés par la technique SAR. Ces techniques ne sont pas applicables dans le cas où le domaine de variation des vitesses radiales relatives des éléments perturbateurs comme les gouttes de pluie ne se dégage pas du domaine des vitesses radiales relatives des éléments de la scène observée qui sont à la même position géographique. Des techniques d'amélioration de la texture des images de scènes SAR maritimes, rurales ou forestière ont également été décrites notamment dans l'ouvrage " Principles and applications of imaging radars» précédemment cité. Ces méthodes sont appliquées sur des images SAR déjà construites. Elles s'appuient sur la construction du spectre 2D de l'image SAR et sur l'élimination d'une partie du spectre qui correspond au speckle, pour ne conserver que les éléments ayant une texture particulière comme le motif de vagues. Ces méthodes ne sont applicables que quand l'image attendue a une texture particulière. Ce n'est pas le cas par exemple, lorsque l'on observe d'une mer belle à peu agitée, à faible angle de rasance en polarisation horizontale. Enfin, dans les situations où le phénomène d'intérêt ne conduit pas à des images ayant une texture particulière, des méthodes de segmentation peuvent être utilisées. Ces méthodes sont basées sur l'identification de lois de densités de probabilités différentes sur le fond de l'image et sur la zone recherchée. Cette segmentation peut aussi être complémentée par des opérations 25 dites de «morphologie mathématique » permettant de lisser l'image résultant de la segmentation en partant de l'hypothèse que les zones recherchées ont une étendue de plusieurs pixels de l'image. Un inconvénient de ces solutions est qu'elles sont appliquées une fois 30 que l'image SAR est créée et que le phénomène parasite a contribué à l'augmentation de tous les niveaux des pixels correspondant à des zones de faible réflectivité que l'on cherche à identifier sur la carte et donc à réduire le contraste entre les zones de faible réflectivité et les zones de plus faible réflectivité. 35 Un but de l'invention est notamment de corriger tout ou partie des inconvénients précités en proposant un procédé de traitement d'imagerie SAR permettant d'atténuer les effets d'échos parasites. A cet effet, l'invention a pour objet un procédé d'élimination d'échos parasites en imagerie SAR comprenant : - une étape de définition des paramètres d'imagerie SAR, - une étape de calcul du spectre du signal reçu, dans chaque case distance, d'une zone d'intérêt, - une étape de filtrage des spectres dans chaque case distance réalisée par un filtre passe bande centré sur la fréquence centrale fix des échos fixes rétrodiffusés par la zone d'intérêt, ladite fréquence centrale fp° étant obtenu par corrélation du spectre mesuré avec la forme théorique du spectre obtenue à partir des caractéristiques du diagramme de l'antenne radar, - une étape de formation d'imagette SAR et - une étape de concaténation des imagettes SAR pour former l'image SAR finale. Selon un mode de mise en oeuvre, en cas de recouvrement de pixels lors de l'étape de concaténation, cette dernière comprend une étape de moyennage de la valeur des pixels de chaque point de l'image SAR finale recouvert par plusieurs imagettes. d'élimination d'échos parasites sur les images SAR tel que décrit par 25 les revendications. De façon avantageuse, la présente invention permet d'améliorer le contraste sur les images SAR en réduisant le niveau des parasites dus à des échos radar naturels, par exemple liés aux particules de pluie ou de 30 déferlement de vagues, qui ont un spectre de vitesse en général plus large que celui des phénomènes ou objets observés qui sont statiques ou faiblement mobiles. D'autres particularités et avantages de la présente invention 35 apparaîtront plus clairement à la lecture de la description ci-après, donnée à titre illustratif et non limitatif, et faite en référence aux dessins annexés, dans lesquels : La figure 1 représente un exemple de plateforme effectuant une imagerie de type SAR d'une cible ; La figure 2 est un exemple de représentation graphique de l'amplitude d'un spectre de signal reçu par un radar de type SAR; La figure 3 est un exemple de représentation graphique de l'amplitude d'un spectre Doppler dans le cas où la rétrodiffusion n'est pas uniforme ; Les figures 4a et 4b sont des exemples de représentation graphique de l'amplitude de spectres Doppler parasités ; La figure 5 illustre le principe de l'invention ; La figure 6 représente un organigramme de principe du procédé selon l'invention ; La figure 7 illustre le représente ; Les figures 8a et 8b sont respectivement des exemples de représentation graphique de l'amplitude d'un spectre " normal » et de celle d'un spectre «pollué ". Le procédé selon l'invention peut être mise en oeuvre par exemple, par des porteurs mobiles de type aéronefs ou satellite embarquant un dispositif de type radar à synthèse d'ouverture (ou SAR) pour la surveillance d'une zone d'intérêt. La cartographie radar à haute résolution à partir de plateformes 10 aéroportées ou spatiales est réalisée en général par la technique d'imagerie SAR. En référence à la figure 1, le principe de la cartographie SAR est rappelé brièvement. Elle est adaptée à la cartographie de points fixes. L'imagerie radar SAR utilise la fréquence Doppler des cibles pour les séparer angulairement. Si on considère une cible mobile M illuminée par le lobe d'antenne d'un radar à synthèse d'ouverture, cette cible a une fréquence Doppler fp qui dépend de son angle 0 par rapport au vecteur vitesse Va de la plateforme 10 telle que : .ii(o) fD = (Equation 1) Pour une plateforme 10 ayant un vecteur vitesse Va horizontal et une 5 visée rasante, on peut simplifier l'équation précédente qui devient : fD= 2.va O. (Equation 1 bis ) OU f 2.va p = . sin a. (avec a=i-c/2 - 0) (Equation 1 ter) Où: 2 représente la longueur d'onde du radar 10 et ii le vecteur unitaire en direction de la cible M. Ainsi, la zone éclairée par le lobe d'antenne, Af D ,lobe a une étendue Doppler qui est directement reliée à l'ouverture à -3dB du lobe d'antenne AG3db : 2a .ù" (0) Afo,lobe = A 03db 15 ou plus simplement : . AfDJobe = 2v sm 03,03db (Equation 2) Le spectre utile à la formation de l'image SAR est ainsi compris dans la bande utile AfD,lobe - Sa fréquence centrale correspond à la fréquence Doppler au centre du lobe d'antenne : 2i; (0) 20 fpo = (Equation 3) Il est à noter qu'en pratique, la fréquence Doppler frn n'est pas connue avec une grande précision à cause des différentes erreurs sur la valeur et la direction du vecteur Va, et sur la direction de l'antenne, du fait notamment des erreurs d'harmonisation. 25 La résolution Doppler sfd peut être obtenue en dérivant l'équation 1. 8fd =11 Ti =2.va .sin 0.80 (Equation 4) Où Ti représente le temps d'intégration cohérente La résolution transverse 5rt peut être donnée par la formule : cYrt = DA' = (Equations 5) 2.va. sin O.Ti Où D représente la distance de la zone à imager Dans le cas où la cible M est mobile, comme par exemple dans le cas de gouttes de pluie se déplaçant sous l'effet du vent ou celui d'un mouton à la surface de la mer, l'équation (1 bis) devient : 2.va 2.Vr fp = . cos O. + (Equation 6) Â, Où Vr représente la vitesse radiale de l'objet. La figure 2 donne un exemple de représentation spectrale de l'allure de l'amplitude du signal reçu par un radar à synthèse d'ouverture, des différents points fixes localisés dans sa cellule de résolution. Cette représentation correspond à un cas simplifié où le coefficient de rétrodiffusion radar est constant. Nous avons vu précédemment que pour un axe cible-plateforme sensiblement perpendiculaire à la trajectoire de la plateforme 10, ce qui correspond à une valeur de l'angle a proche de zéro, il existait une relation entre la fréquence Doppler et l'angle a (cf. équation 1ter). Dans le cas d'un coefficient de rétrodiffusion uniforme, l'allure du centre du spectre du signal reçu de la cellule de résolution correspond à la forme du lobe principal de l'antenne au carré car les signaux reçus de la cellule de résolution radar sont plus ou moins forts selon leur dépointage dans le lobe. Afin de former l'image SAR, après le calcul du spectre du signal reçu dans chaque case distance un filtrage est effectué pour ne conserver que la bande utile de ce spectre. Classiquement cette bande utile correspond à la bande de fréquence pour laquelle l'amplitude du signal est supérieure à l'amplitude maximale moins trois décibels. Un problème est que dans cette bande utile il peut exister un ou plusieurs spectres parasites. A titre d'exemple, la figure 3 représente l'amplitude d'un spectre Doppler dans le cas où la rétrodiffusion dans la cellule de résolution (11) n'est pas uniforme mais présente une valeur plus faible pour un azimut donné. Cette atténuation (30) peut, par exemple, être due à la présence d'une nappe de pétrole, d'hydrométéores ou de speckle. Dans le cas représenté, la rétrodiffusion présente une atténuation (30) ponctuelle de 10 dB de la réflectivité de la mer pour un azimut de 910. Cette échelle d'abscisse comprise entre 87° et 93° correspond à une variation de fréquence Doppler de -350 Hz à +350 Hz. En présence d'hydrométéores et/ou de vagues déferlantes, lorsqu'une plateforme 10 aéroportée ou spatiale effectue une imagerie SAR, plusieurs objets de vitesses radiales Vr différentes peuvent être présents dans la cellule de résolution (11) du radar. Les valeurs de Vr sont en général dispersées autour d'une vitesse radiale Vro (Vitesse radiale moyenne) avec une certaine largeur AVr. De manière schématique, Le spectre Doppler dans la cellule (11) est constitué de la somme du spectre des échos supposés fixes, comme par exemple surface de la mer dans le cas d'une imagerie maritime ou le sol dans un cas terrestre, et des échos mobiles dus par exemple à la pluie et/ou au déferlement des vagues. A titre illustratif, la figure 4a représente, suivant une échelle linéaire, l'amplitude en fonction de l'azimut, d'un exemple de spectre d'une mer présentant une nappe de pétrole 41, d'un exemple de spectre d'échos de pluie 42 ainsi que le spectre résultant de la somme des deux spectres précédents. On peut noter que la puissance des échos de pluie correspond à une puissance de 3dB sous la puissance des échos fixes. La figure 4b représente l'amplitude du spectre résultant de la somme des échos de pluie et des échos de mer présentant une nappe de pétrole, en fonction de l'azimut, avec une échelle logarithmique. Par rapport au cas de la figure 3, on peut noter que pour la partie atténuée, le contraste entre la case spectrale atténuée et les cases spectrales voisines n'est plus que de 3dB au lieu de 10 dB précédemment.35 Afin de réduire la contribution du spectre parasite, l'invention consiste à réduire la taille de la bande utile pour ne conserver qu'une bande de fréquences minimale Af du spectre pour former l'image SAR. Cette réduction de la bande de fréquence peut être assimilée à un filtrage fréquentiel. Une conséquence de ce filtrage est une réduction de la zone image au sol sur l'axe transverse. De façon à combler cette perte d'information dans l'image SAR due au filtrage, le procédé selon l'invention va consister à fractionner le temps d'éclairement de la zone à imager et à former plusieurs sous-images SAR ou imagettes SAR 51 qui seront ensuite juxtaposées afin de former l'image SAR finale. La figure 5 illustre le principe de fractionnement du temps d'éclairement et du traitement en sous images. Cette figure présente une vue schématique d'un exemple de plateforme 10 effectuant une imagerie SAR d'une zone d'intérêt 50. On suppose que la plateforme 10 se déplace à vitesse et route constante pour passer successivement des positions Al à A7. L'antenne du radar à synthèse d'ouverture est pointée à environ 90° de la route empruntée par la plateforme 10. Des analyses Doppler sont réalisées à partir de la réception de Nrec impulsions radar quand la plateforme se déplace sur un segment de droite centré sur la position Al ou A2, ..., ou A7 pour former une imagette SAR 51. La figure 6 représente un organigramme d'étapes possible du procédé selon l'invention. Ce procédé comporte principalement une étape Etpl de définition des paramètres d'imagerie SAR, une étape Etp2 de calcul du spectre du signal reçu dans chaque case distance, une étape Etp3 de filtrage des spectres dans chaque case distance, une étape Etp4 de formation d'imagette SAR 51 et une étape Etp5 de concaténation des imagettes 51 pour former l'image SAR finale. On suppose que la forme d'onde émise par le radar à synthèse d'ouverture a été choisie de façon à obtenir la résolution souhaitée sur l'axe radial et à limiter les ambigüités en distance radiale et en distance transverse. A partir de paramètres d'entrée comme par exemple la vitesse Va du porteur, la distance D de la zone à imager par rapport à la plateforme 10, les paramètres de traitement SAR vont être définis de façon automatique. L'élément d'entrée principal pour définir le traitement SAR sur l'axe 5 transverse est la résolution Ôrt sur cet axe. Suivant l'application radar souhaitée, cette résolution pourra être choisie plus ou moins fine. Le choix de cette résolution conduit à la valeur du temps d'intégration cohérente Ti en utilisant par exemple l'équation 5 citées précédemment. 10 Pendant ce temps Th la plateforme 10 évoluant à une vitesse Va, se déplace d'une distance Lv égale à : Lv=Va.Ti (Equation 7) En référence à la figure 7, si on considère un point M situé à une distance D de la plateforme et à environ 90° de la route de la plateforme, 15 l'azimut de ce point M au milieu de l'intervalle de durée TI va varier d'une valeur égale à: A0Va.Ti = (Equation 8) D En utilisant les équations 5, avec une valeur de 0 d'environ 90° on obtient: 20 M = (Equation 9) 2* D*80 2.8rt Ainsi, une portion Af du spectre Doppler, centrée autour de la fréquence centrale fp°, est suffisante pour réaliser l'image du point M. Cette étendue Doppler Af peut être obtenue en dérivant l'équation (1 bis) : Af =2.Va .sin 0.A0 (Equation 10) 25 Ainsi, pour un angle 0 voisin de 90° on obtient: Af =2.Va .A0 (Equation 11) Ce qui, en utilisant les équations 5, avec une valeur de 0 voisine de 90°, et l'équation 9, conduit à: Va (Equation 12) Af = Va-- D* 80 grt Les valeurs de l'étendue Doppler Af et de la largeur azimutale à traiter dans le lobe antenne AO correspondent aux valeurs minimales permettant de générer l'image d'un point M unique sur l'axe transverse. Dans la technique générale d'imagerie SAR, il est préférable de réaliser, non pas l'image d'un seul point M, mais de l'ensemble des points étendus en transverse autour du point M. Cela permet de générer une image d'une zone d'étendue transverse égale au déplacement du porteur et ainsi 10 d'assurer la formation d'image au fur et à mesure de l'avancement dudit porteur. Par ailleurs, s'il est nécessaire de réaliser un traitement multi vues en azimut de manière à améliorer la résolution radiométrique, les largeurs azimutales AO et les étendues spectrales Af associées peuvent être 15 multipliées par un facteur Nv correspondant au nombre de vues. Il faut également noter que les paramètres des équations définissant les paramètres SAR dépendent de la distance D entre le radar et la cible et qu'il est souvent nécessaire de partager la zone à imager en plusieurs " fauchées » en distance. 20 Pendant la durée où l'avion se déplace sur une longueur de segment telle que définie précédemment, le traitement réalise une analyse spectrale dans chaque case distance correspondant à la fauchée. Cette analyse est réalisée au cours d'une étape Etp2, de façon classique par les méthodes 25 connues de l'homme du métier. De façon connue, une correction des migrations distance est effectuée de façon à tenir compte du fait que la distance parcourue par les échos radar évolue au cour du déplacement de la plateforme 10. La densité spectrale mesurée est ensuite estimée en effectuant une 30 moyenne des spectres sur un nombre prédéterminé N de cases distance. Il s'ensuit une opération de traitement du signal destinée à définir si la densité spectrale calculée est " normale » ou polluée par des phénomènes parasites par exemple d'origine météorologiques. Pour cela, cette opération prend en compte la densité spectrale théorique attendue définie en fonction des caractéristiques du diagramme antenne et la compare avec la densité spectrale mesurée. La définition précise du caractère " normale'> ou non de cette opération dépend de l'application et de l'ampleur des perturbations tolérées. 5 Elle consiste à définir si la densité spectrale diffère trop de l'attendue. Le traitement estime ensuite la fréquence centrale fDo du spectre des échos fixes. Pour cela, l'invention propose une méthode adaptative suivant si le spectre mesuré est considéré comme pollué ou non par des échos parasites. 10 Si la densité spectrale est considérée comme " normale ", la valeur retenue pour la fréquence centrale fDo est obtenue en décalant la densité spectrale attendue sur l'axe des abscisses (axe des fréquences Doppler), de manière à ce que la corrélation entre les densités spectrales mesurée et attendue soit maximale. La fréquence centrale fDo correspond à la fréquence 15 pour laquelle la densité spectrale théorique attendue est maximale. Lorsque la densité spectrale mesurée est considérée comme polluée elle est étalée sur une bande de fréquence beaucoup plus large que dans le cas " non pollué ». Ce phénomène provient du fait que le spectre des échos utiles, comme par exemple ceux du sol ou de la mer, est séparé 20 fréquentiellement des spectres parasites parce que ces parasites ont une vitesse propre qui fait que leur fréquence Doppler est plus importante. La localisation du centre fDo du spectre utile des échos fixes rétrodiffusés par la zone d'intérêt 50, va être effectuée par corrélation avec la forme du spectre théorique obtenue à partir des caractéristiques du 25 diagramme de l'antenne radar. La méthode va consister à rechercher le coté où la pente est la plus raide sur la courbe du spectre mesuré. Ce coté correspond à la bande de fréquences pour laquelle la densité du spectre mesuré est peu perturbée par les effets parasites. La densité spectrale théorique attendue est ensuite décalée sur l'axe des abscisses (axe des 30 fréquences Doppler), de manière à ce que la corrélation entre les densités attendue et mesurée soit maximale sur le côté de la courbe de spectre mesuré ayant la pente la plus raide. La fréquence centrale fDo est ensuite déterminée en recherchant la fréquence pour laquelle l'allure du spectre théorique est maximum. 35 A titre illustratif les figures 8a et 8b présentent des exemples de représentation graphique de l'amplitude d'un spectre " normal » et de celle d'un spectre " pollué ". La figure 8a représente la superposition du spectre mesuré 81 et de l'allure du spectre utile théorique attendu 82 dans le cas où le spectre mesuré n'est pas pollué par un spectre parasite. Dans ce cas les deux courbes sont semblables. Comme vu précédemment la fréquence centrale fbo correspond à la fréquence pour laquelle l'allure du spectre théorique est maximum. La figure 8b correspond au cas où le spectre est pollué. L'allure de la densité spectrale théorique est placée de façon à faire coïncider le coté où la pente de la courbe de la densité spectrale mesurée est la plus raide. Suite au choix de la fréquence centrale fp° et de l'étendue Doppler Af, 15 un filtrage de ces spectres est effectué dans chaque case distance par un filtre passe bande centré sur fp() et de largeur At Ainsi, contrairement à la méthode classique d'imagerie SAR, au lieu de conserver une bande de fréquence égale à la bande utile à -3dB, le procédé mis en oeuvre par l'invention conserve une bande de fréquence plus 20 réduite égale à l'étendue spectrale Af calculée précédemment. Pour chaque sous-éclairement, une imagette SAR 51 est réalisée selon des techniques connues de l'homme du métier, en utilisant le signal filtré au cour d'une étape Etp4. 25 Le traitement comprend une opération de conformité géométrique permettant de placer chaque imagette 51 dans une référence de cartographie terrestre. L'image SAR finale est ensuite formée par concaténation en juxtaposant les différente imagettes 51 formées. 30 Compte-tenu des contraintes sur le fractionnement des temps d'éclairement, il n'a aucun vide entre les différentes imagettes 51, l'image SAR finale peut donc être formée sans perte d'information. Il se peut qu'il y ait des zones de recouvrement entre imagettes 51 35 consécutives. Dans le cas de recouvrement d'imagettes 51, on retrouve la problématique d'une imagerie à multi-vues. Selon des techniques connues de l'homme du métier, ce problème peut être résolu par exemple, en moyennant les valeurs des pixels de chaque point de l'image SAR finale recouvert par plusieurs imagettes 51. De façon avantageuse, le mode d'imagerie SAR défilant tel que défini permet d'obtenir un meilleur contraste entre des zones de sol de faible et de forte réflectivité radar. Le procédé suivant l'invention peut s'appliquer dans le cas de modes SAR de faible ou de moyenne résolution sur l'axe transverse utilisés typiquement pour identifier des zones sur le sol ou la mer de différentes réflectivités radar. Dans ces cas, il est possible d'utiliser un temps d'intégration cohérente inférieur au temps de présence d'un point M du sol à imager dans le lobe principal de l'antenne pendant l'avancement du porteur radar. REVENDICATIONS1. Procédé d'élimination d'échos parasites en imagerie SAR caractérisé en 5 ce qu'il comprend : - une étape Etp1 de définition des paramètres d'imagerie SAR, - une étape Etp2 de calcul du spectre du signal reçu, dans chaque case distance, d'une zone d'intérêt 51, - une étape Etp3 de filtrage des spectres dans chaque case distance 10 réalisée par un filtre passe bande centré sur la fréquence centrale fp() des échos fixes rétrodiffusés par la zone d'intérêt 50, ladite fréquence centrale fp() étant obtenu par corrélation du spectre mesuré avec la forme théorique du spectre obtenue à partir des caractéristiques du diagramme de l'antenne radar, 15 - une étape Etp4 de formation d'imagette SAR 51 et - une étape Etp5 de concaténation des imagettes SAR 51 pour former l'image SAR finale. 2. Procédé selon la revendication précédente caractérisé en ce qu'en cas 20 de recouvrement de pixels lors de l'étape Etp5 de concaténation, cette dernière comprend une étape de moyennage de la valeur des pixels de chaque point de l'image SAR finale recouvert par plusieurs imagettes 51. 25