Brevet FR2897255A1 — publié le 20070817

PROCEDE D'ESTIMATION DU RAYONNEMENT DIFFUSE EN TOMOGRAPHIE PAR RAYONS X DESCRIPTION Le sujet de cette invention est un procédé 5 d'estimation du rayonnement diffusé en tomographie par rayons X. De nombreuses méthodes d'imagerie d'objets sont fondées sur l'atténuation subie par un rayonnement, et qui dépend de la nature de chaque 10 portion des objets. Les rayons se projettent sur un détecteur bidimensionnel placé derrière l'objet. L'image tridimensionnelle est obtenue en procédant à plusieurs irradiations de l'objet à des angles différents d'incidence afin de disposer d'une quantité 15 suffisante de mesures, et en recourant à un algorithme appelé d'inversion pour passer de la série d'images bidimensionnelles prises par le détecteur à l'image tridimensionnelle recherchée. Or chaque pixel du détecteur reçoit, outre 20 un rayon primaire projeté à travers l'objet, du rayonnement diffusé par l'objet et qui trouble la mesure. Le rayonnement diffusé, dont l'intensité est souvent importante, doit être estimé et corrigé pour améliorer la représentation de l'objet. 25 Plusieurs procédés ont été utilisés pour cela. Dans l'un d'entre eux, on recourt à un engin absorbeur du rayonnement composé d'un réseau de petites billes ("beam-stops") séparées de distances connues. Les rayons qui arrivent devant les billes sont entièrement absorbés, les autres passent sans atténuation. Quand on place à la fois l'objet à étudier et l'absorbeur entre la source du rayonnement et le détecteur bidimensionnel, les régions du détecteur qui sont situées derrière les billes forment des taches d'ombre illuminées seulement par le rayonnement diffusé par l'objet. On en déduit, par des interpolations entre les taches d'ombre, l'image du rayonnement diffusé sur toute la surface du détecteur. Cette image est soustraite de l'image ordinaire, prise en l'absence de l'absorbeur, de l'objet étudié. Le procédé donne de bons résultats mais nécessite au moins une autre irradiation de l'objet, qui reçoit donc une dose beaucoup plus importante, nuisible dans le cas courant d'examen d'êtres vivants. Une autre façon de procéder consiste à disposer deux détecteurs l'un derrière l'autre, dont le premier reçoit le rayonnement complet et le second le rayonnement primaire seul, grâce à des collimateurs qui arrêtent le rayonnement diffusé. Chaque acquisition d'une image bidimensionnelle par le détecteur est faite en une seule fois contrairement au procédé précédent, de sorte qu'aucune augmentation de dose d'irradiation n'est subie par l'objet. On estime le rayonnement primaire sur le second détecteur en le mesurant aux endroits qui sont alignés avec une ouverture de collimateur et la source du rayonnement, et en appliquant une correction nécessitée par l'absorption du rayonnement dans le premier détecteur. Une soustraction des mesures du premier détecteur par le rayonnement primaire estimé par le second détecteur aux mêmes endroits y donne une estimation du rayonnement diffusé ; des interpolations entre ces endroits complètent encore une image bidimensionnelle du rayonnement diffusé. L'inconvénient de cette façon de procéder est que l'installation est coûteuse à cause des deux détecteurs. Dans une autre catégorie de procédés, on utilise des calibrations pour estimer le rayonnement diffusé. Des images sont prises en utilisant, à la place de l'objet à étudier, des plaques ou d'autres pièces de forme simple, souvent appelées fantômes de calibration, construites de préférence en des matières ayant des propriétés d'absorption et de diffusion analogues à celle de l'objet. On utilise aussi le réseau de petites billes du premier procédé pour l'appliquer aux fantômes de calibration et obtenir successivement des images du rayonnement primaire et du rayonnement diffusé pour chaque fantôme de calibration. Les mesures faites avec l'objet sont corrigées par une fonctionnelle en utilisant un noyau de convolution défini par les résultats de calibration, en utilisant de préférence les images d'un fantôme choisi, supposé ressembler le plus à l'objet en raison d'une similitude de leurs mesures. Ce procédé est assez simple mais implique des approximations importantes. Un exemple est donné dans le document US 2005/0078787 Al. Dans une autre famille de procédé, on travaille sur des voxels de l'objet et non sur des projections du rayonnement. On peut citer une méthode de Monte Carlo par itérations où, après la reconstruction d'une première image tridimensionnelle de l'objet en utilisant les mesures brutes, le rayonnement diffusé correspondant est simulé, puis soustrait de chacune des mesures pour obtenir une première estimation des projections de rayonnement primaire. Une deuxième image tridimensionnelle de l'objet est reconstruite, et son rayonnement diffusé est de nouveau simulé et soustrait des mesures. Les images tridimensionnelles convergent vers l'image réelle de l'objet en répétant ce procédé, mais les temps de calcul sont excessifs. Dans un autre procédé de cette seconde famille, le rayonnement diffusé par chacun des voxels d'une première image tridimensionnelle de l'objet est évalué par une méthode analytique, puis corrigé par un facteur, et une deuxième image de l'objet est encore construite en soustrayant le rayonnement diffusé estimé des mesures. Contrairement au procédé précédent, aucune nouvelle image de l'objet n'est recherchée, le traitement analytique étant supposé donner une convergence immédiate ; mais le temps de calcul reste important pour des résultats moins bons que ceux du procédé précédent. Un exemple de ce procédé est donné par le document US 2005/0185753 Al. L'invention proposée ici implique une combinaison originale, comprenant à la fois l'application d'une fonctionnelle de convolution typique des méthodes d'analyse projection par projection, tout en recourant temporairement à une décomposition de l'objet en voxels, d'une façon qui permet à la fois d'estimer le rayonnement diffusé avec une meilleure précision que les méthodes utilisant les projections sans nécessiter les temps de calcul importants associés aux méthodes où l'objet est décomposé en voxels. Sous sa forme la plus générale, l'invention est relative à un procédé d'estimation d'un rayonnement diffusé par un objet étudié en utilisant un rayonnement traversant l'objet et un détecteur bidimensionnel recueillant le rayonnement après qu'il a traversé l'objet, comprenant des prises d'images de l'objet par le détecteur en mesurant un rayonnement total, somme du rayonnement diffusé et d'un rayonnement primaire ayant traversé l'objet sur des projections rectilignes ; et aussi des prises d'images de fantômes de calibration traversés successivement par le rayonnement en mesurant encore pour chaque fantôme un rayonnement diffusé par le fantôme et en rayonnement total, somme du rayonnement diffusé par le fantôme et d'un rayonnement primaire ayant traversé le fantôme sur des projections rectilignes ; consistant ensuite à choisir un des fantômes, et à extrapoler le rayonnement diffusé par l'objet du rayonnement diffusé par le fantôme choisi, caractérisé en ce qu'il comprend encore des étapes : de formation d'une image tridimensionnelle de l'objet à partir des images de l'objet prises par le détecteur ; d'attribution de coefficients respectifs de diffusion du rayonnement à chaque portion de l'image tridimensionnelle de détermination, pour les projections du rayonnement à travers l'objet, d'une image de longueurs équivalentes de matériau des fantômes quant à une quantité totale de rayonnement diffusé sur les projections ; et de choix du fantôme d'après une similitude de l'image de longueurs équivalente et d'une épaisseur traversée par le rayonnement de mesure dudit fantôme choisi. De façon favorable, l'extrapolation est faite en appliquant, à chacune des projections, des fonctionnelles au rayonnement total de l'objet et au rayonnement total du fantôme choisi, pour calculer une diffusion du rayonnement par un calcul ; et le procédé comprend encore des étapes consistant à faire un rapport numérique de la fonctionnelle sur l'objet et de la fonctionnelle sur le fantôme choisi, et à multiplier le rayonnement diffusé par le fantôme choisi par le rapport numérique, pour chacune des projections. Dans une réalisation préférée, les fonctionnelles sont de la forme F = (q)acquis x L) *K, où F est le rayonnement diffusé simulé, acquis est le rayonnement total, x est un opérateur de multiplication, L est l'image des longueurs équivalentes de l'objet ou la longueur équivalente du fantôme choisi traversée par la projection, * est un opérateur de convolution bidimensionnelle et K est un noyau de distribution bidimensionnelle du rayonnement diffusé par la projection. L'invention sera décrite maintenant en référence aux figures qui en illustrent certaines réalisations et que voici . - la figure 1 illustre l'acquisition de mesures sur l'objet à étudier, - la figure 2 illustre le mode de calibration, - la figure 3a et 3b donnent la représentation d'un objet effectivement étudié en vues de dessus et de côté, - la figure 4a, 4b, 4c et 4d représentent les résultats des reconstructions obtenues à partir de projections de mesures brutes et de projections corrigées en soustrayant le rayonnement diffusé par cet objet, estimé de trois façons différentes, - et la figure 5 est un organigramme du procédé. La figure 2 montre que l'appareillage comprend une source 1 et un détecteur bidimensionnel 2. Un objet 8 à étudier est placé entre eux et atténue un rayonnement 5 conique émis par la source 1 avant qu'il n'arrive au détecteur 2. Le rayonnement 5 traverse l'objet 8 par des projections rectilignes. Le détecteur 2 mesure l'atténuation du rayonnement sur chacune des projections. Une image tridimensionnelle de l'objet peut être obtenue en faisant tourner ensemble la source 1 et le détecteur 2 autour de lui et en prenant plusieurs images bidimensionnelles par le détecteur 2 à des incidences différentes, dont la synthèse est faite. Dans une étape préliminaire du procédé illustrée à la figure 1, correspondant à une calibration, l'objet 8 est cependant absent et remplacé dans la même installation par une série de fantômes 7 dont l'un d'entre eux est illustré et qui sont étudiés successivement après avoir été placés dans le rayonnement 5. Il s'agit de plaques ou d'autres structures simples construites en un matériau homogène, dont les propriétés de diffusion et d'absorption du rayonnement sont de préférence analogues à celles d'un des matériaux devant constituer l'objet 8 afin d'avoir une bonne similitude des comportements. Un absorbeur 3 est aussi placé dans le rayonnement 5 ; il est composé d'éléments absorbeurs tels que des billes 4 séparées de distances connues. Le rayonnement 5 se projette sur le détecteur 2 avec des taches d'ombre 6 correspondant aux projections des billes 4 et où on mesure, en l'absence de rayonnement primaire ayant traversé le fantôme 7 d'un trajet rectiligne et qui a été complètement absorbé par les billes 4, le rayonnement diffusé par le fantôme 7 à ces endroits (étape El de la figure 5). Une image complète du rayonnement diffusé est obtenue en interpolant les mesures faites aux taches d'ombre 6 sur le reste de l'étendue du détecteur 2 (étape E2). Un noyau de distribution du rayonnement diffusé peut être déduit d'autres mesures faites sur ce fantôme 7 (étape E3), éventuellement au moyen d'autres techniques connues dans l'art. Ce noyau de distribution exprime l'étalement du rayonnement diffusé le long d'une projection sur la région environnante de l'image. On peut souvent le modéliser par une fonction gaussienne bidimensionnelle en fonction des rayons séparant le point d'arrivée de la projection et chacun des points de la régions environnante, et qui dépend alors d'un seul paramètre : l'écart-type. Une autre mesure du fantôme 7, l'absorbeur 3 étant cette fois enlevé, permet d'en obtenir une autre image, qui exprime la somme du rayonnement primaire et du rayonnement diffusé à chaque pixel du détecteur 2 (étape E4). La soustraction des deux images à chaque pixel du détecteur 2 donne ensuite une image du rayonnement primaire seul (E5). Les images du rayonnement primaire, du rayonnement diffusé et le noyau de distribution restent associés au fantôme 7. Cela est répété pour tous les fantômes 7 de façon à obtenir un catalogue de calibration. Les étapes suivantes du procédé concernent l'objet 8. Une série d'images bidimensionnelles en est prise dans la configuration de la figure 1 à plusieurs incidences (étape E6). Une première reconstitution tomographique, tridimensionnelle de l'objet 8 en voxels est obtenue en utilisant un procédé d'inversion quelconque de l'art comme l'algorithme de Feldkamp (étape E7). Chaque voxel est exprimé par le coefficient d'atténuation qu'il est estimé faire subir au rayonnement, sans tenir compte de ce que le rayonnement diffusé a troublé les mesures. Les matériaux constitutifs de l'objet 8, supposés connus au préalable, sont respectivement associés aux voxels en appliquant des critères de seuil d'après les valeurs connues des coefficients d'atténuation de ces matériaux constitutifs et ceux que le procédé d'inversion a calculé aux voxels (étape E8). Des coefficients de pondération "p", égaux aux rapports de la section efficace de diffusion Compton des matériaux sur celle du matériau de calibrage (des fantômes 7), sont attribués aux voxels d'après le matériau qui a été associé à chacun (étape E9). Pour des applications en anatomie, on pourra retenir l'air, les tissus mous et les os comme matériaux constitutifs. Une image bidimensionnelle d'une longueur traversée équivalente, Lequ, est ensuite obtenue à l'aide de l'image tridimensionnelle des coefficients de pondération "p" en appliquant la formule : L(i) Lequ(i)= Jp(1)xdl, où L(i) est l'épaisseur 1=o réelle de l'objet traversée par le rayon primaire qui atteint le détecteur au pixel noté i (étape E10). Lequ correspond à la longueur que chaque rayon du rayonnement 5 devrait traverser dans le matériau des fantômes 7 pour donner le même rayonnement diffusé qu'à travers son parcours dans l'objet 8. Quand Lequ de chaque projection est obtenue, le fantôme 7 ayant l'épaisseur la plus proche de la moyenne de Lequ sur une zone d'indexation, est indexé à la projection (étape E11). En variante, un fantôme fictif ayant une épaisseur égale à Lequ et dont les propriétés seraient obtenues par interpolation entre deux fantômes 7 réels, mesurés, pourrait être indexé à la projection. Une fonctionnelle F, qui donne le taux de diffusion Compton du premier ordre du rayonnement, est ensuite appliquée à chaque projection mesurée du fantôme 7 indexé (étape E12) et de l'objet 8 (étape E13). Elle est définie comme : F = 4 acquis x Lequ K pour l'objet 8, où L représente l'image des longueurs équivalentes calculées sur l'objet, acquis le rayonnement total parvenu au détecteur 2 pour la projection considérée, et K le noyau de distribution du rayonnement diffusé. La même fonctionnelle est appliquée au fantôme 7 indexé en remplaçant Lequ par l'épaisseur de celui-ci. L'application de la fonctionnelle donne une simulation du pouvoir de diffusion du rayonnement. Le rapport de la fonctionnelle F de l'objet 8 sur celle du fantôme 7 est appelé nappe et noté a (obtenu à l'étape E14). Il est défini pour chacune des projections et donne le rayonnement diffusé relatif de l'objet 8 par rapport à celui du fantôme 7. La nappe a est ensuite multipliée par le diffusé du fantôme 7, tel qu'il a été obtenu au moyen de l'absorbeur 3 et des interpolations des valeurs mesurées aux taches d'ombre 6, pour obtenir une estimation du rayonnement diffusé de l'objet 8 sur chacune des projections (étape E15). L'image du diffusé D de l'objet 8 est ainsi estimée par la formule suivante, appliquée à chacune des projections : / D =Dcalibx /~ acquis X Lequ *K OU Dcalib est l'image el) \ calib X L calib du rayonnement diffusé associée au fantôme 7 d'épaisseur équivalente à l'objet 8, acquis et (1)calib sont les projections du rayonnement total reçu pour l'objet 8 et le fantôme 7, Lequ et Lcalib sont l'image de longueurs équivalentes de l'objet 8 traversé par le rayonnement 5 et l'épaisseur du fantôme 7, et K est le noyau de distribution du rayonnement calibré. Des résultats d'essais sont donnés au moyen des figures suivantes. Les figures 3a et 3b représentent un objet composé d'un cylindre 9 de Plexiglas et muni de deux perçages 10 coaxiaux d'air et pouvant donc être comparé, pour l'absorption du rayonnement, à un thorax et à ses deux poumons. Les figures 4a, 4b, 4c et 4d représentent des résultats d'atténuation du rayonnement à travers le cylindre 9, de reconstructions de coefficients à partir de projections mesurées dans toutes les directions par un rayonnement 5 occupant nécessairement plusieurs positions d'incidence autour du cylindre 9. La reconstruction était faite sur une ligne 17 diamétrale du cylindre 9 donnée sur la courbe 11 commune à toutes les figures 4a et 4b passant entre les perçages 10 et donc se traduit par un plateau 12 correspondant à un matériau homogène. L'objet est homogène donc le profil coefficients d'atténuation attendu est un plateau par un rayonnement 5. La figure 4a représente les coefficients d'atténuations bruts 13, obtenus à partir des projections mesurées, sans corriger l'influence du rayonnement diffusé : le plateau 12 est remplacé par une forme d'auge. A la figure 4b, on donne les résultats obtenus en faisant une seconde mesure en utilisant l'absorbeur 3 pour évaluer le rayonnement diffusé puis le soustraire du rayonnement total (courbe 14) : on voit que le rayonnement diffusé est bien évalué, les coefficients d'atténuation estimés étant partout presque semblables aux coefficients d'atténuation réels, mais au prix d'une irradiation accrue ainsi qu'on l'a signalé. A la figure 4c, on donne le résultat obtenu selon le procédé, commenté plus haut, de US 2005/0078787 Al (courbe 15) : le résultat est meilleur que celui de la figure 4a, mais il n'est pas tout à fait satisfaisant, le rayonnement diffusé étant en général surestimé et des variations des coefficients d'atténuation reconstruits étant présentes dans cet exemple. La figure 4d illustre enfin l'utilisation de l'invention (courbe 16) : les résultats sont aussi bons que dans le cas de la figure 4b, et sans augmentation de la dose reçue. REVENDICATIONS 1) Procédé d'estimation d'un rayonnement diffusé par un objet (8) étudié en utilisant un rayonnement (5) traversant l'objet et un détecteur (2) bidimensionnel recueillant le rayonnement après qu'il a traversé l'objet, comprenant : des prises d'images de l'objet par le détecteur en mesurant un rayonnement total, somme du rayonnement diffusé et d'un rayonnement primaire ayant traversé l'objet sur des projections rectilignes ; et aussi des prises d'images de fantômes de calibration (7) traversés successivement par le rayonnement (5) en mesurant encore pour chaque fantôme un rayonnement diffusé par le fantôme et un rayonnement total, somme du rayonnement diffusé par le fantôme et d'un rayonnement primaire ayant traversé le fantôme sur des projections rectilignes ; consistant ensuite à choisir un des fantômes ; et à extrapoler le rayonnement diffusé par l'objet du rayonnement diffusé par le fantôme choisi ; caractérisé en ce qu'il comprend encore des étapes : de formation d'une image tridimensionnelle de l'objet à partir des images de l'objet prises par le détecteur ; d'attribution de coefficients respectifs de diffusion du rayonnement à chaque portion de l'image tridimensionnelle ; de détermination, pour les projections du rayonnement à travers l'objet, d'une image de longueurs équivalentes de matériau des fantômes quant à une quantité totale de rayonnement diffusé sur les projections ; et de choix du fantôme d'après une similitude de la longueur équivalente et d'une épaisseur traversée par le rayonnement de mesure (5) dudit fantôme choisi. 2) Procédé d'estimation d'un rayonnement diffusé par un objet selon la revendication 1, caractérisé en ce qu'un des fantômes est choisi indépendamment pour chacune des projections. 3) Procédé d'estimation d'un rayonnement diffusé par un objet selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que l'extrapolation est faite en appliquant à chacune des projection des fonctionnelles au rayonnement total de l'objet et au rayonnement total du fantôme choisi, pour calculer une diffusion du rayonnement par une simulation, et en ce que le procédé comprend encore des étapes consistant à faire un rapport numérique de la fonctionnelle sur l'objet et de la fonctionnelle sur le fantôme choisi, et à multiplier le rayonnement diffusé par le fantôme choisi par le rapport numérique, pour chacune des projections. 4) Procédé d'estimation d'un rayonnement diffusé par un objet selon l'une quelconque des revendications 1 ou 2, où l'objet est un être vivant, caractérisé en ce que les coefficients de diffusion des portions de l'image tridimensionnelle sont déterminés en attribuant des compositions d'air, de tissus mous ou d'os auxdites portions. 5) Procédé d'estimation d'un rayonnement diffusé par un objet selon la revendication 3, caractérisé en ce que les fonctionnelles sont de la forme F = (acquis x L) *K, où F est le rayonnement diffusé simulé, acquis est le rayonnement total, x est un opérateur de multiplication, L est l'image de longueurs équivalentes de l'objet ou la longueur du 5fantôme choisi traversée par la projection, * est un opérateur de convolution bidimensionnelle et K est un noyau de distribution bidimensionnelle du rayonnement diffusé par la projection. 10