Brevet FR2753279A1 — publié le 19980313

Procédé de traitement des impulsions délivrées par une qamma caméra et qamma caméra mettant en oeuvre ce procédé. L'invention concerne un procédé de traitement des impulsions délivrées par une gamma caméra et une gamma caméra mettant en oeuvre ce procédé. Elle se rapporte aux caméras à scintillation, ou gamma caméra, par exemple 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. L'utilisation des gamma caméras peut se généraliser à l'étude de tout corps radioactif. Le principe de fonctionnement des gamma caméras est représenté sur la figure 1. Un corps radioactif 1 émet des photons gamma y. Ces photons sont émis d'un point du corps 1, potentiellement, dans toutes les directions. Un scintillateur 2 reçoit ces photons gamma y et transforme au moment de l'impact, chaque photon gamma y reçu en source lumineuse s. Un ensemble de tubes photomultiplicateurs 3 disposé en regard du scintillateur 2 transforme en signaux électriques le rayonnement lumineux issu de la source s dans le scintillateur. Chaque tube photomultiplicateur fournit un signal électrique dépendant de la quantité totale de lumière qu'il a reçue. Des matrices de pondération 4 (cinq en général), constituées par exemple d'impédances de même type, transforment les signaux provenant des différents tubes photomultiplicateurs en cinq signaux x+, x-, y+, y et w. Les signaux x+, x, y+, et y contiennent des informations relatives aux positions en X et Y du barycentre de la répartition de la lumière sur le plan de détection constitué par les faces d'entrée des tubes photomultiplicateurs 3. Le signal w représente l'énergie totale récupérée sur l'ensemble de tubes photomultiplicateurs 3. Un circuit de calcul de position 5 récupère les signaux x+, x-, y+, et y pour les intégrer sur une durée équivalente à la durée d'une scintillation et déterminer la position en X et Y de l'impact. Ce circuit de calcul de position 5 fournit deux signaux x et y proportionnel à la position en X et Y. Un circuit de détection d'impact 6 reçoit le signal w et détermine en fonction de ce signal w représentatif de l'énergie reçue par l'ensemble de tubes 3 un signal de validation v qui indique si l'impact est à prendre en compte ou à ignorer. Un dispositif d'affichage 7 reçoit les signaux x, y et v et affiche ou enregistre, s'il y a lieu, le point d'impact de position X et Y, une image étant constituée par l'accumulation de points d'impact. D'autres dispositifs de mesure, sans matrices de pondération, sont envisageables. Leurs principes est de numériser les signaux à la sortie de chaque tube photomultiplicateur et d'utiliser des microprocesseurs rapides pour effectuer des calculs équivalents. De même, des dispositifs à semi-conducteur peuvent remplacer les tubes photomultiplicateurs. Dans la réalité, le corps radioactif 1 émet un certain nombre de photons gamma y dans un intervalle de temps donné. On parle d'un nombre de coups émis par seconde. Le nombre de coups émis par seconde dépend de l'activité du corps radioactif 1. De plus ces coups sont émis de manière aléatoire depuis n'importe quel point du corps radioactif 1 et dans n'importe quelle direction. En fait pour avoir une image, par projection en deux dimensions, représentative du corps radioactif 1, seul les photons gamma y venant d'une direction, par exemple normale au scintillateur 2, sont intéressants. Pour cela, il est possible d'utiliser un collimateur 8 pour collimater les photons gamma y émis par le corps radioactif 1 sur le scintillateur 2. Ce collimateur pourra être, par exemple, constitué d'une plaque en plomb épaisse de plusieurs centimètres percée d'une multitude de trous orientés par exemple perpendiculairement à sa surface afin de ne laisser passer que les rayonnements gamma normaux à cette surface. Le fait de rajouter un collimateur a pour effet de faire un tri. En effet, seul les photons venant perpendiculairement au collimateur de l'exemple traversent celui- ci. Cela représente un coup sur dix mille environ. Les coups sont d'une durée très brève et l'impact associé à un coup est assimilé à une fonction de Dirac dans le temps. Cependant l'impact sur le scintillateur 2 déclenche une cascade de phénomènes à travers le scintillateur 2, l'ensemble de tubes photomultiplicateurs 3, les matrices de pondération 4 et les connexions avec le circuit de calcul de position 5 et le circuit de détection d'impact 6, et donne finalement naissance à une impulsion relativement étalée. Pour localiser l'impact, on intègre, dans le circuit de calcul 5, sur une durée représentative de l'impulsion étalée, les signaux x+, x-, y+, et y pour pouvoir calculer ce qui correspond à un barycentre de la tache de lumière reçue par l'ensemble de tubes photomultiplicateurs 3. Mais si deux coups sont très rapprochés, les impulsions étalées peuvent se chevaucher, faussant ainsi le calcul de la position de l'impact. Le circuit de détection d'impact 6 a pour rôle de déterminer si un impact est à prendre en compte. La détection se faisant en même temps que le calcul de position d'impact, il est fourni au dispositif d'affichage une information qui valide la position en lui associant un niveau d'énergie. Anciennement, la détection était traitée ainsi: - détection du franchissement par le signal w d'un seuil de bruit, - intégration du signal w à partir de ce seuil détecté, - détection des maximums sur ce signal w, - détection du retour du signal w en dessous du seuil de bruit qui met fin à l'intégration du signal w, - élaboration d'un signal de validation en fonction du résultat de l'intégration du signal w, si par ailleurs le signal w est compris entre un maximum et un minimum autorisé, s'il n'y a eu qu'un seul maximum détecté lorsque le signal w était au dessus du seuil de bruit. En fait, une détection comme celle-là revient à vérifier que l'allure temporelle de l'énergie reçue correspond à un certain gabarit. Ce gabarit permet de vérifier qu'il y a bien une seule impulsion par la détection d'un maximum unique dans une plage d'énergie donnée. De ce fait, les impulsions de faible énergie dues à des diffusions COMPTON de photons gamma sont également exclues. La maîtrise des techniques numériques a permis d'améliorer le système. Dans la demande de brevet EP-A-0 470 909, il est décrit un système qui numérise les signaux x+, x-, y+, y et w. Puis à l'aide d'un filtre numérique, dont la fonction de transfert est calculée pour être l'inverse de la fonction de transfert des éléments qui transforment un impact de photon en un signal électrique, on transforme ces signaux x+, x-, y+, y- et w. Le filtrage numérique a pour effet de réduire la durée impulsionnelle des signaux de manière considérable. Avant filtrage, la durée d'une impulsion est un peu supérieure à la microseconde. Après filtrage, cette durée peut être inférieure à 250 ns (en fait la durée dépend entre autre de la fréquence d'échantillonnage et des paramètres du filtre numérique). Les impulsions étant de durée plus courte, les chevauchements d'impulsions sont moins nombreux qu'auparavant. Dans ces conditions, on peut compter plus de coups par durée unitaire et, en définitive, acquérir les images plus rapidement. La détection s'y fait néanmoins par un système de gabarit. Un système avec collimateur utilisant un filtrage numérique peut être utilisé pour compter jusqu'à 200 000 coups par seconde, ce qui représente une émission de l'ordre de 2 milliards de coups par seconde, ce qui correspond à environ 50 à 60 mCi. Or en médecine nucléaire, rendre le patient aussi radioactif avec les isotopes classiquement utilisés présente certains risques pour le patient. En général les systèmes collimatés sont utilisés avec les patients pour traiter au plus quelques dizaines de milliers de coups par seconde. Par ailleurs, on sait qu'il existe des isotopes émetteurs de positon (ou positron) o les photons gamma sont émis par paire dans des directions opposées, par exemple le fluor 18. Une autre technique pour collimater les photons gamma, dans le cas de positons, consiste à utiliser une gamma caméra avec un deuxième détecteur et un circuit de détection de coïncidence qui vérifie la direction de provenance du photon émis dans l'autre direction. Dans le cas général, la gamma caméra utilise deux détecteurs identiques. Le circuit de détection de coïncidence est relié aux deux détecteurs de la gamma caméra. La figure 2 représente un tel système. Sur la figure 2, une source radio active 10 émet des photons gamma yl et y2 par paire dans des directions opposées. De chaque coté de la source radio active 10, se trouve un scintillateur 11 et 12 pour transformer l'énergie des photons gamma yl et y2 en énergie lumineuse. On observe qu'un tel système n'utilise pas de collimateurs. L'énergie lumineuse émise par chacun des deux scintillateurs 11 et 12 est respectivement reçue par deux groupes de tubes photomultiplicateurs 13 et 14. Les groupes de tubes photomultiplicateurs 13 et 14 sont respectivement reliés à des ensembles de matrices de pondération 15 et 16. Les ensembles 15 et 16 transforment chacun l'énergie récupérée depuis les groupes de tubes photomultiplicateurs 13 et 14 en respectivement cinq signaux. Quatre signaux servent à la localisation de l'impact du photon gamma y1 ou y2 sur le scintillateur 11 ou 12, un signal étant représentatif de l'énergie totale dégagée par le photon gamma y1 ou y2. Des circuits respectifs de calcul de position 17 et 18 récupèrent chacun les quatre signaux, en provenance des matrices de pondération 15 et 16, qui sont représentatifs de la position de l'impact pour calculer les coordonnées du point d'impact et les transmettre à un dispositif d'affichage et de mémorisation 19. Un autre circuit de calcul permet de localiser l'origine des deux rayonnements gamma y1 et y2. Des variantes sont possibles au niveau du détecteur; les signaux peuvent être numérisés à la sortie des tubes photomultiplicateurs, rendant les matrices superflues; les tubes photomultiplicateurs peuvent également être remplacés par des dispositifs équivalents. Des circuits de détection d'impact 20 et 21 récupèrent le signal représentatif de l'énergie totale en provenance respectivement des ensembles de matrices 15 et 16. Ces circuits de détection d'impact 20 et 21 vont d'une part vérifier que le coup est valide, c'est à dire conforme au gabarit, et d'autre part vont calculer l'intégrale de puissance pour chaque coup valide. Un signal de validation, si le coup est valide, va ensuite être transmis au dispositif d'affichage 19 et à un détecteur de coïncidence 22. Si deux événements valides arrivent en même temps, c'est à dire dans une fenêtre temporelle par exemple de l'ordre de la dizaine de nanosecondes, le détecteur de coïncidence 22 émet un signal binaire pour signaler au dispositif d'affichage 19 qu'il s'agit bien d'une coïncidence. Le système d'affichage 19 peut reconstruire les trajectoires des photons et ensuite des images. Avec un patient, les émissions de photon gamma se font dans toutes les directions, quelque fois en même temps. En conditions cliniques, les coïncidences détectables ne représentent en général que 1% des coups reçus par les scintillateurs 11 et 12. Les 99% des coups reçus sont simples, c'est à dire qu'un seul des deux photons émis touche un détecteur de la gamma caméra. Dans la pratique, les gamma caméras utilisant la coïncidence, mesurent quelques centaines de milliers de coups par seconde, ce qui revient à détecter et à enregistrer quelques milliers de coïncidence par seconde. L'avantage par rapport à un appareil à collimateur réside essentiellement dans la réduction de dose radioactive à utiliser dans le patient. La dose injectée dans le patient est 100 fois moins radio- active pour un temps d'acquisition équivalent. Un article de Mankoff, Muehllehner et Karp, paru dans Phys. Med. Biol., 1989, Vol. 34, N 4, pages 437-456 et intitulé "The high count rate performance of a two- dimensionally position-sensitive detector for positron emission tomography" dévoile un comptage rapide utilisant la coïncidence. Dans cet article, il est question de compter deux millions de coups par seconde en utilisant un traitement numérique permettant de traiter un coup sur moins de 250 ns. Néanmoins à deux millions de coups émis par seconde, il n'est possible de récupérer que 62% des coups qui sont pourtant émis de manière expérimentale pour produire une coïncidence pour chaque coup émis. Les autres 38% correspondent à des coups non séparés temporellement par la gamma caméra utilisée. Si on augmente le taux d'émission des coups, le pourcentage des coups utilisables serait encore inférieur. Si l'on transpose le dispositif de l'article ci- dessus à un usage médical avec patient, le nombre de coups reçus et exploitables en coïncidence doit être divisé par cent. Quoi qu'il en soit, l'augmentation du nombre de coup émis par le patient tend à se faire chevaucher les coups reçus et donc à diminuer le nombre de coups valides au delà d'un certain nombre de coups émis de l'ordre de deux millions. Les examens utilisant des gamma caméras sont relativement longs car ils nécessitent l'acquisition d'images nécessitant chacune un nombre de coups important (pour qu'une image soit bien définie, plusieurs millions de coups peuvent être nécessaires). Un facteur permettant de diminuer le temps d'examen est d'augmenter le nombre de coups à prendre en compte par seconde. Or actuellement la technologie limite le nombre de coups à accepter. En résumé, l'état de l'art est le suivant: - la gamma caméra de la figure 1 utilisant un collimateur, ne permet pas de mesurer plus de quelques dizaines de milliers de coups par seconde sur un patient bien que technologiquement elle est tout à fait capable de prendre en compte 200 000 coups par seconde, elle nécessite également l'injection d'isotope de forte activité dans le patient, - la gamma caméra de la figure 2 utilisant un système de détection par coïncidence permet de réduire fortement les doses de produits radioactifs à faire absorber par le patient en diminuant légèrement le nombre de coup pris en compte (doses 100 fois moins importantes à nombre de coups égal), néanmoins comme on peut le voir dans l'article cité, avec des activités élevées, de l'ordre de 2 millions de coups, plus de 30 des coups utilisables sont perdus; un tel système permet de récupérer environ 14 000 coups sur 2 millions, de plus la résolution spatiale est en générale meilleure qu'avec les systèmes collimatés. L'invention a pour but d'augmenter le nombre de coups acceptables pour une gamma caméra à détection de coïncidence sans nécessiter ni de modification technologique de grande envergure, ni surtout d'augmentation de dose de radio-activité à injecter au patient. L'invention combine la détection par coïncidence avec une compression des signaux reçus, par exemple par déconvolution, pour en raccourcir la durée des impulsions correspondant à l'impact d'un photon gamma sur un détecteur, la détection de coïncidence se faisant sur les signaux avant déconvolution, sans vérification de la validité du coup à l'aide d'un gabarit. Il en résulte une certaine erreur du fait de tests manquant qui peut être compensée partiellement par un seuil d'énergie minimum. Mais du fait que les coups sont mieux séparés temporellement, on peut se satisfaire de la simple détection de simultanéité. Ainsi l'invention propose un procédé de traitement des signaux délivrés par une gamma caméra comprenant deux détecteurs de part et d'autre d'un corps radioactif, chacun desdits détecteurs produisant des signaux électriques analogiques constitués d'impulsions relatives à des impacts de photons gamma sur lesdits détecteurs, caractérisé en ce que pour chaque détecteur: - on produit au moins un signal comprimé, chaque signal comprimé étant représentatif d'un des signaux électriques sur lequel la durée des impulsions a été raccourcie; - on calcule les coordonnées du point d'impact à partir des signaux comprimés; et en ce que parallèlement: - on détecte les impulsions sur au moins un signal électrique analogique de chaque détecteur; - on produit un signal de coïncidence si une impulsion est détectée, dans une fenêtre temporelle associée à la simultanéité, sur au moins deux des signaux électriques analogiques provenant respectivement des deux détecteurs. Un tel procédé utilise, dans certains cas, des données fausses (coïncidences fortuites et chevauchement d'impulsions), ce qui va à l'encontre de ce qui était recherché jusqu'à présent. Dans l'invention, on a découvert que ce procédé permet de récupérer des données qui présentent un faible taux d'erreur qui est tout à fait acceptable pour de la tomographie. Bien sur, des données complètement erronées sont prises en compte; elles se traduisent par un apport de bruit ou de flou sur l'image récupérée. Le bruit qui correspond aux coïncidences fortuites augmente linéairement avec l'augmentation de l'activité du corps. Il est possible d'augmenter l'activité du patient sans réduire le taux de coïncidences exploitables, une augmentation de bruit étant à prévoir sur l'image. L'invention a aussi pour objet une gamma caméra comportant deux détecteurs produisant chacun des signaux électriques analogiques relatifs aux impacts de photon gamma sur lesdits détecteurs, caractérisé en ce que ladite gamma caméra comporte: des moyens de compression pour comprimer des impulsions sur les signaux électriques analogiques et fournir des signaux comprimés, des moyens pour détecter la présence d'impulsions sur au moins deux des signaux électriques analogiques associés respectivement aux deux détecteurs, des moyens de coïncidence pour détecter si les impulsions des au moins deux signaux électriques analogiques se produisent dans une fenêtre temporelle définie. Pour le procédé ou la gamma caméra, des améliorations sont possibles. Préférentiellement, on utilise des données numérisées afin de simplifier la réalisation de la gamma caméra. Un filtrage numérique ayant une fonction de transfert correspondant à l'inverse de la fonction de transfert de l'un des détecteurs est utilisé. De plus, plutôt que de déclencher les calculs de position ainsi que les intégrations de signaux lorsque l'on considère que le coup reçu est valide, on préférera une architecture qui permet de déclencher une opération à chaque cycle d'horloge réalisant ainsi plusieurs opérations simultanément mais décalées d'un cycle d'horloge, indépendamment de la détection de coïncidence. L'invention sera mieux comprise et des avantages apparaîtront mieux à la lecture de la description, non limitative, qui va suivre à lire conjointement aux dessins annexés parmi lesquels: - la figure 1 représente schématiquement une gamma caméra à collimateur selon l'art antérieur, - la figure 2 représente schématiquement une gamma caméra à coïncidence selon l'art antérieur, - la figure 3 représente schématiquement un mode de réalisation préférée d'une gamma caméra selon l'invention, - les figures 4 à 6 représentent schématiquement des réalisations préférées de circuits utilisés dans le mode de réalisation préféré de la gamma caméra de la figure 3. Les figures 1 et 2 ayant été décrites précédemment, on ne les décrira pas plus en détails. La figure 3 représente un mode de réalisation non limitatif de l'invention. On peut voir sur cette figure deux scintillateurs 23 et 24 placés de chaque coté d'un corps radioactif 25, émetteur de positons, qui émet des photons gamma y' et y" par paires qui sont émises dans des directions opposées. Dans la pratique, 1% des paires ont les deux photons qui atteignent les deux scintillateurs. De plus les photons y' et y" sont émis aléatoirement en lieu, direction et fréquence suivant une loi qui dépend de la densité de radioactivité du corps radioactif. Lorsque un photon y' ou y" arrive en contact avec un scintillateur 23 ou 24, le scintillateur transforme l'énergie du photon en énergie lumineuse qui part dans toutes les directions. L'énergie lumineuse émise par le scintillateur 23 ou 24 est récupérée par un groupe de tubes photomultiplicateurs 26 ou 27, chaque scintillateur 23 ou 24 ayant un groupe de tubes photomultiplicateurs 26 ou 27 qui lui est associé. Chacun des tubes photomultiplicateurs de chaque groupe 26 ou 27 fournit un signal électrique proportionnel à la quantité de lumière reçue du scintillateur qui est elle-même dépendante de la position du tube à considérer par rapport à l'impact du photon y' ou y" sur le scintillateur 23 ou 24. Dans le mode de réalisation préféré, on utilise deux ensembles de matrices de pondération 28 et 29 constituées par exemple de résistances qui sont respectivement connectés aux deux groupes de tubes photomultiplicateurs 26 et 27. Chaque ensemble de matrices de pondération 28 ou 29 est constitué de cinq matrices pondérant différemment les signaux issus de chaque tube photomultiplicateur du groupe de tubes photomultiplicateurs 26 ou 27 associé à l'ensemble de matrices 28 ou 29. Les pondérations effectuées par chacun des ensembles de matrices 28 ou 29 fournissent cinq signaux x+, x, y+, y-, et w dépendant directement de la quantité de lumière instantanée dégagée par le scintillateur 23 ou 24. Les signaux x+, x-, y+, et y sont représentatifs de la répartition de la tache de lumière due à l'impact du photon y' ou y" sur le scintillateur 23 ou 24 suivant les axes X et Y du scintillateur 23 ou 24 à considérer. Le signal w est représentatif de l'énergie lumineuse totale dégagée par le scintillateur 23 ou 24 après l'impact du photon y' ou y" L'exemple préféré de réalisation des détecteurs a été choisi car il est le plus connu à l'heure actuelle et de ce fait est plus simple à comprendre pour l'homme de I métier. Bien évidemment, il existe d'autres systèmes auxquels l'invention est applicable qui sont plus ou moins performant suivant leurs états d'évolution. Il existe par exemple des systèmes sans matrices o une numérisation des signaux se fait directement à la sortie de chaque tube photomultiplicateur, l'invention nécessite une dérivation de la sortie d'un ou plusieurs tubes photomultiplicateurs afin de récupérer les signaux analogiques. Il existe également des systèmes o les tubes photomultiplicateurs sont remplacés par des dispositifs à semi-conducteur équivalent. Par la suite, on appellera également détecteur ou tête de détection le dispositif constitué par exemple d'un scintillateur 23 (ou 24), d'un groupe de tubes photomultiplicateurs 26 (ou 27) et d'un ensemble de matrices de pondération 28 (ou 29). Deux circuits de calcul de position 30 et 31 sont respectivement connectés aux deux ensembles de matrices de pondération 28 et 29. Chacun des circuits de calcul de position 30 ou 31 reçoit les signaux x+, x-, y+, et y. Le but de ces circuits de calcul 30 et 31 est de déterminer la position de l'impact d'un photon en fonction des quatres signaux reçus. De nombreuses réalisations de ce genre de circuit sont connues, aussi bien analogiques que numériques et pourraient être utilisées avec l'invention. On préférera cependant un circuit de traitement numérique. On décrira ultérieurement le circuit de la figure 4 qui est un exemple de réalisation possible d'un de ces circuits de calcul de position 30 et 31. Chacun des circuits de calcul 30 et 31 fournit deux informations X et Y. Dans notre exemple, les informations X et Y sont supportées par deux bus de 16 bits, cependant ces informations peuvent très bien être codées sur un autre nombre de bit o même être analogiques. Dans le cas o les signaux sont numérisés directement à la sortie de chaque tube photomultiplicateur, un bus peut suffire pour que tous les tubes photomultiplicateurs puissent communiquer avec des circuits de calcul qui seront bien évidement adaptés. De plus le calcul de l'énergie du coup reçu peut être également effectué dans ces mêmes circuits de calcul. Deux convertisseurs analogique/numérique 32 et 33 sont respectivement connectés aux ensembles de matrices de pondération 28 et 29 pour transformer les signaux w en informations numériques d'énergie codée sur 16 bits dans notre exemple. Deux circuits de déconvolution 34 et 35 sont respectivement connectés aux convertisseurs analogique/numérique 32 et 33. Ces circuits de déconvolution 34 et 35 sont en fait des filtres numériques d'une structure classique qui ont pour fonction de transfert l'inverse de la fonction de transfert des éléments traduisant l'impact d'un photon y' (ou y") sur le scintillateur 23 (ou 24) en un signal électrique w (que l'on appellera par la suite fonction de transfert initiale). Cette fonction de transfert peut être obtenue par exemple en réalisant la transformée de Fourier inverse de l'inverse de la transformée de Fourier de la fonction de transfert initiale. Pour cela, on considère que la transformée de Fourier de la fonction de transfert initiale est égale à la transformée de Fourier de la réponse impulsionnelle à un impact d'un photon y' (ou y"). Pour plus de renseignement, on peut se reporter au document EP-A-0 470 909 et plus particulièrement au passage, colonne 7 ligne 11 à colonne 8 ligne 3, qui est considéré comme inclus par référence dans la présente demande, de même pour les figure 3 et 4 de cette demande EP-A-0 470 909. Deux intégrateurs 36 et 37 sont connectés respectivement à un bus de sortie de chacun des circuits de déconvolution 34 et 35. Ces intégrateurs 36 et 37 ont une structure qui leur permet d'additionner plusieurs valeurs d'informations successives sortant des circuits de déconvolution 34 et 35. De plus ces intégrateurs 36 et 37 commencent et finissent à chaque cycle d'horloge une intégrale bornée selon la méthode de Riemans en un certain nombre de cycles d'horloge. Ces deux intégrateurs 36 et 37 fournissent à leurs sorties respectives une donnée, codée par exemple sur 16 bits, qui est représentative de la puissance totale reçue par le groupe de tubes photomultiplicateurs 26 ou 27 durant une durée finie, cette donnée pouvant être, dans certains cas, représentative de l'énergie d'un impact reçu par le détecteur. Deux circuits détecteurs d'impulsions 38 et 39 sont également connectés respectivement aux ensembles de matrices de pondération 28 et 29 pour recevoir les signaux électriques analogiques w. Chacun de ces circuits détecteurs d'impulsions 38 ou 39 fournit un signal binaire représentatif de la présence d'une impulsion sur les signaux électriques analogiques w, en se positionnant dans un premier état ou dans un deuxième état. Dans le cas o l'on utilise un signal provenant d'un seul tube photomultiplicateur (de préférence situé au milieu du détecteur), on récupère le même train d'impulsions mais celles-ci sont d'amplitudes plus réduites et variables en fonction de la distance entre le point d'impact du photon et la position du tube photomultiplicateur. Le détecteur d'impulsion devra être plus sensible. Un circuit de coïncidence 40 récupère les deux signaux binaires issus desdeux circuits détecteurs d'impulsions 38 et 39. Ce circuit de coïncidence 40 positionne un signal binaire dans un premier état s'il y a coïncidence et dans un deuxième état s'il n'y a pas coïncidence. Autrement dit si les deux circuits détecteurs d'impulsions 38 et 39 indiquent simultanément une impulsion dans une fenêtre temporelle donnée, le circuit de coïncidence 40 indique alors une coïncidence, une simple porte logique à deux entrées adaptées aux états choisis étant suffisante. Le signal en sortie de la porte logique peut être mis en forme, par exemple à l'aide d'un circuit monostable. Un circuit de mémorisation 41 reçoit de la part des circuits de calcul de position 30 et 31 les coordonnées des points d'impact sur chaque scintillateur à chaque période d'horloge ainsi qu'une énergie associée à chacun de ces points. Or la quasi-totalité des points reçus sont faux. Le signal binaire indiquant la coïncidence va déterminer si les points doivent être mémorisés ou pas. Il faut toutefois faire attention à la synchronisation des informations arrivant des différents circuits 30, 31, 36, 37, et 40. En effet, La synchronisation va dépendre de la fréquence d'échantillonnage et de cadencement des éléments numériques, de la précision des circuits de déconvolution, et du nombre de cycles sur lequel on va intégrer les différents signaux. Le circuit de mémorisation 41 peut également faire un tri rapide des impacts valides en utilisant un seuil minimum d'énergie pour supprimer les erreurs dues aux diffusions importantes. Dans le circuit de mémorisation 41, les points valides auxquels sont associés des énergies sont donc mémorisés et éventuellement affichés ou traités. Un traitement d'image, indispensable pour obtenir une image claire, est de reconstruire les trajectoires (par exemple repérées en coordonnées sphériques) en fonction des points d'impact sur chacun des détecteurs. A partir des trajectoires, différentes images sont constituées par projection droite sur des plans perpendiculaires à la direction de la trajectoire et passant par un point central du corps en examen. Il est également possible de faire subir différents traitements à l'image obtenue tels qu'une élimination de bruit, un ajout de contraste, ou une reconstitution en trois dimensions, ou autres traitements d'image divers et variés qu'il n'est pas nécessaire d'expliquer dans le cadre de l'invention. La figure 4 représente une réalisation préférée des circuits de calcul de position 30 et 31. Ainsi le circuit de calcul de position 30 possède quatre entrées pour les signaux x+, x, y+, et y. A chacune des entrées est connectée un convertisseur analogique/numérique 42 à 45 du même type que les convertisseurs 32 et 33. Chacun de ces convertisseurs 42 à 45 est connecté par l'intermédiaire d'un bus à l'entrée d'un des circuits de déconvolution 46 à 49. Chacun de ces circuits de déconvolution 46 à 49 est identique aux circuits de déconvolution 34 et 35. A la sortie de ces circuits de déconvolution 46 à 49, on récupère des informations pour les introduire dans des intégrateurs 50 à 53. Les intégrateurs 50 à 53 étant du même type que les intégrateurs 36 et 37. Un circuit de calcul 54 de la position en X est connecté à la sortie des deux intégrateurs 50 et 51 qui sont reliés aux signaux x+ et x-. Ce circuit de calcul 54 fournit à sa sortie une information par exemple égale à X = (Ex+ - Ex-) / (Zx+ + Ex-), o Ex+ et Ex- sont les informations en sortie des intégrateurs 50 et 51 qui correspondent aux signaux x+ et x-. Un circuit de calcul 55 de la position en Y est connecté à la sortie des deux intégrateurs 52 et 53 qui sont reliés aux signaux y+ et y-. Ce circuit de calcul 55 fournit à sa sortie une information égale à Y = (Ey+ - Zy-) / (Ey+ + Ey-), o Ey+ et Ey- sont les informations en sortie des intégrateurs 52 et 53 qui correspondent aux signaux y+ et y-. D'autres possibilités de calcul sont envisageables pour calculer les coordonnées de l'impact sur le détecteur. De même, différentes possibilités de réalisation des circuits de calcul 54 et 55 sont possibles. Préférentiellement, on utilise des circuits microprogrammés plus connus sous le nom de DSP (de l'anglais Digital Signal Processor). La figure 5 montre le mode de réalisation préféré pour les intégrateurs 36, 37, 50 à 53. L'intégrateur 36 comporte un certain nombre de registres 65 à 68 égal au nombre de périodes d'échantillonnage à prendre en considération pour le calcul de l'intégrale, en pratique ce sera 8 à 10 périodes donc 8 à 10 registres. Entre les registres 65 à 68, on a placé des additionneurs parallèles 69 à 71, soit un additionneur 69 à 71 de moins que de registres 65 à 68. Chacun des additionneurs 69 à 71 reçoit des données en provenance de la sortie du registre précédent 65 à 67 et des données en provenance de l'entrée de l'intégrateur 36 pour fournir un résultat à l'entrée du registre suivant 66 à 68. Le premier registre 65 de la chaîne a son entrée connectée à l'entrée de l'intégrateur 36 et le dernier registre 68 de la chaîne a sa sortie connectée à la sortie de l'intégrateur 36. Un tel montage présente l'avantage de pouvoir commencer une intégration à chaque cycle d'horloge sans avoir besoin de savoir si le calcul est valide ou pas. Afin de rester dans un format de 16 bits, le bit de poids le plus faible en sortie de l'additionneur 69 à 71 n'est pas récupéré dans le registre suivant 66 à 68. Si l'on dispose d'un DSP très performant, on peut envisager d'implémenter les intégrateurs dessus. La figure 6 représente le mode de réalisation préféré pour les circuits détecteurs d'impulsions 38 et 39. L'entrée du circuit détecteur 38 est connectée d'une part à l'entrée d'un circuit retardateur 72 et d'autre part à l'entrée d'un atténuateur inverseur 73. Le circuit retardateur 72 peut être constitué par exemple d'une ligne à retard qui retarde le signal de quelques nanosecondes entre son entrée et une sortie. L'atténuateur inverseur 73 correspond à un amplificateur linéaire ayant une entrée et une sortie qui amplifie le signal d'entrée par un coefficient K, K étant compris entre 0 et -1. Un additionneur analogique 74 ayant deux entrées et une sortie a sa première entrée connectée à la sortie du circuit retardateur 72 et sa deuxième entrée connectée à la sortie de l'atténuateur inverseur 73. Un comparateur 75 ayant deux entrées et une sortie a une première entrée connectée à O V et une deuxième entrée connectée à la sortie de l'additionneur analogique 74. Les éléments 72 à 75 forment un circuit équivalent à un discriminateur à fraction constante. La sortie du comparateur 75 est dans un premier état s'il n'y a pas d'impulsion puis passe dans un deuxième état dès que le signal subit une variation positive de pente supérieure en valeur absolue à K fois la tension précédente divisée par le retard causé par le circuit retardateur 72. La sortie du comparateur 75 revenant dans le premier état si la pente devient plus faible ou s'inverse ou si la pente est maintenue au delà du retard causé par le circuit retardateur 72. De cette façon pratiquement toutes des impulsions sont détectées. Un circuit monostable 76 ayant une entrée et une sortie a son entrée connectée à la sortie du comparateur 75, sa sortie étant confondue avec la sortie du circuit détecteur 38. Le circuit monostable 76 transforme le signal de sortie du comparateur 75 en signal logique de type binaire. Chaque détection d'une impulsion sur le signal d'entrée du circuit 38 entraîne un changement d'état du signal binaire de sortie du monostable 76 pour une durée prédéfinie, par exemple de l'ordre de la dizaine de nanosecondes. Le circuit monostable 76 ouvre une fenêtre temporelle d'une durée prédéfinie. Des variantes peuvent être envisagées en utilisant par exemple des détecteurs de maximum, par exemple constitués d'un circuit dérivateur suivit d'une détection de zéro, le tout mis en forme par un circuit monostable. Ce mode de réalisation préféré est du à de nombreux compromis. De nombreuses réalisations différentes peuvent être faites sans s'écarter de l'invention par exemple traiter des informations sur un nombre de bits différent, ne pas effectuer les troncatures énoncées, utiliser des circuits qui réalisent les opérations en un nombre de cycles différents, ou ne pas utiliser les mêmes niveaux logiques. Les éléments constituant le mode préféré de réalisation de l'invention doivent en fait être considérés de manière fonctionnelle et ne limitent pas l'invention à cette réalisation préférée. Il est également possible de s'éloigner un peu de la réalisation préférée sans s'éloigner de l'objet de l'invention. Certains choix faits pour la réalisation préférée ne sont pas obligatoires. En effet, il est possible d'avoir un dispositif entièrement analogique, ce qui rend plus complexe le calcul du filtre de déconvolution. Dans le mode de réalisation préféré, on utilise quatre signaux de localisation et un signal relatif à l'énergie. Il est tout à fait possible d'utiliser un nombre de signaux différents sans s'éloigner de l'invention, par exemple un total de 4 ou 6 signaux peut être envisagé (il faudra pour cela modifier le nombre de matrices sur chaque détecteur). De plus, nous avons choisi de détecter une partie des fronts montant des impulsions sur le signal représentatif de l'énergie alors qu'il est possible de détecter les impulsions sur un des signaux de localisation ou même sur un signal issu directement d'un tube photomultiplicateur ou équivalent. Il est même possible d'effectuer une redondance de détection en détectant les impulsions sur plusieurs signaux issus de chaque tête de détection. L'utilisation de l'ensemble de matrices dans une tête de détection peut être remplacée par un dispositif équivalent qui effectue un échantillonnage et une conversion numérique à la sortie de chaque tube photomultiplicateur. Des microprocesseurs rapides sont ensuite chargés d'effectuer le calcul de position et d'énergie. Egalement, le filtrage numérique qui a pour but de réaliser une déconvolution peut être remplacé par n'importe quel autre filtrage ou algorithme de traitement du signal qui réalise une séparation temporelle des impulsions par compression des impulsions. REVENDICATIONS 1. Procédé de traitement des signaux délivrés par une gamma caméra comprenant deux détecteurs (23, 26, 28, 24, 27, et 29) de part et d'autre d'un corps radioactif (25), chacun desdits détecteurs (23, 26, 28, 24, 27, et 29) produisant des signaux électriques analogiques (x+, x, y+ y, et w) constitués d'impulsions relatives à des impacts de photons gamma (y' et y") sur lesdits détecteurs (23 et 24), caractérisé en ce que pour chaque détecteur (23, 26, 28, 24, 27, et 29): - on produit au moins un signal comprimé, chaque signal comprimé étant représentatif d'un des signaux électriques analogiques (x+, x, y+, y-, et w) sur lequel la durée des impulsions a été raccourcie; - on calcule les coordonnées du point d'impact à partir des signaux comprimés; et en ce que parallèlement: - on détecte les impulsions sur au moins un signal électrique analogique de chaque détecteur; - on produit un signal de coïncidence si une impulsion est détectée, dans une fenêtre temporelle associée à la simultanéité, sur au moins deux impulsions des signaux électriques analogiques provenant respectivement des deux détecteurs (w). 2. Procédé selon la revendication 1, caractérisé en ce que les signaux électriques (x+, x, y+, y-, et w) issus des détecteurs (23, 26, 28, 24, 27, et 29) sont pondérés avec des ensembles (28 et 29) de matrices de pondération respectivement associés à chacun desdits détecteurs (23, 26, 28, 24, 27, et 29), les signaux électriques (x+, x, y+, y-, et w) étant ainsi représentatifs de la position d'impact des rayons gamma (r' et y") sur les détecteurs et de l'énergie. 3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que les signaux comprimés sont obtenus par un filtrage selon une fonction de transfert sensiblement inverse à la fonction de transfert correspondant à la transformation du rayonnement gamma (y' ou y") en signal électrique (x+, x-, y+, y, et w) par chacun des détecteurs (23, 26, 28, 24, 27, et 29). 4. Procédé selon l'une des revendications 2 ou 3, caractérisé en ce que les signaux électriques analogiques sur lesquels on détecte la présence d'impulsions correspondent aux signaux électriques analogiques représentatifs de l'énergie. 5. Procédé selon la revendication 4, caractérisé en ce que, à chacune des périodes de l'horloge: - on commence une intégration des signaux (x+, x y, y, et w) sur un nombre de périodes prédéterminé, - on calcule les positions d'impact (X et Y) des rayons gamma (y' et y") sur les détecteurs (23, 26, 28, 24, 27, et 29) à partir des signaux comprimés, relatifs à la position d'impact (x+, x, y+ et y), qui ont été intégrés sur ledit nombre de périodes, 6. Procédé selon la revendication 7, caractérisé en ce que le signal de coïncidence valide certaines des positions d'impact (X, Y) calculées. 7. Gamma caméra comportant deux détecteurs (23, 26, 28, 24, 27, et 29) produisant chacun des signaux électriques (x+, x, y+, y, et w) relatifs aux impacts de photon gamma (y' et y") sur lesdits détecteurs (23 et 24), caractérisé en ce que ladite gamma caméra comporte: - des moyens de compression (34, 35, 46 à 49) pour comprimer des impulsions sur les signaux électriques (x+, x-, y+, y-, et w) et fournir des signaux comprimés, - des moyens pour détecter la présence d'impulsions (38 et 39) sur au moins deux des signaux électriques associés respectivement aux deux détecteurs (w), - des moyens de coïncidence (40) pour détecter si les impulsions des au moins deux signaux électriques se produisent dans une fenêtre temporelle définie. 8. Gamma caméra selon la revendication 7, caractérisé en ce que qu'elle comporte des convertisseurs analogique/numérique (32, 33, 42 à 45) pour transformer les signaux électriques (x+, x-, y+, y, et w) en informations numériques et en ce que les moyens de compression (34, 35, 46 à 49) sont constitués de filtres numériques. 9. Gamma caméra selon la revendication 8, caractérisé en ce qu'elle comporte des circuits de calculs (30 et 31) pour donner la position des impact sur chacun des détecteurs (23 et 24), et en ce que ces circuits de calcul (30 et 31) ont une architecture qui leur permet de commencer un calcul et de finir un calcul à chaque période d'une horloge qui synchronise lesdits circuits de calcul (30 et 31).