Brevet FR2820966A1 — publié le 20020823

B 13768.3 FOCI PROCEDE DE RADIOGRAPHIE A DOUBLE ENERGIE, ET DISPOSITIF DE CALIBRATION POUR CE PROCEDE DESCRIPTION Le sujet de cette invention est un procédé de radiographie à double énergie, ainsi qu'un dispositif de calibration approprié à ce procédé. La radiographie à double énergie consiste à exposer un objet ou un patient à un rayonnement de deux énergies différentes, pour lesquelles les propriétés d'atténuation des matériaux constituant l'objet ou le patient ne varient pas de la même facon. A l' aide dun modèle de leurs fonctions d'atténuation, on obtient alors des informations sur la densité et la nature des matériaux traversés. En particulier, dans le cas de l'ostéodensitométrie, la masse d'os traversée peut être calculée en la distinguant de la contribution des tissus mous à l'atténuation du rayonnement. Le fondement théorique de la méthode sera rapidement résumé ci-dessous. Le flux 0 traversant une longueur l d'un matériau de coefficient datténnation lincique,u à partir d'un flux initial de rayonnement 00 est un 0 = 00e1. On appelle m la mesure datténuation, égale à lnt J o Io et I sont les signaux mesurés par un même détecteur sous les flux 00 et 0. Dans le cas d'un objet complexe composé d'un grand nombre de matériaux (d'indice i), chacun contribue à B 13769.3 JCI latténnation d'après sa longueur Li que les rayons traversent. Or chaque matériau peut être exprimé, pour la propriété d'atténuation, comme une combinaison linéaire de deux matériaux de base, d'après la formule: u = kl'Ul + k2lU2, o k1 et k2 sont des coefficients constants, et p1 et u2 les atténuations de ces matériaux de base, et on pourra exprimer les longueurs équivalentes A1 et A2 des matériaux de base traversées par le rayonnement par: tA,= Lj.k, l Au moyen de ces longueurs équivalentes, les matériaux de base pourront représenter un objet même s'il est en réalité de composition beaucoup plus complexe. Si les matériaux de base ont des propriétés d'atténuation analogues à deux matériaux prépondérants de l' objet, la représentation sera réaliste. Dans le cas d'examen sur des êtres vivants, les matériaux de base classiques sont le plexiglas (polyméthacrylate) pour simuler les tissus mous et l'hydroxyapatite pour simuler les tissus osseux. Le système d'équations reliant les mesures aux atténuations,Uj et aux longueurs équivalentes Aj des matériaux de base est linéaire et donc facile à résoudre si le rayonnement est monochromatique; mais il ne l'est pas dans les situations réelles, et les B 13769.3 JCI longueurs équivalentes sont alors données par des modèles mathématiques plus compliqués, tels que celui-ci: A=aO+a,.mBE+a2.mHE+a3.mBE.mHE+a4.m + aS.m 2 2 IA2=bo+b.mBE+b2.mHE+b3.mBE.mHE+b4.mBE+b5.m HE qui donne en général une précision suffisante, et o mBE et mHE désignent les mesures à haute et basse énergie, et les " a " et " b " sont des coefficients qu'il faut calculer auparavant par une calibration. Cette calibration fait appel à un dispositif souvent appelé " fantôme " et qui est composé des matériaux de base évoqués ci-dessus et choisis pour simuler au mieux l'objet de mesure. Ces matériaux sont distribués dans le fantôme de façon à y présenter des régions o les longueurs, traversées par le rayonnement, des matériaux sont distribuées différemment. Un fantôme classique est construit sous la forme d'un escalier pour chacun des matériaux de mesure, et les escaliers sont superposés de facon que leurs marches soient perpendiculaires. En faisant traverser les escaliers par des rayons verticaux, on obtient ainsi toutes les combinaisons souhaitables entre les diverses épaisseurs des deux matériaux. Un autre fantôme est proposé dans le brevet américain 5 493 601 et comprend une série de tubes convergeant vers la source du faisceau et dotés des hauteurs inégalement réparties des deux matériaux. Il a B 13769.3 JCI pour but d'offrir des mesures plus exactes que le précédent, notamment grâce à la convergence des tubes vers la source, qui fait parfaitement coïncider la longueur parcourue par les rayons dans les tubes avec la hauteur totale de ceux-ci, et aussi pour réduire le rayonnement diffusé par le contenu des tubes. Les mesures expriment donc directement l'atténuation du rayonnement à travers le fantôme et permettent de caleuler sans souci les coefficients recherchés. Il n'en va pourtant pas de même pour les mesures réalisoss à travers l'objet à radiographier, pour lesquelles les conditions favorables exposées ci-dessus ne peuvent être réunies et le rayonnement diffusé ne doit pas être nagligé; il est même particulièrement important pour les sources d'un faisceau conique, souvent plus que le rayonnement primaire, qui a accompli un trajet rectiligne de la source au détecteur et est seul utile à la mesure. L' estimation puis la correction de ce rayonnement diffusé nécessitent des traitements spécifiques. C'est ainsi qu'il est possible, parmi bien d'autres méthodes, de le mesurer séparément puis de le soustraire du rayonnement total requ par le détecteur. On dispose pour cela d'un réseau d'éléments absorLants, comme des billes de plomb disposées au milieu d'une grille, entre l'objet radiographié et le réseau de détecteurs. Les rayons devant passer par ces billes sont interceptés totalement, de sorte qu'il ne parvient aux détecteurs situé sur le chemin de ces rayons que le rayonnement diffusé. Des interpolations numériques permettent ensuite d'estimer, avec une précision suffisante, la B 13769.3 JCI répartition du rayonnement diffusé sur l'ensemble du réseau de détecteurs. Ce procédé commode présente malheureusement ltinconvénient de devoir soumettre l'objet à une deuxième irradiation, susceptible d'être mal acceptée pour des êtres vivants. On peut aussi, au contraire, empêcher le rayonnement diffusé de parvenir aux détecteurs en les munissant d'une collimation stricte qui l'intercepte. Les procédés faisant appel à cette technique obligent à recourir à un balayage de la source qui implique un temps d'acquisition long, et le risque que l'objet ne bouge au cours de la mesure. Des méthodes numériques de genres divers ont encore été proposées pour corriger l' influence du rayonnement diffusé, mais elles ne conviennent généralement bien qu'à des objets particuliers. On doit conclure que le brevet cité plus haut ne semble pas proposer un progrès sensible, puisque le rayonnement diffusé n'est écarté des mesures que pour l'étape de calibration. L'objet de l' invention est de corriger l' influence du rayonnement diffusé d'une autre manière, en le laissant apparaître à la calibration, mais de façon analogue à ce qu'il fait à la radiographie de l'objet, ce qui permet de le corriger par un même procédé numérique, sans qu'on ait à craindre des discordances dans la qualité de ce procédé. Pour résumer, l'invention concerne avant tout un dispositif de calibration d'un système de radiographie à double énergie, comprenant un jeu de blocs présentant différentes épaisseurs d'un premier B 13769.3 JCI matériau, caractérisé en ce que les blocs sont munis d'évidements et en ce que le dispositif comprend encore des inserts à combler les évidements et comprenant différentes répartitions de hauteurs, les hauteurs et les épaisseurs étant considéréss dans une direction identique, entre le premier matériau et un second matériau; ainsi qu'un procédé de radiographie à faisceau conique à double énergie, comprenant une estimation d'épaisseurs de matériaux d'un objet radiographique par une combinaison numérique de mesures d'atténuation des énergies, impliquant une calibration de coefficients de la combinaison, caractérisé en ce que la calibration est mesurée par une radiographie d'un dispositif de calibration conforme à l'une quelconque des revendications précédentes, et quun rayonnement diffusé affectant la radiographie du dispositif de calibration est estimé tout en donnant un critère d'estimation utilisé ensuite pour estimer un rayonnement diffusé affectant la radiographie de l'objet. L'invention sera maintenant décrite au moyen des figures: - la figure 1 représente un organigramme de l'invention; et - les figures 2 et 3 représentent une vue de côté et une vue de dessus d'un fantôme de calibration utilisé pour satisfaire aux exigences exposées plus haut. On se reporte à la figure 1. A l'étape E1, une acquisition de calibration est faite par une radiographie à travers le fantôme B 13769.3 JCI décrit plus loin, ce qui donne des mesures brutes à haute et basse énergie. A l'étape E2, un traitement préliminaire d' estimation et de soustraction du rayonnement diffusé de ces mesures est fait pour calculer les mesures corrigées mHe et mBe des formules précédentes. Ensuite, l'étape E3 consiste en une estimation des paramètres du modèle, notamment les coefficients des polynômes a et b à partir des mesures et des épaisseurs connues que le rayonnement traverse dans le fantôme. D'autre part, une acquisition, à travers l'objet de mesure, est faite à l'étape E4, puis un traitement préliminaire est fait à l'étape E5; il est semblable à celui de l'étape E2 et son but est également d'estimer, puis de soustraire le rayonnement diffusé qui affecte les mesures. L'étape finale E6 consiste à exploiter les mesures corrigées effectuées à travers l'objet et d'appliquer le modèle, et avant tout les coefficients a et b calaulés à l'étape E3, afin de déduire les longueurs équivalentes des matériaux de base al et a2 dans l'objet. Le fantôme utilisé est représenté aux figures 2 et 3. Il comprend quelques blocs d'épaisseurs différentes d'un des matériaux de base, notamment le plexiglas. On pourrait utiliser un jeu de blocs séparés, à travers lesquels les mesures seraient faites successivement, mais on peut aussi unir les blocs en formant un solide à profil d'escalier 1 comprenant différentes couches, ici au nombre de quatre, et numérotées de 2 à 5 de bas en haut, les couches finissant, en général, par des faces dressées 6 ayant B 13769.3 JCI une orientation en dépouille, afin de réduire les composantes à haute fréquence du rayonnement diffusé. De plus, des rangées dinserts sont prévues sous les différents blocs, en comblant des évidements opérés ici à travers la couche inférieure 2. On rencontre, plus précisément, quatre rangées de deux inserts 7 chacune, numérotées de 8 à 11 et qui s'étendent respectivement sous la surface supérieure de la couche inférieure 2, sous la deuxième couche 3, sous les couches 3 et 4 et sous les trois couches supérieures 3 à 5. Les inserts 7 comprennent, en général, une portion en hydroxyapatite 12 et une portion en plexiglas 13. Les portions 12 et 13 ont une même hauteur totale, mais des hauteurs respectives différentes dans chaque rangée 8 à 11, de sorte que les rayons originaires de la source 14 et qui passent par chacun des inserts 7, traversent des combinaisons différentes d'épaisseurs des deux matériaux. De plus, il est produit un rayonnement diffusé analogue à celui qui est produit dans un être vivant. On est redevable de cette similitude à la constitution même du fantôme 1, puisque les matériaux qui le composent, outre qu'ils simulent bien les tissus mous et les os, ont des proportions et des répartitions analogues. En particulier, les inserts 7 sont suffisamment éloignés pour ne pas recevoir de rayonnement diffusé originaire des inserts voisins, mais seulement de la matière de base des couches 2 à 5. On décrira, pour la fin de cet exposé, un procédé d'estimation du rayonnement diffusé, utilisable pour les étapes E2 et E5; ce procédé fait cependant partie d'une autre invention et n'est donc donné que B 13769.3 JCI par souci de compléter la description de celle-ci et de prouver l'intérêt du fantôme 1. A l'étape E2, un rayonnement semblable à celui que subira l'objet et de flux 00 est projeté sur le fantôme 1, et un rayonnement total 0r, somme d'un rayonnement primaire 0 et du rayonnement diffusé 0, est capté par des détecteurs placés derrière le fantôme 1 et mesurant donc l'atténuation du rayonnement initial 00 à travers lui. En particulier, on obtient des mesures de rayonnement ayant passé par chacun des inserts 7, aussi bien pour la haute que pour la basse énergie. Les mêmes mesures sont répétées après avoir placé des absorbours, telles des billes de plomb, entre le fantôme et les détecteurs de facon à intercepter le rayonnement primaire et à ne plus mesurer que le rayonnement diffusé par les détecteurs. On en déduit la valeur du rayonnement diffusé 0a en fonction des hauteurs traversées des deux matériaux de base pour chacune des énergies. L'acquisition des mesures à travers l'objet donne aussi deux jeux de valeurs de rayonnement total 0t à haute et basse énergie. Le rayonnement 0 à travers l'objet peut être estimé d'après les mesures à travers le fantôme 1 et la relation entre 0t et 0 obtenue à travers lui pour les deux énergies, en supposant que la même relation s'applique à travers lobjet. Quoique l' invention trouve utilité pour les faisceaux coniques, son emploi n'est pas limité à eux; on peut concevoir que les étapes d' estimation et de B 13769.3 JCI correct10n flu rayonnement diffuse pourraient glee omises dans des situations oL il serait mains imporLanL, puisqu'il affecLerait les mesures gravers le fanL6me 1 at l'objeL pareillemenL, grace leur similitude de structure, eL que son influence disparaitraiL au calcul des longueurs quivalentes Al eL A,. B 13769.3 JC1 REVENDICATIONS 1. Dispositif (1) de calibration d'un système de radiographie à faisceau conique à double énergie, comprenant un jeu de blocs présentant différentes épaisseurs (2 à 5) d'un premier matériau, caractérisé en ce que les blocs sont munis d'évidements et en ce que le dispositif comprend encore des inserts (7) à combler les évidements et comprenant différentes répartitions de hauteurs (12, 13), les hauteurs et les épaisseurs étant considérées dans une direction identique, entre le premier matériau et un second matériau. 2. Dispositif de calibration suivant la revendication 1, caractérisé en ce que les blocs sont assemblés en une forme d'escalier et les inserts sont répartis en rangées (8 à 11) dans une couche inférieure (2) de l'escalier, les rangées étant situées sous des blocs différents. 3. Dispositif de calibration suivant la revendication 2, caractérisé en ce que l'escalier présente des faces dressées (6) en dépouille. 4. Procédé de radiagraphie à faisceau conique à double énergie, comprenant une estimation d'épaisseurs de matériaux d'un objet radiagraphique par une combinaison numérique de mesures d'atténuation des énergies, impliquant une calibration de coefficients (a, b) de la combinaison, caractérisé en ce que la calibration est mesurée par une radiographie d'un dispositif de calibration conforme à l'une quelconque des revendications précédentes, et qu'un rayonnement B 13769.3 JCI diffusé affectant la radiographic du dispositif de calibration est estimé tout en donnant un critère d'estimation utilisé ensuite pour estimer un rayonnement diffusé affectant la radiographic de l'objet.