FR2932298A1 - Procede de traitement d'une image radiologique d'un organe - Google Patents
Procede de traitement d'une image radiologique d'un organe Download PDFInfo
- Publication number
- FR2932298A1 FR2932298A1 FR0853763A FR0853763A FR2932298A1 FR 2932298 A1 FR2932298 A1 FR 2932298A1 FR 0853763 A FR0853763 A FR 0853763A FR 0853763 A FR0853763 A FR 0853763A FR 2932298 A1 FR2932298 A1 FR 2932298A1
- Authority
- FR
- France
- Prior art keywords
- image
- recombined
- gray level
- recombined image
- energy
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000012545 processing Methods 0.000 title claims abstract description 17
- 210000000056 organ Anatomy 0.000 title claims abstract description 13
- 239000002872 contrast media Substances 0.000 claims abstract description 38
- 238000012937 correction Methods 0.000 claims abstract description 6
- 238000005215 recombination Methods 0.000 claims description 22
- 230000006798 recombination Effects 0.000 claims description 22
- 238000001228 spectrum Methods 0.000 claims description 15
- 230000006835 compression Effects 0.000 claims description 10
- 238000007906 compression Methods 0.000 claims description 10
- 239000000203 mixture Substances 0.000 claims description 7
- 238000013178 mathematical model Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 2
- 244000245420 ail Species 0.000 claims 1
- 235000004611 garlic Nutrition 0.000 claims 1
- 210000000481 breast Anatomy 0.000 description 49
- 230000005855 radiation Effects 0.000 description 21
- 210000001519 tissue Anatomy 0.000 description 16
- 230000007423 decrease Effects 0.000 description 8
- 238000009792 diffusion process Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 210000000038 chest Anatomy 0.000 description 4
- 230000003902 lesion Effects 0.000 description 4
- 238000009607 mammography Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000000295 complement effect Effects 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 230000000762 glandular Effects 0.000 description 2
- 229910052740 iodine Inorganic materials 0.000 description 2
- 239000011630 iodine Substances 0.000 description 2
- 230000000873 masking effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000002083 X-ray spectrum Methods 0.000 description 1
- 210000000577 adipose tissue Anatomy 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 229910052788 barium Inorganic materials 0.000 description 1
- DSAJWYNOEDNPEQ-UHFFFAOYSA-N barium atom Chemical compound [Ba] DSAJWYNOEDNPEQ-UHFFFAOYSA-N 0.000 description 1
- 210000000746 body region Anatomy 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000011282 treatment Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
L'invention concerne un procédé de traitement d'une image radiologique d'un organe comprenant des étapes de : - génération d'une image recombinée (R) à partir d'au moins une image basse énergie (L) et 'une image haute énergie (H) préalablement acquises d'une partie de l'organe, l'image recombinée comprenant d'une part une ou plusieurs zones correspondant à la présence d'un agent de contraste dans l'organe et d'autre part une zone d'arrière-plan, - correction de l'image recombinée (R) par compensation (C) des non-uniformités présentes dans la zone d'arrière-plan, de manière à obtenir une image recombinée compensée dans laquelle la zone d'arrière-plan présente un niveau de gris sensiblement uniforme.
Description
DOMAINE DE L'INVENTION L'invention concerne un procédé de traitement d'images 5 radiologiques, et en particulier d'images mammographiques.
ETAT DE LA TECHNIQUE Dans le domaine de la radiologie, et plus particulièrement en 10 mammographie, on connait des techniques d'acquisition d'images appelées double énergie . Ces techniques consistent à injecter un agent de contraste (par exemple à base d'iode) dans des tissus d'une région du corps à imager (par exemple un sein) puis à acquérir successivement deux images appelées image basse énergie et image haute énergie . Les 15 images basse et haute énergie sont acquises en exposant la région du corps à des rayonnements (rayons X) présentant des spectres ou gammes d'énergie différents. Comme illustré sur la figure 1, l'agent de contraste présente en effet un saut K d'énergie (référence 20) correspondant à un accroissement 20 soudain du coefficient d'atténuation des photons se produisant pour une énergie des photons égale à l'énergie de liaison des électrons de la couche K des atomes de l'agent de contraste interagissant avec les photons. Les agents de contraste couramment utilisés tels que l'iode et le baryum présentent par exemple des sauts K respectifs pour l'absorption des rayons 25 X respectivement à 33,2 keV et 37,4 KeV. L'image basse énergie est acquise en exposant la région à imager à un rayonnement dont l'énergie est inférieure au saut K de l'agent de contraste, tandis que l'image haute énergie est acquise en exposant la région à imager à un rayonnement dont l'énergie est supérieure au saut K de l'agent de contraste. 30 Dans l'image basse énergie, l'atténuation des rayonnements par l'agent de contraste injecté est relativement faible et l'image obtenue montre un fort contraste entre les tissus adipeux et glandulaires. En comparaison, l'image haute énergie montre une atténuation plus importante des rayonnements par l'agent de contraste de sorte que l'image haute énergie présente un plus grand contraste entre l'agent de contraste et les tissus environnants. Afin de permettre une meilleure visualisation des zones de l'organe contenant de l'agent de contraste, les images basse et haute énergies sont recombinées pour produire une image fonctionnelle recombinée faisant ressortir les zones de tissus marquées par l'agent de contraste. L'opération de recombinaison vise à obtenir une image fonctionnelle recombinée sur laquelle apparaissent d'une part une ou plusieurs zones dans lesquelles l'agent de contraste est présent, notamment des zones de lésion, et d'autre part une zone d'arrière-plan présentant un niveau de gris uniforme. En mammographie numérique à contraste amélioré double énergie (Dual-energy contrast-enhanced digital mammography û DE CEDM) ou angio-mammographie numérique, le sein d'une patiente est comprimé afin de présenter une épaisseur la plus constante possible. Puis, deux images sont acquises successivement, en soumettant la zone du corps à des rayonnements présentant des spectres ou gammes d'énergie différents. Les images basse énergie et haute énergie obtenues sont recombinées pour obtenir une image fonctionnelle recombinée de l'agent de contraste. De manière similaire, en tomosynthèse du sein numérique à contraste amélioré double énergie (Dual energy contrast-enhanced digital breast tomosynthesis û DE CEDBT) ou angio-tomosynthèse numérique du sein, le sein d'une patiente est comprimé, et des séquences d'images projectives basse énergie et haute énergie sont acquises en déplaçant la source de rayonnement en des positions différentes. Cette opération peut être réalisée en acquérant successivement une image projective basse énergie et une image projective haute énergie à chaque position de la source, ou en acquérant en premier lieu une séquence d'images projectives basse énergie pour différentes positions de la source puis en second lieu une séquence d'images projectives haute énergie pour différentes positions de la source. Pour chaque position angulaire de la source, les images basse énergie et haute énergie correspondantes sont recombinées de manière à obtenir une séquence d'images recombinées. L'opération de recombinaison peut consister en une soustraction sous la forme : xs =log(xh)ûR•log(xi) où xs est le niveau de gris d'un pixel de l'image recombinée, xh est le niveau de gris d'un pixel de l'image haute énergie, xi est le niveau de gris d'un pixel de l'image basse énergie. Dans l'hypothèse où les spectres d'énergie des rayonnements haute et basse énergie utilisés pour l'acquisition des images sont mono-énergétiques, l'ajustement du paramètre R permet en théorie de supprimer (ou soustraire) les données associés aux tissus de l'organe et de ne conserver que les données relatives à l'agent de contraste. Cependant, en pratique, cette technique de recombinaison ne convient pas lorsque les spectres d'énergie des rayonnements utilisés ne sont pas mono-énergétiques. En particulier, la soustraction logarithmique conduit à l'obtention d'une image recombinée présentant dans la zone d'arrière-plan des non-uniformités dues à la répartition des tissus glandulaires et adipeux dans le sein.
Par ailleurs, d'autres facteurs sont également la cause de variations de niveaux de gris apparaissant dans l'image recombinée. En effet, malgré une mise en compression, le sein ne présente pas une épaisseur uniforme. Plus précisément, le sein comprend une partie comprimée (partie intérieure) présentant une épaisseur sensiblement uniforme et une partie non comprimée (partie de bordure), dans laquelle l'épaisseur du sein décroît très rapidement. Les variations d'épaisseur du sein ainsi que les variations de composante de diffusion des tissus à proximité de la bordure du sein créent des variations de niveaux de gris apparaissant dans la zone d'arrière-plan de l'image recombinée. La décroissance d'épaisseur du sein à proximité de la bordure du sein a pour effet d'accroître les valeurs de niveaux de gris dans les images basse énergie et haute énergie lorsqu'on approche la bordure du sein. Cette augmentation des niveaux de gris résulte en une décroissance des valeurs de niveaux de gris dans l'image recombinée. Comme l'image recombinée est une image où les valeurs de niveaux de gris sont proportionnelles à l'épaisseur d'agent de contraste, cette décroissance est interprétée comme une diminution de la présence de l'agent de contraste à proximité de la bordure du sein. A l'inverse, la décroissance de la composante de diffusion à proximité de la bordure du sein a pour effet de diminuer les niveaux de gris dans les images basse énergie et haute énergie lorsqu'on approche la bordure du sein. Cette diminution des niveaux de gris est interprétée dans l'image recombinée comme une augmentation de la présence de l'agent de contraste à proximité de la bordure du sein. D'autres facteurs, tels qu'un effet talon de l'anode de la source du dispositif d'acquisition, qui créent des variations de distribution spectrale dans le faisceau de rayonnement produit par la source, peuvent également causer des variations de niveaux de gris dans la zone d'arrière-plan de l'image recombinée. Ces variations de niveau de gris sont gênantes pour le radiologue car elles rendent difficile la visualisation de la partie de bordure du sein sans changer le niveau de fenêtrage utilisé pour visualiser l'intérieur du sein (i.e. la zone comprimée du sein). En outre, ces variations de niveau de gris risquent d'être interprétées comme des signes de la présence d'agent de contraste et de 25 conduire à des erreurs d'interprétation de l'image.
PRESENTATION DE L'INVENTION Un but de l'invention est de produire une image recombinée 30 présentant une qualité améliorée, c'est-à-dire dans laquelle les non-uniformités de la zone d'arrière-plan ont été atténuées, voire éliminées.
Ce problème est résolu dans le cadre de la présente invention grâce à un procédé de traitement d'une image radiologique d'un organe comprenant des étapes de : - génération d'une image recombinée à partir d'au moins une image basse énergie et une image haute énergie préalablement acquises d'une partie de l'organe, l'image recombinée comprenant d'une part une ou plusieurs zones correspondant à la présence d'un agent de contraste dans l'organe et d'autre part une zone d'arrière-plan, - correction de l'image recombinée par compensation des non-10 uniformités présentes dans la zone d'arrière-plan, de manière à obtenir une image recombinée compensée dans laquelle la zone d'arrière-plan présente un niveau de gris sensiblement uniforme. Selon ce procédé, l'étape de correction est réalisée après l'étape de 15 recombinaison, ce qui permet d'éliminer les non-uniformités dues aux variations de diffusion qui apparaissent du fait de la recombinaison des images basse et haute énergie. On notera que ce résultat ne pourrait être obtenu par une compensation individuelle opérée dans chacune des images basse et haute 20 énergie, avant l'étape de recombinaison. Le procédé permet d'obtenir une image recombinée présentant une bonne qualité esthétique, et permet potentiellement de mieux détecter d'éventuelles lésions et de mieux quantifier la répartition de l'agent de contraste. 25 De préférence, l'étape de correction de l'image recombinée comprend l'application à l'image recombinée d'une image de compensation pour compenser les non-uniformités présentes dans la zone d'arrière-plan. Cette application consiste en une addition ou une multiplication pixel à pixel de l'image recombinée et de l'image de compensation. 30 Cependant, l'étape de correction de l'image recombinée peut comprendre l'application d'autres traitements de compensation, tels que par exemple l'application d'un filtre linéaire ou non linéaire permettant d'éliminer les non-uniformités présentes dans la zone d'arrière-plan.
Dans une mise en oeuvre de l'invention, le procédé comprend une étape de : - génération d'une image recombinée filtrée par extraction de composantes spatiales basses fréquences de l'image recombinée, - génération de l'image de compensation à partir de l'image recombinée filtrée. Le filtrage des basses fréquences découle du constat que les non-uniformités dues aux variations de diffusion dans l'image recombiné correspondent à des composantes basses fréquences tandis que les lésions, dans lesquelles l'agent de contraste s'est accumulé, correspondent à des composantes de plus hautes fréquences. Plus précisément, l'étape de génération d'une image de compensation comprend une sous-étape de : - détermination d'un niveau de gris d'arrière-plan à partir d'un 15 histogramme de l'image recombinée, - calcul de l'image de compensation comme la différence entre une image uniforme présentant le niveau de gris déterminé et l'image compensée filtrée. L'image de compensation est ainsi calculée comme l'image 20 complémentaire de l'image recombinée pour obtenir un niveau de gris uniforme dans tout le sein. La sous-étape de détermination d'un niveau de gris d'arrière-plan comprend la détection à partir d'un histogramme de l'image recombinée d'un pic de pixels présentant le niveau de gris. 25 En outre, dans un mode de réalisation de l'invention, l'étape de génération d'une image recombinée comprend une opération de recombinaison sous une forme polynomiale : y = Lai, j . p(xl )' .(p(xh )' où y est le niveau de gris d'un pixel de l'image recombinée (R), 30 xi est le niveau de gris d'un pixel de l'image basse énergie (L), xh est le niveau de gris d'un pixel de l'image haute énergie (H ), cp est une fonction associant à des valeurs de niveau de gris acquises une épaisseur radiologique, i et j sont des entiers, et ai,] sont des paramètres prédéterminés de recombinaison.
Les paramètres de recombinaison ai,] peuvent être des paramètres prédéterminés, qui dépendent de conditions d'acquisition des images basse énergie et haute énergie. En particulier, les paramètres de recombinaison dépendent d'une épaisseur de compression de l'organe imagé et de spectres énergétiques de rayonnements basse et haute énergie auxquels a été soumis l'organe lors de l'acquisition. En outre, les paramètres de recombinaison sont déterminés au préalable à partir de points de référence, les points de référence étant obtenus expérimentalement par calibrage d'un dispositif d'acquisition à l'aide d'un fantôme simulant plusieurs compositions de tissus et plusieurs épaisseurs d'agent de contraste, ou par calcul à partir d'un modèle mathématique, en faisant varier une composition de tissu et une épaisseur de l'agent de contraste, et en maintenant les conditions d'acquisition constantes.
De préférence, l'étape de génération d'une image recombinée comprend une opération de recombinaison sous une forme polynomiale du second degré du type : y = a00 +a1,0 . p(xj)+a0,1 .(p(xh)+al,1 .(p(xI).(p(xh)+a2,0 .(p(xi )2 +a0,2 .(P(xh )2 où la fonction cp est une fonction logarithmique. 25
PRESENTATION DES FIGURES D'autres caractéristiques et avantages ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative, et doit 30 être lue en regard des figures annexées, parmi lesquelles : - la figure 1 déjà commentée est un diagramme représentant de manière schématique un coefficient d'absorption d'un agent de contraste en fonction de l'énergie du rayonnement, - la figure 2 représente de manière schématique un dispositif 5 d'acquisition et de traitement d'image mammographique conforme à un mode de réalisation de l'invention, - la figure 3 est un diagramme représentant de manière schématique un niveau de gris dans une image recombinée d'un sein mis en compression en fonction de la distance à la cage thoracique, 10 - la figure 4 est un histogramme représentant de manière schématique les niveaux de gris observés dans une image recombinée après segmentation, - la figure 5 représente de manière schématique l'évolution des diagrammes des images au cours de différentes étapes de traitement de 15 l'image recombinée, - la figure 6 représente de manière schématique des étapes d'un procédé de traitement d'image conforme à un mode de mise en oeuvre de l'invention, - la figure 7 représente de manière schématique une table de points 20 de référence de valeurs de niveau de gris, établie selon une échelle logarithmique, pour des conditions d'acquisitions données des images basse et haute énergies.
25 DESCRIPTION DETAILLEE D'UN MODE DE REALISATION La figure 2 représente de manière schématique un dispositif 10 d'acquisition et de traitement d'image mammographique. Le dispositif 10 comprend une source 11 apte à émettre des rayons X 18, un détecteur 12 numérique apte à recevoir et détecter les rayons émis 30 par la source 11, et une unité de traitement 13 apte à commander la source 11 et à recevoir et traiter des images acquises par le détecteur 12. Un sein 15 d'une patiente est disposé entre la source 11 et le détecteur 12.
Le dispositif 10 comprend également une plaque de compression 14 pour mettre en compression le sein 15 afin que le sein présente une partie comprimée d'épaisseur sensiblement constante. Le sein 15 a reçu préalablement une injection d'un agent de contraste. Le sein 15 comprend une zone 16 de tissus dépourvue d'agent de contraste et une ou plusieurs zone(s) 17 de tissus dans lesquelles l'agent de contraste s'est accumulé. L'unité de traitement 13 est apte à commander la source 11 pour faire varier le spectre d'énergie du rayonnement 18 émis par la source 11.
En particulier, l'unité de traitement 13 est apte à commander la source 11 pour qu'elle émette un rayonnement basse énergie, dont l'énergie moyenne est inférieure au saut K de l'agent de contraste ou un rayonnement haute énergie, dont l'énergie moyenne est supérieure au saut K de l'agent de contraste.
Le détecteur 12 est apte à recevoir le rayonnement émis par la source, après que celui-ci ait traversé les tissus du sein 15. Le détecteur 12 génère une image sous la forme de données comprenant un ensemble de pixels et pour chaque pixel un niveau de gris associé. Le niveau de gris est représentatif de l'atténuation des rayons X ayant traversé les tissus du sein.
Le détecteur 12 est apte à transmettre à l'unité de traitement 13 l'image ainsi acquise. L'unité de traitement 13 est programmée pour recevoir les images basse énergie et haute énergie acquises et exécuter un procédé de traitement pour générer une image recombinée à partir des images basse énergie et haute énergie acquises. L'unité de traitement 13 peut également être programmée pour commander l'affichage sur un écran de l'image recombinée générée, afin qu'un radiologue puisse visualiser l'image recombinée. On remarque sur la figure 2 que le sein 15 comprend une partie comprimée (partie intérieure) présentant une épaisseur e sensiblement constante et une partie non comprimée (partie de bordure, la plus éloignée de la cage thoracique), dans laquelle l'épaisseur du sein décroît très rapidement.
La figure 6 représente de manière schématique des étapes d'un procédé de traitement d'image mis en oeuvre par l'unité de traitement 13. Une image basse énergie L et une image haute énergie H du sein de la patiente ont été préalablement acquises au cours d'un procédé d'acquisition.
En outre, on suppose que l'on connait l'épaisseur e de compression du sein et les spectres des rayonnements basse énergie et haute énergie émis par la source lors de l'acquisition des images L et H.
L'unité de traitement 13 exécute les étapes suivantes : Selon une première étape 21, l'image basse énergie L et l'image haute énergie H sont recombinées pour obtenir une image composite recombinée R .
L'image recombinée R est calculée sous la forme :
y = Lai, j . p(xl )' .(p(xh )' i, j où y est le niveau de gris d'un pixel de l'image composite recombinée R , représentatif d'une épaisseur de l'agent de contraste, xi est le niveau de gris d'un pixel de l'image basse énergie L, xh est le niveau de gris d'un pixel de l'image haute énergie H, cp est une fonction associant à des valeurs de niveau de gris acquises une épaisseur radiologique, i et j sont des entiers. Les paramètres ai,] de recombinaison sont des paramètres prédéterminés, associés aux conditions de l'acquisition, telles que notamment l'épaisseur de compression du sein et les spectres énergétiques 25 des rayonnements basse et haute énergie émis par la source. L'unité de traitement sélectionne une série de paramètres ai,] de recombinaison en fonction de l'épaisseur e de compression du sein, des spectres énergétiques des rayonnements basse et haute énergie émis par la source, ainsi que d'un paramètre de correction de l'épaisseur de 30 compression 6 f .
Le paramètre 6 f est obtenu par analyse de l'histogramme de l'image basse énergie. Un procédé d'obtention et d'application d'un paramètre de correction d'épaisseur est décrit par exemple dans le document EP 1 113 392 Al.
Les paramètres ai ,j de recombinaison ont été préalablement déterminés à partir d'une série de points de référence (y,xi,xh) similaire à celle de la figure 7. Les points de référence (y,xi,xh) peuvent être déterminés soit expérimentalement par calibrage du dispositif d'acquisition à l'aide d'un fantôme simulant plusieurs compositions de tissus et plusieurs épaisseurs d'agent de contraste, soit par calcul à partir d'un modèle mathématique. La série de points de référence (y,xi,xh) est établie en faisant varier une composition de tissu et une épaisseur de l'agent de contraste, et en maintenant les paramètres d'acquisition et l'épaisseur du tissu constants.
Dans le cas où un modèle mathématique est utilisé, le modèle simule la génération d'un spectre d'énergie de rayons X étant donné un potentiel (kVp) et des valeurs de paramètres représentatives de la composition matérielle de la source de rayonnement. Par exemple, on suppose que les données dans la table de la figure 7 soient générées en fonction des valeurs suivantes de spectre d'énergie : Mo/Mo 25kV, 100 mAs pour l'acquisition basse énergie et Mo/Cu 49 kV, 160 mAs pour l'acquisition haute énergie. Le modèle simule la génération du spectre de rayons X en recevant en entrée ou en calculant un nombre de photons générés dans les spectres basse et haute énergie. Le modèle simule également l'atténuation des spectres d'énergie des rayons X à travers divers tissus d'épaisseurs variées de la zone imagée (par exemple en supposant une épaisseur de 40 mm), et simule la transformation du spectre d'énergie des rayons X en valeur de niveau de gris des pixels détectés par le détecteur. En utilisant ce modèle mathématique de simulation, les points de référence (y,xi,xh) de la figure 7 peuvent être déterminés.
Les coefficients ai,j de recombinaison sont déterminés au moyen d'une analyse de régression linéaire de la série de points de référence (Y,xl,xh)• Dans une mise en oeuvre préférée de l'invention, l'image 5 recombinée R est calculée sous la forme d'une équation polynomiale du second ordre : y = a0 0 + a1,0 . (p (x/) + a0,1 . (p (xh) + al,1 . (p (xl) . (p (xh) + a2,0 . (P (xl )2 + a0,2 . (P (xh )2 La fonction cp est typiquement une fonction logarithmique. La partie ao 0 + a1,o . cp(xl )+ ao,1 . cp(xh) de cette équation mathématique représente de 10 manière générale un modèle mathématique pour une soustraction logarithmique. D'autres équations polynomiales d'ordre supérieur pourraient également être utilisées. La figure 3 représente de manière schématique un niveau de gris dans l'image recombinée du sein en fonction de la distance à la cage 15 thoracique. On observe que l'image recombinée présente un niveau de gris sensiblement constant dans la partie comprimée du sein la plus proche de la cage thoracique et un niveau de gris qui décroît très rapidement dans la partie de bordure du sein à mesure que l'on s'éloigne de la cage thoracique. 20 En outre, la partie comprimée présente des zones localisées de niveau de gris maximal correspondant à la présence d'agent de contraste. Par ailleurs, on observe également une augmentation de niveau de gris au niveau de la transition entre la partie comprimée et la partie non-comprimée du sein. Cet artéfact résulte d'un effet de coupelle ( cupping artefact ) due à une 25 modification de la diffusion du rayonnement dans la partie de bordure du sein. Selon une deuxième étape 22, l'image recombinée R est segmentée pour ne conserver qu'une zone de l'image dans laquelle apparaît le sein.
A cet effet, une image de masquage est déterminée à partir de l'image basse énergie, puis l'image de masquage est appliquée à l'image recombinée. L'image recombinée segmentée résultante ne contient qu'une zone 5 où apparait l'intérieur du sein. La figure 4 est un histogramme de l'image recombinée segmentée représentant le nombre de pixels dans l'image recombinée segmentée en fonction du niveau de gris. On observe un premier pic pour un niveau de gris nul, correspondant aux pixels de l'image situés en dehors du sein et un 10 deuxième pic pour un niveau de gris G correspondant au niveau de gris de la zone comprimée du sein ne contenant pas d'agent de contraste. Selon une troisième étape 23, l'image recombinée segmentée est filtrée au moyen d'un filtre passe-bas apte à sélectionner les composantes spatiales basses fréquences de l'image. L'image basse fréquence obtenue 15 PB ne contient que les variations progressives d'épaisseur. Cette étape est basée sur l'hypothèse que les variations de niveau de gris dues aux variations de diffusion correspondent à des composantes basses fréquences de l'image tandis que les lésions, dans lesquelles l'agent de contraste s'est accumulé, correspondent à des composantes de plus 20 hautes fréquences. Par conséquent, l'image recombinée filtrée PB contient uniquement les variations de niveaux de gris dues aux variations d'épaisseur et de diffusé dans le sein. Selon une quatrième étape 24, une image de compensation C est générée. 25 L'image de compensation est calculée comme la différence entre une image présentant un niveau de gris uniforme G et l'image recombinée filtrée PB. C=GûPB où G est le niveau de gris associé au plus grand nombre de pixels de 30 l'image recombinée. Le niveau de gris G est déterminé à partir de l'histogramme de l'image recombinée représenté sur la figure 4.
L'image de compensation C est ainsi calculée comme l'image complémentaire de l'image recombinée filtrée PB pour obtenir un niveau de gris égal à G dans tout le sein. Selon une cinquième étape 25, l'image de compensation C est appliquée à l'image recombinée R . Plus précisément, l'image de compensation est additionnée à l'image recombinée, de manière à obtenir une image finale recombinée compensée : Image finale compensée = R+ C La figure 5 montre l'évolution des diagrammes des images au cours 10 de différentes étapes de traitement de l'image recombinée. L'image finale recombinée compensée comprend une ou plusieurs zones de niveau de gris maximal correspondant à la présence d'un agent de contraste et une zone d'arrière-plan présentant un niveau de gris G uniforme, y compris à proximité de la bordure du sein. L'effet de diffusé 15 présent dans l'image recombinée initiale a été supprimé dans l'image recombinée finale.
Claims (11)
- REVENDICATIONS1. Procédé de traitement d'une image radiologique d'un organe (15) comprenant des étapes de : - génération d'une image recombinée (R) à partir d'au moins une image basse énergie (L) et une image haute énergie (H) préalablement acquises d'une partie de l'organe, l'image recombinée comprenant d'une part une ou plusieurs zones correspondant à la présence d'un agent de contraste dans l'organe et d'autre part une zone d'arrière-plan, - correction de l'image recombinée (R) par compensation (C) des non-uniformités présentes dans la zone d'arrière-plan, de manière à obtenir une image recombinée compensée dans laquelle la zone d'arrière-plan présente un niveau de gris sensiblement uniforme.
- 2. Procédé selon la revendication 1, dans lequel l'étape de correction de l'image recombinée comprend l'application à l'image recombinée d'une image de compensation pour compenser les non-uniformités présentes dans la zone d'arrière-plan.
- 3. Procédé selon l'une des revendications 1 ou 2, comprenant une étape de : - génération d'une image recombinée filtrée (PB) par extraction de composantes spatiales basses fréquences de l'image recombinée (R ), - génération de l'image de compensation (C) à partir de l'image recombinée filtrée(PB).
- 4. Procédé selon la revendication 3, dans lequel l'étape de génération d'une image de compensation (C) comprend une sous-étape 30 de : - détermination d'un niveau de gris (G) d'arrière-plan à partir d'un histogramme de l'image recombinée (R), - calcul de l'image de compensation (C) comme la différence entre une image uniforme présentant le niveau de gris G et l'image recombinée filtrée (PB ).
- 5. Procédé selon la revendication 4, dans lequel la sous-étape de détermination d'un niveau de gris (G) d'arrière-plan comprend la détection à partir d'un histogramme de l'image recombinée (R) d'un pic de pixels présentant le niveau de gris (G ).
- 6. Procédé selon l'une des revendications qui précèdent, dans lequel l'étape de génération d'une image recombinée (R) comprend une opération de recombinaison sous la forme : y = Lai, j . p(xl )' .(p(xh )' où y est le niveau de gris d'un pixel de l'image recombinée (R), xi est le niveau de gris d'un pixel de l'image basse énergie (L), xh est le niveau de gris d'un pixel de l'image haute énergie (H), cp est une fonction associant à des valeurs de niveau de gris acquises une épaisseur radiologique, i et j sont des entiers, et ai,] sont des paramètres prédéterminés de recombinaison.
- 7. Procédé selon la revendication 6, dans lequel les paramètres (ai ,j) de recombinaison sont des paramètres prédéterminés, qui dépendent de conditions d'acquisition des images basse énergie (L) et haute énergie 25 (H).
- 8. Procédé selon la revendication 7, dans lequel les paramètres (ai J) de recombinaison dépendent d'une épaisseur (e) de compression de l'organe imagé et de spectres énergétiques de rayonnements basse et 30 haute énergie auxquels a été soumis l'organe lors de l'acquisition.
- 9. Procédé selon l'une des revendications 6 à 8, dans lequel les paramètres de recombinaison (ai J) sont déterminés au préalable à partir de points de référence (y,xl,xh), les points de référence étant obtenus expérimentalement par calibrage d'un dispositif d'acquisition à l'aide d'un fantôme simulant plusieurs compositions de tissus et plusieurs épaisseurs d'agent de contraste, ou par calcul à partir d'un modèle mathématique, en faisant varier une composition de tissu et une épaisseur de l'agent de contraste, et en maintenant les conditions d'acquisition constantes.
- 10. Procédé selon l'une des revendications 6 à 9, dans lequel l'étape de génération d'une image recombinée (R) comprend une opération de recombinaison sous la forme : y = a0 0 + a1,0 . (p(x/) + a0,1 . (p(xh) + ail . (p(xl) . (p(xh )+ a2,0 . (P(xl )2 + a0,2 . (P(xh )2
- 11. Procédé selon l'une des revendications 6 à 10, dans lequel la fonction cp est une fonction logarithmique.15
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0853763A FR2932298B1 (fr) | 2008-06-06 | 2008-06-06 | Procede de traitement d'une image radiologique d'un organe |
US12/478,943 US8165379B2 (en) | 2008-06-06 | 2009-06-05 | Method of processing radiological images, and, in particular, mammographic images |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0853763A FR2932298B1 (fr) | 2008-06-06 | 2008-06-06 | Procede de traitement d'une image radiologique d'un organe |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2932298A1 true FR2932298A1 (fr) | 2009-12-11 |
FR2932298B1 FR2932298B1 (fr) | 2010-07-30 |
Family
ID=40259804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0853763A Active FR2932298B1 (fr) | 2008-06-06 | 2008-06-06 | Procede de traitement d'une image radiologique d'un organe |
Country Status (2)
Country | Link |
---|---|
US (1) | US8165379B2 (fr) |
FR (1) | FR2932298B1 (fr) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2967888A1 (fr) * | 2010-11-26 | 2012-06-01 | Gen Electric | Procede de galactographie et mammographe pour l'execution dudit procede |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102010041920A1 (de) * | 2010-10-04 | 2012-04-05 | Siemens Aktiengesellschaft | Verfahren zur Darstellung einer Konzentration eines Kontrastmittels in einem vorbestimmten Volumenabschnitt mittels Tomosynthese und entsprechendes Tomosynthesegerät |
RU2495623C1 (ru) * | 2012-03-11 | 2013-10-20 | Федеральное государственное бюджетное учреждение науки Институт космических исследований Российской академии наук (ИКИ РАН) | Способ двухэнергетической делительно-разностной маммографии |
JP6251297B2 (ja) | 2013-03-15 | 2017-12-20 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | デュアルエネルギー画像からの残差モード画像の決定 |
KR20160089688A (ko) * | 2015-01-20 | 2016-07-28 | 삼성전자주식회사 | 엑스선 영상 장치, 및 그 제어방법 |
DE102016204828B4 (de) * | 2016-03-23 | 2022-09-15 | Siemens Healthcare Gmbh | Erzeugen einer kombinierten Kontrastmittel- und Gefäßdarstellung von zu untersuchendem Brustgewebe |
JP6667462B2 (ja) * | 2017-02-21 | 2020-03-18 | 富士フイルム株式会社 | エネルギーサブトラクション処理装置、方法およびプログラム |
WO2020059306A1 (fr) * | 2018-09-18 | 2020-03-26 | 富士フイルム株式会社 | Dispositif, méthode et programme de traitement d'images |
EP3662837A1 (fr) | 2019-03-29 | 2020-06-10 | Siemens Healthcare GmbH | Procédé de fourniture de données d'image d'un organe creux |
CN111904446A (zh) * | 2020-09-11 | 2020-11-10 | 深圳先进技术研究院 | 乳腺成像系统及其成像方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4837686A (en) * | 1985-11-15 | 1989-06-06 | Picker International | Substance quantification in animal bodies |
US6125166A (en) * | 1998-01-13 | 2000-09-26 | Fuji Photo Film Co., Ltd. | Method of forming energy subtraction images |
US6173034B1 (en) * | 1999-01-25 | 2001-01-09 | Advanced Optical Technologies, Inc. | Method for improved breast x-ray imaging |
EP1113392A1 (fr) * | 1999-12-28 | 2001-07-04 | Ge Medical Systems Sa | Procédé de compensation de l'épaisseur d'un organe |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5046118A (en) * | 1990-02-06 | 1991-09-03 | Eastman Kodak Company | Tone-scale generation method and apparatus for digital x-ray images |
US6577762B1 (en) * | 1999-10-26 | 2003-06-10 | Xerox Corporation | Background surface thresholding |
US6683934B1 (en) * | 2000-06-05 | 2004-01-27 | General Electric Company | Dual energy x-ray imaging system and method for radiography and mammography |
US20080167552A1 (en) | 2007-01-04 | 2008-07-10 | General Electric Company | System and method of generating an image of a contrast agent injected into an imaged subject |
US9451924B2 (en) * | 2009-12-30 | 2016-09-27 | General Electric Company | Single screen multi-modality imaging displays |
-
2008
- 2008-06-06 FR FR0853763A patent/FR2932298B1/fr active Active
-
2009
- 2009-06-05 US US12/478,943 patent/US8165379B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4837686A (en) * | 1985-11-15 | 1989-06-06 | Picker International | Substance quantification in animal bodies |
US6125166A (en) * | 1998-01-13 | 2000-09-26 | Fuji Photo Film Co., Ltd. | Method of forming energy subtraction images |
US6173034B1 (en) * | 1999-01-25 | 2001-01-09 | Advanced Optical Technologies, Inc. | Method for improved breast x-ray imaging |
EP1113392A1 (fr) * | 1999-12-28 | 2001-07-04 | Ge Medical Systems Sa | Procédé de compensation de l'épaisseur d'un organe |
Non-Patent Citations (3)
Title |
---|
KAPPADATH S C ET AL: "Dual-energy digital mammography for calcification imaging: improvement by post-image processing", PROCEEDINGS OF THE SPIE - THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING SPIE - THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING USA, vol. 5745, no. 1, 20 April 2005 (2005-04-20), pages 1342 - 1350, XP002512149, ISSN: 0277-786X * |
KAPPADATH S CHEENU ET AL: "Dual-energy digital mammography for calcification imaging: noise reduction techniques.", PHYSICS IN MEDICINE AND BIOLOGY 7 OCT 2008, vol. 53, no. 19, 7 October 2008 (2008-10-07), pages 5421 - 5443, XP002512150, ISSN: 0031-9155 * |
KAPPADATH S ET AL: "Dual-energy digital mammography for calcification imaging: Scatter and nonuniformity corrections", MEDICAL PHYSICS, AIP, MELVILLE, NY, US, vol. 32, no. 11, 21 October 2005 (2005-10-21), pages 3395 - 3408, XP012075202, ISSN: 0094-2405 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2967888A1 (fr) * | 2010-11-26 | 2012-06-01 | Gen Electric | Procede de galactographie et mammographe pour l'execution dudit procede |
Also Published As
Publication number | Publication date |
---|---|
FR2932298B1 (fr) | 2010-07-30 |
US8165379B2 (en) | 2012-04-24 |
US20090304253A1 (en) | 2009-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
FR2932298A1 (fr) | Procede de traitement d'une image radiologique d'un organe | |
JP6999576B2 (ja) | 空間とスペクトル情報に基づく複数エネルギーのct画像におけるノイズ制御のためのシステムと方法 | |
CN107530040B (zh) | X射线ct装置、重构运算装置以及x射线ct图像生成方法 | |
CN102254310B (zh) | 用于减少ct图像数据中的图像伪影、尤其是金属伪影的方法 | |
US6891918B2 (en) | Methods and apparatus for acquiring perfusion data | |
US8340373B2 (en) | Quantitative image reconstruction method and system | |
FR2948481A1 (fr) | Procede d'imagerie pour la realisation d'une modelisation en triple energie, et dispositif pour la mise en oeuvre d'un tel procede | |
US9888897B2 (en) | High dynamic range mammography using a restricted dynamic range FFDM | |
RU2434288C1 (ru) | Способ коррекции цифровых изображений | |
KR102003042B1 (ko) | 멀티 에너지 엑스선에 기초한 의료 영상을 합성 및 표시하는 의료 영상 처리 장치 및 방법 | |
FR2803069A1 (fr) | Procede et systeme de compensation de l'epaisseur d'un organe | |
JPH0263439A (ja) | X線画像サブトラクションの方法 | |
US20140219423A1 (en) | Calibration free dual energy radiography method | |
JP7199813B2 (ja) | 放射線撮影装置 | |
CN110636796B (zh) | X射线暗场成像中的射束硬化校正 | |
KR20160089976A (ko) | 엑스선 영상장치 및 엑스선 영상생성방법 | |
US7471768B2 (en) | Systems and methods for estimating presence of a material within a volume of interest using x-ray | |
CN110458762B (zh) | 一种基于可调双因子的ct图像射束硬化伪影校正系统 | |
CN111899188A (zh) | 一种神经网络学习的锥束ct噪声估计与抑制方法 | |
KR101905745B1 (ko) | k 에지 기반의 내부 토모그라피 영상 처리를 위한 장치 및 방법 | |
DE102012222714A1 (de) | Ermittlung eines Mehrfachenergie-Bildes | |
FR2814666A1 (fr) | Procede et appareil d'examen d'un sein par injection d'un produit de contraste | |
Davis et al. | A modelling approach to beam hardening correction | |
Brun et al. | K-edge spectral computed tomography with a photon counting detector and discrete reconstruction | |
JP2670635B2 (ja) | エネルギーサブトラクション画像の表示方法および装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 9 |
|
PLFP | Fee payment |
Year of fee payment: 10 |
|
PLFP | Fee payment |
Year of fee payment: 11 |
|
PLFP | Fee payment |
Year of fee payment: 13 |
|
PLFP | Fee payment |
Year of fee payment: 14 |
|
PLFP | Fee payment |
Year of fee payment: 15 |
|
PLFP | Fee payment |
Year of fee payment: 16 |
|
PLFP | Fee payment |
Year of fee payment: 17 |