FR 2841987 A1 20040109 FR 0208513 A 20020705 PROCEDE D'EXAMEN RADIOLOGIQUE D'UN OBJET DESCRIPTION Cette invention concerne un procédé d'examen radiologique multi-énergie d'un objet. Les procédés radiologiques consistent à faire traverser un objet à étudier par un rayonnement pour en déduire la répartition de différentes catégories de matériaux, absorbant le rayonnement de façons différentes, dans cet objet. Une application très courante est l'ostéodensitométrie, o on recherche la masse et la densité des tissus osseux dans un patient en distinguant ces tissus des tissus mous. Il est usuel d'utiliser un spectre de rayonnement large et à le diviser en bandes mesurées séparément par des canaux de mesure respectifs. Comme les coefficients d'absorption ou d'atténuation du rayonnement par l'une ou l'autre des catégories de tissu sont différents pour chacune des bandes, le problème théorique se ramène à la solution d'un système d'équations dont le nombre est égal à celui des bandes de mesure et qui comprennent chacune deux inconnues (les épaisseurs ou les masses traversées de tissus mous et osseux). La solution du problème devient possible après une calibration obtenue en faisant traverser par le rayonnement divers étalons dotés d'épaisseurs connues de matériaux aux propriétés d'absorption analogues à celles des matériaux de mesure, et notamment le plexiglas et l'hydroxyapatite pour simuler les tissus mous et les tissus osseux. On peut alors, calculer les paramètres mathématiques d'un modèle B 14155.3/JCI reliant les mesures d'atténuation aux épaisseurs des matériaux. On notera qu'en réalité le corps humain contient trois catégories principales de tissus: les tissus osseux, les tissus maigres et les tissus gras, mais qu'on n'en considère que deux en général à cause de difficultés à distinguer ces trois catégories dans les mesures, de sorte qu'on confond volontairement les tissus maigres et gras. D'autres procédés sont ensuite appliqués pour distinguer leurs proportions dans les tissus mous. La largeur du spectre permet de disposer d'un nombre beaucoup plus grand de bandes de mesures qu'il ne serait nécessaire pour obtenir une solution, et de les employer toutes pour obtenir des résultats plus précis en exploitant toute l'information d'absorption obtenue. Dans l'article "Measurement of bone mineral using a multiple-energy x-ray absorptiometry", par J. Swanpalmer, R. Kullenberg, T. Hansson, Phys.Med. Biol, Vol 43, 1997:pp 379-387, o on dispose de 23 bandes de mesures et o on considère les trois catégories de tissus, il est proposé de combiner les groupes de mesures trois par trois de toutes les façons possibles pour obtenir 1771 (23x22x21/2x3) systèmes de trois équations à trois inconnues qui donnent autant de groupes de résultats. Il convient alors d'appliquer un critère de choix. Les auteurs conseillent de choisir comme résultat véridique celui qui est à la valeur moyenne, ou à la valeur médiane, pour le paramètre le plus important, qui peut être la masse osseuse traversée. B 14155.3/JCI Contrairement à ce qu'on a conseillé plus haut, les auteurs considèrent les trois catégories de tissus. Cela ne met pas en question la validité de leur méthode, à condition d'exposer le patient à une intensité d'irradiation beaucoup plus forte pour ramener les incertitudes des mesures à des proportions acceptables. Il serait d'ailleurs possible de modifier simplement leur méthode pour l'appliquer à des mesures sur deux tissus, ou plus généralement deux catégories de matériaux, seulement. Toutefois, il subsiste le défaut que cette méthode ne tient pas compte du bruit sur les mesures et qu'elle conduit donc à des résultats bruités. L'invention a pour objet un perfectionnement de tels procédés par combinaison de résultats nombreux et comprend une amélioration du critère de choix des résultats. Plus précisément, elle concerne sous sa forme la plus générale un procédé d'examen radiologique d'un objet o sont considérées au moins deux catégories de matériaux, comprenant: l'emploi d'un rayonnement à spectre large; des mesures du rayonnement par bandes du spectre; des expressions (M) d'épaisseurs ou de masses des deux catégories de matériaux traversées par le rayonnement, les expressions (M) étant des fonctions d'au moins deux des mesures (Mesk) et de coefficient (A); et l'application d'un critère de choix parmi les expressions (M) pour en déduire une expression (M finale) jugée véridique; caractérisé en ce que le critère de choix comprend un combinaison (f) des expressions avec des facteurs de pondération (a), et un B 14155.3/JCI calcul de facteurs de pondération de manière que la combinaison ait un bruit minimal (variance minimale en langage mathématique) calculé d'après les bruits sur les mesures (variances sur les mesures). L'invention sera maintenant décrite en référence aux figures: - la figure 1 est une vue d'un spectre, - la figure 2 est une répartition de résultats, - et la figure 3 est un organigramme résumant l'invention, auquel on pourra se reporter dans toute la description qui suit. L'atténuation du rayonnement peut être exprimée par une fonction des épaisseurs traversées de chacun des matériaux d'indices x et y, ou de leurs masses M (densités par unité de surface) dans la direction du rayonnement. Le spectre de mesures de la figure 1 est divisé en N bandes notées généralement par les indices i et j. Les atténuations varieront dans chacune des bandes en raison de coefficients d'absorption variables pour les deux matériaux. Si nous appelons mesi ou mesj les mesures pour une bande d'énergie i ou j, les masses traversées M. et My pourront chacune être exprimées aussi, par exemple, par la formule générale M=A1+ A2.mesi+A3mesj+A4.mes 2 +A5.mes 2+AE.mesi.mesj.. les mesures considérées dans cet exemple étant des meusures d'atténuation, on aura pour chaque canal de mesure i (correspondant à une bande du spectre) la relation mesi = ln (n O1) o noi estle nombre de phatons arrivant sur ni l'objet et ni celui des photons ayant traversé l'objet. B 14155.3/JCI Comme la non-linéarité des fonctions M. et My en fonction des mesures est faible en pratique, on peut s'accommoder de cette fonction polynomiale du deuxième degré qui comprend six coefficients A1 à A6. Le degré du polynôme peut être ajusté en fonction du problème; Par exemple, pour l'analyse d'objets constitués de matériaux ayant des numéros atomiques plus élevés que celui des tissus biologiques, comme le contrôle non destructif pour l'examen d'un objet métallique.les mesures considérées dans cet exemple étant des meusures d'atténuation, on aura pour chaque canal de mesure i (correspondant à une bande du spectre) la relation mesi = in (Onio) o noi estle nombre de phatons arrivant sur 1ni l'objet et ni celui des photons ayant traversé l'objet. Ces coefficients peuvent être trouvés dans une étape d'étalonnage à travers des étalons, appelés parfois fantômes ou cales dans l'art, et qui consistent en des pièces formées d'épaisseurs connues et différentes entre elles de matériaux simulant, par leurs propriétés d'atténuation, les matériaux de l'objet qui sera effectivement à mesurer. Chacun des étalons est donc soumis à une irradiation par le rayonnement pendant une longue durée qui permet de réduire l'influence du bruit sur les mesures. Le spectre mesuré pour chacun des étalons donne encore N mesures résultant de la décomposition du spectre en autant de bandes. En combinant maintenant deux séries de mesures i et j prises pour deux bandes et pour chacun des étalons, on recherche les coefficients A pour ajuster les fonctions M, et My aux mesures. Dans le cas présent, o il y a six coefficients A pour chacune B 14155.3/JCI des deux fonctions, et o on considère les combinaisons de deux bandes de mesure, les mesures pour chacune des bandes devront, pour fournir une solution unique, porter sur six étalons. Des étalons plus nombreux pourront aussi être utilisés pour améliorer la précision sur la solution. Une minimisation de fonction d'erreur sera alors appliquée. Cette détermination des coefficients A est répétée pour diverses combinaisons de paires de mesures. Il en était encore ainsi dans l'article antérieur mentionné plus haut; cependant, on a constaté qu'il était inutile d'effectuer toutes les combinaisons, au nombre de (2) pour exploiter complètement les mesures et que (N-1) combinaison étaient en réalité suffisantes pour recueillir toute l'information. Une façon préférée de procéder consiste à choisir au départ la bande de mesures qui est la moins bruitée (par exemple celle qui a le signal n - le nombre de photons à la réception- le plus important) et à l'associer successivement à chacune des autres bandes de mesures pour les combinaisons. On obtient finalement (N-1) estimations des deux fonctions M, et My, qu'on note MM2,,AN-1 pour chacune de ces deux fonctions. A ce stade du procédé, les fonctions M. et My, représentatives des longueurs traversées de deux matériaux représentatifs des tissus osseux et des tissus mous, pourront être converties en fonction Mu, Mv et Mw représentatives des longueurs équivalentes traversées des tissus osseux, des tissus maigres et des B 14155.3/JCI tissus gras en combinant linéairement Mx et My de trois façons différentes déterminées par l'expérience. Ce procédé de conversion est indépendant de l'invention, déjà connu dans l'art. Comme il n'y a pas de raison de préférer l'une ou l'autre de ces estimations M, un critère de choix doit être appliqué pour obtenir l'estimation M finale qui sera jugée véridique. Dans l'article antérieur, une des expressions obtenues était directement sélectionnée d'après un critère de classement (le résultat médian) ou de moyenne des valeurs prises par les expressions pour un des résultats. Dans l'invention, les expressions M seront combinées, par exemple linéairement, d'après la formule M finale = (al M l)+(a2 M 2)+...+ (a N-1 M N -1) tout en minimisant le bruit; les coefficients a,, etc. ont une somme égale à l'unité (al + a2 +....+ a N-1 = 1). Pour chaque canal de mesure, le bruit susr le nombre de photons suit une loi statistique de Poisson, dont le résultat est indépendant pour chacun des canaux. La matrice de covariance des N-1 résultats peut être exprimée d'après la formule N-I _ aMj 1 lk- amesk amesk Nk La variance sur la combinaison linéaire donnant M finale s'exprime par la formule f =(al,..., aN 1) r.t(al...,aN1); cette quantité f atteint une valeur optimale quand sa dérivée selon toutes les variables B 14155.3/JCI est nulle, c'est-à-dire que l'influence du bruit est [af /a, = 0 I af/aa2 =0 minimisée, soit: * [af/aaN I =0 La résolution numérique de ce système fournit les coefficients ai, a2, etc et M finale, c'est-à-dire les masses traversées des deux catégories de matériaux. Ce sont ces opérations qui sont accomplies dans l'invention; il est important de remarquer que les mesures (mesk) employées dans la formule donnant rij et ensuite sont les mesures effectuées à travers l'objet à étudier lui-même, mais pas celles qui ont été faites à l'étalonnage pour déterminer les coefficients A des fonctions M. L'invention permet d'échapper au dilemme de procédés antérieurs dans lesquels on renonçait à une partie des énergies de mesure ou au contraire on acceptait des mesures peu distinctes: elle utilise en effet tout le spectre, mais en le divisant en bandes assez nombreuses pour que la mesure de chacune puisse être comparée de façon utile à d'autres mesures, faites sur des bandes éloignées. Elle permet ainsi, entre autres avantages, d'étudier aussi bien des organismes maigres que gras. Enfin, l'invention peut être généralisée à un nombre de matériaux plus grand que deux, ce qui peut être intéressant notamment pour les procédés d'imagerie à produit de contraste, o trois variables doivent être B 14155.3/JCI considérées ou pour contrôle des bagages (recherche d'explosifs); et elle peut être appliquée en considérant des combinaisons autres que linéaires des fonctions M. Des application de l'invention sont l'ostéodensitométrie - obtention de la densité de masse osseuse - obtention de la composition corporelle (masse grasse, masse maigre) le contrôle agro-alimentaire, par exemple la détection de fragments d'os dans de la viande ou la détection de morceaux de verres dans des plats cuisinés, le contrôle des bagages recherche d'explosifs, de produits illicites (armes, aliments, drogue...) B 14155.3/JCI REVENDICATIONS 1) Procédé d'examen radiologique d'un objet o sont considérées au moins deux catégories de matériaux, comprenant: l'emploi d'un rayonnement à spectre large; des mesures du rayonnement par bandes du spectre; des expressions (M) d'épaisseurs ou de masses des deux catégories de matériaux traversées par le rayonnement, les expressions (M) étant des fonctions d'au moins deux des mesures (mesk) et de coefficients (A); et l'application d'un critère de choix parmi les expressions (M) pour en déduire une expression (M finale) jugée véridique; caractérisé en ce que le critère de choix comprend une combinaison (f) des expressions avec des facteurs de pondération (a), et un calcul des facteurs de pondération de manière que la combinaison ait un bruit minimal calculé d'après un bruit sur les mesures. 2) Procédé selon la revendication 1, caractérisé en ce que la combinaison (f) des expressions est linéaire. 3) Procédé selon l'une quelconque des revendications 1 ou 2, caractérisé en ce que la variation de la combinaison est calculée avec une matrice de covariance (rf,) des expressions (M). 4) Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que les expressions (M) sont en nombre égal à celui des bandes, moins un, et établies toujours avec une des B 14155.3/JCI bandes (io) et, respectivement, chacune des autres bandes. ) Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que les coefficients (A) sont déterminés à une étape préalable de calibration. 6) Procédé selon l'une quelconque des revendications précédsente, caractérisé en ce qu'il est appliqué à l'ostéodensitométrie. 7) Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce qu'il est appliqué à des contrôles agro-alimentaires. 8) Procédé selon l'une quelconque des revendications 1 à 5, caractérisé en ce qu'il est appliqué au contrôle de bagages. B 14155.3/JCI