Brevet FR2680878A1 — publié le 19930305

Procédé pour localiser des évènements de scintillation multiples L'invention concerne un procédé pour localiser des évènements de scintillation multiples dans une caméra gamma comportant un certain nombre de photomultiplicateurs, dont les signaux de sortie sont soumis à un processus de reconnaissance de formes, en fonction d'un signal de libération ou d'autorisation. L'un des facteurs limitatifs des capacités potentielles d'une caméra gamma ou d'une caméra Anger est la capacité de comptage pouvant être obtenue Dans les systèmes actuels, la limite se situe de façon typique à environ 000 évènements par seconde La raison principale en est que, dans sa forme de réalisation actuelle, la caméra gamma ne peut analyser qu'un seul évènement en l'espace de 1 à 2 ps et l'identifier d'une manière utile pour l'image obtenue en médecine nucléaire Le choix d'évènements valables s'effectue par exemple par l'intermédiaire d'une fenêtre ou gamme d'énergies, c'est-à-dire de l'amplitude d'impulsions à la sortie d'un détecteur. Si deux ou plusieurs quanta de rayonnement gamma pénètrent simultanément ou presque simultanément dans le détecteur, leur énergie totale peut dépasser la gamme ou plage prévue pour des évènements individuels Cependant, les25 quanta sont dispersés dans les caméras connues actuellement, de sorte que la cadence de comptage est fortement réduite Si l'intervalle de temps entre deux évènements individuels est suffisamment grand, des circuits spéciaux permettent de déterminer ce qu'on appelle pile-up ou un évènement double. Au mieux, on peut alors encore analyser de façon approximative et utiliser le premier évènement Le second est dispersé de sorte qu'on obtient à nouveau une réduction importante de la cadence de comptage Si l'intervalle de temps d'un évènement double est trop faible, les deux évènements doivent être rejetés. Dans EP-A-0 155 463, on indique un montage analogue pour l'analyse d'évènements multiples A cet effet, il est prévu un réseau servant à établir la communication entre les photomultiplicateurs de la caméra gamma, mais à l'aide duquel on peut identifier et séparer les évènements multiples,15 lorsque leur intervalle est suffisamment important pour que des groupes de signaux séparés dans l'espace puissent être formés Même dans le cas d'un intervalle de temps important entre les évènements doubles, les groupes de signaux peuvent se perturber réciproquement de sorte que l'analyse et la20 localisation des deux évènements individuels associés sont faussées. Un procédé du type indiqué plus haut est décrit dans l'article de Milster et consorts: "A Full Field Modular Gamma-Camera" paru dans The Journal of Nuclear Medicine,25 Vol 31, N 04, Avril 1990, pages 632 à 639 Dans ce document, on indique une caméra gamma modulaire, qui délivre des informations images utiles concernant l'ensemble de la surface de cristal La localisation s'effectue après numérisation des signaux de sortie des photomultiplicateurs, au moyen de la fonction de probabilité maximale Chaque module est indépendant du point de vue optique et électronique de sorte que plusieurs modules peuvent être combinés pour former un système de caméra présentant une cadence de comptage accrue Cependant cet article n'indique35 pas de localiser également des évènements multiples pour accroître la cadence de comptage. L'invention a pour but d'agencer une caméra gamma de manière à pouvoir séparer, sur la base des signaux de sortie des photomultiplicateurs de la caméra, des évènements multiples possédant des énergies qui diffèrent d'une manière quelconque, c'est-à-dire des évènements qui sont constitués par des évènements individuels possédant l'énergie primaire correcte ou même des évènements individuels apparaissant conjointement et possédant une énergie primaire correcte et des évènements parasites, c'est-à-dire de manière que des évènements multiples puissent être subdivisés en les évènements individuels. Le problème est résolu grâce au fait que lors du processus de reconnaissance des formes, les signaux de sortie sont comparés à des ensembles multiples de signaux de comparaison, qui comprennent respectivement les valeurs escomptées des signaux de sortie des photomultiplicateurs, qui sont fournies par plusieurs évènements de scintillation de comparaison ayant des origines dont les emplacements sont connus, et que les origines sont enregistrées en tant qu'emplacements des évènements de scintillation multiples, qui appartiennent à l'ensemble multiple de signaux de comparaison, qui présente la valeur maximale de similitude par rapport aux signaux de sortie. Grâce à ce procédé, on peut identifier et localiser des évènements multiples qui apparaissent aussi bien simultanément que presque simultanément Bien que l'intervalle de temps entre les évènements individuels constituant un évènement multiple soit nettement inférieur à l'intervalle de temps nécessaire pour identifier un évènement individuel, ces évènements individuels peuvent encore être localisés Par conséquent, la cadence de comptage de la caméra gamma est accrue De ce fait, on peut réduire par exemple la durée de mesure nécessaire pour la formation de l'image ou bien, avec ce procédé, on peut également obtenir des signaux d'image utile de régions d'examen qui présentent une activité accrue En outre, en plus de la localisation, on peut également déterminer l 'énergie des événements individuels Ceci permet de localiser et de déterminer l'énergie d'évènements primaires et d'évènements parasites apparaissant simultanément Ceci peut s'appliquer également dans le cas d'évènements possédant des énergies primaires différentes, et par exemple on peut exécuter des études sur deux isotopes, ou bien on peut utiliser des isotopes possédant plusieurs énergies. Dans un mode de mise en oeuvre avantageux, les ensembles de signaux de comparaison d'un évènement respectif de scintillation de comparaison sont mémorisés sous la forme d'un ensemble individuel de signaux de comparaison dans une mémoire et les ensembles multiples de signaux de comparaison sont formés, à partir de là, dans une unité de calcul Ce processus de reconnaissance des formes requiert une faible place de mémoire. Un autre mode de mise en oeuvre avantageux se caractérise par le fait que chaque ensemble multiple de signaux de comparaison est une combinaison linéaire de plusieurs ensembles individuels de signaux de comparaison et que, lors du processus de reconnaissance de formes, des facteurs sont modifiés dans la combinaison linéaire. Un mode de mise en oeuvre se caractérise par le fait que le signal de libération est délivré lorsque l'énergie de la somme de tous les signaux de sortie est située dans une gamme d'énergies, dont les limites sont situées autour de la valeur escomptée de l'énergie de plusieurs événements primaires de scintillation. Un autre mode de mise en oeuvre avantageux se caractérise par le fait que les signaux de sortie sont intégrés avant d'être soumis au processus de reconnaissance des formes L'intégrale des signaux de sortie présente une dispersion plus faible que l'amplitude des impulsions de signaux De ce fait, lors du processus de reconnaissance des formes, l'influence des variations statistiques de l'amplitude des signaux est réduite. Dans une autre variante particulièrement avantageuse de l'invention, la durée pendant laquelle les signaux de sortie sont intégrés, est modifiée en fonction du fait que les évènements de scintillation apparaissent simultanément ou peu de temps les uns après les autres, en se chevauchant dans le temps Une réduction de la durée d'intégration dans le cas d'évènements de scintillation apparaissant simultanément permet de réduire la durée de traitement et par conséquent d'accroître plus encore la cadence de comptage. Conformément à un mode de mise en oeuvre préféré, les ensembles multiples de signaux de comparaison de plusieurs évènements de scintillation de comparaison sont mémorisés dans une mémoire Les ensembles multiples de signaux de comparaison, qui sont nécessaires pour le processus de reconnaissance des formes, sont interrogés20 directement à partir de la mémoire afin d'être comparés aux signaux de sortie De ce fait, la durée de traitement dans le calculateur pour la localisation est faible. Dans une variante avantageuse, le signal de libération est délivré lorsque l'énergie de la somme de tous les signaux de sortie est située dans une gamme d'énergie, dont les limites sont situées autour de la valeur escomptée de l'énergie de plusieurs évènements primaires de scintillation Ceci permet de réduire la durée de calcul. On va expliquer ci-après des exemples de mise en oeuvre de l'invention en référence à onze figures. La figure 1 représente un schéma-bloc d'une première forme de réalisation d'une caméra gamma servant à localiser des évènements multiples de scintillation. La figure 2 représente une vue en plan d'un cristal à scintillation comportant des photomultiplicateurs. La figure 3 représente des valeurs escomptées des signaux de sortie d'une série de photomultiplicateurs. La figure 4 représente l'agencement d'une unité de signaux de comparaison pour des ensembles de signaux de comparaison. La figure 5 représente un schéma-bloc d'une seconde forme de réalisation d'une caméra gamma pour la localisation d'évènements multiples de scintillation. La figure 6 représente un autre agencement de l'unité de comparaison. La figure 7 représente un signal somme d'un évènement double, dont les évènements apparaissent simultanément. La figure 8 représente un signal somme d'un évènement double, dont les évènements apparaissent peu de temps l'un après l'autre. La figure 9 représente un évènement double différentié selon la figure 8. La figure 10 représente un signal somme intégré. La figure 11 représente des fonctions d'images de points d'évènements doubles en fonction de l'intervalle entre les évènements de scintillation. Le schéma-bloc de la figure 1 représente, sous la forme d'une vue d'ensemble, l'agencement d'une caméra gamma, à l'aide de laquelle on peut localiser des évènements multiples de scintillation à l'aide d'un procédé de reconnaissance des formes On a représenté ici schématiquement, sous la forme d'une vue en élévation latérale, un cristal à scintillation 2, qui est constitué par30 exemple par du iodure de sodium (Na J) Le cristal à scintillation 2 est accouplé optiquement par l'intermédiaire d'une fenêtre 4, constituée par exemple par du verre Pyrex, à un nombre N de photomultiplicateurs PM. On peut voir sur la figure 2 une vue en plan de la disposition du cristal à scintillation 2, de la fenêtre optique 4 et des photomultiplicateurs PM Vingt photomultiplicateurs PM 1 à PM 20 possédant un diamètre de ,8 mm, sont répartis selon un réseau rectangulaire 4 x 5 sur la surface de la fenêtre 4 Le cristal à scintillation 2 et la fenêtre optique 4 possèdent une surface rectangulaire de 254 x 203 mm 2 L'épaisseur du cristal à scintillation 2 est égal à 12,7 mm et l'épaisseur de la fenêtre 4 est égale à 6,4 mm Le cristal à scintillation 2 est métallisé Sur la figure 1, chaque photomultiplicateur PM est raccordé à un canal de traitement de signaux Un canal de traitement de signaux comprend un préamplificateur 6, une ligne à retard 8, un intégrateur 10 et un convertisseur analogique/numérique 12 Les canaux de traitement des signaux raccordent les photomultiplicateurs PM à un processeur de vraisemblance 14, dans lequel est exécuté le processus de reconnaissance de formes. Chaque photomultiplicateur PM est en outre raccordé à respectivement une entrée d'un additionneur 16 La sortie de l'additionneur 16 est raccordée à un analyseur 22 des formes d'impulsions. L'analyseur 22 des formes d'impulsions forme des signaux de commande qui sont envoyés aux intégrateurs 10 par l'intermédiaire d'entrées de commande 28 et du processeur de vraisemblance 14. Le processeur de vraisemblance 14 est raccordé à une unité 30 délivrant des signaux de comparaison dans laquelle, comme cela sera encore décrit plus loin de façon détaillée, sont formés et/ou mémorisés des ensembles multiples de signaux de comparaison Dans la forme de30 réalisation représentée sur la figure 1, l'unité 30 délivrant des signaux de comparaison comprend une première mémoire partielle 32, qui contient des ensembles de signaux de comparaison d'évènements individuels, et une unité de calcul 34, qui forme des ensembles multiples de signaux de comparaison d'évènements multiples, à partir des ensembles de signaux de comparaison mémorisés dans la première mémoire partielle 32. La sortie du processeur de vraisemblance 14 est raccordée par l'intermédiaire d'un seuil de similitude 36 à un calculateur central 38, qui, à partir des évènements localisés, établit une image de scintillation et la représente sur un moniteur 40 Le calculateur central 38 tient compte du fait que seuls des évènements apparaissant avec l'énergie primaire correcte doivent contribuer à former l'image de scintillation. Pour l'obtention des ensembles de signaux de comparaison mémorisés dans la mémoire partielle 32, l'ensemble de la surface 2 du cristal est exploré avec une source radioactive collimatée De cette manière, on obtient, en tous les emplacements xi explorés, des évènements de scintillation de comparaison, qui conduisent à l'apparition de signaux de sortie de photomultiplicateurs PM Les emplacements xi sont par exemple les coordonnées d'un système de coordonnées rectangulaire Étant donné que ces signaux sont soumis à une dispersion statistique, des valeurs escomptées des signaux de sortie sont mémorisées, c'est-à- dire que pour l'emplacement xi pour chaque photomultiplicateur, une valeur moyenne Se (xi) avec (e = 1, N, i = 1,, K) est formée à partir d'un grand nombre d'évènements et est enregistrée dans l'ensemble de signaux de comparaison Grâce à un choix correspondant d'activités de la source et d'une gamme d'énergies pour les quanta primaires individuels, on est certain que les valeurs escomptées correspondent de façon précise à l'absorption de quanta individuels au niveau des emplacements xi. A titre d'exemple d'exploration de l'ensemble de la surface 2 de cristal, la figure 3 indique des valeurs escomptées, qui sont obtenues au moyen d'une exploration le long d'un trajet d'exploration 41 Le trajet d'exploration 41 passe par les centres des photomultiplicateurs PN 9 à PM 12. Sur le diagramme de la figure 3, on a porté en abscisses les points d'exploration xi, qui sont séparés ici par une distance de 4 mm En ordonnées, on a marqué, en unités quelconques, la valeur escomptée Se(xi) des signaux de sortie des photomultiplicateurs individuels PM La valeur escomptée du photomultiplicateur PM 9 est indiquée au moyen de S (xi). Pour le point d'exploration 10, il s'avère par exemple que la valeur escomptée 59 ( 10) du photomultiplicateur est légèrement supérieure à 700, tandis que la valeur escomptée 51 o( 10) du photomultiplicateur PM 10 est égale approximativement à 280. Les valeurs escomptées 51 ( 10) et 512 ( 10) des photomultiplicateurs PM 1, et PM 12 sont inférieures à 100 et varient seulement faiblement dans la zone entourant le point d'exploration 10 Au-dessous des abscisses du diagramme, on a également représenté la position des photomultiplicateur PM, et PM 1, en rapport avec les emplacements d'exploration xi On peut voir que la valeur escomptée Se(xi) du photomultiplicateur PM est maximale lorsque l'emplacement d'exploration xi se situe au centre du photomultiplicateur correspondant. Les valeurs escomptées, ainsi déterminées, d'évènements individuels de scintillation de comparaison avec une origine dont l'emplacement est connu sont mémorisées dans la mémoire 32 comme indiqué sur la figure 4 Le premier ensemble de signaux de comparaison est constitué par les valeurs escomptées des signaux de sortie Si à SN de tous les photomultiplicateurs PM 1 à PMN, lorsque l'évènement de scintillation de comparaison est produit à l'emplacement x 1 dans le cristal à scintillation 2 Le second ensemble de30 signaux de comparaison contient des valeurs escomptées des signaux de sortie si à SNN lorsque l'évènement de scintillation de comparaison est produit à l'emplacement x 2. Lorsque l'ensemble de la surface du cristal à scintillation 2 est exploré en k points, la totalité des k ensembles de signaux de comparaison, qui sont dûs à des évènements primaires, sont mémorisés dans la mémoire 32. Les ensembles multiples de signaux de comparaison pour des événements multiples de scintillation peuvent être formés à partir des ensembles de signaux de comparaison pour5 des événements individuels qui sont mémorisés dans la mémoire 32 A cet effet, on se référera par exemple pour un événement double qui est constitué par deux événements individuels primaires, à la figure 3 L'ensemble de signaux de comparaison pour un événement double est constitué par la somme des valeurs escomptées des événements individuels produits aux emplacements correspondants Ainsi, la valeur escomptée du signal de sortie du photomultiplicateur P dans le cas d'un événement double aux emplacements 5 et 10 conformément à la figure 3 est égale à environ 1400.15 Cependant, non seulement des événements multiples de scintillation, qui proviennent d'événements primaires, sont localisés, mais également des événements multiples de scintillation, qui se composent par exemple d'un événement primaire et d'un événement parasite, dont l'énergie est inférieure à l'énergie d'un événement primaire Le cas indiqué en dernier se présente fréquemment dans le domaine de la pratique clinique et n'a pu être utilisé, avec des procédés connus jusqu'alors, que dans des conditions limites déterminées pour la détermination d'un scintigramme Des événements parasites ou d'une manière plus générale des événements présentant des énergies différentes, sont pris en compte par le fait qu'un ensemble multiple de signaux de comparaison est une combinaison linéaire de plusieurs ensembles individuels de signaux de comparaison Si par exemple un événement double doit être localisé, l'ensemble multiple de signaux de comparaison est formé à partir de deux ensembles individuels de signaux de comparaison, comme représenté sur la figure 4 Les facteurs a 1, a 2 a,,, an dans le cas de m événements déterminent l'énergie des35 événements individuels, dont est constitué l'événement il multiple considéré. En principe, la localisation se déroule de telle sorte que, par exemple, lors de la localisation d'un évènement double, la fonction de vraisemblance est rendue maximale non seulement par rapport aux emplacements xi, Xj, mais également par rapport au facteur a 1, a 2 Dans le cas d'évènements multiples, le nombre des paramètres à modifier augmente de façon correspondante Le résultat est une détermination d'emplacements et une détermination d'énergies.10 Pour accélérer la localisation, on n'utilise pas toutes les combinaisons imaginables xi, xj et a 1 et a 2 pour le calcul de la fonction de vraisemblance maximale, mais on se rapproche du maximum selon un procédé adaptatif avec une largeur de pas de plus en plus faible.15 La caméra gamma de la figure 1 localise aussi bien des évènements simples que des évènements doubles De façon usuelle, un signal somme Sges' dont l'allure de principe dans le temps est représentée sur les figures 7 et 8, est formé à partir des signaux de sortie des photomultiplicateurs PM,20 dans l'additionneur 16 La valeur maximale de l'amplitude est, dans le cas d'évènements individuels, une mesure de l'énergie du quantum de rayonnement gamma absorbée dans le cristal à scintillation 2 Cette valeur maximale est inférieure à un écart statistique Sur la figure 7, une courbe 43 représente par exemple le signal somme d'un évènement individuel de scintillation et une courbe 44 représente le signal somme de deux évènements individuels apparaissant simultanément ou d'un évènement double Deux évènements individuels peuvent cependant apparaître également dans le temps l'un après l'autre, c'est-à-dire presque simultanément, comme cela est représenté sur la figure 8 Les courbes 43 représentent à nouveau le signal somme des évènements individuels Le signal somme de l'évènement double correspondant est représenté par la courbe 46 et est obtenu à partir d'une addition des valeurs des deux courbes 43 Les évènements doubles représentés sur les figures 7 et 8 sont rejetés dans les caméras gamma actuelles, et dans le meilleur des cas, seul l'évènement, qui intervient en premier, peut être localisé sur la figure 8. Le déroulement de la localisation se déroule ici comme cela va être décrit ci-après L'analyseur 22 des formes d'impulsions délivre un signal lorsque le signal somme Sges dépasse un seuil qui est supérieur au bruit dans l'obscurité des photomultiplicateurs PM Un événement apparaît Les intégrateurs 10 et 28 sont ouverts pendant un intervalle de temps d'environ 1 à 2 ps Si pendant cet intervalle de temps d'intégration un autre événement apparaît, l'analyseur 22 des formes d'impulsions délivre à nouveau un signal, à la suite de quoi la durée d'intégration est prolongée de façon correspondante d'environ 1 à 2 ps, etc Si aucun évènement ne se produit pendant le dernier intervalle de temps de prolongement, l'intégration s'arrête et la localisation d'un évènement, intervenant m fois, se produit conformément au signal, qui apparaît m fois, de l'analyseur 22 des formes d'impulsions, pendant la durée d'intégration On indique que m peut être aussi égal à un. Si le maximum atteint de la fonction de vraisemblance maximale est inférieur à un seuil, un autre événement peut être masqué dans l'événement simple ou l'événement multiple localisé, cet autre événement n'ayant pas pu être détecté par l'analyseur 22 des formes d'impulsions, dans le cas d'un intervalle trop faible par rapport à l'événement précédent La localisation d'un événement apparaissant m+ 1 fois est alors déclenchée Si la localisation de l'événement, qui apparaît m+ 1 fois, ne conduit pas à un résultat satisfaisant, on essaie de localiser un événement apparaissant m+ 2 fois, etc Par conséquent, on localise les événements multiples qui peuvent se succéder dans le temps en étant suffisamment proches pour que l'analyseur 22 des formes d'impulsions n'identifie qu'un événement individuel Dans ce cas, la localisation est erronée Dans le cas d'un événement double, la localisation d'un événement individuel conduit à un emplacement situé entre les deux événements effectifs et notamment à une valeur très faible de la fonction de vraisemblance La valeur de similitude formée à partir de la fonction de vraisemblance est alors inférieure à une valeur de seuil, ce qui provoque le déclenchement du processus de reconnaissance des formes pour des événements doubles La localisation ultérieure en tant qu'événement double conduit alors la plupart du temps aux emplacements corrects et aux énergies correctes des événements individuels Le démarrage répété du processus de reconnaissance des formes avec des ensembles de signaux de comparaison, qui sont produits par un nombre, accru de un, d'événements de scintillation de comparaison, s'effectue toujours lorsque la valeur de similitude tombe au-dessous de la valeur de seuil Ce processus se poursuit jusqu'à ce que le procédé de reconnaissance des formes, avec l'ensemble multiple maximal possible de signaux de comparaison, ait été utilisé pour la localisation. Ci-après, en se référant à la figure 5, on va décrire une variante à l'aide de laquelle on localise des événements doubles, qui se composent de deux événements primaires Comme représenté sur la figure 6, les valeurs escomptées pour de tels événements doubles sont indiquées pour toutes les positions et pour tous les intervalles possibles des événements individuels primaires formant les événements doubles et sont mémorisées dans une seconde mémoire partielle 42 Le nombre des ensembles multiples des signaux de comparaison, mémorisés dans la mémoire 42, pour des événements doubles, est approximativement égal au carré du nombre des événements individuels mémorisés dans la mémoire 32. La sortie de l'additionneur est en outre envoyée, par rapport à la figure 1, à un intégrateur 20 par l'intermédiaire d'une ligne à retard 18 La sortie de l'intégrateur 20 est raccordée à deux discriminateurs monocanaux 24 et 26 Les sorties des discriminateurs monocanaux 24,26 sont raccordées au processeur de vraisemblance 14 Les intégrateurs 10 et 20 ainsi que les discriminateurs monocanaux 24 et 26 possèdent des entrées de commande 28, qui sont raccordées à la sortie de l'analyseur 22 des formes d'impulsions. Dans l'analyseur 22 des formes d'impulsions, la forme dans le temps de ce signal somme se maintient. L'analyseur 22 des formes d'impulsions identifie, par exemple par différentiation du signal somme représenté sur la figure 8, des événements individuels apparaissant d'une manière approchée La différentiation d'un tel évènement double fournit l'allure de signal représentée sur la figure 9 Si le signal somme différentié possède cette allure de caractéristiques pour un évènement double formé d'évènements apparaissant presque simultanément, un signal de commande est délivré sur la sortie de l'analyseur 22 des formes d'impulsions Si l'analyseur 22 des formes d'impulsions ne délivre aucun signal, parce que le signal somme ne possède pas la structure de pics doubles représentée sur la figure 8, aussi bien les signaux de sortie que le signal somme sont intégrés pendant un intervalle de temps standard prédéterminé de l'ordre de 1 à 2 ps dans les intégrateurs 10 et 20 Le signal somme intégré est représenté sur la figure 10 La valeur maximale E correspond à l'énergie d'un évènement individuel primaire de scintillation ou d'un évènement multiple primaire de scintillation C'est pourquoi on évalue avant tout les signaux sommes intégrés, étant donné qu'ils présentent une dispersion statistique moins grande que les signaux non intégrés Pour accroître la précision, les signaux devant être intégrés sont envoyés par l'intermédiaire des lignes à retard 8 et 18 de manière que le traitement des35 signaux soit déclenché uniquement lorsque l'allure, dans le temps, a été établie dans l'analyseur 22 des formes d'impulsions. Le signal somme intégré est alors envoyé aux deux discriminateurs monocanaux 24 et 26 Le premier discriminateur monocanal 24 examine le signal somme pour déterminer s'il se situe dans une plage d'énergies Eo AE, Eo étant caractéristique de l'énergie totale de quanta individuels primaires Le second discriminateur monocanal 26 examine le signal somme pour déterminer s'il se situe dans10 une plage d'énergies 2 E AE, qui est caractéristique de l'énergie totale de deux quanta individuels primaires, c'est- à-dire d'un évènement double. Si le signal somme ne tombe dans aucune des deux fenêtres, il s'agit d'un évènement multiple, dans lequel intervient au moins un évènement parasite De tels évènements multiples peuvent être localisés au moyen de la caméra gamma décrite en référence à la figure 1 Les intégrateurs 10 et 20 peuvent être ramenés à l'état initial, par exemple, au moyen d'un autre signal.20Si le discriminateur monocanal 24 répond, c'est qu'on est en présence d'un événement individuel primaire tout d'abord valable Ensuite, le convertisseur analogique/ numérique 12 est déclenché et le processeur de vraisemblance 14 est informé du fait qu'une comparaison à des ensembles individuels de signaux de comparaison doit être exécutée pour des évènements individuels du type de ceux mémorisés dans la mémoire 32 L'emplacement de l'ensemble de signaux de comparaison, qui possède la valeur de similitude maximale par rapport aux signaux de sortie, est considéré en tant30 qu'emplacement de l'évènement individuel devant être analysé. Cette valeur de similitude est en outre envoyée à l'unité à seuil de similitude 36 Si le résultat dépasse un seuil de similitude déterminé, l'évènement accepté comme étant un évènement individuel primaire est envoyé au calculateur 38. Si la valeur de similitude est inférieure au seuil, une localisation est déclenchée d'une manière correspondante, comme cela est décrit en référence à la figure 1 Il peut s'agir, par exemple, d'un événement double, qui est situé seulement de façon aléatoire dans la plage d'énergies, par exemple de deux événements simultanés dispersés, ou d'un événement d'une nature physique entièrement différente. Si, au contraire, le discriminateur monocanal 26 répond, un événement double tout d'abord valable est présent. Comme précédemment, le convertisseur analogique/numérique 1210 est déclenché, mais le processeur de vraisemblance 14 est maintenant informé qu'il faut exécuter une comparaison à des ensembles de signaux de comparaison pour des événements doubles, du type de ceux qui sont soit mémorisés dans la mémoire 42, soit calculés, conformément à la figure 4, à15 partir des ensembles de signaux de comparaison pour des événements individuels dans la mémoire 32 Cependant, ici les facteurs ai et a 2 sont positionnés à un " 1 ", étant donné que tout d'abord il n'y a à localiser aucun événement dispersé. L'emplacement (xi,xj) de l'ensemble de signaux de comparaison, qui possède la valeur de similitude maximale par rapport aux signaux de sortie, est considéré comme étant l'emplacement de l'événement double devant être analysé et qui est constitué par deux événements individuels à l'emplacement xi et à l'emplacement xj Comme dans le cas de l'analyse d'événements individuels, on peut en outre utiliser un seuil de similitude pour supprimer des événements d'une nature physique différente, qui n'interviennent seulement que de façon accidentelle dans la gamme ou plage d'énergies entourant 2 EO, ou faire démarrer une localisation comme cela30 est décrit en référence à la figure 1. S'il est établi dans l'analyseur 22 des formes d'impulsions que la forme du signal somme s'écarte de celle prévue pour des événements individuels ou pour des événements multiples simultanés et possède par exemple la structure à35 deux pics représentée sur la figure 8, c'est qu'il s'agit, avec une grande probabilité, d'un évènement double formé de deux évènements individuels apparaissant en étant séparés par un bref intervalle de temps Si on trouve une telle indication d'évènements doubles, l'analyseur 22 des formes d'impulsions informe les indicateurs 10 et 20 pour qu'il réalise une intégration dans un intervalle de temps prolongé. La valeur finale E du signal somme intégré est indépendante de l'intervalle de temps entre les évènements individuels constituant l'évènement double, en supposant que la durée d'intégration est prolongée de manière que l'intégration soit approximativement complète sur les deux évènements. Naturellement, pour chaque canal de signaux on peut également utiliser deux intégrateurs possédant des durées d'intégration différentes, mais fixes, et entre lesquelles une commutation peut être exécutée en fonction du résultat de l'analyseur 22 des formes d'impulsions Si on a établi l'existence d'une structure à deux pics dans le signal somme, le second discriminateur monocanal 26 détermine alors si l'amplitude du signal somme intégré se situe dans une gamme entourant le double de la valeur d'un évènement individuel 2 EM Si ce n'est pas le cas, comme cela a été décrit en référence à la figure 1, le processus de reconnaissance de formes est déclenché Si, cependant, le second discriminateur monocanal 26 répond, un évènement tout d'abord valable apparaît dans la gamme d'énergies d'un évènement double primaire A nouveau la conversion analogique/numérique des signaux de sortie intégrés est alors déclenchée et le processeur de vraisemblance 15 est informé du fait qu'il faut exécuter maintenant une analyse d'évènements doubles primaires, dont les ensembles de signaux de comparaison sont mémorisés dans la seconde mémoire partielle 42 ou, comme cela a été expliqué en référence à la figure 4, sont respectivement calculés de façon effective, les facteurs étant réglés à UN et n'étant pas modifiés L'évènement double ainsi analysé peut être alors envoyé directement au calculateur central 38 Si la sécurité de l'identification de formes doubles doit être accrue, la valeur de similitude, qui se trouve dans le processeur de vraisemblance 14, de l'évènement double est interrogée de manière à déterminer si elle est supérieure à5 un seuil de similitude, et éventuellement la procédure de reconnaissance de formes conforme à la figure 1 est déclenchée. On voit que les deux procédés de localisation décrits en référence aux figures 1 et 5 peuvent être avantageusement combinés La commande, décrite en référence à la figure 5, du procédé de reconnaissance de formes dans la gamme d'énergies réduit la durée de calcul Il n'est pas nécessaire de réaliser dans tous les cas une modification au moyen du facteur des combinaisons linéaires d'ensembles individuels de signaux de comparaison Si les évènements multiples sont constitués uniquement par des évènements primaires, une modification de l'emplacement de l'évènement est alors suffisante. Sur la figure 11, on a représenté la fonction de points d'images d'évènements primaires doubles, en fonction de la distance de ces évènements En abscisses, on a indiqué l'emplacement des évènements et en ordonnées, le nombre, déterminé au moyen d'une formation de la valeur moyenne, des quanta comptés Les courbes 50 représentent la fonction de points d'image d'une caméra gamma de la figure 1 pour deux évènements individuels simultanés ou presque simultanés séparés par une distance de 5 cm Les courbes 52, 54, 56 représentent de façon correspondante la fonction de points d'image d'évènements doubles séparés par une distance de 4 cm, 3 cm et 2 cm La distance des évènements doubles correspondant aux courbes 58 et 60 est inférieure à 2 cm On peut voir que des évènements qui sont séparés par une distance de 4 cm et de 5 cm peuvent être complètement dissociés en évènements individuels par la caméra gamma décrite De ce fait, la cadence de comptage est fortement 19 2680878 accrue par rapport à la caméra gamma usuelle De même, on peut encore séparer de façon conditionnelle des évènements doubles qui sont séparés par une distance de 2 cm et de 3 cm. C'est seulement au-dessous de 3 cm que l'on ne peut plus dissocier l'évènement individuel constituant des évènements doubles de sorte qu'à partir de là, on n'obtient aucune information concernant l'image. REVENDICATIONS 1 Procédé pour localiser des évènements de scintillation multiples dans une caméra gamma comportant un certain nombre de photomultiplicateurs (PM), dont les signaux de sortie sont soumis à un processus de reconnaissance des formes, en fonction d'un signal d'autorisation, caractérisé par le fait que lors du processus de reconnaissance des formes, les signaux de sortie sont comparés à des ensembles multiples de signaux de comparaison, qui comprennent respectivement les valeurs escomptées des signaux de sortie des photomultiplicateurs, qui sont fournies par plusieurs évènements de scintillation de comparaison ayant des origines dont les emplacements sont connus, et que les origines sont enregistrées en tant qu'emplacements des évènements de15 scintillation multiples, qui appartiennent à l'ensemble multiple de signaux de comparaison, qui présente la valeur maximale de similitude par rapport aux signaux de sortie. 2 Procédé suivant la revendication 1, caractérisé par le fait que les ensembles de signaux de comparaison d'un événement de scintillation de comparaison sont mémorisés sous la forme d'un ensemble individuel de signaux de comparaison dans une mémoire ( 32) et que les ensembles multiples de signaux de comparaison sont formés, à partir de là, dans une unité de calcul. 3 Procédé suivant la revendication 1, caractérisé par le fait que chaque ensemble multiple de signaux de comparaison est une combinaison linéaire de plusieurs ensembles individuels de signaux de comparaison et que, lors du processus de reconnaissance des formes, des facteurs sont modifiés dans la combinaison linéaire. 4 Procédé suivant l'une quelconque des revendications 1 à 3, caractérisé par le fait que le signal d'autorisation est délivré lorsque l'énergie de la somme de tous les signaux de sortie est située dans une gamme ou plage d'énergies, dont les limites sont situées autour de la valeur10 escomptée de l'énergie de plusieurs évènements primaires de scintillation. Procédé suivant l'une quelconque des revendications 1 à 4, caractérisé par le fait que les signaux de sortie sont intégrés avant d'être soumis au processus de reconnaissance des formes. 6 Procédé suivant la revendication 5, caractérisé par le fait que la durée, pendant laquelle les signaux de sortie sont intégrés, est modifiée en fonction du fait que les évènements de scintillation apparaissent simultanément ou20 peu de temps les uns après les autres, en se chevauchant dans le temps. 7 Procédé suivant l'une quelconque des revendications 1 à 6, caractérisé par le fait que le nombre des évènements de scintillation qui apparaissent en25 chevauchement, détermine la multiplicité des ensembles multiples de signaux de comparaison. 8 Procédé suivant l'une quelconque des revendications 5 à 7, caractérisé par le fait que les signaux de sortie intégrés sont numérisés avant d'être soumis au30 processus de reconnaissance des formes. 9 Procédé suivant la revendication 1, caractérisé par le fait que les ensembles multiples de signaux de comparaison de plusieurs évènements de scintillation de comparaison sont mémorisés dans une mémoire ( 42).35 10 Procédé suivant l'une quelconque des revendications 1 à 9, caractérisé par le fait que les évènements de scintillation de comparaison sont des évènements primaires. 11 Procédé suivant l'une quelconque des revendications 1 à 10, caractérisé par le fait que le processus de reconnaissance des formes est commandé en fonction du fait que l'énergie de la somme de tous les signaux de sortie se situe dans une gamme ou plage d'énergies, dont les limites sont situées autour de la valeur10 escomptée de l'énergie de plusieurs évènements de scintillation. 12 Procédé suivant l'une quelconque des revendications 1 à 11, caractérisé par le fait que les emplacements sont enregistrés uniquement lorsque la valeur de15 similitude dépasse une valeur de seuil. 13 Procédé suivant la revendication 12, caractérisé par le fait que le processus de reconnaissance des formes est à nouveau déclenché avec des ensembles multiples de signaux de comparaison, dont la quantité est accrue de UN, lorsque la valeur de similitude dépasse la valeur de seuil. 14 Procédé suivant l'une quelconque des revendications 1 à 13, caractérisé par le fait qu'un processus de reconnaissance des formes pour des évènements doubles est déclenché lorsqu'après une localisation d'évènements individuels lors d'un processus de reconnaissance des formes pour des événements individuels, une valeur de similitude formée lors du processus de reconnaissance des formes tombe au-dessous d'une valeur de seuil.