Brevet FR2606160A1 — publié le 19880506

PROCEDE DE RECONSTRUCTION D'IMAGES DE COUPES PARALLELES SUCCESSIVES D'UN OBJET CONTENANT DES SOURCES EMETTRICES DE RAYONNEMENT GAMMA DESCRIPTION La pressente invention concerne un procédé de reconstruction d'images de coupes parallèles successives d'un objet contenant des sources émettrices de rayonnement gamma, à l'aide d'une gamma-camera å détecteur plan parallèle aux plans de coupe. Elle s'applique à la reconstitution images de coupes parallèles d'un objet, et notamment, à la reconstitution d'images tomographiques d'une ou plusieurs coupes parallèles successives d'un organe du corps humain, contenant des sources émettrices de rayonnement gamma. On sait qu';l est possible de reconstruire les images de coupes parallèles successives d'un objet ou d'un organe contenant des sources émettrices de rayonnement gamma, å l'aide d'une gamma-caméra à détecteur plan fixe, associe å un masque plan, parallèle au détecteur ; ce masque comporte un assemblage de motifs formés chacun par une mosaïque codée d'éléments transparents et opaques aux rayonnements produits par les sources contenues dans ('objet ou l'organe. Le principe de base utilise dans l'obtention d'une image à partir d'une mosaique codée d'éléments transparents et opaques, est simple. Une source ponctuelle de rayonnementS provoque l'apparition d'une image de chaque élément transparent sur le détecteur. La dimension de l'image dépend de la distance qu; sépare la source du détecteur, tandis que l'emplacement de l'image sur Le détecteur dépend du déplacement latéral de la source par rapport au détecteur. Si chaque motif du masque comporte N éléments transparents, on obtient sur le détecteur N images d'une source occupant des positions predéterminees dans un reperde fixe de référence.Les positions des N éléments transparents du masque sont codées ; il en résulte que pour une source ayant une position prédéterminée dans un plan parallèle à celui du détecteur et pour un éloignement prédéterminé de ce plan par rapport au détecteur, N images correspondant respectivement à ces N éléments transparents, se forment sur Le détecteur, en des emplacements codés dépendant de la position de la source dans son plan. La dimension de chaque image dépend de l'éloignement du pLan contenant la source par rapport au détecteur. Les informations codées concernant les positions codées des images, les dimensions de ces images, ainsi que les intensités des rayonnements reçus par le détecteur aux emplacements de ces images sont enregistrées dans une mémoire. Ces informations enregistrées ne permettent pas de reconstruire directement les images de différentes coupes d'un objet, dans des plans parallèles au détecteur. En effet, les images des éléments transparents du masque, obtenues grâce aux rayonnements issus de différentes sources d'un même plan ou de plans différents, se chevauchent. Des images correctes des coupes d'un objet ou d'un organe, ne peuvent être obtenues que grâce à des méthodes de reconstruction qui permettent de déterminer la localisation de chaque source, ainsi que l'intensité du rayonnement émis par celle-ci ; ces méthodes utilisent des traitements des informations enregistrées en mémoire, qui sont de deux types : traitement par déconvolution ou traitement par corrélation. Le principe de reconstruction d'images de coupes d'un objet à l'aide d'une gamma-caméra à détecteur plan, associée à un masque à éléments transparents et opaques, par traitement de convolution ou de corrélation des signaux fournis par cette gamma-catiéra, est décrit dans un article (1) intitulé "Coded aperture imaging with uniformly redendant arrays" par E.E.FENIMORE et T.M. CANNON- revue "Apptied Optics"- Volume 13, nt3, February 1, 1978, pages 337 à 347. Un autre article (2) des mêmes auteurs, intitulé "Tomographical imaging using uniformly redendant arrays" dans la revue "Applied Optics"- Volume 18, n07, April 1, 1979, pages 1052 à 1057 décrit plus particulièrement l'application à la tomographie, du principe de l'article (1) précédent. Un autre article (3) intitulé "An evaluation of techniques for stationary coded a aperture three-dimensional imaging in nuclear medicine" de J.S. FLEMING et B.A.GODDARD, paru dans la revue "Nuclear instruments and methods in physics research" 221-1984- North Holland, pages 242 à 246, décrit une méthode de reconstruction d'images de coupes d'un objet ou d'un organe par traitement de déconvolution des informations issues du détecteur d'une gammacaméra associée à un masque. D'une façon générale, les traitements effectués selon les procédés connus sont longs et compliqués et un problème non résolu, quel que soit le procédé de traitement utilisé, est le bruit de fond important qui apparat dans la reconstruction de chaque image ; l'augmentation du rapport du signal (utile à la reconstruction d'images), au bruit (notamment crée par les sources environnantes de celle qui fournit un signal utile), n'est pas facile à obtenir. Le premier article (1) montre notamment qu'un traitement par déconvolution fournit une mauvaise reconstruction dimages car le bruit reste prédominant. Une meilleure reconstruction d'image est obtenue par un traitement de corrélation qui diminue sensiblement le bruit. Ceci est d'ailleurs confirmé dans les deux autres articles (2), (3). L'amélioration du rapport signal/bruit peut être obtenue, comme le montrent notamment les articles (1) et (2) en disposant les éléments transparents d'un masque selon une répartition aléatoire. Cette répartition qui améliore la valeur du rapport signal/bruit ne permet cependant pas d'obtenir des images extrêmement satisfaisantes, le traitement des informations par corrélation ne parvenant pas à accrottre de façon significative, la valeur du rapport signal/bruit. Dans ces dispositifs connus, Les éléments transparents sont nombreux et occupent cinquante pour cent de la surface du masque. La présente invention a pour but de remédier aux inconvénients des procédés connus de reconstruction d'images de coupes d'un objet, et notamment de fournir un procédé rapide permettant de reduire très sensiblement l'effet du bruit de fond sur les images reconstruites. Ce but est atteint, notamment en utilisant un masque dans lequel les éléments transparents de chaque motif se répartissent d'une façon pseudo-aléatoire prédéterminé, et en effectuant un traitement itératif particulier des valeurs de comptages fournies par la gammacaméra. L'invention concerne un procédé de reconstruction d'images de coupes parallèles successives d'un objet contenant des sources émettrices de rayonnements gamma, à l'aide d'une gamma-caméra à détecteur plan parallèle aux plans de coupes, cette caméra ayant une pluralité de sorties repérées chacune par un code et correspondant respectivement aux différentes coordonnées d'impacts possibles et à l'énergie des rayonnements gamma reçus par le détecteur, cette énergie est mesurée par une valeur numérique fournie sur la sortie correspondant à un impact, le procédé consistant à interposer entre l'objet et le détecteur, parallèlement à ce détecteur, un masque plan formé par assemblage de motifs identiques, de forme rectangulaire ou carrée et comportant chacun une mosaïque d'éléments transparents dans une surface opaque aux rayonnements, ces éléments transparents occupant dans le motif des positions prédéterminées repérées selon une matrice de codage, la reconstitution des images des coupes étant effectuée par un traitement des valeurs de comptages enregistrées dans une mémoire, ces valeurs étant sélectionnées selon des matrices de décodage déduites des matrices de codage et correspondant aux images sur le détecteur d'au moins un motif du masque éclairé par lesdites sources, les dimensions et les emplacements de ces images sur le détecteur dépendant du motif choisi dans le masque, des coordonnées des sources dans chaque plan, par rapport à un repère de référence, ainsi que l'éloignement de chaque plan par rapport au détecteur, le procédé consistant à effectuer les opérations suivantes :: - on sélectionne une portion du détecteur pour chaque éloignement de plan de coupe Pi par rapport au détecteur, - pour chaque plan de coupe Pi, on déplace la matrice de décodage pour repérer les sources éventuellement présentes dans ce plan, ce déplacement correspondant à un balayage de toutes les sources repérées par leurs coordonnées (i, j) dans le plan Pi, et par des valeurs de comptages mesurées, les sources repérées étant celles qui sont situées sur des droites passant par le centre géométrique du détecteur, - à partir des valeurs de comptages correspondantes, on reconstruit les activités des sources de chaque plan, caractérisé en ce que l'on effectue alors les traitements suivants - pour chaque plan Pi d'éloignement prédéterminé, on effectue une comparaison entre L'activité de chaque source de coordonnées (i, j) ainsi reconstruite dans le plan Pi, et les activités reconstruites respectives pour des sources voisines situées dans des plans voisins Pi+1 et Pi-1, ayant les mêmes coordonnées (i, j) que la source considérée dans le plan Pi, - on détermine si les écarts entre l'activité de chaque source de coordonnées (i, j) dans chaque plan Pi et les activités respectives des sources de mêmes coordonnées (i, j) dans les plans Pi-l et Pi+l voisins dépassent un seuilprédéterminé, ce dépassement indiquant le repérage d'une source de coordonnées (i, j) dans le plan Pi, - on attribue une activité à chaque source repérée, - on calcule pour chaque source repérée la valeur de comptages que produirait cette source à partir de l'activité qui lui a été attribuée, - on soustrait, pour chaque source repérée, la valeur de comptages effectivement mesurée précédemment pour la source de mêmes coordonnées (i, j) dans le même plan, de la valeur de comptages calculée pour la source repérée, de manière à obtenir une valeur de comptages corrigée pour la source repérée, - on détermine L'activité de chaque source repérée à partir de la valeur de comptages corrigée correspondante et on recommence les traitements précédents selon un mode itératif jusqu'à ce que la valeur corrigée obtenue pour chaque source soit invariable. Selon une autre caractéristique, chacun desdits éléments transparents d'un motif est repéré par un numéro de ligne et un numéro de colonne correspondant respectivement aux deux nombres d'une paire de nombres d'une série bidimensionnelle équivalente à une série monodimensionnellepseudo-aléatoire. Les caractéristiques et avantages de l'invention ressortiront mieux de la description qui va suivre, donnée en référence aux dessins annexés dans lesquels - la figure 1 est un schéma qui permet de mieux comprendre l'utilisation d'un masque interposé entre La caméra et l'objet, dans le procédé de l'invention, - la figure 2 est un tableau permettant de déterminer la répartition des éléments transparents dans chaque motif du masque utilisé dans l'invention, - la figure 3 représente schématiquement l'un des motifs du masque, - la figure 4 représente schématiquement un système permettant de mettre en oeuvre le procédé de L'invention. La figure 1 permet de mieux comprendre le principe d'un procédé de reconstitution d'images de coupe d'un objet, utilisant une gamma-caméra associée à un masque comprenant des motifs à éléments transparents tels qu'utilisés par exemple dans le procédé de CANNON et FENIMORE mentionné plus haut, ou dans le procédé de l'invention Sur cette figure, on a représenté schématiquement L'écran à pourtour circulaire d'un détecteur 1, d'une gammacaméra. Les rayonnements gamma sont émis par des sources situées dans un objet ou. un organe 2. Le détecteur 1 est un détecteur plan sans collimateur. Les coordonnées des points d'impact des rayonnements sur ce détecteur sont repérées par rapport à un repère de référence X, Y, dont l'origine est par exemple confondue avec le centre géométrique O du détecteur.Les sources appartiennent par exemple à différents plans de coupe P1, P2, P3, de L'organe ou de l'objet 2, parallèles au plan P du détecteur 1. Les éloignements de ces plans par rapport au détecteur sont repérés sur L'axe OZ perpendiculaire aux axes XY du repère de référence. Les coordonnées des sources sont repérées par rapport au repère de référence OXY, et par les éloignements respectifs Z1, Z2, Z3 des plans P1, P2, P3 qui contiennent ces sources. Ces éloignements sont repérés sur l'axe OZ, par exemple par rapport au plan P du détecteur 1. Un masque M, contenu dans un plan P M est interposé entre le détecteur 1 et l'objet ou l'organe 2, dont on veut reconstituer les images des coupes, dans les plans P1, P2, P3 par exemple. Le plan P est parallèle au plan du M détecteur. Ce masque comporte un assemblage de motifs tels que le motif M1, formés chacun par une mosa;que d'éléments transparents 3 et opaques 4 aux rayonnements émis par les sources.Pour simplifier la figure, on na représenté qu'un seul motif de L'assemblage de motifs constituant le masque M. De la même manière, on n'a représenté que quelques éléments transparents et opaques de la mosaique constituant chaque motif. Chaque motif, tel que M1 par exemple, éclairé par une source de rayonnement gamma, a une image sur le détecteur 1. Selon l'assemblage des motifs dans le masque M, les images de ces motifs sont organisées selon une matrice ayant une forme bien définie, dépendant de l'assemblage des motifs du masque, et spécifique à chaque source. L'étendue de chaque image d'un motif dépend de l'éloignement de la source produisant cette image. Les dimensions des images d'un motif, ainsi que l'intensité des rayonnements reçus par le détecteur pour chaque image, dépendent donc des coordonnées de la source correspondante par rapport au repère de référence, et de l'éLoignement du plan de la source par rapport au détecteur. La dimension plane de L'objet est limitée par les distances objet-masque et masque-détecteur.Par contre, avec un masque constitué de nombreux motifs de base, il est possible de choisir dans le plan du détecteur une portion de ce détecteur à une projection d'un motif quelconque du masque correspondant pour une distance prédéterminée du plan de coupe contenant des sources dans l'objet à reconstruire. Si l'on prend une source S1 dans un plan P1 parallèle au plan du détecteur, cette source étant située à l'intersection de l'axe OZ passant par le centre O du détecteur et perpendiculaire au plan P de ce détecteur, l'image sur le détecteur du motif M1 du masque M éclairé par la source S1 est représentée en 11. L'étendue de cette image représente la portion utile du détecteur, utilisée pour reconstruire les images des sources contenues dans le plan P1 situé à une distance prédéterminée du détecteur. Le champ maximal de la coupe C1 de l'objet qu'il est possible d'observer pour le plan P1 s'obtient en projetant à partir du centre O du détecteur, le même motif M1 du masque, sur le plan P1. Pour une autre source S2 du plan P1, qui n'est pas située sur l'axe OZ, on utilise la même portion I1 du détecteur, en projetant sur celui-ci un autre motif M2 du masque M. Si une source S3 est située sur l'axe OZ, mais dans un plan P2 plus éloigné du détecteur que le plan P1, l'image I3 du motif M1 éclairé par cette source S3 couvre une portion plus restreinte du détecteur, que l'image I1 produite par la source S1. Au contraire, si une source est située sur l'axe OZ, dans un plan plus rapproché du détecteur que le plan P1, l'image du motif Ml éclairé par cette source couvre une portion plus étendue du détecteur que L'image Il produite par la source S1. Si à partir du centre O du détecteur, on trace une droite quelconque coupant par exemple trois plans de coupe consécutifs, très rapprochés, respectivement en trois points, ces points ont sensiblement les mêmes coordonées (i, j) par rapport au repère (X, Y) de référence ; seule la distance z de ces points repérée sur OZ varie par rapport au détecteur. On recherche les maximums d'intensité des sources, par comptages, pour des sources qui sont situées dans des plans de coupes parallèles et voisins et qui sont sur une même droite passant par le centre géométrique du détecteur, ces sources ayant les mêmes coordonnées (i,j) dans chaque plan. A partir du centre O du détecteur et en s'appuyant sur un motif du masque, comme indiqué plus haut, on définit une enveloppe de l'objet. Toutes les sources doivent être situées dans cette enveloppe, mais pour chaque plan de coupe, la portion de détecteur est choisie de façon à correspondre à un motif quelconque du masque. Une gamma-caméra équipée d'un tel type de détecteur et utilisant un masque constitué par un assemblage de motifs, fournit sur des sorties, des valeurs numériques relatives aux énergies des rayonnements gamma reçus par le détecteur, à travers les éléments transparents du masque. Chaque sortie est codée en position et permet donc aussi de repérer les coordonnées, par rapport à un repère de référence, de l'impact d'un rayonnement sur le détecteur. Les valeurs numériques pour chaque position d'impact peuvent être enregistrées avant.de subir un traitement permettant la reconstitution des images des sources de chaque plan de coupe. Les valeurs traitées sont celles qui correspondent aux images des motifs, pour chaque éloignement de plan de coupe et pour toutes les coordonnées de sources émettrices qui peuvent être situées dans chaque plan. De manière connue, selon CANNON et FENIMORE, on décode pour tous les plans de coupes et pour toutes les coordonnées de sources émettrices qui peuvent être situées dans ces plans, toutes Les valeurs de comptages correspondant à toutes les images de tous les motifs du masque. Ce décodage est réalisé grâce à des matrices de décodage des images des motifs sur le détecteur. Ces matrices de décodage sont obtenues à partir des matrices de codage. Ce décodage permet de sélectionner dans les valeurs de comptages, celles qui correspondent aux images des motifs éclairés par les différentes sources de l'objet, dans les différents plans de coupes.Le traitement de ces valeurs qui dépendent de la position de chaque source dans un plan de coupe et de l'éloignement de ce plan, permet de reconstruire les sources contenues dans chaque plan de coupe et donc L'image de l'objet pour chaque plan de coupe. Selon L'invention, les éléments transparents de chaque motif du masque sont repérés par un numéro de ligne et un numéro de colonne. -Les motifs du masque sont de forme rectangulaire ou carrée, de même que les éléments transparents de chaque motif. Les numéros de ligne et de colonne de chaque élément transparent correspondent respectivement à deux nombres d'une succession de paires de nombres d'une série bidimensionnelle. Les paires de nombres de cette série sont obtenues à partir des nombres successifs d'une série pseudo-aléatoire monodimensionnelle, prédéterminée. Le tableau de la figure 2 donne dans la première colonne, un exemple de valeurs d'une série pseudo-aléatoire monodimensionnelle connue, de type q-naire, tabulée de O à 650, utlisée dans le procédé de l'invention, pour former une série pseudo-aléatoire bidimensionnelle. Pour construire une série pseudo-aléatoire bidimensionnelle A(i, j) à partir d'une série monodimensionnelle A(x), on écrit que le contenu binaire de A(i,j)=A(x) avec i=mod x et j=mod x. r s La série retenue est dite q-naire et a pour règles de formation : n n-l n-2 q -1 q -1 q -1 i m p = ---- ; p = = ------- ; l = -------- et q = k q-1 q-1 q-1 k est un nombre premier et i est un nombre entier, m représente le nombre d'éléments d'un motif du masque, p le nombre de trous ou éléments transparents, I la corrélation de la série par elle-meme. On prend m=651, q=5 et n=3. Les paires de valeurs de la série bidimensionnelle correspondant à chaque valeur de la série monodimensionnelle, sont données dans les deux colonnes suivantes i et j du tableau de la figure 2. C'est ainsi qu'à la valeur 1 de la série monodimensionnelle, correspondent les paires de valeurs i=l et j=l de la série bidimensionnelle. A la valeur 47 correspondent les valeurs i=16 et j=5 etc. En fait, chaque valeur i de la série bidimensionnelle correspond à la division modulo 31 de la valeur correspondante de la série monodimensionnelle. Chaque valeur J de la série bidimensionnelle correspond à la division modulo 21 de la valeur correspondante de la série monodimensionnelle. Les valeurs de i et J représentent les numéros de lignes et de colonnes à l'intersection desquelles, selon l'invention, on place les éléments transparents d'un motif, de forme rectangulaire ou carrée. Les autres parties du motif sont opaques. Un motif utilisé dans La mise en oeuvre du procédé de l'invention, est repésenté sur la figure 3. Les éléments transparents sont représentés par des surfaces carrées ou rectangulaires, de couleur noire. Les numéros i et J de lignes et de colonnes sont indiqués sur la figure. Pour constituer le motif de la figure 3, on n'a utilisé que les valeurs i comprises entre O et 20 et J comprises entre 0 et 30 dans le tableau de la figure 2. La dimension du masque est choisie de sorte que la plus grande projection d'un motif, s'inscrive dans le diamètre du détecteur. Des calculs simples montrent que le rapport S/s du sional au bruit. est de la forme Dans cette expression C est la valeur totale de comptages du détecteur, k est le rapport entre l'activité d'une source dans l'objet à L'activité totale des sources de l'objet (valeur du comptage pour un temps prédéterminé), m représente le nombre total d'éléments opaques d'un motif du masque, et p est le nombre d'éléments transparents (qui sont en fait des fenêtres ou trous).Cette relation montre qu'à comptage total constant, le rapport S/B augmente lorsque le nombre de trous diminue. Le masque de L'invention dans lequel le nombre d'éléments transparents est faible, fournit un rapport S/B qui est ç fois plus élevé que les masques dont la transparence est voisine de 50X, par exemple du type de ceux utilisés pour FENIMORE et CANNON et n'exige pas une homogénéité de réponse de la caméra, meilleure que le pourcent. Pour reconstruire les images de différentes coupes de L'objet, on opère selon l'invention, de la façon suivante, pour l'une de ces images : les valeurs de comptages fournies par la caméra sont enregistrées dans une mémoire ; ces valeurs de comptages correspondent aux rayonnements reçus par le détecteur, à travers les éléments transparents de l'un ou l'autre des motifs du masque, lorsque celui-ci est éclairé par les différentes sources des différents plans de coupe. Pour le traitement de ces valeurs enregistrées, on opère de la façon suivante. Les valeurs de comptages sont sélectionnées selon Les matrices de décodage déduites des matrices de codage des positions des éléments transparents d'un motif. Ces matrices de décodage correspondent aux images sur Le détecteur, d'au moins un motif du masque éclairé par Les sources contenues dans l'objet. Comme indiqué plus haut, les dimensions et les emplacements de ces images sur le détecteur dépendent des coordonnées des sources dans chaque plan de coupe, ainsi que de l'éloignement de chaque plan par rapport au détecteur. Dans le procédé de l'invention, on choisit par exemple des matrices de décodage correspondant aux images des motifs qui sont centrées sur l'axe du détecteur, ces images occupant une portion du détecteur, au centre de celui-ci. Il en résulte que pour obtenir des images occupant cette position centrale, il est nécessaire de considérer que l'un ou l'autre des motifs du masque, sont éclairés par les sources, les matrices de décodage étant choisies en conséquence. Le procédé consiste à effectuer les opérations suivantes : - comme indiqué plus haut, on sélectionne une portion du détecteur pour chaque éLoignement de plan de coupe Pi par rapport au détecteur, - pour chaque plan de coupe Pi, on déplace la matrice de décodage pour repérer Les sources évertuellement présentes dans ce plan Pi, ce déplacement correspondant à un balayage de toutes les sources repérées par leurs coordonnées (i,j) dans le plan Pi, et par des valeurs numériques mesurees qui sont enregistrées en mémoire, les sources repérées étant celles qui sont situées sur des droites passant par le centre géométrique du détecteur, - à partir des valeurs numériques correspondantes, on reconstruit les activités des sources de chaque plan.Les valeurs de ces activités sont enregistrées en mémoire, - pour chaque plan Pi d'éloignement predéterminé, on effectue une comparaison entre l'activité de chaque source de coordonnées (i,j) ainsi reconstruite dans le plan Pi, et les activités reconstruites respectives pour des sources voisines situées dans des plans voisins Pi+l et Pi-1, ayant les mêmes coordonnées (i,j) que la source considérée dans le plan Pi, - on détermine si les écarts entre l'activité de chaque source de coordonnées (i,j) dans chaque plan Pi et les activités respectives des sources de mêmes coordonnées (i,j) dans les plans Pi-l et Pi+l voisins dépassent un seuil prédéterminé, ce dépassement indiquant le repérage d'une source de coordonnées (i,j) dans le plan Pi, - on attribue une activité à chaque source ainsi repérée, cette activité ainsi que les coordonnées de la source sont enregistrées en mémoire, - on calcule pour chaque source repérée la valeur de comptages que produirait cette source à partir de L'activité qui lui a été attribuée. Cette valeur calculée est enregistrée en mémoire, - on soustrait, pour chaque source repérée, La valeur de comptages effectivement mesurée précédemment au cours du balayage (pour la source de mêmes coordonnées (i,j) dans le même plan), de la valeur de comptages calculée pour la source repérée, de manière à obtenir une valeur de comptages corrigée pour la source repérée.Cette valeur corrigée est enregistree en mémoire, - enfin, on détermine 13activité de chaque source repérée à partir de la valeur de comptages corrigée correspondante et on recommence les traitements précédents selon un mode itératif jusqu'à ce que la valeur corrigée obtenue pour chaque source soit invariable. C'est à partir des valeurs corrigées invariables que sont obtenues les images des coupes de l'objet. Le traitement itératif cesse en fait lorsque le rapport signal/bruit atteint une valeur suffisante, c'est-à-dire lorsque le bruit de fond est devenu peu significatif. La figure 4 représente schématiquement un système de reconstruction d'images de coupes parallèles successives d'un objet ou d'un organe contenant des sources émettrices de rayonnement gamma et permettant de mettre en oeuvre le procédé de l'invention. Ce système comprend une gamma-caméra 5, située en regard de l'objet ou de l'organe 2 dont on veut reconstituer des images de coupes parallèles successives. Cette gamma-caméra comprend notamment un détecteur 1 (sans collimateur) qui n'est pas représenté en détail sur la figure, mais qui comprend de manière connue, un scintillateur associé à des photomultiplicateurs. Ces photomultipLicateurs sont reliés å des moyens d'alimentation, d'amplification et de codage 6, qui fournissent sur une pluralité de sortie 7, une pluralité de valeurs codées de l'énergie des rayonnements reçus par le détecteur. Ces valeurs sont codées de manière qu'à chaque valeur correspondent des informations relatives aux coordonnées du point d'impact de chaque rayonnement sur le détecteur. Ces coordonnées sont repérées dans le repère de référence mentionné plus haut.On a aussi représenté sur cette figure le masque M comprenant un assemblage de motifs tels que MI, décrit plus haut. Ce masque est interposé entre le détecteur 1 et L'objet ou L'organe 2. Dans ce système, la gamma-caméra ne comporte pas de collimateur mais est munie d'un limiteur 8 de champ d'observation. Les valeurs codées fournies par les sorties 7 de la caméra sont enregistrées dans une mémoire 9. Les sorties de cette mémoire sont reliées à des moyens de décodage 10 des valeurs codées mémorisées. Ces moyens de décodage sont reliés à une unité de traitement 11 qui leur applique les matrices de décodage mentionnées plus haut. Les moyens de traitement reçoivent les valeurs ainsi sélectionnées selon ces matrices.Les valeurs décodées peuvent être enregistrées dans un ensemble de mémorisation 12, relié aux moyens de traitement 11, avant d'être traitées pour être ensuite transmises à des moyens de visualisation 13 des images reconstituées des différents plans de coupes de l'objet ou de l'organe 2. Les moyens de mémorisation 12 permettent aussi d'enregistrer les programmes définissant les matrices de décodage des valeurs codées et mémorisées dans La mémoire 9. Les moyens de mémorisation 12 peuvent aussi enregistrer les programmes de traitement des valeurs décodées ; ces programmes ne sont pas décrits ici en détail. Ils mettent en oeuvre un organigramme conforme aux étapes du procédé décrit plus haut. Le dispositif donne un rapport signal sur bruit et une résolution d'autant meilleurs que la source est plus superficielle. Une amélioration sensible de la détection des sources profondes, consiste à prendre sur le détecteur deux projections d'un motif du masque aussi écartées que possible. On obtient deux avantages : l'un dû à L'obliquité de rayonnements gamma participant à la reconstruction d'image, l'autre lié à la statistique de comptage en augmentant la surface de détection utilisée. REVENDICATIONS 1. Procédé de reconstruction d'images de coupes parallèles successives d'un objet (2) contenant des sources émettrices de rayonnements gamma, à L'aide d'une gamma-camera (5) à détecteur plan (1) parallèle aux plans de coupes, cette caméra ayant une pluralité de sorties (7) repérées chacune par un code et correspondant respectivement aux différentes coordonnées d'impacts possibles et à L'énergie des rayonnements gamma reçus par le détecteur, cette énergie étant mesurée par une valeur de comptages fournie par la sortie correspondant à un impact, le procédé consistant à interposer entre l'objet et le détecteur, parallèlement à ce détecteur, un masque (M) plan formé par assemblage de motifs (M1) identiques, de forme rectangulaire ou carrée et comportant chacun une mosaique d'éléments transparents dans une surface opaque aux rayonnements, ces éléments transparents occupant dans le motif des positions prédéterminées repérées selon une matrice de codage, la reconstruction des images des coupes étant effectuée par un traitement des valeurs de comptages enregistrées dans une mémoire (9), ces valeurs étant sélectionnées selon des matrices de décodage déduites des matrices de codage et correspondant aux images sur le détecteur (1) d'au moins un motif (M1) du masque éclairé par lesdites sources, les dimensions et les emplacements de ces images sur Le détecteur (1) dépendant du motif choisi dans le masque, des coordonnées des sources dans chaque plan, par rapport à un repère de référence, ainsi que de L'éloignement de chaque plan par rapport au détecteur, le procédé consistant à effectuer les opérations suivantes :: - on sélectionne une portion du détecteur pour chaque éloignement de plan de coupe Pi par rapport au détecteur, - pour chaque plan de coupe Pi, on déplace la matrice de décodage pour repérer les sources éventuellement présentes dans ce plan, ce déplacement correspondant à un balayage de toutes les sources repérées par leurs coordonnées (i, j) dans le plan Pi, et par des valeurs de comptages mesurées, les sources repérées etant celles qui sont situées sur des droites passant par le centre géométrique du détecteur, à à partir des valeurs de comptages correspondantes, on reconstruit les activités des sources de chaque plan, caractérisé en ce que l'on effectue alors les traitements suivants :: - pour chaque plan Pi d'éloignement predéterminé, on effectue une comparaison entre L'activité de chaque source de coordonnées (i, j) ainsi reconstruite dans le plan Pi, et les activites reconstruites respectives pour des sources voisines situées dans des plans voisins Pi+l et Pi-i, ayant les mêmes coordonnées (i, j) que la source considérée dans le plan Pi, - on détermine si les écarts entre L'activité de chaque source de coordonnées (i, j) dans chaque plan Pi et les activités respectives des sources de mêmes coordonnées (i, j) dans les plans Pi-l et Pi+l voisins dépassent un seuil prédéterminé, ce dépassement indiquant le repérage d'une source de coordonnées (i, j) dans le plan Pi, - on attribue une activité à chaque source repérée, - on calcule pour chaque source repérée la valeur de comptages que produirait cette source à partir de l'activité qui lui a été attribuée, - on soustrait, pour chaque source repérée, la valeur de comptages effectivement mesurée précédemment pour la source de mêmes coordonnées (i, j) dans le même plan, de la valeur de comptages calculée pour la source repérée, de manière à obtenir une valeur de comptages corrigée pour la source repérée, - on détermine L'activité de chaque source repérée à partir de la valeur de comptages corrigée correspondante et on recommence les traitements précédents seLon un mode itératif jusqu a ce que la valeur corrigée obtenue pour chaque source soit invariable. 2. Procédé selon la revendication 1, caractérisé en ce que chacun desdits éléments transparents d'un motif (M1) est repéré par un numéro de ligne et un numéro de colonne correspondant respectivement aux deux nombres d'une paire de nombres choisie dans une succession de paires de nombres d'une série bidimensionnelle équivalente à une série monodimensionnelle pseudo-aléatoire.