Brevet FR2714736A1 — publié le 19950707

SYSTEME DE LOCALISATION D'EVENEXENTS DANS UNE CANERA A RAYONS GAMMA La présente invention se rapporte aux caméras à rayons gamma et plus particulièrement à des systèmes et à des procédés pour localiser des événements dans lesquels un événement est l'impact d'un rayonnement gamma sur un scintillateur en provoquant une scintillation. La plupart des caméras à rayons gamma utilisent actuellement ce que l'on connait sous le nom de mathématiques d'Anger ou calculs d'Anger, tels que les présente le brevet US numéro 3 011 057, pour déterminer la localisation d'un événement au niveau de ce que l'on peut appeler une première approximation. Ensuite, on utilise des facteurs de correction pour corriger la localisation déterminée par le calcul d'Anger. Les facteurs de correction comprennent, par exemple des facteurs de correction de décalage, obtenus en utilisant un champ submergeant en liaison avec un écran aux rayons gamma comportant des perforations dans des localisations connues. Les facteurs de correction sont calculés pour déterminer la distance entre les localisations connues des perforations et les localisations des images des perforations. Ces calculs sont utilisés pour corriger les localisations des événements tels qu'ils se produisent lors de l'utilisation clinique de la caméra à rayons gamma. En pratique il est raisonnable de recalibrer les facteurs de correction périodiquement, c'est-à-dire, d'effectuer périodiquement le traitement de calibrage du décalage par champ submergeant parce que les facteurs de correction de décalage obtenus antérieurement ne sont plus précis en raison de la dérive des composants à base de matériel. Par conséquent, les facteurs de correction qui sont supposés corrects au commencement de la période, deviennent éventuellement relativement incorrects à moins que les circuits électroniques ne soient "réaccordés". En plus d'être incommode, le calibrage de correction de décalage prend du temps et nécessite des équipements. En outre, la correction de décalage actuelle, bien que faite en cours d'utilisation, nécessite toujours une mémoire pour mémoriser les facteurs de correction et un module d'exécution de calculs mathématiques pour effectuer des calculs en utilisant les facteurs de correction. Par conséquent, le besoin de recalibrage nécessite davantage de matériel et/ou de logiciel dans les caméras à rayons gamma. Depuis l'adoption largement répandue des systèmes de calcul de caméra d'Anger, on a imaginé d'autres moyens et d'autres méthodes pour obtenir un procédé plus précis pour localiser de manière précise des événements à mesure qu'ils se produisent. On a suggéré des procédés statistiques pour localiser les événements sur la caméra à rayons gamma à des fins de formation d'image. Cependant, dans le passé ces procédés statistiques ont nécessité de vastes mémoires, des calculs intensifs et un équipement supplémentaire pour les calculs réels de sorte que les procédés statistiques se sont avérés peu pratiques. Sous son aspect le plus large, la présente invention propose une caméra à rayons gamma possédant un système nouveau et amélioré de localisation d'événements, incorporé à l'intérieur de ladite caméra à rayons gamma, le système de localisation d'événement comprenant: un cristal scintillateur qui scintille en réponse à l'impact, sur lui, d'un rayonnement; plusieurs tubes photomultiplicateurs (parfois appelés PM dans la suite) pour produire des signaux électriques en réponse aux scintillations (également appelées événements lumineux); ledit système comportant des PM avec des courbes de sortie qui sont approximativement symétriques radialement en réponse à des événements lumineux détectés par les PM; des circuits de calculateur mis en oeuvre en réponse auxdits signaux électriques pour déterminer la localisation non corrigée des événements détectés et l'énergie de l'événement détecté; lesdits circuits de calculateur pour déterminer la localisation non corrigée et l'énergie de l'événement comprenant en outre: 1) un estimateur pour estimer la localisation en X, Y de l'événement; 2) un calculateur de valeur attendue utilisant la localisation en X, Y estimée de l'événement par rapport aux localisations de chacun des PM pour déterminer de façon analytique une valeur attendue de la sortie d'énergie de chacun des PM en réponse audit événement, en se basant sur la localisation de chacun des PM par rapport à l'événement; 3) un définisseur de localisation en X, Y, mis en oeuvre en réponse à ladite valeur attendue déterminée, pour établir une fonction d'optimisation de localisation d'événement pour définir la localisation en X, Y de l'événement par rapport à la localisation de chacun des PM; 4) et une boucle d'itération utilisant la fonction d'optimisation de localisation d'événement pour déterminer de façon itérative les coordonnées corrigées Xc, Yc de l'événement. Les caractéristiques et avantages de l'invention ressortiront d'ailleurs de la description qui va suivre à titre d'exemple en se référant aux dessins annexés, dans lesquels: la figure 1 est une vue partielle en perspective d'une tête de caméra à rayons gamma en même temps qu'un schéma fonctionnel montrant les circuits de calcul et de traitement d'image qui ont pour résultat l'image affichée; la figure 2 est une vue plus détaillée montrant le circuit de calcul; et la figure 3 est un mode de réalisation d'une boucle d'itération utilisée pour la correction de localisation dans le bloc correcteur de la figure 1. Pour mettre en oeuvre le système on doit prévoir plusieurs tables de consultation fixes. On utilise les tables de consultation pour montrer des éléments mesurés comme les localisations en X et Y du centre de chacun des photomultiplicateurs; globalement désignés par Xi, Yi. Les tables de consultation contiennent également des éléments calculés comme la valeur attendue ai de la sortie d'énergie de chaque photomultiplicateur PMi et la différentielle du logarithme naturel de la valeur attendue indiquée par 8(u) qui est la différentielle du logarithme naturel de la valeur attendue a(u) par rapport à du. La figure 1 montre un système de caméra à rayons gamma de base. Le système de formation d'image en médecine nucléaire 11 de la figure 1 comprend une tête de caméra à rayons gamma 12 incluant un cristal à scintillation 13 et plusieurs tubes photomultiplicateurs (PM) comme le PM 14. Des circuits de calcul 16 sont prévus, qui déterminent la localisation en coordonnées non corrigées de chaque événement Xu, Yu. Un événement est défini comme une scintillation indiquée en 17, par exemple, provoquée par l'impact sur un cristal 13 d'un photon gamma se déplaçant suivant une trajectoire 18 représentée par le trait interrompu. Le circuit de calcul 16 détermine également l'énergie totale mesurée Zu de l'événement. Ainsi, le circuit de calcul délivre, au niveau de sa sortie 19, des valeurs non corrigées de Xu, Yu et Zu. Comme on le sait bien, les coordonnées de localisation Xu, Yu obtenues du circuit de calcul 16 sont influencées par la localisation de l'événement par rapport aux PM. Par conséquent, les coordonnées sont erronées comme l'est l'énergie mesurée Zu. Les caméras à rayons gamma modernes corrigent les coordonnées déterminées par le circuit de calcul d'une manière similaire à celle imaginée par le brevet US numéro 3 745 345. L'énergie mesurée est corrigée en fonction de la localisation de l'événement d'une manière imaginée en premier dans le brevet US numéro 4 095 108. Comme on l'a indiqué au commencement de la présente description, les corrections de localisation de la technique antérieure ont comme défaut, entre autres choses, qu'elles nécessitent un calibrage périodique. Le circuit de correction 21 élimine les défauts des circuits de correction de localisation de la technique antérieure. Le circuit de correction 21 détermine de façon analytique la valeur de l'énergie attendue de chacun des photomultiplicateurs en se basant sur la localisation de l'événement par rapport à chacun des PM. On établit une fonction d'optimisation de localisation qui définit les localisations en X, Y des événements. Les localisations en X, Y définies sont soumises à une boucle d'itération pour déterminer finalement les coordonnées corrigées Xc, Yc- Ces coordonnées corrigées sont utilisées par le processeur d'image 22 pour fournir l'affichage sur un écran 24. L'information mesurée acquise à partir de la caméra comprend le contenu d'énergie mesurée (Zi) de l'événement au droit de chaque PM. Plus précisément Zi est l'énergie à la sortie de chaque PM qui est une fonction de l'entrée d'énergie provenant de l'événement et de la localisation du photomultiplicateur par rapport à la localisation de l'événement. La caméra à rayons gamma mesure la sortie d'énergie Zi de chaque tube photomultiplicateur, comme l'indique le bloc 26 de la figure 2, qui indique la mesure Zi s'appliquant aux PM 1, 2,... n. Les mesures de Zi sont additionnées, c'est-à-dire que l'on obtient la sommation des sorties d'énergie des photomultiplicateurs, laquelle est l'énergie mesurée totale de l'événement. La proportion par tube photomultiplicateur de l'énergie mesurée par rapport à l'énergie mesurée totale est désignée par ni, comme l'indique le bloc 27. Le bloc 28 montre un exemple de calcul des coordonnées en X, Y des événements. Ce calcul est le calcul d'Anger connu. Le calcul détermine les coordonnées Xu et Yu non corrigées de chaque événement. Par exemple, Xu est la sommation de chaque Xi et ni, et Yu est la sommation de chaque Yi et ni, o Xi et Yi sont les coordonnées des localisations de chaque PM qui sont obtenues à partir des tables de consultation, c'est-à- dire de la table de consultation no 1. Ce qui vient d'être décrit jusqu'à maintenant se rapportant aux figures 1 et 2, à partir du bloc 12 jusqu'au bloc 28, pourrait être le calcul d'Anger normal effectué par toutes les caméras à rayons gamma du type Anger. La plupart des caméras à rayons gamma actuellement en utilisation sont des caméras du type Anger. La base mathématique réelle du système de correction de localisation décrit ici suppose une symétrie radiale des courbes de photomultiplicateur. Par définition (Xi, Yi) sont les coordonnées du centre du PM numéroté i, et ni est la sortie d'énergie normalisée du PM numéroté i. La courbe de PM symétrique radialement ai est définie comme la valeur attendue de l'énergie ni, soit ai = E (ni). En supposant que ai soit symétrique radialement, il s'ensuit que: ai (X, Y) = ai (ui) o: ui2 = (X-Xi)2 + (X-Yi)2 La localisation de l'événement est estimée en supposant que M(%l, 12,... a1, a2,...) soit une fonction d'optimisation (de maximisation) de localisation d'événement, c'est-à-dire que les conditions extrêmes: aM aM -0; -= 0 ax aY définissent la localisation la plus probable (X, Y) d'un événement. A partir de la définition de M, et de ai = a(ui) on obtient: aM aM axi aUi 0 = = - * ' ax aXi aui ax Puisque: aui = 2 (X-Xi) ax on peut réécrire l'équation de localisation d'événement comme: aM 0 = - =2 Z (X-Xi) * Wi; ax soit: aM 0 = = 2 Z (X-Xi) Wi; aY o: aM d Wi = (ai) (lnai) aai dui Finalement: X = zXiWi / zWi et d'une manière similaire: Y = EYiWi / zWi Les deux dernières équations indiquent le procédé itératif simple suivant, pour obtenir l'estimation la plus probable pour la localisation d'événement corrigée (Xc, Yc): (a) On prépare les tables de consultation pour a(u) et B(u), o: B(u) = _ (ln aU) au (b) Etant donné un événement qui produit la sortie de PM {Zi}; on calcule les coordonnées de localisation d'événement non corrigées Xu, Yu; par exemple, on utilise "l'arithmétique d'Anger": Xu = Z Xi qi; Yu = z Yi Ti; qi = Zi / Z Zi (c) On utilise les coordonnées de localisation non corrigées (Xu, Yu) pour calculer ui: ui2 = (Xu - Xi)2 + (Yu - Yi)2 (d) En utilisant les tables de consultation pour ai et Bi, on calcule: aM Wi = (di a) Bi aai (e) On calcule l'estimation suivante de (X, Y): X = Z XiWi / Z Wi; Y = Z YiWi / Z Wi (f) On répète les étapes (c) à (e) jusqu'à ce que l'on obtienne une convergence. Le procédé peut être facilement appliqué à la vraisemblance logarithmique maximale: Pour la vraisemblance logarithmique maximale: a i M = ln [n eai] ni! = Z [Ni ln ai - ai - ln (i!)] Ainsi: aM i aiai et par conséquent: d Wi = (Ni ai) (ln ai) du = (ni - ai) Bi La figure 3 représente en détail l'équipement nécessaire pour le procédé de localisation d'événement ci-dessus. On y voit que les valeurs Xu et Yu non corrigées reçues du circuit de calculateur 16 sont maintenant utilisées pour obtenir ui2; c'est-à-dire la localisation vectorielle de l'événement par rapport à chaque tube photomultiplicateur. Plus particulièrement, cette valeur ui2 est égale à (X-Xi)2 + (Y-Yi)2 Comme le montre la figure 3, le circuit de calcul 16 fournit les coordonnées d'événement non corrigées Xu, Yu et l'énergie normalisée, c'est-à-dire, Zi / Z Zi, soit ni. La table de consultation I indiquée au bloc 31 est utilisée pour fournir les valeurs connues des coordonnées Xi, Yi ce qui est la localisation de chacun des PM. La table de consultation II indiquée au bloc 32 fournit la valeur attendue a(u) de l'énergie ni pour chaque PM par rapport à la localisation non corrigée Xu, Yu de l'événement. La table III indiquée en 33 donne les valeurs de la différentielle du logarithme naturel de la courbe symétrique radialement ai pour chaque PM. Le bloc 34 calcule la distance (u) du centre du PMi à l'événement en utilisant les coordonnées non corrigées Xu, Yu de l'événement. Le bloc 36 indique que la valeur attendue a de l'énergie de chaque PM et Bi, le logarithme naturel des valeurs attendues de chaque PM, sont obtenues, respectivement, des tables II et III. En utilisant la valeur de ui2 obtenue du bloc 34, le bloc 37 permet d'obtenir pour chaque PM un poids d'énergie Wi. Le poids d'énergie Wi est le produit de la différentielle de la fonction d'optimisation de localisation d'événement M par rapport à la valeur attendue de l'énergie ni pour le PMi, multipliée par la valeur attendue de l'énergie pour le PMi, et multipliée en outre par la différentielle du logarithme naturel de la valeur attendue de l'énergie in ai du PMi par rapport à la distance du PMi de l'événement. Le bloc 38 montre que les poids d'énergie Wi pour chaque PM sont tous ajoutés ensemble pour fournir un facteur de normalisation W. Les premières coordonnées corrigées par itération (Xcl, Ycl) sont alors fournies par le bloc 39 en utilisant les distances des PM par rapport aux événements Xi, Yi, multipliées par les poids d'énergie des PMi et divisées par le facteur de normalisation W. On utilise alors les coordonnées Xcl et Ycl à la place des coordonnées non corrigées Xu, Yu provenant du circuit de calcul 16 pour l'itération suivante. Lorsque les nouvelles coordonnées corrigées Xcn, Ycn (o n = 2, 3, 4,...) produisent une convergence, comme le détermine le bloc 41, alors les valeurs de Xcn et Ycn sont les valeurs corrigées qui sont utilisées. Le bloc 42 représente une variante du système. Le bloc 37 définit le cas général de la valeur d'énergie Wi en se basant sur une fonction générale d'optimisation de localisation d'événement M (1, q2, ---, a1, a2). Dans un cas plus spécifique o l'on utilise une approche par vraisemblance logarithmique maximale, alors sont applicables les équations du paragraphe f) ci-dessus. On doit noter que, bien que le calcul d'Anger soit montré comme étant utilisé pour obtenir une première approximation des localisations en X et Y de l'événement, n'importe quel autre procédé pourrait être utilisé pour obtenir la première approximation de la localisation. Cependant, le calcul d'Anger est préféré puisqu'il est facilement disponible et bien connu. On détermine alors, comme le montre le bloc 34, la valeur ui, la distance de l'événement par rapport au PM, comme étant l'hypothénuse correspondant à la racine carrée de la somme des carrés de la quantité (X - Xi) et de la quantité (Y - Yi). A partir des valeurs de ui obtenues du bloc 34 on détermine les valeurs de ai, la valeur attendue de la sortie d'énergie de chaque photomultiplicateur PMi pour cette distance ui. On peut par conséquent produire une table de a(ui) qui est universelle pour chaque caméra à rayons gamma. Ainsi, on obtient ai à partir de la table de consultation produite en y entrant par les valeurs de ui calculées. Une autre table de consultation (la table III) donne la valeur de B(u) qui est égale à la différentielle du logarithme naturel de a(u) pour chaque tube photomultiplicateur. A partir des valeurs attendues d'énergie ai, on peut calculer une évaluation d'énergie Wi par PM, o {Wi = (i - ai) Bi} pour tous les photomultiplicateurs de l'ensemble. La valeur W est la valeur de la sommation de tous les Wi pour chaque tube photomultiplicateur. En utilisant les valeurs de W et Wi et les valeurs connues Xi, Yi on obtient les valeurs de coordonnées corrigées Xc et Yc en utilisant les équations: Xcn = E XiWi / W, et Ycn = E YiWi / W. Le bloc 37 ou le bloc 42 donne les valeurs de Wi. Les valeurs de Wi sont utilisées au bloc 38 pour obtenir la valeur de W. La valeur de Wi est également utilisée au bloc 39 en même temps que les valeurs de Xi, Yi et W pour obtenir les valeurs de coordonnées d'événement corrigées Xcn et Ycn On effectue une autre itération comme l'indique la boucle d'itération 25, s'il est répondu négativement à la question posée au bloc 41, concernant la fin des itérations, et si la réponse est affirmative on détermine immédiatement Xc et Yc- On utilise le retour de boucle d'itération 43 lorsque la réponse est non. Puis, les nouveaux Xcn et Ycn corrigés sont utilisés à la place des derniers Xc et Yc corrigés. Ceci fournit une autre réponse pour Xcn et Ycn qui est plus proche de la réponse vraie. Le nombre d'itérations effectuées dépend du but de l'image. Ainsi, on peut déterminer que c'est la fin des itérations lorsque l'amélioration se réduit à une fraction insignifiante. On peut utiliser d'autres critères à l'intérieur de la portée de l'invention. En pratique, trois ou quatre itérations sont plus que suffisantes. Les équations présentées sous forme de blocs fournissent un système très efficace pour obtenir de manière précise, les localisations en X et Y corrigées des événements sans qu'il soit nécessaire d'utiliser une quantité excessive de mémoire et un matériel de mise en oeuvre excessif. Une autre manière pratique d'obtenir de façon itérative des estimations améliorées de localisation d'événement à partir des équations: E (X - Xi) a Wi = 0; Z (Y - Y1) a Wi = 0 réalise un procédé itératif du type Newton-Raphson. En supposant que (X, Y) soit la localisation d'événement connue antérieurement et que (X + AX, Y + AY) soit la localisation d'événement améliorée, les équations de localisation d'événement se liront maintenant comme Z (AX + X - Xi) À Wi (X + AX, Y + AY) = O Z (AY + Y - Yi) ' Wi (X + AX, Y + AY) = O Wi dépend de (X, Y) seulement par l'intermédiaire de u2 = (X - Xi)2 + (y - Y1)2 par conséquent jusqu'au premier ordre en (AX, AY), on obtient awi Wi (X+AY,Y+AY) = Wi(X,Y) + 2 [(X-Xi)OAX=(Y-Yi)*AY] aui On substitue maintenant cette approximation dans les équations de localisation d'événement et seuls les termes allant jusqu'aux termes du premier ordre inclusivement en (tX, AY) sont retenus. Par conséquent, on obtient un ensemble de deux équations linéaires couplées en (AX, AY): sx + Sxx A + + Sxy * AY = O Sy + SXY 5X + Syy AY = 0 o: Sx = E (X - Xi) + Wi Sy = Z (Y - Yi) + Wi 2 aWi Sxx = Z [Wi + 2 * (X Xi)2 *] aui aWi Sxy = Z2 * (X - Xi) (Y - Y1) - Dui 2 wi Syy = E [Wi + 2 * (Y - Yi] * aui La solution pour (AX, AY) est: Sx * Syy - Sy * Sxy Sxx Syy - Sxy2 Sy * Sxx - Sx * Sxy AY = Sxx Syy Sxy En général, le terme Sxy de croisement (X,Y) est faible; cest-à-dire ISxyl possibilité de négliger le terme Sxy. Ainsi, en abandonnant les termes Sxy, la solution pour l'amélioration de la localisation d'événement (AX, AY) se réduit à: -Sx AX = Sxx -Sy AY = Syy L'estimation de localisation améliorée (Xcn, Ycn) est maintenant donnée par: X1 = X + AX Y1 = Y + AY O (X, Y) sont les coordonnées de localisation d'événement connues antérieurement. REVENDICATIONS 1. Caméra à rayons gamma comprenant: un cristal à scintillation (13) qui fournit un événement sous la forme d'une scintillation (17) en réponse à un impact d'un rayonnement (18) sur le cristal (13); plusieurs tubes photomultiplicateurs (PM) (14), ayant approximativement des courbes de PM symétriques radialement, montés sur ledit cristal (13) et mis en oeuvre en réponse à ladite scintillation (17) pour produire des signaux électriques; caractérisée en ce qu'elle comprend en outre: un calculateur de distances (34) pour calculer les distances, en coordonnées et suivant l'hypoténuse, de chacun des PM par rapport à chacune desdites scintillations (17), les coordonnées X, Y des scintillations (17) étant estimées et les coordonnées Xi, Yi des PM (14) étant connues; une source de valeur attendue d'énergie (36) fournissant une valeur attendue d'énergie (ai) pour chaque PM (14) en se basant sur la distance suivant l'hypoténuse du PM (14) par rapport à l'événement telle qu'elle est déterminée par le calculateur de distances (34); un calculateur de pondération (37) pour déterminer un facteur à appliquer à la sommation des coordonnées Xi, Yi des PM (14) pour fournir des coordonnées de localisation en X, Y corrigées (Xc, Yc) des scintillations (17); et des circuits d'itération pour utiliser de manière répétée, dans le calculateur de distance (34), les dernières coordonnées de localisation en X, Y corrigées des scintillations (17), pour déterminer des facteurs qui donnent des coordonnées en X, Y, pour la localisation de la scintillation (17), qui sont davantage corrigées que les dernières coordonnées de localisation en X, Y corrigées. 2. Caméra à rayons gamma selon la revendication 1, caractérisée en ce que des tables de consultation (31, 32, 33) fournissent les données de système utilisées pour la détermination des facteurs. 3. Caméra à rayons gamma selon la revendication 2, caractérisée en ce que lesdites tables de consultation comprennent: une première table (31) qui fournit les positions en coordonnées Xi, Yi de chacun desdits photomultiplicateurs (14); une deuxième table (32) qui fournit une valeur attendue de l'énergie normalisée pour chacun des photomultiplicateurs (14); et une troisième table (33) qui fournit la différentielle par rapport à la distance suivant l'hypoténuse du logarithme naturel de la valeur attendue de l'énergie normalisée pour chacun des photomultiplicateurs (14). 4. Caméra à rayons gamma selon la revendication 1, caractérisée en ce que la distance entre chacun des photomultiplicateurs (14) et l'événement est déterminée par l'équation: Ui2= (X - Xi)2 + (Y - Y1)2 o: U est la distance entre l'événement et le centre du tube photomultiplicateur (14); Xi et Yi sont les coordonnées au niveau du centre du tube photomultiplicateur (14); et X, Y sont les coordonnées estimées de l'événement. 5. Caméra à rayons gamma selon la revendication 1, caractérisée en ce que le facteur est constitué d'un poids de correction Wi par tube photomultiplicateur divisé par la sommation des poids de correction Wi. 6. Caméra à rayons gamma selon la revendication , caractérisée en ce que lesdites localisations en X, Y corrigées, Xc, Yc sont obtenues à partir des équations: Z Xiwi Xc = w YWi Z Yiwi Yc = w o Wi est un facteur de correction pour chacun des PM (14) et W est la sommation Z Wi. 7. Caméra à rayons gamma selon la revendication 6, caractérisée en ce que le facteur de correction Wi par tube photomultiplicateur est égal à la valeur attendue d'énergie ai multipliée par la différentielle d'une fonction d'optimisation de localisation d'événement M multipliée par la différentielle du logarithme naturel de la valeur attendue par rapport à la distance suivant l'hypoténuse U entre le tube photomultiplicateur (14) et la scintillation (17). 8. Caméra à rayons gamma selon la revendication , caractérisée en ce que le facteur de correction Wi par tube photomultiplicateur est égal à l'énergie par tube photomultiplicateur Zi divisée par la sommation Z Zi moins la valeur attendue ai de l'énergie normalisée par tube photomultiplicateur multipliée par la différentielle du logarithme naturel de la valeur attendue ai par rapport à la distance suivant l'hypoténuse U entre le centre du tube photomultiplicateur (14) et la scintillation (17). 9. Caméra à rayons gamma comprenant: un cristal à scintillation (13) qui fournit un événement qui est une scintillation (17) en réponse à un impact d'un rayonnement (18) sur le cristal (13); plusieurs tubes photomultiplicateurs (14), ayant approximativement des courbes de PM symétriques radialement, couplés audit cristal (13) et mis en oeuvre en réponse audit événement pour produire des signaux électriques; une approximation de localisation d'événement mise en oeuvre en réponse auxdits signaux électriques pour déterminer les localisations en X, Y approchées de l'événement et l'énergie Z de l'événement, X, Y étant des coordonnées orthogonales imaginaires du cristal à scintillation (13); un dispositif analytique utilisant les localisations en X, Y approchées de l'événement (17) par rapport aux localisations des tubes photomultiplicateurs (14) pour déterminer de manière analytique une valeur attendue ai de la sortie d'énergie de chacun des tubes photomultiplicateurs (14) en réponse audit événement, en se basant sur les localisations approchées des tubes photomultiplicateurs (14) par rapport à la localisation de l'événement; un calculateur de facteurs de correction (21) qui optimisera le caractère correct des localisations en X, Y approchées des événements, lorsque l'on utilise lesdits facteurs de correction pour modifier, de manière itérative, la sommation des coordonnées de localisation de PM. 10. Procédé de détermination de localisations corrigées d'événements dans des caméras à rayons gamma comprenant les étapes: de scintillation en réponse à l'impact sur un cristal (13) d'un rayonnement (18) pour donner une scintillation (17); de production de signaux électriques en réponse à ladite scintillation (17); caractérisé en ce qu'il comprend en outre les étapes: de calcul de distance, en coordonnées et suivant l'hypoténuse, de chacun de tubes photomultiplicateurs (14) par rapport à chacune desdites scintillations (17), les coordonnées X, Y des scintillations (17) étant estimées et les coordonnées Xi, Yi des PM (14) étant connues; de délivrance d'une valeur attendue d'énergie (ai) pour chaque PM (14) en se basant sur la distance suivant l'hypoténuse du PM (14) par rapport à l'événement, telle qu'elle est déterminée par le calculateur de distances (34); de détermination d'un facteur à appliquer à la sommation des coordonnées Xi, Yi des PM (14) pour fournir des coordonnées de localisation en X, Y corrigées (Xc, Yc) des scintillations (17); et d'utilisation de manière répétée des dernières coordonnées de localisation en X, Y corrigées des scintillations (17), pour déterminer des facteurs qui donnent des coordonnées en X, Y corrigées, pour la localisation de la scintillation (17), qui sont davantage corrigées que les dernières coordonnées de localisation en X, Y corrigées. 11. Procédé selon la revendication 10, caractérisé en ce qu'il comprend l'étape d'obtention de données de système, pour utilisation à la détermination des facteurs, à partir de tables de consultation (31, 32, 33). 12. Procédé selon la revendication 11, caractérisé en ce que lesdites tables de consultation comprennent: une première table (31) qui fournit les positions en coordonnées Xi, Yi de chacun desdits photomultiplicateurs (14); une deuxième table (32) qui fournit une valeur attendue de l'énergie normalisée pour chacun des photomultiplicateurs (14); et une troisième table (33) qui fournit la différentielle par rapport à la distance suivant l'hypoténuse du logarithme naturel de la valeur attendue de l'énergie normalisée pour chacun des photomultiplicateurs (14). 13. Procédé selon la revendication 10, caractérisé en ce qu'il comprend l'étape de détermination de la distance entre chacun des photomultiplicateurs (14) et l'événement en utilisant l'équation: Ui2 = (X - Xi)2 + (Y - Y1)2 o: U est la distance entre l'événement et le centre du tube photomultiplicateur (14); Xi et Yi sont les coordonnées au niveau du centre du tube photomultiplicateur (14); et X, Y sont les coordonnées estimées de l'événement. 14. Procédé selon la revendication 10, caractérisé en ce que le facteur est constitué d'un poids de correction Wi par tube photomultiplicateur divisé par la sommation des poids de correction Wi. 15. Procédé selon la revendication 14, caractérisé en ce qu'il comprend l'étape d'obtention desdites localisations en X, Y corrigées, Xc, Yc à partir des équations: E XiWi Xc = W z YiWi Yc = W o Wi est un facteur de correction pour chacun des PM (14) et W est la sommation Z Wi. 16. Procédé selon la revendication 15, caractérisé en ce qu'il comprend l'étape d'obtention du facteur de correction Wi par tube photomultiplicateur comme étant égal à la valeur attendue d'énergie ai multipliée par la différentielle d'une fonction d'optimisation de localisation d'événement M multipliée par la différentielle du logarithme naturel de la valeur attendue par rapport à la distance suivant l'hypoténuse U entre le tube photomultiplicateur (14) et la scintillation (17). 17. Procédé selon la revendication 14, caractérisé en ce qu'il comprend l'étape d'obtention du facteur de correction Wi par tube photomultiplicateur comme étant égal à l'énergie par tube photomultiplicateur Zi divisée par la sommation S Zi moins la valeur attendue ai de l'énergie normalisée par tube photomultiplicateur multipliée par la différentielle du logarithme naturel de la valeur attendue ai par rapport à la distance suivant l'hypoténuse U entre le centre du tube photomultiplicateur (14) et la scintillation (17). 18. Procédé d'obtention de coordonnées de localisation correcte d'une scintillation dans une caméra à rayons gamma, comprenant les étapes: de création d'un événement qui est une scintillation (17) en réponse à l'impact d'un rayonnement sur un cristal à scintillation (13); de couplage de plusieurs tubes photomultiplicateurs (14), ayant approximativement des courbes de PM symétriques radialement, audit cristal (13) et de mise en oeuvre de ceux-ci en réponse audit événement pour produire des signaux électriques; de détermination de localisations en X, Y approchées de l'événement en réponse aux signaux électriques, X, Y étant des coordonnées orthogonales imaginaires du cristal à scintillation (13); de détermination de façon analytique d'une valeur attendue ai de la sortie d'énergie de chacun des tubes photomultiplicateurs (14) en réponse audit événement, en se basant sur les localisations approchées des tubes photomultiplicateurs (14) par rapport à la localisation de l'événement; d'optimisation du caractère correct des localisations en X, Y approchées des événements, lorsque l'on utilise lesdits facteurs de correction pour modifier, de manière itérative, la sommation des coordonnées de localisation de PM.