FR 2801776 A1 20010608 FR 9915275 A 19991203 PROCEDE D'UTILISATION D'UN SYSTEME OSTEODENSITOMETRIE, PAR RAYONNEMENT X BI-ENERGIE, A FAISCEAU CONIQUE DESCRIPTION Domaine technique La présente invention concerne un procedé d'utilisation d'un système d'ostéodensitométrie par rayons X. On rappelle que l'os téodensitométrie par rayons X une technique de mesure de masses et de densités osseuses à partir d'acquisitions radiographiques effectuées à une pluralité d'énergies. On utilise en général deux énergies que 1 appelle respectivement "haute énergie" et "basse énergie". Diverses configurations d'une chaîne de mesure (source et détecteur de rayons X) sont envisageables pour effectuer les mesures. On distingue trois familles de systèmes selon types de rayonnements utilisés - les systèmes à pinceau de rayonnement ("pencil beam systems") qui utilisent une source rayons X collimatée par un trou et un monodétecteur de rayons X qui est également collimaté, - les systèmes à faisceau en éventail ("fan beam systems") qui utilisent une source de rayons X collimatée par une fente et un détecteur linéaire de rayons X, et - les systèmes à faisceau conique ("cone beam systems") qui utilisent détecteur bidimensionnel de rayons X. L'invention concerne plus particulierement la mise en oeuvre de systèmes d'ostéodensitométrie par rayonnement X biénergie à faisceau conique. Etat de la technique antérieure Les principes méthodologiques de l'ostéodensitométrie par rayonnement X biénergie et les principales solutions techniques actuellement utilisées sont connus par les deux documents suivants auxquels on se reportera [1] 'Technical Principles of Dual Energy X-Ray Absorptiometry", G.M. Blake et I. Fogelman, Seminars in Nuclear Medicine, Vol XXVII, n 3, juillet 1997, pages 210 à 2 et [2] The Evaluation of Osteoporosis : Dual Energy X-Ray Absorptiometry and Ultrasound in Clinical Practice", Second Edition, G.M. Blake, H.W. Wahner et I. Fogelman, Martin Dunitz Editor, 1999, ISBN 1-85317-472-6. On se reportera plus particulièrement aux chapitres 3, 4 et 5 du document [2] où sont décrits les principes de mesure de densités osseuses à deux énergies et les systèmes connus pour faire de telles mesures. Pour réaliser le diagnostic et suivi thérapeutique de l'ostéoporose en des sites pertinents (colonne vertébrale, hanches, avant-bras corps entier), les systèmes commercialement disponibles sont type "pencil beam" et "fan beam". Ces systèmes utilisent des collimateurs qui limitent les rayonnements parasites, en particulier les rayonnements diffusés ("scattered"). Par ailleurs, la géométrie de systèmes d'acquisition limite les déformations géométriques liées à la dimension du détecteur, qui resulteraient, en particulier, de la conicité du faisceau de rayons X. De plus, ces systèmes commercialement disponibles utilisent un balayage en effectuant successivement, pour chaque position de la chaîne d'acquisition, une mesure à haute énergie et une mesure à basse énergie. Cela assure une parfaite cohérence entre les zones du patient qui sont observées lors des deux mesures. Les systèmes à faisceau conique commercialement disponibles ne sont prévus que pour faire des mesures densité osseuse sur des zones périphériques telles que les phalanges, les mains, les avant-bras et les calcanéums. Par rapport aux deux autres familles de systèmes, la mise en oeuvre de ces systèmes à faisceau conique est plus délicate. En effet, le rayonnement diffusé engendré par 'interaction entre le rayonnement X incident et la partie étudiée du corps est important et doit être traité. De plus, la conicité rend la mesure dépendante la position de la zone anatomique considérée dans le champ image. En outre, durant le temps de lecture du capteur bi-dimensionnel, le patient peut avoir bougé entre la mesure à haute énergie et la mesure à basse énergie. si les svstèmes du genre "cone beam" commercialement disponibles sont limités à l'analyse d'extrémités du corps humain telles que les mains, les avant-bras et les calcanéums, c'est parce que ces zones sont peu épaisses et ont de faibles dimensions. En conséquence, les rayonnements parasites sont limités et le positionnement peut être consolidé, par exemple au moyen d'un système de fixation pour un avant-bras ou un calcanéum. En revanche, la mise en oeuvre de tels systèmes à faisceau conique pour des zones anatomiques de dimensions plus importantes, en particulier la colonne vertébrale ou les hanches, zones essentiellement utilisées pour le diagnostic et le suivi de l'ostéoporose, nécessite de prendre en compte les phénomènes parasites ainsi que l'importante conicité du rayonnement X, et de s'adapter aux mouvements potentiels du patient. On connait déjà des systèmes d'os téodensitométrie de zones bidimensionnelles, comprenant une source apte à fournir des flux de rayons X à au moins deux énergies, un écran-convertisseur apte à transformer les rayons X en photons de lumière visible, un système optique de reprise d'image et une caméra CCD, par les documents suivants auxquels on se reportera [3] brevet US-5,150,394 du 22 septembre 1992 "Dual-Energy System for Quantitative Radiographic Imaging", (Andrew Karellas), et [4] Demande internationale publiée le 14 novembre 1996, n de publication WO 96/35372, "A System for Quantitative Radiographic Imaging", (Andrew Karellas). Exposé de l'invention La présente invention a pour objet un procédé d'utilisation d'un système d'ostéodensitométrie genre de ceux qui sont décrits dans les documents [3] [4] ce procédé permettant d'augmenter la précision la reproductibilité des mesures de densité osseuse dans zones anatomiques d'un patient, à partir de radiographies bidimensionnelles de ces zones anatomiques, acquises à une pluralité d'énergies. Pour ce faire, la présente invention utilise mise en correspondance ("matching") entre une image acquise à haute énergie et une image acquise à basse énergie et en outre, de préférence, un moyen d'aide positionnement du patient. De façon précise, la présente invention a pour et un procédé d'utilisation d'un système d'ostéodensitométrie par rayons X à au moins deux énergies, à faisceau conique, ce système comprenant une source de rayons X apte à fournir un faisceau conique de rayons X à au moins une première énergie, appelée haute énergie, et une deuxième énergie, appelée basse énergie et inférieure à la première énergie, un détecteur bidimensionnel de rayons X et des moyens électroniques de traitement des images fournies par ce détecteur, ce procédé étant caractérisé en ce qu'on acquiert l'image à basse énergie d'une zone anatomique d'un patient, on acquiert l'image à haute énergie cette zone anatomique et l'on met en correspondance images respectivement acquises à basse énergie et à haute énergie, avant de construire la carte des densités osseuses de la zone anatomique. Selon un mode de mise en oeuvre particulier du procédé objet de l'invention, pour mettre en correspondance les images respectivement acquises à haute et basse énergies . - on extrait respectivement de ces images les ensembles de contours des zones osseuses de zone anatomique, - on recherche la transformation plane optimale permettant de mettre l'ensemble de contours relatif à l'image acquise à haute énergie correspondance avec l'ensemble de contours relatif à l'image acquise à basse énergie inversement, selon une variante), - on utilise cette transformation pour amener l'image acquise à haute énergie dans référentiel de l'image acquise à basse énergie (ou inversement, selon la variante ci-dessus) - combine ces deux images pour déterminer carte de mesure des densités osseuses. Selon un mode de mise en oeuvre préféré du procédé objet de 'invention, avant de faire les acquisitions à haute et basse energies, on fait un cliché de radioscopie, c'est '-dire un cliché à faible dose de rayons X, de la zone anatomique du patient pour aider à positionner cette zone dans le système. De préférence, lorsque le patient est soumis à un premier examen, on utilise ce cliché à faible dose de rayons X pour rétroagir sur la mécanique du système afin de positionner la zone anatomique par rapport à un référentiel préétabli. De préférence, lorsque le patient est soumis à un examen de suivi, on utilise le cliché à faible dose de rayons X pour placer la zone anatomique dans une position identique à celle qu'elle occupait lors de l'examen précédent. De façon plus explicite, dans le mode de réalisation préféré avant de faire les acquisitions à haute et basse énergies, on fait un cliché, à faible dose de rayons X à une seule énergie, de la zone anatomique du patient. A partir d'une détection des contours osseux sur ce cliché, on déduit les paramètres géométriques permettant soit de centrer la zone anatomique dans champ d'acquisition en cas de premier examen, soit de positionner la zone anatomique par rapport au champ d'acquisition de la même façon que dans l'examen précédent si ce n'est pas le premier examen. Dans les deux cas, ce positionnement est effectué par déplacement du patient par rapport au système source-detecteur, ou du système source- détecteur par rapport au patient, par commande manuelle ou automatique. On connaît déjà la possibilité d'utiliser cliché préalable à l'examen (également appelé prescan scout data ) pour centrer le patient dans domaine de la tomographie par le brevet US 5457724 Automatic field of view and patient centering determination from prescan scout data du 10 octobre 1995. Dans ce brevet, avant de reconstruire une coupe tomographique d'un patient, on acquiert deux projections monodimensionnelles (éventail ou beam ) à 0 et 90 de cette coupe. On détecte les points correspondant aux bords du patient dans les deux projections, et on en déduit la position du centre de la zone et la taille d'acquisition tomographique. Ces paramètres sont donnés à l'opérateur et peuvent être utilises pour déplacer le patient afin de mieux le centrer pour l'acquisition tomographique. Le but d'obtenir la meilleure qualité d'image possible, systèmes tomographiques étant étudiés pour que l'atténuation maximale se trouve au centre de la zone d'acquisition, et les corrections de durcissement de spectre étant effectuées en fonction de la taille du champ acquisition. Dans la présente invention il s'agit de radiologie (avec un détecteur bidimensionnel) et non pas de tomographie (avec un détecteur fan-beam en rotation). L'image finale est une projection bidimensionnelle et pas une coupe reconstruite. De plus, dans le mode de réalisation préféré on utilise un seul cliché préalable et non pas deux clichés à 90 l'un de l'autre. En outre, un but de l'invention est la reproductibilité de la mesure de masse osseuse calculée à partir de l'image, et pas la qualité de 1 image elle- même. De plus, dans les systèmes tomographiques le recentrage du patient n'est pas automatique (il l'est en hauteur mais pas en largeur, car le déplacement latéral de la table n'est pas prévu n'est pas nécessaire). reproductibilité est la propriété de l'appareil de mesure de donner la même mesure pour différents examens sur le même patient (supposé de densité osseuse constante) et le même site anatomique. Dans le cas d'un patient dont la masse osseuse varie dans le temps, sous l'effet d'une maladie ou d'un traitement par exemple, cette propriété de reproductibilité permet de quantifier ces variations de la masse osseuse. plus, dans la présente invention on peut utiliser ce prescan scout data pour positionner automatiquement des caches permettant de limiter la zone irradiation à la zone osseuse, ce qui n'est pas possible en tomographie sous peine de projections tronquées. De cette façon la dose d'irradiation reçue par le patient est minimisée. Rappelons que la mise en correspondance d'images ( image matching ) est utilisée dans plusieurs domaines du traitement d'images tels que, par exemple, la vision stéréoscopique, l'analyse de séquences d'images et la fusion d'images de modalités ou de natures différentes. Différentes approches ont été développées, par exemple la mise en correspondance entre les nappes de niveaux de gris ou l'extraction de primitives et la mise en correspondance entre les primitives. A ce sujet, on se reportera au document suivant "Overview of Image Matching Techniques", C. Heipke Proceedings of OEEPE Workshop on the Application of Digital Photogrammetric Workstation, Kôlb 0. Editor, OEEEPE Official publications, n 33, pages 173 à 189, 1996. En ce qui concerne le positionnement du patient, sur les systèmes de type "pencil beam" ou " beam", un premier positionnement du patient s'effectue à l'aide d'un pointeur laser qui repère la zone examen à partir d'observations morphologiques externes. Ensuite, le balayage débute. Si le patient bien positionné l'examen se poursuit mais si à l'observation sur l'écran des premières lignes acquises le positionnement n'est pas bon, l'opérateur arrête tout et effectue un nouveau positionnement puis relance l'examen. A ce sujet, on se reportera au document [2] pages 198 ' 200 pour ce qui concerne la colonne vertébrale et aux pages 265 à 267 pour ce qui concerne la hanche. étude récente a montré que, pour des systèmes du genre "pencil beam", le repositionnement avait lieu dans des cas et que, pour environ 10% des examens, il fallait repositionner jusqu'à trois fois le patient. ce sujet, on se reportera au document suivant [6] Insights, vol. 10, n 1, mars 1999, pages 10 et (revue éditée par la Société Hologic), 'Independent survey reveals surprising, disappointing results for Lunar users". Sur les systèmes du genre "cone beam" connus, utilisés pour des examens de zones périphériques, le positionnement du patient est assuré par un système mécanique d'aide au positionnement, par exemple une poignée pour l'avant-bras et une cuvette formée pour le talon. Brève description des dessins La présente invention sera mieux comprise à la lecture de la description d'exemples de réalisation donnés ci-après, à titre purement indicatif et nullement limitatif, en faisant référence dessins annexés sur lesquels - la figure 1 est une vue schématique d'un système d'ostéodensitométrie osseuse par rayonnement X à deux énergies, à faisceau conique, qui est utilisable pour la mise en oeuvre de l'invention, - les figures 2 et 3 sont des organigrammes de procédures qui sont utilisées dans des modes de mise en #uvre particuliers de l'invention, - la figure 4 représente le diagramme fonctionnel d'un mode de réalisation avantageux d' procédé de correction des défauts d'une image qui est utilisable dans la présente invention, et - la figure 5 représente le diagramme fonctionnel d'une variante de ce mode de réalisation avantageux du procédé de correction des défauts d'une image, qui est également utilisable dans la présente invention. Exposé détaillé de modes de réalisation particuliers On voit sur la figure 1 un système d'os téodensitométrie osseuse 1 qui comprend une source de rayons X 1a apte à envoyer un faisceau conique 1b de rayons X vers le corps d'un patient 1c en cours d'examen. Cette source 1a est apte à émettre rayonnements X correspondant respectivement à deux niveaux distincts d'énergie. Ces deux niveaux sont utilisés pour obtenir deux images distinctes patient. Un filtre amovible 1d est interposable entre la source 1a et le patient 1c et sert à améliorer les qualités spectrales du faisceau. système 1 comprend aussi un détecteur bidimensionnel 2 qui est très schématiquement représenté en coupe transversale sur la figure 1 et destine à détecter les rayons X émis par la source et ayant traversé le patient 1c. ce détecteur 2 est parallèle à un plan défini par deux directions orthogonales x et y. Ce patient est placé sur un support approprié 2a, par exemple lit, qui est transparent aux rayons X. Dans l'exemple de la figure 1 la source 1a (munie de l'éventuel filtre 1d) est placée au-dessus du patient reposant sur support tandis que le détecteur est placé en-dessous de ce support. Des moyens non représentés sont prévus pour déplacer le support 2a par rapport à la source 1a et au détecteur 2, qui sont alors fixes, ou sont prévus pour déplacer la source 1a et le détecteur 2 par rapport au support 2a qui est alors fixe, ces déplacements ayant lieu parallèlement aux directions x et y. Dans l'invention on peut utiliser tout type de détecteur bidimensionnel, par exemple un écran sensible aux rayons X apte à fournir directement un signal électronique représentatif de l'image acquise par le détecteur sous forme de pixels. Au lieu cela on peut utiliser un écran scintillateur prévu pour recevoir les rayons X ayant traversé le patient et pour convertir ces rayons X en lumière visible. Cette dernière est alors envoyée, par l'intermédiaire d'un miroir, à un capteur CCD muni d'un objectif comprenant un réseau de pixels photosensibles. Sur la figure 1, on voit aussi un dispositif 3 du genre controleur de CCD ou analogue qui lit, pixel par pixel, la représentation d'image fournie par le détecteur et qui numérise cette représentation. La représentation ainsi numérisée est stockée dans une mémoire 3a. Un ordinateur 3b est prévu pour traiter images ainsi mémorisées. Un dispositif d'affichage 3c, comprenant par exemple un tube à rayons cathodiques, est prevu pour afficher les images avant ou après ce traitement. Un tel système est utilisable pour mettre en oeuvre un procédé conforme à l'invention, selon lequel, en vue d'obtenir une bonne précision et bonne reproductibilité de la mesure de densité osseuse - on réalise une mise en correspondance entre les images respectivement acquises à haute energie et à basse énergie et - préférence, on utilise en outre un premier cliché, à faible dose, de type scopie, pour aider à positionner le patient dans le système ostéodensitométrie. Considérons d'abord la mise en correspondance entre l'acquisition à haute énergie et l'acquisition à basse énergie. Afin de pouvoir construire une carte de densité de l'os, il est indispensable que, pour chacune de ces acquisitions, les rayons X aient eu le même parcours dans la zone anatomique. Pour un système du genre "cone beam", le patient est susceptible de bouger entre l'acquisition à haute énergie et l'acquisition à basse énergie. Une façon d'éviter cela est de réaliser ces deux acquisitions dans un court intervalle de temps, inférieur à 200 ms, comme cela se fait sur les systèmes du genre "pencil beam" et "fan beam". Mais, à la différence de ces systèmes, pour le système du genre "cone beam", il faut traiter des images bidimensionnelles c'est-à-dire des images comportant une grande quantité de données, nécessitant un temps de lecture important. Cela nécessite des technologies très sophistiquées et très onéreuses. L'invention permet de ne pas être soumis à cette contrainte de rapidité entre les deux acquisitions, traitant a posteriori les acquisitions à haute et basse énergies. si les deux acquisitions sont réalisées quelques secondes, par exemple environ 10 secondes, la zone anatomique est susceptible d'avoir subi un léger déplacement. En utilisant les acquisitions à haute et basse énergies telles quelles, on induit une erreur dans la mesure de densité osseuse. Pour remédier à cet inconvénient, on propose de mettre en correspondance les acquisitions à haute et basse énergies avant de construire la carte des densités osseuses. Considérons ensuite l'utilisation d'un cliché de scopie pour le positionnement du patient. Etant donné que, pour un système du genre "cone beam", on dispose d'un capteur bidimensionnel qui permet en une seule acquisition, d'avoir une vision globale de la zone analysée, on propose de réaliser un cliché à faible dose (cliché de scopie) avant les acquisitions à haute et basse énergies pour aider au positionnement du patient. ce patient est soumis à son premier examen, on utilise ce cliché de scopie pour rétroagir sur la mécanique du système (c'est-à-dire commander la mécanique du dispositif d'acquisition d'images de façon qu'il se positionne correctement par rapport au patient ou inversement) afin de positionner la zone anatomique rapport à un référentiel préétabli. Si l'examen est un examen de suivi du patient, on utilise le cliché de scopie pour mettre la zone anatomique dans une position identique à celle qu'elle occupait lors du précédent examen (lors duquel on avait fait également un cliché de scopie). Ce type de démarche permet d'augmenter sensiblement la reproductibilité des mesures effectuées avec des systèmes de type "cone beam". En raison de la conicité du faisceau, la mesure depend de la position de la zone anatomique dans ce faisceau. En utilisant le cliché de scopie pour positionner de façon identique la zone anatomique du patient par rapport à un référentiel donné ou pour assurer la cohérence entre deux examens, obtient une bonne reproductibilité de l'examen. A titre d'exemple, pour la mise en correspondance entre les acquisitions à haute et basse énergies, procède de la façon suivante, illustrée par 1 organigramme de la figure 2 1.1. on effectue ces acquisitions à haute et basse énergies (étapes E1 et E2) ; 1.2. on extrait respectivement de ces deux acquisitions les deux ensembles (n 1 et n 2) de contours des zones osseuses (étapes E3 et E4) ; 1.3. on cherche la meilleure (au sens des moindres carrés par exemple) transformation géométrique plane (composition d'homothétie, translation et rotation) permettant de mettre en correspondance les deux ensembles de contours (étape E5) ; 1.4. utilise cette transformation pour amener l'image acquise à haute énergie dans le référentiel de l'image acquise à basse énergie (étapes E6 et E7) ; 1.5. combine ces deux images pour déterminer la carte de mesure des densités osseuses (étape E8). En ce qui concerne l'utilisation d' cliché de radioscopie (en abrégé scopie) pour le positionnement, un exemple de mise en oeuvre pour un patient soumis à son premier examen consiste à 2.1. réaliser un cliché à faible dose ; 2.2. extraire de cette acquisition les contours des zones osseuses ; 2.3. repérer des points caractéristiques dans la carte de contour, par exemple identification de vertèbres ou repérage points caractéristiques sur le col du fémur ; 2.4. construire la fonction de type translation qui permet de placer au mieux ces points par rapport à une position standard définie au préalable ; 2.5. rétroagir sur la mécanique de positionnement pour se ramener à la position standard ; 2.6. réaliser les acquisitions à haute et basse énergies. Dans le cas d'un patient soumis à examen de suivi thérapeutique, l'étape de construction de fonction 2.4 est remplacée par les deux étapes suivantes 2.4.1 récupérer les positions des points caractéristiques dans les acquisitions d'un examen précédent du patient ; .4.2 construire la fonction de type translation qui permet de mettre au mieux en correspondance les points caractéristiques du cliché de scopie par rapport à leur position dans un examen antérieur. Tout ceci est précisé par l'organigramme de la figure 3 . - Étape F1 : on positionne grossièrement le patient pour observer la zone anatomique d'intéret (voir 2.1) - Étape F2 : on extrait les contours de structure osseuse (voir 2.2) - Étape F3 : on identifie les points caractéristiques (voir 2.3) - Étape F4 : on se demande s'il s'agit premier examen - Si ' on va à l'étape F5 dans laquelle on identifie la transformation géométrique amenant les points caractéristiques en position standard (voir 2.4) puis on à l'étape F6 dans laquelle on applique un mouvement à la mécanique pour amener ces points vers la position standard (voir 2.5) - Si non on va à l'étape F7 dans laquelle on récupère la position des points caractéristiques dans l'examen antérieur (voir 2.4.1) puis on va à l'étape F8 dans laquelle on identifie la transformation géométrique amenant les points caractéristiques courants vers ceux de l'examen antérieur (voir 2.4.2) puis on va a l'étape F9 dans laquelle on applique un mouvement à la mécanique pour amener le patient vers une position correspondant à un bon placement des points caractéristiques (voir 2.5). On donne dans ce qui suit des précisions sur traitement des images selon un mode de réalisation préféré du procédé objet de l'invention. Les étapes a) à e) qui suivent constituent étapes préparatoires hors ligne . a) La caractérisation du décalage ("offset") en chacun des pixels qui composent une image est réalisée par acquisition d'images de noir (images sans rayonnement X), notées NOIR, et la caractérisation en gain s'obtient par 1 acquisition d'images, à plein flux à chaque énergie de fonctionnement c'est-à-dire avec le flux maximum de rayonnement X, sans objet ou avec un objet d'épaisseur uniforme, notées PFL. Afin d'obtenir des caractéristiques peu affectées par le bruit d'acquisition, realise un nombre important d'acquisitions dont fait la moyenne. b) La carte des pixels en défaut est obtenue l'analyse d'une acquisition à plein flux d'un objet d'épaisseur uniforme. c) La caractérisation des distorsions géométriques est obtenue par l'acquisition d'images grilles avec des motifs répartis régulièrement ; à partir de la description geométrique de ces grilles et de leurs images construit la fonction de distorsion géométrique du système. Afin de separer cet effet de celui de la conicité, il faut placer les grilles au contact du détecteur (référence 2 dans l'exemple de la figure 1). d) Le calibrage géométrique du système réalisé l'acquisition d'une grille de géométrie connue que l'on place à des positions parfaitement identifiées par rapport au detecteur. Cela permet de remonter à une position de la source de rayons X. e) L'étalonnage "en densité" est le même pour tous les pixels de détection (pixels obtenus avec le détecteur 2 dans l'exemple de la figure 1), les différences de gain ayant été corrigées préalablement. On utilise un fantôme composé d'un certain nombre d'éléments, chaque élément étant composé d'un mélange de deux matériaux de base, d'épaisseurs connues. Les matériaux de base peuvent être un matériau n 1 simulant l'os et un matériau n 2 simulant les tissus mous. On s'assure mécaniquement que les mesures à haute et basse énergie sont acquises perpendiculairement à chaque élément, sans avoir à effectuer de correction de conicité. Un dispositif de collimation (non représenté sur la figure 1), que l'on dispose par exemple entre la source de rayonnement et le fantôme ou entre le fantôme et le détecteur, garantit que ces mesures sont exemptes de rayonnement diffusé. A partir de ces mesures d'étalonnage, on estime les fonctions qui transforment un couple (mesure à basse énergie, mesure à haute énergie) en un couple (épaisseur équivalente de matériau de base n 1, épaisseur équivalente du matériau de base n 2). À ces étapes préparatoires succèdent des étapes de traitement de l'image. On calcule une image corrigée par la formule suivante image corrigée=(Acquis-NOIR)/(PFL-NOIR) ' "Acquis" représente successivement l'acquisition à basse énergie et l'acquisition à haute énergie. Et les pixels en défaut isolés sont corrigés par calcul de moyenne locale ou par un calcul de médiane. Pour les paquets de pixels on utilise une approche fondée sur un procédé pour corriger les défauts d'une image provenant d'un détecteur rayons X (ou y) de type matriciel (définition d'un degré de confiance aux différents pixels et mise en oeuvre d'une convolution normalisée), ce procédé (invention de R. Guillemaud - demande de brevet français n EN 15044 du 30 novembre 1998) étant expliqué à la fin de la présente description. A partir de la fonction de distorsion identifiée dans l'une des étapes préparatoires "hors ligne", c est-à-dire l'étape c), on corrige la position des pixels observés et, par interpolation, on reforme une image échantillonnée sur une grille à pas ("pitch") régulier. La mise en correspondance entre l'acquisition à basse énergie et l'acquisition à haute énergie est effectuée par extraction d'un ensemble de contours dans chacune des deux acquisitions (contours liés à l'interface tissus/os) et par la détermination des parametres de rotation et de translation permettant de passer d'un ensemble de contours à l'autre. Pour le rayonnement diffusé, on utilise le procédé de correction du flux diffusé dans des images de radiographie numérique qui est décrit dans demande internationale publiée le 20 août 1998, n de publication WO 98/36380 (M. Darboux et J.M. Dinten), et qui utilise un modèle tridimensionnel zones anatomiques, en décrivant celles-ci comme des zones composées de deux matériaux : les tissus mous (valeur correspondant à une proportion de graisse et de tissus maigre chez un individu standard) et les os. Les dimensions de ce modèle tridimensionnel, qui sont parallèles au plan de projection, sont données par l'analyse des projections. Les cotes en profondeur résultent d'une estimation de l'épaisseur la zone anatomique (dans une zone sans os) et de données a priori fournies par les traités d'anatomies (par exemple, position de la colonne vertébrale rapport au pour une examen en posterio-anterior). Ce procédé connu consiste ensuite en une approximation des équations de Klein et Nishina en utilisant cette description a priori et conduit à l'expression du rayonnement diffusé comme une fonction non linéaire du rayonnement direct. Pour plus de détails, on se reportera à cette demande internationale WO 98/36380. Pour la détermination de l'image d'os, on applique fonctions qui ont été calculées lors de l'étape "hors ligne" afin de calculer une image " et une image "tissus mous". Les tissus mous, avec une proportion variable de graisse, apparaissent sur les deux images à la fois et biaisent l'image 'os. Une segmentation permet de déterminer les zones ne contenant pas d'os ; sur ces zones on peut calculer une ligne de base" qui correspond à la quantité tissus gras apparaissant comme de l'os. Cette ligne de base peut être modélisée et extrapolée dans zones ne contenant pas d'os puis soustraite de l'image d'os. Pour la détermination de la densité osseuse, la conicité induit un biais sur ces mesures, conduisant à relation entre reproductibilité et position de la zone osseuse dans l'acquisition. Afin de limiter cette sensibilité, on corrige l'effet de conicite sur l'aire mesurée en utilisant les paramètres géométriques du système d'acquisition déterminés dans l'étape "hors ligne". La mesure de densité osseuse n'est, quant à elle, pas affectée (ou très faiblement affectée en raison de la discrétisation) par l'effet de conicité. On explique maintenant le procédé de correction défauts d'une image provenant d'un détecteur de rayons X (ou y) de type matriciel. Pour éviter toute confusion avec la présente invention, ce procédé est appelé "procédé PCD" dans ce suit. Ce procédé PCD est un procédé de correction des défauts d'images X ou y, fondé sur la construction d'une carte de confiance regroupant les indices de confiance (continus, compris entre 0 et 1) de tous les pixels du détecteur. De façon plus précise, ce procède PCD est un procédé de correction des défauts d'images (relatifs à des pixels en défaut ou à des pixels non-mesurés par le détecteur) issus d'un détecteur de rayons X ou y de type matriciel, consistant en une étape calibrage ou calibration du détecteur, puis en étape de correction des pixels en défaut de l'image à corriger. Le procédé PCD se caractérise par le fait que l'étape de calibration consiste en les opérations suivantes a) - acquisition d'une image noire ; b) - acquisition d'au moins une image d'un objet connu ; c) - détermination d'au moins une image d'offset et de gain moyen ; d) - acquisition d'au moins une image l'objet à mesurer avec une atténuation uniforme lorsque le détecteur est de type à rayons X et une émission uniforme lorsqu'il est de type à rayons Y ; e) - détermination d'une première moyenne et d'un premier écart-type pour toute l'image de l'objet à mesurer ; f) - réalisation d'une première carte de confiance en fonction du niveau de gris des pixels de l'image ainsi que de l'écart-type et de la moyenne déterminés en e) ; et en que l'étape de correction consiste en l'acquisition d'une image de l'objet à mesurer et en une correction de l'image de l'objet à partir de la carte de confiance et de la dernière image de l'objet acquise. Selon une variante de l'invention, l'etape de calibration consiste, après l'opération f) de réalisation d'une première carte de confiance en les opérations suivantes g) - sur l'image de l'objet connu, détermination gradient de niveau de gris ; h) - détermination d'une seconde moyenne et d'un second écart-type, puis d'une seconde carte de confiance ; et i) - détermination d'une troisième carte de confiance par produit, pixel à pixel, de la première et la seconde cartes de confiance. Avantageusement, la première et la seconde cartes de confiance sont constituées de l'ensemble indices de confiance attribués à chaque pixel de l'image, l'indice de confiance de chaque pixel étant déterminé en fonction du niveau de gris de ce pixel, ainsi que de l'écart-type et de la moyenne de toute l'image de l'objet à mesurer. Selon un mode de réalisation préféré du procédé PCD, l'étape b) d'acquisition d'image d'un objet connu est réalisée avec une atténuation uniforme. Selon une variante du procédé PCD, dans laquelle les défauts d'images résultent de l'absence de mesure sur au moins une zone de l'image, l'étape de calibration consiste, après l'opération f) ou l'opération i), à - agrandir l'image de l'objet à mesurer ; - ajouter, à la place de la zone de l'image sans mesure, des pixels ayant des indices de confiance nuls ; et - déterminer une quatrième carte de confiance, à partir de la première ou de la troisieme carte de confiance, en tenant compte de ces indices de confiance nuls ; l'étape de correction consiste à effectuer une nouvelle acquisition d'image de l'objet ' mesurer, à lui appliquer l'agrandissement et à corriger l'image de l'obj à partir de la quatrième carte de confiance et de la dernière image de l'objet acquise. donne maintenant une description détaillée de modes de réalisation particuliers du procédé PCD. Ce procédé PCD a pour but de corriger les pixels en défaut ou les pixels manquants d'images provenant d'un détecteur de rayons X ou y de type matriciel constitué de pixels. Plus précisément, ce procédé consiste à construire une carte de confiance à partir d'indices de confiance continus, compris entre 0 et 1, et représentant le niveau de confiance attribué à chaque pixel du détecteur. I1 consiste ensuite à utiliser cette carte de confiance pour corriger les défauts du détecteur matriciel en utilisant un filtrage, tel que le filtrage de KNUTSSON, décrit ultérieurement, ou un filtrage par équations aux dérivés partielles (EDP). Selon le procédé PCD, la détection pixels en défaut n'est pas binaire, mais continue, qui permet traitement plus précis au niveau des frontières des zones de pixels en défaut. Le procédé PCD comporte deux étapes principales, à savoir - une étape de calibration du détecteur qui consiste à effectuer l'acquisition d'une ou plusieurs images d'un objet de référence connu ; puis à partir de images, à détecter les pixels en défaut, puis à déterminer la carte de confiance regroupant l'ensemble des indices de confiance de tous les pixels de l'image ; et - une étape de correction qui consiste à partir de la carte de confiance, à corriger l'image de l'objet à mesurer au moyen d'un filtrage, tel que celui proposé par KNUTSSON. D'une façon générale, la détection des pixels en défaut se fait grâce à l'acquisition d'une plusieurs images d'un objet de référence connu qui, de préférence, a une atténuation uniforme. Dans le cas où il n'y a aucun pixel défaut sur le détecteur, l'image obtenue est uniforme, au bruit près crée par la génération des photons X ou -y, la détection et la conversion de ces rayons X ou y en valeurs numériques. Dans le cas où il y a des pixels en défaut sur le détecteur, la détection de ces pixels en défaut se fait à partir d'une analyse statistique de l'image de calibration obtenue et de l'image de gradient correspondante, car un pixel en défaut a soit une valeur aberrante, soit un fort gradient. Lorsque les pixels en défaut sont détectés, on attribue à chacun de pixels un indice de confiance qui varie de façon continue entre 0 et 1. Le cas de détection de défauts définis par un écart à la moyenne d'une image uniforme le cas de détection de défaut par un pixel à fort gradient pour acquisition a priori uniforme peuvent être combinés ; la combinaison de ces deux criteres se fait grâce au produit des deux cartes de confiance effectués chacune sur un critère différent. Dans ce procédé PCD, on fait l'hypothese justifiée la distribution des pixels en défaut est stable dans le temps. Les défauts vont donc être visibles sur toutes les images acquises sur le dispositif. En effet, correction de pixels en défaut est spécifique au détecteur et non à l'image observée. Elle comble les manques d'information en utilisant mieux les informations fournies par les pixels voisins et en assurant une continuité de l'information. Les images ainsi corrigées sont plus facilement exploitables visuellement par un opérateur dans un ectif, notamment de diagnostic. Sur la figure 4, on a représenté schématiquement diagramme fonctionnel des étapes principales du procédé PCD. De façon plus précise, procédé PCD comporte une étape 10 d'acquisition une ou de plusieurs images noires, c'est-à-dire des images effectuées dans le noir, sans émission rayonnements X. Cette acquisition d'images noires permet de mesurer le bruit de l'électronique du dispositif de lecture associé au détecteur. Le procédé PCD consiste ensuite en une etape 11 d'acquisition d'une ou de plusieurs images d' objet de référence connu. Cette acquisition d'image de référence se fait, de préférence, avec une atténuation uniforme lorsqu'il s'agit d'un détecteur à rayons X. L étape suivante du procédé PCD, référencée 12, consiste à calculer des images d'offset et de gain moyens Moy N et Moy_G, avec où Gn représente les images de l'objet de reference connu, Nn représente les images noires, et où t le nombre d'images acquises aussi bien dans l'étape 10 que dans l'étape 11. Le procédé PCD se poursuit par une étape 13 qui consiste à effectuer l'acquisition d'une image Im de l'objet à mesurer, appelée acquisition image nouvelle . Cette acquisition s'effectue pour une atténuation uniforme, dans le cas d'un détecteur de rayons X. Cette image Im est alors corrigée de 1 offset et du gain de la façon suivante où est une constante choisie par l'utilisateur, et qui a uniquement pour but de définir la dynamique de niveau de gris pour coder la valeur des pixels de l'image corrigée. Toutefois, la correction d'offset et de gain peut se faire selon des méthodes différentes de celle qui été décrite précédemment et, par exemple, en utilisant une série d'objets d'atténuations différentes. Une étape 14 consiste ensuite à déterminer l'écart-type (ou la variance) et la moyenne sur l'image corrigée Im cor. Cette moyenne et cette variance sont calculées sur l'image corrigée de l'offset. La moyenne est déterminée par Moy = 1/nombre-pixel lxyyIm-cor(x, y) et variance par Var = l /nombre-pixel 1xzy(Im-cor(x, y) - Moy)2 . Ces valeurs de moyenne et de variance peuvent éventuellement être affinées en supprimant l'effet points aberrants. Pour cela, on recalcule la variance et la moyenne de l'image uniquement pour les pixels pour lesquels les différences de la valeur avec moyenne Moy est inférieure d'un pourcentage x % à la racine carrée de la variance, x étant défini par l'opérateur. procédé PCD se poursuit par une étape 15 consiste à déterminer l'indice de confiance de chaque pixel détecteur. En d'autres termes, un indice confiance compris de façon continue entre 0 et 1 est calcule pour chaque pixel de l'image en fonction de la valeur du niveau de gris du pixel, de la variance et de la moyenne déterminées pour l'image corrigée. Plus précisément, l'indice de confiance se détermine par - le niveau de gris du pixel, déterminé selon la normalisation des niveaux de gris de l'ensemble des pixels de façon à obtenir une distribution normale - la transformation de la valeur par la fonction F de répartition de la loi normale pour obtenir une distance dite normale dist norm(x, y) = F(norm(x, y)) - le calcul de l'indice de confiance du pixel C1(x, y) = 2(1 - dist norm(x, y)). Pour un pixel dont le niveau de gris de départ est proche de la moyenne Moy, l'indice de confiance vaut 1 ; plus la valeur du niveau de gris du pixel s'éloigne de la valeur moyenne Moy, plus l'indice de confiance se rapproche de 0. Le procédé PCD se poursuit par un test 16 qui consiste ' vérifier si tous les indices de confiance de tous les pixels du détecteur ont bien été déterminés. Si ce n'est pas le cas, on continu la détermination des indices confiance, à l'étape 15. Si c'est le cas le procédé PCD se poursuit par l'étape 17 qui consiste à réaliser une première carte de confiance C1. Cette carte de confiance regroupe tous les indices de confiance déterminés dans l'étape 15. Le procédé PCD se poursuit par une étape 18 qui consiste à déterminer, sur l'image initiale Im, déterminée à l'étape 13, le gradient de niveau de gris de l'image. Cette étape consiste, de façon plus précise - à calculer l'image de gradient grad x(x, y) suivant la direction x de l'image corrigée en 13, par exemple par convolution de l'image Im_cor par un noyau ; - à calculer une image de gradient grad_y(x y) selon la direction Y de l'image Im_cor, par exemple par convolution de l'image Im_cor par un noyau ; - à calculer le module du gradient à partir de l'expression procédé PCD se poursuit par étape 19 de calcul de la variance et de la moyenne de l'image de gradient obtenue à l'étape 18. variance et la moyenne de l'image de gradient se déterminent, respectivement, par Moy_ grad =1 / nombre- pixel i x 1y Mod_ grad(x, Var_ grad = 1 / nombre- pixel lx 1y (Mod_ grad(x, y) Moy_ grad)2 . valeurs de variance et de moyenne peuvent éventuellement être améliorées en supprimant l'effet des points aberrants, comme cela a déj été expliqué pour calcul de la première variance et de la première moyenne à l'étape 14. Le procédé PCD se poursuit par une étape 20) qui consiste à déterminer l'indice de confiance de chaque pixel du détecteur, à partir de son niveau de gris et des moyenne et variance déterminées l'image de gradient. Autrement dit, l'indice de confiance de chaque pixel est déterminé à partir de - la normalisation des niveaux de gris de l'ensemble des pixels, afin d'obtenir une distribution normale - la transformation de la valeur par fonction F de répartition de la loi normale pour obtenir une distance normale dist norm_grad(x, y) = F(norm_grad(x, y)) - calcul de l'indice de confiance du pixel 10 par C2 y) = 2(1 - dist norm_grad(x, y)). Le procédé PCD comporte alors un test 21 qui consiste a vérifier que l'indice de confiance de chaque pixel du détecteur a bien été déterminé. Tant que ce n'est pas le cas, l'étape 20 est réinitialisée pour chacun des pixels. Lorsque c'est le cas, le procédé PCD se poursuit par une étape 22 qui consiste à réaliser une seconde carte de confiance C2 qui est constituée de l'ensemble de tous les indices de confiance de tous les pixels du détecteur déterminés à l'étape 20. Le procédé PCD se poursuit ensuite par une étape 23 qui consiste en une multiplication de première et de la seconde cartes de confiance, de façon à obtenir une troisième carte de confiance C3 (étape 24). Ce sont les indices de confiance de cette troisieme carte qui sont utilisés pour corriger les pixels défaut de l'image de l'objet à mesurer (étape 25) en effet, chaque pixel en défaut est remplacé par l'indice de confiance qui lui correspond sur la carte de confiance C3. La détermination de cette troisième carte de confiance, appelée carte de confiance finale, est déterminée à partir de l'équation C_final (x, y) = C1 (x, y) .C2 (x, y). La calibration du détecteur est effectuée une seule fois. Ensuite, on peut effectuer l'acquisition de plusieurs images sur le détecteur et les corriger, sans avoir à recalibrer le détecteur. L'étape 25 de correction consiste d'abord à effectuer l'acquisition d'une image de l'objet (25a) appelée image courante. L'étape 25b de correction de l'image peut être réalisée moyen d'un filtrage, tel que celui proposé par KNUTSSON, en utilisant la carte de confiance finale calculée dans l'étape 24 de calibration. Ce filtrage décrit, notamment, dans le document de H. KNUTSSON et . WESTING, Normalized and Differential Convolution, Method for Interpolation and Filtering of Incomplete and Uncertain Data , Proceedings CVPR 93, pp. 523, New York 1993. Dans le cas le plus simple (à l'ordre 0), la méthode peut être résumée de façon suivante - on considère une image I qui, par exemple, est l'image acquise sur le détecteur X ; - par ailleurs, on considère que l'on dispose d'une image de confiance C les pixels de l'image ; - on peut alors construire une nouvelle image filt_I de la façon suivante filt_I = (I.C)*a/C*a où - . et / sont les opérateurs multiplication et division d'image pixel à pixel ; - * est l'opérateur de convolution ; a est le noyau de convolution suivant image filt_I ainsi obtenue est une image pour laquelle - sur les zones qui ont un indice de confiance uniforme, l'effet du traitement est équivalent un filtrage ; - pour les pixels dont la confiance est inférieure à ses voisins, l'effet du traitement est de interpolation douce avec une nouvelle valeur dépendra en majorité des pixels voisins de forte confiance. La taille et la forme du noyau de convolution à adapter à la taille (diamètre minimum) des défauts que l'on veut corriger. Le filtrage, selon la méthode de KNUTSSON, permet de corriger les défauts d'image en utilisant au mieux l'information fournie par les pixels corrects de l'image. Par ailleurs,. dans le procédé PCD, afin d'éviter d'introduire un effet de filtrage, même sur les zones correctes où la confiance sur pixels est bonne, le procédé PCD propose de combiner l'image initiale et l'image corrigée de façon a créer une nouvelle image corrigée pondérée par les indices de confiance. Cela se traduit par l'équation Corr_I = Im.C_final + (1 - C_final).filt_I, filt_I est l'image filtrée par méthode de KNUTSSON. La correction des pixels en défaut partir de la troisième carte de confiance peut être réalisée par autres méthodes de filtrage, comme méthode de filtrage par équations aux dérivés partielles, décrite dans le document de R. DE RICHES et FAUGERAS, intitulé Les équations aux dérivés partielles en traitement des images et vision par ordinateur , publié dans la revue Traitement du signal , volume 13, n 6, pages 552 à 577, Février 1996. Une variante du procédé PCD concerne le cas où les défauts de l'image sont dus à l'absence mesure sur ou plusieurs zones du détecteur. Dans ce cas, le procédé PCD consiste à effectuer étapes supplémentaires après l'étape 24 de détermination de la troisième carte de confiance. Ces etapes sont représentées sur la figure 5. Dans cette variante, le procédé PCD comporte une etape 241 qui consiste à agrandir l'image de l'objet de façon à introduire les zones de pixels sans mesure (ou pixels manquants) où le détecteur n'a pas pu faire de mesure. En effet, l'image acquise est de plus petite taille que l'image qui aurait été acquise avec détecteur apte à détecter sur toute sa surface. Une étape 242 consiste ensuite à ajouter, dans ces zones, des pixels d'indices de confiance nuls. Une étape 243 consiste ensuite à réaliser une quatrième carte de confiance qui prend en compte les indices de confiance nuls des pixels de la zone sans mesure. Cette quatrième carte de confiance est alors considérée comme la carte finale. Le détecteur étant calibré, une étape 244 consiste ensuite à effectuer une nouvelle acquisition d'image l'obj à mesurer, ou image courante. Cette dernière est agrandie (étape 245) puis corrigée. C'est à partir de cette image courante agrandie, et en appliquant quatrieme carte de confiance, que l'image est corrigée, dans l'étape 25 décrite précédemment. En d'autres termes, cette variante consiste à introduire les zones sans information, c'est-à-dire les pixels pour lesquels le détecteur n'a pas pu faire de mesure, dans l'image de confiance finale établie à l'étape 24 ; cela se fait en augmentant la taille de l'image de confiance. Pour cela, on coupe l'image au niveau des zones manquantes et on ajoute des pixels à indices de confiance nuls pour remplir les vides de l'image. REVENDICATIONS 1. Procédé d'utilisation d'un systeme (1) d ostéodensitométrie par rayons X à au moins deux énergies, à faisceau conique, ce système comprenant une source (1a) de rayons X apte à fournir un faisceau conique de rayons X à au moins une première énergie, appelée haute énergie, et une deuxième énergie, appelée basse énergie et inférieure à la première énergie, un detecteur bidimensionnel (2) de rayons X et moyens électroniques (3a, 3b) de traitement des images fournies par ce détecteur, ce procédé étant caractérisé en ce qu'on acquiert l'image à basse énergie d'une zone anatomique d'un patient (1c), on acquiert l'image à haute énergie de cette zone anatomique et l'on met en correspondance ces images respectivement acquises à basse énergie et à haute énergie, avant de construire la carte des densités osseuses de la zone anatomique. 2. Procédé selon la revendication 1, dans lequel, pour mettre en correspondance les images respectivement acquises à haute et basse énergies - on extrait respectivement de ces images les ensembles de contours des zones osseuses de la zone anatomique, - on recherche la transformation plane optimale permettant de mettre l'ensemble de contours relatif à l'image acquise à haute énergie en correspondance avec l'ensemble de contours relatif à l'image acquise à basse énergie, - on utilise cette transformation pour amener l'image acquise à haute énergie dans le référentiel de l'image acquise à basse énergie, et - on combine ces deux images pour déterminer la carte de mesure des densités osseuses. 3. Procédé selon la revendication 1, dans lequel, pour mettre en correspondance les images respectivement acquises à haute et basse énergies - on extrait respectivement de ces images les ensembles de contours des zones osseuses de la zone anatomique, - on recherche la transformation plane optimale permettant de mettre l'ensemble de contours relatif à l'image acquise à basse énergie en correspondance avec l'ensemble de contours relatif à l'image acquise à haute énergie, on utilise cette transformation pour amener l'image acquise à basse énergie dans le référentiel de l'image acquise à haute énergie, et on combine ces deux images pour déterminer la carte de mesure des densités osseuses. 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel, avant de faire les acquisitions à haute et basse énergies on fait un cliche, à faible dose de rayons X, de la zone anatomique du patient (1c) pour aider positionner cette zone dans le système. 5. Procédé selon la revendication 4, dans lequel, lorsque le patient (1c) est soumis à un premier examen, on utilise le cliché à faible dose de rayons X pour rétroagir sur la mécanique du système (1) afin de positionner la zone anatomique par rapport à référentiel préétabli. . Procédé selon la revendication 4, dans lequel lorsque le patient (1c) est soumis à un examen suivi on utilise le cliché à faible dose de rayons X pour placer la zone anatomique dans une position identique à celle qu'elle occupait lors de l'examen précédent.