Brevet FR3147415A1 — publié le 20241004

Procédé (100) de mesure d’une aberration dans une image acquise par une matrice de détecteurs, et une pluralité de matrices de filtres spectraux, comprenant les étapes suivantes : acquisition (101) d’une source ponctuelle connue; détermination (102) d’une valeur initiale (VIP) des paramètres d’un modèle d’une réponse impulsionnelle; échantillonnage (103), calcul (104) d’un coût fonction d’une distance entre la réponse impulsionnelle calculée et la réponse impulsionnelle observée ; calcul (105) d’un coût global (CG) pour l’ensemble des bandes spectrales; estimation (102bis) d’une valeur (VEP) de chaque paramètre; échantillonnage (103bis); calcul (104bis) du coût; calcul (105bis) du coût global (NCG) ; comparaison (106) entre le coût global (NCG) et le coût global (CG);répéter (107). Figure 1 Procédé de mesure et/ou de compensation d’une aberration dans une image acquise par un instrument optique comprenant une matrice de détecteurs multispectraux. La présente invention concerne le domaine de l’étalonnage d’un instrument optique comprenant un capteur matriciel, en particulier un capteur optique matriciel multi spectral, par observation d’une source ponctuelle connue, nécessaire notamment pour la compensation d’une ou plusieurs aberrations optiques de l’instrument optique dans une image acquise par l’instrument optique. Il est connu de compenser une ou plusieurs aberrations optiques de l’instrument optique dans une image acquise par un instrument optique comprenant une matrice de détecteurs, à l’aide d’une modélisation paramétrique d’une réponse impulsionnelle de l’instrument, en prenant en compte les différents contributeurs à la formation de l’image par l’instrument, i.e. une combinaison optique, un détecteur, un déplacement de l’instrument durant le temps de pose, etc... Les procédés connus ne traitent cependant pas le cas d’un instrument optique comprenant une matrice de détecteurs multispectraux. L’invention a donc pour but de proposer une solution à tout ou partie de ces problèmes. A cet effet, la présente invention concerne un procédé de mesure d’une aberration dans une image acquise par un instrument optique comprenant une matrice de détecteurs multispectraux, la matrice de détecteurs multispectraux comprenant une matrice de détecteurs, et une pluralité de matrices de filtres spectraux, chaque matrice de filtre spectral étant associée à une bande spectrale, le procédé comprenant les étapes suivantes : - acquisition par la matrice de détecteurs de l’instrument optique d’une image d’une source ponctuelle connue, l’image comprenant, pour chaque matrice de filtre spectral, une partie d’une tâche image de la source ponctuelle acquise dans la bande spectrale de ladite matrice de filtre spectral ; - pour chaque matrice de filtre spectral, détermination d’une valeur initiale de chaque paramètre d’un ensemble de paramètres d’un modèle paramétrique d’une réponse impulsionnelle modélisée de l’instrument optique, l’ensemble de paramètres comprenant un sous-ensemble de paramètres commun à la pluralité de matrices de filtres spectraux, et un sous-ensemble de paramètres spécifiques à chaque matrice de filtre spectral; - pour chaque matrice de filtre spectral, échantillonnage d’une réponse impulsionnelle calculée sur la base de la valeur initiale des paramètres appliquée au modèle paramétrique, l’échantillonnage étant réalisé selon un masque d’échantillonnage associé à la matrice de filtre spectral, - pour chaque matrice de filtre spectral, calcul d’une valeur d’un coût dans la bande spectrale, la valeur du coût étant fonction d’une distance entre la réponse impulsionnelle calculée et la réponse impulsionnelle observée pour chacun des points du masque d’échantillonnage, - calcul d’une valeur d’un coût global pour l’ensemble des bandes spectrales en fonction de la valeur du coût dans chaque bande spectrale, - pour chaque matrice de filtre spectral, estimation d’une valeur estimée de chaque paramètre de l’ensemble de paramètres, la valeur estimée étant fonction d’une variation estimée à partir de la valeur du coût global, la valeur estimée étant égale à la variation estimée ajoutée à la valeur initiale des paramètres; - pour chaque matrice de filtre spectral, échantillonnage d’une nouvelle réponse impulsionnelle calculée sur la base de la valeur estimée des paramètres appliquée au modèle paramétrique, l’échantillonnage étant réalisé selon le masque d’échantillonnage associé à la matrice de filtre spectral, - pour chaque matrice de filtre spectral, calcul d’une nouvelle valeur du coût dans ladite bande spectrale, la nouvelle valeur du coût étant fonction de la distance entre la nouvelle réponse impulsionnelle calculée et la réponse impulsionnelle observée pour chacun des points du masque d’échantillonnage, - calcul (105bis) d’une nouvelle valeur du coût global pour l’ensemble des bandes spectrales en fonction de la nouvelle valeur du coût dans chaque bande spectrale. - comparaison entre la nouvelle valeur du coût global et la valeur du coût global pour obtenir un résultat de la comparaison; - déterminer que la valeur initiale des paramètres est égale la valeur estimée des paramètres et que la valeur du coût global est égale à la nouvelle valeur du coût global, et pour chaque matrice de filtre spectral, répéter les étapes comprises entre l’étape d’estimation d’une valeur estimée des paramètres, jusqu’à l’étape de comparaison entre la nouvelle valeur du coût global et la valeur du coût global jusqu’à ce qu’un critère d’arrêt soit satisfait par le résultat de la comparaison, de sorte qu’à chaque itération de l’étape de répétition, au cours de l’étape d’estimation la variation est estimée à partir d’une valeur initiale des paramètres égale à la valeur estimée des paramètres estimée au cours de l’itération précédente, et à partir d’une valeur du coût global égale à une nouvelle du coût global calculée au cours de l’itération précédente; - déterminer une mesure de l’aberration sur la base de la valeur estimée des paramètres de la réponse impulsionnelle. Selon un mode de mise en œuvre, l’invention comprend une ou plusieurs des caractéristiques suivantes, seules ou en combinaison techniquement acceptable. Selon un mode mise en œuvre, la source ponctuelle est une étoile observée hors de l’atmosphère terrestre, une réponse spectrale de l’étoile observée étant connue. Selon un mode mise en œuvre, la valeur initiale des paramètres est égale à 0 pour chaque paramètre. Selon un mode de mise en œuvre, chaque matrice de filtre spectral est superposée à la matrice de détecteurs, une position de la matrice de filtre spectral relativement à la matrice de détecteurs étant décalée selon au moins l’une parmi les deux dimensions de la matrice de détecteurs, de sorte que chaque point d’échantillonnage du masque d’échantillonnage associé à une matrice de filtre spectral est décalé par rapport à un point correspondant du masque d’échantillonnage associé à une autre matrice de filtre spectral. Selon ces dispositions, plus le nombre de matrices de filtres spectraux est important plus l’estimation de la valeur des paramètres communs du modèle paramétrique est robuste, et meilleure est la qualité de la compensation de l’aberration obtenue grâce au modèle. Selon un mode mise en œuvre, l’étape d’acquisition comprend l’acquisition de plusieurs sources ponctuelles connues et/ou plusieurs acquisitions d’une même source ponctuelle connue, l’ensemble de paramètres du modèle paramétrique comprenant un sous-ensemble de paramètres spécifiques à chaque acquisition, le modèle paramétrique de réponse impulsionnelle étant associé aux différentes réponses impulsionnelles observées dans la bande spectrale de chaque matrice de filtre spectral, et pour chaque source ponctuelle et/ou pour chaque acquisition d’une même source ponctuelle, et dans lequel les étapes suivantes sont mises en œuvre pour chaque matrice de filtre spectral, et pour chaque source ponctuelle, et/ou pour chaque acquisition d’une même source ponctuelle : - étape de détermination de la valeur initiale des paramètres, et étape d’estimation de la valeur estimée des paramètres, - étapes d’échantillonnage de la réponse impulsionnelle calculée sur la base de la valeur initiale des paramètres, respectivement de la nouvelle réponse impulsionnelle calculée sur la base de la valeur estimée des paramètres, - étapes de calcul de la valeur du coût, respectivement de la nouvelle valeur du coût, et dans lequel la valeur du coût global, respectivement la nouvelle valeur du coût global, est calculée pour l’ensemble des bandes spectrales et pour l’ensemble des sources ponctuelles. Selon un mode de mise en œuvre, parmi la pluralité d’acquisitions successives d’une source ponctuelle, au moins un paramètre du modèle paramétrique prend une valeur au cours d’une acquisition et une autre valeur différente au cours d’au moins une autre acquisition, la valeur et l’autre valeur dudit au moins un paramètre étant inconnues lors de l’étape d’acquisition. La variation appliquée entre la valeur et l’autre valeur n’est pas nécessairement connue. Si la variation du paramètre est connue, l’estimation de la valeur des paramètres du modèle paramétrique sera plus robuste. Selon un mode de mise en œuvre, la pluralité de matrices de filtres spectraux comprend au moins trois matrices de filtres spectraux, configurées pour filtrer chacune dans une bande spectrale complémentaire des autres matrices. Selon un mode de mise en œuvre, la pluralité de matrices de filtres spectraux est du type matrice de Bayer. Selon un mode de mise en œuvre, à l’étape d’estimation de la valeur estimée de chaque paramètre, la variation est estimée à partir d’une matrice jacobienne du modèle paramétrique de la réponse impulsionnelle, de sorte que la répétition, des étapes comprises entre l’étape d’estimation d’une valeur estimée des paramètres, et l’étape de comparaison entre la nouvelle valeur du coût global et la valeur du coût global, jusqu’à ce qu’un critère d’arrêt soit satisfait par le résultat de la comparaison, correspond à une mise en œuvre d’un algorithme itératif du type Levenberg-Marquardt. Selon un mode de mise en œuvre, une amplitude de la variation estimée est contrainte, de sorte que la variation estimée est comprise entre une valeur minimum et une valeur maximum. Selon un mode de mise en œuvre, la valeur du coût global, respectivement la nouvelle valeur du coût global, est égale à une somme de la valeur des coûts, respectivement de la nouvelle valeur des coûts, sur l’ensemble des bandes spectrales. Selon un mode de mise en œuvre, le résultat de la comparaison est égal à une différence entre la valeur du nouveau coût global et la valeur du coût global, le critère d’arrêt prédéterminé étant satisfait lorsque le résultat de la comparaison est inférieur à un seuil. Selon un mode de mise en œuvre, l’ensemble des paramètres comprend au moins l’un parmi les paramètres d’une réponse impulsionnelle représentative d’une partie optique de l’instrument optique, les paramètres d’une réponse impulsionnelle représentative de l’intégration des charges sur au moins un détecteur élémentaire de la matrice de détecteurs, les paramètres d’une réponse impulsionnelle représentative d’un filé, i.e. d’un mouvement linéaire, de la matrice de détecteurs pendant un temps de pose, les paramètres d’une réponse impulsionnelle représentative d’un bougé, i.e. d’un mouvement non linéaire, de la matrice de détecteurs pendant le temps de pose, les paramètres représentatifs des coordonnées d’une image géométrique de la source ponctuelle, ladite image géométrique étant ponctuelle.. Selon un mode de mise en œuvre, le procédé comprend une étape dans laquellela mesure de l’aberration est utilisée pour compenser l’aberration par un traitement d’amélioration de l’image acquise par l’instrument optique, le traitement d’amélioration étant fonction de la mesure de l’aberration, le traitement comprenant par exemple une application d’un filtre de compensation calculé à partir de la mesure de l’aberration, ou le traitement comprenant par exemple une restauration de l’image fonction de la mesure de l’aberration Selon un mode de mise en œuvre, le traitement est une commande de l’instrument optique pour actionner une optique active de l’instrument optique, la commande étant basée sur la valeur estimée des paramètres du modèle paramétrique. Selon un mode de mise en œuvre, la commande est l’application d’une consigne de correction des aberrations optiques à des moyens de correction, par exemple en agissant sur un ou plusieurs degrés de libertés d’un miroir secondaire de la combinaison optique ou sur un miroir déformable de la combinaison optique. A titre d’exemple, la refocalisation peut consister à translater le miroir secondaire le long d’un axe optique d’une quantité déterminée à partir de la mesure de l’aberration selon le procédé. Selon un mode de mise en œuvre, la commande de l’instrument optique est une commande de focalisation. Selon un aspect l’invention concerne également un programme d’ordinateur comprenant un ensemble d’instructions configuré pour mettre en œuvre le procédé selon l’un des modes de mise en œuvre décrit ci-avant lorsqu’il est exécuté par un processeur d’une unité de traitement des images et de contrôle, d’un instrument optique comprenant une matrice de détecteurs multispectraux, la matrice de détecteurs multispectraux comprenant une matrice de détecteurs, et une pluralité de matrices de filtres spectraux, chaque matrice de filtre spectral étant associée à une bande spectrale. Selon un autre aspect l’invention concerne également un support lisible par un ordinateur comprenant le programme d’ordinateur décrit ci-avant. Pour sa bonne compréhension, un mode de réalisation et/ou de mise en oeuvre de l’invention est décrit en référence aux dessins ci-annexés représentant, à titre d’exemple non limitatif, une forme de réalisation ou de mise en œuvre respectivement d’un dispositif et/ou d’un procédé selon l’invention. Les mêmes références sur les dessins désignent des éléments similaires ou des éléments dont les fonctions sont similaires. est une représentation schématique de la séquence des étapes du procédé selon un exemple de mise en œuvre de l’invention. Le procédé selon l’invention est basé sur une optimisation combinée de paramètres d’un modèle d’une réponse impulsionnelle modélisée à partir d’un ensemble d’acquisition d’une source ponctuelle connue, par exemple une étoile dans le cas d’observations hors de l’atmosphère terrestre. A partir d’une modélisation fine de la chaîne d’acquisition, i.e de l’instrument optique qui comprend notamment l’optique, le détecteur, les filtres multispectraux, et leur agencement, les paramètres du modèle sont affinés de manière itérative. Les paramètres comprennent des paramètres globaux relatifs à un état global de la chaîne d’acquisition, et des paramètres plus spécifiques à chaque bande spectrale. Selon un exemple de réalisation de la matrice de détecteurs multispectraux, la pluralité de matrices de filtres spectraux comprend au moins trois matrices de filtres spectraux, configurées pour filtrer chacune dans une bande spectrale complémentaire des autres matrices. L’homme du métier comprendra que chaque filtre n’est pas strictement limité à une longueur d’onde, i.e. chaque filtre a une transmission faible mais non nulle dans des longueurs d’ondes voisines autour de la longueur d’onde centrale du filtre. Par ailleurs, le procédé selon l’invention s’applique également à une combinaison comprenant une bande panchromatique, et un ou plusieurs filtres mono-spectraux. Selon un exemple particulier de réalisation, la pluralité de matrices de filtres spectraux est du type matrice de Bayer. Une valeur des paramètres du modèle est estimée en recherchant la valeur qui minimise une fonction de coût, fonction d’un écart entre la réponse impulsionnelle modélisée et la réponse impulsionnelle observée . La recherche de la valeur des paramètres qui minimise la fonction de coût est réalisée en utilisant un algorithme itératif, à partir d’une valeur initiale de chaque paramètre de l’ensemble de paramètres, avec un critère d’arrêt de l’algorithme itératif. L’algorithme itératif est par exemple un algorithme de Levenberg Marquardt. Ainsi, pour chaque observation d’une source ponctuelle et pour chaque bande spectrale, est estimée une valeur des paramètres, pour chaque paramètre de l’ensemble des paramètres considérés, qui minimise ladite fonction de coût. La valeur des paramètres ainsi obtenues est utilisée pour calculer, à partir du modèle paramétrique, une valeur de la réponse impulsionnelle modélisée pour l’ensemble des pixels de la matrice de détecteurs. Pour chaque bande spectrale, la réponse impulsionnelle modélisée sur la base de la valeur estimée des paramètres est ensuite échantillonnée selon un masque d’échantillonnage associé à la matrice de filtre spectral. Un écart entre la réponse impulsionnelle modélisée et la réponse impulsionnelle observée dans ladite bande spectrale est calculé sur la base d’une différence entre la réponse impulsionnelle modélisée et la réponse impulsionnelle observée pour chacun des points du masque d’échantillonnage pour ladite bande spectrale. Autrement dit, pour calculer l’écart, la réponse impulsionnelle modélisée est échantillonnée selon le masque d’échantillonnage de chaque filtre spectral, et la différence entre la réponse impulsionnelle modélisée et la réponse impulsionnelle observée n’est calculée que pour ces points du masque. Chaque point n’est analysé qu’une fois. Un coût, fonction de l’écart ainsi calculé, est déterminé pour chaque bande spectrale, et un coût global, pour l’ensemble des bandes spectrales considérées, est calculé sur la base de l’ensemble des écarts calculés pour chaque bande spectrale. Les étapes ci-dessus sont répétées jusqu’à ce qu’une différence entre la valeur du coût global et la valeur suivante du coût global soit inférieur à un seuil prédéterminé, autrement dit jusqu’à ce que la valeur du coût global s’approche suffisamment d’une valeur limite indéterminée. Le procédé selon un exemple de mise en œuvre de l’invention, tel que décrit ci-avant, permet d’améliorer une précision de détermination des paramètres physiques de l’instrument optique, notamment de la combinaison optique ; il permet également d’étalonner les dispositifs d’actionnement et de correction de l’instrument en simplifiant considérablement les conditions opérationnelles : la mise en œuvre du procédé nécessite uniquement quelques images acquises par l’instrument de points sources connus, par exemple des étoiles. Contrairement aux méthodes actuelles (de refocalisation, d’étalonnage absolue, de mesure de la FTM) qui nécessitent soit l’observation de cibles variées, soit des contraintes opérationnelles fortes (besoin impératif de défocalisation volontaire de l’instrument pour acquisitions à différents points de fonctionnement), la méthode proposée ne nécessite que quelques acquisitions sur des étoiles connues : il n’y a pas besoin d’avoir plusieurs points de focalisation, et si plusieurs points de focalisation sont utilisés, il n’est pas nécessaire de connaître précisément la variation appliquée entre ces points de focalisation. Dans le cas d’une pupille optique symétrique, il y a ambiguité de signe sur la partie symétrique des aberrations optiques. Par exemple, dans le cas du defocus, on ne sait pas dans quel sens déplacer le miroir secondaire. Cette ambiguité est levée si on dispose d’acquisitions pour lesquelles on a volontairement introduit un defocus connu, ou si la pupille n’est pas symétrique , ce qui est le cas usuel en raison de la présence de « bras » radiaux dûs aux supports physiques dudit miroir secondaire. Le procédé permet par exemple de compenser une distorsion dans une image acquise par un instrument optique, l’instrument optique comprenant notamment une combinaison optique et une matrice de détecteurs multispectraux, la combinaison optique étant agencée pour former une image d’une scène observée sur la matrice de détecteurs multispectraux placée dans un plan focal de la combinaison optique. Le traitement de compensation de la distorsion est, selon un exemple de mise en œuvre, une transformation de l’image par application d’une matrice de transformation géométrique, ladite matrice de transformation géométrique étant une matrice de passage de coordonnées angulaires d’un point de la scène observée à des coordonnées dans un référentiel cible. La matrice de transformation est déterminée à partir de paramètres qui permettent de modéliser une réponse impulsionnelle de l’instrument, de manière à compenser le ou les différents facteurs de distorsion introduits par les différents composants de l’instrument optique. Le procédé permet par exemple également compenser la distorsion par une commande d’actionnement d’une optique active de l’instrument optique, en vue d’arriver à état d’aberration optique cible après actuation. Ceci est applicable à la refocalisation de l’instrument optique et participe ainsi à l’étalonnage d’une optique active. Le procédé permet par exemple également l’estimation de paramètres géométriques (position/rotation) de miroirs de la combinaison optique, et plus généralement de tout paramètre intervenant dans le modèle optique ; Le procédé permet par exemple également la mesure de la FTM pour calibrer les algorithmes de restauration, ou le calcul du coefficient d’étalonnage absolue permettant le passage des valeurs de compte numériques en luminance équivalente. Le procédé permet par exemple également une ré-estimation d’une FTM détecteur qui traduit l’effet d’intégration des charges sur un détecteur élémentaire (diffusion des charges élargissant la PSF) pour améliorer le modèle de la chaîne d’acquisition. Le procédé permet par exemple également une estimation du filé (mouvement linéaire) pour estimer les dérives de pointage du satellite ; Le procédé permet par exemple également une estimation du bougé (reliquat non linaire des déplacements) pour estimer les vibrations. Le procédé 100 selon un exemple de mise en œuvre de l’invention va maintenant être décrit en référence à la . Le procédé comprend une étape d’acquisition 101, par la matrice de détecteurs de l’instrument optique, d’une image d’une ou de plusieurs sources ponctuelles connues ; notamment, le spectre d’émission et la magnitude apparente de chacune de ces étoiles est connu. Selon un exemple de mise en œuvre l’étape d’acquisition 101 peut comprendre également plusieurs acquisitions d’une même source ponctuelle connue. La ou les sources ponctuelles sont par exemple des étoiles, une réponse spectrale de chacune des étoiles observées étant connue. L’image comprend, pour chaque source ponctuelle, et pour chaque matrice de filtre spectral, une tâche image acquise dans la bande spectrale de ladite matrice de filtre spectral ; pour chaque source ponctuelle, l’image de la source ponctuelle est étalée sur un ensemble de pixels adjacents appartenant à des bandes spectrales potentiellement différentes. Il faut distinguer la réponse impulsionnelle instrumentale monochromatique, image d’une source ponctuelle monochromatique de longueur d’onde donnée et d’intégrale unité (distribution de Dirac spatiale et spectrale) et l’image d’une source ponctuelle polychromatique qui sera une intégrale spectrale des réponses impulsionnelles monochromatiques pondérées par la réponse spectrale des étoiles. Autrement dit, la tâche image d’une source ponctuelle diffère d’une bande spectrale à une autre bande spectrale. La tâche image d’une source ponctuelle, dans chaque bande spectrale, correspond à une réponse impulsionnelle de l’instrument optique observée dans ladite bande spectrale. La réponse impulsionnelle de l’instrument optique dans chaque bande spectrale est modélisée, d’une manière connue de l’homme du métier, comme un produit de convolution des réponses impulsionnelles dans ladite bande spectrale des différents éléments de la chaîne d’acquisition constitutive de l’instrument optique considéré. L’instrument optique, comprend en particulier l’optique de l’instrument optique, les détecteurs élémentaires de la matrice de détecteurs, sans oublier le cas échéant les autres contributeurs à la réponse impulsionnelle de l’instrument tels que d’éventuels déplacements pendant le temps de pause (filé et/ou bougé). Le modèle paramétrique de la réponse impulsionnelle est ainsi une fonction d’un ensemble de paramètres physiques caractéristiques des différents contributeurs de la chaîne d’acquisition, i.e. de l’instrument optique. L’ensemble de paramètres comprend des paramètres globaux, tel que par exemple les coefficients des polynomes de Zernike, caractéristiques d’un état global de la combinaison optique, et des paramètres spécifiques à certains éléments de la chaîne d’acquisition de l’instrument optique, notamment spécifiques à chaque bande spectrale et à chaque acquisition. Ainsi, par exemple, l’ensemble des paramètres peut comprendre les paramètres d’une réponse impulsionnelle représentative d’une pupille de l’instrument optique, et/ou les paramètres d’une réponse impulsionnelle représentative de l’intégration des charges sur au moins un détecteur élémentaire de la matrice de détecteurs, et/ou les paramètres d’une réponse impulsionnelle représentative d’un filé, i.e. d’un mouvement linéaire, de la matrice de détecteurs pendant un temps de pause, et/ou les paramètres d’une réponse impulsionnelle représentative d’un bougé, i.e. d’un mouvement non linéaire, de la matrice de détecteurs pendant le temps de pause, et/ou les paramètres représentatifs des coordonnées d’une image géométrique de la source ponctuelle, ladite image géométrique étant ponctuelle.. Pour chaque matrice de filtre spectral, optionnellement pour chaque source ponctuelle, l’étape suivante du procédé 100 est une étape de détermination 102 d’une valeur initiale VIP de chaque paramètre de l’ensemble des paramètres du modèle paramétrique de la réponse impulsionnelle modélisée de l’instrument optique, l’ensemble de paramètres comprenant en particulier un sous-ensemble de paramètres commun à la pluralité de matrices de filtres spectraux, et un sous-ensemble de paramètres spécifiques à chaque matrice de filtre spectral. Selon un exemple de mise en œuvre du procédé 100, l’étape de détermination 102 d’une valeur initiale des paramètres du modèle paramétrique, comprend le choix d’une valeur égale à 0, par exemple, pour chaque paramètre. Selon un autre exemple de mise en œuvre, l’étape de détermination 102 d’une valeur initiale comprend le choix, pour chaque paramètre, d’une valeur égale à une valeur déterminée antérieurement, par exemple par une calibration précédente. L’étape suivante du procédé 100, après l’étape de détermination 102 d’une valeur initiale des paramètres du modèle paramétrique, à calculer, pour chaque matrice de filtre spectral, et le cas échéant pour chaque source ponctuelle, une réponse impulsionnelle, sur la base de la valeur initiale VIP des paramètres, puis à échantillonner 103 la réponse impulsionnelle calculée, selon un masque d’échantillonnage associé à la matrice de filtre spectral: autrement dit, la valeur d’un échantillon de la réponse impulsionnelle est calculée, sur la base de la valeur des paramètres déterminés précédemment, en un point du masque d’échantillonnage, ledit masque d’échantillonnage comprenant un ensemble de points répartis régulièrement sur la matrice de filtre spectral considérée, par exemple au point central de chaque filtre élémentaire de la matrice de filtre spectral. Selon un exemple particulier, chaque matrice de filtre spectral est superposée à la matrice de détecteurs, et une position de la matrice de filtre spectral relativement à la matrice de détecteurs est décalée selon au moins l’une parmi les deux dimensions de la matrice de détecteurs, de sorte que chaque point d’échantillonnage du masque d’échantillonnage associé à une matrice de filtre spectral est décalé par rapport à un point correspondant du masque d’échantillonnage associé à une autre matrice de filtre spectral. Pour chaque matrice de filtre spectral, et le cas échéant pour chaque source ponctuelle, une valeur d’un coût, dans la bande spectrale considérée, fonction d’une distance entre la réponse impulsionnelle calculée, et la réponse impulsionnelle observée, est ensuite calculée, au cours d’une étape de calcul 104 d’une valeur d’un coût. Ce coût est calculé, par exemple, sur la base d’une différence entre la réponse impulsionnelle modélisée et la réponse impulsionnelle observée pour chacun des points du masque d’échantillonnage. Par exemple, le coût est égal à une somme des carrés des différences, pour chacun des points du masque d’échantillonnage, optionnellement pondéré par un coefficient de normalisation. Une valeur d’un coût global CG pour l’ensemble des bandes spectrales, le cas échéant pour chaque source ponctuelle, est ensuite calculé au cours d’une étape de calcul 105 d’une valeur du coût global CG. Par exemple, le coût global est égal à une somme des coûts précédemment calculés pour chaque bande spectrale. Pour chaque matrice de filtre spectral, et le cas échéant pour chaque source ponctuelle, une valeur estimée de chaque paramètre VEP de l’ensemble de paramètres est estimée au cours d’une étape d’estimation 102bis, la valeur estimée de chaque paramètre VEP étant estimée à partir de la valeur initiale VIP des paramètres, déterminée à l’étape 102, et d’une variation VAR; ladite variation VAR étant estimée en fonction de la valeur du coût global CG, la valeur estimée VEP étant égale à la variation estimée VAR ajoutée à la valeur initiale des paramètres VIP. Pour chaque matrice de filtre spectral, et le cas échéant pour chaque source ponctuelle, les étapes 103, 104, et 105, précédemment décrites, sont mises en œuvre à nouveau avec la valeur estimée des paramètres VEP. Ces étapes sont représentées sur la , respectivement avec les références 103bis, 104bis, 105bis, de sorte que l’étape 105bis est une étape de calcul d’une nouvelle valeur du coût global NCG pour l’ensemble des bandes spectrales, le cas échéant pour chaque source ponctuelle, pour la valeur estimée des paramètres VEP. La nouvelle valeur du coût global NCG calculé à l’étape 105bis pour la valeur estimée des paramètres VEP est comparé, au cours d’une étape de comparaison 106, à la valeur du coût global CG calculée à l’étape 105 pour la valeur initiale des paramètres VIP déterminée à l’étape 102. Si le résultat de la comparaison RC satisfait à un critère d’arrêt CRA alors la valeur estimée des paramètres VEP est retenue pour la modélisation de la réponse impulsionnelle de l’instrument optique. A titre d’exemple, le résultat de la comparaison RC est une différence entre la nouvelle valeur du coût global NCG et la valeur du coût global CG, et le critère d’arrêt CRA est satisfait si ladite différence est inférieure à un seuil prédéterminé. Si le résultat de la comparaison RC ne satisfait pas au critère d’arrêt CRA alors les étapes du procédé 100 sont répétées 107, depuis l’étape d’estimation 102bis d’une valeur des paramètres VEP, en passant par l’étape 104bis de calcul d’une nouvelle valeur du coût par bande spectrale, et par l’étape 105bis de calcul d’une nouvelle valeur du coût global NCG pour l’ensemble des bandes spectrales, jusqu’à l’étape de comparaison 106 entre la valeur du coût global CG, calculée à l’itération précédente de la répétition 107, et la nouvelle valeur du coût global NCG calculée à l’itération courante de la répétition 107, jusqu’à ce que le critère d’arrêt CRA soit satisfait par le résultat de la comparaison RC. A chaque nouvelle itération, l’étape d’estimation 102bis d’une valeur estimée des paramètres VEP est réalisée à partir d’une valeur des paramètres égale à la valeur estimée des paramètres VEP obtenue à l’étape d’estimation 102bis de l’itération précédente. A cet effet, à chaque nouvelle itération, au cours d’une étape de détermination 107bis, la valeur initiale des paramètres VIP devient égale à valeur estimée des paramètres VEP de l’itération précédente, et la valeur du coût global CG devient égale à la nouvelle valeur du coût global NCG, pour chaque matrice de filtre spectral. Selon un exemple de mise en œuvre, à chaque nouvelle itération, au cours de l’étape d’estimation 102bis d’une valeur estimée des paramètres VEP, la variation VAR de la valeur estimée des paramètres VEP estimée à l’itération précédente, est estimée à partir d’une matrice jacobienne du modèle paramétrique de la réponse impulsionnelle, ladite matrice jacobienne étant déterminée par rapport à la valeur estimée des paramètres VEP estimée à l’itération précédente, de sorte que la répétition 107, des étapes comprises entre l’étape d’estimation 102bis d’une valeur estimée des paramètres VEP, et l’étape de comparaison 106 entre la nouvelle valeur du coût global NCG et la valeur du coût global CG, jusqu’à ce qu’un critère d’arrêt CRA soit satisfait par le résultat de la comparaison RC, correspond à une mise en œuvre d’un algorithme itératif du type Levenberg-Marquardt. Selon un exemple particulier de mise en œuvre, une amplitude de la variation VAR estimée est contrainte, de sorte que la variation VAR estimée est comprise entre une valeur minimum prédéterminée et une valeur maximum prédéterminée. Selon ces dispositions, la valeur estimée des paramètres VEP, converge progressivement, au fur et à mesure des itérations successives, vers une valeur estimée des paramètres VEP de plus en plus proche d’une valeur limite. Ainsi la réponse impulsionnelle modélisée selon la valeur des paramètres VEP estimée à l’issue de l’étape de répétition 107bis, 107 est représentative de la réponse impulsionnelle observée, de sorte qu’une mesure de l’aberration de l’instrument optique peut être déterminée, au cours d’une étape de détermination de la mesure de l’aberration 108, sur la base valeur estimée des paramètres ; optionnellement cette mesure de l’aberration peut être utilisée pour compenser l’aberration, par un traitement utilisant la réponse impulsionnelle ainsi modélisée sur la base de la valeur estimée des paramètres. Selon un exemple de mise en œuvre, le traitement pour la compensation de l’aberration est un traitement d’amélioration de l’image acquise par l’instrument optique, par exemple par application d’un filtre de compensation de l’image calculé à partir de la mesure de l’aberration, i.e. de la valeur estimée des paramètres du modèle paramétrique, ou par application d’un traitement de restauration de l’image, fonction de la mesure de l’aberration, i.e. de la valeur estimée des paramètres du modèle paramétrique. Selon un autre exemple de mise en œuvre, le traitement pour la compensation de l’aberration est une commande de l’instrument optique pour actionner une optique active de l’instrument optique, la commande étant basée sur la valeur estimée des paramètres du modèle paramétrique. Optionnellement le traitement pour la compensation de l’aberration peut combiner une commande de l’instrument optique pour actionner une optique active de l’instrument optique, et un traitement d’amélioration de l’image. Il est à noter que plus le nombre de matrices de filtres spectraux est important plus l’estimation de la valeur des paramètres communs du modèle paramétrique est robuste, et meilleure est la qualité de la compensation de l’aberration, obtenue grâce au modèle. Comme cela a été indiqué, l’ensemble de paramètres comprend des paramètres globaux caractéristiques d’un état global de la combinaison optique, et des paramètres spécifiques à certains éléments de la chaîne d’acquisition de l’instrument optique, notamment spécifiques à chaque bande spectrale. Selon un exemple particulier, l’ensemble des paramètres du modèle paramétrique comprend un sous-ensemble de paramètres spécifiques à chaque nouvelle acquisition d’une source ponctuelle connue, le modèle paramétrique de réponse impulsionnelle étant représentatif des différentes réponses impulsionnelles observées dans la bande spectrale de chaque matrice de filtre spectral, le cas échéant pour chaque source ponctuelle et/ou pour chaque acquisition d’une même source ponctuelle. Les étapes suivantes sont mises en œuvre pour chaque matrice de filtre spectral, le cas échéant pour chaque source ponctuelle, et/ou pour chaque acquisition d’une même source ponctuelle : - étape de détermination 102 de la valeur initiale des paramètres VIP, et étape d’estimation 102bis de la valeur estimée des paramètres VEP, - étapes d’échantillonnage 103 , 103bis de la réponse impulsionnelle calculée sur la base de la valeur initiale des paramètres VIP, respectivement de la nouvelle réponse impulsionnelle calculée sur la base de la valeur estimée des paramètres VEP, - étapes de calcul 104, 104bis de la valeur du coût, respectivement de la nouvelle valeur du coût, et dans lequel la valeur du coût global, respectivement la nouvelle valeur du coût global, est calculée 105, 105bis pour l’ensemble des bandes spectrales et pour l’ensemble des sources ponctuelles. Selon un exemple particulier de mise en œuvre, parmi la pluralité d’acquisitions successives d’une source ponctuelle, au moins un paramètre du modèle paramétrique prend une valeur au cours d’une acquisition et une autre valeur différente au cours d’au moins une autre acquisition, la valeur et l’autre valeur dudit au moins un paramètre étant inconnues lors de l’étape d’acquisition. La variation appliquée entre la valeur et l’autre valeur n’est pas nécessairement connue. Si la variation du paramètre est connue, l’estimation de la valeur des paramètres du modèle paramétrique sera plus robuste. Selon un aspect, l’invention concerne un programme d’ordinateur comprenant un ensemble d’instructions configuré pour mettre en œuvre le procédé 100 selon l’un des modes de mise en œuvre décrits ci-avant, lorsqu’il est exécuté par un processeur d’une unité de traitement des images et de contrôle d’un instrument optique comprenant une matrice de détecteurs multispectraux, la matrice de détecteurs multispectraux comprenant une matrice de détecteurs, et une pluralité de matrices de filtres spectraux, chaque matrice de filtre spectral étant associée à une bande spectrale. Selon un autre aspect, l’invention concerne un support lisible par un ordinateur comprenant le programme d’ordinateur décrit ci-avant. Procédé (100) de mesure d’une aberration dans une image acquise par un instrument optique comprenant une matrice de détecteurs multispectraux, la matrice de détecteurs multispectraux comprenant une matrice de détecteurs, et une pluralité de matrices de filtres spectraux, chaque matrice de filtre spectral étant associée à une bande spectrale, le procédé (100) comprenant les étapes suivantes : - acquisition (101) par la matrice de détecteurs de l’instrument optique d’une image d’une source ponctuelle connue, l’image comprenant, pour chaque matrice de filtre spectral, une partie d’une tâche image de la source ponctuelle acquise dans la bande spectrale de ladite matrice de filtre spectral ; - pour chaque matrice de filtre spectral, détermination (102) d’une valeur initiale (VIP) de chaque paramètre d’un ensemble de paramètres d’un modèle paramétrique d’une réponse impulsionnelle modélisée de l’instrument optique, l’ensemble de paramètres comprenant un sous-ensemble de paramètres commun à la pluralité de matrices de filtres spectraux, et un sous-ensemble de paramètres spécifiques à chaque matrice de filtre spectral; - pour chaque matrice de filtre spectral, échantillonnage (103) d’une réponse impulsionnelle calculée sur la base de la valeur initiale des paramètres (VIP) appliquée au modèle paramétrique, l’échantillonnage étant réalisé selon un masque d’échantillonnage associé à la matrice de filtre spectral, - pour chaque matrice de filtre spectral, calcul (104) d’une valeur d’un coût dans la bande spectrale, la valeur du coût étant fonction d’une distance entre la réponse impulsionnelle calculée et la réponse impulsionnelle observée pour chacun des points du masque d’échantillonnage, - calcul (105) d’une valeur d’un coût global (CG) pour l’ensemble des bandes spectrales en fonction de la valeur du coût dans chaque bande spectrale, - pour chaque matrice de filtre spectral, estimation (102bis) d’une valeur estimée (VEP) de chaque paramètre de l’ensemble de paramètres, la valeur estimée (VEP) étant fonction d’une variation (VAR) estimée à partir de la valeur du coût global (CG), la valeur estimée (VEP) étant égale à la variation estimée (VAR) ajoutée à la valeur initiale des paramètres (VIP); - pour chaque matrice de filtre spectral, échantillonnage (103bis) d’une nouvelle réponse impulsionnelle calculée sur la base de la valeur estimée (VEP) des paramètres appliquée au modèle paramétrique, l’échantillonnage étant réalisé selon le masque d’échantillonnage associé à la matrice de filtre spectral, - pour chaque matrice de filtre spectral, calcul (104bis) d’une nouvelle valeur du coût dans ladite bande spectrale, la nouvelle valeur du coût étant fonction de la distance entre la nouvelle réponse impulsionnelle calculée et la réponse impulsionnelle observée pour chacun des points du masque d’échantillonnage, - calcul (105bis) d’une nouvelle valeur du coût global (NCG) pour l’ensemble des bandes spectrales en fonction de la nouvelle valeur du coût dans chaque bande spectrale. - comparaison (106) entre la nouvelle valeur du coût global (NCG) et la valeur du coût global (CG) pour obtenir un résultat de la comparaison (RC); - déterminer (107bis) que la valeur initiale des paramètres (VIP) est égale la valeur estimée des paramètres (VEP) et que la valeur du coût global (CG) est égale à la nouvelle valeur du coût global (NCG), et pour chaque matrice de filtre spectral, répéter (107) les étapes comprises entre l’étape d’estimation (102bis) d’une valeur estimée des paramètres (VEP), jusqu’à l’étape de comparaison (106) ) entre la nouvelle valeur du coût global (NCG) et la valeur du coût global (CG) jusqu’à ce qu’un critère d’arrêt (CRA) soit satisfait par le résultat de la comparaison (RC), de sorte qu’à chaque itération de l’étape de répétition (107), au cours de l’étape d’estimation (102bis) la variation (VAR) est estimée à partir d’une valeur initiale des paramètres (VIP) égale à la valeur estimée des paramètres (VEP) estimée au cours de l’itération précédente, et à partir d’une valeur du coût global (CG) égale à une nouvelle du coût global calculée au cours de l’itération précédente; - déterminer (108) une mesure de l’aberration sur la base de la valeur estimée des paramètres (VEP) de la réponse impulsionnelle. . Procédé selon la revendication 1, dans lequel l’étape d’acquisition (101) comprend l’acquisition de plusieurs sources ponctuelles connues et/ou plusieurs acquisitions d’une même source ponctuelle connue, l’ensemble de paramètres du modèle paramétrique comprenant un sous-ensemble de paramètres spécifiques à chaque acquisition, le modèle paramétrique de réponse impulsionnelle étant associé aux différentes réponses impulsionnelles observées dans la bande spectrale de chaque matrice de filtre spectral, et pour chaque source ponctuelle et/ou pour chaque acquisition d’une même source ponctuelle, et dans lequel les étapes suivantes sont mises en œuvre pour chaque matrice de filtre spectral, et pour chaque source ponctuelle, et/ou pour chaque acquisition d’une même source ponctuelle : - étape de détermination (102) de la valeur initiale des paramètres (VIP), et étape d’estimation (102bis) de la valeur estimée des paramètres (VEP), - étapes d’échantillonnage (103 , 103bis) de la réponse impulsionnelle calculée sur la base de la valeur initiale des paramètres (VIP), respectivement de la nouvelle réponse impulsionnelle calculée sur la base de la valeur estimée des paramètres (VEP), - étapes de calcul (104, 104bis) de la valeur du coût, respectivement de la nouvelle valeur du coût, et dans lequel la valeur du coût global, respectivement la nouvelle valeur du coût global, est calculée (105, 105bis) pour l’ensemble des bandes spectrales et pour l’ensemble des sources ponctuelles. Procédé selon la revendication 1 ou 2, dans lequel la pluralité de matrices de filtres spectraux comprend au moins trois matrices de filtres spectraux, configurées pour filtrer chacune dans une bande spectrale complémentaire des autres matrices. Procédé selon l’une des revendications 1 à 3, dans lequel, à l’étape d’estimation 102bis de la valeur estimée (VEP) de chaque paramètre, la variation (VAR) est estimée à partir d’une matrice jacobienne du modèle paramétrique de la réponse impulsionnelle, de sorte que la répétition (107), des étapes comprises entre l’étape d’estimation (102bis) d’une valeur estimée des paramètres (VEP), et l’étape de comparaison (106) entre la nouvelle valeur du coût global (NCG) et la valeur du coût global (CG), jusqu’à ce qu’un critère d’arrêt (CRA) soit satisfait par le résultat de la comparaison (RC), correspond à une mise en œuvre d’un algorithme itératif du type Levenberg-Marquardt. Procédé selon l’une des revendications 1 à 4, dans lequel la valeur du coût global, respectivement la nouvelle valeur du coût global, est égale à une somme de la valeur des coûts, respectivement de la nouvelle valeur des coûts, sur l’ensemble des bandes spectrales. Procédé selon l’une des revendications 1 à 5, dans lequel le résultat de la comparaison est égal à une différence entre la valeur du nouveau coût global et la valeur du coût global, le critère d’arrêt prédéterminé étant satisfait lorsque le résultat de la comparaison est inférieur à un seuil. Procédé selon l’une des revendications 1 à 6, dans lequel l’ensemble des paramètres comprend au moins l’un parmi les paramètres d’une réponse impulsionnelle représentative d’une partie optique de l’instrument optique, les paramètres d’une réponse impulsionnelle représentative de l’intégration des charges sur au moins un détecteur élémentaire de la matrice de détecteurs, les paramètres d’une réponse impulsionnelle représentative d’un filé, i.e. d’un mouvement linéaire, de la matrice de détecteurs pendant un temps de pose, les paramètres d’une réponse impulsionnelle représentative d’un bougé, i.e. d’un mouvement non linéaire, de la matrice de détecteurs pendant le temps de pose, les paramètres représentatifs des coordonnées d’une image géométrique de la source ponctuelle, ladite image géométrique étant ponctuelle.. Procédé selon l’une des revendications 1 à 7, dans lequel la mesure de l’aberration est utilisée pour compenser l’aberration par un traitement d’amélioration de l’image acquise par l’instrument optique, le traitement d’amélioration étant fonction de la mesure de l’aberration, le traitement comprenant par exemple une application d’un filtre de compensation calculé à partir de la mesure de l’aberration, ou le traitement comprenant par exemple une restauration de l’image fonction de la mesure de l’aberration. Procédé selon l’une des revendications 1 à 8, dans lequel la mesure de l’aberration est utilisée pour compenser l’aberration par une commande de l’instrument optique pour actionner une optique active de l’instrument optique, la commande étant basée sur la mesure de l’aberration. Programme d’ordinateur comprenant un ensemble d’instructions configuré pour mettre en œuvre le procédé 100 selon l’une des revendications 1 à 9 lorsqu’il est exécuté par un processeur d’une unité de traitement des images et de contrôle, d’un instrument optique comprenant une matrice de détecteurs multispectraux, la matrice de détecteurs multispectraux comprenant une matrice de détecteurs, et une pluralité de matrices de filtres spectraux, chaque matrice de filtre spectral étant associée à une bande spectrale. Support lisible par un ordinateur comprenant un programme d’ordinateur selon la revendication 10.