Brevet FR3055727A1 — publié le 20180309

Procédé et dispositif de caractérisation des aberrations d'un système optique La présente invention concerne un procédé de caractérisation d'au moins une aberration d'un système optique d'un dispositif d'acquisition d'images, apte à acquérir au moins deux images selon un champ de prise de vue, la deuxième image présentant une aberration différentielle par rapport à la première. L'invention concerne également un dispositif et un programme d'ordinateur associés. L'invention se situe dans le domaine de l'optique active, consistant à mesurer et à compenser des défauts optiques, également appelés aberrations, dans des dispositifs d'acquisition d'images. L'invention s'applique notamment dans le domaine de l'imagerie spatiale. Dans ce domaine, les dispositifs d'acquisition d'images comportent des télescopes munis de miroirs. Afin d'obtenir une bonne résolution de prise d'image il est fondamental de réaliser un télescope comportant le moins possible d'aberrations afin d'atteindre la limite physique imposée par la diffraction. Ceci est difficile à réaliser pour de multiples raisons : les bougés relatifs des miroirs constituant la combinaison optique, la qualité de réalisation des surfaces optiques des miroirs, la stabilité mécanique dans le temps vis-à-vis des contraintes thermoélastiques et de gravité induisent des aberrations optiques. La mesure et la compensation des aberrations optiques permettent ainsi de réduire fortement la complexité de réalisation de tels systèmes optiques. La caractérisation et la compensation des aberrations optiques, notamment de défauts de focalisation, est en particulier utile dans le domaine de l'acquisition d'images par balayage de fauchée (en anglais « pushbroom »), réalisé par la fixation de barrettes de détecteurs d'images sur un satellite et un défilement du satellite au-dessus de la scène captée dans le champ de prise de vue. Enfin, dans le domaine de l'acquisition d'image grand public, les appareils photographiques sont munis d'objectifs à lentilles, dont la longueur focale est ajustée par déplacement des lentilles. Dans ce domaine également, il est utile de disposer d'un système automatique de caractérisation des aberrations optiques, notamment d'un défaut de focalisation. Il existe divers systèmes permettant de réaliser un réglage de focalisation, appelés systèmes d'autofocus. Il est également utile de proposer, pour certains dispositifs d'acquisition d'image grand public, une méthode alternative d'estimation automatique d'aberrations comme par exemple de la défocalisation. A cet effet, l'invention propose un procédé de caractérisation d'au moins une aberration optique d'un système optique d'un dispositif d'acquisition d'images, le système optique ayant une fonction de transfert optique associée dépendant des aberrations, le dispositif d'acquisition étant apte à acquérir au moins deux images dans un champ de 3055727 2 prise vue de manière à introduire une aberration différentielle entre les deux images, chaque image étant définie par un signal numérique d'image. Ce procédé comporte des étapes de : -obtention d'une première image et d'une deuxième image d'une même zone du 5 champ de prise de vue, la deuxième image étant acquise avec une aberration différentielle par rapport à la première image, chaque acquisition d'image ayant une fonction de transfert optique associée, -utilisation d'un modèle convolutif de formation d'images pour les images acquises utilisant une fonction de transfert optique dépendant d'au moins une aberration optique, 10 -calcul d'au moins un paramètre caractérisant au moins une aberration optique par minimisation d'une fonctionnelle, définie en utilisant le modèle convolutif de formation d'images, exprimée en fonction des première et deuxième images acquises, d'une modélisation des fonctions de transfert optiques associées aux deux acquisitions dépendant de la ou des aberrations optique à caractériser ainsi que de l'aberration 15 différentielle. Avantageusement, le procédé de l'invention permet d'estimer les aberrations optiques d'un dispositif d'acquisition d'images à partir de deux ou plusieurs images présentant une différence d'aberration connue ou estimée. Avantageusement, le procédé de l'invention s'applique pour des images pour lesquelles l'impact de l'aberration 20 supplémentaire est relativement faible et imperceptible pour un utilisateur. Une telle acquisition d'images peut être réalisée, dans certains cas, sans modification du dispositif d'acquisition d'images. Le procédé selon l'invention peut présenter une ou plusieurs des caractéristiques ci-dessous. 25 Le système optique à caractériser comprend au moins un ensemble de barrettes de détecteurs disposées en quinconce, et dans lequel lesdites première et deuxième images appartiennent à des zones de recouvrement entre barrettes de détecteurs adjacentes. Le système optique comprend une pluralité d'ensembles de barrettes de 30 détecteurs, chaque ensemble de barrettes de détecteurs ayant une bande spectrale d'acquisition associée, le procédé étant appliqué sur l'ensemble des zones de recouvrement des différentes bandes spectrales correspondant à une même zone du champ de vue, afin de déterminer des aberrations relatives à chacune de ces zones du champ de vue. 35 Le procédé est appliqué sur au moins une zone de recouvrement de chaque bande spectrale, de manière à caractériser des aberrations optiques par bande spectrale. 3055727 3 La fonctionnelle dépend en outre d'une longueur d'onde centroïde associée à la bande spectrale d'acquisition, et la longueur d'onde centroïde est estimée lors de la minimisation de ladite fonctionnelle. La fonctionnelle dépend en outre d'un décalage spatial entre la première image et 5 la deuxième images acquises, ledit décalage spatial étant estimé lors de la minimisation de ladite fonctionnelle. Dans lequel l'étape de calcul d'au moins un paramètre caractérisant au moins une aberration optique met en oeuvre une modélisation de ladite fonction de transfert optique en fonction d'au moins un coefficient d'une décomposition dans une base de polynômes 10 de Zernike d'une erreur de front d'onde associée au système optique du dispositif d'acquisition. La minimisation de la fonctionnelle est effectuée par une méthode d'estimation itérative. Le procédé comprend une modélisation de chacune des fonctions de transfert 15 optique associées aux acquisitions d'image par un polynôme de degré inférieur ou égal à deux, défini par des variables représentatives des aberrations optiques à caractériser et une minimisation de ladite fonctionnelle par une méthode non itérative. La valeur d'aberration différentielle entre la première et la deuxième image acquise est estimée lors de la minimisation de ladite fonctionnelle. 20 L'aberration différentielle entre la première image et la deuxième image est une défocalisation relative entre la première image et la deuxième image. Selon un autre aspect, l'invention concerne un dispositif de caractérisation d'au moins une aberration optique d'un système optique d'un dispositif d'acquisition d'images, le système optique ayant une fonction de transfert optique associée dépendant des 25 aberrations, le dispositif d'acquisition étant apte à acquérir au moins deux images dans un champ de prise vue de manière à introduire une aberration différentielle entre les deux images, chaque image étant définie par un signal numérique d'image. Ce dispositif comporte : -module d'obtention d'une première image et d'une deuxième image d'une même 30 zone du champ de prise de vue, acquises par le dispositif d'acquisition d'images, la deuxième image étant acquise avec une aberration différentielle par rapport à la première image, chaque acquisition d'image ayant une fonction de transfert optique associée, -un processeur adapté à mettre en oeuvre un calcul d'au moins un paramètre caractérisant au moins une aberration optique par minimisation d'une fonctionnelle 35 exprimée en fonction des première et deuxième images acquises, utilisant un modèle 3055727 4 convolutif de formation d'images pour les images acquises utilisant une fonction de transfert optique dépendant d'au moins une aberration optique, ladite fonctionnelle étant également exprimée en fonction d'une modélisation des fonctions de transfert optiques associées aux deux acquisitions dépendant de la ou des 5 aberrations optique à caractériser ainsi que de l'aberration différentielle. Selon un autre aspect, l'invention concerne un programme d'ordinateur comportant des instructions logicielles qui, lorsqu'elles sont mises en oeuvre par un dispositif programmable, mettent en oeuvre un procédé de caractérisation d'au moins une aberration optique d'un système optique d'un dispositif d'acquisition d'images tel que 10 brièvement décrit ci-dessus. 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 : 15 -la figure 1 illustre schématiquement un dispositif d'acquisition d'images et un dispositif de caractérisation des aberrations optiques selon un mode de réalisation ; -la figure 2 illustre schématiquement un dispositif d'acquisition d'images par satellite ; -la figure 3 est une vue en coupe schématisée d'un dispositif d'acquisition 20 d'images à barrettes de détecteurs ; -la figure 4 est un logigramme des principales étapes d'un procédé de de caractérisation d'une aberration optique selon un premier mode de réalisation de l'invention ; -la figure 5 est un logigramme des principales étapes d'un procédé de de 25 caractérisation d'une aberration optique selon un deuxième mode de réalisation de l'invention ; -la figure 6 est un synoptique des principaux modules d'un dispositif programmable apte à mettre en oeuvre le procédé de l'invention. 30 Le procédé de l'invention s'applique dans le domaine de l'imagerie satellite mais n'est pas limité à ce domaine. Il sera décrit ci-après dans son application à l'imagerie satellite pour l'observation de scène complexe, par exemple dans le contexte de l'observation de la terre par satellite, étant entendu que le procédé est applicable pour l'acquisition de tout autre type d'image. 35 La figure 1 illustre schématiquement un dispositif d'acquisition d'images 10 utilisé dans un mode de réalisation de l'invention. 3055727 5 Le dispositif d'acquisition d'images 10 comporte un système optique 12, comportant une optique de formation d'image 14, par exemple une ou plusieurs lentilles et un ou plusieurs détecteurs 16 aptes à capter des ondes électromagnétiques formant un signal lumineux représentatif d'une scène, par exemple un paysage, et à le transformer 5 en signal électrique. Le système optique est apte à capter un ou plusieurs signaux lumineux, chacun appartenant à une bande spectrale définie par un intervalle de longueurs d'onde donné. Par exemple, le système optique 12 comporte une pluralité de matrices monodimensionnelles de détecteurs, également appelées barrettes de détecteurs. 10 De préférence, le système optique est apte à capter plusieurs signaux appartenant à des bandes spectrales différentes. Il s'agit dans ce cas d'un système optique multispectral. Le système optique 12 est associé à un module de numérisation 18, permettant de transformer le signal électrique en un signal numérique d'image, échantillonné sur une 15 grille d'échantillonnage bidimensionnelle. De manière connue, un signal numérique d'image est représenté par une matrice de pixels, chaque pixel ayant une position spatiale dans la grille d'échantillonnage et une valeur de radiométrie (ou intensité) associée. Une telle représentation de l'image numérique est appelée représentation dans le domaine spatial. 20 Le dispositif d'acquisition 10 comporte, de manière optionnelle, un module 19 de décalage de plan focal permettant d'introduire une défocalisation relative entre images acquises successivement par le système optique 12. En variante, le dispositif d'acquisition 10 comporte un module 19 permettant d'introduire une aberration différentielle autre qu'une défocalisation relative. 25 Dans une variante, le dispositif d'acquisition d'images 10 ne comporte pas de module 19 spécifique, mais le procédé de l'invention exploite le fait que chaque bande spectrale est acquise par plusieurs barrettes physiques disposées en quinconce et formant un ensemble de barrettes, chaque barrette ayant un numéro associé, telles que les images acquises par les barrettes de numéro pair, dites barrettes paires, présentent 30 une aberration supplémentaire par rapport aux images acquises par les barrettes de numéro impair, dites barrettes impaires, comme dans l'exemple illustré ci-après en référence à la figure 3. Dans cette variante, on exploite les zones de recouvrement ZR entre barrettes paires et impaires. Cela permet d'avoir plusieurs mesures en un même endroit du champ de prise de vue par exploitation des zones de recouvrement ZR 35 homologues de chacune des bandes spectrales et ainsi d'améliorer la performance de mesure par moyennage et filtrage en l'absence de chromatisme, ou de mesurer ce 3055727 6 chromatisme s'il existe. La mesure dans l'ensemble des zones de recouvrement ZR donne accès aux aberrations optiques dans tout le champ de vue du système optique. Deux ensembles de zones de recouvrement, notées respectivement ZR, et ZR2, sont représentés à la figure 3. 5 Le dispositif d'acquisition d'images 10 est connecté, par tout moyen de connexion approprié, par exemple via un bus de transmission de données, à un dispositif de traitement 20, muni d'un ou plusieurs processeurs, et apte à mettre en oeuvre des programmes comportant des instructions logicielles aptes à mettre en oeuvre le procédé de caractérisation d'aberrations optiques selon l'invention. 10 Le dispositif de traitement 20 comporte un module 22 de réception d'images numériques du dispositif d'acquisition d'images 10 et un module 24 de traitement des images numériques reçues. Le dispositif de traitement 20 comporte également un module 26 d'estimation des aberrations optiques recherchées par minimisation d'une fonctionnelle dépendant de deux images acquise avec une aberration différentielle. 15 Enfin, le dispositif de traitement 20 comporte un module 28 d'utilisation du ou des paramètres caractérisant les aberrations optiques calculés par le module de calcul 26 pour une correction calculatoire du signal d'image acquis pour compenser au mieux l'impact des aberrations sur l'image ou pour l'obtention d'une instruction de réglage de la mise au point du dispositif d'acquisition d'image 10, qui est transmise par exemple au 20 module 19 de ce dispositif. La figure 2 illustre schématiquement un dispositif d'acquisition d'images dans le contexte de l'acquisition d'images numériques par satellite. Un satellite non représenté, ayant une trajectoire de déplacement T, est muni d'un dispositif d'acquisition d'images comportant un système optique 30. 25 Le système optique 30 comporte une pluralité de barrettes de détecteurs 32a, 32b, 32c, 32d. L'acquisition d'un paysage observé se fait via un télescope 34. Chaque détecteur d'une barrette de détecteurs réalise l'acquisition d'une portion de la scène (également appelée paysage) observée, enregistrée dans une matrice d'image numérique I sous 30 forme d'un échantillon d'image 36 ou pixel dans les zones d'acquisition simples, pour lesquelles une seule barrette de détecteurs reçoit un flux lumineux. Chaque pixel d'image a une valeur de radiométrie associée. Le dispositif d'acquisition d'image 30 présente, dans les zones de recouvrement inter-barrettes, une aberration différentielle, par exemple une défocalisation relative. 35 Avantageusement, cette aberration différentielle est volontairement introduite et connue. 3055727 7 La figure 3 illustre schématiquement une vue en coupe des barrettes de détecteurs dans un dispositif d'acquisition d'image tel que décrit ci-dessus. Dans le cas de l'imagerie satellite, un dispositif d'acquisition d'images comporte plusieurs ensembles de barrettes de détecteurs, notés respectivement 401, 402, 40,, 40N-1, 5 40N, correspondant à une acquisition d'images appartenant à N bandes spectrales. Chaque ensemble de barrettes 40, comporte une ou plusieurs zones de recouvrement inter-barrettes, comme illustré à la figure 3 où on observe les zones ZR, et ZR2. En pratique, dans un mode de réalisation, les zones de recouvrement ont une 10 largeur supérieure à un nombre donné de pixels, par exemple 50 pixels. Par exemple, dans un mode de réalisation, les ensembles de barrettes de détecteurs 40, correspondent chacun à une voie multi-spectrale, associée à une bande spectrale étroite donnée, en mode multi-spectral, et l'ensemble 40N correspond au mode de prise de vue panchromatique, appelé voie panchromatique, correspondant à une 15 bande spectrale large et ayant une résolution spatiale plus élevée que les autres voies multi-spectrales. De manière classique, un rapport de résolution de 4 existe entre la voie panchromatique et chaque voie multi-spectrale. Dans un premier mode de réalisation, il existe un différentiel, également appelé 20 décalage, natif de focalisation entre les plans focaux de l'ensemble 40N et de l'un quelconque des ensembles de barrettes 40,, et ce décalage est estimé et utilisé dans le procédé de caractérisation d'aberrations optiques comme expliqué en détail ci-après. De même, dans ce mode de réalisation, un décalage de focalisation entre les plans focaux de deux ensembles de barrettes 40, et 40,, dans les zones de recouvrement 25 inter-barrettes, est estimé et utilisé dans le procédé de caractérisation d'aberrations optiques comme expliqué en détail ci-après. Dans un deuxième mode de réalisation, un décalage de focalisation prédéterminé et maîtrisé est introduit, par construction, entre des ensembles de barrettes choisis, de manière à construire un système optique apte à réaliser un décalage de focalisation 30 prédéterminé. Dans ce mode de réalisation, ce décalage de focalisation prédéterminé est utilisé dans le procédé de caractérisation d'aberration optique. Avantageusement, l'introduction d'une aberration différentielle de défocalisation pouvant aller jusqu'à 200nm entre bandes multi-spectrales a un impact négligeable sur la qualité des images en raison du rapport de résolution des voies multi-spectrales avec la voie panchromatique. 35 Plus généralement, une valeur d'aberration différentielle autre qu'une défocalisation relative est volontairement introduite entre ensembles de barrettes choisis. 3055727 8 Selon un mode de réalisation, la détermination des aberrations optiques exploite uniquement les zones de recouvrement de chaque bande spectrale, en exploitant l'aberration différentielle existant entre images acquises par les barrettes paires et impaires. Avantageusement, il est possible alors de caractériser des aberrations optiques 5 par bande spectrale, qu'on appelle également aberrations chromatiques. De plus, la détermination des aberrations optiques est effectuée par zone de recouvrement, ce qui permet avantageusement de caractériser les aberrations optiques selon leur position spatiale dans le champ de prise de vue. Avantageusement, l'exploitation des zones de recouvrement de chaque bande 10 spectrale permet de caractériser les aberrations optiques avec une localisation spatiale et spectrale dans le champ de prise de vue du dispositif d'acquisition d'images. De manière connue, tout système optique est caractérisé par une fonction de transfert optique, qui est la transformation de Fourier de la réponse impulsionnelle optique ou PSF (de l'anglais « Point Spread Function »). 15 La fonction de transfert optique est une fonction complexe des fréquences spatiales, et l'application d'une transformation de Fourier à un signal d'image numérique défini dans un domaine spatial bidimensionnel résulte en un signal d'image transformé, défini dans le domaine-fréquentiel, chaque point ayant un module et une phase associés. La fonction de transfert optique caractérise le flou introduit par le système optique, 20 donc l'atténuation des fréquences spatiales. La fonction de transfert optique idéale, appelée également fonction de transfert optique de diffraction caractérise le flou introduit en l'absence de défauts dans la réalisation de la combinaison optique. Dans la pratique, des aberrations optiques sont introduites par le système optique, et par conséquent la fonction de transfert optique réelle est différente de la 25 fonction de transfert optique de diffraction. On modélise une image numérique I issue d'un signal d'image S représentatif d'une scène captée, formée par un dispositif d'acquisition d'images comportant un système optique de réponse impulsionnelle optique donnée, comme la convolution du signal d'image S observé avec la réponse impulsionnelle optique, puis l'échantillonnage 30 sur une grille bidimensionnelle et l'adjonction de bruit. Il est à noter que l'échantillonnage s'accompagne en général d'un repliement de spectre si la condition de Shannon n'est pas satisfaite. La réponse impulsionnelle h d'un dispositif d'acquisition d'images peut s'écrire comme un produit de convolution de deux termes : 35 h = hop, * hdet (EQ1) 3055727 9 Où l'opérateur * désigne le produit de convolution, hop, est la réponse impulsionnelle optique, et hdé, est la réponse impulsionnelle du détecteur, associée à l'intégration de photons par un détecteur élémentaire et à l'éventuel déplacement du détecteur faisant partie du dispositif d'acquisition d'images pendant le temps d'intégration. 5 L'image numérique 1 s'écrit donc, en utilisant un modèle convolutif d'acquisition d'images : S * hdet * hopt b (EQ2) où b est un terme de bruit. La réponse impulsionnelle optique et la FTM optique qui est sa transformée de 10 Fourier, sont entièrement déterminées par la connaissance de l'erreur front d'onde (WFE) et du masque (P) d'occultation de la pupille de l'objetif du système d'acquisition d'images, par exemple du téléscope dans le cas de l'acquisition d'images par satellite.. On cherche à résoudre le problème de la diversité de phase consistant à estimer l'erreur de front d'onde en utilisant au moins deux images acquises sur un même paysage 15 avec une aberration optique différentielle. La résolution du problème de diversité de phase repose sur une modélisation, selon le modèle convolutif d'acquisition d'images, du couple d'images observées, notées 11 et 12, en fonction des aberrations recherchées, elles-mêmes représentées sous la forme d'une combinaison de polynômes de Zernike, les inconnues étant les poids a, de cette 20 combinaison définis ci-après. Les relations suivantes sont appliquées : {TF (S)x FTM d' X FTM opt(WFE)tp +TF (bi) =TF (11) (EQ3) 'TF (S)x FTM d' X FTMopt+,,a,err(WFE)x RP2IiLep +TF (b2) TF (12 ) (EQ4) 25 Où TF() est la transformation de Fourier, FTMoe désigne la fonction de transfert optique (transformation de Fourier de hopt) et FTMdét est la transformation de Fourier de hdét RP2,1 est une fonction fréquentielle de type rampe de phase linéaire représentant dans le domaine de Fourier le décalage spatial relatif éventuellement non nul entre les deux images. 30 La valeur Aaberr est associée à l'aberration différentielle supplémentaire, introduite ou estimée, entre les deux images 11 et 12. Les termes ID, et b2 désignent le bruit correspondant aux images respectives 1,, 12. 3055727 10 L'expression {}rep indique la présence éventuelle de repliement dans le domaine fréquentiel, appelé également repliement spectral. Le phénomène de repliement spectral est connu. L'erreur de front d'onde WFE s'exprime dans la base orthonormée des polynômes 5 de Zernike par : WFE =WFEa_ priori +la i (EQ5) Où S représente l'ensemble d'indices de la base de polynômes de Zernike correspondant aux aberrations optiques recherchées, Z, est le polynôme de Zernike d'indice i, et a, des coefficients de décomposition (ou poids) correspondants. A chaque 10 polynôme de Zernike est associé un type d'aberration optique, par exemple Z4 représente l'erreur de focalisation, communément appelée défocus, Z5 et Z6 deux types d'astigmatisme etc. Le nombre d'indices dans S est quelconque, par exemple compris entre 1 et 40. Chaque valeur de coefficient ai de la décomposition (EQ5) est représentative de 15 l'aberration correspondante dans l'image considérée. Par exemple, la valeur du coefficient a4 est représentative de l'erreur de focalisation, elle représente l'erreur quadratique moyenne due à la défocalisation. Pour calculer les fonctions de transfert optiques (FTIVlopt) associées aux erreurs de front d'ondes (WFE et WFE + Aaberr), des deux images 1, et 12 on se base sur la relation 20 (EQ6) suivante : ( 2 ( 12z WFE'\ Pe FTM op, = TF TF-1 (EQ6) P étant le masque d'occultation et X la longueur d'onde centroïde de la bande spectrale considérée (moyenne pondérée par le gabarit spectral de la bande spectrale et 25 la réponse spectrale du paysage). Il est à noter que pour une scène complexe, par exemple un paysage observé par imagerie satellitaire, le signal d'image S est inconnu. Il en est de même du module de la fonction de transfert FTMdé, du détecteur. Dans ses divers modes de réalisation, l'invention propose d'utiliser une acquisition 30 d'un couple d'images numériques 1,, 12 représentatives d'une même portion de scène S, et présentant une aberration différentielle de valeur Aaberr, par exemple une différence de focalisation Aaberr, qui est déterminée par construction ou par estimation. 3055727 11 Pour la première image numérique 11, la relation (EQ3), en considérant le bruit négligeable, s'écrit : TTF (S) - FTM d' - FTM opt} rep - TF (I 1) (EQ7) Le terme WFE a été omis pour simplifier l'écriture des équations. 5 Considérant que la valeur d'aberration différentielle Aaberr s'applique sur la deuxième image 12 par rapport à la première image 11, la relation (EQ4) s'écrit : {TF (S) - FTM d' - FTM ,, opt+,,aberr} repTF (I 2) (EQ8) Il est à noter néanmoins que pour une deuxième image 12 présentant une aberration différentielle de valeur Aaberr importante, par exemple 200nm de 10 défocalisation, une quantité importante de flou est introduite dans cette deuxième image, et par conséquent cette deuxième image comprend peu de hautes fréquences spatiales et le repliement de spectre très faible. Avantageusement, l'invention propose de minimiser une fonctionnelle Fonc, obtenue à partir des relations (EQ7) et (EQ8), qui ne nécessite pas la connaissance du 15 signal S, ni de la fonction de transfert spécifique aux détecteurs FTMdét. La fonctionnelle est obtenue en éliminant entre les deux équations (EQ7) et (EQ8) le terme commun TF (S) - FTM de,. Dans un mode de réalisation, la fonctionnelle Fonc à minimiser est donnée par : 21,2. ff TF (12) FTM ot -TF(h) - p11) FTM - ' , Fane = opt+,_.ct,err df df y (EQ9) i 2 i FTM optl +FTM opt+Actberr + E 2 20 Où Gr 2 est la variance du bruit radiométrique b moyen de chacune des deux images, fx et fy représentent les fréquences spatiales, e un terme de correction évitant les instabilités numériques pour les fréquences rendant le dénominateur proche de 0. Par exemple, e peut prendre la valeur 10-12. 25 La variance cr 2 est par exemple estimée en appliquant un modèle de bruit du type : 62 = ldiA 2 ± B .Moy , où Moy est la moyenne de l'image 11, et A et B sont des valeurs réelles dépendant des caractéristiques instrumentales du dispositif d'acquisition d'images. L'utilisation du paramètre cr 2 est optionnelle puisque le critère d'arrêt de la 30 minimisation de la fonctionnelle porte sur la valeur relative de la fonctionnelle entre deux itérations. Néanmoins, la disponibilité d'un tel modèle permet avantageusement d'éliminer 3055727 12 les images trop bruitées et fournit un indice de confiance sur l'estimation, une bonne estimation correspondant à une valeur de fonctionnelle finale proche de 1. Avantageusement, l'expression (EQ9) de la fonctionnelle Fonc à minimiser ne comprend pas les termes TF(S) et FTMdet qui ne sont pas connus. 5 Plusieurs modes de réalisation sont envisagés pour le calcul des aberrations optiques à partir de la minimisation de la fonctionnelle Fonc. La figure 4 est un logigramme des principales étapes d'un procédé de caractérisation d'aberrations optiques selon un premier mode de réalisation de l'invention. Lors d'une première étape 50, des images, respectivement une première image 11 10 et une deuxième image 12, présentant une aberration différentielle et correspondant à une même scène, sont obtenues à partir d'un dispositif d'acquisition d'image. Comme déjà expliqué ci-dessus, dans un mode de réalisation cette aberration différentielle est présente dans le dispositif d'acquisition d'images, par exemple dans les zones de recouvrement inter-barrettes dans le cas d'un dispositif d'acquisition d'images 15 tel que décrit à la figure 3. Dans un autre mode de réalisation une aberration différentielle est introduite volontairement, et la valeur Aaberr correspondante est par conséquent connue. Dans cet autre mode, une valeur Aaberr d'aberration différentielle plus importante peut être introduite, mais la deuxième image 12 obtenue n'est dans ce cas pas exploitable 20 par une application utilisatrice, sauf dans le cas multispectral mentionné ci-dessus. Par exemple, lorsque l'aberration différentielle est une défocalisation, la deuxième image 12 présente une quantité de flou importante. La valeur d'aberration différentielle Aaberr est obtenue à l'étape 52. Dans le mode de réalisation dans lequel l'aberration différentielle est native dans 25 le dispositif d'acquisition d'images, à l'étape 52 on obtient une valeur Aaberr initiale d'aberration différentielle, une valeur affinée étant obtenue à l'étape d'estimation ultérieure décrite ci-dessous. Une transformation de Fourier est ensuite appliquée lors d'une étape de transformation 56 aux signaux d'image, pour obtenir un premier signal d'image transformé 30 TF(Ii)et un deuxième signal d'image transformé TF(I2) . L'étape de transformation 56 est suivie d'une étape de modélisation 58 du module de la fonction de transfert optique. Comme expliqué ci-dessus, on cherche ensuite à minimiser une fonctionnelle Fonc, qui dépend de TF(/1), de TF(I2) et des des fonctions de transfert optique 35 associées aux deux acquisitions et respectivement notées FTM op, et FTM opt+r,",err 3055727 13 Une fonction de transfert optique s'écrit, de manière connue, comme l'autocorrélation de la fonction pupillaire complexe Pc, relative à la pupille de l'objectif du système optique d'acquisition d'images. FTM op, = autocorr(Pa(x, y)) (EQ10) 5 En un point (x,y), la fonction pupillaire complexe s'écrit en fonction du module P(x,y) et de la phase v(x, y) : Pc(X, Y) = P(x, y)eiex,y) (EQ11) Le terme P(x,y) est le masque binaire qui vaut 1 sur la pupille et 0 ailleurs. Le terme de phase v(x, y) , appelé également phase aberrante, représente la 10 différence au point (x,y) entre la phase de l'onde lumineuse réelle sur le plan pupillaire et la phase d'une onde sphérique idéale. On caractérise l'erreur de front d'onde WFE(x,y) par la formule : ex, y) = -27c WFE(x, y) (EQ12) la longueur d'onde centroïde de la bande spectrale considérée. L'erreur de front 15 d'onde s'exprime dans la base orthonormée des polynômes de Zernike selon la formule (EQ5) rappelle ci-dessous WFE(x, y) = WFE a priori (X, y) ±IaiZ i(x, y) (EQ5) Chaque valeur de coefficient ai de la décomposition est représentative de l'aberration correspondante. Ainsi, pour des valeurs des coefficients ai données, on obtient l'erreur de front 20 d'onde selon la formule (EQ5), et en combinant les formules (EQ10), (EQ11) et (EQ12) on en déduit la fonction de transfert optique FTMopt. Dans un mode de réalisation, l'étape 60 met en oeuvre une estimation itérative, en faisant évoluer les valeurs respectives des coefficients a, jusqu'à minimiser la fonctionnelle Fonc. Par exemple, on utilise l'algorithme de Levenberg-Marquardt, ou tout 25 autre algorithme de minimisation itérative. Le critère d'arrêt des itérations porte sur la valeur relative de la fonctionnelle entre deux itérations : la convergence est réputée atteinte lorsque cette valeur relative est inférieure à un seuil prédéterminé. Dans un mode de réalisation, lors de l'étape d'estimation itérative, on estime 30 également la longueur d'onde centroïde X et le décalage spatial entre les deux images 11 et 12, représenté par le terme RP211 dans l'équation (EQ4), qui sont également des inconnues, en partant de valeurs initiales prédéterminées. 3055727 14 En outre, lorsque l'aberration différentielle est native dans le dispositif d'acquisition d'images, la valeur Aaberr de l'aberration différentielle est également estimée lors de l'étape d'estimation itérative à l'étape 60. La figure 5 est un logigramme des principales étapes d'un procédé de 5 caractérisation d'aberrations optiques selon un deuxième mode de réalisation de l'invention. Les étapes 62 à 68 sont analogues aux étapes 50 à 56 du premier mode de réalisation et ne sont pas décrites plus en détail ici. L'étape d'application d'une transformation de Fourier 68 est suivie d'une étape de 10 calcul des modules des fonctions de transfert optique respectifs notés FTM op, et FTM opt+,,cwerr basé sur une simplification des formules applicable lorsque la phase aberrante 9(x,y) associée à la représentation de la fonction pupillaire complexe est de faible amplitude. Dans ce cas, un développement limité, par exemple au second ordre, est 15 appliqué : e1O 1+ jv- v2 (EQ13) 2 La fonction de transfert optique peut alors s'écrire sous la forme d'un développement limité: 20 FTM api FTM (WFEa pria )+Iceicli+Ifluakai (EQ14) ieL k,leL Où FTMop, (WFEa priori ) est la fonction de transfert optique correspondant à la connaissance a priori de la WFE du système optique, les termes {ai} sont les coefficients de la décomposition de l'erreur de front d'onde sur la base des polynômes de Zernike, L 25 est le sous-ensemble des indices des polynômes de Zernike permettant de représenter la phase aberrante et les termes lai 1, Ifikil sont des fonctions fréquentielles. Ainsi, on modélise les fonctions de transfert optique associées aux deux acquisitions par deux polynômes de degré inférieur ou égal à 2 dont les variables a,, sont représentatives des aberrations recherchées. 30 Dans un cas particulier, lorsque la seule aberration optique recherchée est un défaut de défocalisation, seul le coefficient a4 correspondantau polynôme de Zernike Z4 est non nul, et la formule (EQ14) se réduit à : 3055727 15 FTM opt FTM opt (WFE prion)± a4 a4 (EQ15) L'expression (EQ14) ou l'expression (EQ15) est utilisée dans l'expression de la fonctionnelle à minimiser. L'étape 70 est suivie d'une étape 72 d'estimation des valeurs de a4 et de P4 5 permettant de minimiser la fonctionnelle Fonc, et par conséquent on obtient a4, valeur caractérisant l'aberration de défocalisation. La valeur de a4 obtenue est ensuite exploitée dans une étape ultérieure d'utilisation, par exemple pour corriger le réglage de la focalisation de prise de vue, ou pour déconvoluer finement l'image 12 par traitement algorithmique. 10 Avantageusement, l'utilisation de l'expression (EQ14) simplifie le calcul de la fonctionnelle et permet une minimisation non itérative, donc un allégement de la charge calculatoire. Les divers modes de réalisation du procédé de caractérisation d'aberrations optiques de l'invention sont mis en oeuvre par un dispositif programmable de type 15 ordinateur, station de travail, dont les principaux blocs fonctionnels sont illustrés à la figure 8. Un dispositif programmable 100 apte à mettre en oeuvre les procédés de l'invention comprend une unité centrale de traitement 108, ou CPU, apte à exécuter des instructions de programme informatique lorsque le dispositif 100 est mis sous tension. 20 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 100 comporte également des moyens de stockage d'informations 110, 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 25 caractérisation d'aberrations optiques d'un système optique selon l'invention. Le dispositif 100 comporte des moyens de commande 104 permettant de mettre à jour des paramètres et de recevoir des commandes d'un opérateur. Lorsque le dispositif programmable 100 est un dispositif embarqué, les moyens de commande 104 comprennent un dispositif de télécommunication permettant de recevoir des commandes 30 et des valeurs de paramètres à distance. Alternativement et de manière optionnelle, les moyens de commande 104 sont des moyens de saisie de commandes d'un opérateur, par exemple un clavier. De manière optionnelle, le dispositif programmable 100 comprend un écran 102 et un moyen supplémentaire de pointage 106, tel une souris. 35 Les divers blocs fonctionnels du dispositif 100 décrits ci-dessus sont connectés via un bus de communication 112. 3055727 16 Dans un mode de réalisation, le dispositif 100 est un ASIC, un FPGA ou un processeur. Un système de caractérisation des aberrations d'un système optique d'un dispositif d'acquisition d'images selon l'invention comporte un dispositif d'acquisition d'images, tel le 5 dispositif d'acquisition d'images 10 de la figure 1 et un dispositif apte à mettre en oeuvre le procédé de caractérisation des aberrations d'un système optique tel que décrit ci-dessus. Avantageusement, le procédé de l'invention décrit ci-dessus permet d'obtenir la caractérisation des aberrations optiques avec une valeur d'aberration différentielle, par exemple un décalage de focalisation natif ou introduit volontairement. 10 REVENDICATIONS1.- Procédé de caractérisation d'au moins une aberration optique d'un système optique d'un dispositif d'acquisition d'images, le système optique ayant une fonction de transfert optique associée dépendant des aberrations, le dispositif d'acquisition étant apte à acquérir au moins deux images dans un champ de prise vue de manière à introduire une aberration différentielle entre les deux images, chaque image étant définie par un signal numérique d'image, caractérisé en ce qu'il comporte des étapes de : -obtention (50, 62) d'une première image et d'une deuxième image d'une même zone du champ de prise de vue, la deuxième image étant acquise avec une aberration différentielle par rapport à la première image, chaque acquisition d'image ayant une fonction de transfert optique associée, -utilisation d'un modèle convolutif de formation d'images pour les images acquises utilisant une fonction de transfert optique dépendant d'au moins une aberration optique (Aaberr), -calcul (58, 60, 70, 72) d'au moins un paramètre caractérisant au moins une aberration optique par minimisation d'une fonctionnelle, définie en utilisant le modèle convolutif de formation d'images, exprimée en fonction des première et deuxième images acquises, d'une modélisation des fonctions de transfert optiques associées aux deux acquisitions dépendant de la ou des aberrations optique à caractériser ainsi que de l'aberration différentielle (Aaberr). 2.- Procédé selon la revendication 1, dans lequel le système optique à caractériser comprend au moins un ensemble de barrettes de détecteurs (401, 402,..., 40N) disposées en quinconce, et dans lequel lesdites première et deuxième images appartiennent à des zones de recouvrement entre barrettes de détecteurs adjacentes. 3.- Procédé selon la revendication 2, dans lequel le système optique comprend une pluralité d'ensembles de barrettes de détecteurs (401, 402,..., 40N), chaque ensemble de barrettes de détecteurs ayant une bande spectrale d'acquisition associée, le procédé étant appliqué sur l'ensemble des zones de recouvrement des différentes bandes spectrales correspondant à une même zone du champ de vue, afin de déterminer des aberrations relatives à chacune de ces zones du champ de vue. 3055727 18 4.- Procédé selon la revendication 2, appliqué sur au moins une zone de recouvrement de chaque bande spectrale, de manière à caractériser des aberrations optiques par bande spectrale. 5 5.- Procédé selon l'une quelconque des revendications 2 à 4, dans lequel ladite fonctionnelle dépend en outre d'une longueur d'onde centroïde associée à la bande spectrale d'acquisition, et dans lequel ladite longueur d'onde centroïde est estimée lors de la minimisation de ladite fonctionnelle. 10 6.- Procédé selon l'une quelconque des revendications 1 à 5, dans lequel ladite fonctionnelle dépend en outre d'un décalage spatial entre la première image et la deuxième images acquises, ledit décalage spatial étant estimé lors de la minimisation de ladite fonctionnelle. 15 7.- Procédé selon l'une des revendications 1 à 6, dans lequel l'étape de calcul d'au moins un paramètre caractérisant au moins une aberration optique met en oeuvre une modélisation de ladite fonction de transfert optique en fonction d'au moins un coefficient d'une décomposition dans une base de polynômes de Zernike d'une erreur de front d'onde associée au système optique du dispositif d'acquisition. 20 8.- Procédé selon l'une quelconque des revendications 1 à 7, dans lequel la minimisation de la fonctionnelle est effectuée par une méthode d'estimation itérative. 9.- Procédé selon l'une quelconque des revendications 1 à 8, comprenant une 25 modélisation de chacune des fonctions de transfert optique associées aux acquisitions d'image par un polynôme de degré inférieur ou égal à deux, défini par des variables représentatives des aberrations optiques à caractériser et une minimisation de ladite fonctionnelle par une méthode non itérative. 30 10.- Procédé selon l'une quelconque des revendications 1 à 9, dans lequel une valeur de l'aberration différentielle entre la première et la deuxième image acquise est estimée lors de la minimisation de ladite fonctionnelle. 11.- Procédé selon l'une quelconque des revendications 1 à 10, dans lequel ladite aberration différentielle entre la première image et la deuxième image est une défocalisation relative entre la première image et la deuxième image. 3055727 19 12.- Programme d'ordinateur comportant des instructions logicielles qui, lorsqu'elles sont mises en oeuvre par un dispositif programmable, mettent en oeuvre un procédé de caractérisation d'aberrations optiques d'un dispositif d'acquisition d'images 5 selon les revendications 1 à 11. 13.- Dispositif de caractérisation d'au moins une aberration optique d'un système optique d'un dispositif d'acquisition d'images, le système optique ayant une fonction de transfert optique associée dépendant des aberrations, le dispositif d'acquisition étant apte 10 à acquérir au moins deux images dans un champ de prise vue de manière à introduire une aberration différentielle entre les deux images, chaque image étant définie par un signal numérique d'image, caractérisé en ce qu'il comporte : -module d'obtention d'une première image et d'une deuxième image d'une même zone du champ de prise de vue, acquises par le dispositif d'acquisition d'images, la 15 deuxième image étant acquise avec une aberration différentielle par rapport à la première image, chaque acquisition d'image ayant une fonction de transfert optique associée, -un processeur adapté à mettre en oeuvre un calcul d'au moins un paramètre caractérisant au moins une aberration optique par minimisation d'une fonctionnelle exprimée en fonction des première et deuxième images acquises, utilisant un modèle 20 convolutif de formation d'images pour les images acquises utilisant une fonction de transfert optique dépendant d'au moins une aberration optique, ladite fonctionnelle étant également exprimée en fonction d'une modélisation des fonctions de transfert optiques associées aux deux acquisitions dépendant de la ou des aberrations optique à caractériser ainsi que de l'aberration différentielle. 25