FR 3064075 A1 20180921 FR 1752153 A 20170316 1 MÉTHODE D'ESTIMATION DE LA DOSE DÉLIVRÉE PAR UN SYSTÈME DE RADIOLOGIE DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne de manière générale le domaine de la radiologie et plus particulièrement celui de la dosimétrie en radiologie. ÉTAT DE LA TECHNIQUE ANTÉRIEURE L'accroissement du nombre d'examens radiologiques (radiographie conventionnelle, scanographie, radiologie interventionnelle) subis par les patients est devenu un sujet majeur de santé publique. En effet, l'exposition aux rayonnements ionisants induit des dépôts d'énergie dans la matière et conduit à des effets radiobiologiques qui peuvent avoir des conséquences graves, de nature déterministe (radiodermite par exemple) ou stochastique (cancer par exemple). Il est par conséquent essentiel de pouvoir assurer une évaluation précise ainsi qu'un suivi des doses d'énergie délivrées aux patients au cours du temps. A ce jour, il existe principalement deux solutions pour estimer une dose délivrée en radiologie. La première solution est de mesurer le produit dose surface ou PDS par une chambre d'ionisation en sortie du tube à rayons X avec détermination de la surface d'entrée de cette chambre. Toutefois, cette mesure ne prend en compte que les paramètres physiques du tube et non la morphologie du patient. Une seconde solution est de calculer la dose à l'entrée (du patient) en fonction des paramètres physiques du tube et de la géométrie d'exposition. Toutefois, ce calcul est généralement très imprécis, ne tient pas compte de la morphologie du patient et ne permet pas de prendre en compte des variations de géométrie au cours du temps (scanographie ou radiologie interventionnelle par exemple). 3064075 2 On présentera dans un premier temps les grandeurs dosimétriques de référence. De manière générale, l'estimation des doses en radiologie fait appel à une grandeur dosimétrique caractéristique du faisceau, dénommée Kerma (Kinetic Energy 5 Released per unit Mass). Le kerma est défini par : K = dE tr (1) dm où dE, est l'énergie totale transférée (sous forme d'énergie cinétique) des 10 photons X aux particules chargées (électrons) dans un volume de matière de masse dm. Le kerma s'exprime en J.kg-l.s-1 ou, de manière équivalente en Gy. s-1. Le kerma est généralement constitué de deux composantes : le kerma de collision (K,01) et le kerma de radiation ( K rad ) mais dans la gamme d'énergie des photons X pour la radiologie le kerma de radiation peut être négligé par rapport au kerma de collision. 15 En radiologie conventionnelle, le kerma dans l'air peut être exprimé à partir des paramètres d'exposition du tube à rayons X. Plus précisément, le kerma dans l'air, en un point P donné, Ka (P), est proportionnel à la charge du tube utilisée, mAs , (c'est-à-dire le produit de l'intensité du courant traversant le tube par le temps d'exposition) et au carré de la tension accélératrice (1(17' ) des électrons dans le tube. Il est par ailleurs inversement 20 proportionnel au carré de la distance (dFp) du point en question au foyer de la source et à un facteur de filtrage ( F ) du tube, soit : Kair(P) oc mAs .(k17)2 , F 25 La dose absorbée représente l'énergie déposée localement par unité de masse dm. Elle est définie par : (2) 3064075 de D= dm où de est l'énergie déposée par rayonnement ionisant (direct ou indirect) dans un élément de matière de masse dm. Elle s'exprime en J.kg-1 ou, de manière équivalente, 5 en Gy. Comme indiqué précédemment, le kerma K représente l'énergie transférée alors que la dose absorbée, D, représente l'énergie absorbée. En règle générale, l'énergie n'est pas absorbée à l'endroit où elle est déposée en raison du parcours des électrons secondaires dans le matériau. Toutefois, dans le domaine de la radiologie, les photons X 10 étant de basse énergie (énergie inférieure à 150 keV), l'équilibre électronique s'établit à une faible profondeur dans le matériau et il est communément admis que D K. Un certain nombre de grandeurs dosimétriques, dites dérivées, sont obtenues à partir des grandeurs dosimétriques de référence. Il s'agit notamment du débit de kerma 15 dans l'air normalisé, noté nK, . Cette grandeur, encore appelée rendement du tube, représente le kerma par unité de charge électrique passant dans le tube et est donc définie par : nK = Kair (4) mAs 20 Où nKair s'exprime en Gy.(mAs)1 . On trouvera une description détaillée des grandeurs dosimétriques et de leurs relations dans le rapport de la Société Française de Physique Médicale, SFPM n°30, décembre 2014, intitulé « Dosimétrie des explorations diagnostiques en radiologie ». 25 Différents types de dosimètres ont été envisagés dans la littérature pour estimer la dose absorbée lors d'un examen radiologique. En particulier, la demande FR-A-2849697 propose d'utiliser un réseau de fibres optiques croisées couplées à un photomultiplicateur multi-canaux pour mesurer une dose lors d'une procédure radiologique interventionnelle. 3 (3) 3064075 4 Toutefois, ce dispositif est complexe et coûteux, et ne permet pas de mesurer avec justesse la dose absorbée par le patient. Le but de la présente invention est par conséquent de proposer une méthode de mesure qui soit simple, peu coûteuse et permette d'évaluer la dose absorbée par le patient 5 avec une marge d'erreur sensiblement inférieure à celle de l'état de la technique. EXPOSÉ DE L'INVENTION La présente invention est définie par une méthode d'estimation de dose délivrée à un patient par un système de radiologie équipé d'une source de rayons X et d'un capteur 10 plan, selon laquelle : (a) on acquiert une image radiologique du patient au moyen du capteur plan; (b) on détermine le niveau de gris, NG , du pixel au centre de l'image radiologique, correspondant à l'intersection de l'image avec l'axe du faisceau de rayons X; (c) on calcule le kerma Ke à l'entrée du patient à partir dudit niveau de gris, NG , 15 de la charge électrique traversant le tube de rayons X de la source, mAs , pendant le temps d'exposition du patient, de l'épaisseur 5 du corps du patient selon l'axe du faisceau de rayons X et de la distance di séparant la source de rayons X du patient. Le kerma à l'entrée du patient peut alors être calculé au moyen de l'expression : ( d'f di ) ( NG 20 Ke = B a+bln -c mAs rnAs où B est un facteur de rétrodiffusion prédéterminé, dref est une distance de référence à laquelle est effectuée la calibration de dosimétrie, a,b,c sont des coefficients d'une loi de calibration du débit de kerma dans l'air à ladite distance de référence et y est le coefficient 25 d'atténuation des rayons X dans l'eau. 3064075 5 Avantageusement, lors d'une première phase de calibration, on mesure, pour une pluralité de valeurs d'ampérage du tube et/ou de la tension accélératrice et/ou du temps d'exposition, d'une part le débit de kerma dans l'air, flKajr , à ladite distance de référence dref au moyen d'un dosimètre absolu et, d'autre part, le niveau de gris du pixel au centre 5 de l'image radiologique, normalisé par la charge électrique traversant le tube de rayons X, les coefficients a ,b ,c étant déterminés comme ceux vérifiant au mieux la relation nK air = a + b ln (NGn - c) pour ladite pluralité de valeurs. Dans une seconde phase de calibration, on interpose ensuite des fantômes d'eau de différentes épaisseurs x entre la source de rayons X et le capteur plan, et l'on détermine 10 pour chacune de ces épaisseurs le niveau de gris du pixel central de l'image radiologique, NG(x) , ainsi que le niveau de gris de ce pixel en absence de fantôme, NGn(0) , pour une même charge électrique traversant le tube de rayons X, puis l'on calcule les rapports logarithmiques ln NGn(x) pour les différentes épaisseurs, et enfin on en déduit le NGn(0) coefficient d'atténuation par régression linéaire sur ces rapports logarithmiques. 15 Selon une variante, préalablement à l'étape (b), on élimine de l'image radiologique le bruit de fond en lui soustrayant une image radiologique acquise en absence d'irradiation, pour obtenir une image radiologique débruitée. De sucroît, préalablement à l'étape (b) on corrige avantageusement les variations de gain pixel par pixel à partir de la dite image débruitée, de manière à obtenir une image 20 débruitée d'intensité uniforme pour une irradiation uniforme du capteur plan. BRÈVE DESCRIPTION DES DESSINS D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture d'un mode de réalisation préférentiel de l'invention en référence aux figures jointes parmi 25 lesquelles : La Fig. lA représente de manière de manière schématique une première configuration d'un système de radiologie lors d'une première étape de calibration de la méthode de dosimétrie selon l'invention ; 3064075 6 La Fig. 1B représente de manière schématique une deuxième configuration d'un système de radiologie lors d'une seconde étape de calibration de la méthode de dosimétrie selon l'invention ; La Fig. 1C représente de manière schématique une troisième configuration d'un 5 système de radiologie lors de la méthode de dosimétrie selon l'invention ; La Fig. 2 représente de manière schématique un organigramme de la méthode de de dosimétrie pour un système de radiologie, selon un mode de réalisation de l'invention. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS 10 L'idée à la base de l'invention est d'évaluer la dose délivrée à un patient à partir du niveau de gris au centre de l'image radiologique. Par centre de l'image radiologique, on entend le point d'intersection de l'axe du faisceau de rayons X avec le plan de l'image. L'image radiologique est obtenue au moyen d'un capteur plan capable de détecter un rayonnement X incident à faible énergie et de l'échantillonner à la fois spatialement et 15 en intensité. On distingue les capteurs plans à conversion directe et les capteurs plans à conversion indirecte. Les capteurs plans à conversion directe utilisent un matériau photoconducteur (par exemple le sélénium amorphe) pour convertir directement le rayonnement X incident en charges électriques qui sont collectées par une matrice de 20 transistors TFT. Les capteurs plans à conversion indirecte utilisent l'association d'un scintillateur et d'une matrice active de photodiodes (en silicium amorphe), le scintillateur convertissant les photons X en photons optiques et la matrice active convertissant les photons optiques en charges électriques. Les charges électriques ainsi générées sont collectées par une matrice de transistors TFT comme dans le capteur à conversion directe. 25 Le système d'imagerie est composé d'une source de rayonnement basse énergie (keV) et d'un capteur plan, associé avec un système d'acquisition. La source et le capteur plan sont montés sur deux bras robotiques (non représentés) et peuvent être mis en vis-à-vis. 3064075 7 Les bras robotiques peuvent réaliser une rotation (généralement de 3600) autour d'un point fixe appelé isocentre. Le système d'acquisition reçoit l'image radiologique brute du capteur plan, autrement dit les valeurs d'intensité des différents pixels. Il effectue tout d'abord une 5 correction des inhomogénéités d'offset du capteur plan, pixel par pixel. Plus précisément, il soustrait à chaque pixel l'intensité détectée en absence d'irradiation dite intensité DF (Dark Field). En d'autres termes, on soustrait à l'image radiologique brute une image radiologique obtenue sans irradiation de manière à obtenir une image débruitée. Il corrige ensuite, les inhomogénéités de gain, pixel par pixel, dues aux variations des paramètres 10 des amplificateurs dans la matrice de détection, de sorte que l'image débruitée soit d'intensité uniforme pour une irradiation uniforme du capteur plan. Un tel système d'acquisition et les corrections de gain qu'il applique sont bien connus de l'état de la technique. La méthode de dosimétrie selon la présente invention suppose d'effectuer 15 préalablement des étapes de calibration qui sont décrites ci-après en relation avec les Figs. 1A et 1B. La Fig. 1A représente de manière schématique une première configuration d'un système de radiologie lors d'une première étape de calibration de la méthode de 20 dosimétrie. On a désigné la source de rayons X par 110 et le capteur plan par 120. La source et le capteur plan sont montés en vis-à-vis sur un bras robotisé. On a illustré en 130 une table sur laquelle peut être allongé le patient durant l'intervention. On comprendra cependant que cette table est facultative et qu'elle peut être en outre remplacée par un 25 autre type de support permettant de positionner le patient. On supposera dans la suite que la table présente une faible atténuation aux rayons X basse énergie (keV), de sorte que son influence pourra être négligée en première approximation. Lors de la première étape de calibration, les mesures sont effectuées « à feu nu », c'est-à-dire sans atténuation entre la source et le capteur plan, autrement dit en absence de patient ou de fantôme de calibration. 3064075 8 On a pu montrer que le débit de kerma dans l'air, mesuré à une distance de référence, dref , de la source, pouvait s'exprimer sous la forme : nKair =a+bln(NGn-c) (5) 5 où NGn est le niveau de gris du pixel au centre, n, de l'image radiologique (fournie par le système d'acquisition), normalisé par la charge électrique passant dans le tube, soit : NG = NG mAs (6) 10 où NG est le niveau de gris au centre de l'image radiologique. Les coefficients a,b,c figurant dans la loi de calibration (5) sont des constantes qui ne dépendent que des caractéristiques du système de radiologie (source, capteur plan). Les coefficients a,b,c sont calculés à partir d'une pluralité de mesures de niveaux 15 de gris normalisés (NGn ) pour différents débits de kerma (nKair ), en minimisant l'erreur quadratique entre les mesures et la courbe représentant la loi de calibration (curve fitting). Les différents débits de kerma normalisés peuvent être obtenus par exemple en faisant varier l'ampérage du tube, la tension accélératrice et/ou le temps d'exposition. On comprendra que la loi de calibration (5) permet de relier un niveau de gris de 20 l'image radiologique et à un débit de kerma mesuré dans l'air, à une valeur de référence. La mesure de débit de kerma à la distance de référence peut être réalisée par exemple au moyen d'un dispositif de mesure de dose tel que le kit NOMEXTm de la société PTW, ayant bénéficié d'un étalonnage primaire. 25 La Fig. 1B représente de manière schématique une seconde configuration d'un système de radiologie lors d'une seconde étape de calibration de la méthode de dosimétrie. 3064075 9 Cette seconde configuration se distingue de la première par l'interposition d'un fantôme d'eau, 150, entre la source et le capteur plan, les autres éléments du système de radiologie restant identiques. Si l'on note x l'épaisseur du fantôme dans l'axe du faisceau, l'intensité /(x) du 5 faisceau atteignant le capteur plan au centre de l'image radiologique, après avoir traversé le fantôme d'eau, est donnée par : 1(x)= loe-tix (7) 10 où /./ est le coefficient d'atténuation du faisceau dans l'eau et /0 est l'intensité du faisceau au centre de l'image radiologique à feu nu (en absence de fantôme). On peut donc en déduire : NGn(x)= NGn(0)e-' (8) 15 où NGn(x) et NGn(0) sont les niveaux de gris normalisés du pixel au centre de l'image radiologique, n, respectivement en présence de fantôme d'eau et à feu nu. En utilisant des fantômes d'eau de différentes épaisseurs x et en calculant pour ces différentes épaisseurs les rapports logarithmiques ln NGn(x) , on peut en déduire par NGn (0) 20 régression linéaire le coefficient d'atténuation, y . Alternativement, on pourra effectuer la NG(x) , régression linéaire sur les rapports logarithmiques ln a condition de mesurer les NG(0) niveaux de gris NG(x) et NG(0) pour une même charge électrique traversant le tube de rayons X. 25 La Fig. 1C représente de manière schématique une troisième configuration du système de radiologie lors de la mesure de dose délivrée à un patient, conformément à un mode de réalisation de la méthode de dosimétrie selon l'invention. 3064075 10 La troisième configuration est identique à la seconde configuration à la différence près que le patient 170 est disposé entre la source de rayons X et le capteur plan, en lieu et place du fantôme d'eau. Le corps humain étant essentiellement constitué d'eau, on peut considérer en 5 première approximation que le coefficient d'atténuation des rayons X dans le patient est égal au coefficient d'atténuation /./ obtenu dans la seconde étape de calibration. A partir du niveau de gris du pixel au centre, n, de l'image radiologique du patient, on peut déterminer le débit de kerma dans l'air, mesuré à la distance de référence par : 10 nK air = a + b ln (NGn(0)-c)= a+ b ln (NGn(e)e's -c) (9) où NGn(S) est le niveau de gris normalisé du pixel en question, 6 étant l'épaisseur du corps du patient dans l'axe du faisceau. 15 Le kerma dans l'air à la distance de référence de la source est donc donné par : Kair =[a+b ln (NGn(S)e's -c)imAs (10) où mAs est la charge électrique traversant le tube pendant le temps d'exposition. 20 On peut en déduire le kerma, Ki, au point de contact Pi entre le corps du patient et la table, sur l'axe du faisceau. Si l'on note di la distance entre la source de rayons X et le point de contact Pi , on a, d'après la relation (2) : 25 Ki = K. ( dr f d Le kerma à l'entrée du patient, ou dose délivrée au patient, est la grandeur dosimétrique en radiologie conventionnelle. Elle prend en compte la dose, Ki, déposée 3064075 11 par le faisceau incident à la surface du patient mais aussi la dose déposée par le rayonnement rétrodiffusé à la surface du patient. On admet généralement que le kerma à l'entrée du patient, Ke , est le produit du kerma dû au faisceau incident par un facteur de rétrodiffusion, B, soit : Ke = K i.B (12) En définitive, le kerma à l'entrée du patient peut s'exprimer sous la forme suivante : ( d'f di ) Ke =B a +b ln ( NG mAs (13) -c mAs ) _ On rappelle que les coefficients a,b,c sont obtenus par la première étape de calibration et que le coefficient d'atténuation est obtenu par la seconde étape de 15 calibration. La charge électrique mAs passant par le tube est donnée comme le produit de l'ampérage traversant le tube et du temps d'exposition. L'épaisseur 6 pourra être estimée à partir de la corpulence du patient. Alternativement, elle pourra être estimée au moyen d'un dispositif de mesure de distance (équipé d'un capteur à ultrasons ou de lumière IR) intégré dans le capteur plan. Dans ce 20 cas, si l'on note du, la distance entre la source et le détecteur (capteur plan) et z la distance entre le capteur plan et le corps du patient selon l'axe du faisceau, mesurée par le dispositif en question, on aura simplement : Ô= - -z (14) 25 Le facteur de rétrodiffusion est généralement fixé à une valeur moyenne de référence (typiquement B est choisi égal à 1.35 environ ), ou bien tiré de tables obtenues par des simulations Monte-Carlo pour des qualités de faisceaux données. 5 10 3064075 12 Le niveau de gris du pixel au centre de l'image radiologique, NG , est fourni par le capteur plan ou, lorsque celui-ci existe, par le système d'acquisition. La Fig. 2 représente de manière schématique un organigramme de la méthode de 5 de dosimétrie pour un système de radiologie, selon un mode de réalisation de l'invention. Les étapes 210 à 225 sont des étapes de calibration proprement dites et effectuées préalablement à la mesure de dose proprement dite. En 210, dans une configuration du système à feu nu, on effectue des mesures de niveau de gris du pixel au centre, n, de l'image radiologique, ce pour différentes valeurs 10 d'ampérage du tube, de tension accélératrice et/ou différents temps d'exposition, et l'on en déduit les niveaux de gris normalisés correspondants, nNG . Corrélativement, on mesure, au moyen d'un dosimètre absolu, le kerma dans l'air à une distance de référence, dref , de la source de rayons X, et l'on en déduit le débit de kerma nK = Ka' , dans les mêmes conditions.ajrmAs 15 A l'étape 215, on déduit des mesures de nK air et nNG , les coefficients de la loi de calibration (5). A l'étape 220, on mesure le niveau de gris normalisé du pixel au centre, n, de l'image radiologique, NGn , pour des fantômes d'eau de différentes épaisseurs x interposés entre la source de rayons X et le capteur plan, selon la seconde configuration du 20 système de radiologie. A l'étape 225, on déduit des mesures NGn précédentes, le coefficient d'atténuation dans l'eau, u, des rayons X émis par la source. Dans la suite le patient est interposé entre la source et le capteur plan. A l'étape 230, on mesure la distance di de la source au patient le long de l'axe du 25 faisceau. A l'étape 240, on estime ou l'on mesure l'épaisseur 6 du patient le long de l'axe du faisceau. Enfin, à l'étape 250, on estime la dose délivrée au patient à partir du niveau de gris du pixel au centre de l'image radiologique, NG , au moyen de l'expression (13). 3064075 13 On a pu montrer que les doses mesurées à partir du niveau de gris NG étaient très proches des valeurs réelles, avec une marge d'erreur d'au plus 5%-10% environ, soit sensiblement moins que les marges d'erreur généralement constatées pour les méthodes de dosimétrie de l'état de la technique (30%). 5 De manière plus générale, l'homme du métier comprendra qu'après la phase de calibration, la méthode exposée plus haut permet de convertir une image radiologique, en niveaux de gris, fournie par le système d'acquisition, en image en dose. Pour ce faire, on peut remplacer dans l'expression (13), l'épaisseur Ô selon l'axe du faisceau par l'épaisseur 5/cos0 où 0 est l'angle que fait le rayon joignant la source et le pixel avec l'axe du 10 faisceau. L'image en dose ainsi calculée donne en chaque point la dose délivrée au patient. REVENDICATIONS1. Méthode d'estimation de dose délivrée à un patient par un système de radiologie équipé d'une source de rayons X et d'un capteur plan, caractérisée en ce que : (a) on acquiert une image radiologique du patient au moyen du capteur plan; (b) on détermine le niveau de gris, NG , du pixel au centre de l'image radiologique, correspondant à l'intersection de l'image avec l'axe du faisceau de rayons X; (c) on calcule le kerma Ke à l'entrée du patient à partir dudit niveau de gris, NG , de la charge électrique traversant le tube de rayons X de la source, mAs , pendant le temps d'exposition du patient, de l'épaisseur 6 du corps du patient selon l'axe du faisceau de rayons X et de la distance di séparant la source de rayons X du patient. 2. Méthode d'estimation de dose délivrée à un patient selon la revendication 1, caractérisée en ce que le kerma à l'entrée du patient est calculé au moyen de l'expression : ( d'f di ) Ke =B a+b ln ( NG mAs -c rnAs )_ où B est un facteur de rétrodiffusion prédéterminé, dref est une distance de référence à laquelle est effectuée la calibration de dosimétrie, a,b,c sont des coefficients d'une loi de calibration du débit de kerma dans l'air à ladite distance de référence et y est le coefficient d'atténuation des rayons X dans l'eau. 3. Méthode d'estimation de dose délivrée à un patient selon la revendication 2, caractérisée en ce que, lors d'une première phase de calibration, on mesure, pour une pluralité de valeurs d'ampérage du tube et/ou de la tension accélératrice et/ou du temps d'exposition, d'une part le débit de kerma dans l'air, nKair, à ladite distance de référence 3064075 15 dref au moyen d'un dosimètre absolu et, d'autre part, le niveau de gris du pixel au centre de l'image radiologique, normalisé par la charge électrique traversant le tube de rayons X, les coefficients a ,b ,c étant déterminés comme ceux vérifiant au mieux la relation nK air = a + b ln (NGn - c) pour ladite pluralité de valeurs. 5 4. Méthode d'estimation de dose selon la revendication 2 ou 3, caractérisée en ce que, dans une seconde phase de calibration, on interpose des fantômes d'eau de différentes épaisseurs x entre la source de rayons X et le capteur plan, que l'on détermine pour chacune de ces épaisseurs le niveau de gris du pixel central de l'image radiologique, 10 NG(x) , ainsi que le niveau de gris de ce pixel en absence de fantôme, NGn(0) , pour une même charge électrique traversant le tube de rayons X, que l'on calcule les rapports logarithmiques ln NGn(x) pour les différentes épaisseurs, et que l'on en déduit le NGn(0) coefficient d'atténuation par régression linéaire sur ces rapports logarithmiques. 15 5. Méthode d'estimation de dose selon l'une des revendications précédentes, caractérisée en ce que, préalablement à l'étape (b), on élimine de l'image radiologique le bruit de fond en lui soustrayant une image radiologique acquise en absence d'irradiation, pour obtenir une image radiologique débruitée. 20 6. Méthode d'estimation de dose selon la revendication 5, caractérisée en ce que, préalablement à l'étape (b), on corrige les variations de gain pixel par pixel à partir de la dite image débruitée, de manière à obtenir une image débruitée d'intensité uniforme pour une irradiation uniforme du capteur plan. 25