Brevet FR2615979A1 — publié le 19881202

PROCEDE DE CORRECTION DE DISTORSION GEOMETRIQUE D'UNE IMAGE, NOTAMMENT UNE IMAGE PRODUITE PAR UNE GAMMA CAMERA La présente invention a pour objet un procédé de correction de distorsion géométrique d'une image, notamment d'une image produite par une gammacaméra. Elle trouve son application particulièrement dans le domaine médical o les gamma-caméras sont utilisées comme aide au diagnostic. Elle se rapporte aux caméras à scintillation Cou gamma-caméra) de type ANGER dont le brevet américain 3.011.057 décrit le fonctionnement dans ses principes et ses moyens de réalisation. Ces gamma-caméras sont destinées à détecter et à visualiser des photons émis par des corps radioactifs. Les gamma-caméras sont utilisées en médecine nucléaire pour visualiser dans un organe la répartition de molécules marquées par un'isotope radioactif que l'on a injecté à un patient. Une gamma-caméra comprend généralement un collimateur pour focaliser les photons gamma émis par le 2C patient, un cristal scintillateur pour transformer les photons gamma en photons lumineux -ou scintillations, et un réseau de tubes photo- multiplicateurs qui transforment chacune des scintillations en des impulsions électriques dites 2 5 contributions électriques de tubes. Elles comportent en outre des circuits électroniques pour produire à partir des contributions électriques fournies par les tubes photomultiplicateurs, des signaux de coordonnées X et Y du lieu o s'est produite la scintillation, ainsi qu'un signal de validation Z quand l'énergie W de la scintillation appartient à une bande d'énergie prédéterminées. Cette chaine de détection est généralement suivie d'un ensemble de visualisation comportant généralement un oscilloscope cathodique commandé par les signaux de coordonnées X, Y, et par Z pour visualiser par un point lumineux sur l'écran le 1c point d'impact du photon gamma sur le cristal. Cet impact est aussi appelé événement d'image. L'ensemble de visualisation peut éventuellement comporter un dispositif photographique pour Former une image de l'organe observé en intégrant un grand nombre de points lumineux produits sur l'écran cathodique. Il peut par ailleurs comprendre un dispositif de traitement numérique des images. En particulier l'ensemble de visualisation-peut être adapté à la présentation 2C de tomographies de l'organe observé. Pour atteindre ce but on acquiert plusieurs images de cet organe selon une pluralitéd'orientations d'observation de la gamma-caméra par rapport à cet organe. Par des traitements du signal, analogues à ceux rencontrés dans les tomodensitomêtres,.on peut reconstituer des images de coupes des organes examinés. Entre autres qualités une gamma-caméra doit posséder une bonne résolution spatiale, c'est à C dire la capacité de distinguer des petites sources radioactives rapprochées, une bonne réponse en taux de comptage, c'est à dire la capacité de traiter un grand nombre d'événements par unité de temps, et une qualité d'image indépendante de l'énergie de l'isotope considéré. La résolution spatiale dépend de la précision du calcul des coordonnées X et Y de chacun des évènements d'images. La qualité de l'élaboration de ces coordonnées dépend essentiellement des lois phgsiques régissant le Fonctionnement des différentes parties de la gamma- caméra. Ainsi l'interaction d'un photon gamma avec le cristal donne naissance à une scintillation lumineuse dont l'intensité décroit exponentiellement avec le temps. La constante de temps de cette décroissance est caractéristique du cristal scintillateur utilisé. Par exemple, pour un cristal d'iodure de sodium activé au thallium, NaI CT1), elle est de l'ordre de 250 nanosecondes. Cette scintillation est vue par plusieurs tubes photomultiplica- teurs simultanément. Les photons lumineux composant cette scintillation arrachent des photoélectrons aux photocathodes des tubes photomultiplicateurs. Le nombre de photoElectrons arrachés obéit, pour une scintillation donnée, à la loi statistique de Poisson. Ceci signiFie que la contribution électrique d'un tube photomultiplicateur recevant une scintillation a une amplitude dont la valeur suit une distribution statistique de Poisson et dont la valeur est Fonction de l'énergie des photons lumineux incidents. Par ailleurs, à énergie constante, l cette contribution électrique est une Fonction sensiblement gaussienne, de la distance séparant le centre de ce tube photomultiplicateur du lieu o s'est produit la scintillation. Si la scintillation se produit à l'aplomb du centre de ce tube, la contribution électrique est maximum. Plus le lieu de la scintillation est éloigné du 355 centre du tube, plus la contribution électrique 1 est faible. A titre d'exemple, si' une scintillation se produit à l'aplomb d'une paroi du tube la contribution électrique de celui-ci est environ réduite de moitié par rapport à la contribution électrique maximum. Une scintillation est vue par plusieurs tubes photomultiplicateurs simultanément, génàralement six à dix tubes. Aussi, la détermination de l'emplacement de cette scintillation sur le cristal, elle-même représentative du lieu d'émission du photon gamma d'excitation Cet donc de l'évènement d'image), peut être obtenue en calculant l'emplacement du barUcentre des contributions électriques délivres par l'ensemble des tubes photomultiplicateurs excités par cette scintillation. Ce calcul s'effectue simplement, selon ANGER, en injectant les contributions électriques au travers d'un jeu de matrices de résistances dont les valeurs des résistances sont fonction des positions des tubes photomulti- plicateurs auxquels elles sont raccordées. Les positions de ces tubes sont repérées par rapport à des axes cartésiens de référence dont le point d'intersection est généralement situé au centre du réseau de tubes. Dans chaque matrice, il g a autant de résistances qu'il W a de tubes photomultiplicateurs dans le réseau de tubes. Chacune des résistances est relive d'une part à la sortie d'un tube photomultiplicateur différent et d'autre part à un point commun qui constitue la sortie de la matrice. Ces résistances effectuent ainsi une pondération des contributions électriques de chacun des tubes photomulti- plicateurs qui les alimentent. Un des problèmes présentés par les détecteurs de gamma-caméras du twpe ANBER est qu'ils présentent des distorsions géométriques de l'image liées à la structure de captation de la lumière: cristal scintillateur - tubes photomultipli- cateurs - matrices de bargcentration. Les progrès de la médecine nucléaire, et en particulier les améliorations Faites pour recueillir davantage d'informations et des informations meilleures à 1C partir des gamma-caméras, par exemple pour la détection de petites tumeurs, conduisent a une absence de linéarité inhérente à la conception et à la réalisation des caméras. Il en résulte une distorsion spatiale des images gênante. Ces distorsions peuvent être à l'origine de défauts d'uniformité de densité importantes. A titre d'exemple une contraction de O,Lmm du rayon de l'image d'une surface circulaire de lcm de raWon conduit à un défaut d'uniformité de B%. Ces 2C défauts d'uniformité peuvent être particulièrement gênants en tomographie o l'effet d'amplification des dèfauts d'uniformité dû aux procèdès de reconstruction peut être supérieur à un facteur 10. On a tenté, dans l'état de la technique, de corriger les effets de cette distorsion. Le principe général de la correction est le suivant. On effectue la mesure d'une image d'une mire règulière places entre une source radioactive 2C d'émission uniforme et la gamma-camèra. Si la gammacaméra Était parfaite, l'image représentèe devrait correspondre à la répartition régulière des trous ou des fentes de la mire. La mesure effectuée montre en fait des irrégularités provenant des effets de distorsion spatiale. On peut cependant utiliser dans l'image ainsi acquise la connaissance de la distorsion, mesurée en comparant l'image acquise -à la répartition théorique à laquelle on aurait dû aboutir, pour effectuer les corrections des images ultérieurement acquises avec la même gamma-caméra. Ainsi, dans le brevet américain 37Y53B5, la mire, le fantôme d'acquisition, comportent une série de trous de diamètre 3 millimètres espacés les uns des autres de 2q millimètres. Quatre acquisitions sont faites en décalant à chaque fois en X et Y, la mire de 12 millimètres. Pour chaque point de mesure une paire de coefficients de correction est calculée et est stockée. Des points intermédiaires sont déterminés par interpolation Jusqu'à 6B x 6B pour tout le champ image. En fonctionnement normal, l'image distordue est acquise et est stockée. En fin d'acquisition, un programme de correction ?C redistribue les événements d'image de la matrice d'image distordue dans une matrice d'image corrigée obtenue grâce aux facteurs de correction. Cette méthode présente deux inconvénients. D'une part la redistribution des événements d'image en 2 fonction de la surface recouverte conduit à un élargissement -du contenu d'une maille sur quatre mailles. I1 en résulte automatiquement une perte de la résolution spatiale de l'image corrigée: elle devient moins précise. D'autre part, à cause de la distorsion, le nombre de coups à redistribuer n'est pas exact. En effet, il n'W a aucune raison pour que la surface d'une maille élémentaire donnée dans l'image acquise et distordue soit équivalente à la È)5 surface d'une maille élémentaire dans l'image corrigée. Dans un autre état de la technique représenté par le brevet français 2.412.B5E, le fantôme d'acquisition est une mire de fentes de 3 millimètres espacées de 15 millimètres. Quatre acquisitions sont également faites en les décalant de 7,5 mm en X et en Y. Pour chaque coordonnée X et Y de l'image distordue, un programme calcule les nouvelles coordonnées U et U de l'image non distordue. Ce programme ne calcule pas de coefficient de correction mais directement les nouvelles coordonnées. Ces nouvelles valeurs U et U sont stockées dans une matrice 6. x Ei. En fonctionnement de correction, X et Y sont codés sur 12 bits. Les six bits de poids fort de X et de Y adressent la mémoire des U et U. Puis une interpolation linéaire est effectuée en utilisant les six bits de poids faible de X et Y. On utilise ainsi le décalage d'un évènement d'image à placer par rapport aux quatre angles de la maille 2CL élémentaire de l'image corrigée dans laquelle cet évènement d'image doit finalement se trouver. Compte tenu d'un champ image d'environ iOO millimètres de diamètre, une codification sur douze bits des coordonnées des évènements d'image 2devrait conduire normalement à une précision d'image de l'ordre de 0,1 millimètre. Or, l'expérience montre que l'on n'aboutit pas à cette précision de l'image corrigée. Le document 2.q12.85S préconise bien de remplacer l'interpolation linéaire par une interpolation non linéaire, par exemple d'un type cubique, mais ce document ajoute que cependant l'interpolation linéaire donne une précision suffisante pour calculer les coordonnées vraies. Ceci signifie en fait, et l'expérience le vérifie, que l'interpolation non linéaire ne donne pas de meilleurs résultats. On le verra par la suite, la précision attendue de 0,1 millimètre dans l'image redressée, ne peut pas ainsi être atteinte. Un troisième état de la technique, constitué par le brevet européen 0021.366, préconise d'utiliser un fantôme d'acquisition, qui est une mire de Fentes, de 1 millimètre, espacées de 15 millimètres. Les données sont stockées dans une matrice de 256 x 256 sur 12 bits. Un calculateur élabore une matrice de correction de 6* x Bq de coefficients de correction. Là encore, les six bits de poids fort adressent la mémoire des coefficients de correction et permettent d'aboutir dans une maille élementaire de correction à laquelle sont en outre associés huit coefficients correspondant à une correction plus fine d'interpolation qui utilise les six bits de 2 poids faible des coordonnées calculées. Ultérieurement, une analyse d'uniformité est Faite en phase de calibration et de correction. Le calculateur calcule la densité Fi d'événements d'image par surface élémentaire Ai de l'image distordue acquise. Puis il calcule une densité F'i correspondant à une surface élémentaire de l'image corrigée de surface A'i obtenue après application des coefficients de correction. Si la correction de distorsion était parfaite F'i devrait être uniforme. Le brevet européen cité préconise d'utiliser le gradient de cette densité corrigée pour modifier itérativement les coefficients de correction de distorsion. Ce bçevet met dDnc en évidence ce qui est apparu à la lecture des deux précédents, à savoir qu'il faut absolument tenir compte de la densité d'informations en fonction de la distorsion. La solution préconisée dans ce dErnier document présente en outre un inconvénient important du fait qu'elle introduit des défauts de densité de l'image, inhérents à d'autres raisons que la distorsion d'image, dans la correction de distorsion d'image. Uns fois appliquée cette technique, il en résulte que la précision attendue ne peut pas être atteinte puisque le perfectionnement tend à répartir sur plusieurs mailles voisines des défauts qui apparaissent dans une maille, sans par ailleurs en discriminer les raisons. On s'est alors posé la question de savoir pourquoi ces différentes techniques qui, compte tenu du nombre de bits utilisé pour coder les coordonnées des évènements d'images, devaient normalement aboutir à une précision d'image attendue n'U aboutissaient pas. En particulier on s'est posé la question de savoir pourquoi, selon le deuxième état de la technique citée, la correction par approximation par polWnôme cubique n'était pas meilleure que l'interpolation linéaire. On a alors découvert que le choix d'un maillage ES x EU pour tout un champ image de diamètre OO millimètres, définissait en fait une maille trop grande, qui fournissait de ce fait une base d'interpolation présentant dès le départ un échantillonnage insuffisant de la fonction de distorsion et qu'on pouvait déduire du nombre de bits sur lequel étaient codées les coordonnées des événements d'image. D'autre part on s'est apperçu que les calculs intermédiaires d'interpolation devaient être faits avec une précision supérieure au 1/10' de millimètre si on voulait avoir en fin de calcul la précision souhaitee, c'est à dire le 1/10' de millimètre pour un défaut d'uniformité résiduel de l'ordre de 2. En définitive, on a déterminé qu'il W avait une relation entre le nombre minimum de mailles élémentaires dans le maillage utilisé pour établir la matrice de correction d'image et la précision attendue de la correction de distorsion spatiale; 1C cette relation est fonction du gradient de distorsion spatiale. On s'est rendu compte ainsi qu'il était nécessaire de faire un calcul d'interpolation plus précis que la précision minimale souhaitée. Aussi l'invention a pour objet un prDcédé de correction de distorsion spatiale de l'image produite par une gamma-caméra, cette gamma-caméra fournissant des signaux de sortie qui représentent, dans un champ image, des coordonnées 2C d'événements d'image, le procédé comprenant les étapes suivantes: - on mesure les coordonnées des événements d'image obtenus par l'interposition d'une mire test qui correspond à une distribution théorique uniforme de ces événements d'images, - on calcule par approximation cubique les déplacements de coordonnées des points d'un maillage théorique qui correspond à l'image test, - on en déduit des corrections de distorsion à affecter à des coordonnées d'événements d'image, - et on les applique à ces coordonnées, caractérisé en ce que on détermine un nombre minimum de mailles élementaires dans le maillage correspondant au champ image, en fonction d'une précision attendue de la correction de distorsion 5patible. 2615'979 1 1 L'invention sera mieux comprise à la lecture de la description qui suit et à l'examen des Figures qui l'accompagnent. Celles-ci ne sont données qu'à titre indicatif et nullement limitatif de l'invention. Les figures montrent: - figure 1: une mire test utilisable pour mettre en oeuvre le procédé de l'invention; - figure 2: un profil des signaux détectés pour une ligne de l'image de la mire; - 1C - Figure 3: la comparaison d'un réseau de droites parallèles équidistantes correspondant à l'image théorique à obtenir, à l'image distordue effectivement obtenue; - Figure q: le calcul par approximation cubique 'des dépIacements de coordonnées des points d'un maillage de l'image test; - figure 5: un procédé particulier de calcul de la maille élémentaire selon l'invention; - Figure 6: la méthode d'interpolation 2C utilisée dans l'invention pour corriger une image La figure 1 représente une mire test à utiliser pour mettre en oeuvre l'invention. La mire 1 comporte des fentes telles que 2, dont la largeur est de l'ordre d'un millimètre et qui sont espacées les unes des autres par des bandes 3, d'environ douze millimètres de large. Compte tenu d'un champ image de diamètre MOO millimètres, il W a trente trois Fentes. La fente centrale est legêrement décalée par rapport au centre de la mire de 3 millimètres de manière à ce que si elle est retournée de 180 degrés, on puisse définir un autre pas de fentes entrelacé avec le premier. On suppose que cette mire est placée entre une gamma-caméra non représentée (de champ image égal) et une source de rayonnement de raWonris gemma de 2615Z79 densité d'émission à peu prés uniforme dans le champ image. L'image de cette mire est analysêe selon un certain nombre de lignes. Par exemple à l'ordonnée Yi d'une ligne, on relève, au bout d'une certaine durée, l'histogramme H des évènements d'images en fonction de l'abscisse X des points de la ligne. L'histograr-ne H représenté sur la figure 2 a l'allure d'une succession de crêtes approximativement espacées les unes des autres du pas 3 de la mire 1. Si on appelle m le maximum de ces crêtes d'histogramme, on remarque que pour les fréquences 0,35 m la largeur des crêtes est de l'ordre de quatre millimètres. Il devient en conséquence inutile de rechercher à faire l'image d'une mire de pas plus fin: on risquerait de voir s'interpénétrer les crêtes relatives à des fentes adjacentes. Pour chacune de ces crêtes on calcule l'abscisse Xc du barycentre de la crête. Compte 2C tenu du fait que l'on a utilisé trente trois Fentes on obtiendra ainsi trente trois valeurs de Xc. Après retournement de la mire, on peut disposer de trente trois autres valeurs de Xc en principe interlacêes avec les trente trois précédentes. Le calcul des abscisses des barycentres est d'un type connu. De préférence il est effectué par programme de traitement numérique dans un ordinateur. On calcule ainsi pour chacune des lignes d'ordonnée Yi de l'image test les BE abscisses des barycentres. En pratique on effectue ce calcul pour 512 lignes en utilisant les S bits de poids fort des ordonnées des événements d'image pris en compte. La figure 3 représente sous la forme de séries de points les alignements de ces ú615979 barycentres en fonction de leur ligne d'ordonnée Yi. Il est évident qu'ils correspondent à des fentes dont les axes centraux al a2 etc... sont équidistants les uns des autres. La détermination de la position des axes ai par rapport au barycentre Xc Yi est effectuée en tenant compte de cette équidistance d'une part et en appliquant une méthode des moindres carrés d'autre part. Selon cette méthode la somme des carrés des variations d'abscisse DX est minimisée pour un jeu d'axes ai parallèles et équidistants les uns des autres. Le programme qui effectue la détermination du réseau des droites ai effectue de plus une analgse ligne par ligne pour vérifier s'il ne manque pas de barycentre. S'il en manque, ils sont déterminés par interpolation linéaire entre les barycentres voisins correspondant à la même fente ai. Le programme effectue ensuite une interpolation horizontale par un ensemble de PC polynômes cubiques de manière à aboutir à 512 points à partir des S abscisses de barycentres Xc mesurées. La figure L représente schématiquement cette interpolation et montre l'effet de l'invention. Cette figure représente 2- les déplacements DX d'abscisse des barwcentres à l'endroit des abscisses calculées Xci de ces barycentres. Par un ensemble de polynômes cubiques on est capable de reconstituer la courbe ' sur laquelle on détermine les 512 valeurs DX correspondant maintenant à 512 valeurs de Xc. Par opposition à l'état de la technique cité, on a découvert dans l'invention, que pour atteindre la précision attendue (dans le cas particulier de douze bits et de 0,1 mm) il était nécèssaire, compte tenu de l'ordre de grandeur du gradient de distorsion généralement obtenu sur les caméras de ANGER, d'utiliser une maille élémentaire du maillage de correction qui soit de l'ordre du 1/12B du champ image. Dans l'état de la technique cité un maillage 6B X 6S était utilisé. Comme on va le voir par la suite, un maillage EU X 6E conduit à des d@Fauts résiduels importants dans les zones d'image à Fort gradient de distorsion. Ce maillage correspond intuitivement et approximativement au nombre de Fentes de la mire utilisee. L'approximation cubique qui était ainsi Faite dans l'état de la technique tendait à remplacer les mesures d'abscisse de barycentre, par exemple les BE de l'exemple évoqué, par soixante quatre abscisses de correction qui, pour certaines étaient relativement proches des valeurs mesurées puisque 66 est très peu différent de Es. Puis, pour mener à bien la totalité de la 2C correction, dans l'état de la technique cité, on utilisait dans une maille élémentaire ainsi définie, une interpolation linéaire. Mais ceci signifiait que les abscisses corrigées, obtenues par interpolation linéaire, étaient entachées dés l'origine d'une erreur qui dans certains cas était nulle Cpar' exemple à l'endroit du point 5) et qui dans d'autres cas était très importante Cpar exemple l'erreur 6 à l'endroit de l'abscisse 7). Cette erreur pouvant être C supérieure à la précision attendue, le calcul d'interpolation ultérieurement effectué était illusoire car il portait sur de mauvaises bases. Par contre, dans l'invention, on choisit une maille élémentaire plus fine, 126 X 128, de telle 7 Fagon que la courbe q soit approchée précisément, et de telle façon que des interpolations linéaires sur une maille plus fine conduisent à des erreurs structurelles 8 inférieures à la précision attendue. Dans le cas présent on a découvert qu'il était nécessaire de choisir un maillage égal, ou plus fin que 12E X 12e8. Pour effectuer l'approximation cubique les équations aux limites sont telles que le polgnôme cubique passe par S, 10 et 11, et que la 1 drivée 12 du polUnôme en 10 soit la bissectrice des segments S - 10 et 10 - 11. Le mode de calcul des corrections d'abscisse par un polgnôme cubique donne des valeurs de DX en fonction de Yi, DX cYi), qui se suivent assez continuement en fonction de l'abscisse X. Par contre, d'une ligne Yi à une autre ligne immédiatement adjacente les valeurs de DX à une même abscisse sont relativement discontinues. Le programme peut alors eFFectuer un lissage vertical à coefficient 2C variable en fonction de la variation des pentes dans le sens vertical. Sur les bords de l'image, les Fentes sont extrapolées par des paraboles. De plus, les fentes trâs raccourcies aux-bords du champ, à droite ou à gauche, sont extrapolées par des droites. Il en résulte une nette amélioration de la correction de distorsion sur les bords de l'image. Ceci n'est pas trs' utile en imagerie directe C les praticiens s'arrangent toujours pour placer l'organe à visualiser au milieu du champ 3C image). Celà est cependant très utile en tomographie d'organes étendus o l'on réduit de ce fait les artefacts en anneaux. Rendu à ce stade, seulement la moitié du travail a été accompli. Il convient ensuite de tourner la mire 1 de 90' Cpuis encore une fois de 1BO0) et de recommencer toute l'opération-pour calculer les corrections d'ordonnée DYCXi) en fonction de Y. Ce travail est accompli comme précédemment pour X. On obtient finalement deux matrices de coefficients de correction de 512 X 512. Compte tenu de ce qui précède on revient à un maillage de 122 X 122 en faisant la mowenne, comme indiquée sur la figure 5, des corrections d'abscisse à appliquer à chacun des seize éléments 1C d'images adjacents du maillage 512 X 512 et appartenant à une même maille 13 du maillage 122 X 12E. Au moment o une image est acquise on peut décider de ne faire la correction qu'avec la matrice des 512 coefficients de correction, sans interpolation-supplémentaire. L'imaged'uniformité obtenue présente un moire qui est dû au recouvrement ou à l'espace entre des pixels voisins car il n'W a pas d'interpolation. Il est 2C donc nàcàssaire d'effectuer une interpolation pour éviter ce moiré. Cette interpolation a pour principe le repérage d'une maille élémentaire de l'image distordue dans laquelle se trouve un évènement d'image à replacer dans l'image corrigge. Soit en A B C D, figure 6, les quatre points situés dans les coins de la maille élementaire distordue en question. Les coefficients de correction qui permettent de définir la maille élementaire corrigée sont 3C matérialisés sur la figure 6 par les segments en tirets 1q à 17. Si on appelle DX et DY les deux coefficients de correction appliqués à chacun des points A à D et relatifs respectivement aux segments 14 à 17, on peut calculer, connaissant les coordonnées x W de l'événement d'image P dans la maille élémentaire A B C D, les corrections de distorsion d'image, DX et DY, à appliquer aux coordonnées X Y de ce point P dans le champ image. Les formules de correction sont les suivantes: DX-Cl-x-g+xg)DXa+Cx-xg)DXb+C(-xg)DXc+Cxy) DXd DY-(1-x-g+xg)DYa+Cx-xg)gYb+Cg-xg)DYc+Cxg)DYd o a b c et d se rapportent aux points A B C B. On a pu montrer que l'obtention de résultats DX et DY sur 12 bits pouvait être atteinte en calculant les coefficients DX et DY de correction sur onze bits et en calculant les coefficients d'interpolation tels que (1 - x - g + xg) sur B bits seulement. En simplifiant l'explication on peut dire que l'obtention à coup sûr du dixième de millimètre C 12 bits pour un champ image de diamètre 00 millimétres) nécessite d' ajouter au moins deux bits supplémentaires Cil faut être aussi précis que le quart de l'écart le plus fin à C montrer). Avec le bit de signe cela conduit à devoir manipuler 35 bits. Cependant, dans l'invention, utilement, on ne mesure que des déplacements de coordonnées et non pas des coordonnées corrigées. En estimant que le déplacement maximum rencontré serait de l'ordre de plus ou moins 2U millimètres (ce qui serait vraiment exceptionnel et nécessiterait en fait un rerèglage de la gamma-caméra) on a pu se rendre compte qu'on pouvait gagner ainsi un facteur C OO/2q - 16, ce qui est équivalent à quatre bits. Il en résulte qu'il n'est plus nêcêssaire que de manipuler 15 - M - 11 bits (un bit de signe, B bits de dynamique et 2 bits fractionnaires). Et en calculant les coefficients d'interpolation sur 8 bits seulement, à cause d'un raisonnement similaire, on peut obtenir des résultats DX et DY sur 12 bits. En utilisant pour effectuer une opération d'interpolation des sommateurs accumulateurs travaillant sur seize bits on peut facilement mettre en oeuvre la précision demandée. Cette dernière explication montre clairement que la précision est en fait obtenue sans augmenter la puissance de calcul mise en oeuvre par rapport à ce qui était décrit dans l'état de la technique. REVENDICATIONS 1 Procédé de correction de distorsion géométrique d'une image, notamment d'une image produite par une gamma-caméra, cette gamma-caméra fournissant des signaux (Xi, Yi) de sortie qui représentent, dans un champ image, des cDordonnées d'évènements d'images, le procéde comprenant les étapes suivantes: - on mesure les coordonnées des événements d'image CXc, Yi) obtenus par l'interposition d'une mire Cl) test qui correspond à une distribution théorique uniforme des évènements d'images, - on calcule (figure Y) par approximation cubique les déplacements (DX) de coordonnées des points d'un maillage théorique qui correspond à l'image test, - on en déduit des corrections de distorsion à affecter à des coordonnées d'évènements d'image - et on les applique à ces coordonnées, caractérisé en ce que - on détermine un nombre minimum de mailles élémentaires dans le maillage correspondant au champ image en fonction d'une précision attendue de la correction de distorsion spatiale. 2 Procèd, selon la revendication 1 caractérisé en ce que pour déduire on effectue une interpolation Cfigure 6) linéaire des déplacements CDXa, DXb) relatifs aux coordonnées des points CA,B,C,D) du maillage qui correspondent à un voisinage direct du lieu d'un évènement V d'image. 3 Procédé selon la revendication 1 ou 2 caractérisé en ce que pour calculer les déplacements par approximation cubique on les calcule pour une maille de base plus fine (512) que la maille élémentaire correspondant au dit nombre minimum (128) de mailles élémentaires et en ce qu'on effectue une mowenne des déplacements calculés, ceci afin d'augmenter la précision (fig S). Li Procèdè selon la revendication 3 caractérisé en ce que le champ image mesure environ 00 millimètres, en ce que la précision 1C attendue est de 0,1 mm, et en ce que le nombre minimum de mailles élémentaires dans le maillage vaut 128 X 128. Procèdè selon la revendication i caractérisé en ce que les déplacements relatifs (CDXa) sont exprimés sur au moins 11 bits, de manière à déduire des corrections de distorsion sur 12 bits à affecter à des coordonnées d'événements d'image. 6 Procèdé selon l'une quelconque des 2. revendications 3 à 5 caractérisé en ce que les expressions des déplacements comportent 2 bits fractionnaires pour permettre d'exprimer le 1/i du 1/10 de millimètre de précision pour obtenir une précision finale voulue du 1/10 de millimètre.