FR2918487A1 - Procede de traitement d'images en tomosynthese pour une detection d'opacites - Google Patents

Abstract

La présente invention a pour objet un procédé de traitement d'images de projections (P1-Pn) en tomosynthèse pour une détection (26, 29) sélective d'opacités. Le procédé de l'invention permet de détecter automatiquement des opacités, dans des séries de données tridimensionnelles. L'invention met en oeuvre un algorithme permettant d'accroître les données à prendre en compte dans la détection d'opacités en rassemblant toutes les informations disponibles sur lesdites opacités rendant ainsi les critères de sélection des opacités plus sélectives. Ce rassemblement est effectué par une collecte de l'ensemble des données étendues dans plusieurs coupes (40-44) consécutives.

Description

1 Procédé de traitement d'images en tomosynthèse pour une détection

d'opacités

Domaine de l'invention La présente invention a pour objet un procédé de traitement d'images de projections en tomosynthèse pour une détection sélective d'opacités. La présente invention trouve des applications particulièrement avantageuses, mais non exclusives, dans le domaine de l'imagerie médicale et plus particulièrement celui de la mammographie. La présente invention concerne également une station de revue d'images médicales comportant un logiciel de mise en oeuvre d'un tel procédé de traitement d'images. Etat de la technique La mammographie est aujourd'hui largement utilisée pour la détection de lésions et la prévention du cancer du sein. Les signes à détecter sur des images de mammographie sont des signes radiologiques associés à ces lésions. Ces signes sont principalement soit des dépôts de calcium soit des opacités. Les opacités sont des zones denses, où les rayons X sont plus fortement absorbés que dans des zones adjacentes. Une opacité comporte souvent un noyau dense possiblement avec des spicules convergeant vers le noyau. Souvent, la forme du noyau peut être grossièrement décrite par une boule ou un ellipsoïde sur l'image radiographique. Aujourd'hui, il reste difficile, même pour des praticiens, de détecter visuellement les opacités sur des images radiographiques standard. En effet, bien que les opacités soient des éléments de grande taille, de 0.5 centimètre à 3 centimètres de diamètre environ, leur forme n'est pas bien définie et leur contraste est souvent faible. En outre, Il arrive que certaines zones denses ne soient pas repérées par le radiologue, sur des images de mammographies standard. Ce phénomène a plusieurs causes. Notamment, les images de mammographies standard étant le résultat de projections, elles représentent des structures superposées qui perturbent la visibilité des structures du sein, entraînant parfois, soit une interprétation faussement positive lorsqu'une superposition de ces structures ressemble à une lésion, soit une interprétation faussement négative, lorsque ces structures obscurcissent la visibilité d'une lésion.

Pour aider à la résolution de ces problèmes d'interprétations

2 faussement positives ou faussement négatives, il existe dans l'état de la technique de nouveaux dispositifs de mammographie produisant une image tridimensionnelle du sein de la patiente. Ces dispositifs permettent d'acquérir plusieurs projections d'un objet sous différents angles puis de reconstruire la répartition tridimensionnelle de cet objet à l'aide d'un algorithme de reconstruction tomographique. Ce procédé utilisant des projections acquises sur une excursion angulaire réduite est appelé tomosynthèse. Cependant, ces nouveaux dispositifs de mammographie en tomosynthèse présentent des inconvénients. En effet, dans de tels dispositifs, on réalise une reconstruction en volume numérique de 50 à 80 coupes, pour un sein moyen. De ce fait, la quantité d'information à produire est très importante. De même, l'accès à une information présentant un intérêt clinique prend beaucoup plus de temps du fait que cette information est recherchée séquentiellement dans le volume, coupe par coupe.

En outre, si l'espace entre les coupes est d'un millimètre, l'information concernant la présence d'une opacité peut être dispersée sur 5 à 30 coupes consécutives, en fonction de la taille de l'opacité à détecter. Exposé de l'invention L'invention a justement pour but de remédier aux inconvénients des 20 techniques exposées précédemment. Pour cela, l'invention propose un procédé de traitement d'images permettant de détecter automatiquement des opacités, dans des séries de données tridimensionnelles. Le temps passé pour localiser ces opacités est fortement diminué. 25 Aussi, le système peut éviter un éventuel oubli de lésion(s) de la part du radiologue, lors de son examen. Avec ce procédé de traitement d'images, les opacités sont plus faciles à identifier, permettant ainsi un examen plus efficace. L'invention propose une nouvelle approche de la mammographie, qui 30 exploite avantageusement les techniques de traitement numérique de l'image radiologique, pour améliorer la lisibilité de l'information. Ce procédé de traitement d'images permet également de faire face à une quantité importante de données en développant de nouvelles stratégies aptes à diminuer le temps de revue des données et à simplifier l'accès à 35 l'information d'intérêt clinique.

3 Pour ce faire, l'invention met en oeuvre un algorithme permettant de prendre en compte des données réparties sur plusieurs coupes pour la détection d'opacités. Avec cet algorithme, le procédé de détection des opacités de l'invention rassemble toutes les informations disponibles sur lesdites opacités avant d'appliquer des critères de sélection rendant la détection des opacités plus sélective. Pour effectuer cette collecte de données, l'algorithme réunit l'ensemble des données disponibles concernant la présence d'opacités étendues dans plusieurs coupes consécutives. Pour ce faire, les coupes sont associées en des blocs de coupes. L'épaisseur de bloc peut dépendre de la taille de l'opacité à détecter ou de l'élongation d'une particule composée de voxels détectés comme éléments d'opacité potentiels. Au lieu d'appliquer à chaque coupe les critères de sélection d'une opacité comme dans l'état de la technique, l'invention les applique à chaque bloc de coupes. Cette application à chaque bloc permet d'augmenter la précision dans la sélection des opacités. L'invention met également en oeuvre un algorithme de calcul du contraste à appliquer aux coupes associées ou aux coupes reconstruites ou aux projections, en fonction des différents modes de réalisation de l'invention. Le contraste calculé en tout élément ou particule d'éléments du volume est comparé à un seuil de contraste prédéfini d'opacités. Lorsque le contraste calculé est supérieur au seuil de contraste alors l'élément de volume correspondant est validé comme un élément de d'opacité et reçoit un attribut d'opacité.

La présente invention propose également un algorithme de visualisation d'opacités dans des séries de données tridimensionnelles. L'algorithme de l'invention permet d'indiquer la position des opacités sur l'image tridimensionnelle par une colorisation ou un rehaussement de leur intensité et/ou, sur option, les opacités identifiées sont marquées et/ou entourées par un contour. Plus précisément, l'invention a pour objet un procédé de traitement d'images produites par un dispositif à rayons X, dans lequel, - on réalise plusieurs projections radiographiques d'un corps, - on reconstruit un volume numérique en coupes du corps, à partir de l'ensemble des projections,

4 caractérisé en ce que - on associe des coupes consécutives du corps en des blocs de coupes en fonction d'un paramètre d'épaisseur du bloc, - une épaisseur étant caractérisée par une taille d'opacités à analyser, - on détermine un nombre suffisant d'épaisseurs, pour couvrir une gamme de tailles des opacités, le nombre correspondant au nombre des tailles, - on applique à chaque bloc de coupes un algorithme de détection d'opacités, -on donne à chacun des voxels du bloc, correspondant à une opacité détectée, un attribut représentatif d'opacité, - on visualise les voxels ayant un attribut d'opacité dans un volume numérique de présentation, produit à partir des projections. L'invention concerne également un dispositif à rayons X apte à mettre en oeuvre ledit procédé de traitement d'images en tomosynthèse, de l'invention, pour une détection d'opacités. Brève description des dessins L'invention sera mieux comprise à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent. Celles-ci ne sont présentées qu'à titre indicatif et nullement limitatif de l'invention. La figure 1 montre une représentation schématique d'un dispositif de rayonnement X, en particulier un mammographe, muni des moyens perfectionnés de l'invention. La figure 2 montre une représentation schématique des coupes d'un 25 volume numérique associées en bloc de coupes. La figure 3 montre une illustration de moyens mettant en oeuvre un premier mode de réalisation du procédé de l'invention, dans une première stratégie de détection. La figure 4 montre une illustration de moyens mettant en oeuvre un 30 deuxième mode de réalisation du procédé de l'invention, dans la première stratégie de détection. La figure 5 montre une illustration de moyens mettant en oeuvre un troisième mode de réalisation du procédé de l'invention, dans la première stratégie de détection. 35 La figure 6 montre une illustration de moyens mettant en oeuvre un quatrième mode de réalisation du procédé de l'invention, dans une deuxième stratégie de détection. Les figures 7 et 8 montrent une illustration de moyens mettant en oeuvre respectivement un cinquième et un sixième mode de réalisation du 5 procédé de l'invention, dans une troisième stratégie de détection. La figure 9 montre une illustration de moyens mettant en oeuvre un septième mode de réalisation du procédé de l'invention, dans une quatrième stratégie de détection. Description détaillée de modes de réalisation de l'invention La figure 1 montre un dispositif de rayonnement X, particulièrement un mammographe, conforme à l'invention. Ce dispositif 1 de rayonnement X comporte une colonne verticale 2. A cette colonne verticale 2 est articulé un bras 3 portant un tube 4 émetteur de rayonnement X. Un détecteur 5 porté par la colonne 2 est apte à détecter les rayons X émis par le tube 4. Ce bras 3 et le détecteur 5 peuvent être orientés à la verticale, à l'horizontale ou en oblique. Le tube 4 est muni d'un foyer 6 qui est le foyer émetteur des rayons X. Ce foyer 6 émet un faisceau 7 de rayons X suivant une direction d'émission D. Le bras 3 est articulé à la colonne verticale 2 de telle sorte qu'il permet de déplacer le tube 4 suivant une trajectoire en forme d'arc de cercle tout en laissant le détecteur 5 immobile pendant une acquisition. D'autres arrangements sont possibles permettant au tube 4 de se déplacer dans un plan ou dans une portion de sphère. Le tube 4 peut alors occuper différentes positions réparties en basculement entre deux positions extrêmes. Ces deux positions sont par exemple symétriques l'une de l'autre, par rapport à la perpendiculaire au plan du détecteur. Le détecteur 5 est dans un exemple préféré un détecteur numérique. Un plateau porte-sein 8 est superposé au détecteur 5. Le détecteur 5 détecte les rayons X ayant traversé le sein de la patiente et le plateau porte sein 8.

En outre, pour des raisons tant d'immobilisation du sein que de qualité d'image ou d'énergie de rayonnement X délivrée au sein de la patiente, il est nécessaire de comprimer le sein de la patiente, au moment de la radiographie. Différents efforts de compression peuvent être appliqués. Ces efforts sont appliqués via une pelote 9 de compression qui comprime le sein sur le plateau porte-sein 8, en fonction du type d'examen que l'on souhaite

6 faire. A cet effet le système comporte une pelote 9 coulissante et susceptible de venir comprimer le sein perpendiculairement au détecteur, manuellement ou d'une manière motorisée. La pelote 9 est formée en un matériau transparent aux rayons X, notamment en matière plastique.

Alors que la pelote 9, le sein de la patiente, le plateau 8 et le détecteur 5 sont fixes, le tube à rayon X 4 peut occuper diverses positions dans l'espace par rapport à cet ensemble. Dans une variante, le détecteur 5 peut être mobile et occuper diverses positions autour du sein en même temps que le tube à rayons X 4. Dans ce cas, le détecteur 5 n'est plus solidaire au plateau porte sein 8. Le détecteur 5 peut avoir une forme plane ou courbe. Il peut être déplacé en rotation et/ou en translation. En réponse à la réception d'un faisceau 7 qui traverse le sein de la patiente, le détecteur 5 émet des signaux électriques correspondant à l'énergie des rayons reçus. Ces signaux électriques peuvent ensuite être transmis à un circuit logique de commande 10 par l'intermédiaire d'un bus externe 11. Ces signaux électriques permettent à cette logique de commande 10 de produire une image correspondant au sein analysé. Ces images peuvent être visualisées à l'aide d'un écran de cette logique de commande 10 ou imprimées. Dans un exemple, le circuit logique de commande 10 comporte un microprocesseur 12, une mémoire 13 de programme, une mémoire 14 de données, un écran 15 de visualisation, un clavier 16 et une interface 17 d'entrée sortie. Le microprocesseur 12, la mémoire 13, la mémoire 14, l'écran 15, le clavier 16 et l'interface 17 d'entrée sortie sont interconnectés par un bus interne 18. Dans la pratique, lorsque l'on fait réaliser une action par un tel circuit de commande, cette action est réalisée par le microprocesseur 12 de ce circuit. Le microprocesseur reçoit dans son registre des instructions lues dans la mémoire 13 de programme. Le circuit logique de commande 10 est souvent réalisée sous forme de circuit intégré. La mémoire 13 de programme est divisée en plusieurs zones, chaque zone correspondant à des codes instructions pour réaliser une fonction du circuit 10. La mémoire 13 comporte, selon des variantes de stratégies de l'invention, un programme enregistré dans une zone 19 comportant des

7 codes instructions. Ces codes instructions servent à établir une trajectoire du tube 4 et à commander l'acquisition d'une multiplicité d'images de projection le long de cette trajectoire. Une zone 20 comporte des codes instructions pour réaliser un pré-traitement aux images de projection dans le but de supprimer des artefacts dus à la géométrie d'acquisition et au détecteur. Pour toutes les stratégies, la mémoire 13 comporte une zone 21 comportant des codes instructions pour appliquer un algorithme de reconstruction, en particulier tomographique, à l'ensemble des images de projection acquises à différents angles afin d'obtenir un volume numérique reconstruit. Ce volume numérique comporte par exemple 50 à 80 images de coupes, les coupes étant espacées de préférence d'un millimètre. Cette zone 21 comporte également des codes instructions pour appliquer à l'ensemble des images de projection acquises un deuxième algorithme de reconstruction davantage orienté vers un objectif de visualisation. Ce deuxième algorithme est de préférence un algorithme de reconstruction algébrique. Bien que ce deuxième algorithme soit utile, il serait possible d'utiliser deux fois le premier algorithme. Selon une première et une deuxième stratégie, la mémoire 13 comporte une zone 22 comportant des codes instructions pour associer des coupes consécutives du volume numérique reconstruit en bloc de coupes en fonction d'un paramètre d'épaisseur du bloc, typiquement un nombre de coupes dans ce bloc. Dans un exemple, ce paramètre d'épaisseur est caractérisé comme étant supérieur ou égal à une taille d'opacités à analyser. Une zone 23 comporte des codes instructions pour déterminer un nombre suffisant de paramètres d'épaisseur de bloc de coupes. Ce nombre suffisant est tel qu'il permet de couvrir une gamme des tailles possibles des opacités. Par exemple si ce nombre vaut trois, ce choix revient à considérer qu'il y a trois tailles d'opacités à considérer : par exemple les tailles correspondent à 6 coupes, 7 coupes et 8 coupes. Les différentes tailles ne sont toutefois pas nécessairement en progression régulières. Une zone 24 comporte des codes instructions pour déterminer un filtre bidimensionnel à appliquer à chaque coupe reconstruite ou à chaque projection acquise ou à chaque bloc de coupes, en fonction des différents modes de réalisation de l'invention.

Une zone 25 comporte des codes instructions pour appliquer un

8 algorithme de calcul du contraste normalisé en chaque élément de chaque coupe reconstruite ou de chaque bloc de coupes, en fonction des différents modes de réalisation de l'invention. La zone 25 peut être décomposé en deux parties, une partie comportant un programme de calcul de contraste, une autre partie comportant un programme de normalisation de ce contraste. Une zone 26 comporte des codes instructions pour appliquer un algorithme de détection ou d'identification, d'élément d'opacité, en fonction du contraste normalisé, de chaque élément de chaque coupe reconstruite ou de chaque bloc de coupes, en fonction des différents modes de réalisation de l'invention. Les codes instructions de la zone 26 attribuent un attribut d'opacité aux éléments d'opacité. Ils attribuent à tous les autres éléments un attribut de non opacité. Un attribut d'opacité ou de non opacité est un paramètre, au moins temporaire et traité par le circuit. Concernant une troisième stratégie, les attributs d'opacité ou de non opacité, des éléments du volume sont provisoires. Une zone 27 comporte des codes instructions pour appliquer un algorithme de composantes connexes tridimensionnel. Cet algorithme est appliqué sur un volume binaire composé d'un ensemble d'éléments à 1, pour des éléments ayant un attribut d'opacité, et à 0 ailleurs.

Cet algorithme permet de relier entre eux des éléments d'opacité qui sont détectés par l'algorithme de détection et d'affecter un attribut d'opacité provisoire aux éléments intermédiaires. Une zone 28 comporte des codes instructions pour calculer un contraste d'ensemble d'un ensemble d'éléments. A chaque élément de cet ensemble, on affecte ensuite un attribut de contraste égal au contraste d'ensemble dudit ensemble. Une zone 29 comporte des codes instructions pour comparer le contraste d'ensemble de chaque élément à un seuil prédéfini. Lorsque le contraste d'ensemble est supérieur au seuil, on attribue audit élément correspondant un attribut d'opacité définitif. Selon, une quatrième stratégie, une zone 30 comporte des codes instructions pour déterminer un filtre tridimensionnel à appliquer au volume reconstruit. Par la suite, on verra que le filtrage est appliqué au volume entier et est suivi d'une normalisation et d'un seuillage.

Une zone 31 comporte des codes instructions pour calculer un

9 contraste normalisé en chaque élément du volume. Une zone 32 comporte des codes instructions pour appliquer un algorithme de détection d'éléments représentatif d'opacité, en fonction du contraste normalisé de chaque élément du volume. Les codes instructions de la zone 32 attribuent un attribut d'opacité aux éléments représentatifs d'opacité. Pour toutes les stratégies, la mémoire 13 comporte une zone 33 comportant des codes instructions pour coloriser des éléments d'opacités en leur attribuant une information de couleur, dans le volume numérique de présentation. Une zone 34 comporte des codes instructions pour rehausser le signal des éléments d'opacités dans le volume numérique de présentation. Une zone 35 comporte des codes instructions pour afficher le volume de présentation comportant les éléments d'opacités rehaussés ou colorisés, sur l'écran 15 de visualisation. En mode de fonctionnement, le sein et donc le détecteur 5 sont irradiés pendant des poses consécutives. Pour ces poses, le foyer 6 du tube à rayons X occupe dans l'espace des positions qui peuvent être fixes, de préférence régulièrement réparties. Dans un exemple, et bien que cela ne puisse pas être considéré comme une limitation de l'invention, on prévoit que l'exploration angulaire maximale soit ainsi de 60 degrés, de plus ou moins 30 degrés par rapport à une direction centrale d'irradiation, en général verticale par rapport au détecteur pour un mammographe. Au cours de cette exploration, on procède à l'émission d'un certain nombre de faisceaux 7, par exemple neuf, onze, treize, quinze ou autre, en fonction de la précision souhaitée. Dans l'exemple de la figure 1, le nombre de faisceaux 7 est de neuf. Chaque faisceau 7 représente une projection radiographique. La figure 2 montre schématiquement un exemple d'un bloc 50 de coupes 40 à 44 du volume V numérique reconstruit. Ce bloc 50 se développe selon une épaisseur définie par un paramètre 60 d'épaisseur. Ce paramètre 60 d'épaisseur est caractérisé par une taille d'opacité à détecter. Plus la taille des opacités à détecter est grande plus on choisira une épaisseur grande, et réciproquement plus petite dans le cas contraire. L'information concernant la présence d'une opacité est éparpillée sur les coupes 40 à 44 du volume.

10 Pour réunir l'ensemble de cette information, le circuit logique de commande associe l'ensemble des coupes consécutives comportant chacune une partie de l'information. Ainsi, le circuit logique de commande associe des coupes 40 à 44 consécutives du volume V en des blocs 50 de coupes en fonction du paramètre d'épaisseur. Le circuit logique de commande détermine un nombre suffisant d'épaisseurs de blocs de coupes pour couvrir une gamme de tailles d'opacités. Par exemple, elle pourra déterminer trois épaisseurs de bloc 50 pour correspondre à trois tailles d'opacité.

Les blocs 50 de coupes sont formés par une association de coupes 40 à 44. Les éléments du volume des coupes 40 à 44 du bloc 50 reçoivent un attribut résultant d'une opération mathématique. Cette opération mathématique peut être notamment une somme, un maximum, un minimum, une moyenne du niveau de gris selon un axe perpendiculaire au détecteur.

De préférence, cette opération mathématique peut être une combinaison logique ou algébrique de niveaux de gris dans les coupes, suivant une direction d'émission des rayons X de la source au détecteur. Dans l'exemple de la figure 2, le paramètre d'épaisseur est égal à 5 millimètres. Comme les coupes 40 à 44 sont espacées d'un millimètre, alors le circuit logique de commande associe cinq coupes 40 à 44 consécutives en blocs 50 de coupes. Le circuit logique de commande applique à chaque bloc 50 de coupes, l'algorithme de détection des voxels d'opacités. Les figures 3 à 9 montrent différentes illustrations de moyens mettant en oeuvre le procédé de traitement d'images de l'invention, permettant de localiser des éléments susceptibles de constituer des opacités, dans le volume du sein reconstruit, en tenant compte de la dispersion des données d'opacité sur plusieurs coupes. Le procédé de l'invention étant appliqué à un volume numérique, dans ce cas, l'élément du volume est un voxel. Les figures 3 à 5 montrent une première stratégie de détection d'opacités utilisant la méthode de collecte de données décrite à la figure 2. La figure 6 montre une deuxième stratégie de détection d'opacités utilisant la méthode de collecte de données décrite à la figure 2. Les figures 7 et 8 montrent une troisième stratégie de détection d'opacités. La figure 9 montre une quatrième stratégie de détection d'opacités.

La figure 3 montre un premier mode de réalisation de l'invention. Sur

11 la figure 3, le tube 4 à rayons X émet, à l'étape 100, des intensités de rayonnement X traversant le sein de la patiente pour une multiplicité de projections P1 à Pn, selon une trajectoire préalablement déterminée. Le détecteur 5 acquiert les données images 11 à ln de rayons X représentant respectivement les projections P1 à Pn. Le circuit logique de commande 10 traite les données images de rayons X 11 à ln. A l'étape 101, le circuit logique de commande 10 applique un algorithme de reconstruction. Dans un exemple préféré, l'algorithme de reconstruction utilisé est celui connu sous le nom de Simple Back Projection en anglais ou rétroprojection simple. Cet algorithme permet de reconstruire le volume dans différents plans de coupes parallèles au détecteur. En particulier on propose une tomosynthèse du sein. Cette technique de tomosynthèse permet, à partir d'un faible nombre de projections bidimensionnelles ou données images, réparties sur un domaine angulaire restreint et acquises sur un détecteur numérique, de reconstruire le volume tridimensionnel du sein étudié. A ce stade, l'algorithme peut appliquer un algorithme de rehaussement du bord du sein afin de compenser la chute brutale de signal due à une décroissance rapide de l'épaisseur du sein sur son bord.

A l'étape 102, le circuit logique de commande 10 explore une taille xm d'opacité à détecter dans une gamme de tailles préalablement définie. Cette taille xm varie par exemple de 5 millimètres à 3 centimètres. En pratique, les opacités sont considérées comme sphériques et donc la taille correspond à son diamètre. Puis, le circuit logique de commande 10 détermine la valeur du paramètre d'épaisseur en fonction de la taille xm de l'opacité à détecter. Par exemple, une table de correspondance est mémorisée dans la mémoire 14 de données. Cette table permet de fixer les tailles xm à considérer et le nombre de coupes correspondant à traiter. Le circuit logique de commande 10 effectue une association de coupes en bloc de coupes pour la valeur déterminée du paramètre d'épaisseur. Cette association des coupes en bloc de coupes est répétée autant de fois que nécessaire pour couvrir la gamme complète de tailles xm d'opacités. Dans une variante, la valeur du paramètre d'épaisseur peut être fixée à 1, ce qui revient à traiter les coupes reconstruites indépendamment les 35 unes des autres.

12 A l'étape 103, le circuit logique de commande 10 détermine un moyen de calcul du contraste en chaque voxel de chaque bloc de coupes. Ce moyen de calcul du contraste peut être effectué par tout type de filtrage existant permettant de calculer un contraste dans une image. Dans un exemple préféré, le moyen de calcul du contraste est mis en oeuvre par un filtre d'ondelette bidimensionnel de type chapeau mexicain, de paramètre d'échelle s. Le circuit logique de commande 10 choisit l'échelle suivant la taille du de l'opacité à détecter. Ainsi, le paramètre d'échelle du filtre dépend de la taille xm de l'opacité à détecter.

Le circuit logique de commande 10 réalise une convolution de chaque bloc de coupes, déterminée à l'étape 102, avec le filtre d'ondelette, pour chaque taille xm d'opacité à détecter. Lors, de la convolution de chaque bloc de coupe avec rondelette, le circuit logique de commande calcule à chaque voxel une moyenne des niveaux de gris voisins pondérés par la valeur absolue des coefficients négatifs du filtre et une moyenne des niveaux de gris voisins pondérés par des coefficients positifs du filtre. Enfin, la logique de commande calcule le contraste AI de chaque voxel du bloc de coupe en faisant la différence de ces deux moyennes.

L'utilisation du filtre d'ondelette permet d'effectuer un calcul de contraste linéaire en faisant une différence des valeurs moyennes de niveau de gris. A l'étape 104, le circuit logique de commande 10 normalise le contraste AI calculé en chaque voxel de chaque bloc de coupes. Dans un exemple préféré, cette normalisation s'effectue en fonction d'une intensité de fond lb et de la taille xm de l'opacité à détecter. L'intensité de fond lb peut être estimée par la moyenne des niveaux de gris dans un environnement immédiat de l'opacité à détecter que l'on suppose centrée sur le voxel d'intérêt. Cette normalisation définit une différence d'atténuation 4 . Cette différence d'atténuation 4 est égale à O,u = DI lb . xm Le circuit logique de commande 10 peut utiliser d'autres méthodes de calcul et de normalisation de contraste. Cette normalisation n'est pas impérative. A l'étape 105, le circuit logique de commande 10 applique un algorithme d'identification des voxels d'opacités sur tous les blocs de coupes filtrés du volume numérique. Pour chaque bloc de coupes dont l'épaisseur dépend de la taille xm de l'opacité à détecter, lecircuit logique de commande 10 attribue à chaque voxel un attribut d'opacité ou de non opacité. Pour cet algorithme d'identification, le circuit logique de commande 10 vérifie si le contraste calculé de chaque voxel est suffisant pour être considéré comme une opacité. Dans un exemple préféré, le circuit logique de commande 10 effectue un seuillage de chaque bloc de coupes en comparant le contraste normalisé de chaque voxel à un seuil O,un,;n de contraste normalisé préalablement défini. Le circuit logique de commande 10 attribue un attribut d'opacité aux voxels de chaque bloc de coupes dont le contraste est supérieur à ce seuil O,un,;n de contraste. Le seuil O,un,;nest fixé soit par expérimentation, soit il est laissé à l'initiative d'un praticien qui peut le régler à sa volonté. Les étapes 102 à 104 et l'algorithme 105 de détection des éléments d'opacités sont répétées autant de fois que nécessaire pour couvrir la gamme complète de tailles xm d'opacités. Ainsi, seuls les éléments pour lesquels 4 > O,un,;n pour au moins une taille xm d'opacité à détecter sont considérés comme élément d'opacité. A l'étape 106, le circuit logique de commande 10 applique un deuxième algorithme de reconstruction, à partir des données images 11 à ln. Ce deuxième algorithme de reconstruction est davantage orienté vers un objectif de visualisation comme la technique de reconstruction algébrique. Il fournit un volume de présentation avec des voxels ayant des même coordonnées que ceux produit par l'algorithme de reconstruction de l'étape 101 A l'étape 107, le circuit logique de commande 10 applique un algorithme de visualisation des voxels issus de l'algorithme de détection 105, en leur donnant un rehaussement correspondant à leur attribut d'opacité obtenu à l'étape 105, afin de faciliter l'analyse du praticien. Cet algorithme 107 de visualisation peut appliquer un processus de rehaussement de l'intensité des voxels ayant un attribut d'opacité ou peut leur attribuer une information de couleur, ainsi qu'à leur environnement immédiat. La première étape du processus 107 de rehaussement consiste à appliquer un algorithme de composantes connexes 3D sur un volume binaire composé d'un ensemble de voxel avec un attribut à 1 pour des positions d'opacité et à 0 ailleurs. Ceci mène à un ensemble de particules tridimensionnelles formées de voxels d'opacité connectés. Puis, pour chaque voxel du volume de présentation obtenu à l'étape 107 appartenant à une particule et à son environnement immédiat, l'intensité du voxel est augmentée par rapport à l'intensité de l'environnement immédiat de la particule. Dans un mode de réalisation préféré, le circuit logique de commande calcule la moyenne de niveaux gris Mp dans l'environnement immédiat de la particule. L'environnement immédiat d'une particule est défini comme l'ensemble des éléments situés à une distance d, préalablement définie, de tout élément de la particule. Cette moyenne de niveaux de gris Mp estime l'intensité de fond de la particule. Pour chaque voxel (xi, yj, zk) abrègé en (i, j, k) de la particule et de son environnement immédiat, le volume rehaussé noté DCA(i, j, k) est obtenu en amplifiant la différence entre l'intensité S(i, j,k) du voxel de la particule et l'intensité de fond Mp de la particule selon l'équation suivante : d2 DCA (i,j,k)=M+(1+A.e r2).[S(i,j,k)-M] où A influence directement la puissance du rehaussement et d est la distance entre le voxel (i,j,k) et la particule avec d=0 pour chaque voxel à l'intérieur de la particule. Puisqu'il est difficile d'identifier des découpes d'opacité, le circuit logique de commande applique un rehaussement dont la force diminue avec la distance à la particule. Le paramètre i contrôle cette diminution du facteur de rehaussement avec la distance. Dans une variante, la valeur des voxels ayant un attribut d'opacité peut être forcée à la valeur maximale (respectivement minimale) de la dynamique d'affichage permettant ainsi de ramener au premier plan les opacités dans une vue tridimensionnelle lorsque celle-ci utilise l'opérateur de l'intensité maximale (respectivement minimale) plus connu sous le nom anglais de Maximum Intensity Pixel (respectivement Minimum Intensity Pixel). Dans une variante, les opacités sont indiquées sur le volume de présentation par des marqueurs crées par le circuit logique de commande. Ces marqueurs peuvent être une annotation graphique notamment, un triangle, un carré, voire une lettre ou un chiffre, se situant au centre de gravité des opacités.

15 Dans une variante, le circuit logique de commande peut avoir pour conséquence de produire une surface englobante autour de chaque particule (. Pendant la visualisation des coupes, un contour résultant de la coupe de cette surface englobante peut être matérialisé par un trait visible sur la coupe regardée par le praticien, autour de l'opacité. Ceci peut avoir un intérêt si on veut que l'attention du praticien soit plus spécifiquement dirigée sur les lésions détectées. Si on veut aider le radiologue dans sa tâche de détermination de l'extension des lésions, ce contour peut épouser les limites des lésions.

Dans le cas où la visualisation des résultats se fait par colorisation des voxels ayant un attribut d'opacité, une couleur leur est attribuée en fonction de la quantité S(i, j,k) ùMp et de la distance d. Le volume de présentation final sur lequel les particules d'opacités ont été rehaussées ou colorisées est affiché sur l'écran de visualisation en vu d'être visualisé par un praticien. Le radiologue obtient ainsi une vue d'ensemble immédiate de la cartographie de répartition des opacités dans le volume. La figure 4 montre une illustration de moyens mettant en oeuvre un deuxième mode de réalisation du procédé de l'invention. L'étape 200 de la figure 4 correspond à l'étape 100 de la figure 3 où le détecteur 5 acquiert les données images 11 à ln de rayons X, émis par le tube 4, représentant respectivement les projections P1 à Pn. L'étape 201 de la figure 4 correspond à l'étape 101 de la figure 3 où le circuit logique de commande applique l'algorithme de reconstruction aux données images 11 à ln pour obtenir un volume numérique du sein. A l'étape 202, le circuit logique de commande applique à chaque coupe du volume le moyen de calcul du contraste de l'étape 103 de la figure 3. A l'étape 203, le circuit logique de commande effectue une association de coupes filtrées en bloc de coupes filtrées pour une valeur déterminée du paramètre d'épaisseur afin d'obtenir en tout voxel de chaque bloc de coupes un contraste. L'étape 204 de la figure 4 correspond à l'étape 104 bis de la figure 3 où le circuit logique de commande normalise le contraste AI calculé en chaque voxel de chaque bloc de coupes filtrées. Après l'étape 204, le circuit logique de commande applique à chaque bloc de coupes, l'étape 205 correspondante à l'étape 105 de détection

16 d'opacités de la figure 3. Après l'étape 205, le circuit logique de commande applique les étapes 206 et 207 correspondantes respectivement aux étapes 106 et 107 de la figure 3. Les étapes 202 à 205 sont répétées autant de fois que nécessaire pour toutes les tailles de la gamme de tailles d'opacités. Avec ces deux modes de réalisation de l'invention, figures 3 et 4, le circuit logique de commande doit convoluer chaque coupe ou chaque bloc de coupes avec le moyen de calcul du contraste. Ceci est informatiquement exigeant sur le plan des calculs, même pour des coupes espacées de 1 millimètre. L'invention propose pour résoudre ce problème un troisième mode de réalisation comportant un algorithme plus rapide qui dépend seulement du nombre de projections et pas du nombre de coupes ni du nombre de bloc de coupes, comme illustré à la figure 5. Le bon fonctionnement de cet algorithme nécessite une géométrie d'acquisition dans laquelle le facteur d'agrandissement varie peu en fonction de l'angulation du tube. La figure 5 montre ainsi un troisième mode de réalisation du procédé de l'invention. L'étape 300 de la figure 5 correspond à l'étape 100 de la figure 3 où le détecteur 5 acquiert les données images 11 à ln de rayons X, émis par le tube 4, représentant respectivement les projections P1 à Pn. A chaque image brute 11 à ln, le circuit logique de commande applique à l'étape 301 un moyen de calcul du contraste comme précédemment décrit à l'étape 103 de la figure 3. A l'étape 302, le circuit logique de commande applique un algorithme 25 de reconstruction d'un volume numérique en coupes à partir des images brutes 11 à ln filtrées. A l'étape 303, le circuit logique de commande effectue une association des coupes filtrées obtenues en bloc de coupes filtrées, en fonction de la valeur du paramètre d'épaisseur. L'étape 304 de la figure 5 correspond à 30 l'étape 104 de la figure 3 de normalisation de contraste. A l'étape 305, le circuit logique de commande applique à chaque bloc de coupes l'algorithme de détection d'opacités de l'étape 105 de la figure 3. Après l'étape 305, le circuit logique de commande 10 applique les étapes 306 et 307 correspondantes respectivement aux étapes 106 et 107 35 de la figure 3.

17 Les étapes 301 à 305 sont répétées autant de fois que nécessaire pour toutes les tailles de la gamme de tailles d'opacités. Comme le contraste est calculé à l'étape 301, projection par projection, le contraste est calculé plus rapidement qu'à l'étape 202.

Typiquement avec neuf projections, le calcul du contraste est effectué neuf fois. Alors que selon la figure 3, il devrait être calculé autant de fois qu'il y a de taille d'opacité à tester, et, pour chaque taille autant de bloc de coupes qu'il est possible d'associer, sachant par ailleurs qu'il y a beaucoup de coupes (50 à 80). De même, selon la figure 4, il devrait être calculé autant de fois qu'il y a de coupes, sachant par ailleurs qu'il y a beaucoup de coupes (50 à 80). Le filtrage dans ces cas serait lui-même plus long car il porte sur des coupes ou des blocs de coupes. La figure 6 montre une deuxième stratégie de détection d'opacités, un quatrième mode de réalisation, utilisant la méthode de collecte de données 15 décrite à la figure 2. L'étape 400 de la figure 6 correspond à l'étape 100 de la figure 3 où le détecteur 5 acquiert les données images 11 à ln de rayons X, émis par le tube 4, représentant respectivement les projections P1 à Pn. A l'étape 401, le circuit logique de commande 10 applique un 20 algorithme de reconstruction par rétroprojection d'un volume numérique en coupe à partir des images brutes 11 à ln filtrées. A l'étape 402, le circuit logique de commande 10 applique à chaque coupe le moyen de calcul du contraste de l'étape 103 de la figure 3. L'étape 403 de la figure 6 correspond à l'étape 104 de la figure 3 où le circuit logique 25 de commande 10 normalise le contraste AI calculé en chaque voxel de chaque coupe. A l'étape 404, le circuit logique de commande effectue une association des coupes de contrastes normalisés en bloc de coupes, en fonction de la valeur du paramètre d'épaisseur. A l'étape 405, le circuit logique de 30 commande applique à chaque bloc de coupes l'algorithme de détection d'opacités de l'étape 105 de la figure 3. Après l'étape 405, le circuit logique de commande applique les étapes 406 et 40 à 447 correspondantes respectivement aux étapes 106 et 107 de la figure 3. 35 Les étapes 402 à 405 sont répétées autant de fois que nécessaires

18 pour toutes les tailles de la gamme de tailles d'opacités. Par rapport à la figure 3, la normalisation 403 est effectuée après le calcul de contraste. Dans une variante de cette deuxième stratégie, l'invention peut aussi mettre en oeuvre un mode de réalisation comportant un algorithme plus rapide, qui dépend lui aussi comme à la figure 5 seulement du nombre de projections et pas du nombre de coupes, en supposant une géométrie d'acquisition dans laquelle le facteur d'agrandissement varie peu en fonction de l'angulation du tube. Dans ce cas, à chaque image brute 11 à ln, le circuit logique de commande applique un moyen de calcul du contraste comme précédemment décrit à l'étape 301 de la figure 5. Puis le circuit logique de commande applique un algorithme de reconstruction d'un volume numérique en coupes à partir des images brutes 11 à ln filtrées. A partir de ces coupes filtrées le circuit logique de commande applique les étapes 403 à 407 de la figure 6.

Dans cette variante, les étapes de calcul de contraste sur les projections jusqu'à l'étape 405 de détection d'opacités sont répétées autant de fois que nécessaires pour toutes les tailles de la gamme de tailles d'opacités. La figure 7 montre une troisième stratégie, un cinquième mode de réalisation, du procédé de l'invention. L'étape 500 de la figure 7 correspond à l'étape 100 de la figure 3 où le détecteur 5 acquiert les données images 11 à In de rayons X, émis par le tube 4, représentant respectivement les projections P1 à Pn. A l'étape 501, le circuit logique de commande applique un algorithme 25 de reconstruction d'un volume numérique en coupes à partir des images brutes 11 à ln filtrées. A l'étape 502, le circuit logique de commande applique à chaque coupe le moyen de calcul du contraste de l'étape 103 de la figure 3. L'étape 503 de la figure 7 correspond à l'étape 104 de la figure 3. 30 A l'étape 504, le circuit logique de commande applique à chaque coupe reconstruite l'algorithme de détection d'opacités de l'étape 105 de la figure 3. Le circuit logique de commande attribue à chacun des voxels de la coupe correspondant à une opacité détectée un attribut représentatif d'opacité provisoire. La logique de commande crée un volume binaire où 35 tous les voxels sont placés à 0, excepté les voxels ayant un attribut d'opacité

19 provisoire qui sont placés à 1. A l'étape 505, le circuit logique de commande applique un algorithme de composantes connexes 3D au volume filtré et seuillé. L'algorithme de composantes connexes 3D permet d'extraire des ensembles de voxels reliés et dotés d'un attribut d'opacité provisoire. A l'étape 506, le circuit logique de commande calcule un contraste d'ensemble pour chaque ensemble de voxels ainsi reliés dans le volume et dotés de l'attribut d'opacité provisoire. Le circuit logique de commande affecte à chaque voxel d'un ensemble de voxels la valeur du contraste d'ensemble de cet ensemble. Lorsque la valeur de ce contraste d'ensemble est supérieure à un seuil prédéfini, le circuit logique de commande attribue à chaque voxel de l'ensemble un attribut d'opacité définitif. Après l'étape 506, le circuit logique de commande applique l'étape 507 correspondante à l'étape 107 de la figure 3. A cette étape, on visualise les voxels pourvus d'un attribut d'opacité définitif dans un volume numérique de présentation produit à partir des projections. Les étapes 501 à 506 sont répétées autant de fois que nécessaire pour toutes les tailles de la gamme de tailles d'opacités. Dans une variante de la stratégie illustrée à la figure 7, l'invention peut mettre en oeuvre un mode de réalisation comportant un algorithme plus rapide qui dépend seulement du nombre de projections et pas du nombre de coupes en supposant une géométrie d'acquisition dans laquelle le facteur d'agrandissement varie peu en fonction de l'angulation du tube. Dans ce cas, à chaque image brute 11 à ln, le circuit logique de commande 10 applique un moyen de calcul du contraste comme précédemment décrit à l'étape 301 de la figure 5. Puis le circuit logique de commande 10 applique un algorithme de reconstruction d'un volume numérique en coupes à partir des images brutes 11 à ln filtrées. A partir de ces coupes filtrées le circuit logique de commande applique les étapes 504 à 508 de la figure 7. La figure 8 montre une variante de la figure 7. Dans l'exemple de la figure 8, les étapes 600 à 605 correspondent respectivement aux étapes 500 à 505 de la figure 7. A l'étape 606, le circuit logique de commande repère, pour chaque voxel doté d'un attribut d'opacité provisoire, les coupes du volume numérique qui contiennent ce voxel. Le circuit logique de commande

20 détermine une coupe privilégiée C parmi les coupes repérées du volume pour chaque voxel ayant un attribut d'opacité provisoire. Il mesure une élongation L de chaque voxel doté de l'attribut d'opacité provisoire dans les coupes contenant ce voxel. Cette élongation est fournie par le plus grand nombre de Feret. Le nombre de Feret permet de mesurer l'étalement dudit voxel le long d'un ensemble de directions échantillonnant l'intégralité des directions possibles. Puis, il détermine pour chaque voxel doté de l'attribut d'opacité provisoire un bloc de coupes, ce bloc de coupes étant formé par les coupes allant de L/2 coupes précédant la coupe privilégiée C à L/2 coupes suivant la coupe privilégiée. Pour chaque voxel doté de l'attribut d'opacité provisoire, le circuit logique de commande calcule la valeur contraste d'ensemble à lui affecter. Dans un mode de réalisation préféré, la coupe privilégiée est la coupe dans laquelle l'élongation de l'ensemble de voxels, dotés de l'attribut d'opacité provisoire, est maximale. Dans une variante, la coupe privilégiée est la coupe comportant le voxel, de l'ensemble de voxels dotés de l'attribut d'opacité provisoire, ayant le plus fort contraste. L'étape 607 de visualisation de la figure 8 correspond à l'étape de visualisation 507 de la figure 7.

Cette troisième stratégie illustrée aux figures 7 et 8 est moins sensible aux artéfacts de reconstruction selon l'axe perpendiculaire au détecteur dus à un angle d'acquisition limité. La quatrième stratégie de détection est mise en oeuvre en appliquant directement au volume reconstruit un filtre tridimensionnel. Cette quatrième stratégie décrite à la figure 9 permet de calculer un indicateur de présence d'opacité directement sur le volume. La figure 8 montre un septième mode de réalisation du procédé de l'invention. L'étape 700 de la figure 9 correspond à l'étape 100 de la figure 3 où le détecteur 5 acquiert les données images 11 à ln de rayons X, émis par le tube 4, représentant respectivement les projections P1 à Pn.

A l'étape 701, le circuit logique de commande applique un algorithme de reconstruction d'un volume numérique en coupe à partir des images brutes 11 à ln. A l'étape 702, le circuit logique de commande applique à chaque coupe un filtre différentiel tridimensionnel dont le paramètre d'échelle dépend 35 d'une taille d'opacité à rechercher. Ce filtre différentiel tridimensionnel peut

21 être un filtre passe bande tridimensionnel ou un opérateur différentiel comme le Hessian. Ce filtre a un paramètre d'échelle défini en fonction de la taille de l'opacité à détecter. Le circuit logique de commande considère la réponse des filtres comme valeur de contraste en chaque voxel.

L'étape 703 de la figure 9 correspond à l'étape 104 de la figure 3 de normalisation des contrastes. A l'étape 704 le circuit logique de commande applique au volume filtré l'algorithme d'identification d'opacité tel que décrit à l'étape 105 de la figure 3. Les étapes 702 à 704 sont répétées autant de fois que nécessaire pour toutes les tailles de la gamme de tailles d'opacités. Après l'étape 704, le circuit logique de commande applique les étapes 705 et 706 correspondantes respectivement aux étapes 106 et 107 de la figure 3.

Claims (16)

REVENDICATIONS

1. - Procédé de traitement d'images produites par un dispositif (1) à rayons X, dans lequel, - on réalise plusieurs projections (19) radiographiques d'un corps, - on reconstruit (21) un volume numérique (V) en coupes du corps, à partir de l'ensemble des projections, caractérisé en ce que - on associe (22) des coupes (40-44) consécutives du corps en des blocs (50) de coupes en fonction d'un paramètre d'épaisseur (24) du bloc, - une épaisseur (60) étant caractérisée par une taille d'opacités à analyser, - on détermine (23) un nombre suffisant d'épaisseurs, pour couvrir une gamme de tailles des opacités, le nombre correspondant au nombre des tailles, - on applique à chaque bloc de coupes un algorithme (26) de détection d'opacités, - on donne à chacun des voxels du bloc correspondant à une opacité détectée un attribut représentatif d'opacité, - on visualise (35) les voxels ayant un attribut d'opacité dans un volume numérique de présentation produit à partir des projections.

2 - Procédé selon la revendication 1, caractérisé en ce que lors de l'association des coupes en blocs de coupes, les voxels du volume des coupes du bloc reçoivent un attribut résultant d'une opération mathématique qui peut être notamment une somme, un maximum, un minimum, une moyenne du niveau de gris mesuré pour chaque voxel du volume selon un axe perpendiculaire au détecteur ayant des même coordonnées (xi,yi) dans toutes les coupes du bloc l'attribut mesuré leur est ensuite affecté à tous.

3 - Procédé selon la revendication 2, caractérisé en ce que l'opération mathématique est une combinaison logique ou algébrique de niveaux de gris dans les coupes, suivant une direction d'émission des rayons X de la source au détecteur.

4 - Procédé selon l'une quelconque des revendications 1 à 3, caractérisé en ce que l'algorithme de détection d'opacités comporte les 35 étapes suivantes : 23 - on calcule une valeur de contraste (24) pour chaque voxel d'un bloc de coupes, - on définit préalablement un seuil de contraste, -lorsque la valeur de contraste du voxel est supérieure au seuil de contraste, on attribue audit voxel un attribut représentatif d'opacité.

5 - Procédé selon la revendication 4, caractérisé en ce que le calcul du contraste en chaque voxel de chaque bloc comporte les étapes suivantes : -on effectue une convolution (103) de chaque bloc de coupes avec un filtre de paramètre d'échelle s, - le paramètre d'échelle étant défini en fonction de la taille de l'opacité à détecter, - on détermine la valeur du contraste de chaque voxel du bloc à partir du bloc filtré.

6 - Procédé selon la revendication 4, caractérisé en ce que le calcul 15 du contraste en chaque voxel de chaque bloc comporte les étapes suivantes : - on effectue une convolution (202) de chaque coupe avec un filtre de paramètre d'échelle s, - le paramètre d'échelle étant défini en fonction de la taille de l'opacité 20 à détecter, - on associe (203) les coupes filtrées en des blocs de coupes filtrées, - on détermine la valeur du contraste de chaque voxel de chaque bloc à partir du bloc filtré.

7 - Procédé selon la revendication 4, caractérisé en ce que le calcul 25 du contraste de chaque voxel de chaque bloc comporte les étapes suivantes: - on effectue une convolution (301) de chaque projection avec un filtre de paramètre d'échelle s, - le paramètre d'échelle étant défini en fonction de la taille de l'opacité à détecter, 30 - on reconstruit (302) un volume numérique de coupes filtrées à partir des projections filtrées, - on associe (303) les coupes filtrées en des blocs de coupes filtrées, - on détermine la valeur du contraste de chaque voxel de chaque bloc à partir du bloc filtré._ 35

8 - Procédé de traitement d'images produites par un dispositif (1) à 24 rayons X, dans lequel, - on réalise plusieurs projections (19) radiographiques d'un corps, - on reconstruit (21) un volume numérique en coupes du corps, à partir de l'ensemble des projections, - on calcule une valeur de contraste pour chaque voxel, caractérisé en ce que - on applique un algorithme de détection (26) d'opacités à chaque coupe, - on donne à chacun des voxels de la coupe correspondant à une opacité détectée un attribut représentatif d'opacité provisoire, - on forme des ensembles de voxels (27) dotés de l'attribut d'opacité provisoire et reliés dans le volume, - pour un ensemble de voxels ainsi reliés dans le volume, on calcule un contraste d'ensemble (28), - ce contraste d'ensemble étant notamment une moyenne, une somme ou une valeur maximale des contrastes des voxels de l'ensemble, - on affecte à chaque voxel d'un ensemble de voxels le contraste d'ensemble de cet ensemble, - lorsque ce contraste d'ensemble est supérieur à un seuil prédéfini, on attribue à chaque voxel de l'ensemble un attribut d'opacité définitif (29), - on visualise (35) les voxels pourvus d'un attribut d'opacité définitif dans un volume numérique de présentation produit à partir des projections.

9 - Procédé selon la revendication 8, caractérisé en ce que le calcul du contraste d'ensemble comporte les étapes suivantes : - pour chaque voxel de l'ensemble, on repère les coupes du volume numérique qui contiennent ce voxel, - on détermine une coupe privilégiée C parmi les coupes repérées du volume, pour chaque voxel doté d'un attribut d'opacité provisoire, - on mesure une élongation L de l'intersection de l'ensemble avec la coupe privilégiée C, - on détermine un paquet de coupes, ce paquet de coupes étant formé par les coupes allant de L/2 coupes précédant la coupe privilégiée à L/2 coupes suivant la coupe privilégiée, - le contraste d'ensemble est notamment une valeur moyenne, une 35 somme ou une valeur maximale des contrastes des dits voxels de l'ensemble 25 situés dans les coupes du paquet de coupes.

10 - Procédé selon la revendication 9, caractérisé en ce que la coupe privilégiée est la coupe dans laquelle l'élongation de l'ensemble de voxels, est maximale.

11 - Procédé selon la revendication 9, caractérisé en ce que la coupe privilégiée est la coupe comportant le voxel de l'ensemble ayant le plus fort contraste.

12 - Procédé selon l'une quelconque des revendications 4 à 11, caractérisé en ce que le filtre de calcul de contraste est un filtre d'ondelette 10 de type chapeau mexicain de paramètre d'échelle s, - on détermine un nombre suffisant d'échelles s pour couvrir une gamme de tailles des opacités.

13 - Procédé selon l'une quelconque des revendications 1 à 12, caractérisé en ce que 15 - à titre de présentation, on attribue une information de couleur (33) aux voxels ayant un attribut d'opacité dans le volume numérique de présentation, ou - on rehausse (34) l'intensité des voxels ayant un attribut d'opacité dans le volume numérique de présentation, 20 -on visualise le volume numérique de présentation avec les voxels colorisés ou rehaussés.

14 - Procédé selon la revendication 13, caractérisé en ce que le rehaussement des voxels dans le volume numérique de présentation comporte les étapes suivantes : 25 - on détermine des particules de voxels ayant un attribut d'opacité et reliés dans le volume, - pour chaque particule, on amplifie la différence entre l'intensité de chaque voxel de la particule et l'intensité moyenne des voxels environnants ladite particule. 30

15 - Procédé de traitement d'images produites par un dispositif (1) à rayons X, dans lequel, - on réalise plusieurs projections (19) radiographiques d'un corps, - on reconstruit (21) un volume numérique en coupes du corps, à partir de l'ensemble des projections, 35 caractérisé en ce que- on calcule (31) le contraste de chaque voxel du volume en convoluant le volume numérique avec un banc de filtres (30) différentiels tridimensionnels ayant différentes échelles correspondantes à différentes tailles d'opacités, - on définit préalablement un seuil de contraste, -lorsque le contraste calculé pour le voxel est supérieur au seuil de contraste pour au moins une échelle, on attribue audit voxel un attribut représentatif d'opacité, - on visualise (35) les voxels ayant un attribut d'opacité dans un volume numérique de présentation produit à partir des projections.

16. - Dispositif (1) de rayons X caractérisé en ce qu'il comporte un microprocesseur (12) apte à mettre en oeuvre un procédé de traitement d'images selon l'une quelconque des revendications 1 à 15.