FR 3153958 A1 20250411 FR 2310710 A 20231006 Méthode de correction d'un défaut de vignettage dans une séquence d'images à partir de tables d'étalonnage Un aspect de l’invention concerne une méthode mise en œuvre par ordinateur de correction d'un défaut de vignettage dans une séquence d'images acquise par une caméra pour une distance focale F', la séquence d'image comprenant une pluralité de sous-séquences d'images, la méthode étant réalisée à partir d'un ensemble de tables d'étalonnage correspondant au défaut de vignettage et associées chacune à une distance focale propre, au moins une distance focale propre associée à une table d'étalonnage étant différente de F', la méthode comprenant une pluralité d'itérations. Figure à publier avec l’abrégé : Figure 2 Méthode de correction d'un défaut de vignettage dans une séquence d'images à partir de tables d'étalonnage DOMAINE TECHNIQUE DE L’INVENTION Le domaine technique de l’invention est celui de la correction de vignettage dans des séquences d'images. La présente invention concerne une méthode de correction d'un défaut de vignettage dans une séquence d'images. ARRIERE-PLAN TECHNOLOGIQUE DE L’INVENTION Les caméras infrarouges ou détectant la lumière visible comprises dans des viseurs ou aéronefs par exemple comprennent un système de lentilles de verre et/ou des miroirs, la position des lentilles et des miroirs étant variable pour ajuster la distance focale du système et grossir les objets observés. En particulier, le mouvement des lentilles permet de faire varier continument le grossissement et notamment d'ajuster dynamiquement le plan de focalisation. Cependant, le mouvement des miroirs ou des lentilles compris dans les caméras peut contribuer à des défauts de vignettage variables dans le temps, ces défauts affectant la qualité des images produites en temps réel. En particulier, les défauts de vignettage comprennent une composante majoritairement radiale et correspondent à un assombrissement des périphéries d'une image, la lumière extérieure étant atténuée sur les bords de l'image. Le vignettage correspond à une perte progressive d'éclairement en périphérie d'une image, en allant du centre de cette image vers sa périphérie et réduit donc la qualité des images acquises. La FIG. 1 est un exemple de système optique 100 comportant un zoom optique 101 infrarouge refroidi, un détecteur 102 et un écran froid 103. Dans ce cas, l'écran froid 103 est la pupille de sortie du système optique. Le zoom optique 101 comprend notamment une lentille frontale 1011 qui est la pupille d'entrée du système optique et pour laquelle on cherche à obtenir une conjugaison avec l'écran froid 103. La conjugaison entre les deux pupilles n'est pas parfaite et, pour des raisons d’encombrement du zoom optique 101, il est parfois nécessaire de limiter le diamètre utile, autrement dit l'ouverture, de la lentille frontale à la zone de passage Z des rayons sur l’axe optique. Par conséquent, plus l'objet s'éloigne de l'axe et moins l'image reçoit de rayons lumineux: une partie des rayons V arrivant sur lentille frontale et non comprise dans la zone Z est coupée et n'atteint pas le détecteur 102, ce qui est à l'origine du vignettage. Il existe des méthodes de correction de vignettage de séquences d'images acquises, à partir de tables d'étalonnage correspondant à la forme du vignettage, les tables d'étalonnage étant également appelées tables de calibration. Les tables de calibration sont chacune associées à un défaut de vignettage pour des paramètres tels qu'une distance focale donnée, un temps d'intégration donné et une température donnée. Cependant, les tables de calibration ne comprennent pas de paramètres relatifs à la scène filmée par la caméra et ne permettent pas de combler l'absence de rayons à l'origine du vignettage. Il est possible de corriger le vignettage grâce à une corrélation entre les hautes fréquences d'une image acquise en temps réel et les hautes fréquences des tables de calibration. L’estimation des coefficients à partir des hautes fréquences est par exemple décrite par exemple dans les brevets EP2127359B1 et FR3118558. Cependant, l’utilisation des hautes fréquences des images pour déterminer et corriger le vignettage peut introduire des défauts supérieurs aux défauts radiaux originaux. Il est possible d'estimer la composante radiale des défauts de vignettage par des traitements statistiques sur des cercles concentriques puis en l’appliquant de manière radiale sur l’ensemble de l’image comme décrit dans le brevet US9730958B2. Cependant, le résiduel obtenu après une correction par une estimation radiale ou elliptique du défaut reste conséquent. Il existe ainsi un besoin de corriger le vignettage d'images acquises en temps réel de manière optimale. L’invention offre une solution aux problèmes évoqués précédemment, en considérant le défaut de vignettage comme un bruit spatial additif radial et à basse/moyenne fréquence. En particulier, l'invention met en avant ce défaut dans des séquences d'images à corriger en écartant les hautes fréquences avec un filtre radial. Des tables d'étalonnage caractérisant le défaut de vignettage et réalisées en usine ou par modélisations obtenues pour des conditions proches des images à corriger sont projetées sur les images filtrées à corriger, pour quantifier la présence des images de défaut dans les images à traiter. Enfin, les tables d'étalonnage sont soustraites aux images à corriger selon leurs présences dans les images à corriger. Un aspect de l’invention concerne une méthode, mise en œuvre par ordinateur, de correction d'un défaut de vignettage dans une séquence d'images acquise par une caméra pour une distance focale F', la séquence d'image étant ordonnée temporellement et comprenant une pluralité de sous-séquences d'images ordonnées temporellement, la méthode comportant un ensemble d'itérations, chaque itération utilisant une sous-séquence d'images de la pluralité de sous-séquences d'images et une table d'étalonnage, chaque itération comprenant: Réduire des dimensions d'une première image de la séquence d'images pour obtenir une première image réduite, la première image étant corrigée du défaut de vignettage à partir d'une correction du défaut de vignettage déterminée à l'itération précédente ou à partir d'une correction initiale; Réduire des dimensions de la table d'étalonnage pour obtenir une table d'étalonnage réduite; Appliquer un filtre radial passe-haut sur la première image réduite et sur la table d'étalonnage réduite pour obtenir une première image réduite filtrée et une table d'étalonnage réduite filtrée; Déterminer un coefficient de corrélation entre la première image réduite filtrée et la table d'étalonnage réduite filtrée; Déterminer une correction du défaut de vignettage à partir de la table d'étalonnage, du coefficient de corrélation associé à la table d'étalonnage et de la correction du défaut de vignettage obtenue à l'itération précédente; Corriger chaque autre image de la sous-séquence d'images à partir de la correction du défaut déterminée pour obtenir chaque autre image corrigée du défaut de vignettage. Grâce à l’invention, il est possible d'exprimer et quantifier le défaut optique de vignettage compris dans une image en une combinaison linéaire de tables d'étalonnage, la combinaison linéaire étant mise à jour à chaque itération, et de corriger ainsi au moins une autre image de la sous-séquence d'images en soustrayant le défaut de vignettage. En particulier, le coefficient de corrélation déterminé à chaque étape dépend d'une image corrigée précédemment et donc des coefficients de corrélations déterminés précédemment, ce qui permet de faire converger la méthode. En particulier, lorsque le défaut optique à corriger comprend une forte composante radiale à basse ou moyenne fréquence, l'étape de réduction des dimensions de la première image permet d'adapter l'étape de filtrage aux dimensions caractéristiques du défaut de vignettage et de supprimer les motifs de haute fréquence qui ne correspondent pas au défaut de vignettage à corriger. En outre, le filtre radial met généralement en évidence tous les motifs de haute fréquence et tous les détails de l'image compris dans des cercles concentriques, en éliminant des motifs longitudinaux ou verticaux. Cependant, dans le cadre de l'invention, l'application du filtre radial suite à la réduction de dimensions de la première image et de la table d'étalonnage met en évidence les moyennes fréquences de la première image réduite et de la table d'étalonnage réduite. Ainsi, le filtre radial permet de supprimer tous les motifs non orientés dans la direction du rayon et ne correspondant pas au défaut de vignettage. De plus, la corrélation de chaque table d'étalonnage avec chaque première image d'une sous-séquence d'images permet de quantifier la présence de la table d'étalonnage dans la première image afin d'y soustraire la table d'étalonnage proportionnellement à sa présence. Enfin, la correction est mise à jour en prenant en compte, pour chaque sous-séquence, la correction de la sous-séquence précédente afin d'obtenir une précision de correction élevée et permettre une convergence de la méthode. Outre les caractéristiques qui viennent d’être évoquées dans le paragraphe précédent, la méthode selon un aspect de l’invention peut présenter une ou plusieurs caractéristiques complémentaires parmi les suivantes, considérées individuellement ou selon toutes les combinaisons techniquement possibles. Selon un mode de réalisation, la première image d'une sous-séquence d'images utilisée à une itération donnée correspond à la dernière image d'une sous-séquence d'images utilisée à l'itération précédente. Selon un mode de réalisation, la méthode est réalisée pour un ensemble de K tables d'étalonnage, K étant un entier naturel non nul, et pour chaque itération donnée de l'ensemble d'itérations, la table d'étalonnage utilisée est distincte de la table d'étalonnage utilisée lors de l'itération précédente. Avantageusement, lorsqu'un nombre K de tables d'étalonnage est utilisé la méthode selon l'invention peut être mise en œuvre sur une carte FPGA et la combinaison linéaire des tables permettant de déterminer la correction est réalisée en 2*K premières images des sous-séquences traitées en parallèles. Selon un mode de réalisation, l'ensemble des K tables d'étalonnage est choisi selon une règle prédéfinie parmi M tables d'étalonnage, M étant un entier naturel non nul avec K inférieur à M, la règle prédéfinie étant telle que, pour chaque table d'étalonnage parmi les K tables d'étalonnage: la distance focale associée à la table d'étalonnage correspond à l'une des K distances focales les plus proche de F'; ou L'écart entre la distance focale associée à la table d'étalonnage et F' est inférieur à une distance focale moyenne entre au moins deux tables d'étalonnage successives parmi les K tables d'étalonnage choisies. Selon un mode de réalisation, K est compris entre 2 et 10. En particulier, le nombre de tables d'étalonnage K peut être adapté en fonction du compromis entre la réactivité de la méthode et la précision souhaitée. Plus le nombre de tables d'étalonnage est élevé, plus la correction est précise, mais plus le temps de convergence de la méthode est élevé. Ainsi, K compris entre 2 et 10 est un compromis entre la précision de la correction et la convergence. Selon un mode de réalisation, le défaut de vignettage est un bruit additif spatial, ayant une composante radiale et de basse fréquence. On entend par "basse fréquence" une fréquence de période comprise entre 10 et 100 pixels. Selon un mode de réalisation, la réduction de dimensions de la première image et de la table d'étalonnage est réalisée par un filtre moyenneur. Selon un mode de réalisation, le filtre radial est un filtre gradient ou un filtre Laplacien orienté dans la direction du centre de chaque image de chaque sous-séquence d'images. Avantageusement, le filtre radial permet de supprimer des informations parasites qui ne sont pas dans la direction du rayon du défaut de vignettage radial pour gagner en robustesse Selon un mode de réalisation, le coefficient de corrélation entre la première image réduite filtrée et la table d'étalonnage réduite filtrée est déterminé à partir du rapport: d'un produit scalaire canonique entre la première image réduite filtrée et de la table d'étalonnage réduite filtrée; et d'une norme de la table d'étalonnage réduite filtrée. Avantageusement, le produit scalaire entre les deux images permet de projeter la table d'étalonnage réduite filtrée sur la première image réduite filtrée afin d'obtenir la quantité de table dans la scène. Un autre aspect de l'invention concerne un produit programme d'ordinateur comprenant des instructions qui, lorsque le programme est exécuté par un ordinateur, conduisent celui-ci à mettre en œuvre la méthode selon l'un quelconque des modes de réalisations précédents. Un autre aspect de l'invention concerne un support d'enregistrement lisible par ordinateur comprenant des instructions qui, lorsqu'elles sont exécutées par un ordinateur, conduisent celui-ci à mettre en œuvre la méthode selon l'un quelconque des modes de réalisations précédents. L’invention et ses différentes applications seront mieux comprises à la lecture de la description qui suit et à l’examen des figures qui l’accompagnent. BREVE DESCRIPTION DES FIGURES Les figures sont présentées à titre indicatif et nullement limitatif de l’invention. La FIG. 1 est une représentation schématique d'un zoom optique infrarouge refroidi dans lequel une partie des rayons optiques incidents coupée est à l'origine du vignettage. La FIG. 2 est un schéma synoptique représentant des étapes d'une méthode de correction d'un défaut de vignettage selon l'invention. La FIG. 3 est une représentation de l'orientation d'un filtre gradient radial utilisé dans la méthode de correction selon l’invention. La FIG. 4 est un schéma synoptique représentant des étapes d'un premier mode de réalisation d'une méthode de détermination d'une table d'étalonnage. La FIG. 5 est un schéma synoptique représentant des étapes d'un deuxième mode de réalisation de la méthode de détermination d'une table d'étalonnage. La FIG. 6 est une représentation d'un dispositif configuré pour mettre en œuvre la méthode selon l'invention. DESCRIPTION DETAILLEE L'invention concerne une méthode de correction d'un défaut de vignettage dans une séquence d'images. La séquence d’images comporte une pluralité d’images et est ordonnée temporellement, c’est-à-dire que la pluralité d’images est rangée par instant d’acquisition croissant dans la séquence d’images. Ainsi, si une première image est acquise à un instant t, une deuxième image est acquise à un instant t+1, et une troisième image est acquise à un instant t+2, la première image, la deuxième image et la troisième image seront rangées dans cet ordre dans une séquence d’images les comprenant. La séquence d'images est acquise par une caméra infrarouge ou une caméra visible. La séquence d'images est obtenue à partir d'une caméra située par exemple dans un aéronef ou dans un viseur d'observation situé dans un véhicule terrestre. Ainsi, la caméra permet d'acquérir des séquences d'images qui peuvent être affichées en temps réel à un utilisateur de l'aéronef ou du véhicule terrestre par exemple. La séquence d'images comprend une pluralité de sous-séquences d'images. La séquence d'images comprend de préférence au moins deux sous-séquences. Une sous séquence d'images comprend au moins deux images par exemple et comprend de préférence trois images. La sous-séquence d'images peut notamment comprend quatre, cinq voire dix images. En particulier, la séquence d'images étant acquise pour des instants successifs, les sous-séquences d'images de la séquence d'image se succèdent temporellement. En particulier, les images de chaque sous-séquence d'images se succèdent temporellement. On entend par "caméra infrarouge" une caméra configurée pour mesurer et enregistrer des ondes de chaleur et rayonnements infrarouges émis par un corps ou un objet, afin de générer des séquences d'images dites infrarouges. On entend par "caméra visible" une caméra configurée pour mesurer et enregistrer de la lumière visible, dans le spectre du visible qui s'étend de 300 à 750 nanomètres de longueur d’onde afin de générer des séquences d'images dans le spectre du visible. En particulier, la caméra d'acquisition comprend un objectif de distance focale variable. Dans la suite, on confondra "distance focale" et "focale". On note F' la focale de l'objectif de la caméra lors de l'acquisition de la séquence d'images. En particulier, les séquences d'images acquises par la caméra peuvent comprendre un défaut de vignettage qui réduit la qualité des images et qui est corrigé selon la méthode selon invention. Selon l'invention, le défaut de vignettage est modélisé selon une pluralité d'hypothèses décrites dans la suite. Selon une première hypothèse, le défaut de vignettage est une image constante ou "offset" de l'anglais, modélisable comme un bruit additif spatial et fixe dans le temps. On entend par "bruit fixe dans le temps" un bruit dont le temps caractéristique d'évolution est grand devant le temps de convergence de la méthode selon l'invention, le temps caractéristique d'évolution étant de quelques secondes alors que le temps de convergence de la méthode selon l'invention correspond à quelques images. En particulier, lorsque la caméra acquiert des séquences d'images à une cadence comprise entre 25Hz et 50Hz, on entend par "temps de convergence de quelques images" un temps d'environ quelques centaines de millisecondes. En particulier, le temps caractéristique d'évolution du défaut de vignettage est environ dix fois plus grand au temps de convergence de la méthode selon l'invention. Selon une deuxième hypothèse, le défaut de vignettage comprend un signal de moyenne à basse fréquence spatiale. On entend par moyenne à basse fréquence une fréquence dont la période est comprise entre 10 et 100 pixels. Selon une troisième hypothèse, le défaut de vignettage a une composante spatiale radiale dominante par rapport à une composante horizontale ou une composante verticale. Un défaut de composante spatiale radiale dominante par rapport à une composante horizontale ou une composante verticale est un défaut dont la moyenne des gradients est la plus grande dans la direction des rayons passant par le centre de l’image. Cela se traduit visuellement par des motifs concentriques sur l'image comprenant le défaut de vignettage. En pratique la concentricité n’est pas exacte et la position du centre difficilement caractérisable, d’où l’élimination des corrections radiales basées sur un modèle strict. On entend par "image à composante radiale dominante" une image dont la moyenne des valeurs absolue de l’image filtrée par un filtre radial est au moins cinq fois supérieure à la moyenne des valeurs absolues de ladite image filtrée par un filtre orienté orthogonalement au rayons qui passent par le centre de l’image. La méthode de correction selon le premier aspect de l'invention est réalisée à partir d'un ensemble de K tables d'étalonnage sélectionnées parmi un ensemble de M tables d'étalonnage, avec K et M deux entiers naturels non nuls et K inférieur à M. Par exemple, M est compris entre 1 et 100. Par exemple K est compris entre 1 et 10. Par exemple, K peut être égale à 1. En particulier, K est de préférence égal à 2, 3 ou 4. Chaque table d'étalonnage parmi l'ensemble des M tables d'étalonnage peut par exemple correspondre à une image acquise ou modélisée pour une distance focale F' m de l'objectif de la caméra, avec m compris entre 1 et M. Les tables d'étalonnage comprennent des coefficients caractérisant le défaut de vignettage optique en fonction de paramètres tels que la distance focale de l'objectif de la caméra, la température de la caméra et le temps d'intégration par exemple. En particulier, les tables d'étalonnage ont la même taille que la matrice de pixels du détecteur de la caméra et donc la même taille que les séquences d'images acquises par la caméra. Les tables d'étalonnage sont par exemple réalisées en usine. En particulier, la méthode selon l'invention peut être réalisée à partir de n'importe quel ensemble de tables d'étalonnage correspondant chacune à un défaut de vignettage vérifiant les hypothèses citées précédemment et obtenues pour des focales différentes. Des exemples de tables d'étalonnage sont détaillés plus loin dans la description. En particulier, la méthode de correction selon l'invention est réalisée pour une séquence d'images acquise à une distance focale F' et les tables d'étalonnage de l'ensemble de K tables d'étalonnage permettant de réaliser la méthode sont sélectionnées parmi les M tables d'étalonnage, les K tables d'étalonnage étant sélectionnées selon une condition prédéterminée parmi les M tables d'étalonnage. Selon un mode de réalisation, la condition prédéterminée est que la focale F' k de chaque table d'étalonnage de l'ensemble des K tables d'étalonnage, avec k un entier naturel compris entre 1 et K, est l'une des K focales les plus proches de la focale F' parmi les M focales de l'ensemble des M tables d'étalonnage. En d'autres termes, les tables d'étalonnage associées aux K focales les plus proches de la focale F' parmi les M focales sont sélectionnées. En particulier, les K distances focales les plus proches de F' parmi les M distances focales de l'ensemble des M images de défaut sont déterminées par calcul d'une norme de la différence F' m -F', la norme étant par exemple une norme euclidienne, et par sélection des K normes minimales. Selon un autre mode de réalisation, la condition prédéterminée est que l'écart entre chaque focale F' k et la focale F' est inférieur à un seuil prédéfini correspondant à un rapport d'une distance focale moyenne entre au moins deux tables d'étalonnage successives parmi les K tables d'étalonnage choisies et de la focale F'. Le seuil prédéfini est par exemple inférieur ou égal à 20%. Ainsi, on ne sélectionne pas un nombre prédéterminé de tables d'étalonnage mais on sélectionne toutes les tables d'étalonnage correspondant à une focale F' k respectant la condition prédéterminée que l’écart entre la focale F' k et la focale F' est inférieur à un seuil prédéfini. En particulier, au moins une focale F' k associée à une table d'étalonnage de l'ensemble de K tables d'étalonnage est différente de la focale F'. Ainsi, les K tables d'étalonnage sélectionnées varient selon la focale F' d'acquisition de la séquence d'images. La FIG. 2 montre une représentation schématique de la méthode 200 de correction d'un défaut de vignettage selon l'invention dans une séquence d'images telle que décrite précédemment. La séquence d'images est découpée en une pluralité de sous-séquences d'images comportant chacune une pluralité d’images successives dans la séquence d'images. La pluralité de sous-séquences d'images est donc ordonnée temporellement dans la séquence d’images. Par exemple, la séquence d'images comprend P sous-séquences d'images, P étant un entier naturel non nul, compris entre 1 et 100. Chaque sous-séquence d'images comprend de préférence le même nombre d'images, par exemple égal à deux ou de préférence à trois. Par exemple, chaque séquence d'images est découpée de manière qu’une sous-séquence d'images donnée comprenne la dernière image, de la sous-séquence d'images précédant la sous-séquence d'images donnée dans la séquence d'images et la première image, de la sous-séquence d'images suivant la sous-séquence d'images donnée dans la séquence d'images. Par exemple, en considérant une séquence d’images comprenant sept images SEQ = {IM1, IM2, IM3, IM4, IM5, IM6, IM7}, la séquence d’images est découpée en une première sous-séquence d’images SEQ 1 = {IM1, IM2, IM3}, une deuxième sous-séquence d’images SEQ 2 = {IM3, IM4, IM5}, une troisième sous-séquence d’images SEQ 3 = {IM5, IM6, IM7}. Ainsi, la méthode de correction comprend des étapes 201, 202, 203, 204, 205 et 206 réalisées de manière itérative pour chaque sous-séquence SEQ i d'images d'une séquence d'images acquise par une caméra infrarouge ou une caméra visible telle que décrite précédemment. En particulier, la sous-séquence d'images est acquise pour une focale F' de la caméra et comporte un défaut de vignettage sur chacune de ses images. En particulier, la méthode comprend un ensemble d'itérations, l'ensemble comprenant un nombre d'itérations de préférence égal au nombre de sous-séquences SEQ i . A chaque itération une unique séquence est traitée. Ainsi, considérons une itération donnée i, avec i compris entre 1 et P, l'itération donnée utilisant une i-ème sous-séquence SEQ i d'images de la séquence d'images à corriger et sur laquelle la méthode 200 est appliquée. A chaque itération i, une table d'étalonnage, notée T k , parmi l'ensemble de K tables d'étalonnage est utilisée. Dans un mode de réalisation, K est égal à 1 et l'ensemble de K tables d'étalonnage comprend une seule table d'étalonnage, notée T 1 . Dans un autre mode de réalisation, K est supérieur ou égal à 2 et chaque itération parmi l'ensemble d'itérations utilise une table d'étalonnage distincte. De préférence, K est inférieur au nombre total d'itérations de l'ensemble d'itérations et chaque table d'étalonnage est utilisée au moins une fois. Ainsi, considérons l'ensemble de tables d'étalonnage avec K égal à 2 par exemple. L'ensemble de tables d'étalonnage comprend ainsi la table d'étalonnage T 1 et la table d'étalonnage T 2 . Dans cet exemple, si la table d'étalonnage T 1 est utilisée à l'itération i, alors la table d'étalonnage T 2 est utilisée à l'itération i+1, et la table T 1 est à nouveau utilisée à l'itération i+2. Ainsi, pour l'itération i donnée, la méthode comprend une étape 201 de réduction des dimensions d'une première image IM 1 _ i de la séquence SEQ i d'un facteur de réduction N. En particulier, la première image IM 1 _i est corrigée du défaut de vignettage grâce à une correction du défaut de vignettage, notée Corr( i -1 ) k , déterminée à l'itération précédente, c’est-à-dire l'itération i-1. Si l'itération i donnée est telle que i est égal à 1, alors il n y a pas d'itération précédente alors et la première image IM 1 _i est corrigée du défaut de vignettage à partir d'une correction initiale dont un exemple est donné plus loin dans la description. Le facteur de réduction N est un entier naturel non nul compris entre 1 et 64. On entend par "dimension" d'une image sa taille exprimée en nombre de pixels. En particulier, chaque image de la séquence d'images à corriger comprend L lignes de pixels et C colonnes de pixels, sa dimension étant égale à L*C pixels, L et C sont deux entiers naturels non nuls. Ainsi, il en est de même pour chaque image de chaque sous-séquence d'images. La première image a donc une dimension de L*C pixels avant l'étape 201 de réduction de dimension. Après l'étape 201 de réduction de dimension, le nombre de colonnes de la première image est divisé par N, et le nombre de lignes de la première image est divisé par N. Ainsi, la réduction de la première image IM 1 _i permet d'obtenir une première image réduite IM R 1 _i de dimension * = . La méthode comprend une étape 202 de réduction des dimensions de la table d'étalonnage T k utilisée lors de l'itération i en cours, du facteur de réduction N. La table d'étalonnage T k ayant la même dimension que la première image IM 1 _i , la réduction de la table d'étalonnage T k permet d'obtenir une table d'étalonnage réduite T R k de dimension . Les étapes de réduction de dimensions 201 et 202 sont réalisées grâce à un filtre moyenneur, qui, sur chaque image auquel il est appliqué réalise le moyennage de chaque sous-image de N*N pixels successifs sur ladite image. L'image issue de la réduction comprend ainsi N² fois moins de pixels que l'image avant réduction. La réduction de la dimension de la première image permet avantageusement d'adapter le filtrage aux dimensions caractéristiques du défaut de vignettage, et de supprimer des motifs de hautes fréquences qui, par hypothèse, n’appartiennent pas au défaut de vignettage à corriger. En effet, le défaut de vignettage provient de l'absence de rayons lumineux mais également de reflets dans les lentilles du système optiques, qui créent des motifs de moyennes et basses fréquences. Le facteur de réduction N est adapté à la résolution du capteur de la caméra et à la nature du défaut de vignettage qui peut être caractérisé en usine. Ce facteur peut éventuellement être dynamiquement adapté à la position de la focale ou à d’autres paramètres qui influent sur les dimensions caractéristiques du défaut de vignettage. L'étape 203 de la méthode 200 est une étape d'application d'un filtre radial passe-haut sur la première image réduite IM 1_i et la table d'étalonnage réduite T R k . Le filtre radial passe-haut ne filtre que les basses et moyennes fréquences. En particulier, le filtre radial passe haut ne laisse passer que les fréquences basses et moyennes de période inférieures à 3 pixels sur des images dont les dimensions sont réduites d'un facteur N. Ainsi, suite aux étapes 201 et 202, le filtre radial passe-haut est équivalent à un filtre passe bande qui qui ne laisse passer (dans la première image réduite IM 1_i et la table d'étalonnage réduite T R k ) que les fréquences dont les périodes sont proches de 3*N pixels. Par exemple, le filtre radial est un filtre Laplacien de Gaussienne (LoG), un filtre Gradient orienté dans la direction du centre de la première image réduite ou tout filtre basé sur les différences orientées dans la direction du centre de la première image réduite. Le radial filtre a une dimension (ou une taille) notée F*F', avec F et F' deux entiers naturels non nuls. Un exemple de filtre radial avec F = 3 et F' = 3 est le suivant: Pour chaque image à filtrer et pour un pixel donné de l'image à filtrer on associe un angle qui correspond à l'angle formé entre une première droite passant par le pixel central de l'image et par la position dudit pixel donné et une deuxième droite horizontale (dans le sens des lignes de l'image à filtrer) et passant par le pixel central, le sommet de l'angle étant le pixel central. En particulier, l'angle est mesuré entre la deuxième droite et la première droite, et est positif dans le sens trigonométrique. Ainsi, plusieurs pixels peuvent être associés à un même angle . En particulier, le filtre radial comprend des cercles concentriques, chaque cercle correspondant à un ensemble de pixels distants du même nombre de pixels par rapport au pixel central. En outre, le filtre radial calcule par exemple la différence entre deux pixels voisins diagonaux associés au même angle . Avantageusement, l'application du filtre radial sur la première image réduite permet de mettre en évidence la composante radiale du défaut de vignettage. Le filtre radial permet de supprimer tous les motifs non orientés dans la direction de la composante radiale (ou rayon) pour ne garder que les motifs d'intérêt et met ainsi en évidence tous les motifs de haute fréquence et tous les détails de l'image dans la direction des rayons des cercles concentriques. En effet, la composante radiale du défaut de vignettage ne présente des motifs que dans la direction des rayons des cercles concentriques. Ainsi, un filtre radial est avantageux par rapport à un filtre non orienté. La FIG. 3 est un exemple de l'orientation locale du filtre radial, qui est dans ce cas un filtre gradient. Ainsi, la FIG. 3 représente l'orientation du gradient pour chaque position de pixel de la première image. Les traits représentés représentent la direction du gradient. En référence à nouveau à la FIG. 2 , à l'issue de l'étape 203, une première image réduite filtrée IM R F1_i et une table d'étalonnage réduite filtrée T R Fk est obtenue. L'étape 204 est une étape de détermination d'un coefficient de corrélation noté c ik , à partir de la première image réduite filtrée IM R F1_i et d'de la table d'étalonnage réduite filtrée T R Fk . Le coefficient de corrélation est déterminé à partir d'un produit scalaire canonique, dont l'opérateur est noté . Dans la suite de la rédaction, on note v 1RF le vecteur ligne par exemple comprenant toutes les valeurs de pixels de la première image réduite filtrée IM RF1_i , et v T _ kRF le vecteur ligne comprenant toutes les valeurs de pixels de la table d'étalonnage réduite filtrée T RFk . Ainsi, à l'itération donnée, le coefficient de corrélation entre la première image réduite filtrée IM R F1_i et la table d'étalonnage réduite filtrée T R Fk est égal au terme suivant : En particulier, le terme correspond à la norme au carré de la table d'étalonnage réduite filtrée T RFk . Le terme permet de comparer la première image réduite filtrée et la table d'étalonnage réduite filtrée T RFk , la table d'étalonnage réduite filtrée T RFk étant projetée sur la première image réduite filtrée. Le produit scalaire permet de mesurer la présence de la table d'étalonnage T RFk dans la première image IM 1_i et donc de mesurer la présence du défaut de vignettage associé à la table d'étalonnage T RFk dans la première image IM 1_i En particulier, le coefficient de corrélation c i k peut prendre n’importe quelle valeur réelle. Le coefficient c i k peut ainsi être positif ou négatif par exemple. Il est également possible, en fonction de la normalisation des tables d'étalonnage, d’obtenir un coefficient de corrélation c i k dont l’amplitude est supérieure à 1. Le coefficient de corrélation de corrélation c i k permet ainsi de quantifier les similitudes entre la première image réduite filtrée IM R F1_i et la table d'étalonnage réduite filtrée T R Fk et de déterminer dans quelle mesure le défaut de vignettage de la table d'étalonnage considérée est présent dans la première image. L'étape 205 est une étape de détermination d'une correction Corr i k du défaut de vignettage à partir de la table d'étalonnage T R Fk , du coefficient de corrélation c ik obtenu à l'étape 204 à partir de la table d'étalonnage réduite filtrée T R Fk et de la correction Corr( i -1 ) k du défaut de vignettage déterminée à l'itération précédente, c’est-à-dire l'itération i-1 Si l'itération i donnée est telle que i est égal à 1, alors il n y a pas d'itération précédente alors et la correction Corr( i -1 ) k correspond à la correction initiale . En particulier, la correction Corr ik du défaut de vignettage est par exemple égale à une somme de la correction Corr( i -1 ) k obtenue à l'itération précédente et d'un produit de la table d'étalonnage donnée multipliée par le coefficient de corrélation c ik qui lui est associé. Ainsi, une correction Corr i k du défaut de vignettage correspond à une image de même taille que chaque image d'une sous-séquence d'images acquise par la caméra. Ainsi, . Ainsi, la correction est mise à jour à chaque itération i, pour chaque sous-séquence d'images SEQ i et dépend des coefficients de corrélations déterminés aux itérations précédant l'itération i, afin de permettre la convergence de la méthode 200. Grâce à l'aspect itératif de la méthode, celle-ci converge : le calcul du coefficient de corrélation pour une table d'étalonnage à l'itération i est calculé à partir de la première image corrigée par la correction précédemment déterminée à l'itération (i-1) pour prendre en compte les contributions des tables d'étalonnage utilisées aux itérations précédentes. L'étape 206 est une étape de correction de chaque autre image de la sous-séquence SEQ i d'images donnée à partir de la correction . Par exemple, si la séquence d'images SEQ i est notée : { IM 1_i , IM 2_i , IM 3_i } alors les étapes 201 à 205 sont réalisées sur l'image IM 1_i et l'étape de correction 206 est réalisée sur les deuxième et troisième images IM 2_i et IM 3_i afin d'obtenir des deuxième et troisième images corrigées respectives IM' 2_i et IM' 3_ . En particulier, la deuxième image corrigée IM' 2_i est par exemple obtenue par soustraction de la correction à la deuxième image IM 2_i de la sous-séquence SEQ i d'images. Ainsi, L'image IM' 2 _i correspond ainsi à l'image IM 2 _i après correction du défaut de vignettage selon la méthode 200. Ainsi, En particulier, si les sous-séquences d'images comprennent par exemple une pluralité d'images supérieure à troisième images, la correction est de préférence appliquée à toutes les images de la sous-séquences d'images à l'exception de la première image. A l'itération suivante, c’est-à-dire à l'itération i+1, la sous-séquence SEQ (i+1) comprend la dernière image de la sous-séquence SEQ i comme première image IM 1 _ ( i +1). Ainsi, dans ce mode de réalisation, SEQ ( i+1) est notée : { IM 1_ ( i +1) , IM 2_ ( i +1) , IM 3_ ( i +1) } avec IM 1_(i+1) = . En référence à la correction initiale Corr 0 utilisée dans la méthode 200 selon l'invention, la correction initiale est une matrice de même taille que le détecteur de la caméra et peut être obtenue par une méthode de calibration 1 point de la caméra. La calibration 1 point d’une caméra infrarouge refroidie est réalisée par ajustement d'une table de gain G et d'une table de décalage O (de l'anglais "offset"), l'ajustement étant réalisé à partir d’une image défocalisée ou une image d’un diaphragme (de l'anglais "shutter") de la caméra, servant à obtenir une image supposée uniforme sur le capteur. Le gain G et le décalage O sont des matrices de même taille que le détecteur de la caméra. En particulier, la réponse d'un pixel du détecteur de la caméra à un flux peut se caractériser par une fonction linéaire caractérisée par un coefficient directeur qui représente un gain et une ordonnée à l’origine qui représente un décalage. Chaque pixel possède son propre gain et son propre décalage. Par conséquent, à la réception d'un flux lumineux uniforme, les niveaux de gris générés par les pixels diffèrent, ce qui se traduit par l’apparition de non-uniformités sur l’image générée par la caméra. Afin d'éliminer les non-uniformités des pixels, il est nécessaire que les pixels aient le même gain et le même décalage. Une correction est réalisée par application de la table de correction en gain G et de la table de correction en décalage notée O. La table de gain G est multipliée par le flux reçu par chacun des pixels, et s’exprime sous la forme d’une table de facteurs appliqué à chacun des pixels qui est déterminée en usine avant la mise en service de la caméra par exemple. Après son application, tous les pixels ont la même réponse à une variation de flux donnée et leurs fonctions caractéristiques ont donc le même coefficient directeur. La table de gain peut être déterminée en réalisant la dérivée partielle d'ordre 1 de deux images prises à deux températures de scènes différentes. En particulier, la dérivée partielle est réalisée par rapport au paramètre de luminance qui dépend de la température. La table de décalages O est additionnée au flux reçu par chacun des pixels, et s’exprime donc sous la forme d’une table de coefficients ajoutés à chacun des pixels. La table de décalages O permet de corriger des défauts à la fois intrinsèques au détecteur et dus aux flux froids parasites. La table de décalages O est obtenue à partir de l’image d’une scène acquise devant un diaphragme, notée Idiaphragme, multipliée par la table de gain G. En particulier, la table de décalages O est obtenue à partir de la formule suivante: O = Idiaphragme*G-Moyenne(Idiaphragme*G). Ainsi, la correction initiale est égale à la table O. Concernant la ou les tables d'étalonnage utilisées dans la méthode de correction selon l'invention, plusieurs modes de réalisations de ces tables sont décrites dans la suite. Selon un mode de réalisation, une table d'étalonnage peut résulter de la différence entre deux images acquises à partir d'un même fond uniforme respectivement avec des distances focales différentes de la caméra. Ainsi, afin de construire M images de défaut, il est possible d'acquérir deux séries de M images chacune d'un même fond uniforme, les deux séries étant acquises pour deux focales différentes, les M images de défaut étant obtenues par différences des images de chaque séries deux à deux. Selon un mode de réalisation, les tables d'étalonnage peuvent être celles citées dans le brevet US8645104B2, les figures 1 et 2 dudit brevet US permettant de générer des tables de correction appelées "tables for correction" qui correspondent aux tables d'étalonnage. Selon un mode de réalisation, pour chaque focale F' i , la table d'étalonnage T i associée à la focale F' i peut être déterminée selon une méthode de détermination de la table d'étalonnage, la méthode de détermination ayant un premier et un deuxième modes de réalisations distincts. La FIG. 4 présente un schéma synoptique des étapes du premier mode de réalisation la méthode 300 de détermination de la table d'étalonnage T i associée à la focale F' i . En particulier, chaque table d'étalonnage T i est associée à une température interne d'étalonnage T é talonnage de la caméra et un temps d'intégration d'étalonnage t é talonnage de la caméra. On entend par "température interne de la caméra" la température acquise près des optiques de l'objectif de la caméra et notamment près de la lentille frontale de l'objectif par exemple. Dans la suite, toute valeur de température interne de la caméra est comprise par exemple dans l'intervalle [-20; +70] °C. La méthode 300 comprend une étape 301 d'acquisition par la caméra d'une image I S_i d'un flux thermique émis par une source, pour la focale F' i de la caméra, la température interne d'étalonnage T é talonnage de la caméra et le temps d'intégration d'étalonnage t étalonnage de la caméra qui sont mémorisés par un utilisateur dans une mémoire externe à la caméra par exemple. La source émet par exemple un flux thermique (ou un rayonnement thermique) similaire aux flux thermique émis par un corps noir, c’est-à-dire un flux thermique isotrope, donc le rayonnement ne dépend que de la température et non des directions dans l'espace par exemple. Avantageusement, l'utilisation d'une source émettant un rayonnement isotrope permet d’avoir une scène homogène et de n’observer que les défauts optiques de la caméra. En pratique, la source peut être une plaque chauffante recouverte entièrement d'une peinture noire, et ayant une température uniforme, constante et modulable par un utilisateur. Ainsi, la source peut être assimilée à un corps noir. Le flux thermique émis par la source recouvre la totalité de la pupille de l'objectif de la caméra. La méthode 300 comprend une étape 302 de calibration de l'image I S_i selon la méthode de calibration un point de la caméra. La calibration un point est réalisée à partir de la table de gains G et de la table de décalages O décrits précédemment. Après l'étape de calibration, une image de la source calibrée I S 2 _i est obtenue à partir de l'image de la source I S_i , telle que : L'image I S 2 _i obtenue correspond à l'image de la source I S_i corrigée des défauts du détecteur et comprend des défauts optiques résiduels tels que cités précédemment, c’est-à-dire des défauts optiques résiduels correspondant à des flux thermiques parasites de l'objectif de la caméra relatifs à la focale F' i . La méthode 300 comprend en outre une étape 303 d'obtention d'une troisième image I S3 _i résultant de la soustraction d'une image moyenne de l'image I S 2 _i , notée Moyenne(I S2 _i ), à l'image calibrée I S2 _i . En particulier, Moyenne (I S2 _i ) correspond à une image de même taille que I S2 _i et dont chaque coefficient correspond à la moyenne de toutes les valeurs de pixels de I S2 _i. Ainsi, . Selon un premier mode de réalisation, la table d'étalonnage T i associée à la focale F' i correspond à I S 3_i et peut donc être notée: . Ainsi, la table d'étalonnage T i peut être obtenue à l'issue de l'étape 303. Selon un deuxième mode de réalisation, la méthode 300 comprend une étape 304 optionnelle de filtrage spatial de la troisième image I S 3_i pour obtenir une image filtrée F(I S 3_i ). Le filtrage peut par exemple correspondre à un filtrage spatial de l'image I S 3_i dans le cas où des défauts intrinsèques du détecteur subsistent même après l'étape 301 de calibration, par exemple à cause d'un ou plusieurs pixels du détecteur défectueux. Dans ce deuxième mode de réalisation, la table d'étalonnage T i peut être égale à F(I S 3_i ). Selon un troisième de réalisation, la table d'étalonnage T i est obtenue à l'issue d'étapes optionnelles non représentées de la méthode 300, décrites dans la suite. La méthode 300 peut comprendre une première étape optionnelle de réalisation des étapes 301, 302 et 303 à plusieurs instants t différents et d'obtention d'une pluralité d'images (I S3_it ) t>0 . La méthode 300 peut comprendre une deuxième étape optionnelle de moyennage des images I S3i_t pour obtenir une image moyenne Moyenne(I S3_it ) t>0 . En particulier, l'étape 306 de moyennage correspond à un moyennage de chaque valeur de pixel de chaque image I S3_it ayant la même position dans l'image. Autrement dit, chaque coefficient à une position donnée de l'image Moyenne(I S 3_it ) t >0 correspond à une moyenne des valeurs des pixels de même position de chaque image I S 3_it . La méthode 300 peut comprendre une troisième étape optionnelle d'obtention de la table d'étalonnage T i par soustraction de l'image Moyenne(I S3_it ) à l'image I S3_it . Ainsi, . En particulier, le troisième mode de réalisation permet avantageusement de supprimer un bruit temporel présent dans l'image I S3_i , le bruit temporel dégradant la qualité de l'image en lui donnant un aspect granulaire. Selon un mode de réalisation non représenté, les étapes optionnelles sont appliquées à l'image F(IS3_i) suite à l'étape 304. La méthode 300 d'obtention de la table d'étalonnage T i associée à la focale F' i peut être réalisée pour une pluralité de M focales , ce qui permet d'obtenir M premières tables d'étalonnage associées chacune à une distance focale différente F' i , l'indice i étant un entier naturel. En particulier, chaque table d'étalonnage parmi les M premières tables d'étalonnage est obtenue pour une distance focale différente mais avec la même température interne d'étalonnage T étalonnage de la caméra et le même temps d'intégration d'étalonnage t étalonnage de la caméra. En particulier, chaque table d'étalonnage peut être compressée avant d'être sauvegardée dans une mémoire pour être utilisée. La méthode de compression doit être adaptée à la morphologie spatiale des défauts optique résiduels représentés par chaque table d'étalonnage, afin d'éviter toute dégradation de la qualité de la table d'étalonnage. La FIG. 5 présente un schéma synoptique des étapes du deuxième mode de réalisation la méthode 400 de détermination de la table d'étalonnage T i associée à la focale F' i . La méthode 400 comprend une étape 401 d'acquisition par la caméra d'une première image I1 S_i d'un flux thermique émis par la source définie dans la méthode de détermination de la table d'étalonnage. En particulier, l'étape 401 d'acquisition est réalisée pour la focale F' i de la caméra, à une température interne d’étalonnage T1 etalonnage de la caméra et avec le temps d'intégration d'étalonnage t étalonnage de la caméra qui sont mémorisés par un utilisateur dans une mémoire externe à la caméra par exemple ou dans une mémoire de la caméra. La méthode 400 comprend en outre une étape 402 identique à l'étape 302, l'étape 402 étant appliquée à l'image I1 S_i et permettant d'obtenir l'image calibrée I1 S 2 _i . L'image calibrée I1 S 2 _i correspond à l'image I1 S_i corrigée des défauts du détecteur et comprenant des défauts optiques résiduels tels que décrits précédemment. Les défauts optiques résiduels correspondent à des flux thermiques parasites de l'objectif de la caméra pour la focale F' i . La méthode 400 comprend en outre une étape 403 d'obtention d'une image I1 S3 _i résultant de la soustraction d'une image moyenne de l'image I1 S 2 _i , notée Moyenne(I1 S2 _i ) à l'image calibrée I1 S2 _i . En particulier, Moyenne (I1 S2 _i ) correspond à une image de même taille que I1 S2 _i et dont chaque coefficient correspond à la moyenne de toutes les valeurs de pixels de I1 S2 _i. Ainsi, I1 S 3_i = . La méthode 400 comprend en outre une étape 404 de variation de la température interne de la caméra, pour obtenir une température interne T2 étalonnage de la caméra différente de T1 etalonnage . La variation de température interne de la caméra peut être réalisée en plaçant la caméra dans une étuve réglable en température par exemple. La méthode 400 comprend une étape 405 d'acquisition par la caméra d'une deuxième image I2 S_i de la source de flux thermique définie dans la méthode 200, pour la focale F' i de la caméra, la température interne d'étalonnage T2 etalonnage de la caméra et le temps d'intégration d'étalonnage t étalonnage de la caméra. La méthode 400 comprend en outre une étape 406 identique à l'étape 302, l'étape 406 étant appliquée à l'image I2 S_i et permettant d'obtenir l'image calibrée I2 S 2 _i . L'image calibrée I2 S 2 _i correspond à l'image I2 S_i corrigée des défauts du détecteur et comprenant des défauts optiques résiduels tels que cités précédemment. La méthode 400 comprend en outre une étape 407 d'obtention d'une image I' S3 _i résultant de la soustraction d'une image moyenne de l'image I2 S 2 _i , notée Moyenne(I2 S2 _i ) à l'image calibrée I2 S2 _i . En particulier, Moyenne (I2 S2 _i ) correspond à une image de même taille que I2 S2 _i et dont chaque coefficient correspond à la moyenne de toutes les valeurs de pixels de I2 S2 _i. Ainsi, . La méthode 400 comprend une étape 408 d'obtention de la table d'étalonnage T i à partir de l'image I1 S 3_i , la température d'étalonnage T1 étalonnage , l'image I2 S 3_i et la température d'étalonnage T2 étalonnage . En particulier, la table d'étalonnage T i permet d'estimer la variation des défauts optiques résiduels en fonction de la variation de température interne. Une interpolation est réalisée entre la variation des défauts optiques résiduels et la variation de température interne, l’ordre d’interpolation étant choisi selon des critères prédéterminés. En effet, plus l’ordre d'interpolation est élevé, et plus la qualité des images à corriger des défauts optiques résiduels est améliorée car la relation entre la variation des défauts optiques résiduels et la variation de température interne est plus précise. Cependant, plus l'ordre d'interpolation est élevé plus la méthode 400 est complexe à mettre en œuvre car une pluralité d'acquisition d'images du corps noir est requise. En particulier, la variation des défauts optiques résiduels peut dépendre de manière linéaire et indirecte de la variation de température interne. En effet, la variation des défauts optiques résiduels dépend linéairement de la variation de la luminance intégrée qui dépend de la température. Ainsi, la table d'étalonnage T i est obtenue selon la formule suivante: avec L(T1 étalonnage ) la luminance intégrée générée par la température T1 étalonnage et L(T2 étalonnage ) la luminance intégrée générée par la température T2 étalonnage . Ainsi, la luminance intégrée est liée à la température interne T que la caméra intègre (capture). Le niveau de gris capté par la caméra dépend d’une multitude de sources. Pour les besoins de l’interpolation, une simulation de la variation de niveaux de gris capté par le détecteur en fonction de la variation de température interne est réalisée. La détermination de la luminance en fonction de la température (notée L(T1) et L(T2)) est par exemple obtenue via une méthode de simulation de la caméra infrarouge. En particulier, la luminance intégrée prend en compte une température T fixe d'une scène intégrée par la caméra et traversant les optiques, la luminance interne de la caméra (c’est-à-dire les photons partant de la caméra et atteignant le détecteur) et le courant d'obscurité du détecteur de la caméra. Lorsque la caméra infrarouge est placée devant une scène de température T fixe, la température interne de la caméra est variée afin d'observer l'évolution de la variation de niveau de gris en fonction de la variation de la luminance interne de la caméra. Les luminances intégrées sont obtenues à partir d'une résolution des équations de Planck dans un contexte appliqué à la caméra, c’est-à-dire en prenant en compte la variation de niveau de gris en fonction de la variation de la luminance interne de la caméra. L'utilisation des valeurs de luminances en fonction de la température est avantageuse car la luminance d'un pixel d'une image est modélisée linéairement en fonction de son niveau de gris, pour n'importe quelle température. La méthode 400 d'obtention de la table d'étalonnage T i associée à la focale F' i peut être réalisée pour une pluralité de M focales , ce qui permet d'obtenir M deuxièmes tables d'étalonnage associées chacune à une distance focale différente F' i . En particulier, chaque table d'étalonnage parmi les M premières tables d'étalonnage est obtenue pour une distance focale différente mais avec le même temps d'intégration d'étalonnage t étalonnage de la caméra. Un autre aspect de l'invention concerne un dispositif configuré pour mettre en œuvre la méthode selon l'invention. La FIG. 6 représente le dispositif 500 selon un ou plusieurs modes de réalisation de l’invention. Dans ces modes de réalisation, le dispositif comporte un ordinateur 500, comprenant une mémoire 501 pour stocker images de défaut et/ou les séquences d'images par exemple. L'ordinateur 500 comporte en outre un circuit 502. Ce circuit peut être, par exemple un processeur apte à interpréter des instructions sous la forme de programme informatique, une carte électronique dont les étapes de la méthode de l'invention sont décrites dans le silicium, ou encore une puce électronique programmable comme une puce FPGA (pour « Field-Programmable Gate Array » en anglais). En particulier, les étapes 201 à 205 peuvent être programmées pour être mises en œuvre sur la puce FPGA. L’ordinateur 500 comporte une interface d'entrée 503 pour la réception de séquences d'images acquises par la caméra, et une interface de sortie 504. Enfin, l'ordinateur peut comporter, pour permettre une interaction aisée avec un utilisateur, un écran 505 affichant les séquences d'images acquises par la caméra, les tables d'étalonnage et les images corrigées par exemple. L'ordinateur 500 peut également comporter un clavier 506. Bien entendu, le clavier est facultatif, notamment dans le cadre d'un ordinateur ayant la forme d'une tablette tactile, par exemple. Bien entendu, la présente invention ne se limite pas aux formes de réalisation décrites ci-avant à titre d'exemples ; elle s'étend à d'autres variantes. Méthode (200), mise en œuvre par ordinateur, de correction d'un défaut de vignettage dans une séquence d'images acquise par une caméra pour une distance focale F', la séquence d'images étant ordonnée temporellement et comprenant une pluralité de sous-séquences d'images ordonnées temporellement, la méthode comportant un ensemble d'itérations, chaque itération utilisant une sous-séquence d'images de la pluralité de sous-séquences d'images et une table d'étalonnage, chaque itération comprenant: Réduire (201) des dimensions d'une première image de la séquence d'images pour obtenir une première image réduite, la première image étant corrigée du défaut de vignettage à partir d'une correction du défaut de vignettage déterminée à l'itération précédente ou à partir d'une correction initiale; Réduire (202) des dimensions de la table d'étalonnage pour obtenir une table d'étalonnage réduite; Appliquer (203) un filtre radial passe-haut sur la première image réduite et sur la table d'étalonnage réduite pour obtenir une première image réduite filtrée et une table d'étalonnage réduite filtrée; Déterminer (204) un coefficient de corrélation entre la première image réduite filtrée et la table d'étalonnage réduite filtrée; Déterminer (205) une correction du défaut de vignettage à partir de la table d'étalonnage, du coefficient de corrélation associé à la table d'étalonnage et de la correction du défaut de vignettage obtenue à l'itération précédente; Corriger (206) chaque autre image de la sous-séquence d'images à partir de la correction du défaut déterminée pour obtenir chaque autre image corrigée du défaut de vignettage. Méthode (200) selon la revendication précédente selon laquelle la première image d'une sous-séquence d'images utilisée à une itération donnée correspond à la dernière image d'une sous-séquence d'images utilisée à l'itération précédente. Méthode (200) selon l'une des revendications précédentes selon laquelle la méthode (200) est réalisée pour un ensemble de K tables d'étalonnage, K étant un entier naturel non nul, et pour chaque itération donnée de l'ensemble d'itérations, la table d'étalonnage utilisée est distincte de la table d'étalonnage utilisée lors de l'itération précédente. Méthode (200) selon la revendication précédente selon laquelle l'ensemble des K tables d'étalonnage est choisi selon une règle prédéfinie parmi M tables d'étalonnage, M étant un entier naturel non nul avec K inférieur à M, la règle prédéfinie étant telle que, pour chaque table d'étalonnage parmi les K tables d'étalonnage: la distance focale associée à la table d'étalonnage correspond à l'une des K distances focales les plus proches de F'; ou L'écart entre la distance focale associée à la table d'étalonnage et F' est inférieur à une distance focale moyenne entre au moins deux tables d'étalonnage successives parmi les K tables d'étalonnage choisies. Méthode (200) selon l'une des revendications 3 à 4 selon laquelle K est compris entre 2 et 10. Méthode (200) selon l'une des revendications précédentes selon laquelle le défaut de vignettage est un bruit additif spatial, ayant une composante radiale et de basse fréquence. Méthode (200) selon l'une des revendications précédentes selon laquelle la réduction de dimensions de la première image et de la table d'étalonnage est réalisée par un filtre moyenneur. Méthode (200) selon l'une des revendications précédentes selon laquelle le filtre radial est un filtre gradient ou un filtre Laplacien orienté dans la direction du centre de chaque image de la sous-séquence d'images. Méthode (200) selon l'une des revendications précédentes selon laquelle, le coefficient de corrélation entre la première image réduite filtrée et la table d'étalonnage réduite filtrée est déterminé à partir du rapport: d'un produit scalaire canonique entre la première image réduite filtrée et de la table d'étalonnage réduite filtrée; et d'une norme de la table d'étalonnage réduite filtrée. Produit programme d'ordinateur comprenant des instructions qui, lorsque le programme est exécuté par un ordinateur, conduisent celui-ci à mettre en œuvre la méthode selon l'une des revendications précédentes. Support d'enregistrement lisible par ordinateur comprenant des instructions qui, lorsqu'elles sont exécutées par un ordinateur, conduisent celui-ci à mettre en œuvre la méthode selon l'une des revendications 1 à 9.