Brevet FR3035785A1 — publié le 20161111

L'invention concerne le domaine des procédés d'estimation de la répartition de la densité minérale osseuse dans au moins une partie de squelette d'un individu. ARRIERE PLAN DE L'INVENTION Il est connu des procédés pour caractériser la répartition de la densité minérale osseuse dans le squelette d'individus. Cette caractérisation peut être particulièrement utile pour diagnostiquer un risque d'ostéoporose d'un individu. A cet effet on utilise des images radiologiques du squelette de l'individu réalisées en projetant depuis une source de rayons X, tel qu'un tube à rayons X, des rayons X à différents niveaux d'énergie sur l'individu et en détectant, pour chacun de ces niveaux d'énergie, les atténuations de ces rayons X par les tissus de l'individu. L'analyse de l'atténuation différentielle entre les tissus permet de visualiser sur un plan de projection du squelette de l'individu, une répartition surfacique de la densité minérale osseuse. Un inconvénient de ce procédé connu sous le nom de « Dual X-Ray absortiometry » pour absorptiométrie bi photonique à rayons X, est que lorsque l'appareil radiologique nécessite un temps entre la génération des deux niveaux d'énergie, l'individu doit rester immobile, sans quoi on ne peut pas réaliser d'analyse différentielle entre les deux images obtenues aux deux niveaux d'énergie. Un autre problème de cette méthode est qu'il nécessite un niveau de formation important pour les opérateurs de radiologie ainsi que des installations spécifiques à la génération de rayons X bi énergie. OBJET DE L'INVENTION Un objet de la présente invention est de fournir 3035785 2 un procédé alternatif d'estimation de la répartition de la densité minérale osseuse dans au moins une partie de squelette d'un individu. 5 RESUME DE L'INVENTION Afin de répondre à cet objectif, il est proposé selon l'invention un procédé d'estimation de la répartition de la densité minérale osseuse dans au moins 10 une partie de squelette d'un individu, ce squelette étant composé de tissus osseux, ces tissus osseux étant entourés de tissus mous eux même entourés par une enveloppe externe de l'individu. Ce procédé est essentiellement caractérisé en ce 15 qu'il comprend la génération : - d'une représentation géométrique en trois dimensions de ladite partie du squelette de l'individu; et - d'une représentation géométrique en trois 20 dimensions d'au moins une portion de l'enveloppe externe contenant ladite partie du squelette de l'individu ; - d'au moins une première image radiologique de ladite au moins une partie du squelette, cette première image étant générée par observation d'une projection de 25 rayons X, provenant d'une première source, sur une première surface de détection, cette première image étant générée par projection de rayons X et représentant l'individu observé suivant une première incidence d'observation ; 30 - d'au moins une deuxième image radiologique de ladite au moins une partie de squelette, cette deuxième image étant générée par projection de rayons X et représentant l'individu observé suivant une deuxième incidence d'observation distincte de la première 35 incidence d'observation; puis 3035785 3 à l'aide d'un logiciel prédéterminé, et pour plusieurs rayons X donnés parmi lesdits rayons X projetés, on estime pour chaque rayon X donné : - une valeur d'épaisseur de tissus mous 5 traversés par ce rayon X donné ; - une valeur d'épaisseur de tissus osseux traversés par ce rayon X donné ; et - une valeur d'atténuation globale de ce rayon X donné, entre sa source et sa projection sur 10 la surface de détection où il est projeté ; et - à l'aide de la valeur d'épaisseur de tissus mous traversés par ce rayon X donné, et d'une estimée d'un coefficient d'atténuation des tissus mous traversés par ce rayon X donné, on déduit d'une 15 part une valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus mous qu'il a traversée et d'autre part une valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus osseux qu'il a traversée ; puis 20 - à l'aide de la valeur d'épaisseur de tissus osseux traversée par ce rayon X donné, de la valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus osseux qu'il a traversée, et de la valeur d'atténuation globale de ce rayon X 25 donné, le logiciel estime une valeur représentative de la densité minérale osseuse des tissus osseux traversés par ce rayon X donné. Pour la compréhension de l'invention, le terme 30 rayon X correspond à au moins un rayon X ou plusieurs rayon X parallèles entre eux et formant un faisceau. Toujours pour la compréhension de l'invention, le terme coefficient d'atténuation désigne un coefficient permettant de quantifier l'atténuation de l'intensité 35 d'un rayon X donné par une matière donnée en fonction de 3035785 4 l'épaisseur de cette matière traversée par ce rayon X. Chaque rayon X projeté depuis une source de projection se propage dans l'espace et jusqu'à la surface de détection qui lui correspond en suivant une 5 trajectoire rectiligne. Ainsi chaque rayon projeté traverse les tissus mous et les tissus osseux du squelette de l'individu en suivant une ligne droite. Connaissant la source du rayon et son point de projection, c'est-à-dire son point d'impact sur la 10 surface de détection, on peut déterminer la position dans l'espace de chaque rayon X par rapport aux surfaces de détection et par rapport aux sources, sachant que ces surfaces de détection et sources sont fixes les unes par rapport aux autres au moins pendant le temps nécessaire à 15 la projection et que les positions relatives entre ces surfaces de détection et les sources sont connues, notamment par calibration de l'environnement radiographique. Sachant que chaque rayon X traversant une matière 20 est atténué par cette matière en fonction d'un coefficient d'atténuation propre à cette matière, on constate que chaque rayon X traversant l'enveloppe de l'individu, les tissus mous et/ou tissus osseux de l'individu est atténué par l'ensemble des tissus qu'il a 25 traversé. Comme chaque point de l'image radiologique présente un niveau de gris ou de couleur ou de contraste qui est représentatif / fonction de l'intensité du rayon 30 X arrivant en un point correspondant de la surface de détection, on peut déterminer par observation de chaque point de l'image radiologique, un niveau d'intensité du ou des rayons X qui ont été projetés sur le point correspondant de la surface de détection correspondante. 35 3035785 5 Connaissant l'intensité de chaque rayon X émis par la source donnée qui lui correspond et connaissant son niveau d'intensité à son arrivée sur la surface de détection correspondante, on peut retrouver à partir de 5 chaque point d'une image radiologique, un niveau d'atténuation subie par le ou les rayons X qui sont arrivés / ont été projetés sur ce point depuis la source correspondante. Ainsi, par contraste entre les points de l'image, 10 on obtient les limites projetées des différents tissus mous, des différents tissus osseux et de l'enveloppe de l'individu. Comme les sources des rayons et les surfaces de détection sont fixes dans l'espace et localisées dans 15 l'espace à l'aide d'un repère spatial et comme l'individu conserve la même posture pendant le temps de la projection, on peut en observant les première et deuxième images radiologiques, déduire aux incertitudes quantifiables près, les positions exactes dans ledit 20 référentiel spatial de chaque point de l'individu qui a été traversé : - d'une part par un rayon X ou faisceau de rayons X provenant de la première source et projeté(s) sur la première surface de détection 25 comme illustré sur la première image radiologique; et - d'autre part par un autre rayon X ou faisceau de rayons X utilisé pour générer la deuxième image radiologique. 30 On note que le changement d'incidence d'observation entre les première et deuxième images radiologiques peut être obtenu : - soit en déplaçant, la source de rayons X 35 par rapport à l'individu observé, tout en s'assurant 3035785 qu'il conserve une même posture entre les prises de vues nécessaires à l'obtention des première et deuxième images ; - soit en utilisant deux sources de rayons X 5 distinctes et/ou deux surfaces de détection distinctes. Ainsi, l'analyse des première et deuxième images radiologiques, réalisées selon des incidences d'observation de l'individu distinctes, c'est à dire 10 différentes l'une de l'autre, permet : - d'une part d'avoir une représentation spatiale de l'individu, de son enveloppe et de ses tissus mous et osseux ; et - d'autre part d'avoir une estimation de 15 l'atténuation globale subie par un rayon traversant plusieurs tissus mous et éventuellement plusieurs tissus osseux de l'individu. Toutefois la seule estimation d'une valeur d'atténuation globale du rayon X donné, entre sa source 20 et sa projection sur la surface de détection, n'est pas suffisante pour estimer la densité osseuse des tissus osseux traversés par le rayon X donné. En effet, une part de cette valeur d'atténuation globale du rayon X donné est liée non pas aux tissus osseux, mais à des tissus 25 mous. Le procédé selon l'invention apporte une solution à ce problème puisqu'il permet d'estimer la part d'atténuation d'un rayon X donné qui est liée aux seuls 30 tissus mous traversés par ce rayon X donné et la part d'atténuation de ce rayon X donné qui est liée aux tissus osseux traversés par ce rayon X. En effet, le procédé selon l'invention utilise : - une représentation géométrique en trois 35 dimensions de ladite partie du squelette de 3035785 7 l'individu, et - une représentation géométrique en trois dimensions d'au moins une portion de l'enveloppe externe contenant ladite partie du squelette de 5 l'individu ; et connaissant la position de chaque rayon X donné par rapport à ces représentations géométriques de la partie du squelette et de l'enveloppe externe, on sait calculer à l'aide du logiciel : 10 - d'une part la distance totale parcourue par ce rayon X donné au travers des tissus mous situés dans l'enveloppe de l'individu ; et - d'autre part la distance totale parcourue même rayon X donné au travers des tissus situés dans l'enveloppe de l'individu. par ce 15 osseux Ces deux distances correspondent respectivement - une valeur d'épaisseur de tissus mous 20 traversés par ce rayon X donné ; et - une valeur d'épaisseur de tissus osseux traversés par ce rayon X donné. Il existe plusieurs moyens pour estimer le 25 coefficient d'atténuation des tissus mous de l'individu. On peut par exemple rechercher un rayon X particulier traversant uniquement des tissus mous de l'individu sans traverser de tissus osseux. 30 Pour estimer le coefficient d'atténuation des tissus mous de l'individu, le logiciel recherche un rayon X particulier traversant uniquement des tissus mous de l'individu sans traverser de tissus osseux et à partir de la position de ce rayon X particulier et de la 35 représentation géométrique en trois dimensions de la 3035785 8 portion de l'enveloppe externe, le logiciel estime la valeur d'épaisseur de tissus mous traversés par ce rayon X particulier, et par observation d'un point de projection de ce rayon x particulier, le logiciel 5 détermine une valeur d'atténuation subie par ce rayon X particulier lors de sa traversée des seuls tissus mous de l'individu, puis connaissant ladite valeur d'épaisseur de tissus mous traversés par ce rayon ladite valeur d'atténuation subie 10 particulier, le logiciel calcule d'atténuation associé aux seuls l'individu. X particulier et par ce rayon X le coefficient tissus MOUS de Il s'agit ici d'un coefficient d'atténuation 15 spécifique aux tissus mous de l'individu. Il est bien entendu tout à fait possible, pour améliorer la précision, de calculer ce coefficient pour différents rayons X ne traversant que les tissus mous, dans une région considérée comme homogène, et d'estimer 20 le coefficient d'atténuation comme une moyenne des différentes valeurs obtenues. Puis, pour tout rayon X donné, connaissant : - le coefficient d'atténuation des tissus mous de l'individu ; et 25 - l'épaisseur de tissus mous traversée par ce rayon X donné, le logiciel calcule : - d'une part une valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus mous qu'il a traversée ; et 30 - d'autre part, en utilisant une intensité IO de ce rayon donné à sa source et une intensité I de ce rayon donné telle que détectée sur la surface de détection sur laquelle il est projeté, le logiciel calcule une valeur d'atténuation globale de ce rayon X donné, entre sa 35 source et la surface de détection sur laquelle il est 3035785 9 projeté ; puis - en corrigeant la valeur d'atténuation globale de ce rayon X donné pour tenir compte de ladite valeur d'atténuation de ce rayon X donné par l'épaisseur de 5 tissus mous qu'il a traversée, le logiciel calcule la valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus osseux qu'il a traversée. Ainsi, pour chaque rayon X donné, on peut 10 calculer la valeur d'atténuation du rayon X donné par la seule épaisseur de tissus osseux qu'il a traversée. On peut décider que ladite valeur représentative de la densité minérale osseuse des tissus osseux 15 traversés par ce rayon X donné soit un coefficient d'atténuation du tissu osseux traversé par ce rayon x donné. Comme on le verra en détail ci-après, connaissant 20 d'une part l'épaisseur du tissus osseux traversé par ce rayon X donné et d'autre part le niveau d'atténuation de ce rayon X imputable à la seule traversée de ce tissus osseux, il est possible de calculer une densité minérale osseuse dite surfacique qui correspond aux tissus osseux 25 qui ont été traversés par ce rayon X donné. Connaissant les coefficients d'atténuation des tissus osseux estimés avec les différents rayons traversant la structure osseuse d'intérêt, caractérisant l'image projetée de cette structure, il est possible de calculer la densité 30 minérale osseuse surfacique. Il est à noter qu'en complément de la méthode précitée de détermination du coefficient d'atténuation associé aux seuls tissus mous, on peut aussi utiliser une 35 méthode nécessitant une éprouvette de calibration 3035785 10 présentant plusieurs zones ayant chacune un coefficient d'atténuation propre connu. Cette méthode consiste à placer une ou plusieurs éprouvettes qui présentent chacune au moins un 5 coefficient d'atténuation connu, à côté de l'individu lors de la projection des rayons depuis les sources vers les surfaces de détection. Puis, en étudiant le niveau de gris ou de couleur ou de contraste sur les première et seconde images 10 radiologiques, on peut obtenir une correspondance entre chaque coefficient d'atténuation propre connu d'une éprouvette et un niveau gris ou de couleur ou de contraste correspondant sur l'image. On peut ainsi simplement estimer les atténuations 15 subies par chaque rayon X donné projeté sur une surfaces de détection donnée, en recherchant dans une zone de projection de l'éprouvette visible sur une des première ou deuxième images de radiologie, un niveau de gris, de couleur ou de contraste associée équivalent au niveau de 20 gris, de couleur ou de contraste constaté en un point de l'image correspondant au point de contact entre ce rayon X donné et la surface de détection correspondante. On peut ainsi déterminer l'atténuation d'intensité subie par chaque rayon X donné après avoir traversé les tissus de 25 l'individu. En effet, deux points de l'image présentant une intensité similaire ont une même valeur de produit mu*épaisseur, mu étant le coefficient d'atténuation 30 propre à la matière traversée par un rayon X et l'épaisseur étant l'épaisseur de cette matière traversée par le rayon. Ainsi, connaissant l'épaisseur du milieu traversé et l'épaisseur de l'éprouvette on peut estimer le mu des tissus mous traversés. 35 Par correction, on peut ensuite estimer le 3035785 11 coefficient moyen du tissu dur traversé. Dans un mode de réalisation préférentiel, afin d'obtenir les deux images avec les première et deuxième 5 incidences différentes, on peut faire en sorte que ladite deuxième image radiologique de ladite au moins une partie de squelette soit générée par observation d'une projection de rayons X, provenant d'une deuxième source distincte et éloignée de la première source, sur une 10 deuxième surface de détection également distincte et éloignée de la première surface de détection. Dans un mode particulier de réalisation de l'invention, ladite représentation géométrique en trois 15 dimensions de ladite au moins une portion de l'enveloppe externe contenant ladite partie du squelette de l'individu est déduite à l'aide d'au moins une observation de l'enveloppe réalisée à l'aide de moyens de capture optique. 20 Idéalement, la représentation géométrique en trois dimensions de l'enveloppe de l'individu ainsi que les première et deuxième images radiologiques respectivement positionnées dans un même référentiel 25 spatial et correspondent préférentielle ment à un même instant d'observation de l'individu, ce qui permet de faciliter les calculs nécessaires à la mesure de répartition de densité minérale osseuse dans la partie observée du squelette. 30 Enfin, il est à noter que dans certains cas on peut générer la représentation géométrique en trois dimensions de l'enveloppe de l'individu en utilisant les première et deuxième images radiologiques. 35 Pour cela, le logiciel identifie sur les première sont 3035785 12 et deuxième images radiologiques, des paramètres caractéristiques de l'individu tels que positions spatiales relatives entre des points caractéristiques du squelette et/ou de l'enveloppe externe de l'individu 5 visible sur ces première et deuxième images, ou des dimensions ou des distances entre points caractéristiques du squelette et/ou de l'enveloppe externe visibles sur ces première et deuxième images. Puis, le logiciel génère une représentation 10 géométrique en trois dimensions de l'enveloppe de l'individu, à l'aide des paramètres caractéristiques de l'individu ainsi identifiés et à l'aide d'une base de données statistique préalablement élaborée avec une population d'individus, cette base de données illustrant 15 des relations entre des paramètres caractéristiques d'individus et des représentations escomptées en trois dimensions de l'enveloppe de l'individu. Dans un mode particulier lesdites première et deuxième images correspondent respectivement à deux 20 instants distincts de détection des rayons X projetés sur les première et deuxième surfaces de détection. Par ailleurs, les première et deuxième images peuvent illustrer des projections de l'individu en position debout. 25 BREVE DESCRIPTION DES DESSINS D'autres caractéristiques et avantages de l'invention ressortiront clairement de la description qui en est faite ci-après, à titre indicatif et nullement limitatif, en référence aux dessins annexés, dans 30 lesquels: - La figure 1 présente une vue latérale d'un individu sur laquelle on voit l'enveloppe E de l'individu et un plan d'observation « a » choisi pour observer une partie du squelette de l'individu et de l'enveloppe, 35 ici en coupe transversale ; 3035785 13 - La figure la présente la vue en coupe transversale de la partie de squelette de l'individu selon le plan d'observation « a », dans le cas où les tissus mous de l'individu sont considérés comme 5 relativement homogènes du point de vu de leurs coefficients d'atténuation, on voit sur cette figure des première et deuxième sources Slat et Spa de projection de rayons X et des première et deuxième surfaces planes Pl, P2 (plans de détection) de détection de rayons X, chaque 10 surface de détection présentant une pluralité de points de détection chacun adapté à détecter une intensité de rayons X projeté sur le point de détection ; - La figure 2 illustre la même vue que celle de la figure la mais ici l'individu présente une 15 hétérogénéité H de tissus mous conduisant à la prise en compte de plusieurs coefficients d'atténuation, en l'occurrence un coefficient devant être déterminé pour chaque hétérogénéité de tissus mous pour pouvoir estimer la répartition de la densité minérale osseuse en limitant 20 le bruit généré par les tissus mous qui entourent les os et les hétérogénéités de ces tissus mous. DESCRIPTION DETAILLEE DE L'INVENTION Comme indiqué précédemment, l'invention porte sur un procédé d'estimation de la répartition de la 25 densité minérale osseuse dans au moins une partie V de squelette d'un individu. L'invention porte aussi sur un dispositif pour la mise en oeuvre de ce procédé. Comme on le voit sur les figures la et 2, ce dispositif peut comporter des 30 première et deuxième surfaces de détection Pl, P2 qui sont respectivement planes et sont orientées entre elles d'un angle de 900 avec une incertitude d'angle de plus ou moins 60° et préférentiellement de plus ou moins 50 d'angle et préférentiellement inférieure à 1° d'angle, 35 comme sur les figures la et 2. 3035785 14 Ce dispositif comporte également : - une première source de rayons X, dite source postéro antérieure Spa, pour projeter des rayons vers et sur la première surface de 5 détection Pl en traversant l'individu en allant de la face postérieure de l'individu vers sa face antérieure ; et - une deuxième source de rayons X, dite source latérale Slat, pour projeter des rayons vers 10 et sur la deuxième surface de détection P2 en traversant l'individu entre ses faces latérales. Le dispositif comporte aussi un calculateur et un logiciel permettant de réaliser des calculs ou des estimations ou déductions nécessaires à la mise en oeuvre 15 du procédé selon l'invention, ce dispositif étant ainsi automatisé ou du moins semi automatisé pour permettre une estimation de la répartition surfacique ou volumique de densité minérale osseuse, ci-après DMO, dans au moins une partie du squelette d'un individu. 20 A noter que le dispositif peut être constitué d'une source mobile et de moyens de calibration pour connaitre dans chacune des vues l'environnement radiologique global (position des sources par rapport aux films et position relatives des sources), et si 25 nécessaire d'un repérage sur le sujet permettant de positionner un sujet dans un environnement global. Le dispositif peut aussi être constitué d'une seule source et d'un dispositif de rotation et de calibration permettant de tourner le sujet et de définir 30 la paire d'images en prenant en compte la rotation du sujet pour avoir les deux incidences différentes. Le dispositif selon l'invention peut aussi comporter des moyens de capture optique d'observation de 35 l'enveloppe externe E de l'individu. Ces moyens de 3035785 15 capture optique qui sont reliés au calculateur permettent de générer une représentation géométrique en trois dimensions d'au moins une portion de l'enveloppe externe E contenant ladite partie du squelette de l'individu. 5 Ces moyens de capture optique, non représentés, peuvent être un scanner de corps ou une caméra ou appareil photo, stéréoscopique ou non. Ces moyens de capture optique peuvent éventuellement être associés à des moyens de projection 10 de franges de moiré ou de trame éventuellement générée via des projections de lumières sur l'enveloppe de l'individu. Ces moyens de projection peuvent être par exemple adaptés à projeter les franges ou trames sous forme de lumière laser. 15 Pour faciliter l'analyse des images obtenues par le dispositif de l'invention, les positions des première et deuxième surfaces de détection et les positions des première et deuxième sources de rayons X ainsi qu'éventuellement celles des moyens de capture optique 20 sont localisées dans un même espace à l'aide d'un repère spatial. L'objectif du procédé selon l'invention est de permettre la fourniture d'une estimée de la répartition 25 minérale osseuse DMO dans une partie de squelette d'un individu. La partie, plus particulièrement ciblée du squelette pour l'estimation de la répartition de DMO est essentiellement un corps vertébral ou une vertèbre, ou un 30 col du fémur ou plus généralement le fémur. Ces parties sont ciblées car elles sont particulièrement adaptées à la caractérisation du risque d'ostéoporose ou d'une évolution de densité minérale. 35 Dans un premier temps, pour la mise en oeuvre du 3035785 16 procédé selon l'invention, on génère : - une représentation géométrique en trois dimensions de ladite partie du squelette V de l'individu dont on veut étudier la répartition de DMO; et 5 - une représentation géométrique en trois dimensions d'au moins une portion de l'enveloppe externe E contenant ladite partie du squelette V de l'individu ; et - au moins une première image radiologique de 10 ladite au moins une partie du squelette, cette première image étant générée par observation d'une projection de rayons X, provenant de la première source Spa, sur la première surface de détection Pl; et - au moins une deuxième image radiologique de 15 ladite au moins une partie de squelette, cette deuxième image étant générée par observation d'une projection de rayons X, provenant d'une deuxième source Slat distincte et éloignée de la première source, sur une deuxième surface de détection. 20 Idéalement, ces première et deuxième images correspondent détection des deuxième surfaces 25 entend que l'on un même instant de X projetés sur les première et de détection. Par même instant, on a un décalage temporel entre les respectivement à rayons détections nécessaires à la génération des première et deuxième images qui est inférieur à 1 seconde et préférentiellement inférieur à 0.5 seconde. Idéalement, comme il s'agit d'obtenir une 30 répartition de DMO sur le squelette en vue, par exemple d'une analyse de risque d'ostéoporose, les première et deuxième images sont réalisées pour illustrer des projections de l'individu en position debout. La position debout de l'individu est la plus représentative des 35 risques associés à l'ostéoporose. 3035785 17 Bien qu'il soit possible d'obtenir ladite représentation géométrique en trois dimensions de ladite au moins une portion de l'enveloppe externe contenant ladite partie du squelette de l'individu à partir d'une 5 observation de l'individu réalisée avec les moyens de capture optique, il est aussi possible que cette représentation géométrique en trois dimensions soit déduite à partir desdites première et deuxième images radiologiques. 10 Dans certains cas, il est possible de combiner ces deux techniques, les moyens de capture optique permettant de s'assurer de la précision de représentation de l'enveloppe externe et l'usage des première et deuxième images permettant d'assurer que la 15 représentation de l'enveloppe soit bien positionnée dans le même référentiel spatial que la représentation en trois dimensions de la partie du squelette que l'on veut visualiser. Dans le cas où les moyens de capture optique ne 20 fourniraient qu'une seule vue de l'enveloppe externe et en considérant un modèle 3D de vertèbre à partir de la seule image radiographique estimée, il est alors possible d'estimer les épaisseurs de tissus mous traversés et ainsi évaluer l'atténuation des rayons X par ces tissus 25 mous. La DMO surfacique de la vertèbre peut ainsi être estimée mais avec une précision réduite. On note qu'il existe plusieurs solutions pour déterminer une représentation géométrique en trois 30 dimensions de ladite partie du squelette de l'individu. On notera par exemple les solutions illustrées dans les documents brevets FR2856170 et W02008146069. La publication Humbert L, De Guise JA, Aubert B, 35 Godbout B, Skalli W: 3D reconstruction of the spine from 3035785 18 biplanar X-rays using parametric models based on transversal and longitudinal inferences. Med Eng Phys. 2009 Jul;31(6):681-7 décrit un exemple de méthodologie pour reconstruire un modèle virtuel d'une partie de squelette d'un individu à l'aide d'une radiographie biplane, en générant des images du type desdites première et deuxième images radiologiques. Les deux publications suivantes traitent à la fois de la reconstruction 3D du squelette et de celle de 10 l'enveloppe externe. La publication "Personalized Body Segment Parameters From Biplanar Low-Dose Radiography" de Raphael Dumas, Rachid Aissaoui, Member, IEEE, David Mitton, Wafa 15 Skalli, et Jacques A. de Guise publiée au journal IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, Vol 52, NO. 10 Octobre 2005, décrit aussi une méthodologie pour reconstruire un modèle virtuel d'une partie de squelette d'un individu à l'aide d'une radiographie biplane, en 20 générant des images du type desdites première et deuxième images radiologiques. Cette méthodologie est incorporée par référence car elle donne un exemple d'obtention d'un modèle virtuel représentatif d'une partie de squelette / représentation géométrique en trois dimensions de ladite 25 partie du squelette de l'individu. Ce type de modèle obtenu conformément à la méthode décrite dans cette publication est utilisable pour la mise en oeuvre de la présente invention. 30 De même, la publication "Subject-specific body segment parameters estimation using biplanar X-rays: a feasibility study" de Baptiste Sandoz, Sébastien Laporte, Wafa Skalli and David Mitton CNRS, Arts et Métiers Paris Tech, LBM, 151 bd de l'Hôpital 75013 Paris, France, 35 publiée au journal Computer Methods in Biomechanics and 3035785 19 Biomedical Engineering Vol. 13, No. 6, December 2010, 649-654, par l'éditeur Taylor & Francis, décrit également une méthode permettant, à partir d'images de radiologie du type desdites première et deuxième images générées 5 dans des plans distincts non parallèles, d'obtenir un modèle virtuel représentatif d'une partie de squelette et de l'enveloppe externe. Ce type de modèle obtenu conformément à la méthode décrite dans cette dernière publication est 10 utilisable pour la mise en oeuvre de la présente invention. Après avoir déterminé : - la représentation géométrique en trois 15 dimensions de ladite partie du squelette de l'individu; et - la représentation géométrique en trois dimensions de l'enveloppe externe E ; Le logiciel estime automatiquement et pour chaque 20 rayon X Rx donné : - une valeur d'épaisseur de tissus mous Li, L2 traversés par ce rayon X donné Rx ; - une valeur d'épaisseur de tissus osseux LO traversés par ce rayon X donné Rx ; et 25 - une valeur d'atténuation globale de ce rayon X donné, entre sa source et sa projection sur la surface de détection où il est projeté. La valeur d'atténuation d'un rayon X donné pour 30 un tissu traversé d'épaisseur x se calcule à partir de l'intensité Io du rayon X (aussi appelé faisceau incident) en sortie de sa source Slat ou Spa et de l'intensité I de ce même rayon X en sortie du tissu traversé, selon la formule I - I e , g étant le 35 coefficient d'atténuation du tissu 3035785 20 La valeur d'atténuation globale d'un rayon X donné traversant plusieurs tissus d'épaisseur xi et de coefficient d'atténuation lui se calcule à partir de l'intensité Io du rayon X (aussi appelé faisceau 5 incident) en sortie de sa source Slat ou Spa et de l'intensité I de ce même rayon X telle que détectée sur un point unitaire ou pixel de l'image correspondante, selon la formule I = Io e-Egua Le logiciel exécuté par le calculateur, calcule à 10 l'aide de la valeur d'épaisseur de tissus mous traversés par ce rayon X donné : - d'une part une valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus mous qu'il a traversée ; et 15 - d'autre part une valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus osseux qu'il a traversée. A l'aide de la valeur d'épaisseur de tissus osseux traversée par ce rayon X donné, de la valeur 20 d'atténuation de ce rayon X donné par la seule épaisseur de tissus mous qu'il a traversée, et de la valeur d'atténuation globale de ce rayon X donné, le logiciel estime une valeur représentative de la densité minérale osseuse des tissus osseux V traversés 25 par ce rayon X donné. Il va maintenant être explicité en référence aux figures 1, la et 2, un exemple d'estimation de répartition de DMO pour une vertèbre V. 30 Sur la figure la, les tissus mous, c'est-à-dire les organes autres que les tissus durs du squelette, ont un coefficient d'atténuation de rayons X qui est sensiblement homogène. Au contraire, dans l'exemple de la figure 2, il 35 existe une hétérogénéité de tissus mous qui présente 3035785 21 de coefficient d'atténuation différent du coefficient d'atténuation des autres tissus mous. Cet exemple permet d'expliquer comment le dispositif selon l'invention qui comprend un calculateur et le logiciel 5 décrit ci-avant met en oeuvre le procédé selon l'invention de matière automatisée. Rappel de l'équation d'atténuation des rayons X : Soit Io l'intensité du faisceau incident émis par la source, I l'intensité du faisceau / rayon telle que 10 détectée sur un pixel de l'image, g le coefficient d'atténuation de la matière traversée par le faisceau, et x l'épaisseur traversée de cette matière. / = 4e-gx 1 Si l'on a plusieurs milieux gi de plusieurs 15 épaisseurs différentes xi, alors l'équation d'atténuation est : e-E mixt 2 Exemple d'utilisation des formules 1 et 2 pour 20 l'estimation de la densité minérale osseuse d'une vertèbre (figure 1) : Soit une vertèbre V donnée avec deux sources biplanes permettant d'avoir des images RX postéroantérieure PA et latérale Lat, et un rayon 25 incident passant par cette vertèbre sur la vue Lat, et soient : Li et L2 les épaisseurs de tissus mous traversées de part et d'autre de la vertèbre V, - gmle coefficient d'atténuation des tissus 30 mous (qui peut être différencié lui-même en plusieurs coefficients, par exemple gm1et1-42) - ktvi les différents coefficients d'atténuation dans la vertèbre (milieu hétérogène) au niveau de la ligne définissant la trajectoire du rayon 3035785 22 incident dans la vertèbre, avec )4 les épaisseurs traversées dans cette ligne Nous avons alors : - = e AML2-E 10 Et donc = e -11M(Ll+L2)-E liViXVi Io Et donc l0e-gm(1+1,2) = e Et donc nous pouvons reconstruire l'image « DMO surfacique», qui ne comporte que le signal os, et qui est obtenue pour chaque rayon incident (et donc chaque pixel) en considérant le terme I corrigé (Icorr) I cOrr l0e-gm(L1+L2) Avec cette modification du signal, nous pouvons considérer la démarche connue d'absorptiométrie par rayon 20 X à simple énergie (Single-energy X-ray absorptiometry SXA), décrite par exemple dans la publication de Mabilleau G et al pour un os isolé (Mabilleau G, Mieczkowska A, Libouban H, Simon Y, Audran M, Chappard D. Comparison between quantitative X-ray imaging, . dual 25 energy X-ray absorptiometry and microCT in the assessment of bone mineral density in disuse-induced bone loss. J Musculoskelet Neuronal Interact. 2015 Mar;15(1):42-52). Une calibration peut être faite avec un fantôme de calibration (éprouvette), généralement en lucite (dont 30 le coefficient d'atténuation est proche des tissus mous) et aluminium (dont le coefficient d'atténuation est proche de l'os), puis le signal est redistribué entre 5 10 15 3 4 5 3035785 23 cette plage de valeurs pour rehausser le signal os, ce qui permet au final d'avoir la DMO surfacique, et ceci avec une image en simple énergie. Remarques : 5 1) Calibration : Pour la calibration (figure la ou 2), nous pouvons prendre avantage du fait que nous avons des zones qui ne sont constituées que de tissus mous et pour lesquelles nous connaissons la longueur L d'épaisseur traversée, donc nous pouvons avoir 10 directement le coefficient d'atténuation des tissus mous du sujet considéré 2) Rehaussement de contraste « os » : cette image « os », lorsque la calibration a permis d'éliminer le signal tissus mous et de redistribuer le signal sur la 15 plage entre tissus mous et os, permet de rehausser le contraste et peut donc être utile à la reconstruction 3D 3) Présence d'hétérogénéités : Du fait des variations interindividuelles et de spécificités au moment de la prise de radiographies, nous pouvons avoir 20 des inhomogénéités locales de coefficient d'atténuation, qui perturbent le signal. Si nous avons des radiographies biplanes calibrées (figure 2), il est possible d'estimer sur les radiographies (postéoantérieure et latérale) par exemple, la dimension de l'élément hétérogène (noté avec 25 l'indice h), et corriger son effet de la manière suivante : o 1/ Estimation de Lpah et Llath (à partir des radiographies biplanes calibrées et d'un modèle géométrique de l'hétérogénéité) 30 o 2/ sur la vue postéroantérieure, pour les rayons Rpa1 qui traversent les tissus mous sans hétérogénéité, connaissant les épaisseurs traversées, on peut en déduire g, (coefficient d'atténuation des tissus mous) 3035785 24 o 3/ pour un rayon Rpa2 qui traverse l'hétérogénéité sur une longueur totale Lpa2 dont Llath du milieu h, le signal sur l'image résulte dellm(4a-Lie0 etilhlieh, ce qui permet d'en déduire p, 5 o 4/ pour produire l'image de DMO en vue latérale, nous connaissons l'épaisseur globale traversée « hors os » (L1 + L2), l'épaisseur traversée par l'hétérogénéité Lpah (grâce à la vue PA), les coefficients 10 d'atténuation itm et ith, et donc en reprenant l'équation 2 et le raisonnement déployé dans les équations 3 à 6 nous définie par le terme : arrivons à l'image DMO 15 1 I corr loe-iim(L1+L2-Lpah)-gh(Lpah) La valeur de Icorr permet d'obtenir une valeur dite de DMO surfacique, telle qu'elle est obtenue en général à partir des modalités d'imagerie de référence (DXA correspondant à la méthode d'imagerie par 20 absorptiométrie biphotonique à rayons X), cette valeur de DMO surfacique ne prend pas en compte l'hétérogénéité du tissu vertébral traversé. Il est également possible, à partir de cette image de DMO surfacique, d'estimer une distribution de 25 densité volumique selon la méthode décrite par Travert (thèse de sciences disponible sur https://tel.archivesouvertes.fr/file/index/docid/834740/filename/TRAVERT. pdf , dont le titre est « Estimation du risque de fracture ostéoporotique du rachis thoraco-lombaire par un 30 modèle en élément finis personnalisé » soutenue à l'ENSAM en 2012 : Brièvement, il s'agit de positionner une vertèbre comportant une distribution générique de densité minérale osseuse, et d'ajuster cette distribution à partir de la DMO surfacique données sur l'une des vues, 3035785 25 ou de manière itérative à partir des données de DMO surfacique obtenues sur les deux vues. On peut ainsi, à l'aide du logiciel / programme exécuté par un calculateur et analysant les première et 5 deuxième images générées par projections de rayons X, estimer des positions de croisement de rayons dans l'espace, chaque position de croisement caractérisant le croisement entre un rayon donné provenant de la première source et un autre rayon donné provenant de la deuxième 10 source. On peut associer des valeurs de densité minérale osseuse volumique à certains au moins des points de croisement. REVENDICATIONS1. Procédé d'estimation de la répartition de la densité minérale osseuse dans au moins une partie de squelette d'un individu, ce squelette étant composé de tissus osseux, ces tissus osseux étant entourés de tissus mous eux même entourés par une enveloppe externe (E) de l'individu, le procédé étant caractérisé en ce qu'il comprend la génération : - d'une représentation géométrique en trois dimensions de ladite partie du squelette de l'individu; et - d'une représentation géométrique en trois dimensions d'au moins une portion de l'enveloppe externe (E) contenant ladite partie du squelette de l'individu ; - d'au moins une première image radiologique de ladite au moins une partie du squelette, cette première image étant générée par observation d'une projection de rayons X, provenant d'une première source, sur une première surface de détection, cette première image étant générée par projection de rayons X et représentant l'individu observé suivant une première incidence d'observation ; - d'au moins une deuxième image radiologique de ladite au moins une partie de squelette, cette deuxième image étant générée par projection de rayons X et représentant l'individu observé suivant une deuxième incidence d'observation distincte de la première incidence d'observation; puis à l'aide d'un logiciel prédéterminé, et pour plusieurs rayons X donnés parmi lesdits rayons X projetés, on estime pour chaque rayon X donné : - une valeur d'épaisseur de tissus mous traversés par ce rayon X donné ; - une valeur d'épaisseur de tissus osseux 3035785 27 traversés par ce rayon X donné ; et - une valeur d'atténuation globale de ce rayon X donné, entre sa sourcé et sa projection sur la surface de détection où il est projeté ; et 5 à l'aide de la valeur d'épaisseur de tissus mous traversés par ce rayon X donné, et d'une estimée d'un coefficient d'atténuation des tissus mous traversés par ce rayon X donné, on déduit d'une part une valeur d'atténuation de ce rayon X donné par 10 l'épaisseur de tissus mous qu'il a traversée et d'autre part une valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus osseux qu'il a traversée ; puis à l'aide de la valeur d'épaisseur de tissus 15 osseux traversée par ce rayon X donné, de la valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus osseux qu'il a traversée et de la valeur d'atténuation globale de ce rayon X donné, le logiciel estime une valeur représentative de la densité 20 minérale osseuse des tissus osseux traversés par ce rayon X donné. 2. Procédé d'estimation selon la revendication 1, dans lequel ladite deuxième image radiologique de ladite 25 au moins une partie de squelette est générée par observation d'une projection de rayons X, provenant d'une deuxième source distincte et éloignée de la première source, sur une deuxième surface de détection. 30 3. Procédé d'estimation selon la revendication 2, dans lequel les première et deuxième surfaces de détection sont respectivement planes et sont orientées entre elles d'un angle de 90° à plus ou moins 60° et préférentiellement à plus ou moins 5° d'angle. 35 3035785 28 4. Procédé d'estimation selon l'une quelconque des revendications 1 à 3, dans lequel pour déduire ladite valeur d'atténuation dudit rayon X donné, par l'épaisseur de tissus mous qu'il a traversée, le logiciel estime un 5 coefficient d'atténuation des tissus mous de l'individu. 5. Procédé d'estimation selon la revendication 4, dans lequel pour estimer le coefficient d'atténuation des tissus mous de l'individu, le logiciel recherche un rayon 10 X particulier traversant uniquement des tissus mous de l'individu sans traverser de tissus osseux et à partir de la position de ce rayon X particulier et de la représentation géométrique en trois dimensions de la portion de l'enveloppe externe (E), le logiciel estime la 15 valeur d'épaisseur de tissus mous traversés par ce rayon X particulier, et par observation d'un point de projection de ce rayon X particulier, le logiciel détermine une valeur d'atténuation subie par ce rayon X particulier lors de sa traversée des seuls tissus mous de 20 l'individu, puis connaissant ladite valeur d'épaisseur de tissus mous traversés par ce rayon X particulier et ladite valeur d'atténuation subie par ce rayon X particulier, le logiciel calcule le coefficient d'atténuation associé aux seuls tissus mous de 25 l'individu. 6. Procédé d'estimation selon la revendication 5, dans lequel connaissant : - le coefficient d'atténuation des tissus mous de 30 l'individu ; - l'épaisseur de tissus mous traversée par ce rayon X donné ; - une intensité (I0) de ce rayon donné à sa source ; et 35 - une intensité (I) de ce rayon donné telle que 3035785 29 détectée sur la surface de détection (P1, P2) sur laquelle il est projeté, le logiciel calcule : - d'une part une valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus mous qu'il a traversée; 5 et - d'autre part, une valeur d'atténuation globale de ce rayon X donné, entre sa source et la surface de détection sur laquelle il est projeté ; puis - en corrigeant cette valeur d'atténuation 10 globale de ce rayon X donné à l'aide de ladite valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus mous qu'il a traversée, le logiciel calcule la valeur d'atténuation de ce rayon X donné par l'épaisseur de tissus osseux qu'il a traversée. 15 7. Procédé selon l'une quelconque des revendications 1 à 6, dans lequel ladite valeur représentative de la densité minérale osseuse des tissus osseux traversés par ce rayon X donné est un coefficient 20 d'atténuation du tissu osseux traversé par ce rayon x donné. 8. Procédé selon l'une quelconque des revendications 1 à 7, dans lequel lesdites première et 25 deuxième images correspondent respectivement à un même instant de détection des rayons X projetés sur les première et deuxième surfaces de détection. 9. Procédé selon l'une quelconque des 30 revendications 1 à 7, dans lequel lesdites première et deuxième images correspondent respectivement à deux instants distincts de détection des rayons X projetés sur les première et deuxième surfaces de détection. 35 10. Procédé selon l'une quelconque des 3035785 30 revendications 1 à 9, dans lequel les première et deuxième images illustrent des projections de l'individu en position debout. 11. Procédé selon l'une quelconque des revendications 1 à 10, dans lequel ladite représentation géométrique en trois dimensions de ladite au moins une portion de l'enveloppe externe contenant ladite partie du squelette de l'individu est déduite à partir desdites 10 première et deuxième images radiologiques. 12. Procédé selon l'une quelconque des revendications 1 à 11, dans lequel ladite représentation géométrique en trois dimensions de ladite au moins une 15 portion de l'enveloppe externe contenant ladite partie du squelette de l'individu est déduite à l'aide d'une observation de l'enveloppe réalisée à l'aide de moyens de capture optique. 20 13. Dispositif comportant un calculateur et un logiciel adaptés à réaliser des calculs ou des estimations nécessaires à la mise en oeuvre du procédé selon l'une quelconque des revendications 1 à 12.