FR2998160A1 - Procede de traitement d'images radiologiques en double energie - Google Patents
Procede de traitement d'images radiologiques en double energie Download PDFInfo
- Publication number
- FR2998160A1 FR2998160A1 FR1260967A FR1260967A FR2998160A1 FR 2998160 A1 FR2998160 A1 FR 2998160A1 FR 1260967 A FR1260967 A FR 1260967A FR 1260967 A FR1260967 A FR 1260967A FR 2998160 A1 FR2998160 A1 FR 2998160A1
- Authority
- FR
- France
- Prior art keywords
- images
- image
- region
- interest
- recombined
- 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.)
- Withdrawn
Links
Classifications
-
- 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/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
-
- 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/46—Arrangements for interfacing with the operator or the patient
- A61B6/461—Displaying means of special interest
- A61B6/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- 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/48—Diagnostic techniques
- A61B6/482—Diagnostic techniques involving multiple energy imaging
-
- 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/502—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of breast, i.e. mammography
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5217—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- 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/54—Control of apparatus or devices for radiation diagnosis
- A61B6/548—Remote control of the apparatus or devices
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Animal Behavior & Ethology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Dentistry (AREA)
- Human Computer Interaction (AREA)
- Physiology (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Primary Health Care (AREA)
- Epidemiology (AREA)
- Databases & Information Systems (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Procédé de traitement d'images radiologiques d'une région d'intérêt (56) d'un patient, le procédé comprenant les étapes suivantes : - obtention d'un ensemble d'images de la région d'intérêt (56) comprenant au moins deux images (I1, I2), acquis au moyen d'un système d'imagerie médicale, chaque image de l'ensemble d'images étant acquise selon une énergie de rayons X, - détermination (602) d'au moins une image recombinée (I3) par recombinaison d'une pluralité d'images formant un sous-ensemble d'images de l'ensemble d'images pour chaque recombinaison, - segmentation (603) d'au moins une image parmi l'ensemble d'images, pour détecter et isoler au moins une zone (101, 102) de la région d'intérêt (56), - fusion avec (604) l'au moins une image recombinée (I3) consistant à fusionner l'au moins une zone (101, 102) de la région d'intérêt (56) issue de la segmentation (603) afin de visualiser sur l'image recombinée l'au moins une zone.
Description
Procédé de traitement d'images radiologiques en double énergie DOMAINE DE L'INVENTION L'invention concerne l'imagerie médicale, et, en particulier, l'imagerie multi-énergie.
ETAT DE L'ART L'imagerie multi-énergie consiste à acquérir des images d'une même anatomie avec des rayons X d'énergies différentes.
En recombinant les images acquises selon différentes énergies, des images recombinées exploitant les propriétés d'atténuation des différents matériaux imagés peuvent être obtenues. Des images avec des énergies différentes de régions dans lesquelles a été injecté un produit de contraste peuvent fournir des informations fonctionnelles et permettent dans certains cas d'améliorer la détection de développements vasculaires anormaux et de lésions. L'imagerie multi-énergie est notamment utilisée en mammographie et permet dans certains cas d'améliorer la détection et caractérisation des cancers du sein. L'imagerie double énergie est un cas particulier d'imagerie multiénergie. En imagerie double énergie du sein par rayons X avec injection de produit de contraste (dual-energy contrast-enhanced spectral x-ray imaging of the breast en terminologie anglo-saxonne), des paires d'images en basse énergie et en haute énergie du sein sont acquises après injection d'un produit de contraste. 1 L'imagerie multi-énergie permet d'obtenir des informations morphologiques et fonctionnelles des tissus et des lésions dans la région du sein.
L'information morphologique peut être fournie par les images basse énergie - c'est-à-dire l'énergie utilisée en mammographie conventionnelle - tandis que l'information fonctionnelle est fournie par les images recombinées, ici les images double énergie, possiblement améliorées par un produit de contraste. Ces images recombinées sont produites en recombinant les images d'énergies différentes, par exemple basse énergie et haute énergie. Un avantage de l'imagerie multi-énergie est qu'elle permet de fournir simultanément des informations morphologiques et fonctionnelles pour permettre d'intégrer les informations des deux types et ainsi améliorer le diagnostic du cancer. Tel que connu, au cours d'une intervention, les images recombinées et en basse énergie sont visualisées côte à côte. Un problème est que, lorsqu'elles sont présentées séparément, le lien entre les informations issues des deux images doit être effectué mentalement par le radiologue. Ceci n'est pas sans problème car cela peut conduire à des interprétations erronées : la correspondance entre les deux images est en particulier difficile et peut donc conduire à des diagnostics erronés.
PRESENTATION GENERALE DE L'INVENTION L'invention permet de fusionner les informations issues d'images de différentes énergies afin d'en faciliter leur interprétation ultérieure. 2 A cet effet, l'invention propose un procédé de traitement d'images radiologiques d'une région d'intérêt d'un patient, le procédé comprenant les étapes suivantes : - obtention d'un ensemble d'images de la région d'intérêt comprenant au moins deux images de la région d'intérêt, l'ensemble d'images étant acquis au moyen d'un système d'imagerie médicale, chaque image de l'ensemble d'images étant acquise selon une énergie de rayons X, - détermination d'au moins une image recombinée par recombinaison d'une pluralité d'images formant un sous- ensemble d'images de l'ensemble d'images pour chaque recombinaison, - segmentation d'au moins une image parmi l'ensemble d'images, pour détecter et isoler au moins une zone de la région d'intérêt, - fusion avec l'au moins une image recombinée consistant à fusionner l'au moins une zone de la région d'intérêt issue de la segmentation afin de visualiser sur l'image recombinée ladite zone.
L'invention est avantageusement complétée par les caractéristiques suivantes, prises seules ou en une quelconque de leurs combinaisons : - les images du sous-ensemble d'image sont acquises selon des énergies différentes ; - l'image recombinée est une image double énergie ; - les images segmentées fusionnées avec chaque image recombinée sont issues du sous-ensemble d'image associé à cette image recombinée ; - l'étape de fusion comprend une étape de mise à l'échelle des contrastes de la zone dans l'image recombinée à partir des 3 différences de contraste au niveau d'un voisinage de la zone détectée ; - on obtient une première série d'images de projection et une deuxième série d'images de projection, les images de chaque série correspondant à des projections selon différentes positions angulaires, afin de construire une image tridimensionnelle fusionnée ; - après l'étape de segmentation, une étape de reconstruction d'une image tridimensionnelle de la zone détectée et isolée ; - avant l'étape de segmentation, une étape de reconstruction d'une image tridimensionnelle, l'étape de segmentation étant appliquée à l'image tridimensionnelle ; - l'ensemble d'images est issu d'images préalablement acquises et enregistrées. L'invention concerne également un système comprenant une unité de traitement comprenant des moyens pour la mise en oeuvre du procédé de traitement selon l'invention.
Et elle concerne aussi un produit programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé selon l'invention, lorsque celui-ci est mis en oeuvre sur un ordinateur Les avantages de l'invention sont multiples L'invention permet une présentation d'informations morphologiques et fonctionnelles réunies et intègre donc les forces des deux types d'informations. L'invention permet en particulier de présenter ces informations, leur localisation dans la région d'intérêt et les unes par rapport aux autres avec une grande précision de position. 4 L'invention pouvant être réalisée par un produit programme d'ordinateur et pouvant être mise en oeuvre à l'aide de systèmes d'imagerie existants, sa mise en place est aisée et peu coûteuse. En particulier pour les techniques d'imagerie par mammographie, l'invention permet de réunir les informations concernant les lésions mises en valeur par un produit de contraste et les microcalcifications dans une seule image. Ces éléments constituant des indicateurs de la présence ou de l'absence de cancer, l'invention permet une analyse plus aisée et plus précise d'images radiologiques. Ainsi, l'invention permet de faciliter et d'améliorer le diagnostic du cancer. L'invention peut s'appliquer à la mammographie mais également à d'autres techniques d'imagerie en double énergie.
PRESENTATION DES FIGURES D'autres caractéristiques et avantages de l'invention apparaitront lors de la description ci-après d'un mode de réalisation. Aux dessins annexés : - les figures 1 et 2 illustrent schématiquement deux exemples de systèmes d'imagerie médicale selon l'invention, - les figures 3 et 4 illustrent schématiquement des étapes d'un exemple de procédé selon l'invention. DESCRIPTION DE L'INVENTION Système d'imagerie médicale 5 - Système d'imagerie médical pour l'acquisition d'images radiologiques La figure 1 illustre schématiquement un système 100 d'imagerie médicale pour l'acquisition d'images radiologiques. Le système 100 d'imagerie médicale comprend un support 1 destiné à recevoir un patient 10 à examiner, une source 2 destinée à émettre un faisceau 3 de rayons X, un détecteur 4 disposé en face de la source 2 et configuré pour détecter les rayons X émis par la source 2, une unité de commande 6, une unité de stockage 7 et une unité d'affichage 8. La source 2 de rayons X et le détecteur 4 sont reliés par un bras 5 en forme de C. Un tel bras 5 est plus communément appelé arceau. Le bras 5 peut être orienté selon trois degrés de liberté. Le détecteur 4 peut être un capteur d'image à semi-conducteurs comprenant, par exemple, du phosphore d'iodure de césium (scintillateur) sur une matrice de transistor/photodiode en silicium amorphe. D'autres détecteurs adéquats sont : un capteur CCD, détecteur numérique direct qui convertit directement les rayons X en signaux numériques. Le détecteur 4 illustré sur la figure 1 est plan et définit une surface plane d'image. D'autres géométries peuvent bien entendu convenir. L'unité de commande 6 est connectée à l'arceau 5 par connexion filaire ou sans fil. L'unité de commande 6 permet de commander l'acquisition en fixant plusieurs paramètres tels que la dose de radiation à émettre par la source à rayons X et le positionnement angulaire du bras 5.
L'unité de commande 6 permet de commander la position du bras 5, c'est- à-dire la position de la source 2 par rapport au détecteur 4. L'unité de commande 6 peut comprendre un dispositif de lecture (non représenté) par exemple un lecteur de disquettes un lecteur de CD-ROM, DVD-ROM, ou des ports de connexion pour lire les instructions du procédé de traitement d'un support d'instructions (non montré), comme une disquette, un CD-ROM, DVD-ROM, ou clé USB ou de manière plus 6 générale par tout support de mémoire amovible ou encore via une connexion réseau. L'unité de stockage 7 est connectée à l'unité de commande 6 pour l'enregistrement des paramètres et des images acquises. Il est possible de prévoir que l'unité de stockage 7 est située à l'intérieur de l'unité de commande 6 ou à l'extérieur. L'unité de stockage 7 peut être formée par un disque dur ou SSD, ou tout autre moyen de stockage amovible et réinscriptible (clés USB, cartes mémoires etc.). L'unité de stockage 7 peut être une mémoire ROM/RAM de l'unité de commande 6, une clé USB, une carte mémoire, une mémoire d'un serveur central. L'unité d'affichage 8 est connectée à l'unité de commande 6 pour l'affichage des images acquises et/ou d'informations sur les paramètres de commande de l'acquisition.
L'unité d'affichage 8 peut être par exemple un écran d'ordinateur, un moniteur, un écran plat, un écran plasma ou tout autre type de dispositif d'affichage de type connu. Une telle unité d'affichage 8 permet à un praticien de contrôler l'acquisition des images radiologiques.
Le système 100 d'imagerie médicale est couplé à un système de traitement 200. Le système de traitement 200 comprend une unité de calcul 9 et unité de stockage 11. Le système de traitement 200 reçoit des images acquises et stockées dans l'unité de stockage 7 du système 100 d'imagerie médicale à partir desquelles il effectue un certain nombre de traitement (voir ci-après). La transmission des données de l'unité de stockage 7 du système 100 d'imagerie médicale vers l'unité de calcul 9 du système de traitement 30 200 peut être faite à travers un réseau informatique interne ou externe ou 7 à l'aide de tout support mémoire physique adéquat tel que disquettes, CD-ROM, DVD-ROM, disque dure externe, clé USB, carte SD, etc. L'unité de calcul 9 est par exemple un/des ordinateur(s), un/des processeur(s), un/des microcontrôleur(s), un/des micro-ordinateur(s), un/des automate(s) programmable(s), un/des circuit(s) intégré(s) spécifique(s) d'application, d'autres circuits programmables, ou d'autres dispositifs qui incluent un ordinateur tel qu'une station de travail. En variante, l'unité de calcul 9 peut comprendre un dispositif de lecture (non représenté) par exemple un lecteur de disquettes, un lecteur de CD-ROM ou DVD-ROM, ou des ports de connexion pour lire les instructions du procédé de traitement d'un support d'instructions (non montré), comme une disquette, un CD-ROM, un DVD-ROM ou une clé USB ou de manière plus générale par tout support de mémoire amovible ou encore via une connexion réseau.
En outre, le système de traitement comprend une unité de stockage 11 pour le stockage des données générées par l'unité de calcul 9. L'unité de calcul 9 peut être connectée à l'unité d'affichage 8 ou bien à une autre unité d'affichage (non représentée). - Système d'imagerie médicale pour l'acquisition d'images par mammographie La figure 2 illustre schématiquement un système 100 d'imagerie médicale pour l'acquisition d'images par mammographie. Ce système 100 d'imagerie médical diffère du système d'imagerie médicale décrit précédemment en ce que l'unité de commande 6 est connectée à un mammographe 50 au lieu de l'arceau de la figure 1. Le mammographe 50 comprend une source 51 de rayons X apte à émettre un faisceau 52 de rayons X vers un support 53 comprenant un bloc inférieur 54, sur laquelle repose un sein 56 d'une patiente, et une plaque supérieure 55, dite pelote de compression. La plaque supérieure 8 est mobile en translation verticale pour permettre de comprimer le sein 56 de la patiente contre le bloc inférieur 54. Le bloc inférieur 54 comprend en outre un détecteur 57, dont la surface de détection 58 est tournée vers le faisceau 52, directement sous le sein 56. Le faisceau 52 de rayons X émis par la source 51 rencontre le sein 56 de la patiente, le détecteur 57 captant alors les rayons X transmis par le sein 56 afin de réaliser une image mammographique.
Exemple de procédé de traitement d'images radiologiques En référence à la figure 3, il est décrit un procédé de traitement d'images radiologiques conforme à un mode de réalisation de l'invention.
Les étapes du procédé sont représentées schématiquement à la figure 4. Le procédé concerne le traitement d'images radiologiques d'une région d'intérêt 56 d'un patient 10. Un produit de contraste peut avoir été préalablement injecté dans la région d'intérêt 56.
Un tel procédé est mis en oeuvre dans une unité de calcul d'un système 100 d'imagerie médicale. Le système 100 d'imagerie médicale correspond par exemple à l'un des systèmes décrits ci-avant et représenté sur les figures 1 et 2. Dans l'exemple illustré par la figure 3, les images sont obtenues par un système 100 d'imagerie médical pour l'acquisition d'images par mammographie. Par image on entend une représentation pouvant être bidimensionnelle comme tridimensionnelle. Une image peut ainsi être une représentation tridimensionnelle d'un volume. Une telle image peut comprendre une série de représentations selon des modalités fixées. - Première étape : obtention d'images de la région d'intérêt 9 Le procédé comprend une première étape 601 d'obtention d'un ensemble d'images de la région d'intérêt 56 comprenant au moins deux images et 12 de la région d'intérêt 56.
L'ensemble d'images, et en particulier les deux images et 12 sont acquises au moyen du système 100 d'imagerie médicale. L'ensemble d'imgaes peut être issu d'images préalablement acquises et enregistrées. La première image et la deuxième image 12 peuvent être issues d'images préalablement acquises et enregistrées.
L'imagerie médicale multi-énergie consiste à acquérir des images d'une même anatomie avec des rayons X d'énergies différentes. Un tel protocole d'imagerie permet d'exploiter les propriétés d'atténuation des différents matériaux imagés : tissus humains, outils utilisés en radiologie interventionnelle, etc. La première étape 601 comme l'ensemble du procédé sont ici décrits en particulier par rapport à deux images. L'ensemble du procédé est aisément généralisable à l'imagerie multi-énergie impliquant plus de deux images. L'ensemble du procédé est également aisément généralisable à plus de deux images ou à des séries d'images. L'acquisition des images étant en multi-énergie, la première image est acquise selon une première énergie de rayons X, la deuxième image 12 étant acquise selon une deuxième énergie de rayons X différente de la première énergie de rayons X. Par exemple, la première image est une image basse énergie et la deuxième image 12 est une image haute énergie. Chaque image de l'ensemble d'images étant acquise selon une énergie de rayons X De manière préférée, on acquiert les deux images et 12 aux deux énergies différentes de telle manière qu'il n'y ait pas de mouvement significatif de la région d'intérêt entre ces deux images et 12. 10 En imagerie double énergie du sein par rayons X avec injection de produit de contraste, la basse énergie est généralement comprise entre 17 et 22 keV, et la haute énergie est généralement comprise entre 32 et 35 keV.
Compte tenu que les deux images sont issues d'acquisitions selon différentes énergie de rayons X, la première et la deuxième image permettent de visualiser des éléments différents. En imagerie multi-énergie on caractérise un matériau par la variabilité de son absorbance en fonction de l'énergie des rayons émis. Alvarez , Macovski, Lehmann et al. : « Generalized image combinations in dual KVP digital radiography » ; L. A. Lehmann, R. E. Alvarez, A. Macovski, W. R. Brody, N. J. PeIc, S. J. Riederer, and A. L. Hall, Med. Phys. 8, 659 (1981), 01:10.1118/1.595025) ont démontré que l'atténuation linéaire p d'un matériau peut être exprimée comme une combinaison linéaire de deux fonctions dépendant de l'énergie E : p = a, f,(E) + ap fp(E) Les deux constantes a, et ap caractérisent le matériau, et dépendent notamment du numéro atomique Z, et de la masse du matériau. On choisira donc selon l'effet recherché des énergies atténuant et/ou préservant des tissus humains différents. En imagerie double énergie, l'image haute énergie permet de visualiser certains éléments représentés par des zones de contraste particulier. L'image basse énergie permet de visualiser d'autres éléments représentés par des zones de contraste particulier. Les zones de contraste ne sont pas nécessairement les mêmes entre les deux images car elles ont été réalisées selon des énergies différentes mettant en valeur des matériaux différents. 11 La première image 11 de la région d'intérêt 56 représente différents éléments. Par exemple, en référence à la figure 3, sur la première image 11 sont visualisées différentes zones 101 102, 103, 104. Les zones 101, 102 peuvent être des zones de microcalcifications, les zones 103, 104 peuvent être des masses. La deuxième image 12 de la région d'intérêt 56 représente d'autres éléments. En référence à la figure 3, sur la deuxième image 12 sont visualisées 10 différentes zones 201, 203 et 204. Les zones 201, 203 et 204 peuvent être des masses. Alternativement au mode de réalisation présenté sur la figure 3, le système 100 d'imagerie médicale peut être un système d'imagerie 15 médicale pour l'acquisition d'images par tomosynthèse. Il peut s'agir d'un système 100 d'imagerie médicale tel que représenté sur la figure 1. Un tel système permet d'obtenir une image tridimensionnelle d'une région d'intérêt sous la forme d'une série de coupes successives. Le procédé de traitement d'images radiologiques est alors appliqué 20 à une pluralité de couples constitués chacun d'une première image de projection acquise selon une première énergie et d'une deuxième image de projection acquise selon une deuxième énergie. Chaque couple correspond à une projection 2D selon une orientation, par exemple une orientation angulaire, repérée par rapport à la perpendiculaire au support 25 1. L'orientation dite nulle est définie comme la plus proche de la perpendiculaire à la plaque 1. A partir des ensembles d'images de projection acquises selon chaque énergie, on peut reconstruire une image tridimensionnelle ou représentation tridimensionnelle, c'est-à-dire un volume reconstruit, pour 30 chacun des deux niveaux d'énergie. 12 - Deuxième étape : détermination d'une image recombinée Le procédé comprend une deuxième étape 602 de détermination d'au moins une image 13 multi-énergie recombinée par recombinaison d'une pluralité d'images formant un sous-ensemble d'images de l'ensemble d'images pour chaque recombinaison, par exemple par recombinaison de la première image acquise selon la première énergie et de la deuxième image 12 acquise selon la deuxième énergie. Cette recombinaison est par exemple de la forme : 13= a log +13 log 12 où a et R sont des constantes. Les constantes peuvent être déterminées en fonction des caractéristiques physiques d'un matériau constituant d'éléments que l'on veut faire disparaitre ou mettre en valeur. L'image 13 recombinée peut être une image spécifique de tout matériau, par exemple un produit de contraste, qui pourrait être présent dans le sein imagé. En imagerie du sein, cette recombinaison peut comprendre une fonction des images acquises, par exemple la première image et la deuxième image 12, dont la forme est déterminée afin que l'image 13 recombinée présente des informations relatives par exemple à l'épaisseur d'éléments du sein ou à la composition du sein. D'autres méthodes de recombinaison sont possibles. On peut par 25 exemple se référer au document Sylvie Puong : « Imagerie du sein multi spectrale avec produit de contraste » ; thèse de doctorat ; 2008, qui décrit plusieurs méthodes de recombinaison en imagerie multi-énergie. Cette étape 602 de détermination par recombinaison n'implique pas 30 une dégradation de l'image 13 recombinée. En revanche elle permet une suppression d'éléments. 13 Par exemple, sur la figure 3, elle permet de supprimer les masses 104 et les masses 204. Les éléments correspondant à des informations morphologiques, telles que les microcalcifications 101 et 102 sont supprimés.
Les éléments 201 et 303 de l'image 13 recombinée correspondent à des informations fonctionnelles. Il peut s'agir par exemple d'éléments résultant la diffusion d'un agent de contraste dans des lésions. Les éléments 303 représentent les zones où les masses 103 et les masses 203 sont superposées. Les masses 201 sont également 10 conservées. L'image multi-énergie peut être obtenue par recombinaison de plus de deux images. Ces images peuvent être acquises selon plus de deux énergies différentes 15 - Troisième étape : segmentation pour détecter et isoler une zone de la région d'intérêt Le procédé comprend une troisième étape 603 de segmentation d'au moins une image parmi l'ensemble d'images, pour détecter et isoler 20 au moins une zone de la région d'intérêt 56. Par exemple il peut s'agir d'une segmentation sur la première image 11 ou la deuxième image 12 ou les deux pour détecter et isoler au moins une zone de la région d'intérêt 56. La segmentation 603 peut comprendre une étape de sélection d'au 25 moins une zone de la région d'intérêt 56 présente sur la première image 11 ou la deuxième image 12. La zone peut représenter au moins une microcalcification. Par exemple, sur la figure 3, elle représente les microcalcifications 101 et 102. 30 Dans le cadre de l'analyse d'images mammographiques, i I est important de conserver la morphologie, les contrastes de l'échelle de gris 14 et l'arrangement spatial des microcalcifications 101 et 102. En effet, les microcalcifications 101 et 102 constituent d'importants discriminants de la présence ou de l'absence de cancer du sein.
Au cours de l'étape 603, par exemple, les microcalcifications 101 et 102 sont automatiquement détectées et segmentées dans la première image basse énergie à l'aide d'un outil de CAD (« Computer aided detection/diagnosis », diagnostic/détection assisté(e) par ordinateur en terminologie anglo-saxonne). L'outil de CAD repose sur des algorithmes tels que des méthodes de filtrage et de seuillage. D'autres méthodes relevant de la morphologie mathématique, des réseaux neuronaux, des modèles stochastiques ou des approches reposant sur les contours peuvent être utilisées dans la mise en oeuvre de l'outil de CAD. Un exemple d'utilisation d'outil de CAD est donné dans Sylvain Bernard, Serge Muller, Jon Onativia, « Computer-Aided Microcalcification Detection on Digital Breast Tomosynthesis Data: A Preliminary Evaluation », Digital Mammography, Studies in Fuzziness and Soft Computing Volume 210, pp 293-323.
A l'issue de l'étape 603 de segmentation, on obtient une image 14 segmentée représentant la zone détectée et isolée. Cette image 14 a préférentiellement les mêmes dimensions que les première et deuxième images et 12. Sur la figure 3, l'image 14 segmentée représente les microcalcifications 101 et 102. Un voisinage 401, respectivement 402, des microcalcifications 101, respectivement 102, est également représenté. L'image 14 segmentée comprend ainsi des informations relatives aux microcalcifications 101 et 102, à leur position spatiale, mais également des informations sur un voisinage des microcalcifications 101 et 102, notamment sur les 15 contrastes d'échelle de gris au voisinage des microcalcifications 101 et 102. Dans le cas d'images issues d'un système 100 d'imagerie médicale pour l'acquisition d'images par CESM (« contrast enhanced spectral mammography », mammographie spectrale à contraste amélioré en terminologie anglo-saxonne), comme peuvent l'être les images de la figure 3, la détection et la segmentation de l'étape 603 peuvent être réalisées sur la première image en basse énergie uniquement, la deuxième image 12 en haute énergie uniquement, ou les deux en combinaison. Dans le cas d'images issues d'un système 100 d'imagerie médicale pour l'acquisition d'images par CE-DBT (« contrast enhanced digital breast tomography », tomosynthèse mammaire numérique à contraste amélioré en terminologie anglo-saxonne), des images de projection selon différentes positions angulaires sont enregistrées, puis utilisées pour former une image reconstruite. On peut ainsi obtenir une première, respectivement deuxième, série d'images respectivement 12, de projection. A partir de plusieurs de ces images de projection respectivement 12, on peut obtenir un premier, respectivement deuxième, représentation tridimensionnelle ou volume reconstruit J1, respectivement J2. Chaque volume reconstruit J1 ou J2 peut comprendre des images de coupes successives de la région d'intérêt 56.
La détection et la segmentation peuvent alors, selon une première option, être réalisées sur la première série d'images de projection ou sur la deuxième série d'images 12 de projection uniquement, ou, selon une deuxième option, sur la première série d'images de projection en combinaison avec la deuxième série d'images 12 de projection.
Alternativement la détection et la segmentation peuvent, selon une troisième option, être réalisées sur le premier volume reconstruit J1 ou sur 16 2 99 8160 le deuxième volume reconstruit J2 uniquement, ou selon une quatrième option, sur le premier volume reconstruit J1 en combinaison avec le deuxième volume reconstruit J2. Une autre possibilité consiste, selon une cinquième option, à 5 réaliser la détection et la segmentation en utilisant : - la première série d'images de projection ou la deuxième série d'image 12 de projection d'une part, et - le premier volume reconstruit J1, par exemple associé à une énergie basse, ou le deuxième volume reconstruit J2, par 10 exemple associé à une énergie haute, d'autre part. Selon une sixième option : - la première série d'images de projection et la deuxième série d'image 12 de projection sont utilisées en combinaison avec - le premier volume reconstruit J1 ou le deuxième volume 15 reconstruit J2. Une septième option consiste à réaliser la détection et la segmentation en utilisant : - la première série d'images de projection et la deuxième série d'images 12 de projection d'un part, et 20 - le premier volume reconstruit J1 et le deuxième volume reconstruit J2 d'autre part. Une huitième option consiste à utiliser : - la première série d'images de projection ou la deuxième série d'images 12 de projection d'un part, et 25 - le premier volume reconstruit J1 et le deuxième volume reconstruit J2 d'autre part. - Quatrième étape : fusion avec l'image recombinée 30 Le procédé comprend une quatrième étape 604 de fusion de l'au moins une image 13 recombinée obtenue à l'issue de l'étape 602. La 17 quatrième étape 604 de fusion consiste à fusionner en une image 15 fusionnée la zone de la région d'intérêt 56 issue de la troisième étape 603 de segmentation afin de visualiser sur l'image 13 recombinée ladite au moins une zone. Ainsi des données fournies par la zone segmentée et des données fournies par l'image 13 recombinée sont réunies afin de produire un ensemble unique. En référence à la figure 3, on obtient une image 15 fusionnée résultant de la fusion de l'image 13 recombinée et des microcalcifications 101 et 102 issues de l'image 14 segmentée. L'étape 604 de correction comprend une étape de mise à l'échelle des contrastes de la zone de correction dans l'image 13 recombinée à partir des différences de signal au voisinage de la zone détectée.. Il peut s'agir de mise à l'échelle des contraste à partir des différences de signal au voisinage de la zone détectée, dans la première image et/ou la deuxième image 12. En référence à la figure 3, les voisinages 401 et 402 représentés sur l'image 14 segmentée contiennent les informations permettant la mise à l'échelle des contrastes de la zone de correction comprenant les microcalcifications 101 et 102. L'image 13 recombinée et l'image 14 segmentée partageant les mêmes modalités d'imagerie, tels que le format, la taille ou le positionnement des éléments, il est possible d'obtenir une image 15 fusionnée présentant une grande précision de position. Dans le cas illustré par la figure 3 d'images issues de systèmes 100 d'imagerie médicale pour l'acquisition d'images par CESM, les microcalcifications 101 et 102 segmentées de l'image 14 segmentée sont 18 ainsi fusionnées avec l'image 13 recombinée pour obtenir l'image 15 fusionnée. Dans le cas d'images issues de systèmes 100 d'imagerie médicale pour l'acquisition d'images par CE-DBT, la zone segmentée selon la première ou la deuxième option et l'image 13 recombinée de projection sont fusionnées. Un volume fusionné reconstruit J5 est alors obtenu à partir de plusieurs images fusionnées de projection en utilisant un algorithme de reconstruction tomographique.
Alternativement, des images 14 segmentées comprenant la zone segmentée selon la première ou la deuxième option sont utilisées pour obtenir un volume reconstruit de la zone segmentée. Il est alors possible de réaliser un traitement des éventuels artefacts de du volume reconstruit de la zone segmentée. Le volume reconstruit de la zone segmentée est ensuite fusionné avec un volume reconstruit puis recombiné ou un volume reconstruit à partir d'images de projection recombinées. Selon une autre possibilité, la zone segmentée selon la troisième option ou la quatrième option est fusionnée avec un volume reconstruit puis recombiné ou un volume reconstruit à partir d'images de projection 20 recombinées. Selon une possibilité supplémentaire, la zone segmentée selon la cinquième option, la sixième option, la septième option ou la huitième option est fusionnée avec une série d'images de projection recombinées. Un volume fusionné reconstruit est alors obtenu à partir des images ainsi 25 fusionnées de projection en utilisant un algorithme de reconstruction tomographique. Alternativement la zone segmentée selon la cinquième option, la sixième option, la septième option ou la huitième option est fusionnée avec un volume reconstruit puis recombiné ou un volume reconstruit à 30 partir d'images de projection recombinées.
19 Selon un exemple de mode de réalisation les images segmentées fusionnées avec chaque image recombinée sont issues du sous-ensemble d'image associé à cette image recombinée - Cinquième étape : affichage de l'image fusionnée Le procédé comprend une cinquième étape 605 d'affichage de l'image 15 fusionnée sur une unité d'affichage. Il peut s'agir de l'unité d'affichage 8 du système 100 d'imagerie médicale ou d'une autre unité 10 d'affichage. L'image 15 fusionnée présente réunies les informations morphologiques et fonctionnelles relatives à l'imagerie multi-énergie. Sur la figure 3, l'image 15 permet de différencier les 15 microcalcifications 101 qui sont incluses dans une masse 501 fusionnée correspondant à la masse 201 de l'image 13 recombinée, des microcalcifications 102 qui ne le sont pas. L'information morphologique correspondant aux microcalcifications 101 est ainsi fusionnée avec l'information fonctionnelle correspondant à la masse 501 fusionnée.
20 Produit programme d'ordinateur Le procédé de traitement ci-dessus décrit peut être exécuté par un 25 produit programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé lorsque celui-ci est mis en oeuvre sur un ordinateur. 20
Claims (3)
- REVENDICATIONS1. Procédé de traitement d'images radiologiques d'une région d'intérêt (56) d'un patient (10), le procédé comprenant les étapes suivantes : - obtention (601) d'un ensemble d'images de la région d'intérêt (56) comprenant au moins deux images (11, 12) de la région d'intérêt (56), l'ensemble d'images (11, 12) étant acquis au moyen d'un système (100) d'imagerie médicale, chaque image de l'ensemble d'images étant acquise selon une énergie de rayons X, - détermination (602) d'au moins une image recombinée (13) par recombinaison d'une pluralité d'images formant un sous-ensemble d'images de l'ensemble d'images pour chaque recombinaison, - segmentation (603) d'au moins une image parmi l'ensemble d'images, pour détecter et isoler au moins une zone (101, 102) de la région d'intérêt (56), - fusion avec (604) l'au moins une image recombinée (13) consistant à fusionner l'au moins une zone (101, 102) de la région d'intérêt (56) issue de la segmentation (603) afin de visualiser sur l'image recombinée ladite au moins une zone.
- 2. Procédé selon la revendication précédente, dans lequel les images (11, 12) du sous-ensemble d'image sont acquises selon des énergies 25 différentes.
- 3. Procédé selon l'une quelconque des revendications précédentes, dans lequel les images segmentées fusionnées avec chaque image recombinée sont issues du sous-ensemble d'image associé à cette image 30 recombinée. 21. Procédé selon l'une quelconque des revendications précédentes, dans lequel l'étape de fusion (604) comprend une étape de mise à l'échelle des contrastes de la zone (101, 102) dans l'image recombinée à partir des différences de contraste au niveau d'un voisinage (401, 402) de la zone (101, 102) détectée. 5. Procédé selon l'une quelconque des revendications précédentes, dans lequel on obtient une première série d'images (11) de projection et une deuxième série d'images (12) de projection, les images de chaque série correspondant à des projections selon différentes positions angulaires, afin de construire une image tridimensionnelle fusionnée (J5). 6. Procédé selon la revendication précédente, comprenant en outre après l'étape de segmentation (603), une étape de reconstruction d'une image tridimensionnelle de la zone (101, 102) détectée et isolée. 7. Procédé selon la revendication 5, comprenant en outre, avant l'étape de segmentation (603), une étape de reconstruction d'une image tridimensionnelle (Ji, J2), l'étape de segmentation (603) étant appliquée à l'image tridimensionnelle (Ji, J2). 8. Procédé selon l'une des revendications précédentes, dans lequel l'ensemble d'images est issu d'images préalablement acquises et enregistrées. 9. Système d'imagerie médicale comprenant une unité de traitement (200) comprenant des moyens pour la mise en oeuvre d'un procédé selon l'une des revendications précédentes. 10. Produit programme d'ordinateur comprenant des instructions de code de programme pour l'exécution des étapes du procédé selon l'une 22des revendications 1 à 8, lorsque celui-ci est mis en oeuvre sur un ordinateur. 23
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1260967A FR2998160A1 (fr) | 2012-11-19 | 2012-11-19 | Procede de traitement d'images radiologiques en double energie |
| US14/084,253 US20140140604A1 (en) | 2012-11-19 | 2013-11-19 | Method for processing dual-energy radiological images |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR1260967A FR2998160A1 (fr) | 2012-11-19 | 2012-11-19 | Procede de traitement d'images radiologiques en double energie |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| FR2998160A1 true FR2998160A1 (fr) | 2014-05-23 |
Family
ID=47754672
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR1260967A Withdrawn FR2998160A1 (fr) | 2012-11-19 | 2012-11-19 | Procede de traitement d'images radiologiques en double energie |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US20140140604A1 (fr) |
| FR (1) | FR2998160A1 (fr) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110738637A (zh) * | 2019-09-19 | 2020-01-31 | 华中科技大学 | 一种乳腺癌病理切片的自动分类方法及系统 |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2013111813A1 (fr) * | 2012-01-27 | 2013-08-01 | 株式会社 東芝 | Dispositif médical de traitement d'image |
| JP6105903B2 (ja) * | 2012-11-09 | 2017-03-29 | キヤノン株式会社 | 画像処理装置、画像処理方法、放射線撮影システム及びプログラム |
| WO2016186746A1 (fr) * | 2015-05-18 | 2016-11-24 | General Electric Company | Procédés et systèmes pour segmentation automatique |
| US10881358B2 (en) * | 2015-10-05 | 2021-01-05 | Koninklijke Philips N.V. | Tomosynthesis apparatus and method for characterizing a lesion in a breast |
| DE102016211766A1 (de) * | 2016-06-29 | 2018-01-18 | Siemens Healthcare Gmbh | Erzeugung einer Bildsequenz |
| US10096106B2 (en) * | 2016-11-10 | 2018-10-09 | General Electric Company | Combined medical imaging |
| CN112770674B (zh) | 2018-09-28 | 2024-08-06 | 豪洛捷公司 | 通过高密度元素抑制来合成乳腺组织的图像生成的系统和方法 |
| JP7430720B2 (ja) | 2018-11-25 | 2024-02-13 | ホロジック, インコーポレイテッド | マルチモダリティハンギングプロトコル |
| KR20210148132A (ko) | 2019-03-29 | 2021-12-07 | 홀로직, 인크. | 스닙-트리거링된 디지털 영상 보고서 생성 |
| WO2022064957A1 (fr) * | 2020-09-28 | 2022-03-31 | 富士フイルム株式会社 | Dispositif de traitement d'image, procédé de traitement d'image et programme de traitement d'image |
| US12186119B2 (en) | 2021-10-05 | 2025-01-07 | Hologic, Inc. | Interactive model interface for image selection in medical imaging systems |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040101104A1 (en) * | 2002-11-27 | 2004-05-27 | Avinash Gopal B. | Method and apparatus for soft-tissue volume visualization |
| US20090052612A1 (en) * | 2007-08-22 | 2009-02-26 | Xiaoye Wu | System and method of optimizing a monochromatic representation of basis material decomposed ct images |
| US20090103679A1 (en) * | 2007-10-19 | 2009-04-23 | General Electric Company | Image guided acquisition of quantitative dual energy data |
| US20120014504A1 (en) * | 2010-07-19 | 2012-01-19 | Samsung Electronics Co., Ltd. | Apparatus and method for imaging breast |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5389324B2 (ja) * | 2006-12-18 | 2014-01-15 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | X線断層撮影装置 |
-
2012
- 2012-11-19 FR FR1260967A patent/FR2998160A1/fr not_active Withdrawn
-
2013
- 2013-11-19 US US14/084,253 patent/US20140140604A1/en not_active Abandoned
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040101104A1 (en) * | 2002-11-27 | 2004-05-27 | Avinash Gopal B. | Method and apparatus for soft-tissue volume visualization |
| US20090052612A1 (en) * | 2007-08-22 | 2009-02-26 | Xiaoye Wu | System and method of optimizing a monochromatic representation of basis material decomposed ct images |
| US20090103679A1 (en) * | 2007-10-19 | 2009-04-23 | General Electric Company | Image guided acquisition of quantitative dual energy data |
| US20120014504A1 (en) * | 2010-07-19 | 2012-01-19 | Samsung Electronics Co., Ltd. | Apparatus and method for imaging breast |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110738637A (zh) * | 2019-09-19 | 2020-01-31 | 华中科技大学 | 一种乳腺癌病理切片的自动分类方法及系统 |
| CN110738637B (zh) * | 2019-09-19 | 2022-10-14 | 华中科技大学 | 一种乳腺癌病理切片的自动分类系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| US20140140604A1 (en) | 2014-05-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR2998160A1 (fr) | Procede de traitement d'images radiologiques en double energie | |
| Vedantham et al. | Digital breast tomosynthesis: state of the art | |
| US9842415B2 (en) | Method for processing tomosynthesis acquisitions in order to obtain a representation of the contents of an organ | |
| Sechopoulos | A review of breast tomosynthesis. Part II. Image reconstruction, processing and analysis, and advanced applications | |
| FR2963976A1 (fr) | Procede de traitement d'images pour la determination de zones suspectes dans une matrice tissulaire, et son utilisation pour la navigation 3d a travers la matrice tissulaire | |
| US9449403B2 (en) | Out of plane artifact reduction in digital breast tomosynthesis and CT | |
| US8903152B2 (en) | Methods and systems for enhanced tomographic imaging | |
| CN103040479B (zh) | 可能的灌注缺陷的确定 | |
| JP2015518392A (ja) | 光子計数撮像装置で投影データをビニングして画像を生成する方法とシステム{projectiondatabinningandimagegenerationinphotoncountingimagingmodality} | |
| US8611492B2 (en) | Imaging method for rotating a tissue region | |
| CN101082991A (zh) | 利用投影重建对象的方法以及执行该方法的装置 | |
| US20160262713A1 (en) | Creating a resultant image for a specifiable, virtual x-ray quanta energy distribution | |
| FR2919747A1 (fr) | Procede et systeme d'affichage d'images de tomosynthese | |
| US8811695B2 (en) | Methods, apparatus and articles of manufacture to adaptively reconstruct medical diagnostic images | |
| CN102117180A (zh) | 单屏幕多模态成像显示 | |
| FR2897461A1 (fr) | Dispositif de rayonnement x et procede de traitement d'images | |
| EP3821811B1 (fr) | Systèmes et procédés d'imagerie de diffusion cohérente à l'aide d'un détecteur de comptage de photons segmenté pour tomographie assistée par ordinateur | |
| US7564998B2 (en) | Image processing apparatus and method, and program | |
| FR2953964A1 (fr) | Procede de traitement d'images obtenues par tomographie ou tomosynthese a faible nombre de projections | |
| EP2443478B1 (fr) | Fenêtrage à énergie optimale d'événements diffusés en imagerie isotopique | |
| FR2960332A1 (fr) | Procede de traitement d'images radiologiques pour determiner une position 3d d'une aiguille. | |
| CN108065951A (zh) | 在ct成像方法中同时使用不同的造影剂 | |
| EP3597107B1 (fr) | Quantification de graisse à base de topogramme pour un examen de tomographie assistée par ordinateur | |
| Mundinger | 3D supine automated ultrasound (SAUS, ABUS, ABVS) for supplemental screening women with dense breasts | |
| US10507000B2 (en) | Method for simultaneous imaging of functional and morphological X-ray image data of a breast, diagnostic station, computer program product and computer-readable medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PLFP | Fee payment |
Year of fee payment: 4 |
|
| ST | Notification of lapse |
Effective date: 20170731 |