FR 2553898 A1 19850426 FR 8316995 A 19831025 PROCEDE D'ACQUISITION ET DE RECONSTRUCTION D'IMAGE PAR TOMODENSITOMETRIE ET TOMODENSITOMETRE METTANT EN OEUVRE CE PROCEDE L'invention concerne un procédé de reconstruction d'image par tomodensitométrie ayant plus particulièrement pour but d'améliorer la définition de l'image reconstruite -sans alourdir la structure technologique des moyens de mesure de l'absorption de rayonnement pénétrant L'invention concerne également tout tomodensitomètre adapté pour mettre en oeuvre ce procédé. Les tomodensitomètres dits de "troisième génération" se composent essentiellement d'une source de rayonnement pénétrant (des rayons X) émettant un faisceau plat en éventaill d'une rangée courbe 10 de détecteurs places en regard de cette source pour mesurer la fraction non absorbée du rayonnement dans le plan de la coupe à imager, des moyens pour faire tourner la source et les détecteurs dans le plan de cette coupe et un ordinateur pour traiter les données fournies par les détecteurs et reconstruire une image à partir de ces données La définition de l'image reconstruite dépend notamment du nombre des détecteurs embrassant l'éventail du faisceau de rayons X Pour une ouverture de faisceau donné, plus le nombre de détecteurs est élevé et plus la définition est bonne Cependant, l'ensemble multidétecteurs est l'un des organes les plus co Oteux du 20 tomodensitomètre et l'un de ceux qui est le plus difficile à fabriquer avec une bonne fiabilité Par conséquent, l'augmentation du nombre de détecteurs sur un secteur donné correspondant à l'éventail, se heurte à des limites technologiques. L'un des buts de l'invention est d'améliorer sensiblement la 25 définition de l'image sans modifier la structure technologique de l'ensemble source-détecteurs. Un autre but de l'invention est aussi de réduire le coût de construction d'un tomodensitomètre ayant la même définition que les tomodensitomètres actuellement existants, cette réduction de 30 coût étant obtenue en utilisant un ensemble multidétecteurs com- coût étant obtenue en utilisant un ensemble multidétecteurs comportant deux fois moins de détecteurs que les machines existantes. Dans ce dernier cas, non seulement l'ensemble multidétecteurs serait plus facile et moins coûteux à fabriquer, mais la réduction de 5 coût s'étendrait aussi aux moyens électroniques d'acquisition des données. Selon l'invention, ces buts sont atteints grâce à un procédé de reconstruction d'image par tomodensitométrie, du type consistant à faire tourner un ensemble source-détecteurs dans le plan d'une 10 coupe à imager, ladite source émettant un faisceau de rayonnement pénétrant en forme d'éventail et à relever et mémoriser une succession de vues pour différents angles / de rotation de l'ensemble source-détecteurs par rapport à un système d'axes de référence fixe dudit plan, lesdits détecteurs étant agencés suivant un arc de cercle centré sur le foyer de ladite source et chaque vue étant échantillonnée pour des angles Y représentatifs des détecteurs et repérés à partir de la droite joignant le centre de rotation de l'ensemble source-détecteurs et ledit foyer, le décalage angulaire entre deux détecteurs adjacents étant L\ ' constant et le décalage 20 angulaire entre vues étant constant, caractérisé en ce qu'il consiste, le décalage angulaire entre ladite droite et le détecteur le plus proche de celle-ci étant réglé à / et les conditions suivantes étant réalisées: -Y p ( 1 + 0)o(es /'7 e = N, avec N entier positif, en posant e = ( +) o est 25 un coefficient dépendant de la vitesse de rotation de l'ensemble source-détecteurs et du temps d'acquisition d'une vue; à relever et mémoriser deux séries de vues avec un décalage angulaire entre vues de h LS pour la première série et Ap pour la seconde série, au moins certaines des vues des deux séries étant A e décalées entre elles de l 21 ó/, avec E /= 2/e, h étant entier positif; à compléter les vues de la première série par des valeurs échantillonnées de la seconde série dont les coordonnées (e', ') répérées dans l'espace 'l/ 4 e avec O = P 3 + e satisfont à la relation: (O x, x)= (O, + ( 2 k+ I), -t') avec k entier relatif ( x, 'x) étant les coordonnées dans l'espace -, 'e des valeurs échantillonnées recherchées de la première série, et à appliquer un algorithme de reconstruction, connu en soi, comprenant une convolution et une rétro-projection, aux vues complétées de la première série, pour reconstituer une image. L'invention concerne aussi un tomodensitomètre du type com10 portant un ensemble source-détecteurs monté tournant dans un plan de coupe à imager, ladite source étant une source de rayonnement pénétrant en forme d'éventail et lesdits détecteurs étant agencés pour converger vers le foyer de ladite source, lesdits détecteurs étant vus dudit foyer suivant des angles ' respectifs repérés à 15 partir de la droite joignant ledit foyer et le centre de rotation de l'ensemble source-détecteurs, deux détecteurs voisins étant angulairement décalés de 14, des moyens de pilotage pour repérer un certain nombre d'angles/3 de prises de vues de l'ensemble source-détecteurs par rapport à un 20 système d'axe de référence et pour engendrer des signaux de commande desdites prises de vues, des moyens de mémorisation d'informations représentatives desdites vues échantillonnées pour des angles / précités, des moyens de calcul pour appliquer un algorithme de recons25 truction, connu en soi, comprenant une convolution et une rétroprojection à une série de vues, caractérisé en ce que: les conditions suivantes sont satisfaites: = N, avec N entier positif,en posant te = ( 1 +k)>Y o d est / Yeeun coefficient dépendant de la vitesse de rotation de l'ensemble source- détecteurs et du temps d'acquisition d'une vue; un décalage angulaire entre ladite droite et le détecteur le plus proche de celle-ci est prédéterminé à l-Y les moyens de pilotage comportent des moyens pour commander la prise et la mémorisation de deux séries de vues échantillonnées, lesdites séries étant telles qu'au moins certaines vues voisines de séries différentes soient décalées entre elles de /I /; les vues de la première série étant décalées angulairement au pas de h/"et celles de la seconde série étant décalées angulairement au pas de Ap, avec A = 2 ke et h entier positif, et en ce qu'il comporte des moyens de réarrangement pour compléter les vues de la première série par des valeurs échantillonnées prélevées dans la seconde série dont les coordonnées (t,)) dans l'espacez 5,e J avec =p 3 + e satisfont à la relation ( ex,' = ( O ' + ( 2 k+ 1)i, _ ') avec k entier relatif, (e x, x) étant les coordonnées dans l'espace, des valeurs échantillonnées recherchées de la première série. L'invention sera mieux comprise et d'autres avantages de 15 celle-ci apparaîtront mieux à la lumière de la description qui va suivre d'un procédé conforme à son principe et des éléments-constitutifs essentiels d'un tomodensitomètre fonctionnant selon ce procédé, faite en référence aux dessins annexés dans lesquels: la figure I est une représentation schématique des principaux 20 éléments constitutifs d'un tomodensitomètre mettant en oeuvre l'invention; la figure 2 est une représentation dans l'espace l,l J des vues prises par un tomodensitomètre classique fonctionnant en mode pulsé pour la reconstruction d'une image; la figure 3 est une représentation dans le même espace lID Jdes vues prises par un tomodensitomètre du même type dans lequel on a effectué l'une des modifications nécessaires à la mise en oeuvre de l'invention, à savoir un décalage statique entre la source et les détecteurs; la figure 4 est une représentation dans l'espace l>, des vues prises par un tomodensitomètre fonctionnant selon le processus de l'invention; la figure 5 est une représentation dans l'espace Yldu résultat d'une transformation, faisant partie de l'invention, appliquée à une vue; la figure 6 illustre le résultat de l'application de cette transformation à certaines vues de la figure 4; la figure 7 illustre dans l'espacer C,lle résultat de la même transformation appliquée selon une variante de l'invention; la figure 8 est une représentation dans l'espacefz,X d'une vue, dans le cas d'un tomodensitomètre fonctionnant en mode continu. La suite de la description, en référence aux figures I à 7, se place dans le cas particulier d'un tomodensitomètre fonctionnant en mode pulsé Ce genre de tomodensitomètre est bien connu Dans la 10 plupart des cas, le mode pulsé se concrétise par le fait que le faisceau en éventail n'est émis que pendant un intervalle de temps très bref au moment o l'ensemble source-détecteurs se trouve à un angles prédéterminé Dans ce cas, il suffit que la lecture de tous les détecteurs soit achevée lorsque l'ensemble source-détecteurs parvient à l'angle/ suivant Cependant, une autre façon de réaliser ledit mode pulsé peut consister à générer en permanence le faisceau en éventail mais à "lire" tous les détecteurs simultanément à chaque angle/ Dans le cas d'un fonctionnement en mode pulsé et comme on le vérifiera plus loin, le coefficient Se mentionné cidessus est 20 égal à Y et Le =m XEn se reportant à la figure 1, on a représenté les éléments constitutifs essentiels d'un tomodensitomètre Cet appareil comporte un ensemble source-détecteurs mobile en rotation dans un plan d'une coupe à imager Le centre de rotation est O et le plan de la 25 coupe sera considéré comme étant celui de la figure I On associe à ce plan de coupe un système d'axes de référence fixe XOY, orthonormé La source est symbolisée sur le dessin par son foyer F De façon bien connue, il s'agit d'une source de rayonnement pénétrant et plus particulièrement des rayons X, émettant un faisceau 11 en 30 forme d'éventail Dans la suite du texte, on parlera simplement de la source F Le boîtier 13 renfermant tous les détecteurs 14 a une forme courbe pour que ces détecteurs 14 soient agencés suivant un arc de cercle CC' géométriquement centré sur le foyer F Le faisceau en éventail 11 embrasse l'arc de cercle CC' occupé par les faces d'entrée des détecteurs 14 Ainsi, l'image reconstituée sera à l'intérieur du cercle 16 de centre 0, inscrit dans les limites de l'angle CFC' Classiquement, on peut utiliser un ensemble comportant I 024 détecteurs 14 Le boîtier 13 et la source F sont solidaires d'un même support tournant (centre de rotation O) et fixes l'un par rapport à l'autre de sorte qu'on peut repérer la rotation de l'ensemble source-détecteurs par l'angle /3 entre 'axe OY et la droite OF Les valeurs d'atténuation linéraires du faisceau de rayons X pour des angles 19 prédéterminés sont relevées et mémorisées 10 dans une mémoire M d'un ordinateur Pour chaque angle /l prédéterminé, un ordre de lecture des détecteurs 14 est émis, par exemple au moyen d'un codeur incrémental 15 agencé pour lire des informations d'angles solidaires de l'ensemble tournant sourcedétecteurs et pour délivrer (liaison de pilotage 17) des ordres de 15 commande de lecture de l'état de tous les détecteurs 14 L'ensemble des valeurs relevées pour un angle f donné est appelé "vue". On peut donc considérer que chaque vue est échantillonnée de par la structure même de l'ensemble de détecteurs, pour des angles Y représentatifs des détecteurs Chaque angle Y indiquant la 20 position d'un détecteur 14 est repéré à partir de la droite OF De ce fait on appelera IY le décalage angulaire constant entre les centres des faces d'entrée de deux détecteurs 14 adjacents. Si on considère, pour chaque angle Y d'une vue, l'axe O O perpendiculaire en P au rayon d'angle Y frappant un détecteur 25 quelconque, on démontre facilement que l'angle O entre l'axe OX du repère fixe et l'axe OX est tel que: e=p +y ceci dans tous les cas de figures, les angles, p et ' étant orientés, ou encore Y = 9 p D'après cette dernière relation, on voit que les différentes vues (pour des angles p différents) ont une représentation extrêmement simple dans l'espace Le J Comme le montre la figure 2, ces vues sont portées par des droites parallèles à la première bissectrice et décalées de Arle long de l'axe des On peut porter sur ces droites des points représentant les coordonnées, dans l'espace lO, Yl des valeurs échan5 tillonnées effectivement acquises à la lecture des différentes vues. Ainsi, si on suppose que le nombre de détecteurs est pair, soit 2 n, et que la droite OF passe exactement entre les deux détecteurs du centre, 'n et Yn + ceux-ci apparaissent décalés de I / de part et d'autre de l'axe Y Le but de l'invention est d'obtenir pour chaque vue, des valeurs échantillonnées intermédiaires dont les coordonnées dans l'espace le,Yr sont figurés par des étoiles* sur la figure 2 Autrement dit, le but de l'invention est d'obtenir des vues comportant 4 N valeurs échantillonnées avec un ensemble de 2 N détecteurs seulement et d'améliorer ainsi la définition de l'image 15 reconstruite. Avant de définir l'invention sous une formulation plus générale, on va décrire complètement une solution particulière en fixant à priori un certain nombre de paramètres et on expliquera ensuite comment on peut faire varier certains d'entre eux pour définir un 20 ensemble de solutions possibles. Selon une possibilité de l'invention, on provoque avant tout un décalage angulaire permanent et prédéterminé entre l'ensemble des détecteurs 14 (c'est-à-dire le boîtier 13) et la droite OF Selon l'exemple choisi, le décalage angulaire entre la droite OF et le 25 détecteur le plus proche est choisi égal à Cela signifie que l'ensemble des détecteurs est décalé dans le sens inverse des aiguilles d'une montre en considérant la figure 1, de façon que le centre de la face d'entrée du détecteur Yn soit décalé angulairement de AY par rapport à OF et que le détecteur Yn+ I soit 30 décale de 3 par rapport à cette même droite. La figure 3 montre la configuration des vues dans l'espace lY, 'l après ce décalage. Selon l'exemple décrit, on relève et mémorise deux séries Vi et Wi de vues et la machine est agencée pour que: = N; N étant entier positif lx= 2 LY Dans le cas o chaque série de vues est répartie sur un tour de l'ensemble source-détecteurs, alors N représente aussi le nombre de vues des séries -On suppose pour l'exemple que N est pair Dans ces conditions, on décale les vues Wi de la seconde série de + par rapport aux vues Vi de la première série. La figure 4 montre la représentation des vues Vi et Wi dans 10 l'espace lo,îj. Soit la transformation: T(O,Y)=(O -, Y) On peut démontrer que la fonction d'absorption d'une vue est invariante par cette transformation T Or cette transformation T 15 peut se décomposer dans l'espace lt Aï en une translation -_( suivied'une symétrie par rapport à O Cette transformation est illustrée en soi à la figure 5 Si on applique ladite transformation à la série des vues Wi, on constate (figure 6) que les valeurs échantillonnées des vues W i ainsi transformées prennent automa20 corone tiquement les coordonnées des valeurs échantillonnées intermédiaires recherchées, dans la série des vues Vi Autrement dit, l'application de cette transformation aux valeurs échantillonnées des vues de la seconde série a pour conséquence de compléter les vues de la première série Chaque vue Vi complétée est donc deux fois 25 plus échantillonnées que la vue Vi d'origine. On peut alors en déduire une étape de réarrangement des données traduisant cette transformation Selon ce réarrangement, on complète les vues de la première série par des valeurs échantillonnées de la seconde série dont les coordonnées (O', X ') dans 30 l'espace l 0, J satisfont à la relation (O x, X x) = (O + ( 2 k + 1),)i ') k étant un entier relatif, (Ox, x) étant les coordonnées, dans l'espace f, Yl des valeurs échantillonnées recherchées de la première série. L'expression ( 2 k + l)l Ttraduit le fait que les vues réellement acquises Wi de la deuxième série peuvent être théoriquement I prolongées dans l'espace L Pe yj de - à + En raison de la périodicité des vues au- delà d'un tour complet, qui peut se traduire ainsi: Wi Wj i = j modulo N Autrement dit, même si la translation 9 ' + ( 2 k + 1)'rne désigne pas une vue Wi réellement acquise, il suffit de rechercher k pour que: Px = 1 + ( 2 k + l)Tsoit vérifiée avec O ' désignant une vue 10 réellement acquise de la seconde série W I Une fois le réarrangement des données, tel que décrit cidessus, effectué, il ne reste plus qu'à appliquer un algorithme de reconstruction connu en soi, comprenant une convolution et une rétroprojection aux vues complétées de la première série pour 15 reconstruire une image de définition améliorée En se reportant à nouveau à la figure 1, on a décrit sous forme de schéma-bloc fonctionnel les parties essentielles de l'ordinateur associé aux détecteurs 14 Les valeurs échantillonnées lues aux sorties des détecteurs 14 sont transformées en informations numériques grâce à un 20 convertisseur analogique-numérique 18 puis traitées dans un préprocesseur 19 effectuant les opérations habituelles de calibration et de transformation logarithmique Les informations ainsi transformées sont les valeurs d'atténuation linéaires et ces informations numériques sont mémorisées dans la mémoire M Celle-ci est divisée 25 en une mémoire A recevant dans un premier temps des informations numériques représentatives des vues Vi de la première série et en une mémoire B recevant des informations numériques représentatives des vues Wi de la seconde série Selon ce qui précède, cette dernière se présente sous forme d'une matrice d'unités de mémori30 sation comprenant 2 N colonnes (nombre de détecteurs par vues) et N lignes (nombre de vues Wi) dans le cas o l'acquisition de la seconde série de vues W se fait sur un tour complet et avec /Wj N * Par ailleurs, pour des raisons qui apparaîtront clairement plus loin, la mémoire A se présente sous forme d'une matrice d'unités de mémorisation comprenant 4 N colonnes et N lignes (nombre de vues Vi) puisque l'exemple décrit correspond aussi à une acquisition des vues Vi sur un tour complet et avec V-N Le remplissage de la mémoire A par les informations représentatives des vues Vi se fait à raison d'une unité de mémorisation sur 2; ces unités sont représentées en grisé sur la figure 1 La mémoire A est associée à des moyens de lecture 20 qui organisent le transfert des informations vers un processeur 21, connu en soi, et dont le rôle est d'appliquer l'algorithme de reconstruction mentionné plus haut Cependant la lecture d'une ligne de la mémoire A n'est commandée que lorsque toutes les unités de mémorisation de cette 15 ligne qui n'ont pas été "remplies" pendant la phase de mémorisation des vues Vi, auront reçu des informations numériques lues dans la mémoire B. L'algorithme qui va être décrit maintenant et qui découle de la relation ( O x, i' x) = ( &' + ( 2 Ik + 1)', Y ') permet d'organiser le 20 transfert des informations depuis la mémoire B vers la mémoire A en vue de la compléter Cette relation a été exprimée jusqu'ici dans l'espace l'g Il pour faciliter la démonstration car les vues V ou W i I I y ont une représentation particulièrement simple Cependant, les 25 paramètres pris en compte pour l'adressage des mémoires au moment de la lecture des vues Viet Wi sont en fait: È: la position angulaire de l'ensemble source-détecteurs, et Y: représentatif des détecteurs. On montre que la relation simple: O = 3 +/ permet de transcrire la relation qui précède dans l'espace 3, de la façon suivante: (( 2 k + 1) + x + 2 x 'x)P' ( x) et ( Yx) étant les coordonnées dans l'espace j,' des valeurs échantillonnées manquantes des vues Vi et (Pl) et (g') étant les coordonnées des valeurs échantillonnées des vues Wi Pour simplifier l'algorithme, on suppose que N est pair, soit N = 2 K = 2-r, avec K entier positif, mais il est à noter que N impair ne constitue pas un obstacle majeur à la définition d'un algorithme semblable dont la mise en oeuvre aboutit au même résultat. Le problème se résume à un simple calcul d'adresse pour aller chercher dans la mémoire B les informations numériques représentatives des valeurs échantillonnées recherchées et pour recopier ces informations dans la mémoire A. iéme Soit la iime vue de la série Vi, celle-ci se trouve mémorisée 10 toute entière dans la iième ligne de la mémoire-matrice A mais les valeurs échantillonnées à compléter ont une adresse de la forme: A (i A, j A = 2 m -1) avec m e 1,2 nl On déduit de la relation précédente que le numéro i B de la ligne de la mémoire-matrice B o se trouve la valeur échantillonnée 15 à transférer à l'adresse A (i A, j A) est i B = (K + m + i A N 1) modulo N avec i Bî,NlJ Autrement dit, si l'expression (K + m + i A N 1) Ll,N, on retranche ou on ajoute N autant de fois que nécessaire pour ramener ladite expression dans l'intervalle l 1,Nl. Une fois la ligne i B déterminée, le numéro de la colonne j B permettant de sélectionner la valeur échantillonnée recherchée est j B = ( 2 N + I m) puis il suffit de transférer l'information contenue à l'adresse B (i B, j B) pour l'inscrire à l'adresse A (i A, j A) Les opérations de calcul d'adresse et de transfert d'information peuvent 25 être pilotées par un processeur spécialisé ou de préférence, par un sous-programme de l'ordinateur Ce processeur ou ce sousprogramme est symbolisé à la figure 1 par le bloc fonctionnel 22 interconnecté entre les mémoires A et B Il pilote aussi les moyens de lecture 20 (liaison fonctionnelle 23) gérant le transfert des infor30 mations depuis la mémoire A vers le processeur de reconstruction 21. L'invention a été décrite ci-dessus en fixant pour simplifier certains paramètres On va maintenant envisager d'autres variantes possibles. Ainsi, on s'est placé dans le cas o 12 i R= 2 t X avec = N Dans ce cas, il y a autant de vues dans les deux séries et, de toute façon, les vues de la seconde série doivent être prises dans ces conditions Il se peut cependant que le nombre de vues de la première série soit supérieur à ce qui est normalement nécessaire pour reconstituer une image Dans ce cas une premiere généralisation de l'invention consistera à relever et à mémoriser deux séries de vues avec un décalage entre vues de h Lr 3 pour la première 10 série Viet 1 pour la seconde série Wi Dans ce cas, il apparaît que seulement certaines vues des deux séries seront décalées entre elles de Y puisque des vues de la seconde série n'ont pas de vues 'voisines" de la première série, à cause de la différence de pas angulaires La figure 7 illustre cette situation, dans le cas o h = 2. 15 On voit que W 2 i + 1 n'a pas de vue "voisine" de la première série décalée de Y 2 D'autre part, la figure 7 montre que des valeurs échantillonnées de la série Wi sont inutiles après la transformation T si h est différent de 1 Pour optimiser le réarrangement, on pourra donc ne relever et/ou traiter qu'une valeur échantillonnée sur h dans 20 chaque vue de la série Wi Ceci compense la contrainte de devoir relever les vues Wi avec un pas angulaire 2 h. Une autre situation dans laquelle seulement certaines vues des deux séries seront décalées de est celle d'une acquisition "optimisée" des vues de chaque série, en moins d'un tour par série On 25 sait en effet que le nombre de vues nécessaire pour reconstruire une image peut être acquis en un peu plus d'un demi-tour de l'ensemble source-détecteurs, à savoir un demi-tour augmenté du secteur angulaire de l'éventail 11 Il existe des algorithmes de reconstruction adaptés à ce mode de prise de données Si on applique le 30 principe de l'invention à un tomodensitomètre utilisant un tel algorithme, les deux séries de vues n'auront qu'un secteur de rotation en commun (la première série étant par exemple prélevée seule sur un premier secteur angulaire puis des vues des deux séries étant prélevées sur un second secteur angulaire commun et la fin de la seconde série étant prélevée seule sur un troisième secteur angulaire) et bien sur, dans un tel cas, le décalage entre vue de 2 tel que défini plus haut, n'a de sens que pour le second secteur angulaire donc pour une partie des vues des deux séries. Par ailleurs, d'autres valeurs de décalage angulaire (des détecteurs par rapport à la source d'une part et entre des vues des deux séries d'autre part) sont possibles. Si on considère l'ensemble des détecteurs de la figure 1, soit 2 N détecteurs, on a analysé en détail le cas d'un décalage des 10 détecteurs tels quef N = -A-4 & avec N pair et un décalage de + entre vues des deux séries, ce qui signifie que la vue Wi sera prise à P = + si la vue Vi a été prise à /3 = O, p étant orienté positivement dans le sens trigonométrique On démontre que le résultat serait le même si le décalage entre vues était de 3 D'autre part, si N est impair, les deux valeurs possibles du décalage angulaire entre les vues sont 3 _ ou / Par ailleurs le décalage angulaire entre les détecteurs et la source F peut aussi être tel que Tn 3 c'est-à-dire n+ + 4 Dans ces conditions, si N est pair, le décalage angulaire entre vues peut être de 2 ou + 3 et si N est impair, le décalage angulaire entre vues peut être de ±5 ou Or un décalage angulaire de 32 par rapport à une vue de la première série équivaut à un décalage de + par rapport à la 25 vue précédente de la même série et un décalage de + 32 équivaut de la même façon à un décalage de 2 par rapport à la vue suivante. La généralité de la solution peut donc s'exprimer de la façon suivante: le décalage angulaire entre la droite OF et le détecteur le plus proche de celle-ci est réglé à a en valeur absolue soit /l le décalage angulaire relatif entre les vues des deux séries Viet W est de A 2 en valeur absolue, soit,r/ Sur la plupart des tomodensitomètres, il existe des moyens de règlage précis de la position de l'ensemble des détecteurs (c'est-àdire le boîtier 13) par rapport à la source Dans la plupart des cas, ce "vernier" est ajusté pendant la phase de mise au point de 5 l'appareil pour annuler aussi précisément que possible tout décalage entre les détecteurs et la source, c'est-à-dire: et 1 = Y n:-2 n+I 2 Les mêmes moyens peuvent être utilisés pour provoquer le décalage statique des détecteurs à / En ce qui concerne le décalage dynamique entre les vues des séries Vi et Wi, on pourra procéder de différentes manières On pourra par exemple acquérir la première série pendant un premier tour et la seconde série pendant un second tour, de l'ensemble source-détecteur Une façon simple d'obtenir le décalage voulu 15 entre le premier et le deuxième tour consiste à utiliser un codeur incrémental dont le pas élémentaire de comptage soit -(ou un sous- multiple de cette valeur), c'est-à dire qu'il fournisse une impulsion exploitable tous les -t Deux vues successives d'une série seraient ainsi déclenchées toutes les 4 h impulsions du codeur 20 incrémental pour la série Vi et toutes les 4 impulsions pour la série Wi et le décalage dynamique à la fin d'un tour se ferait par exemple "en sautant" une impulsion du codeur incrémental avant d'acquérir la première vue de la seconde série Les moyens de traitement logique à associer au codeur incrémental pour réaliser cette fonction de 25 décalage dynamique sont d'une très grande simplicité et leur conception est à la portée de l'homme du métier Si la rapidité des circuits électroniques le permet (ou éventuellement en adoptant une vitesse de rotation plus faible del'ensemble source-détecteurs) il est aussi parfaitement possible de procéder à une acquisition "entre30 lacée" des vues des deux séries Dans le cas particulier décrit cidessus, par exemple, on pourrait acquérir les vues de la façon suivante: -la vue Vi serait acquise à /3 = i 1 -la vue Wi serait acquise à = i + + i 2 la vue Vi+ 1 serait acquise à/ = (i + 1) A la vue Wi + 1 serait acquiseà p = (i + 1)A + 'ê i+ 1 2 et ainsi de suite Bien entendu, dans ce dernier cas, le programme de l'ordinateur serait modifié pour écrire les données alternativement dans les mémoires A et B. On va maintenant généraliser le principe de l'invention à un tomodensitomètre fonctionnant en "mode continu" Ce type de tomodensitomètre est également connu et il se distingue essentiellement du tomodensitomètre fonctionnant en mode pulsé en ce que le faisceau en éventail est généré en permanence et en ce que la 15 lecture des détecteurs à chaque prise de "vue" se fait par scrutation progressive, d'une extrémité à l'autre de l'ensemble multidétecteurs, par exemple depuis le détecteur '1 jusqu'au détecteur Y 2 n. Deux paramètres supplémentaires doivent donc être pris en considération: la vitesse de rotation, v, de l'ensemble source-détecteurs; l'intervalle de temps Z séparant l'acquisition des valeurs échantillonnées correspondant à deux détecteurs voisins, en considérant une vitesse de scrutation constante. La figure 8 illuste la façon dont la représentation d'une vue 25 dans l'espace re, évolue si on prend en compte les paramètres v et Z La vue portée par la droite en trait interrompu correspond au mode pulsé tandis que la même vue portée par la droite en trait plein correspond au mode continu. Les valeurs échantillonnées se déplacent parallèlement à l'axe e puisque O = 3 +; P variant pendant l'acquisition Les deux droites concourrent en Y 1 puisque la scrutation des détecteurs commence par le détecteur Y 1 et la valeur échantillonnée suivante, Y 2 ' se trouve déplacée en ',2 a d'une distance v Z dans l'espace Q Jparallèlement à l'axe des Soit â l'angle entre 1 12 et I '2 a' Du fait que '1 et X 2 sont portés par une droite parallèle à la première bissectrice, la pente de la droite I 2 a est: tga A+v +v soit v tga= c L Les vues sont donc portées par des droites parallèles à la droite d'équation O = ( 1 +d L)'Y. Si on pose (I +o e)>'= e on constate que toutes les opérations et transformations qui ont été décrites en référence au mode pulsé dans l'espace o sont conservées dans l'espace lt, 'elpour le mode continu Le mode pulsé peut même être considéré comme un cas limite du mode continu puisqu'une émission pulsée de la source 15 peut être considérée comme réalisant des prises de "vues" à vitesse nulle (v = O, donc o = O) et une lecture simultanée de tous les détecteurs correspond à T = O, donc Q = O. Dans le cas du mode continu, Y e (gamma équivalent) remplace donc X et A e remplace h sauf en ce qui concerne le décalage de 20 l'ensemble de détecteurs par rapport à la droite OF, qui reste 4 ^/. D'autres modes de scrutation entraînent une variation du coefficient c Par exemple, si on commence la scrutation progressive en partant de l'autre extrémité de l'ensemble de détecteurs, 25 tout ce qui précède est valable avec o = v Z;. REVENDICATIONS 1 Procédé de reconstruction d'image par tomodensitométrie, du type consistant à faire tourner un ensemble source-détecteurs dans le plan d'une coupe à imager, ladite source émettant un faisceau de rayonnement pénétrant en forme d'éventail et à relever 5 et mémoriser une succession de vues pour différents angles O de rotation de l'ensemble sourcedétecteurs par rapport à un système d'axes de référence fixe (XOY) dudit plan, lesdits détecteurs étant agencés suivant un arc de cercle (CC') centré sur le foyer (F) de ladite source et chaque vue étant échantillonnée pour des angles ' 10 représentatifs des détecteurs et repérés à partir de la droite (OF) joignant le centre de rotation (O) de l'ensemble source-détecteurs et ledit foyer (F), le décalage angulaire entre deux détecteurs adjacents étant L constant et le décalage angulaire entre vues étant constant, caractérisé en ce qu'il consiste, le décalage angulaire 15 entre ladite droite (OF) et le détecteur le plus proche de celle-ci étant réglé à / / et les conditions suivantes étant réalisées: = N avec N entier positif en posant X e = ( 1 +)Y o O est un coefficient dépendant de la 20 vitesse de rotation de l'ensemble source-détecteurs et du temps d'acquisition d'une vue; à relever et mémoriser deux séries de vues (Vi, Wi) avec un décalage entre vues de h /3 pour la première série (Vi) et 4/ pour la seconde série (Wi), au moins certaines des vues des deux séries 25 étant décalées entre elles de / A 2 el avec A= 2/= e, h étant entier positif à compléter les vues de la première série par des valeurs échantillonnées de la seconde série dont les coordonnées (e', Y') repérées dans l'espace (,f avec 9 =/ + ', satisfont à la relation: (ex, "x) = ( 6 ' + ( 2 k + 1)I/, ') avec k entier relatif (C Ox, Y x) étant les coordonnées dans l'espace lô &J de valeurs échantillonnées recherchées de la première série; à appliquer un algorithme de reconstruction, connu en soi, 5 comprenant une convolution et une rétro-projection, aux vues complétées de la première série, pour reconstituer une image. 2 Procédé selon la revendication 1, caractérisé en ce qu'on relève et mémorise les vues (Vi) d'une série pendant un tour de l'ensemble sourcedétecteurs et les vues (Wi) de l'autre série pen10 dant le tour suivant et qu'on provoque un décalage de / 2 e/dans le relevé des deux séries, entre les deux tours. 3 Procédé selon la revendication 1, caractérisé en ce qu'on relève et mémorise les vues (Vi, Wi) des deux séries alternativement pendant un même tour. 4 Procédé selon la revendication 2 ou 3, caractérisé en ce que h étant égal à 1, N vues de la première série sont complétées par N vues de la seconde série. Procédé selon la revendication 3 ou 4, caractérisé en ce que h étant supérieur à 1, on ne relève et/ou traite qu'une valeur 20 échantillonnée sur h dans chaque vue de la seconde série. 6 Procédé selon l'une des revendications précédentes pour une tomodensitomètrie en mode pulsé, caractérisé en ce que Ye = 7 Procédé selon l'une des revendications I à 5 pour une tomodensitométrie en mode continu, caractérisé en ce que 25 =-/, é étant la vitesse de rotation de l'ensemble sourcedétecteurs et U étant l'intervalle de temps séparant l'acquisition des valeurs échantillonnées correspondant à deux détecteurs voisins. 8 Tomodensitomètre du type comportant un ensemble sourcedétecteurs (F, 13) monté tournant dans un plan de coupe à imager, 30 ladite source (F) étant une source de rayonnement penètrant en forme d'éventail ( 11) et lesdits détecteurs ( 14) étant agencés pour converger vers le foyer (F) de ladite source, lesdits détecteurs étant vus dudit foyer suivant des angles Y respectifs repérés à partir de la droite joignant ledit foyer et le centre de rotation (O) de l'ensemble source-détecteurs, deux détecteurs voisins étant angulairement décalés de 'Y, des moyens de pilotage ( 15) pour repérer un certain nombre d'angles /P de prises de vues de l'ensemble source-détecteurs par 5 rapport à un système d'axes de référence (XOY) et pour engendrer des signaux de commande desdites prises de vues (Vi, Wi), des moyens de mémorisation (M) d'informations représentatives desdites vues échantillonnées pour des angles P précités, des moyens de calcul ( 21) pour appliquer un algorithme de 10 reconstruction, connu en soi, comprenant une convolution et une rétro- projection à une série de vues, caractérisé en-ce que: les conditions suivantes sont satisfaites: qr = N avec N entier positif en posant Ye = ( 1 +o()' o c est un coefficient dépendant de la vitesse de rotation de l'ensemble source-détecteurs et du temps d'acquisition d'une vue, un décalage angulaire entre ladite droite (OF) et le détecteur ( 14) - le plus proche de celle-ci est prédéterminé à / les moyens de pilotage comportent des moyens ( 15) pour commander la prise et la mémorisation de deux séries (Vi, Wi) de vues échantillonnées, lesdites séries étant telles qu'au moins certaines vues voisines de séries différentes sont décalées entre elles de / Mel; les vues de la première série étant décalées angulairement au pas de h 13 et celles de la secondes série étant décalées angulairement au pas de /, avec ln = 2 Y et h entier positif, et en ce qu'il comporte des moyens de réarrangement (A, B, 22) pour compléter les vues de la première série par des valeurs échantillonnées prélevées dans la seconde série dont les coordonnées (O ', ') 30 dans l'espace),l avec O = J 3 +e satisfont à la relation (ex, x) = ( O O + ( 2 k + 1), È') avec k entier relatif, (O x, Y x) étant les coordonnées dans l'espace le 'ïel de valeurs échantillonnées recherchées de la première série. série. 9 Tomodensitomètre selon la revendication 8, caractérisé en ce que lesdits moyens de pilotage comportent un codeur incrémental ( 15). Tomodensitomètre selon la revendication 9, caractérisé en 5 ce que le pas élémentaire de comptage dudit codeur incrémental est -, ou une sous-multiple de cette valeur. 11 Tomodensitomètre selon l'une des revendications 8 à 10, du type fonctionnant en mode pulsé, caractérisé en ce que 'e = ô. 12 Tomodensitomètre selon l'une des revendications 8 à 10 du 10 type fonctionnant en mode continu, caractérisé en ce que = NZ, v étant la vitesse de rotation de l'ensemble source-détecteurs et T étant l'intervalle de temps séparant l'acquisition des valeurs échantillonnées correspondant à deux détecteurs voisins.