Brevet FR2678460A1 — publié le 19921231

I La présente invention concerne un procédé d'enregistrement et de traitement d'une image d'une scène située dans un plan P donné rapporté à deux directions perpendiculaires x, y, ladite image étant fournie par un instrument de prise de vues présentant, dans la direction u (u = x, y) du plan P, une fonction de transfert de modulation HI (f) de fréquence de coupure f Cu, f désignant la fréquence spatiale. L'invention trouve une application particulièrement avanta- geuse dans le domaine de l'observation de la Terre à partir de satellites artif iciels. Les instruments de prise de vues utilisés dans les systèmes satellitaires d'observation de la Terre comportent, soit un détecteur élémentaire dont l'image dans le plan P parcourt la scène à enregistrer à l'aide d'un dispositif de balayage mécanique couplé au mouvement de translation ou de rotation du satellite, soit une barrette linéaire de détecteurs élémentaires, par exemple du type à couplage de charges (CCD), avec balayage en translation dans la direction perpendiculaire à l'axe de ladite barrette (balayage dit "pushbroom"). De façon usuelle, le procédé conforme au préambule comprend une première étape d'enregistrement de l'image à une fréquence d'échantillonnage feu dans la direction u, égale à la fréquence de coupure f Cu dont l'inverse correspond à la résolution spatiale au niveau du plan P d'observation Cette résolution est en général limitée par les dimensions du (des) détecteur(s) élémentaire(s) utilisé(s). Toutefois, du fait que le dispositif de formation d'images de l'instrument de prise de vues, un télescope par exemple, n'est pas lui- même parfait, la fonction de transfert de modulation H u(f) n'est pas égale à 1 sur tout le spectre de fréquence, ce qui conduit à une dégradation de l'image enregistrée de la scène à observer Il est possible de traiter cette image en lui appliquant un filtre de déconvolution Si on appelle P u(f) et S u(f) les spectres de fréquence respectifs dans la direction u de la scène à observer et de l'image enregistrée, on a la relation SU(f) = H u(f) PU(f) Pour restituer Pu (f), il faut donc appliquer à Su(f) un filtre de déconvolu- tion qui dans le domaine fréquentiel est égal à l'inverse de la fonction de transfert de modulation Pu (f) i Su(f) u u H u Hu(f) Ce qui s'exprime dans le domaine spatial par le produit de convolution p(u) = TF_ 1 l I l s(u) HU(f) TF-1 désignant la transformée de Fourier inverse, et p(u) et s(u) étant les fonctions spatiales représentant respectivement la scène et son image. Cependant la théorie du traitement du signal, notamment le théorème de Shannon, montre que seules les fréquences inférieures ou égales à la moitié de la fréquence d'échantillonnage peuvent être restituées par déconvolution, cette limitation étant due au phénomène général de repliement spectral De plus, lors de la déconvolution, les hautes fréquences, comprises entre f eu/2 et feu, viennent perturber les basses fréquences comprises entre O et feu/2. Il a été proposé dans l'article "Repliement spectral lié à l'utilisation de détecteurs constitués de cellules CCD" de G Jacquemod et al., paru dans la revue "Traitement du Signal" Vol 5, N O 1, 1988, de limiter l'influence du repliement spectral par un suréchantillonnage de l'image, la fréquence d'échantillonnage feu étant prise égale à kf Cu, k étant un nombre entier supérieur à 1 Dans ce cas, il est possible de restituer les fréquences au moins jusqu'à la fréquence de coupure f Cu D'autre part, comme au-delà de f la fonction de transfert de modulation est Cu pratiquement nulle, les fréquences supérieures à f Cu sont sans effet sur les basses fréquences inférieures à f Cu. On dispose donc d'un moyen d'augmenter d'un facteur k la résolution dans la direction u d'un instrument de prise de vues donné ou, selon un autre aspect du procédé, de réaliser un instrument de prise de vues ayant les mêmes performances de résolution mais de moindre encombre- ment dans la mesure o la distance focale du dispositif de formation d'images peut être divisée par k tout en conservant la même ouverture numérique, et donc le même rapport signal à bruit Cette dernière possibilité présente un avantage très important dans le cas de l'application spatiale envisagée ci-dessus puisque ne serait-ce qu'un facteur k = 2 permettrait un gain de 8 en volume et 4 à 5 sur la masse de l'instrument de prise de vues, d'o une diminution appréciable de la charge utile du satellite. Il faut cependant signaler que les avantages précités ne sont pas sans inconvénients En effet, que ce soit dans le cas d'une augmentation de la résolution d'un instrument de prise de vues donné, ou dans celui de l'utilisation d'un dispositif de formation d'images de distance focale réduite, l'opération de déconvolution a pour conséquence une diminution du rapport signal à bruit, ceci étant dû au fait que, compte tenu du spectre de bruit B u(f), le signal radiométrique S u(f) de l'image s'écrit SU(f) = H u(f) P (f) + B u(f) L'application du filtre de déconvolution conduit à une scène restituée: P' (f) fu (P)(f) + U u u Hu(f) Hu(f) La fonction de transfert de modulation H u(f) étant au plus égale à 1, on comprend que l'effet du bruit B u(f) a tendance à augmenter. Aussi, le problème technique à résoudre par l'objet de la présente invention est de proposer un procédé conforme au préambule, comprenant, d'une part, une première étape d'enregistrement de ladite image à une fréquence d'échantillonnage feu de l'image dans la direction u égale à kfcu, k étant un nombre entier supérieur à 1, et, d'autre part, une deuxième étape de traitement de l'image enregistrée consistant à appliquer à ladite image un filtre de déconvolution, procédé qui permettrait, lors de l'opération de déconvolution, de conserver l'essentiel des avantages du procédé de suréchantillonnage tout en ayant la possibilité de toujours maîtriser le niveau de bruit final à une valeur acceptable. La solution au problème technique posé consiste, selon la présente invention, en ce que ledit filtre de déconvolution est égal, dans le domaine fréquentiel et pour chaque direction u, au produit de l'inverse de la fonction de transfert de modulation H u(f) par une fonction de troncature Iu(f) ajustable. Ainsi, la scène reconstituée a pour spectre de fréquence P'(f) = I (f) P (f) + I (f) BU(f) u u u u HU(f) On notera que l'effet de la fonction de troncature I u(f) varie en sens contraire selon que l'on considère la qualité image ou le rapport signal à bruit Si l'on désire limiter sensiblement le bruit final, on choisira une fonction de troncature faible, au détriment de la résolution spatiale représentée par Iu(f) Pu(f), ce qui n'est pas très pénalisant dans le cas d'un instrument de prise de vues donné dont le suréchantillonnage augmente la résolution et réduit le signal Par contre, pour un instrument sous- dimensionné dont le suréchantillonnage permet de maintenir la résolution avec une même ouverture et donc le même signal, on pourra choisir une fonction de troncature I U(f) plus importante Bien entendu, selon l'instrument considéré, la fonction de troncature sera choisie et adaptée, en tenant compte des caractéristiques de la mission à accomplir. La description qui va suivre en regard des dessins annexés, donnés à titre d'exemples non limitatifs fera bien comprendre en quoi consiste l'invention et comment elle peut être réalisée. La figure I est un schéma d'un instrument de prise de vues pour la mise en oeuvre du procédé selon l'invention. La figure 2 montre une fonction de transfert de modulation typique de l'instrument de la figure 1. La figure 3 montre une fonction de troncature applicable en combinaison avec la fonction de transfert de modulation de la figure 2. La figure 4 a montre un mode de réalisation du suréchantillon- nage colonne pour une barrette linéaire de détecteurs. La figure 4 b montre un mode de réalisation du suréchantillon- nage ligne pour une barrette linéaire de détecteurs. La figure 5 est un premier montage de barrettes linéaires réalisant un suréchantillonnage ligne d'ordre 2. La figure 6 est un deuxième montage de barrettes linéaires réalisant un suréchantillonnage ligne d'ordre 2. La figure 1 représente de façon schématique un instrument de prise de vues, installé par exemple à bord d'un satellite artificiel de la Terre, et destiné à fournir une image d'une scène située dans un plan P donné, rapporté à deux directions perpendiculaires x, y Dans l'application envisagée, le plan P correspond au sol. L'instrument 100 de prise de vues est constitué, d'une part, d'un dispositif de formation d'images schématisé par une lentille 110, mais qui pourrait être tout aussi bien un télescope à miroirs par exemple, et, d'autre part, d'une barrette linéaire 120 de détecteurs élémentaires 121 du type CCD, placée au foyer F de la lentille 110 Dans la suite de la description et sans que cela ne soit limitatif, les détecteurs élémentaires 121 sont supposés de forme carrée Comme on peut le voir sur la figure 1, le mode de balayage du plan P est du type "pushbroom" dans lequel l'image d'axe x de la barrette linéaire 120 se déplace dans la direction perpendiculaire y avec une vitesse V égale à la vitesse au sol du satellite. L'instrument 100 de prise de vues présente dans chaque direction u = x, y du plan P une fonction de transfert de modulation H U(f) dont une allure typique est donnée sur la figure 2 Cette fonction de transfert de modulation s'annule pour une fréquence de coupure f Cu correspondant à la résolution de l'instrument 100 Généralement, le pouvoir de résolution est limité par la structure discontinue des détecteurs plus que par la qualité du dispositif de formation d'images Si dx représente la taille des détecteurs élémentaires 121, la résolution au sol Au dans la direction u du plan P est donnée par A u = H dx fo H étant l'altitude du satellite et fo la distance focale du dispositif de formation d'images A titre d'exemple, pour un satellite dont l'orbite est située à 800 km de la Terre, une distance focale de 1 m et des détecteurs élémentaires de 7 pm de côté, on atteint une résolution au sol de 5,6 m. La fréquence de coupure fcu est alors sensiblement égale à I fo fcu = x_ Au Hdx L'image de la scène à observer est d'abord enregistrée à une fréquence d'échantillonnage feu dans la direction u égale à kf Cu, k étant un nombre entier supérieur à 1 Dans le mode de réalisation particulier montré à la figure 2, le coefficient de suréchantillonnage k a été pris égal à 2, de sorte que feu vaut f = 2 f 2 eu cu L u Après cette première étape d'enregistrement, l'image enregistrée Su(f) est soumise à une opération de traitement consistant à lui appliquer un filtre de déconvolution égal, dans le domaine fréquentiel et dans chaque direction u, au produit de l'inverse de la fonction de transfert de modulation Hu(f) par une fonction de troncature Iu(f). La figure 3 montre un exemple d'une telle fonction de troncature analytiquement donnée par 1 si f fl Cette courbe est connue sous le nom de fenêtre de Hanning. Comme cela a été expliqué plus haut, la fonction de troncature Iu(f) est ajustable de manière à obtenir, en fonction de la mission envisagée, le meilleur compromis entre résolution spatiale et rapport signal à bruit. Les figures 4 a et 4 b indiquent respectivement comment, dans le cas d'une barrette linéaire de CCD, un suréchantillonnage par 2 en colonne et en ligne peut être obtenu. En ce qui concerne le suréchantillonnage colonne (Figure 4 a), il suffit de disposer d'un temps d'intégration Ti de la barrette égal au temps At qui met la trace au sol du satellite pour parcourir la distance A y/2 soit: Ti = A t A y Y Au 2 V 2 V Pour réaliser le suréchantillonnage ligne (Figure 4 b), il suffit de disposer de deux barrettes élémentaires parallèles, de même pas dx, décalées l'une par rapport à l'autre de la distance dx/2 de sorte que les images de ces deux barrettes dans le plan P soient décalées de t x/2 = A u/2. La figure 5 donne un premier montage de deux barrettes linéaires réalisant un suréchantillonnage ligne d'ordre 2 Les barrettes 120 ' et 120 " sont montées chacune sur une face d'un prisme respectif 210 ', 210 ", lesdits prismes étant accolés de façon à constituer un diviseur optique 200 de ligne Le flux incident 300 se partage à l'interface 220 des deux prismes formant lame semi-réfléchissante Chaque barrette linéaire 120 ' et 120 " ne reçoit donc que la moitié de l'énergie lumineuse provenant de la scène à observer, d'o une limitation du rapport signal à bruit. Le montage représenté sur la figure 6 permet d'obtenir un suréchantillonnage d'ordre 2 avec un signal multiplié par un facteur 2 par rapport au montage précédent Les deux barrettes linéaires 120 ', 120 ", décalées de dx/2, sont réalisées dans un même composant 400, et, de façon à assurer la synchronisation des images enregistrées séparément par chacune des barrettes, le retard à l'enregistrement de l'une des barrettes par rapport à l'autre doit être égal à k'Tiy k' étant un nombre entier non nul, ce qui impose que la distance d séparant les barrettes soit égale à k'dx/2 Dans ce deuxième montage, chaque ligne de détecteurs reçoit l'intégralité du flux lumineux incident, d'o un signal double de celui correspondant au montage de la figure 5 Dans le cas o le bruit photonique est prépondérant, le rapport signal à bruit est multiplié par un facteur égal à 2. Il faut également noter que le dernier montage peut se généraliser au cas d'un suréchantillonnage d'ordre k supérieur à 2 Il suffit de disposer dans un même composant de k barrettes linéaires décalées de dx/k et distantes de k'dx/k Le nombre k est seulement limité par les possibilités technologiques de réaliser un composant à plusieurs barrettes. Enfin, l'invention a été décrite en mettant en oeuvre une pluralité de détecteurs élémentaires assemblés en barrettes, mais il est bien entendu qu'elle pourrait tout aussi bien être réalisée avec un seul détecteur élémentaire Dans ce cas, l'enregistrement de l'image sur le détecteur est effectué par balayage de l'image dudit détecteur élémentaire dans les deux directions du plan P avec suréchantillonnage sur chacun des balayages. REVENDICATIONS 1 Procédé d'enregistrement et de traitement d'une image d'une scène située dans un plan P donné rapporté à deux directions perpendiculaires x, y, ladite image étant fournie par un instrument ( 100) de prise de vues présentant, dans la direction u (u = x, y) du plan P, une fonction de transfert de modulation Hu (f) de fréquence de coupure f Cu, f désignant la fréquence spatiale, procédé comprenant, d'une part, une première étape d'enregistrement de ladite image à une fréquence d'échantillonnage feu de l'image dans la direction u, égale à kf Cu, k étant un nombre entier supérieur à 1 et, d'autre part, une deuxième étape de traitement de l'image enregistrée consistant à appliquer à ladite image un filtre de déconvolution, caractérisé en ce que ledit filtre de déconvolution est égal, dans le domaine fréquentiel et pour chaque direction u, au produit de l'inverse de la fonction de transfert de modulation H u(f) par une fonction de troncature I u(f) ajustable. 2 Procédé selon la revendication 1, caractérisé en ce que, ledit instrument de prise de vues comportant au moins un détecteur élémentaire ( 121) et la fréquence de coupure f Cu étant sensiblement égale à l'inverse de la dimension Au de l'image dudit détecteur élémentaire dans la direction u du plan P, la fréquence d'échantillonnage feu est égale à kl Au. 3 Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que, ledit instrument de prise de vues comportant un détecteur élémentaire ( 121), l'enregistrement de ladite image est réalisé par balayage de l'image dudit détecteur élémentaire dans les deux directions x, y du plan P. 4 Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que, ledit instrument ( 100) de prise de vues comportant au moins une barrette linéaire ( 120) d'une pluralité de détecteurs élémentaires ( 121) dont les images dans le plan P sont alignées dans la direction x, l'enregistrement de ladite image est réalisé par balayage de l'image de ladite barrette linéaire ( 120) dans la direction y. Instrument de prise de vues pour la mise en oeuvre du procédé selon la revendication 4, caractérisé en ce que ledit instrument ( 100) comporte k barrettes linéaires parallèles ( 120, 120 "), de même pas dx, décalées successivement de la distance dx/k et séparées l'une de l'autre de la distance k'dx/k, k' étant un nombre entier non nul. 6 Instrument de prise de vues pour la mise en oeuvre du procédé selon la revendication 4, caractérisé en ce que, le nombre k étant égal à 2, ledit instrument ( 100) comporte au moins deux barrettes linéaires parallèles ( 120 ', 120 ") de même pas dx, disposées sur un diviseur optique de ligne et décalées l'une par rapport à l'autre de la distance dx/2.