Verrou
Le modèle de Rayleigh, utilisé dans la plupart des approches de segmentation d'images échographiques, s'avère insuffisant pour la modélisation et la segmentation des régions échographiques correspondant au myocarde. Existe-t-il d’autres modèles plus adaptés permettant d’envisager la segmentation d’une plus grande classe d’images échographiques ? Par ailleurs, l’évaluation des algorithmes de segmentation, d’estimation de mouvement et de recalage d’images en imagerie thoracique et cardiovasculaire est un problème difficile du fait de l’absence de référence. Des modèles numériques de géométries et de comportements réalistes (thorax respirant, cœur battant, fluides en écoulement dans des structures vasculaires) constituent une telle référence. Leur développement permettrait d’acquérir simultanément un ensemble de connaissances qui pourraient alors être exploitées dans l’action ‘Analyse spatio-temporelle d’images’. Enfin, est-il possible de prévoir à l’avance l’impact d’une séance d’irradiation sur les tissus d’un patient cancéreux? La simulation numérique est une voix prometteuse en vue d’apporter des réponses.
Développements méthodologiques
Statistiques d’images US (D. Friboulet, O. Bernard Interaction avec l’équipe 3 Imagerie Ultrasonore)
Nous proposons le développement d'un modèle statistique basé sur les K-distributions, permettant la modélisation conjointe des régions sanguines et myocardiques et adapté à la segmentation. Le modèle des K-distributions a été développé initialement dans un objectif de caractérisation tissulaire pour des images échographiques enveloppe. L'application d'un tel modèle dans un cadre de segmentation pose un certain nombre de difficultés : d'une part la consistance des estimateurs des paramètres caractérisant la distribution n'est pas assurée dans la mesure où ces paramètres tendent vers l'infini pour les zones correspondant au sang; d'autre part la complexité de ce modèle rend difficile une implantation efficace pour une tâche de segmentation. Afin de pallier ces difficultés, nous proposons d'étendre le modèle des K-distributions à l'image échographique radiofréquence (RF). Les premiers résultats [BERN-06a, IEEE Trans Ultrason Ferroelectr Freq Control] [BERN-07, IEEE Trans Ultrason Ferroelectr Freq Control] montrent en effet que l'expression du modèle pour les données RF permet d'envisager une approximation efficace des K-distributions et le développement d'estimateurs plus stables.
Multimodality digital model of the thorax, heart and vessels (P. Clarysse, M. Orkisz, I. Magnin, D. Sarrut, Interaction with teams 2 Imagerie Volumique and 8 Imagerie morphologique et fonctionnelle)
Nous avons développé un modèle numérique du thorax respirant et du cœur battant [HADD-07, thèse], [HADD-05c, ITBM-RBM] (Figure 3). Ce modèle intègre une représentation de l’anatomie et du mouvement des différentes structures, les paramètres géométriques et fonctionnels associés et un ensemble d’images issues de la simulation d’acquisitions en différentes modalités, IRM (simulateur IRM SIMRI développé au laboratoire) et TEP (Simulateur TEP SORTEO, CERMEP, Lyon) notamment. Nous cherchons maintenant à inclure une variabilité de formes et de fonctions permettant ensuite de représenter une gamme d'anatomies personnalisées et de comportements normaux et pathologiques. D’autre part, de plus en plus de simulateurs d’images médicales devenant disponibles, une plateforme d’imagerie virtuelle multimodalités est en cours de développement [Sorina-08, MICCAI Workshop] dont un des principaux objectifs est de masquer la difficulté généralement rencontrée dans la mise en œuvre pratique de ces simulateurs. Il s’agit de reproduire dans un contexte de simulation le protocole clinique d’acquisition d’images dans diverses modalités.
The thorax-heart model also has to include a representation of the coronary tree obtained from image segmentation. The segmentation of 3D arborescent structures, such as the coronary tree, requires a reliable detection of bifurcations, as well as a capability of distinguishing the end of a branch from an apprent discontinuity owing to a local narrowing. Nous privilégions l’approche locale multi-échelle avec suivi de la branche mère et détection récursive des bifurcations et des branches filles ORKI-08, Mach Graph Vision [1]. La localité est assurée par un opérateur sphérique. Sa taille varie en permettant l’analyse multi-échelle et, au final, l’estimation du calibre local de la branche. Nous nous intéressons à l’évolution du contenu et de la surface de la sphère à travers les échelles. Dans certaines images, il est possible de trouver un seuil adaptatif permettant de séparer, à l’intérieur de la sphère, les structures d’intérêt du fond. On peut alors réaliser le suivi récursif quasi-automatiquement en analysant les composantes connexes sur la surface de la sphère CARR-05, CAIP [2], CARR-06a, CLEI [3], CARR-06b, CARS [4], CARR-06c, CCR [5], CARR-07, Int J CARS [6]. Une nouvelle technique de segmentation de la lumière artérielle par approche fast-marching a été développée et son auteur, Michael Baltaxe, s’est vu décerner le prix national Colombien Otto de Greiff BALT-07, EMBC [7] et FLOR-07a, CARS [8]. Dans le cas de l’extraction semi-interactive de l'arbre coronaire (un clic par branche), l’approche proposée est basée sur les moments d'inertie et un modèle élastique de l'axe. Avec cette méthode, nous avons été parmi les 8 meilleures équipes en finale du concours ‘Coronary Artery Tracking’ (CAT'08) associé à MICCAI 2008 HERN-08, MICCAI Workshop [9] (en ligne dans Midas Journal http://hdl.handle.net/10380/1401).
Extraction of coronary artery centerlines. Left: comparison between centerlines extracted starting from one click per artery (blue) and centerlines traced by experts (other colors), in four principale coronary branches (data from CAT'08 challenge). Right: details of the comparison in a right coronary artery.
En utilisant les modélisations a priori effectuées à partir des images TDM 4D (cf. section 2.1), nous avons pu proposer deux méthodes de reconstruction compensée en mouvement d'images Cone-Beam (une méthode analytique et une algébrique) [RIT-07, EMBC], [RIT-07, SPIE MI], [RIT-07, ICCR]. Ces méthodes permettent de reconstruire des images sans ou avec très peu d'artéfact de mouvement. D'autre part, un nouveau modèle d'estimation de mouvement directement à partir des projections Cone-Beam est proposé [VAND-08, MICCAI Workshop][VAND-07, ICCR]. Notons que, dans le cadre de ces travaux, un logiciel open-source (vv) a été développé et est proposé à la communauté [SERO-08, MICCAI Workshop].
Vascular digital simulation (G. Courbebaisse, M. Orkisz, Interaction with team 8 Imagerie morphologique et fonctionnelle)
Parmi les techniques utilisées en simulation numérique (éléments finis, méthodes particulaires) la méthode Lattice Boltzmann trouve des applications dans de nombreux domaines. Son adaptation pour la modélisation et la simulation organe ou modalité –dépendante est un challenge pour les prochaines années. Nous étudions sa mise en œuvre dans le cadre de la simulation d'écoulements viscoélastiques et plus spécifiquement de la thrombose dans les anévrismes. A priori, cette méthode qui prend en compte des grandeurs physiques à l’échelle microscopique pour simuler le comportement du système étudié à l’échelle macroscopique est adaptée aux systèmes multi-échelles auquel nous sommes confrontés. Ce travail de recherche s’appuie sur une collaboration entre le SPC (Scientific and Parallel Computing Group) de l'Université de Genève, le LIN (Laboratoire d’Ingénierie Numérique) de l’EPFL et CREATIS-LRMN [COUR-06b, MMM2006] et entre dans le cadre du projet européen COST P19. Les premiers résultats obtenus donnent une estimation du volume ‘thrombosé’ dans un anévrisme [MALA-08, Modern Physics], [NAVA-08, Computer Physics Communications]. A terme, l’objectif est de permettre un diagnostic patient spécifique pour détecter les anévrismes à risque qui nécessiteraient une intervention.
Par ailleurs, le modèle déformable simplexe de cylindre précédemment développé dans l’équipe a été exploité pour la simulation de la pose de stents artériels à partir d’une imagerie angiographique 3D FLOR-07b, ITBM [10].
Simulation Monte-Carlo pour la dosimétrie (D. Sarrut, Interaction avec l’équipe 2 Imagerie Volumique)
Le premier objectif consiste à proposer des représentations d'images spatio-temporelles adaptées aux logiciels de simulation d'interaction rayonnement-matière par méthode Monte-Carlo afin d'en accélérer les temps de calcul tout en préservant le degré de précision requis par les médecins [AYAD-06, AAPM]. Pour des milieux de géométrie complexe et dynamique, cela nécessite l'étude de grilles irrégulières isothétiques et d'algorithmes de navigation de particules au sein d'environnements 3D/4D. Sur le versant dosimétrie, le but est d'obtenir une carte volumique de la dose reçue par patient (ou dose transmise dans le cas de la simulation d'un imageur portal) lors d'une irradiation thérapeutique lorsque le milieu est décrit par une géométrie complexe et dynamique (plusieurs millions de voxels).
Ce travail est intégré dans l’ANR « blanche » SIMCA2 (« SIMulation CAncer CArbone ») dont l’objectif est de proposer une plate-forme permettant de simuler le dépôt de dose issu d'une irradiation par faisceaux d'ions Carbone (Hadronthérapie) à l'intérieur d'un patient virtuel. Nous avons jusqu’ici proposé un algorithme de navigation permettant d'accélérer fortement (jusqu'à un facteur 15) les simulations de distribution de dépôt de dose effectuée avec le code Monte-Carlo Geant4 [SARR-08, Medical Physics]. Cette méthode repose sur une segmentation et une carte de distance adaptées au déplacement des particules simulées. Une partie de ce travail s’intègre également dans le cadre du projet régional d'hadronthérapie ETOILE.
Résultats