Brevet FR3071344A1 — publié le 20190322

L'invention concerne un procédé de détermination d'une fonction d'étalement de point d'un système d'imagerie, un procédé d'amélioration d'image acquise par un système d'imagerie et un dispositif mettant en oeuvre au moins un desdits procédés. Contexte de l'invention Une fonction d'étalement du point (« Point Spread Function (PSF) » en terminologie anglo-saxonne), appelée également réponse impulsionnelle optique (RIO), est une fonction mathématique décrivant une réponse d'un système d'imagerie à une source ponctuelle. La PSF modélise comment un point est imagé par un système d'imagerie et reflète des déformations optiques du système d'imagerie. On considère qu'une image d'un objet issue d'un système optique est le fruit d'une convolution entre un signal issu de l'objet (i.e. une image issue d'un système d'imagerie parfait sans déformation optique) et la PSF modélisant le système d'imagerie. Lorsque cette fonction est connue, il est possible de déconvoluer (i.e. appliquer une convolution inverse) l'image issue du système d'imagerie par la PSF pour obtenir une image améliorée de l'objet. Il existe diverses méthodes pour caractériser un système d'imagerie. On connaît par exemple la méthode dite des bords obliques (« Slanted Edge » en terminologie anglo-saxonne) comprise dans la norme ISO 12333 :2014, Photography -- Electronic still picture imaging -- Resolution and spatial frequency responses. Cette méthode permet de déterminer une fonction de transfert de modulation (« Modulation Transfert Function (MTF) » en terminologie anglo-saxonne) permettant d'évaluer des capacités d'un système d'imagerie à restituer fidèlement un signal. En optique, la MTF est la fonction qui définit un contraste en fonction d'une finesse des détails à restituer, i.e. en fonction de fréquences spatiales. La MTF correspond à un module d'une fonction de transfert optique (« Optical Transfert Function (OTF) » en terminologie anglo-saxonne) qui correspond elle-même à une transformée de Fourier de la PSF. La méthode décrite dans le document US20160080737A1 permet de déterminer une PSF d'un système de lentilles. Cette méthode utilise une mire, sous forme de grille 2 de test, pour déterminer la PSF. Cette grille de test est placée en face du système de lentilles, à une position correspondant à une distance focale dudit système, cette distance focale dépendant de la taille de la grille de test. Cette méthode ne permet donc que de déterminer la PSF correspondant à la position de la grille de test. La PSF obtenue 5 n'est donc pas valide pour toute autre position de la grille de test. Une image d'un objet situé à une autre position que la position pour laquelle a été calculée la PSF ne pourrait pas être corrigée correctement par la PSF. Cette méthode implique de calculer une autre PSF pour cette autre position et d'avoir des moyens pour déterminer la position de l'objet par rapport au système d'imagerie. 10 Il est souhaitable de pallier ces inconvénients de l'état de la technique. Il est notamment souhaitable de proposer un procédé et un dispositif qui permettent de déterminer une PSF « générique », valide pour plusieurs positions d'un objet par rapport au système d'imagerie. Par ailleurs, ce procédé doit être simple à mettre en oeuvre. 15 EXPOSE DE L'INVENTION Selon un premier aspect de l'invention, l'invention concerne un procédé de détermination d'une fonction d'étalement de point d'un système d'imagerie, comprenant : pour chaque position d'une pluralité de positions d'une mire dans un 20 champ optique dudit système d'imagerie : acquérir une image de la mire, dite image réelle ; obtenir une image synthétique de la mire représentant un modèle numérique de la mire recalé sur une zone de l'image réelle correspondant à la mire de manière à ce que ledit modèle coïncide à ladite zone ; estimer la fonction d'étalement de point en utilisant l'image réelle et l'image synthétique ; et, calculer une fonction de transfert 25 optique en appliquant une transformée de Fourier à la fonction d'étalement de point ; calculer une fonction de transfert optique moyenne à partir des fonctions de transfert optique calculées pour chaque position de la pluralité de positions ; et, obtenir une fonction d'étalement de point moyenne en appliquant une transformée de Fourier inverse à la fonction de transfert optique moyenne. 30 Ainsi, on obtient une fonction d'étalement de point valable pour une pluralité de positions d'un objet qui se situerait dans le champ optique du système d'imagerie. Il n'est alors pas nécessaire d'avoir une fonction d'étalement de point dédiée à chaque position possible d'un objet dans le champ optique du système d'imagerie, ni de déterminer la position dudit objet par rapport au système d'imagerie. 3 Selon un mode de réalisation, chaque position de la pluralité de positions se situe dans une zone de travail du système d'imagerie dans laquelle des objets devant être imagés par le système d'imagerie doivent se situer pour pouvoir être imagés par le système d'imagerie. 5 Ainsi, la fonction d'étalement de point est valide pour un objet situé à une position quelconque dans la zone de travail. Selon un mode de réalisation, l'estimation de la fonction d'étalement de point comprend : calculer une image de contour à partir de l'image réelle, dite image de contour réelle, et une image de contour à partir de l'image synthétique, dite image de 10 contour synthétique, une image de contour étant obtenue à partir d'une image originale et représentant exclusivement des contours d'objets contenus dans l'image originale ; déterminer la fonction d'étalement de point minimisant une fonction de coût comprenant un premier terme correspondant à une différence entre l'image de contour réelle et un produit de convolution entre la fonction d'étalement de point à estimer et 15 l'image de contour synthétique. Selon un mode de réalisation, la fonction de coût à minimiser s'écrit de la manière suivante : f (h) = 11h * JO -J/11 + yllhC + P11h111 où f (h) est la fonction de coût à minimiser, h est une matrice représentative de la 20 fonction d'étalement de point recherchée, Il. Il représente une norme L2au carré, Il. Ill représente une norme * est un opérateur de convolution, yIIhII est un premier terme de régularisation utilisé pour accélérer la convergence de la minimisation, le 11 h111 est un deuxième terme de régularisation utilisé pour obtenir une fonction d'étalement de point dérivable sur toute une image acquise par le système d'imagerie, y et le sont des 25 constantes prédéterminées. Selon un mode de réalisation, la mire comprend des objets comprenant des bords ayant de forts gradients spatiaux dans une pluralité de directions différentes. Selon un mode de réalisation, une transformation géométrique permettant de transformer le modèle numérique de la mire pour qu'il coïncide à ladite zone de l'image 30 réelle correspondant à la mire est obtenue pour chaque position de la pluralité de positions. Selon un mode de réalisation, la transformation géométrique est calculée en utilisant des points particuliers prédéfinis de la mire. 4 Selon un mode de réalisation, le procédé comprend une segmentation de la zone de l'image réelle correspondant à la mire, l'image réelle utilisée lors du recalage du modèle géométrique de la mire et lors de l'estimation de la fonction d'étalement de point comprenant exclusivement la zone comprenant la mire, chaque zone de l'image 5 réelle ne correspondant pas à la mire ayant été supprimée. Selon un deuxième aspect de l'invention, l'invention concerne un procédé d'amélioration d'image acquise en utilisant un système d'imagerie, comprenant : obtenir une fonction d'étalement de point du système d'imagerie en appliquant le procédé selon le premier aspect; puis, appliquer une convolution inverse entre la 10 fonction d'étalement de point ainsi obtenue et chaque image acquise en utilisant le système d'imagerie. Selon un troisième aspect de l'invention, l'invention concerne un dispositif de détermination d'une fonction d'étalement de point d'un système d'imagerie, comprenant : des moyens mis en oeuvre pour chaque position d'une pluralité de 15 positions d'une mire dans un champ optique dudit système d'imagerie comprenant : des moyens d'acquisition d'images pour acquérir une image de la mire, dite image réelle ; des moyens d'obtention pour obtenir une image synthétique de la mire représentant un modèle numérique de la mire recalé sur une zone de l'image réelle correspondant à la mire de manière à ce que ledit modèle coïncide à ladite zone ; des moyens d'estimation 20 pour estimer la fonction d'étalement de point en utilisant l'image réelle et l'image synthétique ; et, des moyens de calcul pour calculer une fonction de transfert optique en appliquant une transformée de Fourier à la fonction d'étalement de point ; et, des moyens de calcul pour calculer une fonction de transfert optique moyenne à partir des fonctions de transfert optique calculées pour chaque position de la pluralité de 25 positions ; et, des moyens d'obtention pour obtenir une fonction d'étalement de point moyenne en appliquant une transformée de Fourier inverse à la fonction de transfert optique moyenne. Selon un quatrième aspect de l'invention, l'invention concerne un programme d'ordinateur comprenant des instructions pour mettre en oeuvre, par un dispositif, le 30 procédé selon le premier ou le deuxième aspect, lorsque ledit programme est exécuté par un processeur dudit dispositif Selon un quatrième aspect de l'invention, l'invention concerne des moyens de stockage, stockant un programme d'ordinateur comprenant des instructions pour mettre 5 en oeuvre, par un dispositif, le procédé selon le premier ou le deuxième aspect, lorsque ledit programme est exécuté par un processeur dudit dispositif BREVE DESCRIPTION DES DESSINS 5 Les caractéristiques de l'invention mentionnées ci-dessus, ainsi que d'autres, apparaîtront plus clairement à la lecture de la description suivante d'un exemple de réalisation, ladite description étant faite en relation avec les dessins joints, parmi lesquels : 10 - la Fig. 1 décrit schématiquement un système d'imagerie pour lequel on détermine une PSF en utilisant le procédé selon l'invention ; - la Fig. 2 illustre schématiquement un exemple de mire adaptée au procédé selon 1 ' invention; - la Fig. 3 illustre schématiquement un exemple d'architecture matérielle d'un 15 module de traitement mettant en oeuvre le procédé de détermination d'une PSF selon l'invention ; - la Fig. 4 illustre schématiquement un exemple de procédé de détermination d'une PSF selon l'invention ; et, - la Fig. 5 illustre schématiquement un exemple de procédé d'amélioration 20 d'images comprenant le procédé de détermination d'une PSF selon l'invention. DESCRIPTION DETAILLEE DE DIVERS MODES DE REALISATION La description qui suit détaille plus particulièrement un mode de réalisation de la 25 présente invention dans un contexte où le système d'imagerie comprend un système d'acquisition d'images tel qu'une caméra numérique comprenant un objectif comprenant un réseau de lentilles, un système de mise au point (tel qu'un système d'autofocus) et un capteur d'images tel qu'un capteur d'images CCD (dispositif de transfert de charge, « Charged Couple device » en terminologie anglo-saxonne) ou 30 CMOS (semiconducteur metal-oxyde complementaire, « complementary metal-oxide semiconductor » en terminologie anglo-saxonne). Le procédé est toutefois adapté à d'autres types de systèmes d'imagerie produisant des images pouvant être traitées numériquement, tels que des caméras analogiques, des objectifs, etc. 6 La Fig. 1 décrit schématiquement un système d'imagerie 1 pour lequel on détermine une PSF en utilisant le procédé selon l'invention. Le système d'imagerie 1 de la Fig. 1 comprend une caméra numérique 10. Cette caméra 10 est connectée à un dispositif de traitement 14 mettant en oeuvre le procédé 5 selon l'invention. Devant la caméra 10 se trouve un objet 12, ici un objet plan, comprenant une mire. Dans un mode de réalisation, l'objet 12 est placé de sorte que le centre de la mire soit sur un axe optique de la caméra 10 à une distance Z. La distance Zi est comprise entre une distance minimum 13B et une distance maximum 13A, la zone entre les distances 13B et 13A formant une zone de travail. La zone de travail est 10 définie comme une zone dans laquelle des objets devant être imagés par la caméra 10 doivent se situer pour pouvoir être imagés par la caméra 10. Cette zone de travail pourrait être déterminée par exemple en utilisant un critère permettant de restituer fidèlement certains détails d'un objet. Une mesure d'un degré de flou dans une image pourrait constituer un tel critère. 15 Le système d'imagerie 1 comprend en outre, un système d'éclairage 11 éclairant l'objet 12. Dans un mode de réalisation, la caméra 10 comprend un capteur d'images comprenant des pixels de «2,5 !am », un objectif de focale « 25 mm », et une ouverture fixe de «2.8 », le système d'éclairage produit une lumière infrarouge de longueur 20 d'onde « 860 nm ». La zone de travail est comprise entre « 55 cm » et « 95 cm », une distance optimale de mise au point étant de « 75 cm ». Dans un mode de réalisation, l'objet 12 pourrait être placé de sorte que le centre de la mire soit dans le champ optique de la caméra 10, i.e. à un endroit visible par la caméra 10, mais sans être sur l'axe optique de la caméra 10. 25 La Fig. 2 illustre schématiquement un exemple de mire 120 adaptée au procédé selon 1 ' invention. Dans un mode de réalisation, la mire 120 comprend des objets (des formes géométriques) ayant de forts gradients spatiaux dans une pluralité de directions différentes. Dans l'exemple de la Fig. 2, la mire 120 comprend des portions d'un disque 30 tronquées possédant des bords tranchants disposés radialement. Dans un mode de réalisation, la mire 120 comprend des points particuliers permettant d'estimer une transformation géométrique de la mire due à des déformations optiques du système d'imagerie 1 (imperfection de l'objectif, imperfection du capteur 7 d'images, ...). Dans la mire 120 de la Fig. 2, les points particuliers sont quatre paires de carrés reliés par un de leurs coins. Dans un mode de réalisation, lorsque le système d'imagerie 1 induit des déformations optiques complexes, il est possible d'utiliser une mire comprenant plus de 5 quatre points particuliers, les points particuliers pouvant être des paires de carrés reliés par un de leurs coins ou toute autre forme géométrique adaptée telle que des disques, des croix, des triangles, etc. La Fig. 3 illustre schématiquement un exemple d'architecture matérielle du module de traitement 14. 10 Selon l'exemple d'architecture matérielle représenté à la Fig. 3, le module de traitement 14 comprend alors, reliés par un bus de communication 140 : un processeur ou CPU (« Central Processing Unit » en anglais) 141 ; une mémoire vive RAM (« Random Access Memory » en anglais) 142 ; une mémoire morte ROM (« Read Only Memory » en anglais) 143 ; une unité de stockage telle qu'un disque dur ou un lecteur 15 de support de stockage, tel qu'un lecteur de cartes SD (« Secure Digital » en anglais) 144 ; au moins une interface de communication 145 permettant au module de traitement 14 de communiquer avec la caméra 10 pour recevoir les images acquises par ladite caméra 10 et de contrôler des acquisitions d'images par la caméra 10. Le processeur 141 est capable d'exécuter des instructions chargées dans la RAM 20 142 à partir de la ROM 143, d'une mémoire externe (non représentée), d'un support de stockage (tel qu'une carte SD), ou d'un réseau de communication. Lorsque le module d'analyse 14 est mis sous tension, le processeur 141 est capable de lire de la RAM 142 des instructions et de les exécuter. Ces instructions forment un programme d'ordinateur causant la mise en oeuvre, par le processeur 141, du procédé décrit en relation avec la 25 Fig. 4 ou la Fig. 5. Les procédés décrits en relation avec les Figs. 4 ou 5 peuvent être implémentés sous forme logicielle par exécution d'un ensemble d'instructions par une machine programmable, par exemple un DSP (« Digital Signal Processor » en anglais), un microcontrôleur ou un GPU (processeur graphique, « Graphics Processing Unit » en 30 terminologie anglo-saxonne), ou être implémentés sous forme matérielle par une machine ou un composant dédié, par exemple un FPGA (« Field-Programmable Gate Array » en anglais) ou un ASIC (« Application-Specific Integrated Circuit » en anglais). On note que le module de traitement 14 aurait tout aussi bien pu être compris dans la caméra 10. 8 La Fig. 4 illustre schématiquement un exemple de procédé de détermination d'une PSF selon l'invention. Un objectif du procédé de la Fig. 4 est de déterminer une PSF du système d'imagerie 1 générique dans le sens où cette PSF doit être adaptée à une pluralité de 5 positions d'un objet dans la zone de travail et non à une seule position particulière. Le procédé de la Fig. 4 débute par une étape 41 au cours de laquelle le module de traitement 14 initialise une variable i à zéro. La variable i est un nombre entier positif utilisé par la suite pour parcourir un ensemble prédéfini de distances Zi comprises dans la zone de travail. L'ensemble prédéfini comprend un nombre max de distances Z. 10 Dans une étape 42, le module de traitement 14 initialise une variable z, représentant une distance, à une valeur d'une des distances Zi correspondant à la valeur courante de la variable i. Dans une étape 43, le module de traitement provoque le positionnement de l'objet 12 à la distance z du système d'imagerie 1, c'est-à-dire à une distance z par rapport au 15 bord extérieur de l'objectif de la caméra 10 le plus proche de la mire. Dans un mode de réalisation, l'objet 12 est fixé sur un dispositif coulissant positionné sur un rail horizontal parallèle à l'axe optique de la caméra 10. Le module de traitement 14 contrôle alors le dispositif coulissant pour positionner l'objet 12 à la distance souhaitée de la caméra 10. 20 Au cours de l'étape 43, le module de traitement 14 provoque une acquisition d'une image, dite image réelle, de l'objet 12 par le système d'imagerie 1. Dans une étape 45, le module de traitement 14 obtient une image synthétique de la mire 120 représentant un modèle numérique de la mire 120 recalé sur une zone de l'image réelle correspondant à la mire de manière à ce que ledit modèle numérique 25 coïncide (i.e. se superpose) à ladite zone de l'image réelle. Le modèle numérique de la mire 120, correspond à une image idéale de la mire 120 sans déformation optique et sans convolution par la PSF du système d'imagerie 1. Dans un mode de réalisation, l'objet 12 est positionné perpendiculairement à l'axe optique de telle sorte que le centre de la mire soit sur l'axe optique de la caméra 10. 30 Dans ce cas, le recalage du modèle numérique se fait de manière itérative. On suppose que le module de traitement 14 stocke le modèle numérique de la mire 120 dans une taille prédéfinie correspondant à une taille maximum que pourrait avoir la mire 120 dans l'image réelle. Dans une première itération, le module de traitement utilise le modèle numérique de la mire 120 dans la taille prédéfinie et calcule une différence entre le 9 modèle numérique et la mire 120 dans l'image réelle. Dans chaque itération suivante, le module de traitement 14 réduit la taille de la mire et calcule à nouveau la différence entre le modèle numérique et la mire 120 dans l'image réelle. Les itérations s'arrêtent lorsque le module de traitement 14 a trouvé une différence minimum. Le modèle 5 numérique de la mire 120 dans la taille minimisant la différence est alors considéré comme recalé sur le mire 120 de l'image réelle. Le module de traitement 120 génère alors l'image synthétique à partir du modèle numérique de la mire 120 ainsi recalé. Dans un mode de réalisation, le module de traitement 14 détermine une transformation géométrique permettant de passer du modèle numérique de la mire 120 10 dans la taille prédéfinie, à la mire 120 telle qu'elle apparaît dans l'image réelle. Cette détermination est basée sur une détection des points particuliers dans l'image réelle et une détermination de la transformation géométrique permettant de projeter les points particuliers du modèle numérique sur les points particuliers de l'image réelle. La transformation est ensuite appliquée au modèle numérique de la mire 120 pour obtenir 15 l'image synthétique. La détermination de la transformation géométrique peut utiliser par exemple une méthode décrite dans les publications Goshtasby, Ardeshir, "Piecewise linear mapping functions for image registration," Pattern Recognition, Vol. 19, 1986, pp. 459-466 et Goshtasby, Ardeshir, "Image registration by local approximation methods," Image and Vision Computing, Vol. 6, 1988, pp. 255-261. 20 Dans une mise en oeuvre particulière, l'image réelle est affichée sur un dispositif d'affichage d'un utilisateur et l'utilisateur détoure une position de chaque point particulier. La position de chaque point particulier ainsi détourée est ensuite fournie au module de traitement 14 afin qu'il mette en oeuvre une procédure d'estimation de la transformation géométrique. 25 Dans un mode de réalisation, le module de traitement 14 obtient pour chaque distance Zi, une transformation géométrique prédéterminée. Chaque transformation géométrique a été déterminée en plaçant un objet à la distance Zi du système d'imagerie 1, en acquérant une image dudit objet, puis en déterminant la transformation géométrique permettant de passer d'une image de référence dudit objet à l'image 30 acquise par le système d'imagerie 1. Dans une étape 46, le module de traitement 14 estime la PSF du système d'imagerie 1 en utilisant l'image réelle et l'image synthétique. Dans un mode de réalisation, le module de traitement 14 calcule une image de contour, dite image de contour réelle, à partir de l'image réelle J/ et une image de 10 contour, dite image de contour synthétique, à partir de l'image synthétique JO. Une image de contour est obtenue à partir d'une image originale et représente exclusivement des contours d'objets contenus dans l'image originale. On peut obtenir une image de contour par diverses méthodes telles que par exemple, un filtrage de Canny ou un 5 filtrage de Sobel. Dans une mise en oeuvre particulière, chaque image de contour est obtenue en appliquant un opérateur gradient à chaque image originale respective (i.e. l'image réelle et l'image synthétique). Le module de traitement 14 détermine ensuite la PSF minimisant une fonction de coût f (h) comprenant un premier terme correspondant à une différence entre l'image de contour réelle J/ et un produit de convolution entre la 10 PSF du système d'imagerie 1 à estimer et l'image de contour synthétique O. La fonction de coût f (h) peut s'écrire de la manière suivante : f (h) = 11h * JO -J/11 + yllhC + P11h111 où h est une matrice représentative de la PSF recherchée, Il. II représente une norme L'au carré, 11.11i représente une norme * est un opérateur de convolution. Le 15 terme y II hC est un premier terme de régularisation que le module de traitement 14 utilise optionnellement pour accélérer la convergence de la minimisation. Le terme P11h111 est un deuxième terme de régularisation que le module de traitement 14 utilise optionnellement pour obtenir une PSF dérivable sur toute l'image, c'est-à-dire, ayant une forme lissée. y et le sont des constantes prédéterminées. Dans un mode de 20 réalisation, la constante y est déterminée en effectuant plusieurs minimisations sur une base d'images en faisant varier la constante y à chaque image de la base. A la fin de chaque minimisation, une information représentative d'un niveau de netteté de l'image obtenue par déconvolution par la PSF est calculée et on garde la constante y permettant d'obtenir le meilleur niveau de netteté. La constante le est obtenue de manière similaire. 25 Il existe diverses méthodes de minimisation connues applicables à la minimisation de la fonction de coût f (h). Dans un mode de réalisation, la fonction de coût f (h) est minimisée par une méthode de gradient conjugué ou par une méthode des multiplicateurs de Lagrange. La PSF minimisant la fonction de coût f (h) est notée hi. 30 Dans une étape 47, le module de traitement 14 calcule l'OTF de la PSF hi, notée Hi, en appliquant une transformée de Fourier à la PSF hi. Dans une étape 48, le module de traitement 14 détermine si la variable i est égale au nombre de distances dans l'ensemble prédéfini de distances Z. Si ce n'est pas le cas, 11 la variable i est incrémentée d'une unité par le module de traitement 14 lors d'une étape 49 et le module de traitement 14 retourne à l'étape 42. Si la variable i a atteint le nombre de distances dans l'ensemble prédéfini de distances Zi, le module de traitement 14 passe à une étape 50, au cours de laquelle le 5 module de traitement 14 calcule une OTF moyenne, notée Hmoy, à partir des OTF Hi calculées pour chaque distance Z. Dans un mode de réalisation : vi REVENDICATIONS1) Procédé de détermination d'une fonction d'étalement de point d'un système d'imagerie (1), caractérisé en ce qu'il comprend : pour chaque position d'une pluralité de positions d'une mire (120) dans un champ optique dudit système d'imagerie : acquérir (43) une image de la mire, dite image réelle ; obtenir (45) une image synthétique de la mire représentant un modèle numérique de la mire recalé sur une zone de l'image réelle correspondant à la mire de manière à ce que ledit modèle coïncide à ladite zone ; estimer (46) la fonction d'étalement de point en utilisant l'image réelle et l'image synthétique ; et, calculer (47) une fonction de transfert optique en appliquant une transformée de Fourier à la fonction d'étalement de point ; calculer (50) une fonction de transfert optique moyenne à partir des fonctions de transfert optique calculées pour chaque position de la pluralité de positions ; et, obtenir (51) une fonction d'étalement de point moyenne en appliquant une transformée de Fourier inverse à la fonction de transfert optique moyenne. 2) Procédé selon la revendication 1, caractérisé en ce que chaque position de la pluralité de positions se situe dans une zone de travail du système d'imagerie dans laquelle des objets devant être imagés par le système d'imagerie doivent se situer pour pouvoir être imagés par le système d'imagerie. 3) Procédé selon la revendication 1 ou 2, caractérisé en ce que l'estimation de la fonction d'étalement de point comprend : calculer une image de contour à partir de l'image réelle, dite image de contour réelle, et une image de contour à partir de l'image synthétique, dite image de contour synthétique, une image de contour étant obtenue à partir d'une image originale et représentant exclusivement des contours d'objets contenus dans l'image originale ; déterminer la fonction d'étalement de point minimisant une fonction de coût comprenant un premier terme correspondant à une différence entre l'image de contour 14 réelle et un produit de convolution entre la fonction d'étalement de point à estimer et l'image de contour synthétique. 4) Procédé selon la revendication 3, caractérisé en ce que la fonction de coût à 5 minimiser s'écrit de la manière suivante : f (h) = 11h * JO -J/11 + yllhC + 11h111 où f (h) est la fonction de coût à minimiser, h est une matrice représentative de la fonction d'étalement de point recherchée, Il. Il représente une norme L2au carré, 11.111 représente une norme * est un opérateur de convolution, yIIhII est un premier terme 10 de régularisation utilisé pour accélérer la convergence de la minimisation, p11 h 111 est un deuxième terme de régularisation utilisé pour obtenir une fonction d'étalement de point dérivable sur toute une image acquise par le système d'imagerie, y et le sont des constantes prédéterminées. 15 5) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que la mire comprend des objets comprenant des bords ayant de forts gradients spatiaux dans une pluralité de directions différentes. 6) Procédé selon l'une quelconque des revendications précédentes, caractérisé en 20 ce qu'une transformation géométrique permettant de transformer le modèle numérique de la mire pour qu'il coïncide à ladite zone de l'image réelle correspondant à la mire est obtenue pour chaque position de la pluralité de positions. 7) Procédé selon la revendication 6, caractérisé en ce que la transformation 25 géométrique est calculée en utilisant des points particuliers prédéfinis de la mire. 8) Procédé selon l'une quelconque des revendications précédentes, caractérisé en ce que le procédé comprend une segmentation (44) de la zone de l'image réelle correspondant à la mire, l'image réelle utilisée lors du recalage (45) du modèle 30 géométrique de la mire et lors de l'estimation (46) de la fonction d'étalement de point comprenant exclusivement la zone comprenant la mire, chaque zone de l'image réelle ne correspondant pas à la mire ayant été supprimée. 15 9) Procédé d'amélioration d'image acquise en utilisant un système d'imagerie, caractérisé en ce qu'il comprend : obtenir (51) une fonction d'étalement de point du système d'imagerie en appliquant le procédé selon l'une quelconque des revendications 1 à 8 ; puis, 5 appliquer (53) une convolution inverse entre la fonction d'étalement de point ainsi obtenue et chaque image acquise en utilisant le système d'imagerie. 10) Dispositif de détermination d'une fonction d'étalement de point d'un système d'imagerie (1), caractérisé en ce qu'il comprend : 10 des moyens mis en oeuvre pour chaque position d'une pluralité de positions d'une mire (120) dans un champ optique dudit système d'imagerie comprenant : des moyens d'acquisition d'images (10) pour acquérir une image de la mire, dite image réelle ; des moyens d'obtention (45) pour obtenir une image synthétique de la mire 15 représentant un modèle numérique de la mire recalé sur une zone de l'image réelle correspondant à la mire de manière à ce que ledit modèle coïncide à ladite zone ; des moyens d'estimation (46) pour estimer la fonction d'étalement de point en utilisant l'image réelle et l'image synthétique ; et, des moyens de calcul (47) pour calculer une fonction de transfert optique en 20 appliquant une transformée de Fourier à la fonction d'étalement de point ; et, des moyens de calcul (50) pour calculer une fonction de transfert optique moyenne à partir des fonctions de transfert optique calculées pour chaque position de la pluralité de positions ; et, des moyens d'obtention (51) pour obtenir une fonction d'étalement de point 25 moyenne en appliquant une transformée de Fourier inverse à la fonction de transfert optique moyenne. 11) Programme d'ordinateur, caractérisé en ce qu'il comprend des instructions pour mettre en oeuvre, par un dispositif (14), le procédé selon l'une quelconque des 30 revendications 1 à 9, lorsque ledit programme est exécuté par un processeur (141) dudit dispositif. 12) Moyens de stockage, caractérisés en ce qu'ils stockent un programme d'ordinateur comprenant des instructions pour mettre en oeuvre, par un dispositif (14), 16 le procédé selon l'une quelconque des revendications 1 à 9, lorsque ledit programme est exécuté par un processeur (141) dudit dispositif