Brevet FR2514904A1 — publié le 19830422

La présente invention concerne une caméra à scin tillation capable de corriger les erreurs dans les positions du rayonnement détecté de façon à donner un tomogramme d'une grande précision. Les erreurs qui se manifestent sur les positions du rayonnement détecté sont dues à la caractérlsti- que de détection non linéaire du détecteur de rayonnement. On utilise habituellement les caméras à scintillation pour détecter le rayonnement ou les rayons gamma émis par un corps humain ou autre, dans lequel on a injecté un radio-isotope, dans le but d'obtenir un tomogramme nécessaire pour faire un diagnostic concernant une partie désirée du corps humain. Le dispositif de détection de rayonnement qui est incorporé à la caméra à scintillation présente une caractéristique de détection affectée d'une non linéarité inhérente, et le tomogramme détecté manifeste une distorsion si on ne corrige pas la non linéarité de la caractéristique de détection. On trouvera ci-après la description d'une caméra à scintillation connue capable de corriger effectivement la non linéarité de la caractéristique de détection. Cette caméra à scintillation connue comprend un détecteur de rayonnement destiné à détecter le rayonnement ou les rayons gamma un circuit de calcul de position destiné à calculer1 à partir des signaux de sortie provenant du détecteur de rayonnement, des positions d'émission de rayonnement qu'on appellera ci-après positions de rayonnement ou positions ; un convertisseur analogique-numérique destiné à convertir les signaux de sortie ou de position en signaux numériques, ces signaux de position provenant du circuit de calcul de position ; une mémoire destinée à enregistrer une valeur de correction de référence pour chacun des points d'intersection d'une matrice donnée qui est définie de la manière désirée sur un plan X-Y, ces valeurs de correction de référence étant utilisées pour corriger les signaux de position, qui ont été convertis en signaux numériques par le convertisseur analogique-numérique, afin de donner des signaux de position exacts représentant les positions exactes, alors que les signaux de position con tiennent déjà des composantes d'erreur à cause de la non linéarité de la caractéristique de détection qui est inhérente au détecteur de rayonnement, et ils ne représentent pas des positions correctes ou exactes sur le plan X-Y correspondant aux positions réelles ; un circuit de calcul de valeur de correction destiné à lire dans la mémoire des valeurs de correction de référence pour quatre points d'intersection adjacents aux positions calculées, afin de corriger les positions de rayonnement calculées, qui ont été calculées par le circuit de calcul de position et numérisées par le convertisseur analogique-numérique, pour donner des positions exactes dans un mode de mesure, et destiné à calculer des valeurs de correction pour corriger les positions calculées afin de donner des valeurs exactes, ces quatre points d'intersection contenant des positions calculées dans une zone formée en reliant mutuellement quatre points d'intersection par des lignes droites ; un additionneur destiné à additionner la valeur de correction calculée à la position calculée ; et un dispositif d'affichage destiné à présenter les positions résultantes pour visualiser un tomogramme concernant une partie désirée du corps humain. Du fait que la caméra à scintillation de ce type effectue la correction de non linéarité du circuit de calcul de position en faisant appel aux techniques de l'informatique, il est relativement facile de procéder aux réglages concernant par exemple la variation horaire, et on peut également parvenir à une correction très précise sans dégrader le pouvoir de résolution, la qualité de l'image, etc. Le circuit de calcul de valeur de correction fonctionne de la manière suivante : on suppose qu'une position de rayonnement calculée P(Xp,Yp) convertie en une valeur numérique par le convertisseur analogique-numérique est contenue, comme le montre la figure 1, dans une zone entourée par quatre points d'intersection A(Xi,Y), B(Xi'Yj+1) C(Xi+1'Yj+1) et D(Xi+1 Yj) d'une matrice donnée dans la mémoire. Dans ce cas, le circuit de calcul de valeur de correction lit dans la mémoire les valeurs de correction de référence (Xi,j; Yi,j), (Xi,j+1, Yi,j+1), (Xi+l,j+1; Yi+1,j+1) et (Xi+1) Yj+1,j) aux quatre points d'intersection A, B, C et D. En ce qui concerne la valeur de correction de référence (Xi,j; Yi,j) qu'on envisage à titre d'exemple en se référant à la figure 2, X. j représente une composante X du vecteur de correction et Yi j représente une composante Y de ce vecteur. On peut dire la meme chose pour les autres valeurs de correction de référence. On suppose maintenant que les positions exactes corrigées des quatre points d'intersection A, B, C et D sont A', B', C', D' et que les points d'intersection obtenus en traçant des lignes verticales à partir de la position P(Xp,Yp) vers les côtés respectifs AB, BC, CD et DA, sont les points Q, R, S et T.On établit dans ce cas les relations suivantes exprimées par les équations (1), (2), (3) et (4) BQ : QA = B'Q' : Q'A' ..... (1) BR : RC = B'R' : R'C' .. . (2) CS : SD = C'S' : S'D' ..... (3) DT : TA , D'T' : T'A' ..... (4) Les relations ci-dessus permettent de déterminer un point d'intersection P' entre les lignes R'T' et Q'S'. Le vecteur PP' devient maintenant un vecteur de correction correspondant à une valeur de correction liée à la position P. Lorsqu'on effectue la correction de non linéarité de cette manière dans la caméra à scintillation, la précision de la correction augmente lorsqu'on augmente le nombre de bits de quantification du convertisseur analogiquenumérique. Cependant, lorsque le nombre de bits de quantification du convertisseur analogique-numérique devient de plus en plus grand, le convertisseur analogique-numérique devient complexe et il présente un coût élevé et une faible vitesse de conversion. Il est donc souhaitable d'employer un convertisseur analogique-numérique ayant un nombre de bits de quantification aussi faible que possible, sans pour autant que ce nombre réduit de bits de quantification ait une influence sur la précision de la correction. Lorsqu'on considère de façon théorique le pouvoir de résolution, on voit qu'il n'est pas nécessaire que le convertisseur analogique-numérique ai un nombre de bits de quantification élevé, et on peut obtenir le nombre de bits -de quantification en appliquant le théorème de l'échar.tillon- nage à la fréquence spatiale qui correspond au pouvoir de résolution déterminé par le dispositif lui-même. Cependant, lorsque le nombre de bits de quar.tifica- tion du convertisseur analogique-numérique devient faible, le pouvoir de résolution n'est théoriquement pas diminué, mais le problème suivant se manifeste en pratique. On suppose que les mailles Mi, M2, M3, ... repré sentées sur la figure 3 sont les mailles unités pour la quantification des signaux de position par le convertisseur analogique-numérique. Tous les points contenus dans cha cune des mailles unités M1, lI2, M3, ... sont rassemblés et considérés comme des points caractéristiques MP1, MP2, MP3, ... par la conversion analogique-numérique du conver tisseur analogique-numérique. Le circuit de calcul de valeur de correction calcule des vecteurs de correction CV1, CV2, CV3, ... qui correspondent aux points caractéris tiques MP1, MP2, MP3, ...Du fait que ces vecteurs de correction CV1, CV2, CV3, ... sont obtenus sur la base des points MP1, MP2, MP3, ... qui sont les points caractéristi ques des mailles unités M1, M2, M3, ..., ils sont appliqués à tous les points contenus dans les mailles unités Mi, M2, M3, ... Par conséquent, toutes les positions corrigées sur la base du vecteur de correction CV1 sont contenues dans une région unité M1' qui est déterminée essentiellement par le vecteur de correction CV1. Toutes les positions corrigées sur la base du vecteur de correction CV2 sont comprises dans une région unité M21 qui est déterminée essentiellement par le vecteur de correction CV2. De façon similaire, toutes les positions corrigées sur la base du vecteur de correction CV3 sont contenues dans une région M3' qui est déterminée essen tiellement par le vecteur de correction C5l3. Comme décrit ci-dessus, les régions unités M1', M2', M3', ... dans lesquelles se trouvent les positions corrigées sont limitées aux régions qui correspondent aux régions unités Mi, M2, M3, ... et les positions des régions unités M1', M2', M3', ... elles-memes sont limitées par des vecteurs de cor rection correspondant aux points caractéristiques MP1, MP2, - MP3, ... des mailles unités Par conséquent, dans les régions dans lesquelles des vecteurs de correction mutuelle ment adjacents doivent prendre un écartement mutuel plus grand que leur écartement avant la correction, il apparat une zone qui n'est pas couverte par les régions unités, comme le montrent les régions unités M1 et M2'. Il apparatt ainsi une région dans laquelle il n'y a pas de positicn corrigée. Dans les régions dans lesquelles des vecteurs de correction mutuellement adjacents doivent se rapprocher par rapport à leur situation avant la correction, il apparat une zone dans laquelle les régions unités se chevauchent mutuellement, comme le montrent les régions unités M1' et M3'. I1 apparat ainsi une région dans laquelle les posi tions corrigées se chevauchent mutuellement. Une correction irrégulière de ce type diminue la précision du tomogramme qui est obtenu. Cependant, il n'est pas souhaitable d'employer un convertisseur analogiquenumérique ayant un nombre de bits de quantification plus élevé que le nécessite le pouvoir de résolution, dans le seul but d'éliminer les erreurs qui résultent d'une correc tion irrégulière. L'invention a pour but de faire disparaître les inconvénients mentionnés ci-dessus. Le but de l'invention est donc de réaliser une caméra à scintillation ayant une structure simple et capa ble de supprimer l'irrégularité dans la correction de la non non linéarité, en utilisant un convertisseur analogique- numérique ayant un nombre de bits de quantification aussi faible que possible, pour réduire le coût, et de donner un tomogramme d'une grande précision. Conformément à l'invention, une caméra à scintil lation comprend un détecteur de rayonnement destiné à détecter le rayonnement qui correspond à chaque événement de scintillation ; un circuit de calcul de position de rayonne ment destiné à calculer des positions de rayonnement à par tir des signaux de sortie qu'il reçoit du détecteur de rayonnement ; un convertisseur analogique-numérique destiné à convertir les signaux de sortie de position analogiques qu'il reçoit du circuit de calcul de position de rayonnement en signaux de position numériques ; une première mémoire destinée à enregistrer, en mode d'étalonnage, l'image d'une source de rayonnement de référence comportant un ensemble de trous disposés en une matrice ; une unité arithmétique destinée à rechercher pour chaque point d'intersection d'une matrice la discordance entre une position réelle dans la source de rayonnement de référence et la position correspondante dans l'image qui est enregistrée dans la première mémoire, et à calculer une valeur de correction de référence à partir d'une valeur de discordance pour chaque point d'intersection ; une seconde mémoire destinée à enregistrer la valeur de correction calculée à son point d'intersection correspondant de sa matrice correspondante ; des premier et second générateurs de données aléatoires destinés à appliquer des premières et des secondes données aléatoires indépendantes à chaque événement de scintillation, en mode de mesure ; un circuit de calcul de valeur de correction destiné à recevoir un signal de position prover.ant du convertisseur analogique-numérique et des premières et secondes données aléatoires provenant des premier et second générateurs de données aléatoires, de manière à prendre les données correspondant A-la composante X dans le signal de position en tant que bits de rangs supérieurs, et les premières données aléatoires en tant que bits de rangs inférieurs pour une position X détectée, et également de manière à prendre les données correspondant à la composante Y dans le signal de position en tant que bits de rangs supérieurs, et les secondes données aléatoires en tant que bits de rangs inférieurs pour une position Y détectée dans le mode de mesure, ce circuit de calcul de valeur de référence lisant également, sur la base de ces données, des valeurs de correction de référence dans des points d'intersection qui entourent la position du signal de position dans la matrice de la seconde mémoire, afin de calculer une valeur de correction concernant la position ; et un additionneur destiné à additionner les valeurs de correction calculées aux signaux de position qui proviennent du convertisseur analogiquenumérique. D'autres caractéristiques et avantages de l'invention seront mieux compris à la lecture de la description qui va suivre de modes de réalisation et: en se référant aux dessins annexés sur lesquels La figure 1 montre un mode de correction destiné à la recherche de positions de rayonnement exactes La figure 2 montre un vecteur comportant des composantes dans les directions X et Y dont origine se trcuve à la position A sur la figure i La figure 3 est une représentation utilisée pour expliquer une difficulté qui apparat lorsque des signaux de position détectée sont quantifiés par un convertisseur analogique-numérique La figure 4 est un schéma synoptique montrant la configuration d'un exemple de caméra à scintillation correspondant à l'invention La figure 5 est une coupe montrant une source de rayonnement de référence employée dans la caméra à scintillation qui est représentée sur la figure 4 La figure 6 est une vue en plan montrant la source de rayonnement de référence La figure 7 montre la configuration détaillée d'un circuit de calcul de valeur de correction qui est employé dans la caméra à scintillation représentée sur la figure 4 ; et La figure 8 montre la configuration de caméra à scintillation à laquelle est en outre ajouté un dispositif de sortie numérique. La figure 4 montre la configuration d'un exemple de caméra à scintillation conforme à l'invention. Au cours de son fonctionnement, la caméra à scintillation met en oeuvre un mode d'étalonnage puis un mode de mesure. On va expliquer sommairement le mode d'étalonnage pour lequel on utilise une source de rayonnement de référence qui est représentée en détail sur les figures 5 et 6. La source de rayonnement de référence comprend un élément de source de rayonnement de référence 101, du type plaque, une plaque de plomb 102 espacée par rapport à l'élément de source 101 et une plaque 103 fixée à la manière d'une couche sur la face inférieure de la plaque de plomb 102 afin d'éviter que la plaque de plomb 102 se courDe.La plaque 103 consiste en une matière telle qu'une résine acrylique qui permet le passage du rayonnement ou des rayons gamma que rayonne l'élément de source lOi. La plaque de plomb 12 comporte un ensemble de trous (par exemple 6x6 = 36 trcus) disposés en matrice. Le mode d'étalonnage est mis en oeuvre en utilisant la source de rayonnement de référence ayant une telle configuration et, dans ce cas, on obtient une valeur de correction de référence à chacun des points d'intersection d'une matrice, à partir d'une valeur de discordance entre une position réelle (ou point d'intersection) de chacun des trous dans la matrice de la plaque de plomb 1G2, et la position mesurée correspondante qui est détectée par un détecteur de rayonnement et calculée par un circuit de calcul de position. La valeur de correction de référence à chaque point d'intersection dans la matrice donnée est enregistrée dans une mémoire pour être utilisée pour le calcul d'une valeur de référence relative à un signal de position détecté dans le cycle de mesure qui sera accompli ultérieurement. Pour mettre en oeuvre le mode de mesure, on injecte un radio-isotope dans un corps humain, par exemple et on dispose le scintillateur au-dessus de l'abdomen du corps humain pour détecter le rayonnement ou les rayons gamma qui sont émis par le radic-isotope injecté. On calcule une position de rayonnement exacte en utilisant la valeur de correction de référence qui a été enregistrée dans la mémoire dans le mode d'étalonnage. On va maintenant décrire en détail en se référant à la figure 4 une caméra à scintillation correspondant à l'invention. La référence 1 désigne un scintillateur destiné à convertir le rayonnement ou les rayons gamma qui sont émis par une source de rayonnement en un ensemble de quanta de lumière. La référence 2 désigne un guide de lumière destiné à guider les quanta de lumière vers un tube photomultiplicateur 3. Le scintillateur 1, le guide de lumière 2 et le tube photomultiplicateur 3 forment un détecteur de rayonnement 4. Le détecteur de rayonnement 4 est connecté à un circuit de calcul de position de rayonnement 5 qui a pour fonction de calculer une position de rayonnement P(Xp,Yp) dans des coordonnées X-Y, sur la base du signal de sortie qu'il reçoit du détecteur de rayonnement 4, et de présenter en sortie des signaux de position analogiques Xp et Yp représentant la position P(Xp,Yp). Comme on l'a déjà indiqué, les signaux de position analogiques Xp et Yp contiennent des composantes d'erreur du fait que le détecteur de rayonnement 4 a une caractéristique de détection non linéaire.Le signal de sortie du détecteur de rayonnement 4 est appliqué à l'entrée d'un analyseur de hauteur d'impulsion 6 qui a pour fonction d'analyser l'amplitude du signal de sortie qui provient du détecteur de rayonnement 4 et d'émettre un signal de sortie de suppression d'effacement U lorsque l'amplitude est dans une plage prédéterminée. Le circuit de calcul de position 5 est connecté à des convertisseurs analogique-numérique 11 et 12. Le convertisseur analogique-numérique 11 convertit le signal de position analogique Xp provenant du circuit de calcul 5 en un signal numérique Xp, tandis que le convertisseur analogique-numérique 12 convertit le signal de position analogique Yp en un signal numérique Yp.Dans toute la description, on désigne les signaux numériques tels que Xp et Yp en ajoutant un tréma (') au-dessus de symboles tels que Xp et Yp représentant les signaux de position analogiques. On utilise une première mémoire 13 et une unité arithmétique 15 dans le mode d'étalonnage. Ainsi, dans le mode d'étalon- nage, une source de rayonnement de référence, représentée sur les figures 5 et 6, est placée face au scintillateur 1 et la mémoire 13 reçoit, en tant que signaux d'adresse et par l'intermédiaire de commutateuu 14a et 14b,des signaux de position numériques Xp et Yp qui proviennent des convertisseurs analogique-numérique li et 12, dans ce cas, ainsi qu'un signal de suppression d'effacement qui provient de l'analyseur de hauteur d'impulsion 6 par l'intermédiaire d'un commutateur 14c.De ce fait, la-configuration matricielle de trous disposés en matrice dans la plaque de plomb 102 de la source de rayonnement de référence, soit en d'autres termes l'image d'une source de rayonnement de référence, est enregistrée dans la mEmoire 13, avec des erreurs dues à la caractéristique de détection non linéaire du détecteur de rayonnement 4. L'unité arithmétique 15 lit dans la mémoire 13 les données d'image correspondant aux points d'intersection successifs dans la matrice, et elle calcule des valeurs de correction de référence aux points d'intersection de la matrice donnée, sur la base de ces données.Une seconde mémoire 16 reçoit le signal de suppression d'effacement provenant de l'analyseur de hauteur d'impulsion 6 et les signaux de position numériques Xp et Yp provenant des convertisseurs analogique-numérique 11 et 12, et elle enregistre des valeurs de correction de référence calculées par l'unité arithmétique 15 et concernant les points d'intersection dans la matrice donnée. Des premier et second générateurs de données aléatoires 17 et 18 reçoivent le signal de suppression d'effacement U qui provient de I'analyseur de hauteur d'impulsion 6, de. façon à générer indépendamment des données aléatoires ayant un nombre de bits prédéterminé. On suppose dans ce mode de réalisation que les générateurs de données aléatoires 17 et 18 sont des registres à décalage à cinq bits. Ce sont naturellement des données aléatoires à cinq bits qui sont présentées en sortie dans ce cas. De plus, les données aléatoires sont renouvelées chaque fois que le signal de suppression d'effacement U est reçu. Un circuit de calcul de valeur de correction 19 reçoit les I, signaux de position numériques Xp et Yp et il lit à des adresses (points d'intersection) de la mémoire 16 désignées par les signaux de position, des valeurs de correction de référence enregistrées à ces adresses. Le circuit de calcul de valeur de correction 19 reçoit simultanément des données de sortie aléatoires Xr provenant du générateur de données aléatoires 17, en tant que bits de rangs inférieurs du signal .. de position numérique Xp provenant du convertisseur analogi que-numérique 11, et des données de sortie aléatoires Yr du générateur de données aléatoires 18, en tant que bits de rangs inférieurs du signal de position nurrérique Yp qui pro- vient du convertisseur analogique-numérique 12.Les données comprenant le signal de position Xp du convertisseur analogique-numérique 11 en tant que bits de rangs supérieurs et les données aléatoires de sortie Xr du générateur de données aléatoires 17 en tant que bits de rangs inférieurs sont traitées en tant que signal de position Xp détecté, et des 1? données comprenant le signal de position Yp du convertisseur analogique-numérique 12 en tant que bits de rangs supérieurs et les données aléatoires de sortie Yr du générateur de données aléatoires 18 en tant que bits de rangs inférieurs sont traitées en tant que signal de position Yp détecté, dans le circuit de calcul de valeur de correction 19, ce qui permet de calculer les valeurs de correction Xp et Yp. Un convertisseur numérique-analogique 20 convertit en un signal analogique Xp la valeur de correction calculée Xp qui provient du circuit de calcul 19, et un convertisseur numérique-analogique 21 convertit en un signal analogique Yp la valeur de correction calculée Yp qui provient du circuit de calcul 19. Un additionneur 22 additionne le signal de position analogique Xp qui provient du circuit de calcul de position 5 et la valeur de correction Xp qui provient du convertisseur numérique-analogique 20. Un additionneur 23 additionne le signal de position analogique Yp qui provient du circuit de calcul de position 5 et la valeur de correction de sortie Yp qui provient du convertisseur numériqueanalogique 21. Un oscilloscope 24 reçoit des signaux de sortie qui proviennent des additionneurs 22 et 23, ainsi qu'un signal de suppression d'effacement U qui provient de i'ana- lyseur de hauteur d'impulsion 6 par l'intermédiaire d'un cir cuit de retard 25, afin de présenter ainsi sur son écran un tomogramme dans lequel la non linéarité est corrigée et la distorsion est supprimée- dans la mesure du possible.Le circuit de retard 25 a pour fonction de retarder le signal de suppression d'effacement depuis l'instant auquel les signaux de position Yp et Yp sont émis par le circuit de calcul de position 5, jusqu'à l'instant auquel'les valeurs de correction qui correspondent aux signaux de position Xp et Yp sont calculées dans le circuit de calcul 19 et sont ajoutées aux signaux de position Xp et Yp qui proviennent du circuit de calcul 5. Un circuit de commande destiné à commar.der en synchronisme l'ensemble du dispositif n'est pas représenté sur le schéma synoptique de la figure 4, pour la commodité de l'explication. On va maintenant décrire en détail le mode d'étalonnage. La source de rayonnement de référence qui est représentée sur les figures 5 et 6, par exemple, est placée dans ce cas face au scintillateur 1. Le rayonnement ou les rayons gamma qui sont émis par l'élément de source de rayonnement 101 traversent la plaque de support 103 pour être rayonnés vers le scintillateur 1 à travers les trous disposés de façon matricielle qui sont formés dans la plaque de plomb 102. Lorsque des rayons gamma sont dirigés vers le scintillateur 1 à travers chacun des trous, une scintillation apparalt à chaque point du scintillateur 1 vers lequel un rayon gamma est dirigé, et des quanta de lumière ayant un niveau qui correspond à l'énergie du rayon gamma sont générés. Les quanta de lumière générés sont transmis au tube photomultiplicateur 3 par le guide de lumière 2. Les quanta de lumière sont amplifiés dans le tube photomultiplicateur 3 de façon à avoir une amplitude appropriée.Les quanta de lumière amplifiés dans le tube photomultiplicateur 3, constituant le signal de sortie du détecteur de rayonnement 4, sont appliqués à l'entrée du circuit de calcul de position 5 et de l'analyseur de hauteur sinaidesortneu d'impulsion 6. Comme indiqué précédemmen, e ecteur de rayonnement 4 comprend des composantes d'erreur, à cause de la non linéarité de la caractéristique de détection Q détecteur de rayonnement 4. Par conséquent, des signaux com prenant des composantes d'erreur sont appliqués à l'entrée -du circuit de calcul 5 et de l'analyseur de hauteur d'impul sion 6. Le circuit de calcul 5 calcule une position d'irradiation ou de rayonnement des rayons gamma, sur la base de ce signal d'entrée. Les signaux de position analogi ques Xp et Yp, représentant la position P(Xp,Yp) qui est obtenue en tant que résultant du calcul, sont appliqués à l'entrée des convertisseurs analogique-numérique 11 et 12 pour être convertis en signaux de position numériques Xp et Yp. Il n'est pas nécessaire que les convertisseurs analogi que-numérique 11 et 12 aient un nombre de bits de quantifi- cation extrêmement élevé, et leur nombre de bits de quanti fication peut être égal au nombre le plus faible ne dégradant pas le pouvoir de résolution.Pour la commodité de l'explica- tion1 on suppose dans ce mode de réalisation que chaque con vertisseur analogique-numérique 11 et 12 a un nombre de bits de quantification égal à cinq. L'analyseur de hauteur d'impulsion 6 détermine si l'énergie des quanta de lumière d'entrée, c'est-à-dire l'énergie du rayonnement, est dans une plage prédéterminée, et il émet un signal de suppression d'effacement U lorsqu'il détermine que l'énergie est dans la plage prédéter minée. En mode d'étalonnage, les commutateurs 14a, 14b et 14c sont placés sur la position correspondant au contact I, tandis que les interrupteurs 14d, 14e et 14f demeurent ouverts. Par conséquent, les signaux de position numériques 'I Xp et Yp sont appliqués à la mémoire 13 en tant que signaux d'adresse et ils désignent les adresses correspondantes dans la mémoire 13. Dans le cas où le signal de suppression d'ef facernent U a été appliqué à la mémoire 13 à partir de l'ana lyseur de hauteur d'impulsion 6, un "1" logique est ajouté au contenu des adresses que désignent les signaux de posi tion numériques Xp et Yp. L'image de la configuration matricielle de trous qui est formée dans la plaque de plomb 102 de la source de rayonnement de référence, c'est-à-dire l'image d'une source de rayonnement de référence, est enre gistrée de cette manière dans la mémoire 13, avec des compo santes d'erreur dues à la caractéristIque de détection non linéaire du détecteur de rayonnement 4. Les données d'image qui sont enregistrées dans la mémoire 13 sont lues par l'unité arithmétique 15 qui est constituée par exemple par un micro-ordinateur.L'unité arithmétique 15 mémorise dans sa mémoire interne les positions réelles (ou positions exactes) dans la matrice de trous qui est formée dans la plaque de plomb 102. L'unité arithmétique 15 lit les données d'image dans la mémoire 13 et elle calcule des valeurs de correction aux points d'intersection dans la matrice donnée. Le calcul de valeurs de correction est effectué de manière à rechercher la discordance entre une position exacte (ou un point d'intersection exact) dans la matrice donnée enregistrée précédemment dans la mémoire interne de l'unité arithmétique 15, et la position correspondante (ou le point d'intersection correspondant) dans la matrice qui figure dans les données d'image enregistrées dans la mémoire 13, et de manière à calculer une valeur de correction sur la base de cette valeur de discordance. La valeur de correction de référence calculée à chaque point d'intersection dans la matrice donnée est enregistrée dans la seconde mémoire 16. L'opération d'étalonnage en mode d'étalonnage, décrite ci-dessus, est ainsi terminée. On va maintenant décrire le mode de mesure. En mode de mesure, on détecte la position du rayonnement ou des rayons gamma qui sont émis par un radioisotope injecté par exemple dans un corps humain. Dans ce cas, le scintillateur 1 est placé au-dessus du corps humain et les commutateurs 14a, 14b et 14c sont placés sur les contacts II, tandis que les interrupteurs 14d, 14e et 14f sont fermés. Le détecteur de rayonnement 4 détecte le rayonnement qui correspond à chaque scintillation, et le circuit de calcul 5 calcule la position de la scintillation. Le processus est similaire à celui du mode d'étalonnage, et on ne le décrira donc pas.Du fait que les commutateurs 14a et 14b sont placés sur les positions correspondant aux contacts II, les signaux de position numériques à cinq bits Xp et Yp qui proviennent des convertisseurs 11 et 12 sont appliqués à la seconde mémoire 16 et au circuit de calcul de valeur de correction 19. Du fait que les interrupteurs 14d et 14e sont fermés, les signaux de position analogiques Xp et Yp provenant du circuit de calcul de position 5 sont respectivement appliqués aux additionneurs 22 et 23. Du fait que le commutateur 14c est placé sur la position correspondant au contact II et que l'interrupteur 14f est fermé, le signal de suppression d'effacement U qui provient de l'analyseur de hauteur d'impulsion 6 est appliqué à l'entrée des générateurs de données aléatoires 17, 18 et du circuit de retard 25.Dans le cas où l'analyseur de hauteur d'impulsion 6 applique le signal de suppression d'effacement U à la mémoire 16, des valeurs de correction de référence à quatre points d'intersection sont lues à ces quatre points d'intersection de la matrice, dans la mémoire 16, au moyen du circuit de calcul de valeur de correction 19. Ces quatre points d'intersection dans la matrice présente dans la mémoire 16 entourent la position P(Xp,Yp) qui est représentée par les ,. .. signaux de position Xp et Yp. Simultanément, des données aléatoires à cinq bits Xr qui sont ajoutées, en tant que t' bits de rangs inférieurs, au signal de position Xp et des données aléatoires à cinq bits Yr qui sont ajoutées en tant que bits de rangs inférieurs au signal de position Yf sont respectivement appliquées au circuit de calcul de valeur de correction 19, à partir des premier et second générateurs 17 et 18.Le circuit de calcul de valeur de correction 19 reçoit donc un signal de position à dix bits (Xp+Xr) qui comporte un signal de position à cinq bits Xp en tant que bits de rangs supérieurs et des données aléatoires à cinq bits Xr en tant que bits de rangs inférieurs, et un signal de position à dix bits (Yp+Yr) qui comprend un signal de ll position à cinq bits Yp en tant que bits de rangs supérieurs et des données aléatoires à cinq bits Yr en tant que bits de rangs inférieurs.En mettant en oeuvre les équations (5) et (6), qui sont déterminées à partir des équations (1), (2), (3) et (4) exprimées en relation avec les figures 1 et 2 et qu'on décrira ultérieurement, le circuit de calcul de valeur de correction 19 calcule la valeur de correction de position P(Xp,Yp) en utilisant les valeurs de correction de référence aux quatre points d'intersection qui entourent la position P(Xp,Yp). La valeur de correction calculée Xp est appliquée à l'entrée du convertisseur numérique-analogique 20, convertie en composantes analogiques Xp, puis appliquée à l'entrée de l'additionneur 22. La composante Yp, relative à Y, est appliquée à l'entrée du convertisseur numérique-analogique 21, convertie en un signal analogique Yp, puis applique à l'entrée de l'additionneur 23. Le signal analogique Xp est additionné dans l'additionneur 22 au signal analogique Xp qui provient du circuit de calcul de position 5. Le signal analogique Yp est additionné dans l'additionneur 23 au signal analogique Yp qui-provient du circuit de calcul de position 5. Les résultats des additions accomplies dans les additionneurs 22 et 23 sont appliqués à l'oscilloscope 24 en tant que signaux de position corrigés. Lorsque le signal de suppression d'effacement U est appliqué, par le circuit de retard 25 dans ce cas, l'oscilloscope 24 présente de façon visible sur l'écran 26 une position qui est déterminée dans le système de cordonnées X-Y par les signaux de position corrigés. Comme décrit ci-dessus, l'écran présente un tomogramme d'une grande précision. On va maintenant décrire le circuit de calcul de valeur de correction 19 en se référant au schéma synoptique détaillé qui est représenté sur la figure 7. Le circuit de calcul 19 reçoit le signal de position à cinq bits Xp du convertisseur analogique-numérique 11 en tant que bits de rangs supérieurs (bits de fort poids) de la composante X, soit Xp, au signal de pos~vlcn, et le signal de position à cinq bits Yp du convertisseur analogique-numériquè 12 en tant que bits de rangs supérieurs (bits de fort poids) de la composante Y, soit Yp, du signal deposition. Le circuit de calcul 19 reçoit également les données aléatoires à cinq bits Xr du générateur de données aléatoires 17, en tant que bits de rangs inférieurs de la composante X du signal de position et les données aléatoires à cinq bits Yr du générateur de données aléatoires 18, en tant que bits de rangs inférieurs de la composante Y. De façon plus concrète, les signaux de position Xp sont appliqués à l'entrée d'un circuit de désignation d'adresse 204 dans lequel ces signaux forment un signal d'adresse à dix bits qui comprend cinq bits du signal de position Yp en tant que bits de rangs supérieurs et cinq bits du signal de position Xp en tant que bits de rangs inférieurs. Ce signal d'adresse est appliqué à la mémoire 16 pour désigner le point A représenté sur la figure 1.En synchronisme avec un signal sous forme d'impulsion qui est appliqué de l'extérieur, le circuit de désignation d'adresse présente successivement en sortie un signal d'adresse formé en additionnant "1" au signal d'adresse à dix bits et désignant le point "D", un signal d'adresse formé en ajoutant en outre "32" en code décimal lorsque la matrice de la memoire 16 est une matrice de dimensions 32x32, et désignant le point "C", et un signal d'adresse formé en soustrayant "1" du signal d'adresse qui désigne le point "C", le signal obtenu désignant le point "B". Le circuit de désignation d'adresse 204 applique également un signal de commande C aux registres 2061, 2062, 2063, 2064, 2065, 2066, 2067 et 2068.Les registres 2061 - 2068 mémorisent sous la commande du signal de commande C les valeurs de correction de référence aux points A, B, C et D qui sont lues dans la mémoire 16. De façon plus concrète, les registres 2061 et 2065 mémorisent les valeurs de correction de référence #Yi,j et bai 4Xi J au point A. Les registres 2062 et 2066 mémorisent les valeurs de correction de référence #Yi+1,j et #Xi+1,j au point D. Les registres 2063 et 2067 mémorisent les valeurs de correction de référence Y1+1,j+1 et bXi+l,j+l au point C. Les registres 2064 et 2068 mémorisent les valeurs de correction de référence DYi,j+l et Xi,j+1 au point B. Le générateur de données aléatoires 17 applique les données aléatoires à cinq bits Xr à un registre 2081 qui les mémorise. Comme indiqué précédemment, on utilise les données aléatoires Xr en tant que bits de rangs inférieurs -de la composante X dans le signal de position, et elles désignent la position du point P d'une manière aléatoire lorsqu'on prend le point A pour origine. Le générateur de données aléatoires 18 applique les données aléatoires à .. cinq bits Yr à un registre 2082 qui les mémorise. On utilise de façon similaire les données aléatoires Yr en tant que bits de rangs inférieurs de la composante Y dans le signal de position et elles représentent la position du point P de façon aléatoire lorsqu'on prend pour origine le point A. Le circuit de calcul de dXp, 2101, calcule la valeur de correction X1,p du signal de position Xp sur la t, t, base des données aléatoires Xr et Yr qui proviennent des registres 2081 et 2082, et des valeurs de correction de référence qui proviennent des registres 2065 - 2068. Le circuit de calcul de tYp 2102 calcule de façon similaire la valeur de correction tYp du signal de position Yp, sur la base des données aléatoires Xr et Yr qui proviennent des registres 2081 et 2082, et des valeurs de correction de référence qui- proviennent des registres 206i - 2064. On va maintenant expliquer le fonctionnement. Lorsque le circuit de désignation d'adresse 204 reçoit les signaux de position Xp et Yp, il présente à la mémoire 16 un signal dtadresse à dix bits qui comprend cinq bits du signal Xp en tant que bits de rangs supérieurs et cinq bits du signal Yp en tant que bits de rangs inférieurs, et il désigne le point d'intersection prédéterminé A-dans la mémoire 16, au moyen de ce signal d'adresse. Le circuit de désignation d'adresse 204 applique également au même moment le signal C aux registres 2061 - 2068. Sous la commande de ce signal de commande C, les registres 2061 et 2065 enre gistrent des valeurs de correction Y1,j et hXi ; qui sont lues au point d'intersection A que désigne le signal d'adresse.Lorsque les valeurs de correction dYi,j et hXi,J sont mémorisées, le signal d'adresse auquel la quanti té "1" est additionnée pour désigner le point D est émis par le circuit de désignation d'adresse 202 vers la mémoire 16. Les valeurs de correction de référence ti+1 J et X1+11j au point d'intersection D dans la mémoire 16 sont ainsi mémorisées respectivement dans les registres 2062 et 2066. Le signal d'adresse auquel la quantité 32 seulement a été additionnée, en mode décimal, pour désigner le point C est ensuite émis vers la mémoire 16. Les valeurs de correction de référenceYi+1,j+1 et # Xi+1,j+1 sont lues au point d'intersection C dans la mémoire 16 et elles sont mémorisées dans les registres 2063 et 2067. Le signal d'adresse duquel la quantité "1" est soustraite pour désigner le point B est ensuite émis par le circuit de désignation d'adresse 204 vers la mémoire 16.Les valeurs de correction de référence tj+l et #Xi,j+1 sont alors lues au point d'intersection B dans la mémoire 16 et elles sont mémorisées respectivement dans les registres 2064 et 2068. Conformément à la description précédente, les valeurs de correction de référence qui proviennent de la mémoire 16 sont mémorisées dans les registres 2061 - 2068 tandis que les données Yr qui sont utilisées en tant que bits de rangs inférieurs de la composante Y dans le signal de position d'entrée et les données Xr qui sont utilisées en tant que bits de rangs inférieurs de la composante X dans le signal de position d'entrée, sont appliquées par les générateurs de données aléatoires 17 et 18 aux registres respectifs 2081 et 2082 et elles sont mémorisées dans ces registres. Les valeurs de correction mémorisées dans les registres 2061 et 2064, ainsi que les données aléatoires Xr et Yr mémorisées dans les registres 2081 et 2082 sont appliquées au circuit de calcul de AY"p1 2102. De façon similai- re, les valeurs de correction qui sont mémorisées dans les registres 2065 et 2068, ainsi que les données aleatoires Xr et Yr qui sont mémorisées dans les registres 2081 et 2082 sont appliquées au circuit de calcul de JXp, 2101. Les circuits de calcul 2101 et 2102 calculent les valeurs de correction AXp et êYp relatives aux composantes X et Y du signal de position.On peut obtenir les valeurs de correc tion A,X"p et AYp en résolvant les équations (5) et (6) ci-après dans les circuits de calcul 2101 et 2102. #Xp = 1/a(#Xi,j.(a-Xr).(a-Yr) + (BXi,j j.Xr.(a-Yr) + #Xi,j+1.(a-Xr).Yr +#Xi,j+1.Xr .Yr) ... (5) #Y1p = 1/a(#Yi,j.(a-Xr).(a-Yr) + #Yi+1,j.Xr.(a-Yr) + #Yi,j+1.(a-Xr).Yr + #Yi+1,j.Xr.Yr) ... (6) Dans ces équations, "a" représente le nombre de lignes (ou de colonnes) de la matrice dans la mémoire 16. On peut facilement obtenir (a-Xr) et (a-Yr) par calcul binaire des compléments à 32 de Xr et Yr. On peut aisément réaliser les circuits de calcul 2101 et 2102 au moyen de techniques classiques, de façon à résoudre les équations (5) et (6). Comme il ressort de ce qui précède, l'invention offre une structure dans laquelle deux générateurs de données aléatoires indépendants, 17 et 18, fournissent en sortie deux données aléatoires indépendantes pour chaque événement de scintillation et dans laquelle ces deux données aléatoires indépendantes sont ajoutées en tant que bits de rangs inférieurs, aux signaux de position Xp et Yp qui sont soumis à une conversion analogique-numérique par les convertisseurs analogique-numérique 11 et 12.Par conséquent, un point dans une maille unité de quantification M1, par exemple, soumis à une conversion. analogique-numérique de la manière représentée sur la figure 3, n'est pas reconnu comme un point caractéristique MP1 soumis à une conversion analogique-numérique, mais comme un point dans la maille unité M1 qui est déterminé de façon aléatoire conformément à la valeur des données aléatoires. Du fait que les données aléatoires sont ajoutées, en tant que bits de rangs inférieurs, aux signaux de position de sortie qui proviennent des convertisseurs analogique-numérique 11 et 12, ces données déterminent de façon aléatoire dans la maille unité de quantification M1 le point caractéristique qui est défini par les signaux de position qui proviennent des convertisseurs analogique-numérique 11 et 12.Ainsi, les points situés dans la maille unité Ml ne sont pas reconnus comme un point caractéristique mais sont traités comme un ensemble de points répartis dans la maille unité Mi. On détermine une valeur de - correction portant sur le point reconnu de cette manière. Par conséquent, les positions qui résultent de l'addition de valeurs de correction ainsi obtenues aux signaux de posi tion Xp et Yp ne sont pas concentrées dans la région M1', mais sont réparties de façon régulière, comme dans le cas dans lequel le signal de position détecté est soumis à une conversion analogique-numérique suffisamment fine, c'est-à- dire le cas dans lequel le nombre de bits de quantification de chaque convertisseur analogique-numérique il et 12 est suffisamment augmenté.On peut dire la même chose au sujet des autres mailles unités MP2, MP3, On peut obtenir des valeurs de correction dans la région de positions qui est définie par les données aléatoi res ajoutées, qui n'ont pas de relation directe avec les positions de rayonnement réelles, mais l'étalement des signaux de position sous l'effet des données aléatoires est accompli d'une manière telle que des points étalés en prati que soient tout d'abord regroupés en un point caractéristi que par la conversion analogique-numérique, puis étalés à nouveau dans la région dans laquelle ces points étaient étalés à l'origine. Ainsi, lorsqu'on choisit pour le nombre de bits de quantification de la conversion analogique numérique une valeur qui conserve le pouvoir de résolution désiré , on peut parvenir à un calcul de position d'une grande précision. Conformément à l'invention, on peut donc obtenir un tomogramme d'une grande précision dans lequel la non linéarité est corrigée, même si le nombre de bits de quan tification de chaque convertisseur analogique-numérique est juste suffisant pour conserver le pouvoir de résolution. Dans le cas où des signaux de position obtenus après correction sont présentés en sortie sous la forme de signaux numériques destinés à être appliqués au dispositif informatique ou à un dispositif analogue, il peut exister un additionneur numérique 31 (voir la figure 8) qui est des tiné à ajouter le signal de sortie tXp provenant du circuit de calcul de valeur de correction 19 aux données comprenant le signal de position ip du convertisseur analogique-numérique 11 en tant que bits de rangs supérieurs et les données aléatoires Xr du générateur de données aléatoires 17 en tant que bits de rangs inférieurs, et un additionneur numérique 32 destiné à additionner le signal de sortie bYp provenant du circuit de calcul de valeur de correction 19 aux données qui comprennent le signal de position Yp du convertisseur analogique-numérique 12 en tant que bits de rangs supérieurs, et les données aléatoires Yr en tant que bits de rangs inférieurs, ce qui permet de présenter des signaux additionnés en sortie des additionneurs numériques 31 et 32. On peut remplacer les générateurs de données aléatoires par des moyens destinés à générer des données constituant une approximation de données aléatoires. On peut par exemple échantillonner par le signal de suppression d'effacement U les signaux de sortie de deux générateurs d'ondes triangulaires, et soumettre les valeurs ainsi obtenues à une conversion analogique-numérique. On peut employer dans ce cas deux générateurs d'ondes triangulai res produisant des ondes triangulaires qui ne soient pas synchronisées mutuellement, afin d'obtenir des donnes aléatoires indépendantes relatives aux composantes X et Y, soit Xp et Yp. Dans le mode de réalisation considéré, chaque convertisseur analogique-numérique 11 et 12 comporte cinq bits mais leur nombre de bits de quantification peut être inférieur, à condition de ne pas dégrader le pouvoir de résolution. Le nombre de bits de quantification des convertisseurs analogique-numérique il et 12 et celui des générateurs de données aléatoires 17 et 18 est égal ici à cinq, mais ces éléments peuvent avoir des nombres de bits de quantification différents les uns des autres, en fonction du mode d'utilisation. Bien entendu diverses modifications peuvent être apportées par l'homme de l'art aux dispositifs et procédés qui viennent d'être décrits uniquement à titre d'exemples non limitatifs, sans sortir du cadre de l'invention. REVENrDICATIONS 1. Caméra à scintillation, caractérisée en ce qu'elle comprend : un détecteur de rayonnement (4) destiné à détecter le rayonnement correspondant à chaque événement de scintillation ; un circuit de calcul de position de rayonnement (5) destiné à calculer les positions de rayonnement à partir des signaux de sortie qu'il reçoit du détecteur de rayonnement ; un circuit convertisseur analogique-numérique (11, 12) destiné à convertir les signaux de sortie de position analogiques qu'il reçoit du circuit de calcul de position de rayonnement pour donner des signaux de position numériques ; une première mémoire (13) destinée à enregistrer en mode d'étalonnage l'image d'lue source de rayonnement de référence qui comporte un ensemble de trous disposés en une matrice ; une unité arithmétique (15) destinée à rechercher pour chacun des points d'intersection d'une matrice la discordance entre une position réelle dans la source de rayonnement de référence et la position correspondante dans l'image qui est enregistrée dans la première mémoire, et à calculer une valeur de correction de référence à partir de la valeur de discordance pour chacun des points d'intersection ; une seconde mémoire (16) destinée à enregistrer la valeur de correction calculée au point d'intersection correspondant de sa matrice correspondante ; des premier et second générateurs de données aléatoires (17, 18) destinés à appliquer des premieres et secondes données aléatoires indépendantes à chaque événement de scintillation, en mode de mesure ; Lcci circuit de calcul de valeur de correction (19) qui est destiné à recevoir des signaux de position provenant du circuit convertisseur analogique-numérique et les premières et secondes données aléatoires provenant des preier et second générateurs de données aléatoires, afin de prendre en tant que bits de rangs supérieurs des données correspondant à une composante X dans les signaux de position, et en tant que bits de rangs inférieurs les premières données aléatoires, pour une position X détectée, et afin de prendre également en tant que bits de rangs supérieurs des données correspondant à une composante Y dans les signaux de position, et en tant que bits de rangs inférieurs les secondes données aléatoires, pour une position Y détectée, en mode de mesure, et destiné à lire, sur la base de ces données, des valeurs de correction de référence aux points d'intersection qui entourent la position des signaux de position dans la matrice de la seconde mémoire, afin de calculer une valeur ae correction liée à la position; et un circuit additionneur (22, 23) destiné à additionner les valeurs de correction calculées aux signaux de position qui proviennent du circuit convertisseur analogique-numérique. 2. Caméra à scintillation selon la revendication 1 caractérisée en ce que le circuit de calcul de valeur de correction (19) comprend un circuit de désignation d'adresse (204) qui est destiné à appliquer un signal d'adresse pour lire une valeur de correction de référence à chacun des points d'intersection dans la matrice de la seconde mémoire (16), sur la base de signaux de position provenant du circuit convertisseur analogique-numérique ; des registres (206l - 2068) destinés à enregistrer les valeurs de correction de référence lues ; des premier et second circuits de registre (208 - 2082) destinés à enregistrer les premières et secondes données aléatoires ; et un circuit de calcul (210 - 2102) qui est destiné à calculer les valeurs de correction concernant les signaux de position, en mode de mesure, sur la base des valeurs de correction de référence qui sont enregistrées dans les registres (2061 - 2068) et des premières et secondes données aléatoires qui sont enregistrées dans les premier et second circuits de registre (2101 2102)