Brevet FR2833713A1 — publié le 20030620

La présente invention concerne un procédé d'amélioration de la qualité d'images radar et notamment de la qualité des images obtenues par les radars à ouverture synthétique ou SAR pour"Synthetic Aperture Radar" selon la dénomination anglo-saxonne Les techniques disponibles actuellement permettent d'obtenir des images SAR d'une résolution extrêmement fine-de l'ordre de 50cm à 100Km-mais proposent des images polluées par différents phénomènes parmi lesquels on peut citer : - Les images d'échos mobiles qui présentent un étalement de spectre important sur l'axe doppler. - Les images "fantômes" liées aux phénomènes d'ambiguïté en doppler et en distance qui présentent également un étalement de spectre important sur l'axe doppler et l'axe distance. - Les images doppler parasites issues de l'émergence d'échos forts aux travers de lobes secondaires lointains présents dans le signal imparfaitement focalisé. Ce dernier problème est surtout sensible sur des images de résolution extrêmement fine et n'est pas résolu par les méthodes classiques d'autofocalisation. Pour remédier à ces phénomènes de pollution on utilise généralement des méthodes qui présentent des inconvénients et dont les résultats sont imparfaits : - La suppression des images d'échos mobiles peut être réalisée par extraction de ces échos au moyen de procédés de filtrage spatiotemporel. Ces procédés, de type STAP (SpatioTemporal Adaptative Process), ont notamment pour inconvénients de nécessiter l'utilisation de plusieurs voies de réception et d'apporter une charge de calcul importante. De plus ils n'apportent pas de solution à la pollution liée aux lobes secondaires lointains issus d'échos forts. - L'émergence d'échos liés au niveau de lobes secondaires lointains est habituellement traitée par des méthodes d'autofocalisation dont l'efficacité n'est pas suffisante dans le cas des images de résolution extrêmement fines. Aucune de ces méthodes ne s'avère donc efficace pour supprimer les remontées d'échos de fort niveau a travers les lobes secondaires lointains. Ceci représente un inconvénient majeur dans le cas d'images de résolution extrêmement fine. Un but de l'invention est de pallier les inconvénients précités et en particulier d'apporter une solution à la pollution apportée par les lobes secondaires lointains d'échos de fort niveau. A cet effet, l'invention a pour objet un procédé d'amélioration d'une image radar, comportant une étape de modélisation d'un écho ponctuel en la somme d'un signal focalisé et d'un signal défocalisé, puis une étape de filtrage du signal défocalisé. En particulier pour une case distance D le signal S (m, n) est modélisé comme suit : m D=PL/W-+)] + () expD' (2J+ () +) m m L'étape de filtrage effectue la décomposition spectrale du signal SD (m) et de n signaux obtenus à partir de SD (m) au travers d'une fenêtre glissante. Elle effectue aussi la recherche du signal minimum à partir des représentations spectrales obtenues. Avantageusement, cette décomposition spectrale est réalisée par application d'une FFT au signal SD (m), suivie d'une pondération adaptative destinée à éliminer les lobes secondaires de FFT. Avantageusement le procédé comporte une étape de suréchantillonnage de l'écho reçu. L'invention a notamment pour principaux avantages : - D'apporter une solution globale aux problèmes de pollution liés à la composante du signal résultant des imperfections de l'autofocalisation. - De ne pas nécessiter la mise en place de plusieurs voies de réception et d'être donc peu coûteuse en matériel et en puissance de calcul. D'autres caractéristiques et avantages de l'invention apparaîtront au travers de la description qui suit, faite en regard des dessins annexés qui représentent : La figure 1, l'image SAR d'un écho ponctuel focalisé dans le repère axe doppler 1 axe distance La figure 2, l'illustration du passage du domaine axe temps long/ axe fréquence rapide au domaine axe doppler/axe distance par une FFT suivant deux dimensions. La figure 3, l'exposé des étapes possibles du procédé selon l'invention La figure 4, l'illustration de l'effet de l'étape de filtrage sur un signal ponctuel focalisé modélisé La figure 5, l'illustration de l'effet de l'étape de filtrage sur un cas particulier de signal ponctuel défocalisé modélisé La figure 6, L'illustration des étapes précédentes du procédé associées à une étape de suréchantillonnage et une étape de restitution d'une image doppler en sinus cardinal La figure 7, l'illustration d'une méthode de reconstitution d'un signal temporel extrapolé, à partir du seul lobe principal de son image sur l'axe doppler suréchantillonné. La figure 8, l'illustration de l'application de la méthode précédente au procédé suivant l'invention. La figure 1 présente l'image SAR d'un écho ponctuel immobile dans un plan défini par un axe doppler et un axe distance. Cet écho est caractérisé par sa fréquence doppler f, sa distance r et son amplitude complexe A. On constate que cette image présente un défaut de netteté 11 selon les axes doppler et distance. Dans le cas de la figure 1, cet écho ponctuel immobile étant correctement focalisé, la pollution de l'image est due aux lobes secondaires de FFT. Cette pollution est éliminée par le procédé selon l'invention. La figure 2 illustre une façon d'obtenir cette image. Le signal S (m, n) reçu par le radar, exprimé dans l'espace Temps long 1 Fréquence rapide, est traité par FFT 2D ou transformée de Fourier rapide suivant deux dimensions. L'axe doppler est l'image par FFT de l'axe temps long qui représente la succession des récurrences émises. L'axe distance est pour sa part l'image de l'axe fréquence rapide, qui correspond au spectre du signal reçu pendant une récurrence. Chaque terme de la FFT s'écrit : m n S (m, n) = Aexp (j2' -) exp (j2--) (1) M N Dans l'expression (1) m est la coordonnée du point sur l'axe temps long qui varie de 0 à M-1 par valeur entière ; et n est la coordonnée du point sur l'axe fréquence rapide qui varie de 0 à N-1 également par valeur entière. L'application de cette FFT 2D à un signal ponctuel focalisé donne dans le plan axe doppler 1 axe distance, une image fréquentielle en sinc (sinus cardinal) qui centrée sur une fréquence doppler f et une distance r. Le sinus cardinal est le signal qui correspond à l'image par FFT d'un signal rectangulaire. Cette image fréquentielle a par ailleurs pour défaut de présenter d'importants lobes secondaires de FFT dont on constate l'effet de pollution sur l'image de la figure 1. L'application simple d'une FFT 2D à un écho ponctuel, même bien focalisé, ne donne donc pas une image distance/ doppler ponctuelle. Le procédé selon l'invention complète donc ce traitement par une pondération adaptative qui va fortement atténuer ces lobes secondaires. Si on effectue sur un écho ponctuel focalisé la FFT sur l'axe fréquence rapide, le signal correspondant s'écrit, pour une case distance D : S D (m) = ADexp[j (2ï ( ; + tp D) ] avec mE [0, M -1] (2) M OD)] où Ao est un nombre réel et XD un terme de phase appartenant à l'intervalle [0, 2rc]. Pour un écho ponctuel défocalisé c'est à dire par exemple un écho mobile ou un écho ambigu en doppler ou en distance, la relation (2) est généralement reformulée comme suit : m SD (m) = AD (m) exp [j (2V+ (m) +)]. (3) Ici Ao (m) et (m) sont fonction de m et dépendent de l'écho considéré. Le terme de phase rp (m) constitue le déphasage parasite responsable de la défocalisation de l'écho considéré et des remontées de signal par les lobes secondaires. Selon l'invention, le signal correspondant à la FFT sur l'axe fréquence rapide d'un écho ponctuel est modélisé, pour une case distance D, par la relation suivante : ) = xpD- (2-+] + () exp [7 (2-+) +) (4) SD (M) =AD M +9D)] +AD M (O (M) + (PD) Selon l'invention S D (m) s'écrit donc comme la somme de deux signaux : focal exp [j (2) zf m- (5) (PDW,)] (5) qui correspond à la réponse d'un écho ponctuel focalisé. m RD (m) = xpD- (2'- + (ï) +) (6) qui correspond à la réponse d'un écho défocalisé L'énergie du signal focalisé S focal est de préférence plus forte que celle du D signal défocalisé RD' Dans le cas d'un point focalisé, RD (m) aura une valeur nulle. Les descriptions qui suivent sont basées sur l'exploitation de la modélisation décrite par la relation (4). La figure 3 présente des étapes possibles du procédé selon l'invention. Dans une étape 31 le signal S D (m) est modélisé selon la relation (4) précédente. Dans une étape 32 le signal RD (m), modélisation d'un point defocalisé, est éliminé ou atténué. La figure 4 illustre cette étape 32 par un exemple de filtrage, appliqué sur un signal de type Sfo"" (m) correspondant à un point focalisé. D Dans ce cas RD (m) est nul. Le chronogramme 41 est la représentation sur l'axe temps long des M échantillons du signal analysé S ; cal (m). Le signal étant issu d'un point D focalisé sa période est constante sur l'intervalle des M échantillons. On obtient ensuite par FFT la représentation spectrale 42 de ce signal sur l'axe doppler. Ce spectre est une fonction en sinus cardinal de la fréquence, dont la caractéristique est de présenter des lobes secondaires élevés. On applique alors à ce spectre une pondération adaptative, également appelée"Spatially Variant Apodization", ou"SVA", selon la dénomination anglo-saxone. L'action de cette dernière sur un signal de type sinus cardinal consiste à supprimer les lobes secondaires sans toucher au lobe principal. On obtient ainsi un signal dont la représentation spectrale 43 ne présente plus que le lobe principal. Une méthode de pondération adaptative est notamment décrite dans le document"Spatially Variant Apodization for sidelobe control"de MM. Stankwitz, Dallaire et Finup contenu dans la revue"Proceedings of IEEE 1994 National Radar Conference, Atlanta GA, March 1994". En parallèle, on applique ce traitement par FFT et pondération adaptative à n groupes de M/2 échantillons. Ces groupes sont pris dans une fenêtre glissante 47 parmi l'ensemble des M échantillons, de telle sorte que par glissements successifs la fenêtre ait couvert l'ensemble des M échantillons. On obtient donc n représentations spectrales 44 après FFT, puis n représentations spectrales 45 après pondération adaptative. Ces n représentations spectrales correspondent à n portions du signal SD (m). Du fait du plus faible nombre d'échantillons utilisés, ces représentations spectrales ont une résolution moins bonne que la représentation primaire à M échantillons. On parle de résolution fine 43 et de résolution grossière 45. Ces différentes images sont ensuite utilisées pour une recherche du signal minimum 46. Celui-ci est obtenu en gardant pour chaque fréquence doppler le signal d'amplitude le plus faible de l'ensemble des images utilisées. Dans cette opération on utilise également l'image de résolution fine qui comme on le voit sur la figure, est nécessaire pour maintenir la bonne largeur de lobe. Ceci implique de tenir compte de la différence de gain après filtrage, entre l'image de résolution fine et les autres images. Pour cela il peut être nécessaire de les normaliser entre elles avant la recherche du minimum 46. Le gain obtenu après traitement par FFT, de groupes de M/2 échantillons est inférieur de trois décibels à celui obtenu après traitement de groupes de M échantillons. On multiplie donc par un coefficient deux l'amplitude des spectres des images de résolution grossière. Il apparaît donc que lorsque l'on applique le traitement décrit précédemment à un point focalisé, on peut constater que les n portions de taille M/2 du signal SD (m) ont des périodes constantes et identiques : Leurs spectres d'amplitudes sont identiques et centrés sur la même fréquence. La recherche du minimum restitue donc l'image pleine résolution. La figure 5 illustre l'application de l'exemple de filtrage sur un signal de type RD (m) correspondant à un point défocalisé. Dans ce cas focal (m) est nul. SD Un chronogramme 51 représente, sur l'axe temps long, les M échantillons du signal analysé RD (m). Le signal étant, cette fois-ci, issu d'un point défocalisé, sa période n'est plus constante sur l'intervalle des M échantillons. RD (m) est choisi ici, à titre d'exemple, avec un terme de phase parasite zu (m) de nature quadratique, représentatif par exemple d'un point mobile en translation à vitesse constante. On constate que le chronogramme 51 n'a plus une période constante sur l'intervalle des M échantillons. Après FFT on obtient un spectre 52 quelconque centré sur une fréquence fo auquel on applique la pondération adaptative qui restitue un signal de spectre 53 quelconque. Comme sur la figure 4, on applique en parallèle le traitement, suivi d'une normalisation, à n groupes de M/2 échantillons, pris dans une fenêtre glissante 56 parmi l'ensemble des M échantillons. On obtient donc n représentations spectrales 54, correspondant à n portions du signal RD (m) qui, contrairement à ce qui se passe dans le cas d'un point focalisé, n'ont pas une période constante. Ces n représentations spectrales sont donc centrées sur des fréquences différentes. De même on effectue alors la recherche du signal minimum ce qui a pour effet de supprimer intégralement le signal : le point défocalisé est éliminé 55. Le même processus appliqué au cas où rp (m) est une fonction quelconque de m donne un résultat similaire, à ceci près que l'on a atténuation et non disparition du signal après recherche du minimum. Les figures 4 et 5 illustrent également l'application de l'exemple de filtrage, au cas d'un point quelconque dont le signal correspondant sur l'axe temps long est décrit par la relation (4). Ce signal est exprimé comme la somme d'une composante focalisée et d'une composante défocalisée. Dans ce cas le traitement présenté conduira au maintien de la composante focalisée de l'écho et à la disparition ou au moins à une atténuation sensible de la composante défocalisée. De plus, t'élimination des lobes secondaires de FFT, au cours de l'opération, permet de s'affranchir d'une éventuelle pollution par ces lobes, pollution qui pourrait dégrader l'opération de recherche du minimum. Les étapes 31 et 32, illustrés à la figure 3, permettent avantageusement d'obtenir, dans le système d'axe Doppler/Distance, une image débarrassée des points défocalisés. Cependant, une étape 32 telle que décrite par les figures 4 et 5 ne permet pas le retour par FFT inverse à un signal représenté sur l'axe temps long parS"' (m). En effet ce signal a un spectre en sinus cardinal avec un lobe principal et des lobes secondaires, et non un spectre limité au seul lobe principal. Cette caractéristique peut s'avérer gênante, notamment dans le cas où l'on veut utiliser certains procédés de traitement des images SAR qui s'appliquent normalement à des points focalisés dont la représentation spectrale sur l'axe doppler à la forme d'un sinus cardinal. Le fait de ne pas disposer après filtrage d'un signal de ce type peut altérer l'efficacité du traitement utilisé. Ceci est par exemple le cas des traitements de très grande résolution, communément appelés traitements d'hyper résolution. La figure 6 présente une nouvelle suite d'étapes du procédé selon l'invention pour pallier cet inconvénient. L'étape 31 de la figure 3 est précédée d'une étape 61 de suréchantillonnage sur l'axe doppler, basée sur l'extrapolation du signal SD (m) à au moins 2M échantillons. L'étape 32 de filtrage est par exemple suivie d'une étape 62 de restitution d'une image doppler en sinus cardinal. La figure 7 Illustre un principe basé sur l'introduction d'une étape de suréchantillonnage sur l'axe doppler du spectre du signal ! (). Cette étape est préalable à l'application du traitement décrit précédemment. Elle consiste à effectuer la FFT du signal SD (m) extrapolé sur au moins 2M échantillons 71. L'extrapolation complète les M échantillons initiaux par des échantillons nuls pour obtenir au moins 2M échantillons. Le spectre 72 obtenu après FFT sur 2M points est identique au spectre obtenu précédemment mais échantillonné deux fois plus finement. Après pondération adaptative le signal obtenu se présente sous la forme du lobe principal débarrassé des lobes secondaires. En revanche l'image sur l'axe temps long de ce signal par FFT inverse n'est plus de la forme SD (). C'est un signal extrapolé dont la durée est supérieure à M échantillons. C'est de plus un signal dégradé qui correspond à un signal de la forme SD (m) affecté d'une pondération P (m) dont chaque échantillon SD (i) serait multiplié par un coefficient de pondération P (i). L'expression de la pondération P (m) qui affecte le signal, est donnée par le résultat de l'application à un signal rectangulaire de M échantillons du traitement décrit aux figures 4 et 5. Pour restituer à ce signal sa forme originale, on effectue sa division par P (m). Chaque échantillon SD (i) est divisé par le coefficient P (i) correspondant. Le résultat de cette opération est un signal formé des M échantillons du signal d'origine auxquels s'adjoignent des échantillons latéraux extrapolés. La figure 8 illustre l'application de la méthode décrite au paragraphe précédent. Après être passé par l'étape 61 de la figure 6 et avoir subi une extrapolation, le signal SD (m) est traité par FFT, pondération adaptative et recherche du minimum. Il est alors représenté par le signal 81 sur l'axe doppler suréchantillonné. On lui applique ensuite une FFT inverse qui donne le signal temporel extrapolé et pondéré 74 présenté figure 7 dont on ne retient que les M échantillons centraux. Après application de la pondération inverse de P aux M échantillons retenus on obtient la restitution 82 d'un signal S focal (m) dont D la représentation spectrale sur l'axe doppler a une expression en sinus cardinal. Ces dernières opérations constituent l'étape 62 de la figure 6. Associée au procédé selon l'invention cette méthode permet donc avantageusement de récupérer à partir d'un signal SD (m) quelconque l'image spectrale en sinus cardinal du signal S focal (m) correspondant D Le procédé selon l'invention ainsi que la méthode de restitution qui lui est associée, ont été décrits, à titre d'exemple, avec l'emploi d'un traitement de type FFT, complété par une pondération adaptative. Ce traitement ayant pour but d'effectuer le passage du domaine axe temps long /axe fréquence rapide au domaine axe doppler/axe distance, d'autres types de traitements sont bien sûr possible, pour autant qu'ils aboutissent à la représentation spectrale des signaux. L'étape de filtrage 32 a été décrite avec l'emploi d'un traitement par FFT, suivi d'une pondération adaptative des représentations spectrales. Ce filtrage est bien adapté à un traitement numérique. D'autres solutions peuvent bien sûr être envisagées, notamment avec des filtres analogiques. L'invention a notamment comme avantage de ne nécessiter qu'une seule voie de réception, contrairement aux méthodes de filtrage spatial qui nécessitent au moins, pour l'analyse du signal, une voie Somme et une voie Différence. Corollairement, elle présente comme avantage de nécessiter une charge de calcul moins importante. REVENDICATIONS 1. Procédé d'amélioration d'une image radar caractérisé en ce qu'il comporte au moins une étape (31) de modélisation d'un écho ponctuel (S (m, n)) en un signal focalisé (S focal (m et un signal défocalisé (RD (m)), D puis une deuxième étape (32) de filtrage du signal défocalisé. 2. Procédé selon la revendication 1, caractérisé en ce que l'écho S (m, n) est modélisé, pour une case distance D, selon la relation suivante : SD =exp[j (2+)] + AD (m) exp [7 (2+) +) m m 3. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'étape (32) de filtrage effectue la décomposition spectrale du signal SD (m) et de n signaux obtenus à partir de SD (m) au travers d'une fenêtre glissante (47), effectue le filtrage des lobes secondaires des représentations spectrales obtenues (42 et 44) puis effectue la recherche du signal minimum (46) à partir des représentations spectrales après filtrage (43 et 45). 4. procédé selon la revendication 3, caractérisé en ce que, la décomposition spectrale étant effectuée par FFT, les lobes secondaires de FFT sont éliminés par pondération adaptative. 5. Procédé selon l'une quelconque des revendications 3 ou 4, caractérisé en ce que les n représentations spectrales (44) des signaux obtenus à partir deSD (m) au travers de la fenêtre glissante (47), font l'objet d'une normalisation par rapport à la représentation spectrale de S D (m) (42), normalisation consistant en une multiplication par deux de leur niveau. 6. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte une étape (61) de suréchantillonnage de la représentation spectrale de l'écho sur l'axe doppler, qui précède l'étape (31) de modélisation. 7. Procédé selon la revendication 7, caractérisé en ce que le suréchantillonnage en doppler est réalisé à partir du signal SD (), extrapolé à au moins 2M échantillons par ajout d'au moins M échantillons nuls au M échantillons constituants. 8. Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce qu'il comporte une étape (62) de reconstitution d'un écho ponctuel focalisé (82) à partir de l'image spectrale (81)