Brevet FR2605745A1 — publié le 19880429

-- i -- CAMERA A SCINTILLATION La présente invention a trait au domaine de la mé- decine nucléaire et concerne plus spécialement une caméra à scintillation comprenant un cristal scintillateur éventuelle- ment équipé d'un collimateur et destiné à convertir chaque photon reçu en une scintillation, un guide de lumière pour le couplage dudit cristal à la fenêtre d'entrée d'un jeu de p photodétecteurs destinés à convertir en courant chaque scin- tillation, p voies d'acquisition recevant les signaux de sor- tie desdits photodétecteurs et délivrant p signaux électriques de caractéristiques liées notamment à l'intensité de la scin- tillation et à la distance de cette scintillation à chacun des photodétecteurs, et un calculateur destiné à délivrer les coordonnées xj et yj d'une scintillation j et l'énergie Ej qui est associée à cet événement j. Pour déterminer l'image de fixations radioactives à l'intérieur d'un organe, la médecine fait appel, entre autres moyens, au principe de la scintigraphie. Ce principe consiste à introduire dans l'organisme d'un patient un élément radioactif qui va se fixer plus ou moins sur certains organes selon que ceux-ci sont sains ou malades. La mesure de l'inten- sité de rayonnement gamma émis fournit alors une indication sur la répartition de l'élément radioactif dans l'organisme et constitue donc une aide au diagnostic. Une telle mesure est réalisée à l'aide d'une caméra à scintillation. Dans les caméras à scintillation traditionnelles, par exemple de type Anger (le physicien Anger étant le premier à avoir proposé une caméra à scintillation, dont les principes fondamentaux sont décrits dans le brevet des Etats-Unis d'Amérique n0 3011057), les rayons gamma représentatifs de la distribution radioactive dans le milieu examiné pénètrent, - 2 - après traversée d'un collimateur, dans un cristal scintilla- teur. Les scintillations qui se produisent dans ce cristal sont alors détectées par l'intermédiaire de toute une série de tubes photomultiplicateurs (par exemple 37), après traversée d'un guide de lumière assurant un couplage optique entre le cristal et les tubes. Ces tubes sont répartis devant le bloc optique (cristal + guide de lumière) pour couvrir pratiquement toute sa surface et transformer l'énergie lumineuse de chaque scintillation apparue en un signal électrique mesurable. En sortie de chaque tube photomultiplicateur est alors prévue une voie d'acquisition analogique, opérant successivement une amplification, une intégration et une mise en forme des signaux fournis par le tube. Les sorties Sij de l'ensemble de ces voies d'acquisition sont envoyées vers un calculateur qui fournit par estimation les coordonnées xj et yj d'une scintillation j et son énergie Ej (l'indice i désigne celle des voies d'acquisition qui est concernée). Dans le calculateur, plusieurs types de dispositifs de calcul peuvent être prévus, mais deux d'entre eux essentiellement sont réellement utilisés, à savoir un dispositif de calcul barycentrique à rapport arithmétique, ou un dispositif de calcul barycentrique à rapport logarithmique. Dans un dispositif de calcul barycentrique à rap- port arithmétique, les grandeurs xjj,, Ej sont données par les expressions: xj = (1) zi yj = Yj (2) zi i=p Ej = 1 Gi Sij (3) 1=1 - 3 Dans ces expressions, on a: xj = E Ki Sij (4) i=1 Yj = E Hi Sij (5) i=1 j=p Zj = E ' Ji Sij (6) i=l Les coefficients Gi, Ki, Hi, Ji sont des facteurs de pondération liés à la position de l'axe de chacun des p tubes photomultiplicateurs. Dans un dispositif de calcul barycentrique à rapport logarithmique, les grandeurs xj, yj, Ej sont cette fois données par les expressions: x= Log (7) = J a7 xi yj = Log (8) a y i=p Ej =1 Gi Sij (9) + i=p + avec Xj =E Ki Sij (10) i= i - =p - Xj =i Ki Sij (11) i= i + i=p + Yi E Hi Sij (12) i= I - i-=p =jE Hi Si de l'axe de chacun des p tubes photomultiplicateurs. Dans l'un ou l'autre cas, quelle que soit l'arith- métique utilisée, une modification de calcul des coordonnées a été proposée pour améliorer la résolution spatiale intrinsèque de la caméra. Cette modification, désormais utilisée dans la - 4 - plupart des caméras à scintillation développées dans le do- maine de la médecine nucléaire et décrite par exemple dans le brevet des Etats-Unis d'Amérique US-A- 3 732 419, consiste à introduire un seuil so dans le calcul des signaux composites + - + - Xj, Yj, Zj ou Xj, Xj, Yj, Yj. Plus précisément, et en prenant par exemple le cas d'un dispositif de calcul barycentrique à rapport arithmé- tique, les signaux composites Xj, Yj sont maintenant donnés par les expressions: i=p Xj = 1 Ki (Sij-so) (14) i=1 i=p Yj = E Hi (Sij-so) - (15) i=1 avec, en outre, la convention Sij-so = 0 si Sij est in- férieur à so. Les relations (1), (2), (3) donnant xj, yj, Ej et la relation (6) restent inchangées. Cette introduction d'un seuil est réalisée au moyen d'un circuit à seuil introduit dans chacune des p voies analogiques d'acquisition. Ces p voies possèdent alors deux sorties, l'une pour les signaux Sij-so qui vont permettre le calcul de Xj et Yj, et l'autre pour les signaux Sij utilisés pour le calcul de Zj et Ej. La valeur so du seuil est fixe et a été réglée une fois pour toutes lors de la fabrication, afin d'obtenir la meilleure résolution spatiale à une énergie donnée, par exemple 140 keV. Cependant, comme ce réglage est justement indépen- dant de l'énergie réelle des événements, une dégradation des performances des caméras (notamment distorsions linéaires et pertes de résolution) se produit pour des énergies d'événe- ments trop différentes de celle utilisée lors du réglage à la fabrication. C'est en particulier le cas lorsqu'on effectue des examens multi-isotopiques. - 5 -- Le brevet des Etats-Unis d'Amérique US-A-4 475 042 remédie à cet inconvénient en proposant une caméra à scintil- lation à seuil dynamique, ledit seuil ayant une amplitude qui est une fonction du temps et de l'énergie de l'événement en cours de traitement. Ce seuil, appliqué comme le seuil fixe mentionné préalablement dans les p voies analogiques d'acqui- sition, en sortie des préamplificateurs associés aux tubes photomultiplicateurs, est constitué par la somme de deux si- gnaux. Le premier des deux signaux est une fraction du signal d'énergie qui résulte de la sommation pondérée des signaux fournis par tous les préamplificateurs, sans effet de seuil. Le deuxième signal est une fraction du signal d'énergie inté- gré qui correspond à la valeur intégrée du premier des deux signaux. A ces deux signaux, une composante continue peut encore, éventuellement, être superposée. Le signal de seuil ainsi constitué présente en fonction du temps des fluctuations de type statistique c'est- à-dire liées au nombre de photons relatif à chacun des évé- nement inclus dans l'observation. La valeur moyenne de ces fluctuations est liée à l'énergie de la scintillation, ou évé- nement, en cours de traitement. Ce seuil est retranché de cha- cun des signaux fournis à la sortie des préamplificateurs. L'amélioration ainsi apportée par le document US-A-4 475 042 permet en effet de préserver les qualités d'image lorsqu'on travaille avec plusieurs isotopes simulta- nément. Cette solution présente malgré tout, dans la pra- tique, des inconvénients sensibles: (a) les signaux qui composent le seuil dynamique ne sont obtenus qu'après des déformations et des retards dé- terminés par la fonction de transfert des circuits traversés; (b) le seuil dynamique est un signal qui fluctue en fonction du temps de façon aléatoire, et en outre varie d'un événement à l'autre, et ces fluctuations et variations statistiques d'amplitude relative et de phase différentes pour - 6 - les signaux issus des préamplificateurs et pour le seuil dy- namique conduisent à ce que les signaux délivrés soient par- fois, au lieu d'un signal nul ou d'un signal utile, du bruit qui dégrade à nouveau la résolution spatiale; (c) le seuil dynamique ne convient plus à t4ux de comptage élevé, lorsqu'il y a empilement au moins partiel des événements, car la valeur moyenne du seuil dynamique n'est dans ce cas plus seulement liée à l'énergie du seul événement en cours de traitement et est donc incorrecte; (d) parmi l'ensemble des tubes photomultiplica- teurs concernés par une scintillation déterminée, certains sont situés dans une position intermédiaire, à une distance telle de cette scintillation que l'amplitude maximale des im-- pulsions fournies par les préamplificateurs n'est que deux à trois fois supérieure à l'amplitude moyenne du seuil dyna- mique. Pour les voies d'acquisition associées à ces tubes, les impulsions de sortie sont alors incomplètes puisque le seuil dynamique élimine justement la partie basse terminale de ces impulsions. Le but de l'invention est donc de proposer une caméra à scintillation dans laquelle on a remédié à ces incon- vénients et en particulier au fait que les réalisations anté- rieures à seuil dynamique ne conviennent pas à taux de comptage élevé. A cet effet, la caméra à scintillation selon l'in- vention est caractérisée: (A) en ce que lesdites p voies d'acquisition réalisent l'amplification, le filtrage et l'échantillonnage desdits signaux de sortie des photodétecteurs, puis la conver- sion analogique-numérique des' échantillons obtenus et leur sommation, et délivrent p signaux numériques à l'entrée du calculateur; (B) en ce que le calculateur comprend lui-même: (a) un bus de transfert desdits p signaux numériques; - 7- (b) un étage de sommation numérique comprenant lui-même: - une mémoire-tampon recevant les p signaux de sortie Mij du bus de transfert; - un circuit soustracteur de seuil recevant sur une première entrée lesdits p signaux, retardés par la mé- moire-tampon et sur une deuxième entrée une valeur moj de seuil dit asservi et délivrant p signaux Mij-moj, ces signaux prenant conventionnellement la valeur 0 si Mij est inférieur à moj; - quatre dispositifs de sommation pondérée numé- rique dont les deux premiers reçoivent les p signaux Mij-moj délivrés par le circuit soustracteur de seuil et dont les deux autres reçoivent directement les p signaux de sortie Mij du bus de transfert; - un calculateur de seuil recevant la sortie du quatrième dispositif de sommation pondérée numérique pour l'évaluation du seuil asservi moj selon une relation de pro- portionalité à la valeur Emj associée à chaque événement j; ledit étage de sommation numérique délivrant les signaux Xm,j, Ym,j, Zm,j, Em,j suivants (ki, Hi, Ji, Gi, étant des coefficients de pondération): i=p Xmj = t Ki (Mij-moj) i=1 i=p ym, = Hi (Mij-moj i=1 i=p Zm,j = i= i i Emj = -1 Gi Mij i=1 soit directement en sortie des deux premiers dispositifs de sommation pondérée numérique, pour les deux premiers de ces -- 8 -- signaux, soit par l'intermédiaire de deux circuits de réali- gnement temporel prévus respectivement en sortie des deux autres dispositifs, pour les deux derniers de ces signaux; (c) un étage de traitement des événements in- cluant des circuits de calcul de désempilement et deux divi- seurs et délivrant à partir des signaux Xmj, Ymj, Zmj, Emj les trois signaux de coordonnées et d'énergie xj, yj, Ej relatifs à l'événement j; (C) en ce qu'un étage de détection, séquencement et stockage recevant un signal qui correspond à la somme des p signaux de sortie des photodétecteurs est prévu pour délivrer d'une part les différents signaux d'horloge pour la synchroni- sation des éléments de p voies d'acquisition et des éléments du calculateur et d'autre part des coefficients de correction destinés à l'étage de traitement des événements. Dans deux variantes de réalisation de la caméra ainsi proposée, le calculateur peut n'utiliser qu'une seule des deux grandeurs Z ou E pour le calcul des coordonnées, compte tenu des possibilités de correction ultérieure offertes par les caméras à scintillation actuellement disponibles. Le calculateur ne comprend dans ce cas que trois dispositifs de sommation pondérée numérique dont les deux premiers reçoivent les p signaux Mij-moj délivrés par le circuit soustracteur de seuil et dont le dernier reçoit directement les p signaux de sortie Mij du bus de transfert. Le calculateur de seuil reçoit alors la sortie dudit dernier dispositif de sommation pondérée numérique pour l'évaluation du seuil asservi moj selon une relation de proportionalité à la valeur Emj associée à chaque événement j. Dans la première de ces deux variantes, l'étage de sommation numérique délivre les signaux Xmj, Ym,j, Zm,j suivants (Ki, Hi, Ji, étant des coefficients de pondération): i=p Xm,i = f Ki (Mij-moj) 9- i=p Ym,j = [ Hi (Mij-moj i=1 i=p Zm,j = E Ji Mij i=1 soit directement en sortie des premier et deuxième dispositifs de sommation pondérée numérique, pour les deux premiers de ces signaux, soit par l'intermédiaire d'un circuit de réalignement temporel prévu en sortie du dernier dispositif, pour le dernier de ces signaux. Dans la deuxième de ces variantes, l'étage de som- nation numérique délivre les signaux Xm,j, Ym,j, Em,j suivants (Ki, Hi, Gi, étant des coefficients de pondé- ration): i=p Xm,j = E Ki (Mij-moj) i=1 i=p ym,j = i1 Hi (Mij-moj i=1 i=p Em,j = [ Gi Mij i=1 Ces signaux sont également fournis soit direc- tement en sortie des premier et deuxième dispositifs de somma- tion pondérée numérique, pour les deux premiers de ces si- gnaux, soit par l'intermédiaire d'un circuit de réalignement temporel prévu en sortie du dernier dispositif, pour le der- nier de ces signaux. Dans l'une ou l'autre de ces structures de ca- méra, il est possible, grâce à l'utilisation du seuil asservi - 10 - à la fois à l'énergie des événements et à l'effet d'empilement desdits événements, d'obtenir une résolution spatiale optimale quel que soit le taux de comptage et dans toute la gamme d'énergie des rayonnements utilisés. La caméra à scintillation ainsi obtenue est plus rapide et plus performante que les réalisations antérieures et, entre autres, convient tout particulièrement pour les examens multi-isotopiques. Les particularités et avantages de l'invention apparaîtront maitenant de façon plus détaillée dans la des- cription qui suit et dans les dessins annexés, donnés à titre d'exemples non limitatifs et dans lesquels: - les figures la à lc montrent respectivement la forme des signaux individuels correspondant à des scintilla- tions rapprochées entraînant un empilement partiel (figure la), la forme du signal global résultant de cet empilement (figure lb), et la forme du signal représentant pour chaque événement j la valeur mesurée à l'instant to,j+1 et qui ré-- sulte de la sommation des échantillons pendant l'intervalle de temps $j,j+ l pour la voie d'acquisition i (figure lc); - les figures 2 et 3 montrent les schémas de prin- cipe d'une caméra à scintillation selon l'état de la technique et selon l'invention respectivement; - la figure 4 montre un premier mode de réalisa- tion du calculateur d'une caméra à scintillation conforme à l'invention; la figure 5 montre une exemple de réalisation d'un des circuits de calcul de désempilement de l'étage-de traitement des événements; - la figure 6 montre un exemple de réalisation de l'étage de détection, séquencement et stockage; - la figure 7 montre une variante de réalisation des circuits de calcul de désempilement tels que celui de la figure 5; - 11 - - la figure 8 montre un deuxième mode de réalisa- tion du calculateur d'une caméra à scintillation conforme à l'invention, incluant la variante de réalisation du circuit de calcul de désempilement représentée sur la figure 7; - la figure 9 montre un troisième mode de réalisa- tion du calculateur d'une caméra à scintillation conforme à l'invention, incluant un troisième mode de réalisation de l'étage de traitement des événements; - les figures 10 à 12 montrent, en correspondance aux figures 4, 8 et 9, les modifications du calculateur lors- qu'on n'utilise plus que trois voies de calcul X, Y, Z; - les figures 13 à 15 montrent de même, en corres- pondance aux figures 4, 8 et 9, les modifications du calcula- teur lorsqu'on n'utilise plus que trois voies de calcul X, Y, E. Dans les caméras à scintillation munies de p voies d'acquisition utilisant des circuits convertisseurs-intégra- teurs tels que le dispositif décrit dans la demande de brevet français FRA- 2 540 995 ou celui décrit dans la demande de brevet français FR-A- 2 552 233, les signaux obtenus en sortie des voies d'acquisition sont une mesure de la charge des im- pulsions électriques fournies par les tubes photonultiplica- teurs respectifs en réponse à une scintillation. Par exemple, en l'absence d'effet d'empilement, la voie d'acquisition i (i variant de 1 à p) fournirait un signal Sij en réponse à une scintillation j survenant en un point xj, yj du cristal scintillateur. S'il y a au contraire em- pilement d'événements comme l'illustre à titre d'exemple la figure la qui montre des signaux individuels correspondant à des scintillations rapprochées entraînant un empilement par- tiel, le signal Mij (voir la figure lc qui montre la forme des signaux Mij successifs), fourni par la voie d'acquisi- tion i pour un événement j et mesuré à la fin de l'intervalle de temps 8j,j+1 = to,j+l-to,j pendant lequel cet événe- ment j n'est perturbé ni par les événements précédents ni par - 12 - les événements suivants, n'est plus Sij mais un signal fonction de Sij et des valeurs Sik (k = j-1, j-2,..., j-q+1,j-q) correspondant aux q événements précédents qui par- ticipent aux perturbations de l'événement j. Plus précisément, on peut montrer que, pour cha- cune des voies d'acquisition, le signal mesuré Mij est de la forme théorique suivante: Sij k j-q Mij = =S j-+ Ykj Sik (16) ai k=j- Dans cette relation (16), les coefficients ajet %kj sont des coefficients de correction, prédéterminés à partir de la connaissance de la forme moyenne, en fonction du temps, des impulsions analogiques fournies à l'entrée des voies d'acqui- sition. Les coefficients aj, obtenus par extrapolation, sont seulement fonction de l'intervalle de temps 8j,j+1 de la-me- sure. Ces coefficients aj, supérieurs à 1, tendraient donc vers la valeur 1 si la durée de l'intervalle de mesure Gj,j+1 devenait très grande. Les coefficients Ykj, obtenus par interpolation, sont fonction de l'intervalle de temps ekj séparant chaque événement précédent j-1, j-2, etc... de l'événement j et exprimé par la relation Gkj = to,j-to,k avec k = j-1, j-2,..., q-1, q. Ces coefficients Ykj, inférieurs à 1, tendraient vers la valeur 0 si la durée des intervalles Gkj devenait très grande. Sur la figure la, la charge électrique à détermi- ner, qui serait notée Qj et qui est proportionnelle à la quantité de courant débité pendant la durée ej,j+1 du signal individuel associé à l'événement j, est estimée à partir de la mesure Mij effectuée pendant cette durée et à partir des q charges Qk précédemment estimées. La figure lb montre le signal somme, ou signal composite, résultant de l'empilement - 13 - des signaux individuels (de la figure la) correspondant à plu- sieurs scintillations rapprochées dans le temps. La figure lc montre, on l'a vu, la forme des signaux Mi,j successifs. Lorsqu'on considère le cas de l'absence d'empile- ment, le calculateur présent dans la caméra, qui permet no- tamment, on l'a vu, d'évaluer l'énergie Ej de chaque événe- ment j à partir des signaux fournis par les p voies d'acquisi- tion, effectue ce calcul à l'aide de la somme pondérée sui- vante (relation (3) déjà nommée): i=p Ej = 1 Gi Sij i=1 S'il y a au contraire empilement d'événements, la quantité Em,j est donnée par: i=p Emj = 1 Gi Mij (17) i=1 ou encore, d'après la relation (16): Ej k =j-q Eij = + E Nkj Ek (18) aj kj-1 Pour l'évaluation des coordonnées xj, yj de chaque évé- nement j par le calculateur à partir des signaux fournis par les p voies d'acquisition, et ce indépendamment du fait qu'il y ait ou non un empilement d'événements, on introduit donc, conformément à l'invention, un seuil noté moj et variable d'un événement à l'autre. De la sorte, le calculateur réalise dans un pre- mier temps les sommes pondérées suivantes: i=p Xmj = E Ki(Mij-moj) (19) i=1 i=p ZmJ = i JiMij (21) i = i=( - 14 - avec Mij-moj = 0 si Mij est inférieur ou égal à moj. Dans ces relations (19) à (21) les facteurs de pondération sont par exemple les mêmes que précédemment. On peut montrer que la valeur du seuil moj est avantageusement -sans bien entendu que ce soit le seul choix possible- prise proportionnelle à la valeur Emj obtenue à l'aide des relations (17) et (18). Le facteur de proportiona- lité est appelé fo et est exprimé en V/eV si moj est une tension et si Emj est exprimé en électron-volts, et l'on a donc: moj = fo.Emj Ej k = j- q moj = fo À ( Ej + E = Ykj Ek) (22) 0j k = j-1 = kE1 De cette expression, il résulte qu'en cas d'absence d'empile- ment correspondant aux valeurs aj = 1 et Ikj = 0, le seuil, qualifié d'asservi, prend la valeur soj = fo.Ej, et les signaux fournis par les p voies d'acquisition sont les signaux Sij. Dans ces conditions, (absence d'empilement), les relations (19), (20), (21) s'-écrivent maintenant: i=p Xj = 1 Ki (Sij-soj) (23) 1=1 :l=p Yj = i Hi (Sij-Soj) (24) i=1 1=p Zj= jÉ Ji Sij (25) Les sommes pondérées ne sont pas entachées d'erreurs dues à l'empilement, et le calcul des coordonnées selon les relations (1) et (2) est immédiat. - 15 - Lorqu'au contraire il y a empilement des événe- ments, compte tenu des expressions (16), (22) et (23) à (25) les sommes pondérées (19) à (21) peuvent maintenant être ex- primés de la façon suivante: X k = j-q Xm j = -+E kj Xk (26) aj k=j-1 Y- k =j-q Ymj = kj +Ykk (27) raj k=- Z k = j-q Zmj Zj- + E N- kj Zk (28) =j-1 Après cette opération de sommation pondérée réa- lisée dans un premier temps, une deuxième opération consiste à effectuer un traitement des événements empilés de façon à dé- terminer, à partir des sommes pondérées (26) à (28), les véri- tables coordonnées xj, yj de chaque scintillation. Ce traitement est réalisable selon différents modes de réalisa- tion qui vont être décrits plus loin, au cours de la descrip- tion d'exemples de caméra à scintillation selon l'invention. On rappellera préalablement la structure d'une ca- méra à scintillation de type classique. Une telle caméra, re- présentée sur la figure 2, comprend un cristal scintillateur équipé d'un collimateur 20 et destiné à convertir chaque photon reçu en une scintillation. Ce cristal est couplé par l'intermédiaire d'un guide de lumière 30 à la fenêtre d'entrée d'un jeu de p photodétecteurs constitué ici par des tubes pho- tomultiplicateurs 50. Ces tubes 50 convertissent chaque scin- tillation en un courant qui est alors traité par p voies d'ac- quisition 60 entièrement analogiques. Ces voies d'acquisition réalisent notamment l'amplification, le filtrage, l'inté- gration et la mise en forme des signaux de sortie des tubes photomultiplicateurs 50, et sont suivies d'un calculateur 100 fournissant les coordonnées xj, Yj et l'énergie Ej. - 16 - Selon le mode de réalisation qui va être mainte- nant plus particulièrement décrit, en correspondance à la fi- gure 3 montrant par rapport à la figure 2 les modifications du schéma de principe pour une caméra selon l'invention, les p voies d'acquisition 60 ne sont plus totalement analogiques comme dans la caméra classique décrite. Ces voies fournissent maintenant p signaux numériques Mi,j (i = indice variant de 1 à p) à l'entrée du calculateur 100, et réalisent maintenant successivement l'amplification, le filtrage et l'échantillon- nage des signaux de sortie des tubes photomultiplicateurs 50 puis la conversion analogique-numérique des échantillons obte- nus et la sommation de ces échantillons numériques. La valeur de ces p signaux numériques est liée à celle du courant de sortie des tubes 50 et donc à une fraction de l'intensité de la scintillation initiale, mais différemment selon le taux d'empilement des scintillations.(cette fraction est elle-même liée à la réalisation du bloc optique et en particulier à la distance entre le point de scintillation et l'axe des tubes). S'il n'y avait pas d'empilement, la valeur de chacun de ces signaux serait notée Si,j; l'estimation de ces valeurs en présence d'empilement sera notée Sij. Chacune de ces p voies comprend en série, pour réaliser les fonctions citées ci-dessus, un circuit 61 d'am- plification et de filtrage recevant la sortie du tube 50 cor- respondant, un circuit de réalignement temporel 62, puis un dispositif de conversion et d'intégration 63 assurant succes- sivement l'échantillonnage des signaux de sortie du circuit 62 correspondant, la conversion analogique-numérique des échan- tillons obtenus et la sommation de ceux-ci. Il est en outre prévu un amplificateur analogique de sommation 64 dont les p entrées reçoivent les p sorties des circuits d'amplification et de filtrage 61 et dont la sortie est fournie à un détecteur de début d'impulsion localisé dans l'étage de détection, sé- quencement et stockage 400 décrit plus loin. - 17 - La sortie de chacun des p dispositifs de conver- * sion et d'intégration est envoyée vers le calculateur 100, éventuellement par l'intermédiaire de p mémoires dites FIFO (de l'anglais "First-In, First-Out") permettant de travailler OS ultérieurement à des fréquences plus faibles en régularisant le débit des événements. Cette disposition de mémoires FIFO permet en effet une réduction de la vitesse des calculs ulté- rieurs et le rythme plus lent ainsi obtenu peut être pratique- ment égal au rythme moyen d'arrivée des événements (par exem- ple 2 microsecondes pour un rythme moyen de 500.000 événements par seconde) et non plus au rythme aléatoire d'arrivée des événements (environ 0,2 microseconde dans l'exemple précé- dent). Chacun des dispositifs de conversion et d'intégration 63 est ici équivalent à celui décrit dans la demande de brevet français FR-A-2 552 233, et les opérations effectuées par ces dispositifs sont contrôlées par l'étage de détection, séquen- cement et stockage 400. Le calculateur 100 recevant les p sorties des voies d'acquisition comprendlui-même différents dispositifs de calcul qui vont permettre de déterminer distinctement les coordonnées xj, yj et l'énergie Ej de toute scintilla- - tion j, par exemple à l'aide des relations (1) à (6) dans le cas d'un dispositif de calcul barycentrique à rapport arithmétique. Plus précisément, ce calculateur 100, représenté sur la figure 4 dans le cas d'un dispositif de calcul barycen- trique à rapport arithmétique, est agencé de la façon sui- vante. Il comprend tout d'abord un bus de transfert 150 de. signaux numériques Mi,j disponibles en sortie des p voies d'acquisition. Pour l'une de ces voies, par exemple la voie i, on a représenté les signaux analogiques individuels associés à plusieurs scintillations rapprochées dans le temps (figure la) et le signal résultant de l'empilement de ces signaux individuels (figure lb), l'événement j étant perturbé en amont par plusieurs événements j-1, j-2, etc... Si ej et lk,j - 18 - sont les coefficients de correction respectivement par extrapolation et par interpolation qui, on l'a vu, peuvent être déterminés à partir de la connaissance de la forme moyenne, en fonction du temps, des impulsions correspondant à un événement détecté et à partir de la mesure de la durée ej,j+1 comprise entre toj et to,j+1, si k est suc- cessivement égal à j-1, j-2,..., j-q, et si les signaux Sik représentent pour ces valeurs respectives de k les valeurs corrigées des effets d'empilement-qui seraient fournies par la voie d'acquisition numérique i en réponse aux scintillations j-1, j-2,..., j-q, alors les signaux numériques obtenus en sortie des voies d'acquisition sont de la forme: ki +kj- q Mij = j+ kj-q1 Nkj Si,k (29) ai kj- o Mi,j représente la valeur mesurée à l'instant to,j+1 résultant de la sommation des échantillons pendant l'inter- valle de temps Gj,j+1 pour la voie i. Le calculateur comprend ensuite, en sortie du bus de transfert 150, un étage de sommation numérique 200. Comme indiqué sur la figure 4, qui montre un premier exemple de réalisation du calculateur selon l'invention, cet étage de sommation numérique 200 comprend tout d'abord deux dispositifs de sommation pondérée numérique 223 et 224 qui réalisent les sommations pondérées suivantes: i=p Zm,j = i Ji Mij (30) i=1 i=p 25.Emj = 1 Gi Mij (31) î==i La grandeur Em,j est utilisée par un calculateur de seuil 231 qui peut être un multiplicateur de cette grandeur par le facteur fo (déterminé expérimentalement et mémorisé) et qui délivre sur sa sortie la valeur du seuil asservi moj. - 19 - L'étage 200 comprend également une mémoire tampon 201 qui reçoit les p signaux de sortie du bus de transfert 150 et est destinée à retarder ces p signaux pour laisser s'effec- tuer préalablement les opérations des circuits 224 et 231. La commande de lecture de la mémoire 201 est fournie, comme indi- qué sur la figure 4, par le calculateur de seuil 231, et la sortie de cette mémoire 201 est alors envoyée vers un circuit soustracteur de seuil 211 qui reçoit, sur une première entrée, les p signaux Mij retardés par la mémoire tampon 201 et, sur une deuxième entrée, la valeur du seuil asservi moj. Les si- gnaux de sortie du soustracteur sont les signaux Mij-moj. L'étage 200 comprend enfin deux autres dispositifs de sommation pondérée numérique 221 et 222 qui réalisent res- pectivement les sommations pondérées suivantes: i=p Xm,j = 1 Ki(Mij-moj) (32) i=1 i=p Ym,j = i Hi(Mij-moj) (33) i=1 Chacun de ces quatre dispositifs de sommation pon- dérée numérique 221 à 224 est par exemple du type du multi- plieur-accumulateur TDC 1009 (commercialisé par la société TRW, La Jolla, CA 92038, USA), dont l'une des entrées reçoit la sortie correspondante Mij ou Mij-moj selon le cas, et dont l'autre reçoit les coefficients de pondération sous forme numérique, stockés dans une mémoire auxiliaire. Dans le cas o l'on utilise effectivement ce type de multiplieur-accumula- teur, ladite mémoire auxiliaire, qui doit être synchronisée avec le fonctionnement d'ensemble du calculateur, peut être par exemple incorporée à l'étage de détection, séquencement et stockage 400 décrit plus loin. Les signaux de sortie Xmj, Ymj, Zmj, Emj disponibles en sortie de l'étage de sommation numérique 200, - 20 - provenant des quatre dispositifs de sommation pondérée numé- rique 221 à 224 soit directement pour les deux premiers de ces signaux, soit par l'intermédiaire de circuits de réalignement temporel 233 et 234 pour les deux autres, sont alors fournis à un étage de traitement des événements 500. Cet étage 500 comprend, comme indiqué sur la figure 4, quatre circuits de calcul de désempilement 501 à 504, deux diviseurs 505 et 506 et un circuit de réalignement temporel 507. Les quatre circuits 501 à 504 étant identiques, on ne décrit ici que l'un d'eux, par exemple le circuit 501. Ce circuit, représenté sur la figure 5, comprend un soustracteur 510 recevant sur sa première entrée positive la sortie du dispositif de sommation pondérée numérique correspondant 221 (aux dispositifs 221 à 224 correspondent respectivement les circuits 501 à 504). Le soustracteur 510 est suivi d'un premier multiplieur 511 et d'un registre de stockage 512 dont la sortie est celle du circuit 501. Le soustracteur 510 est également suivi, en parallèle sur les éléments 511 et 512, d'un deuxième multiplieur 513 et d'un deuxième registre de stockage 514. Ces multiplieurs peuvent être remplacés par un seul circuit de multiplication associé à un multiplexeur-démultiplexeur temporel. L'entrée négative du soustracteur 510 est reliée à la sortie du registre de stockage 514. La deuxième entrée du multiplieur 511 est reliée à la sortie d'une mémoire 470 stockant le coefficient aj et celle du multiplieur 513 est reliée à la sortie d'une mémoire 480 stockant le coefficient Nj,k. Les sorties X, Y, Z, E des quatre circuits de calcul de désempilement 501 à 504 sont données, pour l'événement j, par les expressions: k=j-q Xj = aj [Xm,j - k=j-1 k,j Xk] (34) - 21 - k=j-q Yj = [Y,j - k__j Yk] (35) yj = 0aj [Ymi - k=j-1 Yk,j - k=j-q Zj = aj [Zm,j - kEj-1 [j Zk (36) k=-1kZ] - k=j-q Ej = aj[Em,j - kEj- kj Ek (37) k=j-l kj k Les éléments de chaque circuit 501 à 504, par exemple les éléments 510 à 514 du circuit 501, constituent un circuit de calcul de désempilement équivalent à celui qui est décrit dans la demande FR-A-2 552 233 avec les références 120 à 160. Les trois autres circuits 502 à 504 comprennent les mêmes éléments que le circuit 501. La sortie X du circuit de calcul de désempilement 501 est envoyée vers la première entrée du diviseur 505 et la sortie Y du circuit 502 vers la première entrée du diviseur 506. La deuxième entrée de chacun de ces diviseurs est consti- tuée par la sortie Z du circuit de calcul de désempilement 503. Les trois sorties de l'étage de traitement, qui sont aussi celles du calculateur, sont constituées par la sor- tie xj = Xj/Zj du diviseur 505, la sortie xj = Yj/Zj du diviseur 506 et la sortie Ej du circuit de réalignement temporel 507 prévu en sortie du circuit de calcul de désempi- lement 504. A ces éléments est encore associé, dans le calcu- lateur, l'étage de détection, séquencement et stockage 400. Cet étage 400, représenté sur la figure 6, comprend d'abord un détecteur de début d'impulsion 410, qui reçoit la sortie de l'amplificateur analogique de sommation 64 (voir les figures 3 et 4). Ce détecteur 410 est suivi d'un circuit d'horloge 420 qui génère des signaux périodiques assurant la commande - 22 - d'un circuit de séquencement 450 et celle d'un compteur 430 de ces signaux d'horloge. Le nombre ainsi compté est soumis à un circuit de test 440 dont la sortie est envoyée vers le circuit de séquencement 450. Ce dernier assure la synchronisation des opérations effectuées dans les voies d'acquisition, dans l'étage 200 et dans l'étage 500, et valide le contenu d'un registre 460 de stockage de la sortie du compteur 430, registre qui est en parallèle sur le circuit de test 440. En sortie du registre 460 sont prévues les deux mémoires 470 et 480 mentionnées plus haut en référence à la figure 5 et stockant respectivement le coefficient aj et le coefficient lj,k. Les éléments 410 à 480 qui viennent d'être cités constituent un étage de détection, séquencement et stockage similaire à celui qui est décrit dans la demande FR-A-2 552 233. Dans le premier mode de réalisation qui vient d'être décrit, la structure du calculateur a permis de déduire A les sommes pondérées exactes Xj, Yj, Zj d'une part des valeurs calculées Xm,j, Ym,j, Zm,j et d'autre part des coefficients aj et 'k,j, en utilisant les formules (34) à (37). Dans un deuxième mode de réalisation du calculateur, le calcul des coordonnées des événements peut être effectué sans introduire d'extrapolation, selon les expressions suivantes: k=j-q X'j Xm,j - k=j-1 k,j X'k (38) k=j-q Y'j = Ym,j - kE Y-1'k (39) k=j-l k,j k=j-q Z' j = Xm, -j - ú Zk (40) k=-1k, - 23 - avec toujours Xj = X'j/Z'j et yj = Y'j/Z'j. Dans ce mode de réalisation, l'étage de traitement des événements est maintenant référencé 600, et les trois circuits de calcul de désempilement qui reçoivent respectivement les signaux Xm, Ym, Zm, sont, par rapport à la figure 5, modifiés par suppression du multiplieur 511. Ces trois circuits ont en effet la structure re- présentée sur la figure 7 pour l'un quelconque, par exemple le premier de ces trois circuits 601 à 603. Ce circuit 601 com- prend un soustracteur 610 qui reçoit sur sa première entrée la sortie du dispositif de sommation pondérée numérique corres- pondant 221. Le soustracteur 610 est suivi d'une part direc- tement d'un registre de stockage 612, dont la sortie est comme précédemment celle du circuit de calcul de désempilement, et d'autre part, en parallèle, d'un multiplieur 613 puis d'un re- gistre de stockage 614. La sortie de ce registre 614 est ren- voyée vers la deuxième entrée du soustracteur 610, et l'autre entrée du multiplieur 613 est reliée à la sortie de la mémoire 480 prévue dans l'étage 400 pour le stockage du coefficient y. Les deux autres circuits 602 et 603 comprennent des éléments similaires. Le circuit de calcul de désempilement 504 restant inchangé, le calculateur 100 prend maintenant la configuration représentée sur la figure 8. L'énergie E est, comme précédemment, disponible en sortie du circuit de réali- gnement temporel 507. Dans le deuxième mode de réalisation qui vient d'être décrit, la structure du calculateur a permis cette fois de procéder d'abord par détermination de sommes pondérées sans extrapolation puis par calcul des coordonnées comme précédem- ment. Dans un troisième mode de réalisation, une autre va- riante de calculateur est proposée, dans laquelle on détermine d'abord des coordonnées xmj, Ymj dites mesurées, qui sont entachées d'erreurs du fait de l'empilement des événements: xmj = Xm, j/Zm,j (41) Ymj = Ym,j/Zm, j (42) 2605-745 - 24 - A cette détermination succède alors le calcul des coordonnées exactes de l'événement j, au moyen des expressions suivantes: k=j-q xj = xmj - ú r (43) k=j-1 kj (Xk - Xmj) k=j-q Yj = Ymj - k=j-1 kj (Yk - Ymj) (44) Plus précisément, dans ce troisième mode de réali- sation de l'invention, représenté sur la figure 9, le calcula- teur 100 comprend un troisième type d'étage de traitement des événements, référencé 700. Dans cette variante de réalisation, le traitement de désempilement est effectué non plus sur les signaux Xm, Ym, Zm de sortie de l'étage de sommation nu- mérique 200 mais sur les coordonnées xmj et Ymj obtenues en sortie de deux diviseurs 705 et 706. Ce diviseur reçoit les signaux Xm et Zm de sortie de l'étage de sommation numé- rique 200, et le diviseur 706 reçoit les signaux Ym et Zm. Les coordonnées xj, yj correspondant à l'événement j sont alors obtenues en sortie de circuits de calcul de désem- pilement 701 et 702, à partir d'une part de ces coordonnées dites mesurées et d'autre part des coordonnées déjà connues Xk, Yk des événements précédents qui perturbent l'événe- ment j, conformément aux expressions (43) et (44). Les coefficients rk,j, qui sont fonctions des intervalles de temps mesurés ej,j+1, Ok,j et des rapports Ek/Ej sont calculés dans un circuit de calcul additionnel 707 qui reçoit d'une part la sortie E du circuit de calcul de désempilement 504, toujours présent dans l'étage de traite- ment, et d'autre part les coefficients a et y délivrés par l'étage de détection, séquencement et stockage 400. Dans cet exemple, le circuit 707 calcule d'une part, à partir des valeurs de E successivement reçues, les rapports successifs Ek/Ej et d'autre part les produits aj Nk,j à partir desquels sont évalués les coefficients rk,j selon la relation r k,j = "j yk,j Ek/Ej. Ces coefficients - 25 - rkj sont fournis aux circuits de calcul de désempilement 701 et 702, qui ont une configuration similaire à celle des circuits 601 et 602. L'énergie E est, là encore comme précé- demment, disponible en sortie du circuit de réalignement tem- porel 507. Bien entendu, la présente invention n'est pas li- mitée aux exemples de réalisation décrits et représentés, à partir desquels d'autres variantes peuvent encore être propo- sées sans pour cela sortir du cadre de l'invention. Par exemple, dans les trois modes de réalisation décrits, le calcul du seuil dans le circuit 231 peut être ef- fectué à partir du signal Zm en lieu et place du signal Em. Dans ce cas, la connexion d'entrée du calculateur de seuil 231 est alors reliée non plus à la sortie du dispositif 224 mais à celle du dispositif 223. D'autre part, un circuit de réjection en amplitude du nombre des événements à traiter peut être prévu, pour n'effectuer les calculs que sur des évé- nements sélectionnés (à l'aide d'un seuil, d'une fenêtre d'énergie, etc...). Par ailleurs, un multiplexeur temporel peut être prévu pour n'utiliser qu'un diviseur au lieu de deux dans cha- cun des modes de réalisation des figures 4, 8 et 9. Un circuit de multiplexage temporel peut également être prévu pour ré- duire le nombre des circuits de calcul de désempilement et n'utiliser qu'un seul circuit de calcul de désempilement au lieu des quatre circuits 501 à 504 dans le cas du mode de réa- lisation de la figure 4, des quatre circuits 601 à 603 et 504 dans le cas du mode de réalisation de la figure 8, ou des trois circuits 701, 702 et 504 dans le cas du mode de réalisa- tion de la figure 9. L'ensemble des variantes proposées dans ce paragraphe s'applique également aux autres modes de réalisation des figures 10 à 15 mentionnés ci-dessous. Compte tenu des moyens de correction des défauts de linéarité et d'énergie qui sont mis en oeuvre dans les caméras gamma actuellement disponibles à partir des signaux x, - 26 - y, E obtenus en sortie du calculateur, on sait en effet qu'il - est possible d'utiliser indifféremment Z ou E pour le calcul des coordonnées. Dans ce cas on peut ne calculer qu'une seule de ces deux grandeurs et, selon le choix effectué, déduire l'autre grandeur de calculs faisant intervenir des corrections spécifiques de ce choix et exécutées par lesdits moyens. L'étage de sommation numérique 200 ne comprend plus alors que trois dispositifs de sommation pondérée numérique fournissant des signaux Xm, Ym, Zm ou Xm, Ym, Em respectivement, et de même l'étage de traitement des événe- ments ne comprend plus que trois circuits de calcul de désem- pilement. Les figures 10 à 12 montrent, en correspondance aux figures 4, 8, 9, les modifications du calculateur lorsqu'on n'utilise plus que trois voies X, Y, Z, tandis que les fi- gures 13 à 15 montrent, toujours en correspondance aux fi- gures 4, 8, 9, les modifications du calculateur lorsqu'on n'utilise plus que trois voies X, Y, E. D'autre part, on notera également qu'il est possible, pour travailler à des fréquences plus faibles en régularisant le débit des événements, d'introduire en amont de l'étage de sommation numérique des mémoires FIFO à écri- ture et lecture commandées par l'étage de détection, séquen- cement et stockage 400. Enfin, il est équivalent de considérer que le bus 150 est incorporé au calculateur, dont il est l'élément d'en- trée ou d'accès, ou qu'au contraire il lui est associé sans lui être inclus. - 27 - REVENDICATIONS: 1. Caméra à scintillation comprenant un cristal scin- tillateur éventuellement équipé d'un collimateur et destiné à convertir chaque photon reçu en une scintillation, un guide de lumière pour le couplage dudit cristal à la fenêtre d'entrée d'un jeu de p photodétecteurs destinés à convertir en courant chaque scintillation, p voies d'acquisition recevant les si- gnaux de sortie desdits photodétecteurs et délivrant p signaux électriques de caractéristiques liées notamment à l'intensité de la scintillation et à la distance de cette scintillation à chacun des photodétecteurs, et un calculateur destiné à déli- vrer les coordonnées xj et yj d'une scintillation j et l'énergie Ej qui est associée à cet événement j, caracté- risée: (A) en ce que lesdites p voies d'acquisition réa- lisent l'amplification, le filtrage et l'échantillonnage des- dits signaux de sortie des photodétecteurs, puis la conversion analogique-numérique des échantillons obtenus et leur somma- tion, et délivrent p signaux numériques à l'entrée du calcu- lateur; (B) en ce que le calculateur comprend lui-même: (a) un bus de transfert desdits p signaux numériques; (b) un étage de sommation numérique comprenant lui-même: - une mémoire-tampon recevant les p signaux de sortie Mij du bus de transfert; - un circuit soustracteur de seuil recevant sur une première entrée lesdits p signaux, retardés par la mé- moire-tampon, et sur une deuxième entrée une valeur moj de seuil dit asservi et délivrant p signaux Mij-moj, ces signaux prenant conventionnellement la valeur 0 si Mij est inférieur à moj; - quatre dispositifs de sommation pondérée numé- rique dont les deux premiers reçoivent les p signaux - 28 - Mij-moj délivrés par le circuit soustracteur de seuil et dont les deux autres reçoivent directement les p signaux de sortie Mij du bus de transfert; - un calculateur de seuil recevant la sortie du quatrième dispositif de sommation pondérée numérique pour l'évaluation du seuil asservi moj selon une relation de pro- portionalité à la valeur Emj associée à chaque événement j; ledit étage de sommation numérique délivrant les signaux Xm,j, Ym,j, Zm,j, Em,j suivants (Ki, Hi, Ji, Gi, étant des coefficients de pondération): i=p Xm, j = E Ki (Mij-moj) i= i= i=p Ym' ij = E Hi (Mij-moj) i=1 i=p Zmj = 1 Ji Mij i=1 i=p Emj = E Gi Mij i=1 soit directement en sortie des deux premiers dispositifs de sommation pondérée numérique, pour les deux premiers de ces signaux, soit par l'intermédiaire de deux circuits de réa- lignement temporel prévus respectivement en sortie des deux autres dispositifs, pour les deux derniers de ces signaux; (c) un étage de traitement des événements in- cluant des circuits de calcul de désempilement et deux divi- seurs et délivrant à partir des signaux Xmj, Ymj, Zmj, Emj les trois signaux de coordonnées et d'énergie xj, Yj, Ej relatifs à l'événement j; (C) en ce qu'un étage de détection, séquencement et stockage recevant un signal qui correspond à la somme des p signaux de sortie des photodétecteurs est prévu pour délivrer d'une part les différents signaux d'horloge pour la synchro- nisation des éléments des p voies d'acquisition et des éléments - 29 - du calculateur et d'autre part des coefficients de correction destinés à l'étage de traitement des événements. 2. Caméra à scintillation selon la revendication 1, caractérisée en ce que l'étage de traitement des événements comprend quatre circuits de calcul de désempilement, recevant pour chaque événement j respectivement les signaux numériques Xmj, Ymj, Zmj, Emj, et délivrant quatre signaux Xj, Yj, Zj, Ej, et deux diviseurs délivrant respectivement deux signaux xj = Xj/Zj et yj = Yj/Zj, les trois signaux xj, yj, Ej étant constitués par les signaux de sortie respec- tivement du premier diviseur, du deuxième diviseur et d'un circuit de réalignement temporel recevant la sortie du qua- trième circuit de calcul de désempilement, et en ce que les coefficients de correction sont des coefficients de correction par extrapolation et interpolation, a et y respectivement, des- tinés auxdits quatre circuits de calcul de désempilement. 3. Caméra à scintillation selon la revendication 2, caractérisée en ce que chacun des quatre circuits de calcul de désempilement comprend un soustracteur recevant sur sa pre- mière entrée la sortie correspondante de l'étage de sommation numérique, ce soustracteur étant suivi d'une part d'un premier multiplieur et d'un premier registre de stockage, et d'autre part, en parallèle sur ce premier multiplieur et ce premier registre, d'un deuxième multiplieur et d'un deuxième registre de stockage, la sortie de ce deuxième registre étant reliée à la deuxième entrée du soustracteur et les deuxièmes entrées des premier et deuxième multiplieurs étant reliées la pre- mière à la sortie de la mémoire de stockage du coeffi- cient a et la seconde à la sortie de la mémoire de stockage du coefficient y. 4. Caméra à scintillation selon la revendication 3, caractérisée en ce que les multiplieurs sont remplacés par un seul circuit de multiplication associé à un multiplexeur- démultiplexeur temporel. - 30 - 5. Caméra à scintillation selon la revendication 2, caractérisée en ce que chacun des trois premiers circuits de calcul de désempilement comprend un soustracteur recevant sur sa première entrée la sortie correspondante de l'étage de som- mation, suivi d'une part d'un troisième registre de stockage dont la sortie est celle du circuit de calcul de désempilement et d'autre part, en parallèle, d'un troisième multiplieur puis d'un quatrième registre de stockage, la sortie de ce quatrième registre étant reliée à la deuxième entrée du soustracteur et l'autre entrée du troisième multiplieur étant reliée à la sor- tie de la mémoire de stockage du coefficient w. 6. Caméra à scintillation selon la revendication 1, caractérisée en ce que l'étage de traitement des événements comprend trois circuits de calcul de désempilement, deux divi- seurs, un circuit de réalignement temporel, et un circuit de calcul additionnel, les deux diviseurs recevant les signaux Xmjl Ymjl Zmj pour délivrer deux signaux xmj = Xmj/Zmj et Ymj = Ymj/Zmj, les premier et deuxième circuits de calcul de désempilement recevant lesdits signaux xmj, Ymj et délivrant les signaux xj, yj et le troisième recevant le signal Emj et délivrant le signal E, et le circuit de calcul additionnel recevant d'une part ledit signal E et d'autre part lesdits coefficients de correction pour fournir aux premier et deuxième circuits de calcul de désempilement un coefficient de correction additionnel r. 7. Caméra à scintillation comprenant un cristal scintillateur éventuellement équipé d'un collimateur et destiné à convertir chaque photon reçu en une scintil_ lation, un guide de lumière pour le couplage dudit cristal à la fenêtre d'entrée d'un jeu de p photodé- tecteurs destinés à convertir en courant chaque scin- tillation, p voies d'acquisition recevant les signaux de sortie desdits photodétecteurs et délivrant p signaux - 31 - électriques de caractéristiques liées notamment à l'intensité de la scintillation et à la distance de cette scintillation à chacun des photodétecteurs, et un calculateur destiné à déli- vrer les coordonnées Xj et yj d'une scintillation j et l'énergie Ej qui est associée à cet événement j, caractérisée: (A) en ce que lesdites p voies d'acquisition réalisent l'amplification, le filtrage et l'échantillonnage desdits signaux de sortie des photodétecteurs, puis la conver- sion analogique-numérique des échantillons obtenus et leur sommation, et délivrent p signaux numériques à l'entrée du calculateur; (B) en ce que le calculateur comprend lui-même: (a) un bus de transfert desdits p signaux numériques; (b) un étage de sommation numérique compre- nant lui-même: - une mémoire-tampon recevant les p signaux de sortie Mij du bus de transfert; - un circuit soustracteur de seuil recevant sur une première entrée lesdits p signaux, retardés par la mé- moire-tampon, et sur une deuxième entrée une valeur moj de seuil dit asservi et délivrant p signaux Mij-moj, ces signaux prenant conventionnellement la valeur 0 si Mij est inférieur à moj; - trois dispositifs de sommation pondérée numé- rique dont les deux premiers reçoivent les p signaux Mij-moj délivrés par le circuit soustracteur de seuil et dont le dernier reçoit directement les p signaux de sortie Mij du bus de transfert; - un calculateur de seuil recevant la sortie dudit dernier dispositif de sommation pondérée numérique pour l'évaluation du seuil asservi moj selon une relation de proportionalité à la valeur Emj associée à chaque évé- nement j; - 32 - ledit étage de sommation numérique délivrant les signaux Xmj, Ym,j, Zm,j suivants (Ki, Hi, Ji, étant des coefficients de pondération): i=p Xmj E Ki (Sij-m J) j=p Ymj = iE Hi (Mij-moj) i=1 i=p Zmj = Ji Mij i=1i soit directement en sortie des premier et deuxième dispositifs de sommation pondérée numérique, pour les deux premiers de ces signaux, soit par l'intermédiaire d'un circuit de réaligne- ment temporel prévu en sortie du dernier dispositif, pour le dernier de ces signaux; (c) un étage de traitement des événements in- cluant des circuits de calcul de désempilement et deux divi- seurs et délivrant à partir des signaux Xmj, Ymj, Zmj, les trois signaux de coordonnées et d'énergie xj, yj, Ej relatifs à l'événement j; (C) en ce qu'un étage de détection, séquencement et stockage recevant un signal qui correspond à la somme des p signaux de sortie des photodétecteurs est prévu pour délivrer d'une part les différents signaux d'horloge pour la synchronisation des éléments des p voies d'acquisition et des éléments du calculateur et d'autre part des coefficients de correction destinés à l'étage de traitement des événements. 8. Caméra à scintillation comprenant un cristal scintillateur éventuellement équipé d'un collimateur et destiné à convertir chaque photon reçu en une scintillation, - 33 - un guide de lumière pour le couplage dudit cristal à la fenêtre d'entrée d'un jeu de p photodétecteurs destinés à convertir en courant chaque scintillation, p voies d'acqui- sition recevant les signaux de sortie desdits photodétec- os05 teurs et délivrant p signaux électriques de caractéristiques liées notamment à l'intensité de la scintillation et à la distance de cette scintillation à chacun des photodétecteurs, et un calculateur destiné à délivrer les coordonnées xj et Yj d'une scintillation j et l'énergie Ej qui est associée à cet événement j, caractérisée: (A) en ce que lesdites p voies d'acquisition réalisent l'amplification, le filtrage et l'échantillonnage desdits signaux de sortie des photodétecteurs, puis la conver- sion analogique-numérique des échantillons obtenus et leur sommation, et délivrent p signaux numériques à l'entrée du calculateur; (B) en ce que le calculateur comprend lui-même: (a) un bus de transfert desdits p signaux numériques; (b) un étage de sommation numérique comprenant lui- même: une mémoire-tampon recevant les p signaux de sortie Mij du bus de transfert; - un circuit soustracteur de seuil recevant sur une première entrée lesdits p signaux, retardés par la mé- moire-tampon, et sur une deuxième entrée une valeur moj de seuil dit asservi et délivrant p signaux Mij-moj, ces signaux prenant conventionnellement la valeur 0 si Mij est inférieur à moj; - trois dispositifs de sommation pondérée numéri- que dont les deux premiers reçoivent les p signaux Mij-moj délivrés par le circuit soustracteur de seuil et dont - 34 - le dernier reçoit directement les p signaux de sortie Mij du bus de transfert; - un calculateur de seuil recevant la sortie dudit dernier dispositif de sommation pondérée numérique pour l'évaluation du seuil asservi moj selon une relation de proportionalité à la valeur Emj associée à chaque évé- nement j; ledit étage de sommation numérique délivrant les signaux Xm,j, Ym, j, Em,j suivants (Ki, Hi, Gi, étant des coefficients de pondération): i=p Xm,j = Ki (Mij-moj) i=p Ymj = E Hi (Mij-moj Mi= i=p Emj = Gi Mij ji= soit directement en sortie des premier et deuxième dispositifs de sommation pondérée numérique, pour les deux premiers de ces signaux, soit par l'intermédiaire d'un circuit de réa- lignement temporel prévu en sortie du dernier dispositif, pour le dernier de ces signaux; (c) un étage de traitement des événements incluant des circuits de calcul de désempilement et deux diviseurs et délivrant à partir des signaux Xmj, Ymj, Emj les trois signaux de coordonnées et d'énergie Xj, yj, Ej relatifs à l'événement j; (C) en ce qu'un étage de détection, séquencement et stockage recevant un signal qui correspond à la somme des p signaux de sortie des photodétecteurs est prévu pour délivrer d'une part les différents signaux d'horloge pour la synchroni- sation des éléments des p voies d'acquisition et des éléments du calculateur et d'autre part des coefficients de correction destinés à l'étage de traitement des événements. - 35 - 9. Caméra à scintillation selon l'une des revendica- tions 7 et 8, caractérisée en ce que l'étage de traitement des événements comprend trois circuits de calcul de désempilement, recevant respectivement lesdits signaux numériques, Xmj, Ymj, Zmj, ou Xmj, Ymj,. Emj respectivement, et délivrant trois signaux Xj, Yj, Zj, ou Xj, Yj, Ej respectivement, et deux diviseurs délivrant respectivement deux signaux A A ^ A xj = Xj/Zj et yj = Yj/Zj, ou xj = Xj/Ej et yj = Yj/Ej respectivement, les trois signaux xj, yj, Ej étant cons- titués par les signaux de sortie respectivement du premier diviseur, du deuxième diviseur et d'un circuit de réaligne- ment temporel recevant la sortie du dernier circuit de calcul de désempilement, et en ce que les coefficients de correction sont des coefficients de correction par extrapolation et inter- polation, a et X respectivement, destinés aux trois circuits de calcul de désempilement de l'étage de traitement des évé- nements. 10. Caméra à scintillation selon la revendication 9, caractérisée en ce que chacun des trois circuits de calcul de désempilement comprend un soustracteur recevant sur sa pre- mière entrée la sortie correspondante de l'étage de sommation numérique, ce soustracteur étant suivi d'une part d'un premier multiplieur et d'un premier registre de stockage et d'autre part, en parallèle sur ce premier multiplieur et ce premier registre, d'un deuxième multiplieur et d'un deuxième registre de stockage, la sortie de ce deuxième registre étant reliée à la deuxième entrée du soustracteur et les deuxièmes entrées des premier et deuxième multiplieurs étant reliées la première à la sortie de la mémoire de stockage du coefficient a et la seconde à la sortie de la mémoire de stockage du coefficient. 11. Caméra à scintillation selon la revendication 10, caractérisée en ce que les multiplieurs sont remplacés par un seul circuit de multiplication associé à un multiplexeur-dé- multiplexeur temporel. - 36-- 12. Caméra à scintillation selon la revendication 9, caractérisée en ce que chacun des trois circuits de calcul de désempilement comprend un soustracteur recevant sur sa pre- mière entrée la sortie correspondante de l'étage de sommation numérique, ce soustracteur étant suivi d'une part d'un troisième registre de stockage dont la sortie est celle du circuit de calcul de désempilement et d'autre part, en parallèle, d'un troisième multiplieur puis d'un quatrième registre de stockage, la sortie de ce quatrième registre étant reliée à la deuxième entrée du soustracteur et l'autre entrée du troisième multiplieur étant reliée à la sortie de la mémoire de stockage du coefficient y. 13. Caméra à scintillation selon l'une des revendica- tions 7 et 8, caractérisée en ce que l'étage de traitement des événements comprend trois circuits de calcul de désempilement, deux diviseurs, un circuit de réalignement temporel, et un circuit de calcul additionnel, les deux diviseurs recevant les signaux Xmj, Ymj, Zmj pour délivrer deux signaux xmj = Xmj/Zmj et Ymj = Ymj/Zmj, ou respectivement les si- gnaux Xmj, Ymj, Emj pour délivrer deux signaux xmj = Xmj/Emj et Ymj = Ymj/Emj, les deux premiers circuits de calcul de désempilement recevant lesdits signaux xmj, Ymj et délivrant les signaux x,y et le troisième recevant le signal Zm ou Em et délivrant le signal E, et le circuit de calcul additionnel recevant d'une part ledit signal E et d'autre part lesdits coefficients de correction pour fournir aux premier et deuxième circuits de calcul de désempilement un coefficient de correction additionnel r. 14. Caméra à scintillation selon l'une des revendica- tions 1 à 13, caractérisée en-ce que les diviseurs sont rem- placés par un seul circuit de division associé à un multi- plexeur-démultiplexeur temporel. 15. Caméra à scintillation selon l'une des revendica- tions 1 à 14, caractérisée en ce que les circuits de calcul de 26057q5 37 - désempilement sont remplacés par un seul de ces circuits au- quel est associé un multiplexeur-démultiplexeur temporel. 16. Caméra à scintillation selon l'une des revendica- tions 1 à 15, caractérisée en ce qu'est prévu un circuit de réjection en amplitude pour la réduction du nombre d'événe- ments à traiter. 17. Caméra à scintillation selon l'une des revendica- tions 1 à 16, caractérisée en ce que des mémoires dites FIFO sont prévues en sortie de chacune des p voies d'acquisition ou en amont de chacun des dispositifs de sommation pondérée numé- rique. 18. Caméra à scintillation selon l'une des revendica- tions 1 à 17, caractérisé en ce que le bus est associé au calculateur sans lui être inclus.