FR 3138230 A1 20240126 FR 2207574 A 20220722 L’invention concerne un dispositif (5) de formation d’une image haute résolution à partir de couples formés d’une image basse résolution et d’une durée d’exposition comprenant : - un estimateur (17) agencé pour estimer, avec un algorithme Lucas-Kanade, un paramètre de distorsion pour chaque image basse résolution, - une unité de traitement (21) agencé pour obtenir une image haute résolution initiale, - une unité d’alignement (23) agencé pour aligner une image basse résolution de référence en fonction de chaque paramètre de distorsion, - un réseau de neurones convolutif (25) agencé pour générer, pour chaque image basse résolution, un facteur de confiance, - un calculateur (27) agencé pour calculer, pour chaque image basse résolution, un poids de déparasitage. - un module d’optimisation (13) agencé pour minimiser une fonction d’énergie quadratique à deux variables : une variable d’estimation de l’image haute résolution et une variable auxiliaire initialisées par l’image haute résolution initiale. [Fig. 1] Dispositif et procédé de formation d’une image à faible niveau de bruit à partir d’une rafale d’images Le domaine de l’invention se rapporte au traitement d’une rafale d’images pour la formation d’une image présentant un niveau de bruit inférieur à celui des images de la rafale d’images et une résolution spatiale supérieure. Les photographies prises avec un appareil photographique de milieu de gamme ou un téléphone intelligent (aussi désigné par le terme anglophone « smartphone ») présentent généralement un niveau de bruit élevé et une faible résolution spatiale. Ces caractéristiques ont des répercussions sur le niveau de détail de la photographie : un niveau de bruit faible permet de distinguer les détails dans les zones sombres tandis qu’une résolution spatiale élevée permet de zoomer sur la photographie. Par ailleurs, la qualité de la photographie peut être altérée par des artefacts colorés causés par des couleurs vives. Une solution possible, pour un capteur photographique de taille donnée, est de diminuer la taille des pixels pour augmenter la résolution spatiale. La contrepartie est cependant une augmentation du bruit dans les zones sombres causée par le fait que chaque pixel reçoit une plus faible quantité de photons. Ce phénomène est particulièrement problématique avec les téléphones intelligents dont le capteur photographique est généralement de taille réduite. Il est également connu d’utiliser plusieurs photographies d’une même scène – on parle alors d’une rafale d’images (aussi connue sous le terme anglophone « burst ») – pour construire une photographie d’une qualité supérieure, c’est-à-dire un niveau de bruit plus faible et une meilleure résolution spatiale. Les images de la rafale d’images peuvent être exploitées sous la forme de fichiers RAW – aussi appelés images brutes –, c’est-à-dire des fichiers d’images numériques contenant les données brutes pour la formation d’une image visible. Ce principe est notamment décrit dans l’article de B. Wronski et al. (2019) : « Handheld multi-frame super-resolution » (ACM Transactions On Graphics, vol. 38, No. 4, art. 28, p. 1-18) pour augmenter la résolution spatiale. Plus récemment, les inventeurs de la présente invention, B. Lecouat, J. Ponce et J. Mairal, ont proposé, dans une publication (2021) : « Lucas-Kanade reloaded: end-to-end super resolution from raw image burst » (International Conference on Computer Vision (ICCV)), une solution qui permet d’aligner les images brutes avec une précision sous-pixellique, de rester aussi fidèle que possible aux données acquises par le capteur photographique et de concevoir une fonction de régularisation (aussi connue sous le terme anglophone « regularizer ») performante. Cette solution présente néanmoins plusieurs inconvénients. En particulier, elle ne tient pas compte de l’hétérogénéité des images brutes, c’est-à-dire leurs durées d’exposition respectives, et ne prévoit pas non plus de déparasitage (aussi désigné par le terme anglophone « deghosting »). Par ailleurs, l’optimisation proposée pour minimiser une fonction d’énergie et résoudre un problème inverse exige des ressources de calcul importantes et présente un coût de temps élevé. La présente invention vient améliorer la situation. À ce titre, l’invention concerne un dispositif de formation d’une image d’une première résolution spatiale à partir d’une pluralité de couples formés chacun d’une image d’une deuxième résolution spatiale inférieure ou égale à la première résolution spatiale et d’une durée d’exposition. Le dispositif comprend : - un module d’initialisation agencé pour recevoir les images de deuxième résolution spatiale et comprenant : • un estimateur agencé pour estimer, pour chaque image de deuxième résolution spatiale, un paramètre de distorsion en appliquant un algorithme Lucas-Kanade, • une unité de traitement agencée pour appliquer, au moins à une image de deuxième résolution spatiale de référence parmi les images de deuxième résolution spatiale, un dématriçage et, lorsque la deuxième résolution spatiale est inférieure à la première résolution spatiale, une mise à l’échelle pour obtenir une image de première résolution spatiale initiale, et - un module d’optimisation agencé pour recevoir les images de deuxième résolution spatiale, les paramètres de distorsion associés et l’image de première résolution spatiale initiale et pour minimiser une fonction d’énergie quadratique paramétrée par une variable d’estimation de l’image de première résolution spatiale et une variable auxiliaire, la variable d’estimation de l’image de première résolution spatiale et la variable auxiliaire étant initialisées sur la base de l’image de première résolution spatiale initiale. La minimisation est réalisée par itérations successives comprenant chacune la mise à jour de la variable auxiliaire par une descente de gradient en fonction des images de deuxième résolution spatiale et des paramètres de distorsion associés et la mise à jour de la variable d’estimation de l’image de première résolution spatiale en fonction de la variable auxiliaire mise à jour. Le dispositif est caractérisé en ce que le module d’initialisation est agencé pour recevoir en outre les durées d’exposition respectives des images de deuxième résolution spatiale et comprend en outre : • une unité d’alignement agencée pour aligner par interpolation l’image de deuxième résolution spatiale de référence en fonction de chaque paramètre de distorsion, • un réseau de neurones convolutif agencé pour recevoir, en entrée, chaque image de deuxième résolution spatiale et l’image de deuxième résolution spatiale de référence alignée en fonction du paramètre de distorsion associé et pour générer, en sortie, un facteur de confiance, • un calculateur agencé pour calculer, pour chaque image de deuxième résolution spatiale, un poids de déparasitage en fonction de la durée d’exposition correspondante et du facteur de confiance associé, et en ce que le module d’optimisation est agencé pour recevoir en outre les poids de déparasitage respectivement associés aux images de deuxième résolution spatiale et comprend : • une unité de mise à jour agencée pour mettre à jour la variable auxiliaire par descente de gradient en fonction en outre des poids de déparasitage. Dans un ou plusieurs modes de réalisation, le module d’initialisation comprend en outre : • un réseau de neurones convolutif agencé pour recevoir, en entrée, les images de deuxième résolution spatiale et pour générer, en sortie, une carte de caractéristiques pour chaque image de deuxième résolution spatiale, et en ce que l’estimateur est agencé pour appliquer l’algorithme Lucas-Kanade à chaque carte de caractéristiques pour estimer le paramètre de distorsion associé à chaque image de deuxième résolution spatiale. Avantageusement, l’estimateur est agencé pour traiter chaque carte de caractéristiques, pour l’application de l’algorithme Lucas-Kanade, sous la forme d’une pyramide gaussienne. Dans un ou plusieurs modes de réalisation, le module d’optimisation comprend en outre : • un réseau de neurones convolutif agencé pour recevoir, en entrée, la variable auxiliaire mise à jour et pour générer, en sortie, la variable d’estimation de l’image de première résolution spatiale mise à jour. Avantageusement, le réseau de neurones convolutif du module d’optimisation est agencé pour mettre en œuvre un partage de poids d’une itération à l’autre. Par exemple, le calculateur est agencé pour calculer chaque poids de déparasitage comme suit : où : - y k est la k-ème image de deuxième résolution spatiale, - w k est le poids de déparasitage calculé pour l’image de deuxième résolution spatiale y k , - Δt k est la durée d’exposition de l’image de deuxième résolution spatiale y k , - K est le nombre d’images de deuxième résolution spatiale, - m c (·) est un masque binaire remplaçant par zéro toute valeur de pixel inférieure ou égale à un seuil c, - g ρ est le facteur de confiance paramétré par un paramètre ρ du réseau de neurones convolutif, - W k est un opérateur de distorsion caractérisé par le paramètre de distorsion associé à l’image de deuxième résolution spatiale y k , et - y k0 est l’image de deuxième résolution spatiale de référence. Avantageusement, l’unité d’alignement est agencée pour aligner l’image de deuxième résolution spatiale de référence en fonction de chaque paramètre de distorsion par interpolation bilinéaire. Dans un ou plusieurs modes de réalisation, l’unité de traitement est agencé pour dématricer chaque image de deuxième résolution spatiale, aligner par interpolation chaque image de deuxième résolution spatiale dématricée en fonction du paramètre de distorsion associé, former une image de deuxième résolution spatiale par moyennage des images de deuxième résolution spatiale dématricées et alignées et, lorsque la deuxième résolution spatiale est inférieure à la première résolution spatiale, mettre à l’échelle l’image de deuxième résolution spatiale formée pour obtenir l’image de première résolution spatiale initiale. Par exemple, l’unité de mise à jour est agencée pour mettre à jour la variable auxiliaire par descente de gradient comme suit : avec : où : - z est la variable auxiliaire, - x est la variable d’estimation de l’image de première résolution spatiale, - y k est la k-ème image de deuxième résolution spatiale, - K est le nombre d’images de deuxième résolution spatiale, - w k est le poids de déparasitage associé à l’image de deuxième résolution spatiale y k , - Δt k est la durée d’exposition de l’image de deuxième résolution spatiale y k , - W k est un opérateur de distorsion caractérisé par le paramètre de distorsion associé à l’image de deuxième résolution spatiale y k , - C est un masque binaire modélisant une matrice de filtres colorés, - s est le facteur de mise à l’échelle utilisé pour obtenir l’image de première résolution spatiale initiale, - D s est un opérateur de décimation caractérisé par le facteur s de mise à l’échelle, - B est un opérateur de floutage, - δ est un paramètre de la descente de gradient, et - η est un paramètre de la fonction d’énergie quadratique à deux variables. Typiquement, le module d’optimisation est agencé pour implémenter trois itérations successives. Dans un ou plusieurs modes de réalisation, au moins un réseau de neurones convolutif présente une architecture de type U-net. L’invention concerne également un procédé de formation d’une image d’une première résolution spatiale à partir d’une pluralité de couples formés chacun d’une image d’une deuxième résolution spatiale inférieure ou égale à la première résolution spatiale et d’une durée d’exposition. Le procédé est mise en œuvre par le dispositif décrit précédemment et comprend : - une phase d’initialisation comprenant les opérations suivantes : recevoir les images de deuxième résolution spatiale, estimer, pour chaque image de deuxième résolution spatiale, un paramètre de distorsion en appliquant un algorithme Lucas-Kanade, appliquer, au moins à une image de deuxième résolution spatiale de référence parmi les images de deuxième résolution spatiale, un dématriçage et, lorsque la deuxième résolution spatiale est inférieure à la première résolution spatiale, une mise à l’échelle pour obtenir une image de première résolution spatiale initiale, et - une phase d’optimisation comprenant les opérations suivantes : recevoir les images de deuxième résolution spatiale, les paramètres de distorsion associés et l’image de première résolution spatiale initiale, minimiser une fonction d’énergie quadratique paramétrée par une variable d’estimation de l’image de première résolution spatiale et une variable auxiliaire, la variable d’estimation de l’image de première résolution spatiale et la variable auxiliaire étant initialisées sur la base de l’image de première résolution spatiale initiale. La minimisation est réalisée par itérations successives comprenant chacune la mise à jour de la variable auxiliaire par une descente de gradient en fonction des images de deuxième résolution spatiale et des paramètres de distorsion associés et la mise à jour de la variable d’estimation de l’image de première résolution spatiale en fonction de la variable auxiliaire mise à jour. Le procédé est caractérisé en ce que la phase d’initialisation comprend les opérations suivantes : recevoir en outre les durées d’exposition respectives des images de deuxième résolution spatiale, aligner par interpolation l’image de deuxième résolution spatiale de référence en fonction de chaque paramètre de distorsion, utiliser un réseau de neurones convolutif recevant, en entrée, chaque image de deuxième résolution spatiale et l’image de deuxième résolution spatiale de référence alignée en fonction du paramètre de distorsion associé et générant, en sortie, un facteur de confiance, calculer, pour chaque image de deuxième résolution spatiale, un poids de déparasitage en fonction de la durée d’exposition correspondante et du facteur de confiance associé, et en ce que la phase d’optimisation comprend les opérations suivantes : recevoir en outre les poids de déparasitage respectivement associés aux images de deuxième résolution spatiale, mettre à jour la variable auxiliaire par descente de gradient en fonction en outre des poids de déparasitage. Dans un ou plusieurs modes de réalisation, le procédé comprend en outre une phase d’entraînement préalable du dispositif avec un ensemble de données d’apprentissage comprenant une image de première résolution spatiale et une pluralité de couples formés chacun d’une image de deuxième résolution spatiale et d’une durée d’exposition, chaque image de deuxième résolution spatiale étant générée à partir de l’image de première résolution spatiale comme suit : avec : où : - x (1) est l’image de première résolution spatiale de l’ensemble de données d’apprentissage, - y k (1) est la k-ème image de deuxième résolution spatiale de l’ensemble de données d’apprentissage, - Δt k (1) est la durée d’exposition de l’image de deuxième résolution spatiale y k (1) , - W k (1) est un opérateur de distorsion associé à l’image de deuxième résolution spatiale y k (1) , - C (1) est un masque binaire modélisant une matrice de filtres colorés, - D s (1) est un opérateur de décimation caractérisé par un facteur s de mise à l’échelle, - B (1) est un opérateur de floutage, - ε k (1) est un bruit associé à l’image de deuxième résolution spatiale y k (1) . Avantageusement, la phase d’entraînement préalable est supervisée en utilisant plusieurs ensembles de données d’apprentissage comprenant chacun une image de première résolution spatiale et une pluralité de couples formés chacun d’une image de deuxième résolution spatiale et d’une durée d’exposition. La phase d’entraînement préalable supervisée permet de paramétrer le dispositif de manière à minimiser la quantité suivante : où : - est l’image de première résolution spatiale du i-ème ensemble de données d’apprentissage, - est l’image de première résolution spatiale estimée par le dispositif à partir du i-ème ensemble de données d’apprentissage, et - n est le nombre d’ensembles de données d’apprentissage. Enfin, l’invention concerne aussi un programme informatique comprenant des instructions dont l’exécution par au moins un processeur se traduit par la mise en œuvre du procédé décrit précédemment. D’autres caractéristiques, détails et avantages apparaîtront à la lecture de la description détaillée ci-après, et à l’analyse des dessins annexés sur lesquels : illustre un dispositif de formation d’une image d’une première résolution spatiale à partir d’une rafale d’images d’une deuxième résolution spatiale selon l’invention ; illustre un procédé de formation d’une image d’une première résolution spatiale à partir d’une rafale d’images d’une deuxième résolution spatiale selon l’invention ; illustre une phase d’initialisation du procédé de la ; illustre une phase d’optimisation du procédé de la ; et illustre une phase d’entraînement du procédé de la La illustre un système 1 comprenant un capteur photographique 3 et un système 5. Le capteur photographique 3 est agencé pour recevoir un rayonnement électromagnétique – et plus particulièrement la lumière visible – en provenance d’une scène et pour convertir ce rayonnement en un signal électrique. Ce signal électrique est destiné à être numérisé, par exemple avec un convertisseur analogique-numérique, puis amplifié et traité pour obtenir une image numérique. Le capteur photographique 3 constitue une surface photosensible formée par une matrice d’éléments actifs appelés photosites. Chaque photosite peut être assimilé à un capteur élémentaire – de type photodiode – agencé pour convertir la lumière reçue en un courant électrique. Chaque photosite fonctionne selon le principe de l’effet photoélectrique : les photons incidents arrachent des électrons des photosites. Par ailleurs, la matrice de photosites du capteur photographique 3 est couplée avec une matrice de filtres colorés (parfois désigné par l’acronyme anglophone CFA pour « color filter array »), c’est-à-dire une mosaïque de filtres de couleur disposée de sorte que chaque photosite est associé à un filtre de couleur. En effet, les photosites sont sensibles à l’intensité lumineuse et non à la couleur ; la matrice de filtres colorés permet ainsi de séparer les couleurs. Parmi les différentes matrices de filtres colorés connues, la matrice de Bayer est la plus utilisée. La matrice de Bayer est constituée à 50% de filtres verts, à 25% de filtres rouges et à 25% de filtres bleus. L’image analogique directement obtenue par le capteur photographique 3 prend la forme d’une carte d’éclairement lumineux (aussi connue sous les termes anglophones « irradiance map » ou « irradiance image »). Chaque pixel de cette carte – ou image – d’éclairement lumineux correspond à un photosite. Ainsi, la carte d’éclairement lumineux se présente sous la forme d’une matrice dont chaque pixel a une valeur réelle positive caractéristique de l’éclairement lumineux mesuré au niveau du photosite correspondant. Pour générer une image numérique, le capteur photographique 3 est agencé pour appliquer, pixel par pixel, une transformation S à la carte d’éclairement lumineux. La fonction S est une fonction qui, à une valeur réelle positive, associe un entier appartenant à un ensemble P q , c’est-à-dire l’ensemble des entiers compris entre 0 et 2 q -1. On dit alors que le capteur photographique 3 a une profondeur de q bits. En outre, il convient de prendre en compte la durée d’exposition de la carte d’éclairement lumineux considérée. La durée d’exposition – également appelée temps de pose ou vitesse d’obturation – est la durée pendant laquelle le capteur photographique 3 est exposée à la lumière. Chaque pixel de l’image numérique obtenu est calculé comme suit : où : - u est un pixel donné, - y(u) est la valeur du pixel u de l’image numérique, - x(u) est la valeur du pixel u de la carte d’éclairement lumineux, et - Δt est la durée d’exposition. La sortie du capteur photographique 3 est donc une image matricée en noir et blanc dont chaque pixel a une valeur entière caractérisant la quantité de photons collectée par le photosite correspondant. Typiquement, l’image numérique obtenue par le capteur photographique 3 prend la forme d’un fichier RAW. Le dispositif 5 est agencé pour former une image d’une première résolution spatiale à partir d’une pluralité de couples formés chacun d’une image d’une deuxième résolution spatiale et d’une durée d’exposition. Il doit être compris ici que, au sein d’un couple donné, la durée d’exposition correspond au temps de pose de l’image de deuxième résolution spatiale associée. Plus généralement, le dispositif 5 est agencé pour former, à partir d’une rafale d’images, une image dont la qualité – en particulier le niveau de détail – est supérieure à celle des images de la rafale d’images. Par « niveau de détail », on comprend ici que l’image formée par le dispositif 5 a un niveau de bruit inférieur, et optionnellement une résolution spatiale supérieure ou égale à celle des images de la rafale d’images ou une gamme dynamique plus élevée. En photographie, la gamme dynamique désigne le rapport entre l’intensité lumineuse la plus élevée et l’intensité lumineuse la plus faible qu’un capteur photographique peut capturer. La première résolution spatiale est par conséquent supérieure ou égale à la deuxième résolution spatiale. Il doit être noté en particulier que le dispositif 5 peut être utilisé simplement pour baisser le niveau de bruit sans augmenter la résolution spatiale, par exemple pour une rafale d’images d’une scène sombre ou nocturne, auquel cas la première résolution spatiale de l’image formée par le dispositif 5 est égale à la deuxième résolution spatiale des images de la rafale d’images. Par souci de simplification, l’image de première résolution spatiale est appelée « image haute résolution », tandis que les images de deuxième résolution spatiale sont appelées « images basse résolution » dans la suite de la description. Par ailleurs, concernant les formats respectifs de l’image haute résolution et de l’image basse résolution, le dispositif 5 peut être configuré de sorte que les images basse résolution reçues sont des fichiers RAW – ou images brutes – alors que l’image haute résolution formée est une carte – ou image – d’éclairement lumineux. Il convient également de noter que les images basse résolution et leurs durées d’exposition respectives peuvent être reçues par le dispositif 5 en provenance du capteur photographique 3, auquel cas le système 1 forme un système de capture d’une rafale d’images basse résolution et de formation, à partir de cette rafale d’images basse résolution, d’une image haute résolution. Le dispositif 5 peut également être agencé pour recevoir les images basse résolution en provenance d’une pluralité de capteurs photographiques tels que le capteur photographique 3 ou encore recevoir les images basse résolution depuis une mémoire ou une plateforme de stockage quelconque. On comprend ainsi que le dispositif 5 peut être totalement dissocié du capteur photographique 3. Le dispositif 5 est donc agencé pour recevoir un ensemble Y d’images basse résolution y k et un ensemble Δ de durées d’exposition Δt k , et pour former une image haute résolution . Ici, k est un entier naturel compris entre 1 et K, où K est le nombre de couples formés d’une image basse résolution y k et d’une durée d’exposition Δt k associée. Parmi les images basse résolution y k , on distingue par ailleurs une image basse résolution de référence y k0 par rapport à laquelle l’alignement de chacune des autres images basse résolution y k est destiné à être évalué. En général, l’image basse résolution de référence y k0 est l’image basse résolution centrale de la rafale d’images y k . Toutefois, l’image basse résolution de référence y k0 peut également être la première image basse résolution reçue par le dispositif 5. Les images basse résolution y k d’une même scène sont prises en rafale et peuvent donc être chronologiquement ordonnées, en utilisant par exemple un marqueur temporel attaché à chaque image basse résolution y k . On comprend ici que les expressions « image basse résolution centrale » et « première image basse résolution » sont relatives à l’ordre chronologique des images basse résolution y k . Comme illustré sur la , le dispositif 5 comprend une mémoire 7, un processeur 9, un module d’initialisation 11 et un module d’optimisation 13. La mémoire 7 est agencée pour stocker des instructions dont la mise en œuvre, par le processeur 9, se traduit par le fonctionnement du dispositif 5. En particulier, de telles instructions permettent de mettre en œuvre la formation d’une image haute résolution à partir d’une pluralité de couples formés chacun d’une image basse résolution y k et d’une durée d’exposition associée Δt k . La mémoire 7 peut également être agencée pour stocker les images basse résolution reçues par le dispositif 5 et leurs durées d’exposition respectives. La mémoire 7 peut désigner tout support de stockage de données agencé pour recevoir et conserver des données numériques, par exemple un disque dur, un disque à semi-conducteurs (aussi désigné par l’acronyme anglophone SSD pour « solid-state drive ») ou plus généralement tout matériel informatique permettant le stockage de données sur de la mémoire flash. La mémoire 7 peut également être une mémoire vive ou encore un disque magnéto-optique. Une combinaison de plusieurs types de stockages de données peut aussi être envisagée. Le processeur 9 peut être réalisé de toute manière connue, par exemple sous la forme d’un microprocesseur, d’un circuit logique programmable (aussi désigné par l’acronyme anglophone PLD pour « Programmable Logical Device ») ou d’une puce dédiée de type FPGA (acronyme anglophone pour « Field Programmable Gate Array ») ou SoC (acronyme anglophone pour « System on Chip »), d’une grille de ressources informatiques, d’un microcontrôleur ou de toute autre forme propre disposant de la puissance de calcul nécessaire à la mise en œuvre du procédé de formation d’une image haute résolution décrit ci-après. Un ou plusieurs de ces éléments peuvent également être réalisés sous la forme de circuits électroniques spécialisés tel un ASIC (acronyme anglophone pour « Application-Specific Integrated Circuit »). Une combinaison de processeurs et de circuits électroniques peut aussi être envisagée. Un procédé de formation d’une image haute résolution à partir d’une pluralité de couples formés chacun d’une image basse résolution y k et d’une durée d’exposition Δt k va à présent être décrit en référence à la . Comme illustré sur la , ce procédé, mis en œuvre par le dispositif 5, comprend une phase d’entraînement 200, une phase d’initialisation 210 et une phase d’optimisation 220. La phase d’entraînement 200 sera abordée en fin de description. La phase d’initialisation 210 est mise en œuvre par le module d’initialisation 11 et est illustrée plus en détail sur la . Comme illustré sur la , le module d’initialisation 11 comprend un réseau de neurones convolutif 15, un estimateur 17, une unité de déparasitage 19 et une unité de traitement 21. Lors d’une opération 300, le module d’initialisation 11 reçoit une pluralité de couples formés chacun d’une image basse résolution y k et d’une durée d’exposition Δt k associée. L’ensemble Y d’images basse résolution y k et l’ensemble Δ de durées d’exposition Δt k sont par exemple transmis au dispositif 5 par le capteur photographique 3. Lors d’une opération 310, l’estimateur 17 estime, pour chaque image basse résolution y k , un paramètre de distorsion p k en appliquant un algorithme Lucas-Kanade. La méthode Lucas-Kanade est une méthode différentielle connue de l’homme du métier et qui permet d’estimer, pour une image donnée, un flux optique caractérisant un mouvement apparent d’un objet présent dans la scène capturée et causé par un mouvement du capteur photographique utilisé. Cette méthode repose sur le principe que, dans le voisinage d’un pixel donné de l’image, le flux optique est constant. L’équation du flux optique est généralement résolue, au sein du voisinage considéré, par une méthode des moindres carrés. Dans le cadre de l’invention, l’estimateur 17 compare chaque image basse résolution y k à l’image basse résolution de référence y k0 pour l’estimation du paramètre de distorsion p k qui correspond donc à un flux optique entre l’image basse résolution y k et l’image basse résolution de référence y k0 . En d’autres termes, l’estimateur 17 estime le flux optique p k entre chaque image basse résolution y k et l’image basse résolution de référence y k0 . Les paramètres de distorsion p k sont caractéristiques d’un opérateur de distorsion W k (aussi connu sous le terme anglophone « warping operator ») qui traduit les tremblements ou mouvements du capteur photographique lors de l’acquisition de image basse résolution y k . Étant donné que l’algorithme Lucas-Kanade prend, comme image de référence, l’image basse résolution de référence y k0 , on en déduit que l’opérateur de distorsion W k0 correspondant à l’image basse résolution de référence y k0 est l’identité. Dans le mode de réalisation décrit ici, l’opération 310 est mise en œuvre par l’implémentation d’opérations 311 et 312 décrites ci-après – et qui constituent en fait des sous-opérations de l’opération 310. Ce mode de réalisation utilise le réseau de neurones convolutif 15. Ainsi, lors de l’opération 311, le réseau de neurones convolutif 15 reçoit, en entrée, les images basse résolution y k et génère, en sortie, une carte de caractéristiques (aussi désignée par le terme anglophone « feature map ») FM k pour chaque image basse résolution. Lors de l’opération 312, l’estimateur 17 applique l’algorithme Lucas-Kanade à chaque carte de caractéristiques FM k pour estimer le paramètre de distorsion p k associé à chaque image basse résolution y k . En particulier, chaque carte de caractéristiques FM k peut être traitée par l’estimateur 17 sous la forme d’une pyramide gaussienne. En traitement d’image, une pyramide désigne une représentation multi-résolution d’une image dans laquelle l’image est présente en plusieurs résolutions spatiales distinctes – l’image de plus grande résolution spatiale formant la base de la pyramide, tandis que l’image de plus faible résolution spatiale en forme le sommet. Chaque résolution spatiale de l’image constitue ainsi un « étage » de la pyramide. Dans le cas en l’espèce, la pyramide gaussienne comprend par exemple six résolutions spatiales, donc six étages. Il est connu d’appliquer l’algorithme Lucas-Kanade à des images sous forme pyramidale. Un tel algorithme PLK (acronyme anglophone pour « Pyramidal Lucas-Kanade ») consiste alors, pour comparer deux images, à comparer les étages correspondants des pyramides représentant respectivement les images considérées. Typiquement, trois itérations de l’algorithme Lucas-Kanade sont implémentées pour chaque étage. En d’autres termes, dans le cas présent, l’estimateur 17 modélise chaque carte de caractéristiques FM k comme une pyramide gaussienne et applique un algorithme Lucas-Kanade pour déterminer le flux optique entre les résolutions spatiales correspondantes de la pyramide gaussienne modélisant la carte de caractéristiques FM k et de la pyramide gaussienne modélisant la carte de caractéristiques FM k0 . Il en résulte, après comparaison deux à deux des étages respectifs, le paramètre de distorsion p k de l’image basse résolution y k . Le paramètre de distorsion p k correspond donc à un flux optique entre la carte de caractéristiques FM k et la carte de caractéristiques FM k0 de l’image basse résolution de référence y k0 . En variante, l’estimateur 17 applique directement l’algorithme Lucas-Kanade aux images basse résolution y k pour estimer le paramètre de distorsion p k associé à chacune. On comprend donc que le réseau de neurones convolutif 15 est optionnel et que l’estimateur 17 peut estimer les paramètres de distorsion p k sans les cartes de caractéristiques FM k . Les paramètres de distorsion p k sont transmis à l’unité de déparasitage 19 et à l’unité de traitement 21. L’unité de déparasitage 19 est agencée pour générer, pour chaque image basse résolution y k , un poids de déparasitage w k . Le poids de déparasitage w k caractérise à la fois la durée d’exposition Δt k de l’image basse résolution y k et son alignement par rapport à l’image basse résolution de référence y k0 . Par conséquent, l’unité de déparasitage 19 permet de conférer à chaque image basse résolution y k , aux fins de la mise en œuvre de la phase d’optimisation 220 par le module d’optimisation 13, un poids d’autant plus élevé que l’image basse résolution y k considérée a une durée d’exposition Δt k élevée par rapport aux autres images basse résolution y k et un meilleur alignement relativement à l’image basse résolution de référence y k0 . Comme illustré sur la , l’unité de déparasitage 19 comprend une unité d’alignement 23, un réseau de neurones convolutif 25 et un calculateur 27. Lors d’une opération 320, l’unité d’alignement 23 aligne, par interpolation, l’image basse résolution de référence y k0 en fonction de chaque paramètre de distorsion p k . L’interpolation mise en œuvre par l’unité d’alignement 23 est par exemple une interpolation bilinéaire qui présente l’avantage d’être rapide et d’offrir une qualité satisfaisante. Toutefois, d’autres méthodes d’interpolation peuvent être utilisées par l’unité d’alignement 23, notamment une interpolation bicubique ou une interpolation spline. En d’autres termes, l’unité d’alignement 23 réalise, pour chaque paramètre de distorsion p k , un alignement de l’image basse résolution de référence y k0 . L’image basse résolution de référence y k0 alignée en fonction de l’opérateur de distorsion W k associé à l’image basse résolution y k est notée W k y k0 . Lors d’une opération 330, le réseau de neurones convolutif 25 reçoit, en entrée, chaque image basse résolution y k et l’image basse résolution de référence alignée W k y k0 en fonction du paramètre de distorsion p k associé et génère, en sortie, un facteur de confiance g ρ . On comprend ici qu’un facteur de confiance g ρ est généré pour chaque image basse résolution y k ou, plus exactement, pour chaque opérateur de distorsion W k . Lors d’une opération 340, le calculateur 27 calcule, pour chaque image basse résolution y k , un poids de déparasitage w k en fonction de la durée d’exposition Δt k correspondant et du facteur de confiance g ρ associé. Chaque poids de déparasitage w k peut être calculé comme suit où : - w k est le poids de déparasitage calculé pour l’image basse résolution y k , - m c (·) est un masque binaire remplaçant par zéro toute valeur de pixel inférieure ou égale à un seuil c, - g ρ est le facteur de confiance paramétré par un paramètre ρ du réseau de neurones convolutif 25. Le masque binaire m c permet d’exclure les pixels saturés de chaque image basse résolution y k . Il est à noter que, contrairement à la méthode proposée par les inventeurs B. Lecouat, J. Ponce et J. Mairal (2021), la présente invention effectue un déparasitage qui tient compte, le cas échéant, de l’hétérogénéité des images basse résolution. Les poids de déparasitage w k permettent en effet de donner une plus grande importance aux images basse résolution y k présentant une durée d’exposition Δt k plus élevée. Comme mentionné précédemment, les paramètres de distorsion p k respectifs des images basse résolution y k générés par l’estimateur 17 sont également transmis à l’unité de traitement 21. L’unité de traitement 21 est agencée pour former une image haute résolution initiale x 0 aux fins de la mise en œuvre de la phase d’optimisation 220 par le module d’optimisation 13. Cette image haute résolution initiale x 0 est destinée à constituer, lors de la phase d’optimisation 220, un point de départ pour la formation de l’image haute résolution . Lors d’une opération 350, l’unité de traitement 21 applique un dématriçage (plus connu sous le terme anglophone « demosaicking ») et une mise à l’échelle au moins à l’image basse résolution de référence (yk 0 ) pour obtenir une image haute résolution initiale x 0 . Plus exactement, une telle mise à l’échelle pour augmenter la résolution spatiale est appelée conversion ascendante. Le facteur de mise à l’échelle est noté s dans la suite de la description. Une telle conversion ascendante peut être réalisée par interpolation. Là encore, la méthode d’interpolation utilisée par l’unité de traitement 21 peut être une interpolation bilinéaire, une interpolation bicubique ou une interpolation spline. Il est à noter que la mise à l’échelle n’est pas nécessaire dans le cas où l’image formée par le dispositif 5 a la même résolution spatiale que les images de la rafale d’images reçues. Comme expliqué précédemment, un tel cas se produit notamment lorsque la rafale d’images concerne une scène sombre ou nocturne. Avantageusement, l’unité de traitement 21 forme l’image haute résolution initiale x 0 à partir de toutes les images basse résolution y k et non pas seulement de l’image basse résolution de référence y k0 . Dans un tel cas, l’opération 350 est mise en œuvre par l’unité de traitement 21 par l’implémentation d’opérations 351, 352, 353 et 354 décrites ci-après – et qui constituent en fait des sous-opérations de l’opération 350. Ainsi, lors de l’opération 351, l’unité de traitement 21 applique un dématriçage à chaque image basse résolution y k . Lors de l’opération 352, l’unité de traitement 21 aligne par interpolation chaque image basse résolution y k dématricée en fonction du paramètre de distorsion p k associé. De même que pour l’opération 320, plusieurs méthodes d’interpolation peuvent être utilisées par l’unité de traitement 21, notamment l’interpolation bilinéaire mais aussi l’interpolation bicubique ou l’interpolation spline. Lors de l’opération 353, l’unité de traitement 21 forme une image basse résolution par moyennage (aussi connu sous le terme anglophone « averaging ») des images basse résolution dématricées et alignées. Enfin, lors de l’opération 354, l’unité de traitement 21 met à l’échelle l’image basse résolution formée pour obtenir l’image haute résolution initiale x 0 . Là encore, la mise à l’échelle – ou conversion ascendante – n’est mise en œuvre que si le dispositif 5 est agencé pour former une image d’une résolution spatiale supérieure à celle des images de la rafale d’images reçues. En référence à nouveau à la , le procédé se poursuit avec une phase d’optimisation 220 mise en œuvre par le module d’optimisation 13 et illustrée plus en détail sur la . Le module d’optimisation 13 est agencé pour minimiser une fonction d’énergie quadratique à deux variables : une variable d’estimation de l’image de première résolution spatiale x et une variable auxiliaire z. La fonction d’énergie quadratique à deux variables à minimiser est par exemple la suivante : avec : où : - x est la variable d’estimation de l’image haute résolution, - z est la variable auxiliaire, -║·║ F est une norme, - η et λ sont des paramètres de la fonction d’énergie quadratique à deux variables, - Ω(·) est une fonction de régularisation, - C est un masque binaire modélisant une matrice de filtres colorés, - s est le facteur de mise à l’échelle appliqué pour obtenir l’image haute résolution initiale x 0 , - D s est un opérateur de décimation caractérisé par le facteur s de mise à l’échelle, et - B est un opérateur de floutage. Comme illustré sur la , le module d’optimisation 13 comprend une unité de mise à jour 29 et un réseau de neurones convolutif 31. Lors d’une opération 400, le module d’optimisation 13 reçoit les images basse résolution y k , les paramètres de distorsion p k associés, l’image haute résolution initiale x 0 et les poids de déparasitage w k respectivement associés aux images basse résolution y k . La suite de la phase d’optimisation 220 consiste à minimiser la fonction d’énergie quadratique à deux variables par itérations successives selon une méthode de type HQS (acronyme anglophone pour « half-quadratic splitting »). Chaque itération comprend la mise à jour de la variable auxiliaire z et la mise à jour de la variable d’estimation de l’image haute résolution x. Il est à noter que, dans la publication précédente des inventeurs B. Lecouat, J. Ponce et J. Mairal (2021), la fonction d’énergie quadratique à minimiser est une fonction à trois variables. En effet, en plus de la variable auxiliaire z et de la variable d’estimation de l’image haute résolution x, les paramètres de distorsion p k – rassemblés dans un vecteur p – sont également mis à jour selon une approche de Gauss-Newton. Chaque itération est ainsi alourdie d’une opération supplémentaire, ce qui demande des ressources de calcul supplémentaires et représente un surcoût de temps. En effet, dans la présente invention, les paramètres de distorsion p k ne varient pas pendant la phase d’optimisation : on utilise directement les paramètres de distorsion p k estimés par l’estimateur 17. Lors d’une opération 410, le module d’optimisation 13 initialise la variable auxiliaire z et la variable d’estimation de l’image haute résolution x avec l’image haute résolution initiale x 0 . Par ailleurs, une variable locale i, correspondant au nombre d’itérations, est initialisée à 1. Lors d’une opération 420, l’unité de mise à jour 29 met à jour la variable auxiliaire z par une descente de gradient en fonction des images basse résolution y k , des paramètres de distorsion associés p k et des poids de déparasitage w k . Par exemple, l’unité de mise à jour 29 met à jour la variable auxiliaire z par descente de gradient comme suit : où : δ est un paramètre de la descente de gradient. Lors d’une opération 430, le module d’optimisation 13 met à jour la variable d’estimation de l’image haute résolution x en fonction de la variable auxiliaire z mise à jour. Avantageusement, le réseau de neurones convolutif 31 reçoit, en entrée, la variable auxiliaire z mise à jour et génère, en sortie, la variable d’estimation de l’image haute résolution x mise à jour. Par ailleurs, le réseau de neurones convolutif 31 peut mettre en œuvre un partage de poids (plus connu sous le terme anglophone « weight-sharing ») d’une itération à l’autre. Le réseau de neurones convolutif 31 permet ainsi, à chaque itération, de réaliser un débruitage (aussi désigné par le terme anglophone « denoising ») de la variable d’estimation de l’image haute résolution x. Alternativement, il est aussi possible d’utiliser l’opérateur proximal G de la fonction de régularisation Ω pour mettre à jour la variable d’estimation de l’image haute résolution x et donc de ne pas avoir recours à un réseau de neurones convolutif tel que le réseau de neurones convolutif 31. Le recours à l’opérateur proximal G a notamment fait l’objet d’une publication de N. Parikh et S.P. Boyd (2014) : « Proximal A lgorithms » (Fondations and Trends in Optimization, vol. 1, No. 3, p. 127-239). Le débruitage peut ainsi être mis en œuvre par un filtrage bilatéral ou par un opérateur proximal d’une fonction de régularisation telle que la variation totale (parfois désignée par l’acronyme anglophone TVR pour « total variation regularization »). Toutefois, préférentiellement, le réseau de neurones convolutif 31 est sollicité pour mettre à jour la variable d’estimation de l’image haute résolution par un opérateur G ω vérifiant : où : - ω est un paramètre du réseau de neurones convolutif 31, et - γ est un paramètre vérifiant : . Lors d’une opération 440, le module d’optimisation 13 détermine si un nombre prédéterminé d’itérations – ou nombre de boucles – a été atteint. Dans l’exemple illustré sur la , ce nombre d’itérations est fixé à 3. Toutefois, le nombre d’itérations peut être strictement supérieur à 3. Si le nombre prédéterminé d’itérations n’a pas été atteint, la phase d’optimisation 220 continue avec une opération 450 au cours de laquelle la variable locale i est incrémentée d’un pas de 1. En revanche, si le nombre prédéterminé d’itérations a été atteint, la phase d’optimisation 220 se termine avec une opération 460 qui correspond à la sortie de l’image haute résolution estimée , c’est-à-dire la variable d’estimation de l’image haute résolution x mise à jour au cours de la dernière itération. L’image haute résolution estimée est donc l’image formée par le dispositif 5 à partir de la pluralité de couples formés chacun d’une image basse résolution y k et d’une durée d’exposition Δt k . Comme expliqué précédemment, le dispositif 5 a recours à un ou plusieurs réseaux de neurone convolutifs. Typiquement, au moins un réseau de neurones convolutif parmi les réseaux de neurones convolutifs 15, 25 et 31 présente une architecture de type U-net. Un U-net est un réseau de neurones dit « entièrement convolutif » généralement dédié à la vision par ordinateur (plus connue sous l’expression anglophone « Computer Vision ») et plus particulièrement à la segmentation sémantique, notamment dans le domaine biomédical. Un U-net comprend une partie de contraction – aussi appelée « encodeur » – constitué d’un assemblage de couches de convolution, d’unité linéaire rectifiée (« ReLU ») et de max-pooling pour récupérer le contexte d’une image et une partie d’expansion – aussi appelée « décodeur » – qui exploite les informations obtenues par la partie de contraction pour obtenir une localisation précise par convolution transposée et concaténation ascendante. La partie d’expansion permet ainsi de retrouver la taille initiale de l’image analysée. À titre d’exemple, K. Zhang et al. décrivent l’utilisation d’un réseau de neurones convolutif de type U-net dans la publication (2020) : « Deep unfolding network for image-super-resolution » (Proceedings of the Conference on Computer Vision and Pattern Recognition (CVPR), p. 3217-3226). Toutefois, d’autres architectures de réseaux de neurones convolutifs peuvent être utilisées. Par conséquent et comme illustré sur la , il peut être avantageux de mettre en œuvre une phase d’entraînement préalable 200 du dispositif 5 avec un ensemble de données d’apprentissage. Cet ensemble de données d’apprentissage comprend une image haute résolution et une pluralité de couples formés chacun d’une image basse résolution et d’une durée d’exposition. Par souci de simplification, les notations utilisées ci-après pour désigner les données d’apprentissage sont similaires à celles utilisées pour décrire la phase d’initialisation 210 et la phase d’optimisation 220. La phase d’entraînement préalable 200, et plus précisément la génération des données d’apprentissage, est illustrée sur la . La génération des données d’apprentissage est semi-synthétique. En d’autres termes, le principe est de partir d’une image haute résolution x (1) et de lui appliquer des opérations successives de manière à obtenir des images basse résolution y k (1) . Dans l’exemple illustré sur la , l’image haute résolution x (1) est une image de dimensions sh x sw x 3. Ici, la dimension 3 correspond aux canaux RVB (acronyme de « rouge, vert, bleu » ; ou RGB en anglais) : l’image haute résolution x (1) est une image en couleur. Les images basse résolution y k (1) sont ici des images de dimensions h x w. Il y a donc un facteur d’échelle s entre les images brutes y k (1) générées et l’image haute résolution x (1) de départ. Pour chaque image basse résolution y k (1) qu’on souhaite générer, on utilise en entrée une durée d’exposition Δt k (1) et un opérateur de distorsion W k (1) . L’opérateur de distorsion W k (1) permet de simuler des tremblements ou mouvements d’un capteur photographique. Par ailleurs, pour toutes les images basse résolution y k (1) qu’on souhaite générer, on utilise un opérateur de floutage B (1) (ou « blur » en anglais), un opérateur de décimation D s (1) pour réduire la résolution de l’image d’un facteur d’échelle s et un masque binaire C (1) modélisant une matrice de filtres colorés. Un bruit ε k (1) propre à chaque image basse résolution y k (1) est également ajouté à l’image obtenue. Chaque image basse résolution y k (1) est ainsi générée comme suit à partir de l’image haute résolution x : avec : où : - x (1) est l’image haute résolution de l’ensemble de données d’apprentissage, - y k (1) est la k-ème image basse résolution de l’ensemble de données d’apprentissage, - Δt k (1) est la durée d’exposition de l’image basse résolution y k (1) , - W k (1) est l’opérateur de distorsion associé à l’image basse résolution y k (1) , - C (1) est le masque binaire modélisant une matrice de filtres colorés, - D s (1) est l’opérateur de décimation caractérisé par le facteur s de mise à l’échelle, - B (1) est l’opérateur de floutage, et - ε k (1) est le bruit associé à l’image basse résolution y k (1) . Par ailleurs, lorsque l’image haute résolution x (1) est une carte d’éclairement et que les images basse résolution y k (1) générées sont des images brutes, il convient d’appliquer la fonction S décrite précédemment qui, à une valeur réelle positive, associe un entier compris entre 0 et 2 q -1, q étant un entier strictement positif. On a alors l’expression suivante : L’ensemble de données d’apprentissage ainsi généré est ensuite utilisé pour entraîner le dispositif 5. En d’autres termes, le dispositif 5 reçoit les images basse résolution y k (1) et les durées d’exposition Δt k (1) de l’ensemble de données d’apprentissage et leur appliquer la phase d’initialisation 210 et la phase d’optimisation 220 détaillées plus haut. L’image estimée est comparée à l’image x (1) de l’ensemble de données d’apprentissage – qui est donc connue –, ce qui permet de paramétrer le dispositif 5, et notamment les différents réseaux de neurones convolutifs utilisés. Les paramètres à apprendre sont par exemple les poids des couches du réseau de neurones convolutif 15, le paramètre ρ du facteur de confiance g ρ et les poids des couches du réseau de neurones convolutif 25, les paramètres λ et η de la fonction d’énergie quadratique minimisée par le module d’optimisation 13, le paramètre δ de la descente de gradient mise en œuvre par l’unité de mise à jour 29, les paramètres ω et γ et les poids des couches du réseau de neurones convolutif 31. En rassemblant l’ensemble de ces paramètres dans un vecteur θ, on peut modéliser le dispositif 5 par un moteur d’inférence agencé pour appliquer F θ à un ensemble Y (1) d’images basse résolution y k (1) et à un ensemble Δ (1) de durées d’exposition Δt k (1) associées et en tirer une image haute résolution estimée . En d’autres termes : Avantageusement, la phase d’entraînement préalable 200 est supervisée en utilisant plusieurs ensembles de données d’apprentissage comprenant chacun une image haute résolution x (i) et une pluralité de couples formés chacun d’une image basse résolution y k (i) et d’une durée d’exposition Δt k (i) . La phase d’entraînement préalable 200 supervisée permet de paramétrer le dispositif 5 de manière à minimiser la quantité suivante : où : - est l’image haute résolution du i-ème ensemble de données d’apprentissage, - est l’image haute résolution estimée par le dispositif 5 à partir du i-ème ensemble de données d’apprentissage, - Y (i) est l’ensemble des images basse résolution y k (i) du i-ème ensemble de données d’apprentissage, - Δ (i) est l’ensemble des durées d’exposition Δt k (i) du i-ème ensemble de données d’apprentissage, et - n est le nombre d’ensembles de données d’apprentissage. Dispositif (5) de formation d’une image d’une première résolution spatiale ( ) à partir d’une pluralité de couples formés chacun d’une image d’une deuxième résolution spatiale (y k ) inférieure ou égale à la première résolution spatiale et d’une durée d’exposition (Δt k ), ledit dispositif (5) comprenant : - un module d’initialisation (11) agencé pour recevoir les images de deuxième résolution spatiale (y k ) et comprenant : • un estimateur (17) agencé pour estimer, pour chaque image de deuxième résolution spatiale (y k ), un paramètre de distorsion en appliquant un algorithme Lucas-Kanade, • une unité de traitement (21) agencée pour appliquer, au moins à une image de deuxième résolution spatiale de référence (y k0 ) parmi les images de deuxième résolution spatiale (y k ), un dématriçage et, lorsque la deuxième résolution spatiale est inférieure à la première résolution spatiale, une mise à l’échelle pour obtenir une image de première résolution spatiale initiale (x 0 ), et - un module d’optimisation (13) agencé pour recevoir les images de deuxième résolution spatiale (y k ), les paramètres de distorsion (p k ) associés et l’image de première résolution spatiale initiale (x 0 ) et pour minimiser une fonction d’énergie quadratique paramétrée par une variable d’estimation de l’image de première résolution spatiale (x) et une variable auxiliaire (z), ladite variable d’estimation de l’image de première résolution spatiale (x) et ladite variable auxiliaire (z) étant initialisées sur la base de l’image de première résolution spatiale initiale (x 0 ), la minimisation étant réalisée par itérations successives comprenant chacune la mise à jour de la variable auxiliaire (z) par une descente de gradient en fonction des images de deuxième résolution spatiale (y k ) et des paramètres de distorsion (p k ) associés et la mise à jour de la variable d’estimation de l’image de première résolution spatiale (x) en fonction de la variable auxiliaire (z) mise à jour, ledit dispositif (5) étant caractérisé en ce que le module d’initialisation (11) est agencé pour recevoir en outre les durées d’exposition (Δt k ) respectives des images de deuxième résolution spatiale (y k ) et comprend en outre : • une unité d’alignement (23) agencée pour aligner par interpolation l’image de deuxième résolution spatiale de référence (y k0 ) en fonction de chaque paramètre de distorsion (p k ), • un réseau de neurones convolutif (25) agencé pour recevoir, en entrée, chaque image de deuxième résolution spatiale (y k ) et l’image de deuxième résolution spatiale de référence alignée (W k y k0 ) en fonction du paramètre de distorsion (p k ) associé et pour générer, en sortie, un facteur de confiance (g ρ ), • un calculateur (27) agencé pour calculer, pour chaque image de deuxième résolution spatiale (y k ), un poids de déparasitage (w k ) en fonction de la durée d’exposition (Δt k ) correspondante et du facteur de confiance (g ρ ) associé, et en ce que le module d’optimisation (13) est agencé pour recevoir en outre les poids de déparasitage (w k ) respectivement associés aux images de deuxième résolution spatiale (y k ) et comprend : • une unité de mise à jour (29) agencée pour mettre à jour la variable auxiliaire (z) par descente de gradient en fonction en outre des poids de déparasitage (w k ). Dispositif (5) selon la revendication 1, caractérisé en ce que le module d’initialisation (11) comprend en outre : • un réseau de neurones convolutif (15) agencé pour recevoir, en entrée, les images de deuxième résolution spatiale (y k ) et pour générer, en sortie, une carte de caractéristiques (FM k ) pour chaque image de deuxième résolution spatiale (y k ), et en ce que l’estimateur (17) est agencé pour appliquer l’algorithme Lucas-Kanade à chaque carte de caractéristiques (FM k ) pour estimer le paramètre de distorsion (p k ) associé à chaque image de deuxième résolution spatiale (y k ). Dispositif (5) selon la revendication 2, caractérisé en ce que l’estimateur (17) est agencé pour traiter chaque carte de caractéristiques (FM k ), pour l’application de l’algorithme Lucas-Kanade, sous la forme d’une pyramide gaussienne. Dispositif (5) selon l’une des revendications précédentes, caractérisé en ce que le module d’optimisation (13) comprend en outre : • un réseau de neurones convolutif (31) agencé pour recevoir, en entrée, la variable auxiliaire (z) mise à jour et pour générer, en sortie, la variable d’estimation de l’image de première résolution spatiale (x) mise à jour. Dispositif (5) selon la revendication 4, caractérisé en ce que le réseau de neurones convolutif (31) du module d’optimisation (13) est agencé pour mettre en œuvre un partage de poids d’une itération à l’autre. Dispositif (5) selon l’une des revendications précédentes, caractérisé en ce que le calculateur (27) est agencé pour calculer chaque poids de déparasitage (w k ) comme suit : où : - y k est la k-ème image de deuxième résolution spatiale, - w k est le poids de déparasitage calculé pour l’image de deuxième résolution spatiale y k , - Δt k est la durée d’exposition de l’image de deuxième résolution spatiale y k , - K est le nombre d’images de deuxième résolution spatiale, - m c (·) est un masque binaire remplaçant par zéro toute valeur de pixel inférieure ou égale à un seuil c, - g ρ est le facteur de confiance paramétré par un paramètre ρ du réseau de neurones convolutif (25), - W k est un opérateur de distorsion caractérisé par le paramètre de distorsion (p k ) associé à l’image de deuxième résolution spatiale y k , et - y k0 est l’image de deuxième résolution spatiale de référence. Dispositif (5) selon l’une des revendications précédentes, caractérisé en ce que l’unité d’alignement (23) est agencée pour aligner l’image de deuxième résolution spatiale de référence (y k0 ) en fonction de chaque paramètre de distorsion (p k ) par interpolation bilinéaire. Dispositif (5) selon l’une des revendications précédentes, caractérisé en ce que l’unité de traitement (21) est agencé pour dématricer chaque image de deuxième résolution spatiale (y k ), aligner par interpolation chaque image de deuxième résolution spatiale dématricée en fonction du paramètre de distorsion (p k ) associé, former une image de deuxième résolution spatiale par moyennage des images de deuxième résolution spatiale dématricées et alignées et, lorsque la deuxième résolution spatiale est inférieure à la première résolution spatiale, mettre à l’échelle l’image de deuxième résolution spatiale formée pour obtenir l’image de première résolution spatiale initiale (x 0 ). Dispositif (5) selon l’une des revendications précédentes, caractérisé en ce que l’unité de mise à jour (29) est agencée pour mettre à jour la variable auxiliaire (z) par descente de gradient comme suit : avec : où : - z est la variable auxiliaire, - x est la variable d’estimation de l’image de première résolution spatiale, - y k est la k-ème image de deuxième résolution spatiale, - K est le nombre d’images de deuxième résolution spatiale, - w k est le poids de déparasitage associé à l’image de deuxième résolution spatiale y k , - Δt k est la durée d’exposition de l’image de deuxième résolution spatiale y k , - W k est un opérateur de distorsion caractérisé par le paramètre de distorsion (p k ) associé à l’image de deuxième résolution spatiale y k , - C est un masque binaire modélisant une matrice de filtres colorés, - s est le facteur de mise à l’échelle utilisé pour obtenir l’image de première résolution spatiale initiale (x 0 ), - D s est un opérateur de décimation caractérisé par le facteur s de mise à l’échelle, - B est un opérateur de floutage, - δ est un paramètre de la descente de gradient, et - η est un paramètre de la fonction d’énergie quadratique à deux variables. Dispositif (5) selon l’une des revendications précédentes, caractérisé en ce que le module d’optimisation (13) est agencé pour implémenter trois itérations successives. Dispositif (5) selon l’une des revendications précédentes, caractérisé en ce qu’au moins un réseau de neurones convolutif (15, 25, 31) présente une architecture de type U-net. Procédé de formation d’une image d’une première résolution spatiale ( ) à partir d’une pluralité de couples formés chacun d’une image d’une deuxième résolution spatiale (y k ) inférieure ou égale à la première résolution spatiale et d’une durée d’exposition (Δt k ), ledit procédé étant mise en œuvre par le dispositif (5) selon l’une des revendications précédentes et comprenant : - une phase d’initialisation (210) comprenant les opérations suivantes : recevoir (300) les images de deuxième résolution spatiale (y k ), estimer (310), pour chaque image de deuxième résolution spatiale (y k ), un paramètre de distorsion en appliquant un algorithme Lucas-Kanade, appliquer, au moins à une image de deuxième résolution spatiale de référence (y k0 ) parmi les images de deuxième résolution spatiale (y k ), un dématriçage et, lorsque la deuxième résolution spatiale est inférieure à la première résolution spatiale, une mise à l’échelle pour obtenir une image de première résolution spatiale initiale (x 0 ), et - une phase d’optimisation (220) comprenant les opérations suivantes : recevoir (400) les images de deuxième résolution spatiale (y k ), les paramètres de distorsion (p k ) associés et l’image de première résolution spatiale initiale (x 0 ), minimiser une fonction d’énergie quadratique paramétrée par une variable d’estimation de l’image de première résolution spatiale (x) et une variable auxiliaire (z), ladite variable d’estimation de l’image de première résolution spatiale (x) et ladite variable auxiliaire (z) étant initialisées sur la base de l’image de première résolution spatiale initiale (x 0 ), la minimisation étant réalisée par itérations successives comprenant chacune la mise à jour (420) de la variable auxiliaire (z) par une descente de gradient en fonction des images de deuxième résolution spatiale (y k ) et des paramètres de distorsion (p k ) associés et la mise à jour (430) de la variable d’estimation de l’image de première résolution spatiale (x) en fonction de la variable auxiliaire (z) mise à jour, ledit procédé étant caractérisé en ce que la phase d’initialisation (210) comprend les opérations suivantes : recevoir (300) en outre les durées d’exposition (Δt k ) respectives des images de deuxième résolution spatiale (y k ), aligner (320) par interpolation l’image de deuxième résolution spatiale de référence (y k 0 ) en fonction de chaque paramètre de distorsion (p k ), utiliser (330) un réseau de neurones convolutif (25) recevant, en entrée, chaque image de deuxième résolution spatiale (y k ) et l’image de deuxième résolution spatiale de référence alignée (W k y k0 ) en fonction du paramètre de distorsion (p k ) associé et générant, en sortie, un facteur de confiance (g ρ ), calculer (340), pour chaque image de deuxième résolution spatiale (y k ), un poids de déparasitage (w k ) en fonction de la durée d’exposition (Δt k ) correspondante et du facteur de confiance (g ρ ) associé, et en ce que la phase d’optimisation (220) comprend les opérations suivantes : recevoir (400) en outre les poids de déparasitage (w k ) respectivement associés aux images de deuxième résolution spatiale (y k ), mettre à jour (420) la variable auxiliaire (z) par descente de gradient en fonction en outre des poids de déparasitage (w k ). Procédé selon la revendication 12, caractérisé en ce qu’il comprend en outre une phase d’entraînement préalable (200) du dispositif (5) avec un ensemble de données d’apprentissage comprenant une image de première résolution spatiale (x (1) ) et une pluralité de couples formés chacun d’une image de deuxième résolution spatiale (y k (1) ) et d’une durée d’exposition (Δt k (1) ), chaque image de deuxième résolution spatiale (y k (1) ) étant générée à partir de l’image de première résolution spatiale (x (1) ) comme suit : avec : où : - x (1) est l’image de première résolution spatiale de l’ensemble de données d’apprentissage, - y k (1) est la k-ème image de deuxième résolution spatiale de l’ensemble de données d’apprentissage, - Δt k (1) est la durée d’exposition de l’image de deuxième résolution spatiale y k (1) , - W k (1) est un opérateur de distorsion associé à l’image de deuxième résolution spatiale y k (1) , - C (1) est un masque binaire modélisant une matrice de filtres colorés, - D s (1) est un opérateur de décimation caractérisé par un facteur s de mise à l’échelle, - B (1) est un opérateur de floutage, - ε k (1) est un bruit associé à l’image de deuxième résolution spatiale y k (1) . Procédé selon la revendication 13, caractérisé en ce que la phase d’entraînement préalable (200) est supervisée en utilisant plusieurs ensembles de données d’apprentissage comprenant chacun une image de première résolution spatiale (x (i) ) et une pluralité de couples formés chacun d’une image de deuxième résolution spatiale (y k (i) ) et d’une durée d’exposition (Δt k (i) ), ladite phase d’entraînement préalable (200) supervisée permettant de paramétrer ledit dispositif (5) de manière à minimiser la quantité suivante : où : - est l’image de première résolution spatiale du i-ème ensemble de données d’apprentissage, - est l’image de première résolution spatiale estimée par le dispositif (5) à partir du i-ème ensemble de données d’apprentissage, et - n est le nombre d’ensembles de données d’apprentissage. Programme informatique comprenant des instructions dont l’exécution par au moins un processeur (9) se traduit par la mise en œuvre du procédé selon l’une des revendications 12 à 14.