FR 2670038 A1 19920605 FR 9014957 A 19901129 PROCEDE ET DISPOSITIF DE RECONSTRUCTION D'IMAGES TRIDIMENSIONNELLES D'UN OBJET EN UTILISANT DEUX TRAJECTOIRES CIRCULAIRES D'AXE COMMUN DESCRIPTION L'invention se rapporte à un procédé et un dispositif de reconstruction d'images tridimen- sionnel Les d'un objet qui comporte des mesures sur deux trajectoires circulaires d'axe commun, c'est-à-dire que les trajectoires appartiennent à des plans paral Lèles et Leurs centres sont joints par un axe perpendiculaire au plan des trajectoires. Si les plans des trajectoires sont confondus, elles sont concentriques On peut aussi envisager des trajec- toires plus nombreuses. La reconstruction des images concerne un objet déterminé à examiner, et les moyens matériels employés comportent deux réseaux bidimensionnels de capteurs parcourant chacun une des trajectoires ou, de manière équivalente, un seul réseau qui parcourt toutes les trajectoires. Limage de l'objet est définie par des valeurs prises par une fonction sur chacun de ses points. La fonction est d'autre part une propriété d'un rayon- nement (X ou gamma entre autres) de forme conique ayant un point focal et qui passe à travers l'objet; chaque rayon est reçu par un des capteurs du réseau bidimensionnel et représente donc la somme de la fonc- tion sur tous les points de l'objet appartenant à ce rayon Un traitement approprié des sommes sur tous les rayons, et ceci pour un nombre suffisant de mesures suivant des incidences différentes autour de l'objet, permet de reconstituer l'image de l'objet. On se contente en pratique de considérer des nombres finis de rayons et de points de l'objet selon des discrétisations ou des maillages. Cette invention constitue un perfectionnement d'une invention précédente, décrite dans la demande de brevet européen EP-A-0 292 402, mais on peut envisa- ger d'employer l'invention actuelle dans d'autres circonstances ou en utilisant d'autres méthodes mathé- matiques de traitement des données. Les méthodes envisageables utilisent en particulier ce qu'on appelle la transformée de Radon de la fonction à mesurer, qui est définie comme la somme de la fonction sur chacun des plans, appelés plans de Radon, qui passent à travers l'objet considéré ou, mieux encore, la dérivée première de cette trans- formée La contribution des points de ces plans qui n'appartiennent pas à l'objet est considérée comme inexistante, ce qui est valable dans le cas d'un rayon- nement qui traverse un milieu gazeux ambiant sans être atténué Ici encore, une discrétisation est effec- tuée pour n'effectuer les calculs que sur un nombre fini de plans. La dérivée première de la transformée de Radon est définie comme la dérivée de la transformée de Radon en fonction de la variable p définissant la distance du plan considéré à une origine Pour chaque plan, elle correspond à la somme sur ce plan de la dérivée première de la transformée de Radon dans la direction perpendiculaire au plan L'invention décrite dans le brevet EP-A-0 292 402 montre que pour un plan passant par une position du point focal du cône et rencontrant le réseau bidimensionnel de cap- teurs, il est possible de calculer à partir des mesures La valeur exacte de la dérivée première de la trans- formée de Radon sur ce plan Ce calcul, suivant les formules exposées dans le brevet EP-A-0 292 402, met en oeuvre de façon préférentielle une pondération de correction d'éloignement du point focal, deux fil- trages correspondant aux calculs des dérivées premières respectivement suivant Les lignes et Les colonnes du réseau de capteurs, deux sommations suivant la droite d'intersection entre le plan à traiter et le réseau bidimensionnel de capteurs, puis une combinaison Linéaire et une normalisation des résultats Cette sommation met en oeuvre les interpolations nécessaires car ces droites d'intersection passent entre les ran- gées de capteurs ou les coupent. Dans le reste de la demande, on entend par somme la somme pondérée des valeurs mesurées (dans le cas de la transformée de Radon) aussi bien que la combinaison linéaire des sommes des valeurs pondé- rées et filtrées (dans le cas de la dérivée première de la transformée de Radon) obtenues selon les lignes et les colonnes du réseau. La somme de la fonction sur les points des plans de Radon est facile à obtenir à condition que ces plans aient une intersection avec le réseau bi- dimensionnel de capteurs et qu'ils passent par le point focal unique que visent les capteurs Il suffit d'effectuer la somme des valeurs mesurées par chaque capteur situé à l'intersection, avec les interpolations nécessaires car les intersections des plans de Radon passent entre des rangées de capteurs ou les coupent. Une fois que les valeurs de la fonction sur les plans de Radon ont été calculées, il existe des formules d'inversion, exposées dans la demande de brevet men- tionnée ci-dessus, qui permettent de parvenir aux valeurs de la fonction sur Les points du maillage de l'objet. Il faut toutefois revenir sur les conditions qui permettent d'obtenir un nombre de plans de Radon suffisant pour permettre une description satisfaisante de l'objet Chaque plan de Radon peut être défini par ce qu'on appelle son point caractéristique, c'est- à-dire Le point de projection sur ce plan d'une origine choisie arbitrairement Ce point caractéristique, noté C sur la figure 1, peut être défini par ses coordonnées sphériques P, X et Ode rayon, longitude et colatitude respectivement à partir de L'origine 0 Le plan de Radon P passant par le point caractéris- tique C peut être défini par le rayon p et Le vecteur unitaire -nde direction OC. Les valeurs de la fonction sur les plans de Radon ne peuvent effectivement être calculées que pour les plans de Radon qui coupent la trajectoire parcourue par Le point focal du rayonnement Ce point focal est concrètement, dans Le cas d'une fonction d'atténuation, une source ponctuelle de rayons X, gamma, etc C'est la conception qui est décrite dans la demande de brevet européenne mentionnée Les mêmes conditions géométriques existent en tomographie d'émis- sion, quand la fonction à mesurer est l'activité émise par Le corps Le point focal n'a alors pas d'existence physique et correspond simplement au point de conver- gence vers Lequel tous les collimateurs utilisés devant le réseau bidimensionnel sont dirigés. IL faut, pour obtenir un nombre suffisant de plans de Radon, effectuer plusieurs mesures avec des positions différentes du point focal. Considérons une trajectoire circulaire T parcourue par le point focal (ou la source) S, le réseau plan de capteurs Pdét parcourant La même trajec- toire que le point focal S ou éventuellement une tra- jectoire concentrique de rayon différent Supposons encore que L'origine O qui sert à définir les points caractéristiques C soit confondue avec le centre de La trajectoire T Le volume enveloppant les points caractéristiques correspond à un tore To produit par la rotation d'une surface sphérique de diamètre OS autour de La trajectoire T En effet, L'axe de rotation des plans de Radon qui passent par Le point focal S ont leurs points caractéristiques répartis sur la surface sphérique de diamètre OS car L'angle SCO est droit Le tore To est appelé le volume caractéristique des mesures, qui dépend donc de la forme de la trajec- toire T et de sa position par rapport à l'origine O. Il est d'autre part suffisant, pour obtenir une description complète de l'objet par Les procédés utilisant Les sommes de la fonction sur les plans de Radon, de pouvoir disposer des points caractéris- tiques appartenant à un volume caractéristique de L'objet Ce volume caractéristique de l'objet est toujours inclus dans une sphère V centrée sur l'origine O et qui englobe l'objet M On peut donc assurer que si Le volume caractéristique de l'objet est inclus dans Le volume caractéristique des mesures, la recons- truction de l'image de L'objet M sera possible Dans le cas d'une trajectoire circulaire T, cette condition n'est malheureusement pas remp Lie car une zone d'ombre subsiste pour les points caractéristiques de La sphère V qui n'appartiennent pas au tore To et dont les plans de Radon ne coupent pas la trajectoire T. L'objet essentiel de l'invention est de faciliter les reconstructions en les rendant plus rapides si plusieurs chaînes de mesure sont employées, ce qui permet de réduire la durée des examens, en imagerie médicale o le patient doit rester immobi Le 3 u et en contrôle non destructif o on est astreint à des cadences de travail; ou encore, Le procédé permet de traiter de-s cas o chaque chaîne d'occupation est définie par des paramètres propres et en particulier Le cas o les deux centres des trajectoires sont proches mais distincts pour des raisons mécaniques. C'est alors La qualité des images qui est améliorée. Enfin, quand les deux trajectoires circulaires sont de part et d'autre de l'objet, La zone d'ombre peut être réduite et l'image améliorée ici encore La super- position de plusieurs trajectoires permet aussi d'aug- menter le volume caractéristique de mesures et autorise l'examen de plus gros objets Mais on se heurte dans tous les cas à des difficultés de calcul car le regrou- pement des acquisitions faites sur plusieurs trajec- toires nécessite des conversions mettant en oeuvre des opérations mathématiques à cause du paramétrage de chaque série de mesures par des coefficients carac- téristiques propres. De plus, si on utilise des chaînes de mesures différentes et si les sommes de la fonction sont d'abord produites dans des repères séparés, les imprécisions mécaniques sur l'orientation réelle des capteurs, les rayons réels de leur trajectoire et du point focal, etc, se traduisent par des calculs et des interpolations supplémentaires quand les changements de repère sont effectués pour regrouper les résultats des séries de mesures Les informations de reconstruction des images sont donc dégradées. Si par contre on utilise le procédé conforme à l'invention, les paramètres des trajectoires et de séries de mesures, qui intègrent les imprécisions mécaniques et ont été mesurés au cours d'étalonnages préliminaires de la machine, sont utilisés directement dans le repère utilisé pour reconstruire l'objet, dont les images seront alors reconstruites avec une netteté meilleure. L'invention repose en conséquence sur la constatation que, si on définit un repère commun pour paramétrer les trajectoires et les volumes caractéris- tiques de mesure, aucune conversion de coordonnées ni de changement de repère n'est nécessaire pour inté- grer Les résultats provenant de trajectoires diffé- rentes On obtient directement les informations inté- ressantes dans Le repère d'objet dans lequel Le maillage utilisé pour reconstruire les images est défini Or, une telle façon de procéder est facile si les trajectoires sont coaxiales. L'invention est, pour récapituler, relative à un procédé de reconstruction d'images tridimension- nelles d'un objet défini par des valeurs prises par une fonction sur des points de L'objet, la fonction étant une propriété d'un rayonnement conique ayant un point focal et passant à travers l'objet, dans lequel la fonction est calculée par l'intermédiaire de sommes de la fonction sur des plans passant par au moins un point de l'objet et définis dans un repère d'objet, comportant au moins deux séries de mesures, chaque somme de la fonction sur lesdits plans para- métrés dans le repère d'objet étant calculée à partir d'au moins une des séries de mesures, chacune des séries de mesures étant menée avec un réseau bidimen- sionnel de capteurs de rayonnement, orientés vers le point focal, et qui est déplacé autour de l'objet sur une trajectoire circulaire respective, caractérisé en ce que l'objet et les trajectoires ont des centres joints par un axe perpendiculaire aux trajectoires, sont définis par un même repère, en ce que, pour chaque série de mesures, une série de sommes de la fonction est calculée dans ce repère pour les plans sécants ou tangents à une trajectoire associée du point focal, et en ce que les séries de sommes de la fonction sont regroupées pour obtenir les sommes de la fonction dans ledit repère. Avantageusement, le repère d'objet a un axe confondu avec l'axe perpendiculaire aux trajec- toires. Les coordonnées utilisées sont avantageuse- ment sphériques. L'invention est aussi relative aux dispo- sitifs conçus spécialement pour la mise en oeuvre du procédé. Plusieurs variantes de l'invention sont possibles en utilisant des dispositifs légèrement différents et ces procédés vont être décrits en réfé- rence aux figures suivantes annexées à titre illus- tratif et non limitatif: la figure 1, déjà décrite, indique les conditions géométriques liées à des acquisitions à l'aide de trajectoires circulaires; la figure 2 représente une généralisation de ces conditions géométriques; la figure 3 représente un mode d'acquisi- tion de mesures; la figure 4 représente un dispositif appro- prié à l'acquisition des mesures selon la figure 3; et les figures 5 à 7 représentent trois dispo- sitifs plus complexes envisageables en accord avec L'invention. Examinons maintenant le cas de la figure 2, o la trajectoire circulaire référencée par Tl et de centre 01 appartient à un plan qui passe à dis- tance de l'origine O, ce plan étant perpendiculaire à l'axe Z Le volume caractéristique des mesures, c'est-à-dire l'ensemble des points caractéristiques des plans de Radon définis à partir de l'origine 0, est à nouveau un tore Toi engendré en faisant tourner une surface sphérique de diamètre O 51 le long de la trajectoire Tl et autour de l'origine O Le tore Toi a une forme complexe puisqu'il présente une partie d'intersection I au centre, car les coupes méridiennes du tore To 1 forment deux cercles qui sont partiellement confondus Un raisonnement géométrique simple permet de se convaincre que les plans associés aux points situés à l'intérieur de la zone d'intersection I ne coupent pas la trajectoire et que cette partie appar- tient à la zone d'ombre: le volume caractéristique des mesures représenté suivant cette coupe méridienne est celui qui est hachuré. On a représenté en pointillés la superposi- tion d'un autre tore To 2 résultant de la rotation d'un cercle de diamètre 052 autour d'une autre trajec- toire T 2 de centre 52 distinct de l'origine O, qui est ici représentée sensiblement au centre du segment 0102 Le volume caractéristique de mesures correspon- dant est également hachuré On s'aperçoit que bien que ces volumes caractéristiques ne coïncident pas et que la multiplication des trajectoires ayant un axe commun permette de décrire des objets plus volumi- neux, il reste toujours une zone d'ombre autour de l'axe Z Le procédé de l'invention n'est donc pas approprié pour décrire parfaitement les objets examinés à cet endroit même si le volume de la zone d'ombre est réduit En pratique, on effectuera des interpola- tions à l'intérieur de la zone d'ombre conformément aux indications de l'invention antérieure. Un mode possible pour la réalisation du procédé a été représenté sur la figure 3 Les deux trajectoires sont confondues en une trajectoire unique T sur laquelle on effectue N mesures en des points régulièrement répartis Q 1 à Qn définis par des angles i =I(i-1)2 à partir de l'axe X On mesure pour n chaque point Qi les valeurs de la fonction sur les plans de Radon dont les points caractéristiques appar- tiennent à la sphère de diamètre O Qi Dans le cas général, la somme de la fonction sur ces plans peut 1 O être calculée à partir de deux positions du point focal car ces p Lans coupent La trajectoire T en deux points, Ql et Q 5 dans l'exemple considéré En pra- tique, les points W retenus sont cependant des points qui sont définis régulièrement sur un maillage défini dans le repère 0,X,Y,Z par ses coordonnées sphériques régu Lières (rayon P, Longitude 4 et colatitude 0, si bien que le plan de Radon ne coupe pas La trajec- toire exactement aux points Q et que des approximations ou des interpolations sur les angles A/i sont donc nécessaires, suivant la méthode envisagée dans L'inven- tion précédente par exemple. Selon La présente invention, on dispose de deux chaînes de mesures dont chacune effectue Les mesures simultanément sur un sous-ensemble des points Q Les deux façons les plus normales de procéder consistent soit à attribuer les points d'ordre impair à l'un et Les points d'ordre pair à L'autre, soit à partager La trajectoire en deux moitiés égales dont Les points Q respectifs sont affectés à L'une des chaînes de mesures, chaque sous-ensemble comprenant alors une moitié des points Q' IL est encore possible que les sous-ensembles soient plus grands qu'une moi- tié, ce qui rendra possibles des recouvrements Une autre idée consiste à doubler les mesures aux points Q pour améliorer La qualité de L'image par une moyenne statistique. Pour chacun des points du volume caractéris- tique de l'objet, Le regroupement des séries de mesures consiste par exemple soit à considérer alternativement les sommes de fonction associées à chaque série de mesures pour les plans rencontrant Les deux trajec- toires des points focaux, soit à faire Leurs moyennes, soit à ne considérer qu'une somme de fonction associée quand Le plan ne rencontre qu'une des trajectoires. 1 1 Un dispositif approprié est matérialisé sur la figure 4 Il s'agit d'une installation de tomo- graphie médicale o le patient est couché sur une table 1 horizontale Un rail circulaire 2 s'étend dans un plan vertical et transversal et entoure la table 1 Il s'agit en fait d'une crémaillère qui porte deux chariots mobiles qui engrènent sur elle. Ces chariots ont été schématisés et on peut se reporter le cas échéant à l'invention antérieure pour plus de détails techniques Chacun d'eux porte une source de rayons X à proximité d'un écran portant un réseau bidimensionnel de capteurs Les capteurs de l'un des écrans 4 sont focalisés sur la source 3 de l'autre chariot mobile, et ceux de l'autre écran 4 ' sont foca- lisés sur la source 3 ' du premier chariot mobile. Les sources 3 et 3 ' émettent un faisceau conique 5 ou 5 ' L'objet à examiner est placé dans le cône des deux faisceaux 5 et 5 ' En fonctionnement, les deux chariots mobiles se déplacent simultanément en rotation égale de manière à rester diamétralement opposés. Les mouvements sont réglés par l'ordinateur de commande de l'installation Dans le cas d'une installation analogue en tomographie d'émission, les sources 3 et 3 ' seraient omises et le point focal des capteurs de chaque écran 4 ou 4 ' serait avantageusement situé au centre de l'autre écran 4 ou 4 ', pour simplifier les calculs d'acquisition Il serait encore possible de séparer les sources 3, 3 ' des écrans 4, 4 ' de sorte que les faisceaux 5 et 5 ' se coupent à angle droit par exemple. Une conception un peu différente est repré- sentée à la figure 5 Les chariots mobiles sont iden- tiques aux précédents mais placés cette fois-ci sur un rail respectif 2 ou 2 ' Les rails 2 et 2 ' sont coplanaires et concentriques Rien n'est modifié à l'exception des coefficients des relations géométriques de grossissement en fonction des distances des sources 3 et 3 ' et des écrans 4 et 4 ' à l'origine O. L'invention peut être mise en oeuvre (figure 6) avec un rail unique 2 qui porte une source unique 3 diamétralement opposée à un écran 4 portant un réseau bidimensionnel de capteurs Dans cette conception, on fait effectuer deux tours complets à la source 3 et l'écran 4, entre lesquels la table 1 subit une translation perpendiculaire au plan du rail 2 On utilise une telle solution pour accroître le volume caractéristique des mesures conformément à l'illustra- tion de la figure 2. Le même résultat peut être obtenu en dédou- blant le dispositif (figure 7): deux rails circulaires 2 et 2 ' portent chacun une source 3 ou 3 ' et un écran 4 ou 4 ' Les déplacements sur les deux rails 2 et 2 ' sont simultanés, ce qui accélère l'acquisition des mesures. Si les résultats pour un plan de Radon déter- miné sont obtenus à chaque trajectoire, une moyenne ou une imbrication sont effectuées avant de procéder aux opérations de reconstruction de l'image, telles que les opérations de filtrage, de pondération et de rétroprojection décrites dans le brevet antérieur. L'invention peut être mise en oeuvre d'autres manières C'est ainsi que le débattement angulaire sur chaque trajectoire circulaire peut être réduit à une partie de circonférence Un cas avantageux est celui d'une circonférence complète formée par l'union des débattements angulaires. On rappelle en particulier que des trajec- toires plus nombreuses sont possibles sans sortir du cadre de l'invention. Par ailleurs, on a considéré avant tout l'utilisation de ce procédé pour l'inversion de la transformée de Radon proprement dite ou de sa dérivée première On peut encore procéder sensiblement de la même façon avec une transformée de Fourier de tous Les rayons de L'espace de Radon qui ne coupent pas la zone d'ombre, et ceci pour chacun des vo Lumes caractéristiques de mesures Le regroupement des mesures ainsi transformées peut être mené dans ce repère commun. Au lieu de La transformée de Radon ou de sa dérivée première, on peut encore procéder de La même façon avec La transformée de Hi Lbert de ces trans- formées. L'invention peut servir pour Les applications usuelles d'imagerie médicale ou de contrôLe non des- tructif de pièces. REVENDICATIONS 1 Procédé de reconstruction d'images tri- dimensionnel Les d'un objet (M) défini par des valeurs prises par une fonction sur des points de L'objet, La fonction étant une propriété d'un rayonnement conique ( 5, 5 ') ayant un point focal et passant à travers l'objet (M), dans lequel la fonction est calcu- lée par l'intermédiaire de sommes de la fonction sur des plans (P) passant par au moins un point de l'objet et définis dans un repère d'objet, comportant au moins deux séries de mesures, chaque somme de la fonction sur lesdits plans paramétrés dans le repère d'objet étant calculée à partir d'au moins une des séries de mesures, chacune des séries de mesures étant menée avec un réseau bidimensionnel de capteurs de rayonne- ment, orientés vers le point focal, et qui est déplacé autour de l'objet sur une trajectoire circulaire res- pective, caractérisé en ce que l'objet et les trajec- toires ont des centres ( 01, 02) joints par un axe perpendiculaire aux trajectoires, sont définis par un même repère, en ce que, pour chaque série de mesures, une série de sommes de la fonction est calcu- lée dans ce repère pour les plans sécants ou tangents à une trajectoire associée du point focal, et en ce que les séries de sommes de La fonction sont regroupées pour obtenir Les sommes de la fonction dans ledit repère. 2 Procédé de reconstruction d'images tri- dimensionnelles suivant la revendication 1, caractérisé en ce que le repère d'objet a un axe (Z) confondu avec l'axe perpendiculaire aux trajectoires. 3 Procédé de reconstruction d'images tri- dimensionnelles suivant l'une quelconque des revendica- tions 1 ou 2, caractérisé en ce que les coordonnées sont des coordonnées sphériques. 4 Procédé de reconstruction d'images tri- dimensionnelles suivant l'une quelconque des revendica- tions 1 à 3, caractérisé en ce que le rayonnement est un rayonnement d'atténuation émis par une source ponctuelle au point focal (S), tournant autour de l'objet suivant une trajectoire circulaire de source. Procédé de reconstruction d'images tri- dimensionnelles suivant l'une quelconque des revendica- tions 1 à 3, caractérisé en ce que le rayonnement est un rayonnement d'émission provenant de l'objet, le point focal étant un point vers lequel convergent les directions de mesure des capteurs. 6 Dispositif de reconstruction d'images tridimensionnelles suivant l'une quelconque des reven- dications 1 à 5, caractérisé en ce qu'il comprend au moins un réseau bidimensionnel de capteurs, les réseaux étant déplacés autour d'un objet à examiner suivant au moins deux trajectoires circulaires et reliés à des moyens d'acquisition de mesures et de traitement des mesures appropriés à la mise en oeuvre d'un procédé conforme à l'une quelconque des revendica- tions précédentes.