FR2946172A1 - Procede de calibration d'une console de traitement d'images en imagerie tep/tdm - Google Patents
Procede de calibration d'une console de traitement d'images en imagerie tep/tdm Download PDFInfo
- Publication number
- FR2946172A1 FR2946172A1 FR1056677A FR1056677A FR2946172A1 FR 2946172 A1 FR2946172 A1 FR 2946172A1 FR 1056677 A FR1056677 A FR 1056677A FR 1056677 A FR1056677 A FR 1056677A FR 2946172 A1 FR2946172 A1 FR 2946172A1
- Authority
- FR
- France
- Prior art keywords
- contour
- tumor
- image
- light intensity
- pet
- 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
- 206010028980 Neoplasm Diseases 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000002685 pulmonary effect Effects 0.000 title claims abstract description 17
- 238000012545 processing Methods 0.000 title claims description 12
- 238000002600 positron emission tomography Methods 0.000 claims abstract description 52
- 210000004072 lung Anatomy 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims description 20
- 238000003384 imaging method Methods 0.000 claims description 11
- 206010056342 Pulmonary mass Diseases 0.000 claims description 4
- 238000003325 tomography Methods 0.000 claims description 4
- 208000020816 lung neoplasm Diseases 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 238000002591 computed tomography Methods 0.000 abstract description 28
- 238000013170 computed tomography imaging Methods 0.000 abstract description 8
- 239000000700 radioactive tracer Substances 0.000 description 14
- 210000001519 tissue Anatomy 0.000 description 11
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 4
- 238000012879 PET imaging Methods 0.000 description 4
- 239000008103 glucose Substances 0.000 description 4
- 230000005855 radiation Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 3
- 238000002059 diagnostic imaging Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 208000037841 lung tumor Diseases 0.000 description 3
- 210000000056 organ Anatomy 0.000 description 3
- 210000004879 pulmonary tissue Anatomy 0.000 description 3
- 230000002285 radioactive effect Effects 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 210000003484 anatomy Anatomy 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- ZCXUVYAZINUVJD-AHXZWLDOSA-N 2-deoxy-2-((18)F)fluoro-alpha-D-glucose Chemical compound OC[C@H]1O[C@H](O)[C@H]([18F])[C@@H](O)[C@@H]1O ZCXUVYAZINUVJD-AHXZWLDOSA-N 0.000 description 1
- YCKRFDGAMUMZLT-UHFFFAOYSA-N Fluorine atom Chemical compound [F] YCKRFDGAMUMZLT-UHFFFAOYSA-N 0.000 description 1
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000019522 cellular metabolic process Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 239000011737 fluorine Substances 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 230000003211 malignant effect Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 238000002661 proton therapy Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000000241 respiratory effect Effects 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10104—Positron emission tomography [PET]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30061—Lung
- G06T2207/30064—Lung nodule
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine (AREA)
Abstract
L'invention concerne un procédé de calibration d'une console (2) de traitement d'images d'une machine d'imagerie TEP/TDM couplée (1), comprenant les étapes de : a) mesure de la surface d'un nodule cancéreux pulmonaire sur une image TDM; b) récupération d'une image TEP dudit nodule pulmonaire effectuée par ladite machine d'imagerie; c) détermination d'une valeur seuil de luminosité définissant un contour sur l'image TEP, le contour étant tel que la zone de l'image à l'extérieur du contour présente une luminosité inférieure à cette valeur seuil et tel que la surface à l'intérieur du contour soit égale à la surface mesurée sur l'image TDM; d) mesure automatique de paramètres d'intensité lumineuse sur l'image TEP à l'intérieur dudit contour ; -répétition des étapes a) à d) pour différents nodules pulmonaires ; -génération d'une fonction calibrée définissant un seuil d'intensité lumineuse, par approximation mathématique des mesures obtenues aux étapes d).
Description
PROCEDE DE CALIBRATION D'UNE CONSOLE DE TRAITEMENT D'IMAGES EN IMAGERIE TEP/TDM
L'invention concerne les machines d'imagerie médicale, et en particulier 5 le traitement numérique en imagerie médicale destiné à identifier et quantifier des tumeurs pulmonaires.
La tomodensitométrie (TDM), dite aussi scanographie ou simplement scanner, est une technique d'imagerie médicale qui consiste à mesurer 10 l'absorption des rayons X par les tissus puis, par traitement informatique, à numériser et enfin reconstruire des images 2D ou 3D des structures anatomiques. Pour acquérir les données, on emploie la technique d'analyse tomographique ou "par coupes", en soumettant le patient au balayage d'un faisceau de rayons X. L'émetteur de rayons X effectue généralement une 15 rotation autour du patient en même temps que des récepteurs situés en face. Les récepteurs ont pour fonction de mesurer l'intensité des rayons après qu'ils ont été partiellement absorbés durant leur passage à travers le corps du patient. Les données obtenues sont ensuite traitées par ordinateur, ce qui permet de recomposer des vues en coupes bidimensionnelles puis des vues en trois 20 dimensions des organes. Grâce aux tomodensitomètres multidétecteurs (ou multi-barrettes) à acquisition spiralée (déplacement lent de la table d'examen durant l'acquisition), on peut obtenir une exploration très précise d'un large volume du corps humain pour un temps d'acquisition de quelques secondes. L'image obtenue en TDM est satisfaisante pour déterminer l'anatomie du 25 patient. De plus, une telle image présente une résolution et une netteté très satisfaisante. Cependant, une telle image s'avère imprécise pour définir les contours d'une tumeur. L'imagerie TDM est notamment inadaptée pour la détection de nouvelles extensions des tumeurs.
30 La Tomographie par Emission de Positron (TEP) est devenue un outil majeur pour le diagnostic, le bilan d'extension et la prise en charge d'un cancer. La TEP est une imagerie d'émission puisque la radioactivité provient du corps du patient. La TEP est une scintigraphie réalisée après injection d'un traceur faiblement radioactif (le plus souvent du 18F-Fluorodéoxyglucose (18FDG)) par 35 voie intraveineuse. Ce traceur est semblable au glucose et se fixe ainsi au niveau des tissus qui utilisent beaucoup de glucose. Les tissus cancéreux utilisent typiquement une grande quantité de glucose pour puiser l'énergie nécessaire à leur multiplication accélérée, et ont ainsi tendance à fixer une grande quantité de traceur. L'accumulation du traceur dans un tissu permet ainsi 40 de suspecter la présence d'une pathologie maligne. La partie des traceurs non fixée se dirige vers la vessie.
Le fluor présent dans le traceur émet un rayonnement de façon temporaire. Ce rayonnement induit l'émission de deux photons dans l'organisme selon des directions opposées. Une caméra TEP comprend une multitude de détecteurs destinés à détecter ce rayonnement provenant de l'organisme du patient. La détection de la trajectoire de deux photons en coïncidence (puisque arrivant en même temps sur deux détecteurs opposés) par le collimateur de la caméra TEP permet de localiser le lieu de leur émission. On peut donc déterminer la concentration du traceur en chaque point de l'organe. Cette information quantitative est restituée sous la forme d'une image faisant apparaître avec des contrastes de couleurs les différentes zones en fonction de leur concentration en traceur. L'examen se déroule en deux temps. À son arrivée, le patient est installé sur un lit et doit se reposer. Le traceur radioactif est ensuite injecté, dilué dans une solution salée. Le patient reste au repos une heure au moins. Ce temps permet au produit de bien se répartir dans l'organisme et d'être capté par les tissus (comme les tumeurs). Le patient est installé au calme dans une pièce sombre. À partir de l'injection, on évite toutes les activités qui peuvent perturber la diffusion du traceur, afin de ne pas stimuler les muscles, ce qui augmenterait la consommation de glucose et risquerait de modifier les résultats de l'examen.
Avant de s'installer sur le lit d'examen de la machine TEP, le patient se rend aux toilettes afin de vider au maximum le traceur non fixé de la vessie. Lors de l'examen, le patient est allongé sur un lit qui se déplace à l'intérieur d'un anneau détecteur. Seule une partie du corps se trouve à l'intérieur de l'appareil. L'enregistrement des images dure de 20 à 40 minutes selon les machines ou l'indication de l'examen. Ainsi la TEP permet de visualiser les activités du métabolisme des cellules : on parle d'imagerie fonctionnelle par opposition aux techniques d'imagerie dite structurelle comme celles basées sur les rayons X. La TEP permet ainsi de cerner avec précision de nouvelles tumeurs à un stade précoce et donc de suivre l'évolution d'une tumeur. La TEP ne permet par contre pas de localiser précisément la fixation du traceur sur le plan anatomique. La TEP fournit également une image relativement floue. En effet, l'image TEP restituant des différences d'intensité de pixels, cette image ne fournit pas de limite précise définissant les frontières de la tumeur.
Les images saisies en imageries sont restituées à un opérateur par l'intermédiaire d'un écran après avoir subi divers traitements numériques. Une console raccordée à la machine d'imagerie gère la restitution des images et propose différentes fonctions d'analyse pour interpréter les images restituées. Différentes techniques ont notamment été développées pour délimiter automatiquement des tumeurs sur des images TEP, notamment afin de déduire leur volume. Une méthode de contourage usuelle de la tumeur est basée sur des fonctions d'ajustement de seuils d'intensité de l'image. Pour différentes tailles de tumeurs, on calibre une fonction mathématique de fixation du seuil par des expériences réalisées sur des fantômes de calibration. Un fantôme de calibration comprend des sphères de taille calibrée, remplies de radio-isotopes en forte concentration et représentant des tumeurs. Le fantôme de calibration comprend par ailleurs un cylindre rempli d'une solution radioactive diluée représentant des tissus sains. Les sphères sont plongées dans le cylindre. Différentes images sont réalisées, avec différents paramètres d'images, par exemple différentes concentrations en solutions radioactives dans les sphères et dans le cylindre. On teste alors différentes valeurs de seuil de délimitation pour chaque image pour retenir celle qui génère un contour le plus proche des sphères calibrées. On réalise ensuite une interpolation entre les fonctions obtenues pour les différents paramètres d'image pour générer la fonction de seuil calibrée. Cette fonction peut ensuite être utilisée pour la délimitation de tumeurs de patients subissant un examen. Cependant, cette méthode de calibration ne donne pas des résultats suffisamment performants. En pratique, le fantôme utilisé comme modèle pour réaliser la calibration s'avère trop éloigné des conditions réelles d'examen, les tissus sains du patient autour de la tumeur présentant des hétérogénéités, contrairement à la solution radioactive diluée dans laquelle baignent les sphères. La délimitation automatique de la tumeur d'un patient en imagerie TEP s'avère donc peu précise.
Le couplage de l'imagerie TEP et de l'imagerie TDM (imagerie TEP/TDM couplée) permet d'améliorer la résolution spatiale, la sensibilité et la spécificité de l'imagerie fonctionnelle. L'imagerie TEP/TDM couplée est généralement mise en oeuvre sur une machine hybride. En effet, des publications récentes montrent que l'imagerie TEP/TDM réalisée sur une machine hybride est supérieure à l'imagerie seule TDM, ou TEP, ou à la combinaison visuelle des imageries TDM et TEP. Les machines à imagerie TEP/TDM couplée fournissent les images acquises à une console de traitement. Après identification d'une tumeur, la console de traitement croise généralement l'image TEP et l'image TDM pour définir le domaine de cette tumeur. Afin de garantir le traitement de toute la tumeur, on génère une image de la tumeur couvrant à la fois la tumeur identifiée sur l'image TEP et la tumeur identifiée sur l'image TDM. On maximise ainsi le domaine dans lequel la tumeur est susceptible d'être présente. Le domaine maximisé peut ensuite être traité par irradiation. La différence de domaine entre les deux images permet également de déterminer l'évolution de la tumeur.
Cette maximisation présente cependant des inconvénients. En effet, si une telle maximisation garantit l'irradiation de toute la tumeur, cette maximisation induit également l'irradiation de tissus sains inclus dans le domaine maximisé. L'irradiation de tissus sains génère un risque d'initiation d'un nouveau cancer. À cet égard, la maximisation du domaine de la tumeur n'est ni suffisamment objective ni suffisamment précise. L'invention vise à résoudre un ou plusieurs de ces inconvénients. L'invention porte ainsi sur un procédé de calibration d'une console de traitement d'images connectée à une machine d'imagerie tomographie par émission de positron/tomodensitométrie couplée dite TEP/TDM couplée, comprenant les étapes de : a) mesure par traitement numérique de la surface d'un nodule cancéreux pulmonaire sur une image TDM en coupe effectuée par ladite machine d'imagerie, le nodule pulmonaire étant un nodule pulmonaire isolé et non traité ; b) récupération d'une image TEP dudit nodule pulmonaire effectuée par ladite machine d'imagerie selon la même coupe ; c) détermination d'une valeur seuil de luminosité définissant un contour sur l'image TEP, le contour étant tel que la zone de l'image à l'extérieur du contour présente une luminosité inférieure à cette valeur seuil et tel que la surface à l'intérieur du contour soit égale à la surface mesurée sur l'image TDM; d) mesure automatique de paramètres d'intensité lumineuse sur l'image TEP à l'intérieur dudit contour ; - répétition des étapes a) à d) pour différents nodules pulmonaires ; - génération d'une fonction calibrée définissant un seuil d'intensité lumineuse pour contourer une tumeur en fonction desdits paramètres d'intensité lumineuse mesurés, la fonction calibrée étant une approximation mathématique des mesures obtenues aux étapes d). Selon une variante, ledit nodule de l'image de tomodensitométrie est un nodule du lobe supérieur d'un poumon.
Selon une autre variante, lesdits paramètres d'intensité mesurés à l'étape d) comprennent l'intensité lumineuse maximale et l'intensité lumineuse moyenne. Selon encore une variante, l'étape d) comprend également la mesure de l'intensité lumineuse minimale à proximité et à l'extérieur du contour.
Selon encore une autre variante, le procédé comprend la génération automatique d'un masque incluant le poumon comportant la tumeur lors de l'étape a), la mesure de l'intensité lumineuse minimale étant réalisée à l'intérieur de ce masque à l'étape d). Selon une variante, la fonction calibrée générée est fonction du rapport 40 entre Tmoy et Smin, Tmoy étant l'intensité lumineuse moyenne à l'intérieur du contour et Smin étant l'intensité lumineuse minimale à l'extérieur du contour.
Selon encore une variante, la fonction calibrée générée est de la forme : Tseuil= A* e[-B*Lnû)]*Tmax, avec A, B, C et D des paramètres obtenus lors de l'approximation mathématique des mesures des étapes d), Tseuil la valeur seuil de l'intensité lumineuse destinée à former le contour d'une tumeur, Tmax la valeur maximale de l'intensité lumineuse à l'intérieur du contour. Selon une autre variante, une fonction de seuil distincte est calibrée pour les plans sagittal, coronal et transversal. Selon encore une autre variante, la fonction calibrée est générée à partir d'images d'au moins 20 nodules pulmonaires distincts.
Selon encore une variante, la génération de la fonction calibrée comprend : - une première approximation mathématique des mesures obtenues aux étapes d), - l'élimination de plusieurs nodules pulmonaires dont les mesures sont les 15 plus éloignées de la première approximation mathématique, - une deuxième approximation mathématique des mesures restantes pour obtenir la fonction calibrée. D'autres caractéristiques et avantages de l'invention ressortiront clairement de la description qui en est faite ci-après, à titre indicatif et nullement 20 limitatif, en référence aux dessins annexés, dans lesquels : - la figure 1 est une représentation schématique d'une machine d'imagerie TEP/TDM couplée ; - la figure 2 est une image TDM d'un patient au niveau d'une coupe transverse ; 25 -la figure 3 est un agrandissement de l'image de la figure 2 centré sur une tumeur ; - la figure 4 est une image TEP du patient au niveau de la coupe transverse de la figure 2 ; - la figure 5 est un agrandissement d'une image retouchée de la figure 4 ; 30 -les figures 6 et 7 sont des diagrammes illustrant le résultat de différentes mesures effectuées sur les images, ainsi que des courbes approximant ces mesures.
L'inventeur a constaté qu'un nodule pulmonaire isolé, bien limité 35 visuellement et non traité présente un profil et une taille similaires sur une image de tomodensitométrie et sur une image TEP. L'image anatomique du nodule pulmonaire est alors extrêmement proche de son image fonctionnelle. Par conséquent, une calibration de la fonction de seuil de délimitation d'une image TEP basée sur l'utilisation d'une image de tomodensitométrie d'un tel nodule 40 peut être réalisée de façon particulièrement précise et réaliste pour améliorer la détermination du contour des différentes tumeurs de patients ultérieurs.
Une telle calibration s'avère encore plus précise lorsque le nodule est un nodule du lobe supérieur d'un poumon. En effet, on constate que cette zone des lobes supérieurs présente un moindre mouvement respiratoire des patients, ce qui rend les images de tomodensitométrie et de TEP particulièrement précises pour effectuer la calibration. On peut faire une calibration pour les lobes supérieurs des poumons, et une autre calibration pour le reste des tissus pulmonaires.
La première étape du procédé de calibration selon l'invention comprend la sélection d'images TEP et TDM effectuées avec la machine (ou une machine de la même série et pour les mêmes modes d'examens et de reconstruction d'images) pour différentes tumeurs pulmonaires de patients. Chaque image doit restituer un nodule pulmonaire isolé, bien limité visuellement, n'ayant pas fait l'objet d'un traitement. De cette manière, on arrive à conserver une grande identité entre les images TEP et TDM. On désigne généralement une tumeur plus ou moins arrondie présentant une dimension comprise entre 10 et 35 mm par le terme nodule. Un nodule isolé permet d'éviter d'introduire des interférences avec d'autres nodules lors des traitements d'image. Un nodule n'ayant pas fait l'objet d'un traitement permet de conserver une grande identité entre ses images TEP et TDM. Les cellules cancéreuses non traitées poursuivant leur consommation de traceur radioactif. Pour chaque tumeur pulmonaire sélectionnée, on récupère un couple d'images TEP et TDM selon une même coupe. Dans l'exemple illustré, les figures 2 et 3 sont des images TDM d'une coupe d'un nodule pulmonaire. Les figures 4 et 5 sont des images TEP selon la même coupe de ce module pulmonaire. Dans ces figures, le nodule 4 est entouré par un trait discontinu pour faciliter son identification. Durant la calibration, on peut procéder à la superposition de l'image TDM avec l'image TEP pour bien associer le nodule 4 sélectionné sur les deux images. On peut ainsi contourer les éléments anatomiques dans l'image TDM et les superposer sur l'image TEP pour faciliter la lecture de cette image. Pour des facilités de lecture, la forme de la tumeur 4 a été schématisée sur la figure 5 et les niveaux d'intensité lumineuse ont été représentés par des pavés comportant des motifs distincts. L'échelle d'intensité lumineuse croissante est indiquée en légende.
Lors d'une seconde étape, par une fonction de traitement d'image classique, on détermine un contour dans l'image TEP défini par un seuil d'intensité lumineuse Tseuil. A l'extérieur du contour, l'intensité lumineuse est inférieure au seuil et à l'intérieur du contour, la surface est égale à la surface mesurée sur l'image TDM de la tumeur. Le seuil d'intensité lumineuse définissant un tel contour peut par exemple être trouvé manuellement ou calculé automatiquement de façon itérative en testant successivement différentes valeurs de seuil d'intensité lumineuse.
Lors d'une troisième étape, sur l'image TEP, on mesure automatiquement d'une part l'intensité lumineuse maximale Tmax à l'intérieur du contour. Sur la figure 5, on a par exemple pu mesurer une valeur Tmax=20273. L'échelle d'intensité lumineuse utilisée est spécifique à la machine d'imagerie mentionnée ultérieurement. On mesure automatiquement par ailleurs l'intensité lumineuse moyenne Tmoy à l'intérieur du contour. Sur la figure 5, on a par exemple pu mesurer une valeur Tmoy=15209. Les mesures de l'intensité lumineuse maximale et de l'intensité lumineuse moyenne à l'intérieur d'un contour sont des fonctionnalités de traitement usuellement disponibles dans les applications de traitement d'image. L'automatisation de telles mesures d'intensité est aisément mise en oeuvre par l'homme du métier.
Avantageusement, on mesure également un paramètre caractéristique de chaque image TEP d'une coupe de nodule. On mesure par exemple l'intensité lumineuse minimale Smin sur l'image TEP à proximité du contour du nodule 4. Cette intensité lumineuse minimale est représentative du bruit induit dans l'image TEP par les tissus sains du poumon. Cette intensité lumineuse est ainsi représentative de la fixation minimale du traceur par les tissus sains pulmonaires. Sur la figure 5, on a par exemple pu mesurer une valeur Smin=400, correspondant à la zone entourée en pointillés. Pour déterminer Smin, on peut générer un masque sur l'image TDM, dans lequel on supprime par analyse numérique les zones hors des tissus pulmonaires. En appliquant le masque généré sur l'image TEP, la valeur Smin peut être calculée à l'intérieur de ce masque. Smin sera calculée pour chaque image TEP d'une coupe du nodule. Pour chaque coupe, et pour l'ensemble des nodules sélectionnés, on dispose donc des valeurs Smin, Tseuil, Tmax et Tmoy mesurés pour cette coupe. Tseuil peut être exprimé en pourcentage par rapport à l'intensité Tmoy. Pour chaque coupe, on dispose donc d'un couple (Tmoy, Tmax). Lors d'une quatrième étape, à partir de ces couples, on détermine une première équation Tmax=f(Tmoy). Par une analyse statistique, on détermine par exemple une droite représentative de cette fonction, comme illustré à la figure 6.
Pour améliorer la précision de cette équation, on peut éliminer des images considérées comme non-conformes car leur couple (Tmoy, Tmax) s'écarte nettement de la droite représentative. Dans l'exemple illustré à la figure 6, après avoir exclu des images non-conformes, on déduit l'équation suivante : Tmax=1,6381 *Tmoy-1691.9 (1) Dans cet exemple, on a pu déterminer une précision quadratique R2=0,9578, ce qui signifie que cette équation est particulièrement représentative des valeurs mesurées sur les images. Lors d'une cinquième étape, à partir des valeurs Tseuil, Tmax et Smin de chaque image, on détermine une deuxième équation Tseuil=g(Tmax/Smin). On détermine par exemple une courbe représentative de cette fonction, comme illustré à la figure 7. Comme pour la courbe de la figure 6, on peut éliminer des images non conformes car leurs valeurs s'écartent trop de la courbe déterminée. Dans l'exemple illustré à la figure 7, après avoir exclu des images non-conformes, on déduit l'équation suivante : Tseuil=142.13*(Tmax/Smin)-x245 (2) Dans cet exemple, on a pu déterminer une précision quadratique R2=0.8321 Lors d'une sixième étape, en combinant les équations (1) et (2), on déduit 15 l'équation (3) suivante : Tseuil= 142,13 * e[-0,2454*Ln(1,6381 *(Tmoy/Smin)-(1691,9/ Smin))] (3) Cette dernière équation (3) définit ainsi une fonction de seuil calibrée. On pourra définir une fonction de seuil différente pour le plan sagittal, coronal et transversal, ou une même fonction de seuil pour l'ensemble de ces 20 plans.
Les diagrammes des figures 6 et 7 ont été réalisés à partir d'images générées par une machine commercialisée par la société Philips sous la référence commerciale Gemini. Les diagrammes ont été réalisés à partir 25 d'images de tumeurs pulmonaires de 32 patients distincts. Seules les images de 26 patients ont été conservés pour réaliser la calibration de la fonction de seuil, certaines images ayant été exclues pour non-conformité comme détaillé auparavant.
30 Lors de l'utilisation de la console calibrée 2, une série d'images TEP et d'images TDM sont réalisées selon plusieurs coupes sur un patient. Au moins une image en coupe du patient est restituée sur l'écran de la console 2. L'utilisateur sélectionne une tumeur présente sur cette image en coupe. Pour chaque coupe de la tumeur, la console 2 définit un masque autour de la tumeur 35 dans l'image TDM. La console 2 place ce masque sur la tumeur dans l'image TEP correspondante. La console 2 mesure alors Smin dans ce masque. Puis Tmax est mesuré dans la tumeur. Ces mesures sont utilisées pour un calcul itératif du seuil d'intensité lumineuse Tseuil sur la base de l'équation (3). Pour la première itération, on introduit la valeur Tmax à la place de Tmoy 40 dans l'équation (3) pour déterminer une valeur Tseuill. On trace le contour pour cette valeur Tseuill sur l'image TEP et on mesure la valeur Tmoyl à l'intérieur de ce contour. On détermine une valeur Tseuil2 en introduisant les valeurs Tmoyl et Smin dans l'équation (3). On trace le contour pour cette valeur Tseuil2 sur l'image TEP et on mesure la valeur Tmoy2 à l'intérieur de ce contour. On répète ensuite le processus, jusqu'à ce que la valeur d'intensité lumineuse moyenne mesurée converge. On peut alors déterminer la valeur finale de Tseuil. On trace le contour autour de l'image TEP en coupe de la tumeur sur la base de cette valeur Tseuil. De façon connue en soi, on peut réaliser un lissage du contour de la tumeur au moyen d'un filtre numérique. La console 2 délimite ainsi le contour de la tumeur sur chacune des images TEP en coupe, coupes par exemple effectuées dans des plans sagittaux, coronaux et transversaux. De façon connue en soi, par interpolation entre les contours des différentes images en coupe repositionnées dans l'espace, la console 2 génère une représentation en trois dimensions de la tumeur.
La console 2 peut présenter différentes fonctionnalités à l'utilisateur, notamment une fonctionnalité de calcul du volume de la tumeur ou d'indication de ses dimensions suivant différents axes. Le volume de la tumeur permet en général aux praticiens d'affiner le dosage d'un traitement par radiothérapie, d'un traitement par protonthérapie ou d'un traitement par hadronthérapie.
L'image TEP peut donc être utilisée à la fois pour déterminer un volume précis de la tumeur et disposer de précieuses informations fonctionnelles sur l'évolution de cette tumeur. Par ailleurs, l'équation (3) étant calculée à partir de la valeur Smin, elle prend avantageusement en compte un facteur empirique de fixation lié à une partie bien spécifique qui caractérise (et différencie) les tissus pulmonaire des patients examinés. La fonction de seuil calibrée sur une machine et une console données peut bien entendu être utilisée sur d'autres machines et consoles de même type. Bien que le calibrage de la fonction de seuil ait été effectué avec des nodules pulmonaires, la console 2 peut bien entendu être utilisée pour délimiter le contour d'autres types de tumeurs, de dimensions différentes ou présentes sur d'autres organes.
Différents modèles d'équations (3) peuvent être déterminés pour différentes tailles de nodules. On pourra générer une fonction de seuil calibrée 35 par interpolation des différentes équations déterminées. Dans l'exemple décrit, l'équation (3) s'exprime en fonction de l'intensité lumineuse moyenne à l'intérieur d'un contour. Une telle équation s'avère particulièrement avantageuse pour parvenir à une valeur de seuil par un processus itératif, l'intensité lumineuse moyenne à l'intérieur du contour pouvant 40 varier à chaque itération. Cependant, l'équation (3) de la fonction seuil peut également être déterminée autrement, et prendre en compte d'autres paramètres que ceux décrits, sans pour autant sortir du cadre de l'invention. L'équation de la fonction de seuil peut par exemple tenir compte des propriétés du traceur radioactif.
Bien que l'on ait décrit de façon générale l'utilisation des intensités lumineuses dans la calibration de la console 2, dans la plupart des cas, l'intensité lumineuse s'exprimera plus simplement sous la forme d'un niveau de gris.
Claims (10)
- REVENDICATIONS1. Procédé de calibration d'une console (2) de traitement d'images connectée à une machine d'imagerie tomographie par émission de positron/tomodensitométrie couplée dite TEP/TDM couplée (1), comprenant les étapes de : a) mesure par traitement numérique de la surface d'un nodule cancéreux pulmonaire sur une image TDM en coupe effectuée par ladite machine d'imagerie, le nodule pulmonaire étant un nodule pulmonaire isolé et non traité ; b) récupération d'une image TEP dudit nodule pulmonaire effectuée par ladite machine d'imagerie selon la même coupe ; c) détermination d'une valeur seuil de luminosité définissant un contour sur l'image TEP, le contour étant tel que la zone de l'image à l'extérieur du contour présente une luminosité inférieure à cette valeur seuil et tel que la surface à l'intérieur du contour soit égale à la surface mesurée sur l'image TDM ; d) mesure automatique de paramètres d'intensité lumineuse sur l'image TEP à l'intérieur dudit contour ; -répétition des étapes a) à d) pour différents nodules pulmonaires ; -génération d'une fonction calibrée définissant un seuil d'intensité lumineuse pour contourer une tumeur en fonction desdits paramètres d'intensité lumineuse mesurés, la fonction calibrée étant une approximation mathématique des mesures obtenues aux étapes d).
- 2. Procédé de calibration selon la revendication 1, dans lequel ledit nodule de l'image de tomodensitométrie est un nodule du lobe supérieur d'un poumon.
- 3. Procédé de calibration selon la revendication 1 ou 2, dans lequel lesdits paramètres d'intensité mesurés à l'étape d) comprennent l'intensité lumineuse maximale et l'intensité lumineuse moyenne.
- 4. Procédé de calibration selon l'une quelconque des revendications précédentes, dans lequel l'étape d) comprend également la mesure de l'intensité lumineuse minimale à proximité et à l'extérieur du contour.
- 5. Procédé de calibration selon la revendication 4, comprenant la génération automatique d'un masque incluant le poumon comportant la tumeur lors de l'étape a), la mesure de l'intensité lumineuse minimale étant réalisée à l'intérieur de ce masque à l'étape d).
- 6. Procédé de calibration selon la revendication 4 ou 5, dans lequel la fonction calibrée générée est fonction du rapport entre Tmoy et Smin, Tmoy étant l'intensité lumineuse moyenne à l'intérieur du contour et Smin étant l'intensité lumineuse minimale à l'extérieur du contour.
- 7. Procédé de calibration selon la revendication 6, dans lequel la fonction calibrée générée est de la forme : Tseuil= A* e[-B*Ln(C*(Tmoy/Smin)û(D/ Smin))]*Tmax, avec A, B, C et D des paramètres obtenus lors de l'approximation mathématique des mesures des étapes d), Tseuil la valeur seuil de l'intensité lumineuse destinée à former le contour d'une tumeur, Tmax la valeur maximale de l'intensité lumineuse à l'intérieur du contour.
- 8. Procédé de calibration selon l'une quelconque des revendications 15 précédentes, dans lequel une fonction de seuil distincte est calibrée pour les plans sagittal, coronal et transversal.
- 9. Procédé de calibration selon l'une quelconque des revendications précédentes, dans lequel la fonction calibrée est générée à partir d'images 20 d'au moins 20 nodules pulmonaires distincts.
- 10. Procédé de calibration selon la revendication 9, dans lequel la génération de la fonction calibrée comprend : - une première approximation mathématique des mesures obtenues aux 25 étapes d), - l'élimination de plusieurs nodules pulmonaires dont les mesures sont les plus éloignées de la première approximation mathématique, - une deuxième approximation mathématique des mesures restantes pour obtenir la fonction calibrée. 30
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1056677A FR2946172B1 (fr) | 2010-08-19 | 2010-08-19 | Procede de calibration d'une console de traitement d'images en imagerie tep/tdm |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1056677A FR2946172B1 (fr) | 2010-08-19 | 2010-08-19 | Procede de calibration d'une console de traitement d'images en imagerie tep/tdm |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2946172A1 true FR2946172A1 (fr) | 2010-12-03 |
FR2946172B1 FR2946172B1 (fr) | 2013-02-22 |
Family
ID=43027451
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR1056677A Expired - Fee Related FR2946172B1 (fr) | 2010-08-19 | 2010-08-19 | Procede de calibration d'une console de traitement d'images en imagerie tep/tdm |
Country Status (1)
Country | Link |
---|---|
FR (1) | FR2946172B1 (fr) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102488528A (zh) * | 2011-12-07 | 2012-06-13 | 华中科技大学 | 一种层析成像几何参数的校准方法 |
-
2010
- 2010-08-19 FR FR1056677A patent/FR2946172B1/fr not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
BRAMBILLA M ET AL: "Threshold segmentation for PET target volume delineation in radiation treatment planning: The role of target-to-background ratio and target size", MEDICAL PHYSICS, AIP, MELVILLE, NY, US LNKD- DOI:10.1118/1.2870215, vol. 35, no. 4, 7 March 2008 (2008-03-07), pages 1207 - 1213, XP012115975, ISSN: 0094-2405 * |
FORD ERIC ET AL: "Tumor delineation using PET in head and neck cancers: Threshold contouring and lesion volumes", MEDICAL PHYSICS, AIP, MELVILLE, NY, US LNKD- DOI:10.1118/1.2361076, vol. 33, no. 11, 20 October 2006 (2006-10-20), pages 4280 - 4288, XP012091936, ISSN: 0094-2405 * |
MOUSSALLEM MAZEN PHILIPPE: "Optimisation du contourage des tumeurs pulmonaires à partir de l'imagerie PET-CT pour la planification dosimétrique des traitements en radiothérapie", MASTER 2 PHYSIQUE MÉDICALE (À LYON1-GRENOBLE1), pages 1 - 38, XP002563710, Retrieved from the Internet <URL:http://www.creatis.insa-lyon.fr/rio/Rapports_de_stages?action=AttachF ile&do=view&target=Moussalem2008.pdf> [retrieved on 20100105] * |
NESTLE URSULA ET AL: "Comparison of different methods for delineation of 18F-FDG PET-positive tissue for target volume definition in radiotherapy of patients with non-Small cell lung cancer", JOURNAL OF NUCLEAR MEDICINE, SOCIETY OF NUCLEAR MEDICINE, RESTON, VA, US, vol. 46, no. 8, 1 August 2005 (2005-08-01), pages 1342 - 1348, XP002563712, ISSN: 0161-5505 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102488528A (zh) * | 2011-12-07 | 2012-06-13 | 华中科技大学 | 一种层析成像几何参数的校准方法 |
CN102488528B (zh) * | 2011-12-07 | 2013-04-24 | 华中科技大学 | 一种层析成像几何参数的校准方法 |
Also Published As
Publication number | Publication date |
---|---|
FR2946172B1 (fr) | 2013-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP7203852B2 (ja) | 深層学習を使用した低線量petイメージングからの全線量pet画像の推定 | |
JP6999576B2 (ja) | 空間とスペクトル情報に基づく複数エネルギーのct画像におけるノイズ制御のためのシステムと方法 | |
JP2020168352A (ja) | 医用装置及びプログラム | |
JP7059298B2 (ja) | 心筋ct灌流画像合成 | |
EP1233700B1 (fr) | Procede d'utilisation d'un systeme d'osteodensitometrie, par rayonnement x bi-energie | |
US20060025680A1 (en) | Method and apparatus for contrast enhanced medical imaging | |
US8831323B2 (en) | Method and apparatus for measuring activity of a tracer | |
JP7359851B2 (ja) | 陽電子放出断層撮影(pet)のための人工知能(ai)ベースの標準取込み値(suv)補正及び変動評価 | |
FR2849241A1 (fr) | Procede et dispositif d'imagerie radiographique | |
Wein et al. | Self-calibration of geometric and radiometric parameters for cone-beam computed tomography | |
JP5820549B2 (ja) | 灌流イメージング | |
JP2019536538A (ja) | スペクトルctにおける骨及びハードプラークのセグメンテーション | |
US20100046820A1 (en) | Framing of positron emission tomography data to assess activity peak | |
FR2946172A1 (fr) | Procede de calibration d'une console de traitement d'images en imagerie tep/tdm | |
FR2927719A1 (fr) | Procede de traitement d'images obtenues par tomosynthese et dispositif associe | |
FR2868544A1 (fr) | Procede de quantification de la radioactivite de structures vivantes de petites dimensions mettant en oeuvre la tomographie par emission | |
Kruzhilov et al. | Whole-body PET image denoising for reduced acquisition time | |
Lee et al. | Human airway measurement from CT images | |
Debrun et al. | Volume measurements in nuclear medicine gated SPECT and 4D echocardiography: validation using a dynamic cardiac phantom | |
EP4154035B1 (fr) | Procede de traitement des donnees relatives a un examen radiologique d'un patient | |
KR101113747B1 (ko) | 마이크로 ct를 이용한 스텐트 내 재협착 평가 방법 | |
EP3566204B1 (fr) | Procede de determination d'un indice d'heterogeneite de la perfusion myocardique a partir d'images medicales | |
CN113057582B (zh) | 基于mlem的自适应fista初始图像的定量锥束x射线发光计算机断层扫描方法 | |
JP2023108605A (ja) | 方法、医用画像処理装置及びプログラム | |
Dixon et al. | A robust and versatile method for the quantification of myocardial infarct size |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 7 |
|
PLFP | Fee payment |
Year of fee payment: 8 |
|
PLFP | Fee payment |
Year of fee payment: 9 |
|
PLFP | Fee payment |
Year of fee payment: 10 |
|
PLFP | Fee payment |
Year of fee payment: 11 |
|
ST | Notification of lapse |
Effective date: 20220405 |