Brevet FR2694412A1 — publié le 19940204

PROCEDE ET DISPOSITIF DE TRAITEMENT DES SIGNAUX ACQUIS PAR UN SYSTEME DE CARTOGRAPHIE La présente invention concerne un procédé et un disposi- tif de traitement des signaux acquis par un système de cartographie. Bien que nous ne décrirons par la suite que l'application à un radar de cartographie, il est d'ores et déjà bien entendu que le terme "système" s'applique également au sonar de cartographie. Les radars de cartographie sont soit aéroportés, soit embarqués sur satellites, et leur mission est d'acquérir des signaux provenant de la zone que l'on cherche à étudier en envoyant sur cette zone des impulsions à une fréquence de récurrence FR donnée Ces radars sont généralement à visée latérale On utilise alors la fréquence Doppler des signaux acquis pour déterminer la distance et la position des points de la zone. On sait à l'heure actuelle que l'utilisation d'un radar à compression d'impulsions permet d'obtenir une bonne résolution en distance On connaît des radars SAR à visée latérale, SAR étant les initiales anglo- saxonnes de Synthetic Aperture Radar (Radar à ouverture synthétique), dont la particularité en mode classique de prise de vue, est d'avoir toujours le même faisceau d'éclairement en azimut au fur et à mesure que le satellite ou l'avion avance sur sa trajectoire. L'inconvénient de ces radars est que la bande Doppler Bd occupée par les signaux échos est égale à la fréquence de récurrence FR du radar Or, la résolution en azimut est proportionnelle à la quantité vx, v étant la vitesse de déplacement du radar et t, la résolution temporelle, X étant elle-même inversement proportionnelle à la bande Doppler Bd occupée par les signaux échos. La nécessité d'augmenter cette bande Doppler Bd en vue d'avoir une meilleure résolution en azimut a donné naissance au ra- dar SAR utilisant un mode télescope, plus connu sous le nom em- prunté à l'anglais de mode spotlight, dont on peut trouver une description détaillée dans le document "SPOT MODE SAR FOCUSING WITH THE W-K TECHNIQUE" de C Prati, A Monti Guarnieri et F Rocca (IGARSS 91, pages 631, 634) Le principe général du mode spotlight est d'éclairer toujours la même zone au fur et à mesure que le radar avance, et par conséquent, d'incliner le faisceau du radar au fur et à mesure de son déplacement de telle sorte que cette zone soit tou- jours éclairée L'acquisition des signaux échos ne commence que lorsque ce faisceau intercepte la zone complète. Avec ce mode, la bande Doppler Bd qui correspond à l'intervalle de décalage en fréquence Doppler pendant le temps d'éclairement de la zone, est beaucoup plus grande que la bande Doppler l FR/2; + FR/21 du radar, et, consécutivement, la résolu- tion en azimut est meilleure Cependant, l'inconvénient majeur du mode spotlight est que la longueur de la zone imagée se réduit à la largeur à 3 d B du lobe de l'antenne radar. Pour voir une zone plus étendue qu'en mode spotlight, un mode spotlight amélioré, que nous appellerons mode spotlight étendu pour le distinguer du mode spotlight classique vu précédem- ment, est actuellement à l'étude Le principe du mode spotlight étendu est illustré schématiquement sur la figure 1 sur laquelle l'axe XX' représente la trajectoire du radar à vitesse sensiblement constante, et la droite passant par les trois points E 1, E 2, E 3, une portion azimutale de la zone explorée, parallèle à l'axe XX'. L'acquisition des signaux échos selon le mode spotlight étendu s'effectue dès que le faisceau du radar intercepte un point de la zone Comme illustré sur la figure 1, le faisceau id à l'instant td vient de terminer l'acquisition des échos provenant du point E 1, et commence juste l'acquisition des échos provenant du point E 2 A l'instant tc, le faisceau lc du radar éclaire symétriquement de part et d'autre du point E 2 Cet instant tc correspond à la distance minimale du radar au point E 2 et à une fréquence Doppler nulle du signal écho provenant de E 2 Enfin, à l'instant tf, le faisceau if du radar correspond à la fin de l'acquisition des échos provenant de E 2, et au début de l'acquisition de ceux provenant de E 3. Dans ce mode spotlight étendu, la longueur de la zone imagée dépend de la capacité de dépointage de l'antenne en azimut. Il est également nécessaire que la rotation du faisceau en azimut as- sure la couverture désirée en respectant le temps d'éclairement im- posé par la résolution azimutale Tout comme dans le cas du mode spotlight classique, la bande Doppler Bd est plus grande que la bande Doppler l FR/2; + FR/2 l du radar. L'avantage est, comme nous l'avons déjà dit, une amélioration de la résolution en azimut, mais il existe cependant un inconvénient consistant en un phénomène de repliements des fré- quences Doppler dans le plan temps-fréquence Doppler, ces replie- ments étant sources d'ambiguïtés Par conséquent, il est impossible d'appliquer un traitement SAR standard aux signaux échos ainsi ac- quis On sait résoudre actuellement ce problème pour le mode spot- light classique Une solution est par exemple exposée dans le docu- ment "A TOMOGRAPHIC FORMULATION OF SPOTLIGHT MODE SYNTHETIC APERTURE RADAR" de Munson D C, J D O'Brien et W K Jenkins (IEEE, Volume 71, pages 917-925 Août 1983). Par contre, aucune solution n'a été proposée pour le mode spotlight étendu, et c'est ce que la présente invention se pro- pose de faire. Plus précisément, un premier objet de l'invention concerne un procédé de traitement de signaux échos acquis par un système de cartographie de fréquence de récurrence FR, le système fonctionnant en mode spotlight étendu à vision latérale, pour lequel chaque signal écho acquis suit théoriquement dans le plan temps- fréquence Doppler, une loi continue quasi-linéaire et comprise entre deux courbes sensiblement linéaires possédant une même inclinai- son initiale dans ledit plan et définissant une bande oblique, ladite loi étant en pratique discontinue en raison de repliements des fré- quences Doppler dans la bande utile du système l FR/2; + FR/2 l, le procédé étant caractérisé en ce qu'il consiste à effectuer successivement les étapes suivantes: faire un sur-échantillonnage des signaux échos acquis, à une fréquence F proportionnelle à la fréquence de récurrence FR, avec un coefficient de proportionnalité choisi de façon à ce que la- dite loi redevienne continue dans la bande oblique; rendre la bande oblique horizontale dans ledit plan de fa- çon à la superposer à la bande utile l FR/2; + FR/2 l du radar faire un filtrage passe-bande de bande l-FR/2; FR/2 l rendre à la bande oblique son inclinaison initiale. Un second objet de l'invention concerne un dispositif de traitement pour la mise en oeuvre du procédé selon les revendica- tions 1 à 5, caractérisé en ce qu'il comporte: un dispositif de suréchantillonnage à la fréquence F, recevant en entrée les signaux échos acquis par le radar; un premier module de génération d'une rampe de fréquence de pente opposée à l'inclinaison initiale de la bande oblique; un premier module de multiplication effectuant la multiplication des échantillons issus du dispositif de sur- échantillonnage par la rampe de fréquence générée par le premier module de génération d'une rampe; un filtre passe-bande de bande l-FR/2; + FR/2 l filtrant les échantillons issus du premier module de multiplication; un second module de génération d'une rampe de fréquence de pente égale à l'inclinaison initiale de la bande oblique; un second module de multiplication effectuant la multiplication des échantillons issus du filtre passe-bande par la rampe de fréquence générée par le second module de génération d'une rampe. L'invention, ainsi que ses avantages, seront mieux compris au vu de la description suivante, faite en référence aux figures annexées: La figure 1 illustre le principe de l'acquisition des signaux échos en mode spotlight étendu; La figure 2 montre, dans le plan temps-fréquence Doppler, les lois suivies théoriquement par les signaux échos acquis provenant des points E 1, E 2, E 3 de la figure 1, et ce que voit le radar en pratique; La figure 3 donne, dans le plan temps-fréquence Doppler, l'allure des signaux à l'issue de la première étape du procédé de traitement selon l'invention; La figure 4 montre, dans ce même plan, l'allure des signaux après la seconde étape du procédé de traitement selon l'invention; La figure 5 illustre, dans le même plan, l'allure des signaux après la dernière étape du procédé selon l'invention; La figure 6 illustre les courbes de migrations, dans le plan distance-fréquence Doppler, relatives aux trois signaux; La figure 7 est un exemple de réalisation d'un dispositif pour la mise en oeuvre du procédé selon l'invention. La figure 1, illustrant le principe de l'acquisition en mode spotlight étendu, a déjà été décrite dans l'introduction On supposera par la suite que le radar a commencé à émettre à un temps t' antérieur à td et correspondant au début de l'acquisition pour le point E 1, et continue d'émettre jusqu'à un temps t" correspondant à la fin de l'acquisition pour le point E 3. Avec les hypothèses précédentes, on a représenté simultanément sur la figure 2 la loi dans le plan temps-fréquence Doppler suivie théoriquement par les signaux échos acquis et le phénomène de repliement des fréquences dû en pratique à la fréquence de récurrence FR du radar. Compte tenu de la vitesse sensiblement constante de déplacement du radar, du mode d'acquisition qui se déclenche dès qu'un point Ei est intercepté par le faisceau du radar, et de la vitesse de rotation sensiblement constante du faisceau au fur et à mesure que le radar se déplace, il est facile de démontrer que les signaux acquis devraient avoir chacun une loi en fréquence Doppler quasi-linéaire en fonction du temps et de même pente On a représenté ainsi les trois droites continues théoriques référencées 1, 2 et 3 relatives respectivement aux signaux échos en provenance des points E 1, E 2, et E 3, comprises respectivement dans les intervalles de temps lt', tdl, ltd, tfl et ltf, t"l. Pour faciliter la compréhension, on a pris comme conventions: de représenter en trait continu tout ce qui se situe dans la bande utile l FR/2; + FR/2 l du radar, et en traits pointillés, tout ce qui est extérieur à cette bande utile; de distinguer les portions des trois lois relatives aux trois points E 1, E 2, E 3, d'une part, par des traits (continus ou pointillés) d'épaisseurs différentes, et d'autre part, par des références respectives 1, 2 et 3 Ces conventions valent également pour les figures 3 et 4. La propriété caractéristique du mode spotlight étendu est que les lois précédentes sont comprises entre deux courbes 4 a, 4 b, quasi-linéaires, définissant une bande oblique de largeur sensiblement égale à la fréquence de récurrence FR du radar. L'inclinaison de cette bande oblique est causée par le fait que les points E 1, E 2 et E 3 apparaissent sous un angle de dépointage différent à mesure que le radar avance Cette inclinaison est aussi fonction de la vitesse de rotation du faisceau du radar. s Cependant, ces trois droites comprises dans la bande oblique ne sont pas obtenues en pratique si l'on n'effectue aucun traitement des signaux acquis En effet, la fréquence de récurrence FR du radar réalise un échantillonnage naturel des signaux, et il en résulte un repliement des portions de droites en traits pointillés à l'intérieur de la bande utile du radar Ainsi, comme le montre la figure 2, on obtient dans la bande de fréquence l FR/2; + FR/2 l trois portions de droites discontinues pour chaque droite 1, 2, 3, les portions correspondant aux repliements constituant des cibles fantômes. Le procédé de traitement selon l'invention propose une série d'étapes qui permet de reconstituer les trois droites théoriques comme nous allons le décrire ci-après: La première étape du procédé selon l'invention consiste à faire un sur-échantillonnage des signaux acquis à une fréquence F proportionnelle à la fréquence FR de récurrence Le coefficient de proportionnalité est choisi de façon à reconstituer les droites 1, 2 et 3 dans la bande oblique Une fréquence F égale à 4 FR suffit pour le cas de figure traité. Dans le cas général, si la fréquence F est trop importante, on pourra partitionner le traitement en blocs le long de l'axe azimutal, effectuer le traitement pour chaque bloc, puis concaténer les résultats. La figure 3 montre le résultat obtenu dans le plan temps- fréquence à l'issue de cette étape de sur-échantillonnage On peut constater sur cette figure 3 que la bande oblique ( 4 a, 4 b), dans l'intervalle de temps allant de t' à t", ne contient que trois droites continues 1, 2, 3 qui sont les trois lois reconstituées Dans la pratique, on peut procéder au sur-échantillonnage de différentes manières: Si le coefficient de proportionnalité est un nombre entier N, on peut intercaler, entre chaque signal écho acquis (N-1) zéros. Une autre méthode consiste à faire la transformée de Fourier des signaux échos acquis, de rendre périodique de période FR cette transformée sur un intervalle de fréquence l-F/2; +F/2 l, puis de revenir en temporel en prenant la transformée de Fourier inverse Le problème est qu'il reste encore beaucoup d'échos parasites à l'extérieur de cette bande oblique. La seconde étape du procédé selon l'invention consiste à rendre la bande oblique horizontale de façon à la superposer à la bande utile l FR/2; + FR/2 l du radar Pour ce faire, on peut multiplier les échantillons de signaux issus de l'étape de sur- échantillonnage par un signal dont la loi en fréquence en fonction du temps est une droite passant par l'origine du plan temps-fréquence Doppler, c'est-à-dire le point de coordonnées (tc, o), et dont la pente est opposée à celle de la bande utile Le résultat de cette deuxième étape est illustré sur la figure 4: tous les signaux utiles 1, 2, 3 en traits continus se trouvent à l'intérieur de la bande l FR/2; + FR/2 l, entre les deux droites 4 a, 4 b devenues horizontales, et, à l'intérieur de cette bande, les signaux suivent bien une loi quasi- linéaire et continue. L'étape suivante selon l'invention consiste en un filtrage passe-bande de bande égale à l-FR/2; +FR/2 l de manière à éliminer tous les signaux parasites 1, 2, 3 dont seulement quelques uns ont été représentés en pointillés sur la figure 4. Finalement, la dernière étape du procédé selon l'invention consiste à rendre à la bande comprise entre les droites 4 a, 4 b, son inclinaison initiale Ceci peut se faire en multipliant les échantillons après filtrage par un signal dont la loi en fréquence en fonction du temps est une droite passant par l'origine du plan temps-fréquence Doppler, et dont la pente est égale à l'inclinaison initiale de la bande oblique Le résultat de cette étape est illustré sur la figure 5: on peut voir que les lois théoriques 1, 2, 3 ont été reconstituées dans la pratique et qu'il n'y a plus aucune ambiguïté possible Le spectre des signaux de départ a été déplié et les cibles fantômes éliminées. La figure 6 donne les courbes 1 ', 2 ', 3 ' de migration des points respectifs E 1, E 2, E 3 dans le plan distance radar-fréquence. Ces trois courbes 1 ', 2 ', 3 ' sont quasiment superposées et possèdent l'allure générale de ce que l'on obtiendrait avec un radar SAR à visée latérale classique Les échantillons issus de traitement selon le procédé de l'invention peuvent donc être traités par tout processeur SAR classique. La figure 7 donne un exemple non limitatif d'un dispositif pour la mise en oeuvre du procédé selon l'invention: Le signal radar vidéo est tout d'abord traité par un dispositif de sur-échantillonnage 11 à la fréquence F Les échantillons issus du dispositif de sur- échantillonnage 11 sont ensuite multipliés, au moyen d'un premier module de multiplication 22, par des échantillons issus d'un premier module de génération de rampe 33, lequel génère un signal modulé quasi- linéairement en fréquence avec une pente opposée à l'inclinaison initiale de la bande oblique Les échantillons résultant de cette multiplication sont ensuite filtrés par un filtre passe-bande 44 de bande l-FR/2; + FR/21 puis sont multipliés, au moyen d'un second module de multiplication 55 par des échantillons issus d'un second module de génération de rampe 66, lequel génère un signal modulé quasi- linéairement en fréquence avec une pente égale à l'inclinaison initiale de la bande oblique Ceci constitue la fin du traitement selon l'invention, et les échantillons issus du second module de multiplication 55 peuvent être traités par un processeur SAR classique 77. Un grand avantage de la présente invention réside ainsi dans la capacité à transformer des données fournissant une image radar à résolution azimutale plus élevée qu'en SAR classique, sur une zone plus étendue qu'en mode spotlight classique, en des5 échantillons susceptibles d'être traités par un processeur SAR classique conçu pour la visée latérale. REVENDICATIONS 1 Procédé de traitement de signaux échos acquis par un système de cartographie de fréquence de récurrence FR, le système fonctionnant en mode spotlight étendu à vision latérale, pour lequel chaque signal écho acquis suit théoriquement dans le plan temps- fréquence Doppler, une loi continue quasi-linéaire et comprise entre deux courbes sensiblement linéaires possédant une même inclinaison initiale dans ledit plan et définissant une bande oblique, ladite loi étant en pratique discontinue en raison de repliements des fréquences Doppler dans la bande utile du système l FR/2; + FR/21, le procédé étant caractérisé en ce qu'il consiste à effectuer successivement les étapes suivantes faire un sur-échantillonnage des signaux échos acquis, à une fréquence F proportionnelle à la fréquence de récurrence FR, avec un coefficient de proportionnalité choisi de façon à ce que ladite loi redevienne continue dans la bande oblique; rendre la bande oblique horizontale dans ledit plan de façon à la superposer à la bande utile l FR/2; + FR/21 du radar; faire un filtrage passe-bande de bande l FR/2; + FR/2 l; rendre à la bande oblique son inclinaison initiale. 2 Procédé de traitement selon la revendication 1, caractérisé en ce que le coefficient de proportionnalité est un nombre entier N et en ce que l'étape de sur-échantillonnage consiste à intercaler entre des signaux échos acquis adjacents N-1 zéros. 3 Procédé de traitement selon la revendication 1, caractérisé en ce que l'étape de sur-échantillonnage consiste à faire la transformée de Fourier des signaux échos acquis, de la rendre périodique de période FR sur un intervalle de fréquence l-F/2 +F/21, puis de faire la transformée de Fourier inverse. 4 Procédé de traitement selon l'une quelconque des revendications précédentes, caractérisé en ce que, pour rendre la bande oblique horizontale, on multiplie les échantillons issus de l'étape de sur-échantillonnage par un signal dont la loi, dans le plan temps-fréquence Doppler, est une droite quasi-linéaire passant par l'origine du plan de pente opposée à l'inclinaison initiale de la bande oblique. Procédé de traitement selon l'une quelconque des revendications précédentes, caractérisé en ce que, pour rendre à la bande oblique son inclinaison initiale, on multiplie les échantillons restants après ledit filtrage par un signal dont la loi, dans le plan temps-fréquence Doppler est une droite quasi-linéaire passant par l'origine et de pente égale à l'inclinaison initiale de la bande oblique. 6 Dispositif de traitement pour la mise en oeuvre du procédé selon les revendications 1 à 5, caractérisé en ce qu'il comporte: un dispositif de sur-échantillonnage ( 11) à la fréquence F, recevant en entrée les signaux échos acquis par le radar; un premier module de génération d'une rampe de fréquence ( 33) de pente opposée à l'inclinaison initiale de la bande oblique; un premier module de multiplication ( 22) effectuant la multiplication des échantillons issus du dispositif de sur- échantillonnage ( 11) par la rampe de fréquence générée par le premier module de génération d'une rampe ( 33); un filtre passe-bande ( 44) de bande l FR/2; + FR/21 filtrant les échantillons issus du premier module de multiplication ( 22); un second module de génération d'une rampe de fréquence ( 66) de pente égale à l'inclinaison initiale de la bande oblique; un second module de multiplication ( 55) effectuant la multiplication des échantillons issus du filtre passe-bas ( 44) par la rampe de fréquence générée par le second module de génération d'une rampe ( 66).