FR 3015740 A1 20150626 FR 1363273 A 20131220 Procédé, dispositif et système de traitement de correction d'images La présente invention concerne un procédé de traitement de correction d'images acquises par un dispositif d'acquisition d'images en déplacement rectiligne relatif par rapport à une zone observée. Elle concerne également un dispositif de traitement de correction d'images et un système de traitement de correction d'images associés. La présente invention se situe dans le domaine de l'amélioration de la qualité d'images numériques, et trouve des applications notamment dans le domaine de la télédétection d'images, de l'imagerie par satellite et de l'inspection industrielle. Dans le domaine de la télédétection, il est connu d'utiliser, pour la capture d'images, des matrices de capteurs élémentaires, qui peuvent être bi-dimensionnelles ou mono-dimensionnelles, également appelées barrettes, qui sont attachées à un satellite et qui défilent au droit de la zone observée, perpendiculairement à la direction de déplacement du satellite. Ainsi, chaque capteur élémentaire défile au regard du paysage observé, chaque ligne de l'image numérique acquise correspondant à l'acquisition effectuée à un instant donné. L'image numérique obtenue est composée d'une ou plusieurs matrices bidimensionnelles d'échantillons appelés pixels, chaque échantillon ayant une valeur de radiométrie associée. Dans ce mode d'acquisition, chaque colonne d'une matrice d'image numérique obtenue est associée physiquement à un capteur élémentaire, chaque pixel d'une ligne d'image correspondant à une acquisition de signal de ce capteur élémentaire. Chaque matrice d'images correspond à un échantillonnage d'un signal photométrique sur une grille d'échantillonnage bidimensionnelle, de pas d'échantillonnage P, P étant la distance entre barycentres photométriques théoriques de deux pixels successifs dans chaque direction. Ce mode d'acquisition par déplacement rectiligne d'une matrice de capteurs est connu sous le nom d'acquisition par balayage, ou « push-broom » en anglais. Le dispositif d'acquisition d'images est en déplacement relatif rectiligne par rapport à la zone observée. Les capteurs élémentaires intègrent le rayonnement électromagnétique d'une surface observée dans une bande spectrale donnée et le transforment en charges électriques. Chaque capteur est caractérisé par une fonction d'étalement appelée PSF (pour « point spread function ») qui a un module dans le domaine de Fourier noté FTM (pour Fonction de transfert de Modulation »). En cas de saturation, une accumulation de charges électriques excédentaires est constatée au niveau des capteurs élémentaires. Afin d'évacuer les charges excédentaires, il a été proposé d'utiliser des dispositifs d'anti-éblouissement (en anglais « anti-blooming drain »), positionnés latéralement par rapport aux capteurs élémentaires, de manière alternée, tous les deux capteurs élémentaires d'une barrette. Un tel système de placement des dispositifs d'anti- éblouissement latéral est connu sous le nom de système d'anti-éblouissement latéral alterné. Ainsi, les charges excédentaires sont éliminées pour chaque pixel de l'image acquise, chaque pixel étant associé à un capteur élémentaire. Cependant, en l'absence de saturation, la présence d'un tel système d'anti- éblouissement latéral alterné a un impact sur l'accumulation de charges électriques dans un puits de potentiel correspondant au capteur élémentaire, ce qui a pour effet de décaler le barycentre photométrique de chaque pixel par rapport à la grille d'échantillonnage idéale. En effet, si chaque dispositif d'anti-éblouissement latéral est implanté sur la face d'incidence des photons sur les capteurs élémentaires, un tel dispositif d'anti- éblouissement introduit une occultation partielle ou totale de la surface photosensible. Si chaque dispositif d'anti-éblouissement latéral est implanté sur la face opposée à la face d'incidence des photons, il perturbe les champs électriques modifiant la répartition des flux d'électrons créés par les photons incidents. En dehors des dispositifs anti-éblouissement latéral, dans les capteurs de type CMOS (pour « complementary metal-oxide semi-conductor »), diverses fonctions électroniques sont implantées à proximité des capteurs élémentaires, par exemple l'implantation de bus électrique entre deux capteurs. Dans ce cas également, des surfaces photosensibles se trouvent occultées ou décalées, et par conséquent le barycentre photométrique d'au moins certains pixels est décalé par rapport à une grille d'échantillonnage idéale. L'invention a pour objectif d'améliorer la qualité des images numériques acquises dans ces conditions d'acquisition aboutissant à un décalage du barycentre photométrique d'au moins une partie des pixels acquis par rapport à une grille d'échantillonnage initiale qui est la grille d'échantillonnage idéale. A cet effet, l'invention propose, selon un premier aspect, un procédé de traitement de correction d'image numérique acquise par un dispositif d'acquisition d'images comportant au moins une matrice de capteurs élémentaires, en déplacement relatif par rapport à une zone observée, une image acquise étant représentée par au moins une matrice de pixels correspondant à un signal d'image acquis par les capteurs élémentaires échantillonné sur une grille d'échantillonnage initiale, ayant un pas d'échantillonnage horizontal initial et un pas d'échantillonnage vertical initial, chaque pixel ayant un barycentre photométrique associé, le dispositif d'acquisition d'images générant, dans au moins une direction d'irrégularité parmi la direction horizontale et la direction verticale, un décalage des barycentres photométriques par rapport à la grille d'échantillonnage initiale. Le procédé selon l'invention comprend les étapes, appliquées à au moins une portion d'une dite matrice de pixels, de : - extraction de ladite portion d'un premier sous-ensemble de pixels correspondant au signal d'image échantillonné sur une première sous-grille d'échantillonnage et d'un deuxième sous-ensemble de pixels correspondant au signal d'image échantillonné sur une deuxième sous-grille d'échantillonnage, chacune desdites première et deuxième sous-grille d'échantillonnage ayant un pas d'échantillonnage selon la direction d'irrégularité égal au double du pas d'échantillonnage correspondant de la grille d'échantillonnage initiale ; - obtention d'un premier signal d'image à partir dudit premier sous-ensemble de pixels et d'un deuxième signal d'image à partir dudit deuxième sous-ensemble de pixels, - application d'un premier filtre de ré-échantillonnage sur le premier signal d'image et d'un deuxième filtre de ré-échantillonnage sur le deuxième signal d'image, - ré-échantillonnage des signaux d'image filtrés pour obtenir un signal d'image ré-échantillonné sur une grille d'échantillonnage finale ayant un pas d'échantillonnage horizontal égal au pas d'échantillonnage horizontal initial et un pas d'échantillonnage vertical égal au pas d'échantillonnage vertical initial, lesdits premier et deuxième filtre de ré-échantillonnage étant déterminés en fonction dudit décalage des barycentres photométriques par rapport à la grille d'échantillonnage initiale et d'un positionnement souhaité de la grille de ré- échantillonnage finale par rapport à la grille d'échantillonnage initiale. Avantageusement, le procédé de traitement de correction selon l'invention permet d'améliorer la qualité des images numériques acquises par ré-échantillonnage à partir de sous-grilles régulières, permettant de compenser les irrégularités induites par la présence d'un dispositif d'acquisition d'image, en déplacement relatif par rapport à une zone observée, induisant un décalage de barycentre photométrique des pixels. Le procédé selon l'invention peut présenter une ou plusieurs des caractéristiques ci-dessous, prises indépendamment ou selon toutes combinaisons techniquement réalisables. L'image acquise est acquise par une matrice monodimensionnelle de capteurs élémentaires, et ledit premier sous-ensemble de pixels comprend tous les pixels de la portion de matrice de pixels appartenant aux colonnes ou lignes d'indice pair et ledit deuxième sous-ensemble de pixels comprend tous les pixels de la portion d'image traitée appartenant aux colonnes ou lignes d'indice impair de ladite portion de matrice traitée. Le dispositif d'acquisition d'images induit des zones d'occultation de l'acquisition de charges par les capteurs élémentaires auxquels ils sont associés, disposées tous les deux pixels selon la direction d'irrégularité, la détermination desdits premier et deuxième filtre étant fonction de la largeur d'une dite zone d'occultation. Les premier et deuxième filtre sont définis à partir d'un filtre de ré-échantillonnage G par les formules suivantes : Gi (x) = G(x + di ) Et G2 (X) = G(C12 + 1 - 4 . où dl est une distance, selon la direction d'irrégularité, entre les barycentres photométriques de la première sous-grille d'échantillonnage et la grille d'échantillonnage finale, et d2 est une distance, selon la direction d'irrégularité, entre les barycentres photométriques de la deuxième sous-grille d'échantillonnage et la grille d'échantillonnage finale. La détermination des premier et deuxième filtre est effectuée à partir d'un filtre de ré-échantillonnage défini par : G(x) = sinc(zy) - 7r - cotan( 71- - 1+a x sinc2( 71- .x - - où sinc est la fonction sinus 2 , 2 2, cardinal, cotan est la fonction co-tangente, et a est égal à la demi-largeur d'une dite zone d'occultation, x et a étant exprimés en pixels de la grille d'échantillonnage initiale. La détermination des premier et deuxième filtre comprend une détermination de longueur de support de chacun desdits filtres et l'application d'une fenêtre d'apodisation. La détermination desdits premier et deuxième filtre comprend une normalisation des coefficients de chacun desdits premier et deuxième filtre. Le ré-échantillonnage sur la grille d'échantillonnage finale comporte une addition pixel par pixel des signaux d'image filtrés pour obtenir une valeur de pixel du signal d'image ré-échantillonné. L'obtention d'un premier signal d'image consiste à insérer une colonne ou ligne de zéros entre deux colonnes ou lignes consécutives dudit premier sous-ensemble de pixels, et l'obtention d'un deuxième signal d'image consiste à insérer une colonne ou ligne de zéros entre deux colonnes ou lignes consécutives dudit deuxième sous-ensemble de pixels. L'invention concerne, selon un deuxième aspect, un dispositif de traitement de correction d'image numérique acquise par un dispositif d'acquisition d'images en déplacement rectiligne relatif par rapport à une zone observée, comportant au moins une matrice de capteurs élémentaires, une image acquise étant représentée par au moins une matrice de pixels correspondant à un signal d'image acquis par les capteurs élémentaires échantillonné sur une grille d'échantillonnage initiale, ayant un pas d'échantillonnage horizontal initial et un pas d'échantillonnage vertical initial, chaque pixel ayant un barycentre photométrique associé, le dispositif d'acquisition d'images générant, dans au moins une direction d'irrégularité parmi la direction horizontale et la direction verticale, un décalage des barycentres photométriques par rapport à la grille d'échantillonnage initiale. Ce dispositif comprend des moyens, appliquées à au moins une portion d'une dite matrice de pixels, de : - extraction de ladite portion d'un premier sous-ensemble de pixels correspondant au signal d'image échantillonné sur une première sous-grille d'échantillonnage et d'un deuxième sous-ensemble de pixels correspondant au signal d'image échantillonné sur une deuxième sous-grille d'échantillonnage, chacune desdites première et deuxième sous-grille d'échantillonnage ayant un pas d'échantillonnage selon la direction d'irrégularité égal au double du pas d'échantillonnage correspondant de la grille d'échantillonnage initiale ; - obtention d'un premier signal d'image à partir dudit premier sous-ensemble de pixels et d'un deuxième signal d'image à partir dudit deuxième sous-ensemble de pixels ; - application d'un premier filtre de ré-échantillonnage sur le premier signal d'image et d'un deuxième filtre de ré-échantillonnage sur le deuxième signal d'image, - ré-échantillonnage des signaux d'image filtrés pour obtenir un signal d'image ré- échantillonné sur une grille d'échantillonnage finale ayant un pas d'échantillonnage horizontal égal au pas d'échantillonnage horizontal initial et un pas d'échantillonnage vertical égal au pas d'échantillonnage vertical initial. Ce dispositif comporte en outre des moyens de détermination desdits premier et deuxième filtre de ré-échantillonnage en fonction dudit décalage des barycentres photométriques par rapport à la grille d'échantillonnage initiale et d'un positionnement souhaité de la grille de ré-échantillonnage finale par rapport à la grille d'échantillonnage initiale. Le dispositif de traitement de correction d'image numérique présente des avantages analogues aux avantages du procédé associé. Ce dispositif comporte des moyens de mise en oeuvre du procédé de traitement de correction d'image numérique tel que brièvement décrit ci-dessus. Selon un troisième aspect, l'invention concerne un système de traitement de correction d'images, ce système comportant un dispositif d'acquisition d'images en déplacement rectiligne relatif par rapport à une zone observée, comportant au moins une matrice de capteurs élémentaires, apte à acquérir une image numérique représentée par au moins une matrice de pixels correspondant à un signal d'image acquis par les capteurs élémentaires échantillonné sur une grille d'échantillonnage initiale, ayant un pas d'échantillonnage horizontal initial et un pas d'échantillonnage vertical initial, chaque pixel ayant un barycentre photométrique associé, le dispositif d'acquisition d'images générant, dans au moins une direction d'irrégularité parmi la direction horizontale et la direction verticale, un décalage des barycentres photométriques par rapport à la grille d'échantillonnage initiale, et un dispositif de traitement de correction tel que brièvement défini ci-dessus. Le système de traitement de correction d'image numérique présente des avantages analogues aux avantages du procédé de traitement de correction d'image numérique brièvement décrit ci-dessus. Dans un mode de réalisation, chaque capteur élémentaire est apte à convertir un rayonnement électromagnétique en charges électriques, et le dispositif d'acquisition d'images comporte des dispositifs d'anti-éblouissement aptes à évacuer des charges électriques excédentaires, les dispositifs d'anti-éblouissement étant positionnés latéralement selon une direction donnée tous les deux capteurs élémentaires. Selon une caractéristique particulière, chaque capteur élémentaire a une face photosensible orientée vers la zone observée, et chaque dispositif d'anti-éblouissement est positionné sur la face photosensible d'un capteur élémentaire. Selon une autre caractéristique particulière, chaque capteur élémentaire a une face photosensible orientée vers la zone observée, et chaque dispositif d'anti- éblouissement est positionné sur la face opposée à la face photosensible d'un capteur élémentaire. Dans un mode de réalisation, le dispositif d'acquisition d'images est un dispositif de type CMOS. D'autres caractéristiques et avantages de l'invention ressortiront de la description qui en est donnée ci-dessous, à titre indicatif et nullement limitatif, en référence aux figures annexées, parmi lesquelles : -la figure 1 illustre schématiquement le principe d'acquisition d'image par une matrice de capteurs élémentaires dans le cas de l'imagerie satellite ; -la figure 2 illustre une grille d'échantillonnage idéale correspondant à une matrice d'image, en l'absence d'utilisation de système d'anti-éblouissement latéral alterné ; - la figure 3 représente une matrice de pixels correspondant à une grille d'échantillonnage irrégulière en cas d'utilisation d'un système d'anti-éblouissement latéral alterné, - la figure 4 est un schéma illustrant les blocs fonctionnels d'une chaîne de traitement d'images selon un mode de réalisation ; -la figure 5 est un synoptique d'un procédé de traitement de correction d'image selon un mode de réalisation de l'invention; -les figures 6 et 7 représentent schématiquement des matrices de pixels correspondant à des sous-grilles d'échantillonnage ; - la figure 8 représente schématiquement un ensemble de blocs fonctionnels d'un dispositif apte à mettre en oeuvre l'invention. L'invention sera décrite ci-après plus particulièrement dans le domaine de l'acquisition d'image par satellite, par des matrices monodimensionnelles ou barrettes de capteurs élémentaires, dans le mode d'acquisition par le mode dit «push-broom », dans lequel chaque colonne de l'image numérique obtenue est associée physiquement à un capteur de la matrice monodimensionnelle de capteurs élémentaires, et dans lequel la grille d'échantillonnage est régulière, ayant un pas d'échantillonnage horizontal égal au pas d'échantillonnage vertical. Il est cependant à noter que le concept de l'invention s'applique dans d'autres cas d'application, notamment dans le cas de l'acquisition d'images par une matrice monodimensionnelle ou bidimensionnelle de capteurs élémentaires, selon d'autres modes de balayage d'un paysage ou si le pas d'échantillonnage vertical est différent du pas d'échantillonnage horizontal. Par exemple, l'invention s'applique également lorsque la matrice de capteurs élémentaires est fixe, et le paysage acquis défile au droit de la matrice de capteurs élémentaires. Le contexte de l'acquisition d'images numériques par satellite est illustré schématiquement dans la figure 1. Un satellite non représenté ayant une trajectoire de déplacement T comporte un dispositif d'acquisition d'images 2, qui est une matrice monodimensionnelle de capteurs élémentaires 4, qui est apte à acquérir des images via un système optique 6. Chaque capteur élémentaire 4 réalise l'acquisition d'une portion du paysage observé. Le paysage observé défile devant chaque capteur élémentaire 4 qui intègre le flux électromagnétique reçu pendant un temps de pose et le transforme en une charge électrique proportionnelle. Les charges obtenues sont converties en tensions électriques proportionnelles au flux électromagnétique observé. Par exemple, des capteurs CCD (pour « coupled charged device ») sont utilisés en tant que capteurs élémentaires. Une matrice d'image numérique, composée d'échantillons d'image ou pixels 8, chaque pixel ayant une valeur de radiométrie associée, appelée également valeur d'intensité associée, est obtenue à partir des tensions électriques précédemment obtenues. La matrice d'image numérique I est composée de lignes de pixels Li et de colonne de pixels Ci, chaque ligne d'image correspondant à l'acquisition du paysage à un instant donné. Dans l'exemple de la figure 1, chaque colonne de pixels est associée à un capteur élémentaire 4. La matrice I est formée sur une grille d'échantillonnage théorique, appelée grille d'échantillonnage initiale, illustrée à la figure 2. La figure 2 illustre une grille d'échantillonnage 10, formée des barycentres photométriques théoriques 12 (représentés par un symbole `+' sur la figure) de chaque surface d'acquisition 14 d'un rayonnement électromagnétique par un capteur élémentaire 4. Le pas d'échantillonnage horizontal de la grille d'échantillonnage 10, correspondant à la distance selon la direction horizontale entre deux barycentres 12 successifs, est noté P. Le pas d'échantillonnage P est de 1 pixel dans chaque direction dans cet exemple. Plus généralement, pour une grille d'échantillonnage on définit un pas d'échantillonnage vertical et un pas d'échantillonnage horizontal. Pour une grille d'échantillonnage régulière, le pas d'échantillonnage horizontal est égal au pas d'échantillonnage vertical. Dans le cas où un dispositif d'anti-éblouissement latéral, appelé dispositif AEL, est positionné tous les deux capteurs, sur la face photosensible des capteurs, ce dispositif occulte une partie du rayonnement électromagnétique reçu par les capteurs élémentaires 4. Dans d'autres modes d'implantation, par exemple lorsque le dispositif AEL est positionné tous les deux capteurs, sur la face opposée à la face photosensible des capteurs, un effet équivalent à une occultation directe du flux d'électrons recueilli est constaté. La figure 3 illustre une grille d'échantillonnage irrégulière 20 associée à une barrette de capteurs élémentaires 4 lorsqu'un dispositif AEL alterné est positionné tous les n=2 capteurs élémentaires 4. Dans l'exemple de la figure 3, une portion de matrice à six colonnes notées Co à 05 est obtenue par une barrette de six capteurs. Comme illustré sur cette figure, le positionnement des dispositifs AEL entraîne l'apparition de bandes d'occultation Z1 à Z3, chacune de ces bandes correspondant au positionnement d'un dispositif AEL entre deux capteurs voisins, tous les deux capteurs. La largeur ,e des zones d'occultation induites par les dispositifs AEL est constante et prédéterminée. Dans un exemple numérique, lorsque les surfaces d'acquisition 14 ont une dimension de 15 lm par 15 lm (micromètres), la largeur d'une bande AEL est de l'ordre de 5 lm correspondant à ,e =1/3 de pixel d'une grille d'échantillonnage régulière de pas égal à 1 pixel. Ainsi, chaque surface d'acquisition 22 correspondant à un capteur élémentaire donné est diminuée en largeur, et la grille d'échantillonnage 20 présente une irrégularité horizontale. Les barycentres 12 correspondant à une grille d'échantillonnage 10 régulière sont notés par le symbole `+' sur la figure 3. Les barycentres photométriques 24 réels de chaque surface d'acquisition 24 sont décalés par rapport aux barycentres 12, de manière différente selon la parité de l'indice de la colonne. Dans l'exemple de la figure 3, les barycentres photométriques des colonnes d'indice pair C2J, notés par le symbole 'o', sont décalés d'une distance ,e/4 vers la gauche du barycentre théorique 12, alors que les barycentres photométriques des colonnes d'indice pair C2,-0, notés par le symbole 'x', sont décalés d'une distance ,e/4 vers la droite du barycentre théorique 12. L'ensemble des barycentres des colonnes d'indice pair notés 'o' forment une première sous-grille d'échantillonnage GP, dite grille paire, qui est régulière dans chaque direction, de pas d'échantillonnage horizontal 2P et de pas d'échantillonnage vertical P. L'ensemble des barycentres des colonnes d'indice impair notés 'x' forment une deuxième sous-grille d'échantillonnage GI, dite grille impaire, qui est régulière dans chaque direction, de pas d'échantillonnage horizontal 2P et de pas d'échantillonnage vertical P. Cependant, la grille d'échantillonnage 20 formée par l'ensemble des barycentres photométriques 24 est irrégulière, la direction d'irrégularité étant la direction horizontale. L'invention propose de modifier le traitement radiométrique de l'image acquise par un dispositif d'acquisition induisant un décalage des barycentres photométriques des pixels selon au moins une des directions, afin d'améliorer la qualité de l'image numérique résultante. La figure 4 montre un synoptique d'une chaîne de traitement 30 pour l'acquisition d'images dans un mode d'acquisition tel que décrit ci-dessus en référence à la figure 1, faisant partie d'un système de traitement de correction d'images selon l'invention. La chaîne de traitement d'images 30 comporte une unité 32 d'amplification du signal électrique S obtenu en sortie des capteurs 4 connectée en entrée d'un codeur analogique-numérique 34. Une unité de correction radiométrique 36, adaptée pour implémenter un procédé de traitement de correction selon l'invention est connectée à la sortie de l'unité 34 et en entrée d'une unité 38 de traitement de restauration, tel une déconvolution ou un débruitage. Le traitement de correction radiométrique 36 est appliqué avant tout ré- échantillonnage de l'image numérique nécessaire pour une application de l'image numérique acquise. Un ordinogramme des principales étapes d'un procédé de traitement de correction d'images selon l'invention est représenté dans la figure 5. Dans ce mode de réalisation le cas d'un décalage des barycentres photométriques en colonnes adjacentes, tel qu'illustré dans la figure 3, est considéré. On appelle f le signal d'image continu idéal, qui est indépendant de toute représentation d'une image sur une grille d'échantillonnage pour une représentation numérique. Pour chaque ligne d'une portion d'image numérique I à traiter, on dispose des valeurs échantillonnées f(i,j),i=0,...,M ; j=0,1,2...., N, où M est l'indice maximal de ligne et N est l'indice maximal de colonne de la portion d'image à traiter. La grille d'échantillonnage idéale, appelé grille d'échantillonnage initiale, a un pas d'échantillonnage horizontal P. Lors d'une première étape 40, un premier sous-ensemble de pixels comportant les colonnes d'indice pair, est extrait de la grille d'échantillonnage complète de la portion d'image à traiter. Le premier sous-ensemble de pixels correspond à un premier sous-signal d'image, noté f1, échantillonné sur la première sous-grille d'échantillonnage paire GP, de pas d'échantillonnage horizontal P'=2P, Le premier sous-signal d'image f1 est complété par des zéros lors de l'étape suivante 42 pour obtenir un premier signal d'image F1 de pas d'échantillonnage horizontal égal au pas d'échantillonnage initial P. La formule suivante s'applique : F1 (i,2 j) = f (i,2 j) et F1 (i ,2 j +1) = 0 Comme illustré à la figure 6, une matrice de pixels (pair, ou premier ensemble de pixels, de même taille que la portion d'image numérique à traiter est obtenue. La première matrice de pixels correspond au premier signal d'image F1 échantillonné sur une grille d'échantillonnage de même pas d'échantillonnage que la grille initiale dans les deux directions horizontale et verticale. Lors d'une étape 44, un deuxième sous-signal d'image noté f2, correspondant aux colonnes d'indice impair, échantillonné sur la deuxième sous-grille d'échantillonnage impaire GI, de pas d'échantillonnage P'=2P, est extrait de la grille d'échantillonnage complète de la portion d'image à traiter. Le deuxième sous-signal d'image f2 est complété par des zéros lors de l'étape suivante 46 pour obtenir un deuxième signal d'image F2 de pas d'échantillonnage horizontal égal au pas d'échantillonnage initial P. La formule suivante s'applique : F2 (i,2 j +1) = f (i,2 j +1) et F2 (i,2 j) = 0 . Comme illustré à la figure 7, une matrice de pixels 1,1',r, ou deuxième ensemble de pixels, de même taille que la portion d'image numérique à traiter est obtenue. La deuxième matrice de pixels correspond au deuxième signal d'image F2 échantillonné sur une grille d'échantillonnage de même pas d'échantillonnage que la grille initiale dans les deux directions horizontale et verticale. Après l'étape 46, une étape 50 consiste à obtenir deux filtres de ré- échantillonnage G1, G2. Ces filtres sont définis chacun par un ensemble de coefficients, une taille de support de filtre, des valeurs de pondération d'apodisation. Les filtres de ré-échantillonnage G1, G2 sont normalisés : la somme des coefficients de filtres sur leur support de taille finie est égale à 1. Dans le mode de réalisation préféré, la formule mathématique suivante est utilisée pour fournir les coefficients d'un filtre de ré-échantillonnage G : G(x) = sinc(zy) - - cotan 1 +ct x sinc- ( où sinc est la function sinus 2 71--- 2) 2 ) cardinal, cotan est la fonction co-tangente, et a = / 2 est la demi-largeur des zones d'occultation, comme illustré à la figure 3. Avantageusement, la fonction G(x) est la fonction mathématique qui permet de reconstruire, à partir d'une grille d'échantillonnage décalée (correspondant aux symboles 'o' et 'x'), l'image qui aurait été acquise par un dispositif d'acquisition régulier exempt de décalages des barycentre photométrique (correspondant directement aux symboles `+') avec le même niveau de netteté, sous certaines hypothèses concernant le dispositif d'acquisition : image et filtres de supports non limités, FTM identique entre les deux sous-grilles d'échantillonnage, régularité parfaite de l'échantillonnage des sous-grilles, absence de repliement de spectre au pas final P, absence de bruit instrumental d'acquisition. Les propriétés avantageuses de la fonction G(x) sont en grande partie préservées lorsque ces hypothèses ne sont pas parfaitement respectées, dans le cas d'un dispositif d'acquisition réel. Dans ce mode de réalisation, le filtrage est effectué par ligne de pixels de chaque signal image, dans la mesure où les barycentres photométriques sont régulièrement espacés dans la direction verticale, donc aucune correction radiométrique n'est nécessaire au niveau des colonnes pour améliorer la qualité d'image. Un premier filtre de ré-échantillonnage G1(x) et un deuxième filtre de ré- échantillonnage G2(x) sont définis respectivement pour filtrer le premier signal d'image F1 et le deuxième signal d'image F2. Dans le mode de réalisation préféré, les premier et deuxième filtre de ré-échantillonnage sont définis comme suit : ( ct Gi(x)= G x + - (Eq 1) - 2) Et ( G2(x) =G 1+ -a- : (Eq 2) - 2 Cette définition par les formules (Eq 1) et (Eq 2) des filtres de ré-échantillonnage permet de réaliser le recalage de chacun des signaux d'image F1 et F2 sur une grille d'échantillonnage finale qui est une grille médiane, notée GM sur la figure 3, dont les barycentres photométriques par pixel sont situés à égale distance des barycentres des sous-grilles paire et impaire respectives. Avantageusement, la contribution des premier et deuxième signaux d'image est équilibrée lorsque les filtres de ré-échantillonnage sont définis par les formules (Eq 1) et (Eq 2). Les coefficients de chacun des filtres G1 et G2 sont obtenus sur un support fini, de préférence compris entre 16 et 256, et typiquement égal à 128. Une apodisation ou fenêtrage est appliquée, consistant à multiplier les filtres G1 et G2, sur leur support, par une fonction symétrique (centrée au milieu du support) valant 1 au centre et à décroissance douce vers les bords. Classiquement, les fonctions dites de Hann, de Hamming, gaussienne ou en cosinus relevé sont utilisées. Lors d'une première étape de ré-échantillonnage 50 le premier filtre G1 est appliqué au premier signal F1, ligne par ligne, pour obtenir un premier signal d'image filtré, noté F1 * G1 , où le symbole * représente le produit de convolution. Ensuite, lors d'une deuxième étape de ré-échantillonnage 52, le deuxième filtre G2 est appliqué au signal F2, ligne par ligne, pour obtenir un deuxième signal d'image filtré, noté F2 * G2 . Finalement, lors d'une étape de sommation 54, les deux signaux filtrés sont additionnés, pixel par pixel, pour obtenir un signal final ré-échantillonné F, qui correspond au signal d'image initial f reconstitué sur une grille régulière finale GM de pas d'échantillonnage P. F(i, j)= Fi(i, j)* Gi(j)+ F2 (i, j) * G2 ( j) Il est à noter que selon une variante, les étapes 50, 52 et 54 sont effectuées simultanément. Selon une variante du mode de réalisation décrit ci-dessus, lors de l'étape 50, les premier et deuxième filtres de ré-échantillonnage sont définis comme suit : G, (x)= G(x) (Eq 3) Et G2 (X) = G(a +1- x) (Eq 4) Avec la définition donnée par les formules (Eq 3) et (Eq 4) ci-dessus, la grille de ré-échantillonnage choisie est une grille de pas d'échantillonnage horizontal P, positionnée sur les barycentres photométriques 24 notés 'o' sur la figure 3. Dans ce mode de réalisation, les barycentres photométriques des colonnes d'indice pair ne sont pas modifiés, seuls les barycentres photométriques des colonnes d'indice impair sont modifiés. Selon une autre variante, lors de l'étape 50, les premier et deuxième filtres de ré-échantillonnage sont définis comme suit : Gi(x) = G(a +1- x) (Eq 5) Et G2 (X) = G(X) (Eq 6) Avec la définition ci-dessus, la grille de ré-échantillonnage choisie la grille de ré- échantillonnage choisie est une grille de pas d'échantillonnage horizontal P, positionnée sur les barycentres photométriques 24 notés 'x' sur la figure 3. Dans ce mode de réalisation, les barycentres photométriques des colonnes d'indice impair ne sont pas modifiés, seuls les barycentres photométriques des colonnes d'indice pair sont modifiés. D'une manière plus générale, on définit les filtres de ré-échantillonnage G1 et G2 à partir de la formule de filtre de ré-échantillonnage G et de la demi-largeur a des zones d'occultation, pour un premier décalage dl par rapport aux barycentres photométriques de la première sous-grille d'échantillonnage et un décalage d2 par rapport aux barycentres photométriques de la deuxième sous-grille d'échantillonnage. La condition suivante s'applique d1 + d2=a pour compenser la demi-largeur a des zones d'occultation. Les décalages dl et d2 définissent le positionnement de la grille d'échantillonnage finale : (x) = G(x + (Eq 7) Et G2 (X) = G(C12 ± 1 x) (Eq 8) Les procédés de traitement de correction de l'invention sont mis en oeuvre par un dispositif programmable de type ordinateur, station de travail, dont les principaux blocs fonctionnels sont illustrés à la figure 8. Un dispositif programmable 60 apte à mettre en oeuvre les procédés de l'invention comprend une unité centrale de traitement 68, ou CPU, apte à exécuter des instructions de programme informatique lorsque le dispositif 60 est mis sous tension. Dans un mode de réalisation, une unité centrale de traitement multi-processeurs est utilisée, permettant d'effectuer des calculs parallèles. Le dispositif 60 comporte également des moyens de stockage d'informations 70, par exemple des registres, aptes à stocker des instructions de code exécutable permettant la mise en oeuvre de programmes comportant des instructions de code aptes à mettre en oeuvre le procédé de traitement de correction d'images selon l'invention. Le dispositif 60 comporte des moyens de commande 64 permettant de mettre à jour des paramètres et de recevoir des commandes d'un opérateur. Lorsque le dispositif programmable 60 est un dispositif embarqué, les moyens de commande 64 comprennent un dispositif de télécommunication permettant de recevoir des commandes et des valeurs de paramètres à distance. Alternativement et de manière optionnelle, les moyens de commande 64 sont des moyens de saisie de commandes d'un opérateur, par exemple un clavier. De manière optionnelle, le dispositif programmable 60 comprend un écran 62 et un moyen supplémentaire de pointage 66, tel une souris. Les divers blocs fonctionnels du dispositif 60 décrits ci-dessus sont connectés via un bus de communication 72. Un système de traitement de correction d'images selon l'invention comporte un dispositif d'acquisition d'images, tel le dispositif d'acquisition d'images 2 de la figure 1, générant un décalage des barycentres photométriques selon au moins une direction, et un dispositif de traitement de correction mis en oeuvre par un dispositif programmable tel que décrit ci-dessus en référence à la figure 8. Selon un premier mode de réalisation, le dispositif d'acquisition d'images est une matrice de capteurs élémentaires portée par un satellite et munie d'un dispositif AEL positionné tous les deux capteurs sur la face photosensible des capteurs. Selon un deuxième mode de réalisation, le dispositif d'acquisition d'images est une matrice de capteurs élémentaires portée par un satellite et munie d'un dispositif AEL positionné tous les deux capteurs, sur la face opposée à la face photosensible des capteurs, ayant un effet équivalent à une occultation directe du flux d'électrons recueilli est constaté. Selon un troisième mode de réalisation, le dispositif d'acquisition d'images est une matrice de capteurs CMOS comprenant une fonction électronique, par exemple un bus électrique, implantée tous les deux capteurs selon une direction donnée, résultant en une occultation des surfaces photosensibles et en un décalage des barycentres photométriques selon la direction donnée. L'invention s'applique, de manière plus générale, dans tout système comportant un dispositif d'acquisition d'images réalisant une occultation de la surface d'incidence du rayonnement électromagnétique ou une perturbation du flux d'électrons de ce rayonnement électromagnétique. Le principe de l'invention a été décrit ci-dessus dans un exemple dans lequel une irrégularité des barycentres photométriques dans la direction horizontale est présente. De manière analogue, l'invention s'applique lorsque la direction d'irrégularité des barycentres photométriques est la direction verticale. En outre, le procédé de l'invention est également applicable de manière séparable dans la direction horizontale d'une part et dans la direction verticale d'autre part, pour compenser des irrégularités verticales et horizontales, qui surviennent notamment lors de l'utilisation d'un dispositif d'acquisition comportant une matrice bidimensionnelle de capteurs élémentaires.30 REVENDICATIONS1.- Procédé de traitement de correction d'image numérique acquise par un dispositif d'acquisition d'images comportant au moins une matrice de capteurs élémentaires, en déplacement relatif par rapport à une zone observée, une image acquise (I) étant représentée par au moins une matrice de pixels correspondant à un signal d'image acquis par les capteurs élémentaires échantillonné sur une grille d'échantillonnage initiale (10), ayant un pas d'échantillonnage horizontal initial et un pas d'échantillonnage vertical initial, chaque pixel ayant un barycentre photométrique (12) associé, le dispositif d'acquisition d'images générant, dans au moins une direction d'irrégularité parmi la direction horizontale et la direction verticale, un décalage des barycentres photométriques (24) par rapport à la grille d'échantillonnage initiale (10), caractérisé en ce qu'il comprend les étapes, appliquées à au moins une portion d'une dite matrice de pixels, de : - extraction (40, 44) de ladite portion d'un premier sous-ensemble de pixels correspondant au signal d'image échantillonné sur une première sous-grille d'échantillonnage (GP) et d'un deuxième sous-ensemble de pixels correspondant au signal d'image échantillonné sur une deuxième sous-grille d'échantillonnage (GI), chacune desdites première et deuxième sous-grille d'échantillonnage (GP, GI) ayant un pas d'échantillonnage selon la direction d'irrégularité égal au double du pas d'échantillonnage correspondant de la grille d'échantillonnage initiale, - obtention (42, 46) d'un premier signal d'image à partir dudit premier sous- ensemble de pixels et d'un deuxième signal d'image à partir dudit deuxième sous-ensemble de pixels, - application (50, 52, 54) d'un premier filtre de ré-échantillonnage (G1) sur le premier signal d'image et d'un deuxième filtre de ré-échantillonnage (G2) sur le deuxième signal d'image, - ré-échantillonnage (50, 52, 54) des signaux d'image filtrés pour obtenir un signal d'image ré-échantillonné sur une grille d'échantillonnage finale (GM) ayant un pas d'échantillonnage horizontal égal au pas d'échantillonnage horizontal initial et un pas d'échantillonnage vertical égal au pas d'échantillonnage vertical initial, lesdits premier (G1) et deuxième (G2) filtre de ré-échantillonnage étant déterminés (48) en fonction dudit décalage des barycentres photométriques par rapport à la grille d'échantillonnage initiale et d'un positionnement souhaité de la grille de ré- échantillonnage finale par rapport à la grille d'échantillonnage initiale. 2.- Procédé selon la revendication 1, caractérisé en ce que ladite image acquise est acquise par une matrice monodimensionnelle de capteurs élémentaires, et en ce que ledit premier sous-ensemble de pixels comprend tous les pixels de la portion de matrice de pixels appartenant aux colonnes ou lignes d'indice pair et ledit deuxième sous- ensemble de pixels comprend tous les pixels de la portion d'image traitée appartenant aux colonnes ou lignes d'indice impair de ladite portion de matrice traitée. 3.- Procédé selon l'une quelconque des revendications 1 ou 2, caractérisé en que le dispositif d'acquisition d'images induit des zones d'occultation (Z1, Z2, Z3) de l'acquisition de charges par les capteurs élémentaires auxquels ils sont associés, disposées tous les deux pixels selon la direction d'irrégularité, et en ce la détermination (48) desdits premier et deuxième filtre est fonction de la largeur d'une dite zone d'occultation. 4.- Procédé selon la revendication 3, caractérisé en ce lesdits premier et deuxième filtre sont définis à partir d'un filtre de ré-échantillonnage G par les formules suivantes : G1(x) = G(x + di ) Et G2 (X) = G(d2 + 1 - 4 . où dl est une distance, selon la direction d'irrégularité, entre les barycentres photométriques de la première sous-grille d'échantillonnage et la grille d'échantillonnage finale, et d2 est une distance, selon la direction d'irrégularité, entre les barycentres photométriques de la deuxième sous-grille d'échantillonnage et la grille d'échantillonnage finale. 5.- Procédé selon la revendication 4, caractérisé en ce que la détermination des premier et deuxième filtre est effectuée à partir d'un filtre de ré-échantillonnage défini par : G(x) = sinc(2x) - 7 I - - cotan( 7.1. - 1+ a x , ( z .x sinc - - - où sinc est la fonction sinus 2 ) 2 2) cardinal, cotan est la fonction co-tangente, et a est égal à la demi-largeur d'une dite zone d'occultation, x et a étant exprimés en pixels de la grille d'échantillonnage initiale. 6.- Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce que la détermination (48) desdits premier et deuxième filtre comprend une déterminationde longueur de support de chacun desdits filtres et l'application d'une fenêtre d'apodisation. 7.- Procédé selon la revendication 6, caractérisé en ce que la détermination desdits premier et deuxième filtre comprend une normalisation des coefficients de chacun desdits premier et deuxième filtre. 8.- Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le ré-échantillonnage sur la grille d'échantillonnage finale (GM) comporte une addition (54) pixel par pixel des signaux d'image filtrés pour obtenir une valeur de pixel du signal d'image ré-échantillonné. 9.- Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que l'obtention d'un premier signal d'image consiste à insérer une colonne ou ligne de zéros entre deux colonnes ou lignes consécutives dudit premier sous-ensemble de pixels, et en ce l'obtention d'un deuxième signal d'image consiste à insérer une colonne ou ligne de zéros entre deux colonnes ou lignes consécutives dudit deuxième sous-ensemble de pixels. 10.- Dispositif de traitement de correction d'image numérique acquise par un dispositif d'acquisition d'images en déplacement rectiligne relatif par rapport à une zone observée, comportant au moins une matrice de capteurs élémentaires, une image acquise étant représentée par au moins une matrice de pixels correspondant à un signal d'image acquis par les capteurs élémentaires échantillonné sur une grille d'échantillonnage initiale, ayant un pas d'échantillonnage horizontal initial et un pas d'échantillonnage vertical initial, chaque pixel ayant un barycentre photométrique associé, le dispositif d'acquisition d'images générant, dans au moins une direction d'irrégularité parmi la direction horizontale et la direction verticale, un décalage des barycentres photométriques par rapport à la grille d'échantillonnage initiale, caractérisé en ce qu'il comprend des moyens, appliquées à au moins une portion d'une dite matrice de pixels, de : - extraction de ladite portion d'un premier sous-ensemble de pixels correspondant au signal d'image échantillonné sur une première sous-grille d'échantillonnage et d'un deuxième sous-ensemble de pixels correspondant au signal d'image échantillonné sur une deuxième sous-grille d'échantillonnage, chacune desdites première et deuxième sous-grille d'échantillonnage ayant un pas d'échantillonnage selon la directiond'irrégularité égal au double du pas d'échantillonnage correspondant de la grille d'échantillonnage initiale, - obtention d'un premier signal d'image à partir dudit premier sous-ensemble de pixels et d'un deuxième signal d'image à partir dudit deuxième sous-ensemble de pixels, - application d'un premier filtre de ré-échantillonnage sur le premier signal d'image et d'un deuxième filtre de ré-échantillonnage sur le deuxième signal d'image, - ré-échantillonnage des signaux d'image filtrés pour obtenir un signal d'image ré-échantillonné sur une grille d'échantillonnage finale ayant un pas d'échantillonnage horizontal égal au pas d'échantillonnage horizontal initial et un pas d'échantillonnage vertical égal au pas d'échantillonnage vertical initial, le dispositif comportant en outre des moyens de détermination desdits premier et deuxième filtre de ré-échantillonnage en fonction dudit décalage des barycentres photométriques par rapport à la grille d'échantillonnage initiale et d'un positionnement souhaité de la grille de ré-échantillonnage finale par rapport à la grille d'échantillonnage initiale. 11.- Système de traitement de correction d'images caractérisé en ce qu'il comporte un dispositif d'acquisition d'images en déplacement rectiligne relatif par rapport à une zone observée, comportant au moins une matrice de capteurs élémentaires, apte à acquérir une image numérique représentée par au moins une matrice de pixels correspondant à un signal d'image acquis par les capteurs élémentaires échantillonné sur une grille d'échantillonnage initiale, ayant un pas d'échantillonnage horizontal initial et un pas d'échantillonnage vertical initial, chaque pixel ayant un barycentre photométrique associé, le dispositif d'acquisition d'images générant, dans au moins une direction d'irrégularité parmi la direction horizontale et la direction verticale, un décalage des barycentres photométriques par rapport à la grille d'échantillonnage initiale, et un dispositif de traitement de correction conforme à la revendication 10. 12.- Système de traitement de correction d'images selon la revendication 11, caractérisé en ce que chaque capteur élémentaire est apte à convertir un rayonnement électromagnétique en charges électriques, et en ce que le dispositif d'acquisition d'images comporte des dispositifs d'anti-éblouissement aptes à évacuer des charges électriques excédentaires, les dispositifs d'anti-éblouissement étant positionnés latéralement selon une direction donnée tous les deux capteurs élémentaires.35 13.- Système de traitement de correction d'images selon la revendication 12, caractérisé en ce que chaque capteur élémentaire a une face photosensible orientée vers la zone observée, et en ce que chaque dispositif d'anti-éblouissement est positionné sur la face photosensible d'un capteur élémentaire. 14.- Système de traitement de correction d'images selon la revendication 12, caractérisé en ce que chaque capteur élémentaire a une face photosensible orientée vers la zone observée, et en ce que chaque dispositif d'anti-éblouissement est positionné sur la face opposée à la face photosensible d'un capteur élémentaire. 15.- Système de traitement de correction d'images selon la revendication 11, caractérisé en ce que le dispositif d'acquisition d'images est un dispositif de type CMOS. 10