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.