FR 2921177 A1 20090320 FR 0757639 A 20070917 PROCEDE DE TRAITEMENT D'IMAGES ANATOMIQUES EN VOLUME ET SYSTEME D'IMAGERIE METTANT EN îUVRE CE PROCEDE DOMAINE DE L'INVENTION ET ART ANTERIEUR La présente invention est relative à l'imagerie médicale et notamment le traitement d'images anatomiques en volume. Que l'on utilise la radiographie tomographique (scanner CT) ou d'autres techniques d'acquisition d'images en volume, l'acquisition des images peut concerner juste une partie de l'anatomie (coeur, poumon, tête) ou en contenir plusieurs (du cou jusqu'aux jambes, par exemple). Or, selon les parties d'anatomie ou les organes à examiner, les traitements d'images qui doivent être mis en oeuvre sont souvent différents, la performance de chaque traitement étant optimisée pour des images où seule la partie de l'anatomie à examiner est présente. Egalement, il est souhaitable pour le praticien de pouvoir disposer immédiatement d'images correspondant exactement à la partie de l'anatomie qu'il a à examiner, et ce même si l'acquisition s'est faite en couvrant une zone d'anatomie plus grande que celle sur laquelle il souhaite travailler. C'est pourquoi il est apparu souhaitable de pouvoir segmenter et identifier de façon automatique les images des différentes parties d'anatomie et des différents organes qui ont fait l'objet de l'acquisition. Il a déjà été proposé par le brevet US 7 177 453 une méthode d'analyse d'images 3D selon laquelle les images sont partitionnées en différents sous-volumes pour notamment permettre des traitements ultérieurs optimisés. La technique décrite dans ce document concerne toutefois principalement l'analyse de la partie haute du corps. PRESENTATION GENERALE DE L'INVENTION L'invention propose quant à elle un traitement permettant de partitionner les principales parties d'anatomie d'un individu : tête, cou, thorax, poumons, coeur, abdomen, membres inférieurs. Egalement, l'invention propose un traitement permettant d'isoler des vues de certains organes d'une anatomie : coeur, poumon, arche de l'aorte, etc. Plus particulièrement , il est proposé un procédé pour le traitement d'images d'anatomie en volume acquises par un système d'imagerie médicale selon lequel on partitionne les images en sous- volumes correspondant à différentes parties d'anatomie identifiées sur lesdites images, caractérisé en ce que, pour réaliser ce partitionnement, on traite les images acquises pour identifier sur celles-ci les poumons et en ce que, si les poumons sont identifiés, on détermine des coupes supérieure et inférieure délimitant ceux-ci selon l'axe de l'anatomie sur laquelle les images ont été acquises et on met en oeuvre un traitement distinct sur les images correspondant à des zones d'anatomie respectivement au dessus et en dessous de ces deux coupes. Ainsi, il est proposé d'utiliser les poumons en tant que repère physique pour partitionner les images. Notamment, même s'ils ne sont pas contenus dans l'examen, leur absence sur les images donne déjà une information intéressante sur ce qui n'est éventuellement pas présent (par exemple le coeur et le thorax) et donc sur ce qu'il n'est pas nécessaire d'essayer de détecter. Notamment, pour identifier les poumons sur les images, on met en oeuvre un traitement visant à soustraire l'air dans la trachée des images, selon lequel : • on détermine les différentes composantes connexes d'air qui apparaissent sur les images acquises, • on en détermine une excentricité moyenne 3D, • on compare à un seuil les excentricités moyennes 3D ainsi calculées pour ne conserver que celles de ces composantes connexes dont l'excentricité 3D est inférieure à un premier seuil et dont la longueur selon l'axe z est supérieure à un deuxième seuil, les composantes connexes ainsi isolées étant ensuite soustraites des images. Egalement, lorsque des poumons sont détectés, on met en oeuvre un traitement en vue d'identifier et de partitionner les images correspondant au coeur. Toutefois, on compare préalablement la longueur de l'anatomie figurant sur les images à un seuil donné, ledit traitement en vue d'identifier et de partitionner les images correspondant au coeur n'étant mis en oeuvre que lorsque ladite longueur est supérieure audit seuil. Notamment, pour déterminer la coupe correspondant à la limite supérieure du coeur, on détermine au moins trois coupes axiales sur laquelle on identifie l'aorte ascendante et l'aorte descendante. Egalement, lorsque les poumons sont détectés, on détermine un profil axial de surface occupée par les voxels sur les coupes en dessous de la coupe limite inférieure des poumons et on détermine la présence de l'abdomen et de différentes parties des jambes en analysant ce profil. Lorsque les poumons ne sont pas détectés, on compare la longueur de l'anatomie ayant fait l'objet de l'acquisition des images à un seuil et si cette longueur est supérieure audit seuil, on détermine un profil axial de surface occupée par les voxels sur les coupes et on détermine la présence de l'abdomen et de différentes parties des jambes en analysant ce profil. Lorsque l'analyse du profil ne permet pas de déterminer la présence de l'abdomen et de différentes parties des jambes, on détermine un profil axial de l'air apparaissant sur les coupes et on analyse ce profil pour en déduire une information sur la zone de la tête ou du cou présente sur les images. Par ailleurs, lorsque les poumons sont détectés, on détermine un profil axial de l'air apparaissant sur les coupes au-dessus de la coupe limite supérieure des poumons et on analyse ce profil pour en déduire une information sur la zone de la tête ou du cou présente sur les images. Notamment, on compare au préalable à un seuil donné la longueur d'anatomie au-dessus de la coupe supérieure limite des poumons, la détermination du profil axial de l'air n'étant mise en oeuvre que lorsque ladite longueur est supérieure audit seuil. Lorsque les poumons ne sont pas détectés, on compare la longueur de l'anatomie ayant fait l'objet de l'acquisition des images à un seuil et si cette longueur est inférieure audit seuil, on détermine un profil axial de l'air apparaissant sur les coupes et on analyse ce profil pour en déduire une information sur la zone de la tête ou du cou présente sur les images. En outre et indépendamment du procédé présenté ci-dessous, l'invention concerne également un procédé pour le traitement d'images en volume acquises par un système d'imagerie médicale, caractérisé en ce que : • on parcourt les coupes axiales perpendiculaires à l'axe de l'anatomie pour identifier une coupe axiale correspondant à une limite supérieure du coeur, • on détermine sur cette coupe deux composantes connexes correspondant à l'aorte ascendante et à l'aorte descendante, • on détermine les composantes connexes en volume qui correspondent à ces composantes connexes, • on fusionne lesdites composantes connexes en volume, • on détermine sur la composante connexe en volume ainsi obtenue les centres des composantes connexes qui lui correspondent sur une pluralité de coupes axiales, • on détermine un plan optimisant une distance par rapport au nuage de points constitués par ces différents centres. L'invention propose également un système d'imagerie médicale comportant des moyens d'acquisition d'images d'anatomie en volume et des moyens de traitement de ces images aptes à partitionner les images en sous-volumes correspondant à différentes parties d'anatomie identifiées sur lesdites images, caractérisé en ce que lesdits moyens comportent des moyens aptes à mettre en oeuvre un procédé de traitement du type évoqué ci-dessous. Elle propose également un programme mémorisé sur des moyens aptes à être lus par un ordinateur, caractérisé en ce qu'il comporte des données et/ou instructions pour mettre en oeuvre un tel procédé. PRESENTATION DES FIGURES D'autres caractéristiques et avantages de l'invention ressortiront encore de la description qui suit, laquelle est purement illustrative et non limitative, et doit être lue en regard des figures annexées sur lesquelles : • les figures 1 et 2 sont deux représentations schématiques d'un système d'imagerie CT ; • la figure 3A illustre une partie d'un arbre décisionnel pouvant être utilisé dans un mode de réalisation ou de mise en oeuvre possible de l'invention ; • la figure 3B illustre une autre partie de cet arbre décisionnel ; • la figure 4 illustre différentes étapes d'un traitement possible pour la détection des poumons ; • les figures 5A à 5E sont des images en coupe illustrant différentes étapes du traitement de détection de poumons de la figure 4 ; • la figure 6 illustre les étapes d'un exemple de traitement permettant la détection d'une fenêtre autour du coeur ; • les figures 7A à 7G sont des images en coupe illustrant différentes étapes du traitement de détection du coeur de la figure 6 ; • la figure 8 illustre quant à elle un traitement permettant la détection de l'arche aortique ; • les figures 9A à 9F sont des images en coupe illustrant différents résultats issus de la détection de l'arche aortique ; • la figure 10 illustre quant à elle les étapes du traitement permettant la détection de la partie basse du corps et notamment l'abdomen et les jambes ; • la figure 11 illustre la correspondance entre un profil obtenu en mettant en oeuvre le traitement de la figure 10 et différents os des membres inférieurs et de l'abdomen ; • la figure 12 illustre quant à elle les différentes étapes d'un traitement permettant la détection de la partie haute du corps, et notamment de la tête et du cou ; • la figure 13 illustre la correspondance entre un profil obtenu en mettant en oeuvre le traitement de la figure 12 et différents os de la tête et du cou. DESCRIPTION D'UN OU PLUSIEURS MODES DE REALISATION ET DE MISES EN OEUVRE RAPPELS SUR LES SYSTEMES D'IMAGERIE 3D Un exemple de système d'imagerie 10 permettant de mettre en oeuvre l'invention est illustré schématiquement sur les figures 1 et 2. Dans cet exemple, cet appareil d'imagerie est un appareil de tomographie CT, d'autres techniques d'imagerie en volume (imagerie 3D) étant toutefois envisageables : imagerie par résonance magnétique ou imagerie par émission de positron par exemple. De manière classique, cet appareil de tomographie 10 comprend : • une source 12 de rayons X, • un collimateur 14 qui définit les dimensions et la forme du faisceau 16 de rayons X qui traverse une région dans laquelle un sujet tel qu'un patient humain 18 est installé, • un capteur 22, qui reçoit la partie 20 du rayonnement qui traverse l'anatomie 18 sur laquelle on souhaite réaliser l'examen. Le capteur 22 peut être une matrice de détecteurs, chacun de ces détecteurs générant en sortie un signal analogique représentant l'intensité des rayons X incidents. Ces signaux sont acquis et traités pour reconstruire, en temps réel ou presque, une image des détails du sujet. De manière classique, les signaux sont enregistrés selon plusieurs angles autour du patient de manière à collecter plusieurs vues radiographiques. La source 12 est contrôlée par une unité de contrôle 24 qui comprend un dispositif 30 de commande des rayons X et qui commande également, par l'intermédiaire d'un dispositif de contrôle de motorisation 32, d'une part le mouvement de rotation de la source 12 et du capteur 22 sur le portique 52 qui les porte (commande du système de rotation 26) et d'autre part le déplacement linéaire du lit 58 qui porte le patient par rapport au portique 52 (commande du système de déplacement linéaire 28). Par ailleurs, l'unité de contrôle 24 comprend un système d'acquisition de données 34 qui reçoit les signaux analogiques en sortie du capteur 22 et qui les convertit en signaux digitaux pour le traitement ultérieur par un processeur, par exemple un ordinateur 36. L'ordinateur 36 peut comprendre û ou communiquer avec û une mémoire 38 qui peut stocker les données traitées par l'ordinateur 36, ou les données à traiter par l'ordinateur 36. N'importe quel type de dispositif de mémoire accessible par un ordinateur et permettant de stocker la quantité voulue de données et/ou de codes, peut être utilisé. De plus, la mémoire 38 peut comporter un ou plusieurs dispositifs de mémoire, de types semblables ou différents, pouvant être locaux ou distants par rapport au système 10. Les dispositifs de mémoire peuvent stocker des données, des paramètres de traitement et/ou des programmes informatiques pour exécuter les différents processus décrits ici. L'ordinateur 36 est typiquement utilisé pour contrôler l'appareil d'imagerie 10. A cet effet, l'ordinateur 36 est configuré pour recevoir des commandes et des paramètres d'acquisition par un opérateur par l'intermédiaire d'un poste de travail 40, classiquement équipé d'une souris, d'un clavier et/ou d'autres périphériques. Un écran d'affichage 42 couplé au poste de travail permet l'affichage de l'image reconstruite. L'image peut également être imprimée au moyen d'une imprimante 43, qui peut être couplée au poste de travail. En outre, le poste de travail peut également être couplé à un système 44 d'archivage et de communication d'images (PACS ou � Picture Archiving and Communication System � selon la terminologie anglo-saxonne). Le PACS 44 peut être couplé à un système distant 46, afin que d'autres personnes puissent accéder à l'image et aux données d'image à distance. Bien sûr, cet exemple n'est en aucun cas limitatif et chacun des dispositifs présentés ci-dessus peut être couplé à d'autres dispositifs, selon l'implémentation qui est souhaitée. Les données collectées par le capteur 22 subissent typiquement une correction, un prétraitement et/ou un étalonnage dans le système d'acquisition 34 et/ou l'ordinateur 36 pour conditionner les données afin de représenter les intégrales des coefficients d'atténuation des objets analysés le long des rayons allant de la source au détecteur. Les données traitées, communément appelées images de projection, peuvent être utilisées par un algorithme de reconstruction. En tomographie, on acquiert un certain nombre d'images de projection, chacune selon un angle différent par rapport au sujet et/ou au détecteur. Les algorithmes de reconstruction tomographique sont bien connus de l'homme du métier et permettent de former une image tridimensionnelle du volume balayé (images en volume). ARBRE DE DECISION On va maintenant décrire un prétraitement qui peut être mis en oeuvre sur les images acquises par un système d'imagerie 3D. Ce prétraitement partitionne les images acquises grâce au système 10 ù ou tout autre système d'imagerie 3D û en sous-volumes correspondant aux différentes parties d'anatomie du corps du patient ou de la partie du corps du patient qui a fait l'objet de l'acquisition : tête, cou, thorax, abdomen, membres inférieurs. Il met en outre en évidence différents organes ou parties d'organes qu'il peut être souhaitable d'isoler pour les examens : poumons, coeur, arche aortique, 15 etc. A cet effet, ce prétraitement se déroule selon l'arbre de décision 100 qui est illustré sur les figures 3A et 3B. Le principe de cet arbre de décision consiste à commencer à partitionner les images selon qu'elles comprennent ou non les poumons puis à traiter les 20 parties au-dessus et en dessous des poumons. Ainsi, dans une première étape 110, le traitement met en oeuvre sur les images un test de détection des poumons. Si les poumons sont effectivement détectés sur certaines des images, les images acquises sont considérées comme se trouvant au moins au niveau du 25 thorax et l'arbre de décision se poursuit avec les différentes étapes illustrées sur la figure 3A. Le traitement met alors en oeuvre un test 120 sur la longueur de poumons visible sur les images et compare cette longueur à une longueur seuil donnée (7cm par exemple). 30 Lorsque la longueur des poumons est supérieure à cette longueur seuil, le traitement met en oeuvre un algorithme de détection du coeur (test 130). En parallèle au test 120, il est également mis en oeuvre sur les images des tests 150a et 150b permettant de partitionner les images qui se trouvent au-dessus ou en dessous des poumons. Pour les images au-dessus des poumons, le traitement met en oeuvre un test sur la longueur de l'anatomie visible sur les images et au-dessus des poumons, qu'il compare à un seuil donné, par exemple de 10cm. Si cette longueur est supérieure audit seuil, on met alors en oeuvre deux tests 170 et 180 de détection respectivement de la partie tête et de la partie cou. Dans le cas contraire, c'est-à-dire si la longueur de l'anatomie qui est au dessus des poumons et qui est visible sur les images est inférieure audit seuil, le traitement considère que l'on se trouve uniquement en présence de la partie cou. Pour les images en dessous des poumons, le traitement met en oeuvre des tests de détection des jambes et de l'abdomen. Ainsi, si les jambes sont détectées (test 190), le traitement met ensuite en oeuvre un test pour la détection de l'abdomen (200). Si les jambes ne sont pas détectées, le traitement considère que les images se trouvent immédiatement au niveau de l'abdomen. La figure 3B illustre l'autre partie de l'arbre de décision, dans le cas où les poumons ne sont pas détectés lorsque le traitement met en oeuvre le test 110. Dans cette hypothèse, le traitement met alors en oeuvre un test 210 comparant la longueur de la partie d'anatomie qui a fait l'objet de l'acquisition d'images à une longueur seuil, qui est par exemple de 15cm. Si celle-ci est supérieure audit seuil, le traitement met en oeuvre les tests 170 et 180 déjà évoqués en référence à la figure 3A de détection de la partie tête et de la partie cou. Le reste des images est considéré comme faisant partie de l'abdomen. Si par contre la longueur de la partie de l'anatomie qui se trouve sur les images est inférieure audit seuil, la détection met en oeuvre une détection des jambes au moyen du test 200, également déjà évoqué. Ainsi, l'arbre décisionnel tient compte de facteurs anatomiques et d'informations sur l'examen (par exemple la longueur). DETECTION DES POUMONS L'étape 110 de détection des poumons met en oeuvre un traitement permettant la détection de l'air sur les images, lequel est généralement bien défini sur les histogrammes puisqu'il ne présente pas beaucoup d'atténuation (valeur 5 d'intensité la plus petite des histogrammes). Ce traitement vise notamment à identifier les autres organes susceptibles de comprendre de l'air en leur intérieur pour les éliminer, voire également à identifier l'extérieur du corps (l'air entourant le patient faisant très souvent partie des acquisitions). 10 Les étapes de ce traitement sont les suivantes : - Seuillage du volume (étape 111) On applique un seuillage au volume original entre -1023 H.U. et -500 H.U. pour préserver l'air. - Filtrage de taille 3D (étape 112) 15 On enlève les composantes connexes 3D dont la taille est petite (afin d'éliminer des structures telles que des parties du colon). - Détection et séparation de la trachée (étape 113) La trachée relie les poumons avec l'air à l'intérieur de la tête. Elle naît dans le cou et fait suite au larynx. Elle se prolonge ensuite dans le thorax. Elle est en 20 forme de fer à cheval, dont la paroi postérieure est musculaire. La présence de la trachée et de l'air de la tête est assez gênante lorsque l'on souhaite détecter des poumons. Pour permettre de définir le thorax à partir des poumons et ne garder pour le thorax que la partie du corps où les poumons sont présents, on supprime des images l'air correspondant à la trachée. 25 On notera qu'une façon de détecter la trachée est d'utiliser sa forme en fer à cheval, qui, au niveau des poumons, s'approche d'une forme plus circulaire. Ainsi, une mesure du degré de cette circularité est prise en compte grâce à l'excentricité, laquelle est définie à partir des moments centraux : = (x- x)' C@- y)3 y io 30 Les moyennes étant calculées comme suit : , Cly = ~y y N où N représente la taille de la composante, l'excentricité pouvant quant à elle s'exprimer de la façon suivante : '\/(P20 - /102) 2 - 4/1i 1 /120 + /102 Les paramètres de cette formule représentent les moments centraux définis précédemment. Cette expression fournit une mesure pour chaque composante connexe (CC) apparaissant sur les vues en coupe ; plus sa valeur s'approchera de 0, plus on peut considérer que la composante connexe correspondante est circulaire, l'excentricité d'un cercle valant 0. Une composante connexe est une portion d'une image 2D ou 3D dont les éléments (pixels ou voxels) sont tous connectes entre eux. Chaque image segmentée peut se décomposer en une liste de composantes connexes indépendantes. Traitement de détection de la trachée Pour détecter la trachée, on utilise un traitement sur les images en coupe axiale (images 2D), c'est-à-dire sur les images selon des coupes perpendiculaires à l'axe principal de l'anatomie (axe parallèle à la longueur de la table et appelé axe z dans la suite du texte). Les étapes sont les suivantes : - Evaluation de l'excentricité et la taille en 2D On parcourt les coupes axiales et on réalise sur celles-ci un étiquetage des composantes connexes (étape 113a). Puis on calcule l'excentricité de chaque CC identifié par cet étiquetage (étape 113b), pour ne garder que les CCs dont l'excentricité se trouve en dessous d'un seuil (proche de 0) et dont la taille se 2 2 trouve entre n ' rmin et n • rmax (étape 113c). Excentricités moyennes des CCs 3D On calcule ensuite la moyenne de l'excentricité de chaque CC 3D à partir des valeurs des CCs 2D qui en font partie (étape 113d). - Longueur des CCs 3D On garde la CC 3D dont l'excentricité moyenne est la plus petite et dont la longueur (sur l'axe z, c'est-à-dire selon l'axe du corps sur lequel l'acquisition a été réalisée) se trouve au-dessus d'un seuil. Cette composante doit correspondre à une partie de la trachée (étape 113e). Une fois détectée (si elle est présente), elle est soustraite du volume seuillé pour séparer les poumons de l'air dans la tête et le cou. - Sélection des poumons : utilisation de descripteurs sur les CCs 3D (étape 114) Plusieurs descripteurs sont utilisés pour bien choisir la ou les CCs correspondant aux poumons (ils représentent une composante connexe s'ils sont reliés à travers la trachée). Notamment : • Centre de masses : les poumons sont à peu près centrés sur l'image (alors que l'air autour du corps ne l'est pas) • Longueur : une longueur maximale est prise en compte (information anatomique) • Moyenne (p) et écart-type (a) : les deux poumons ont des statistiques similaires. Les deux valeurs sont calculées sur les intensités des voxels de chaque CC 3D • Ratio volume/longueur : les parties résiduelles d'air pouvant rester à l'intérieur du corps sont éliminées grâce à ce ratio (qui est très petit dans ce cas). On fait un premier parcours des composantes pour choisir des candidats en testant la position, l'écart-type (les poumons sont relativement homogènes), la moyenne (elle n'est jamais plus petite que -900 H.U. pour les poumons) et le ratio volume/longueur. Si un seul candidat a été trouvé, il est pris comme poumons. Par contre, si l'on en a choisi plusieurs, on compare si les coordonnées z de ses centres sont proches (les poumons sont situés à la même hauteur) et si les moyennes sont proches aussi afin de garder 2 CCs. Si aucune CC ne vérifie ces conditions, la CC prise est celle dont le ratio volume/longueur est le plus grand. - Choix de la première et de la dernière coupe des poumons (étape 115) Une fois que les poumons ont été identifiés, il faut choisir les valeurs en z de la fenêtre qui leur correspond. Le choix pour la limite inférieure est assez simple : il suffit de garder la valeur en z de la coupe contenant le dernier segment (le dernier voxel) des poumons. Le choix de la limite supérieure est fonction de l'information disponible sur les coupes 2D axiales, sur lesquelles on peut retrouver : • soit une partie de la trachée au-dessus des poumons et les deux poumons ; • soit une partie de la trachée au-dessus des poumons et un seul poumon ; • soit les deux poumons ; • soit un seul poumon. Lorsque la présence de deux poumons est détectée, on fixe comme début pour ceux-ci la première coupe où il apparaît deux composantes connexes 2D. Ceci donne de bons résultats dans les cas où la trachée est présente au-dessus des poumons : le nombre de CCs sera toujours égal à 1 jusqu'à arriver sur la première coupe des poumons. Par ailleurs, lorsqu'un seul poumon est détecté, le critère appliqué est le suivant : une coupe est considérée comme le début des poumons si : • la surface totale des CCs sur la coupe est plus grande qu'un seuil (typiquement plus grande que l'aire maximale de la trachée) ; • le nombre de CCs est plus grand que 1 et la condition précédente n'a pas été vérifiée. Cela permet de bien fixer le début lorsque la trachée est présente au-dessus des poumons et permet surtout d'éviter de fausses détections lorsque une partie seulement des poumons est contenue dans l'examen (les poumons peuvent être divisés en plus de 2 CCs à cause de la présence d'autres organes dans la partie inférieure). Les figures 5A à 5E illustrent différentes étapes du traitement d'identification des poumons. Les figures 5A et 5B montrent la trachée et les poumons superposés au volume. Les figures 5C à 5E montrent les résultats obtenus correspondant aux poumons. LA DETECTION DU COEUR La détection du coeur est basée sur la détection de l'aorte ascendante et l'aorte descendante dans la même coupe (laquelle se trouve juste au début du coeur et en dessous de l'arche aortique) grâce notamment à la forme circulaire qu'elles présentent. Cette coupe est considérée comme le début de la fenêtre cardiaque et permet de couper l'arche et donc d'améliorer les algorithmes cardiaques dont la performance est très impactée par sa présence. Les vaisseaux sont injectés avec un produit de contraste afin de les faire ressortir et les rendre assez visibles par rapport aux autres structures. Les différentes étapes d'un traitement de détection du coeur sont illustrées sur la figure 6, différentes images obtenues en mettant en oeuvre sur une image de coupe (figure 7A) différentes étapes d'un tel traitement étant représentées sur les figures 7B à 7G : Etapes - Filtrage gaussien (étape 131 ù figure 7B) Parfois les examens sont bruités et un seuillage peut générer des trous à l'intérieur des composantes connexes. Le filtre gaussien permet de réduire le bruit, mais il rend également l'image plus floue (ce qui rend les contours plus faibles, par exemple). Le filtre s'exprime : x2+y2 G(x,y) = 1 e 26-2 2~za-2 La valeur de uest choisie pour permettre d'enlever une partie du bruit sans trop dégrader l'image. - Seuillage du volume (étape 132 û figure 7C) Le volume est seuillé entre 120 H.U. et 750 H.U. afin d'isoler les vaisseaux. Néanmoins, d'autres structures sont présentes également. -Soustraction des contours (étape 133) 2921177 - 15 Dans la plupart des cas, le seuillage ne suffit pas pour séparer l'aorte de son entourage. Cependant, malgré la proximité en niveaux de gris, le gradient n'est pas négligeable et peut être utilisé. On extrait donc les contours pour créer un masque que l'on soustrait au volume seuillé. Ainsi, les structures dont 5 l'interface est plus marquée seront séparées les unes des autres. -Filtre de Sobel et seuillage du gradient (étape 133a) : Le gradient (figure 7D) est calculé par un masque de Sobel. Les valeurs du module sont plus petites à l'intérieur des objets car elles ont été atténuées par le filtrage (il y a moins de bruit). Les contours sont également un peu plus épais et 10 moins importants que si le filtrage n'avait pas été appliqué. On choisit un seuil pour le gradient tel que l'on ne garde que les contours afin de créer le masque de soustraction (la figure 7E correspondant à l'image du gradient après seuillage). Soustraction des contours sur le volume (étape 133b ù figures 7F) : 15 Le masque est ensuite soustrait sur l'image seuillée de façon que plusieurs structures soient décollées. Cela permet notamment de séparer l'aorte ascendante de l'artère pulmonaire, lesquelles se trouvent souvent assez collées. -Remplissage 2D (étape 134) En raison du seuillage et de la soustraction des contours, il est possible que 20 des trous apparaissent, donc il faut les reboucher. - Ouverture (étape 135 ù figure 7G) Après la soustraction des contours, il peut y avoir de petits morceaux encore collés à l'aorte. Cela rend la forme des aortes moins circulaire et pourrait alors perturber les mesures de l'excentricité. Une ouverture est appliquée dans le 25 but de produire un effet de lissage autour des structures. En pratique, la taille que l'on a choisie pour l'ouverture est : taille(cm) - Filtrage de taille 2D (étape 136) taille(voxels) = resolution (cm) Cette étape élimine les plus petites composantes qui restent dans l'image risquant d'affecter les étapes suivantes. La taille est fixée (en voxels) comme •r2. à partir du rayon minimal considéré pour l'aorte. resolution - Détection de la CC la plus circulaire pour chaque coupe (étape 137) Une fois ce prétraitement effectué, on analyse les CCs qui restent afin d'extraire et d'évaluer plusieurs descripteurs : • Taille ; • Position ; • Excentricité ; • Moyenne (p) et écart-type (a) (Les deux valeurs sont calculées sur les intensités des voxels de chaque CC) ; • Distance entre eux. On garde pour chaque coupe la composante dont l'excentricité est la plus petite, dont la variance est plus petite qu'un seuil (la densité du produit de contraste doit être relativement homogène) et dont la taille est plus petite 2 que r max (rayon maximal pour l'aorte). resolution - Détection de la deuxième CC qui correspond à l'aorte (étape 138) Après cela on parcourt le reste des CCs et on calcule la distance par rapport à la composante déjà enregistrée ainsi que la pente de la droite les reliant. On garde la deuxième composante pour chaque coupe si : • l'excentricité est en dessous d'un seuil (proche de 0) ; • la distance se trouve dans un intervalle fixe (pris à partir de critères anatomiques) ; • la pente se trouve dans un intervalle fixe (l'orientation du patient peut varier mais pas énormément) ; • les tailles sont similaires ; • les statistiques sont similaires. Cela permet de ne garder que l'aorte ascendante et l'aorte descendante pour chaque coupe, mais il peut y avoir des composantes isolées qui sont sélectionnées (par exemple sur une seule coupe). - Détection de trois coupes consécutives (étape 139) Afin d'éviter de fausses détections, on décide de fixer le début du coeur s'il y a trois coupes consécutives contenant les deux parties de l'aorte. Pour faire cela, on vérifie que les centres des CCs varient très peu entre les coupes, et la première des coupes est prise comme le début du coeur. - Choix de la limite inférieure du coeur (étape 140) La limite inférieure pour la fenêtre du coeur n'est pas critique puisque le fait d'avoir une surinformation en dessous du coeur n'a pas autant d'impact sur les performances des protocoles. Néanmoins, il faut fixer une limite en raison du temps de calcul. Comme le temps et les priorités n'ont pas permis de le faire, la limite est prise à partir de la limite supérieure du coeur (la première coupe du coeur) et en tenant compte des poumons : Z c e1.tr = (zmincoeur+16cm si(zm;ncoeur+l6cm) zmaxpoumons GENERATION AUTOMATIQUE D'UNE VUE DITE � CANDY CANE � (ARCHE DE L'AORTE' La détection de l'arche de l'aorte est d'un grand intérêt pour les médecins puisqu'elle permet de diagnostiquer facilement des maladies liées notamment à la dissection de l'aorte. Comme on le comprendra à la lecture de ce qui suit, le traitement qui est proposé pour détecter l'arche de l'aorte peut se mettre en oeuvre de façon totalement indépendante de l'arbre décisionnel des figures 3A et 3B. Ce traitement utilise l'identification de la première coupe de coeur et des centres de l'aorte ascendante et descendante dans cette coupe. La connaissance de ces centres permet de récupérer l'arche (en tant que composante 3D). Puis on extrait de l'arche une succession de points correspondants à différents centres de l'aorte ascendante et de l'aorte descendante, le nuage de points ainsi constitué étant utilisé pour calculer le plan médian de l'arche de l'aorte, plan dans lequel on obtient la vue dite � candy cane �. Plus précisément, les étapes de ce traitement sont les suivantes. -Détection des poumons (étape 110) Cette étape est facultative, mais elle permet de réduire le temps de calcul et de pouvoir traiter des volumes ayant une grande couverture en z. - Détection de la limite supérieure du coeur (étape 141 correspondant aux étapes 137 à 139) On garde les centres de l'aorte ascendante et descendante, lesquels feront partie des calculs du plan et servent à obtenir d'autres points. - Soustraction du gradient (étape 142) Une fois que l'on a détecté le début du coeur, l'arche (en tant que composante 3D) peut être collée aux artères pulmonaires au- dessus de la coupe que l'on a prise comme début du coeur. Pour éviter cela, on fait un seuillage du gradient moins restrictif (on garde plus de contours) et il est soustrait du volume. Bien que cette étape génère des trous à l'intérieur de l'aorte, ils ne sont pas gênants pour ce qui vient après car on préserve toujours des points sur l'aorte. - Ouverture à partir de la carte de distances 3D (étape 143) Il y a plusieurs ramifications qui partent de l'aorte et qui font donc partie de l'arche (de la CC 3D). Ces points peuvent décaler le plan calculé, donc ils sont enlevés grâce à une ouverture de taille n faite à partir d'une carte de distances 3D. Autrement dit, on calcule la carte de distances et on enlève les voxels dont la distance au fond est plus petite ou égale à n. Ensuite, on applique une dilatation avec le même principe. - Séparation de l'aorte du coeur (étape 144) L'aorte ascendante sort du coeur, donc il faut séparer celui-ci de l'arche. Comme l'aorte ascendante 2D avait déjà été trouvée sur la coupe du début du coeur, on peut la soustraire du volume et ainsi la séparer du coeur. - Récupération de l'arche 3D (étape 145) A part l'arche, il reste encore d'autres structures 3D. On sélectionne l'arche de la façon suivante : • récupération de la CC 3D (composante connexe en volume) contenant le point préalablement calculé sur l'aorte ascendante (étape 145a) : récupération de la partie ascendante de l'arche (ou de toute l'arche) ; • récupération de la CC 3D contenant le point préalablement calculé sur l'aorte descendante (étape 145b) : récupération de la partie descendante de l'arche (ou de toute l'arche) ; • union des deux composantes (étape 145c) (l'arche peut être parfois découpée en 2 CCs à cause des opérations des étapes précédentes). A la fin de cette étape, ce qui a été gardée contient soit toute l'arche (1 CC) soit les deux parties indépendantes (2 CCs). - Extraction de points sur l'arche (étape 146) Les n points sont extraits à partir d'un étiquetage 2D axial de l'arche en récupérant les centres de chaque composante 2D. - Calcul du plan au sens de moindres carrés (étape 147) Le système obtenu est surdéterminé. On minimise les distances des points au plan au sens de moindres carrés. Les figures 9A à 9F illustrent différents résultats issus du calcul de la vue dite � candy cane �. DETECTION DE L'ABDOMEN ET DES MEMBRES INFERIEURS La détection des anatomies se trouvant en dessous des poumons est basée sur la projection de la densité des os tout au long de l'axe z. Ainsi, on analyse la courbe tracée à partir de la surface d'os sur chaque coupe. Comme on le comprend en se référant aux figures 3A et 3B, d'un côté, si les poumons sont détectés (figure 3A), le profil est calculé à partir de la dernière coupe des poumons (partitionnement de l'étape 150b). De l'autre côté, s'ils ne sont pas présents (figure 3B), le profil contiendra au maximum l'abdomen et les jambes. Etapes - Soustraction de la table d'examen (étape 191) Il faut d'abord enlever la table d'examen parce que sa densité est proche de celle des os. Une fonction qui était déjà implémentée est utilisée pour le faire. -Seuillage du volume (étape 192) 292117T 20 Le volume est seuillé à partir de 120 H.U. pour garder les parties d'intérêt. - Remplissage 2D (étape 193) Il nous faut boucher certains trous qui apparaissent sur quelques parties osseuses afin d'éviter d'avoir des discontinuités. Cela peut arriver au niveau de la 5 hanche ou sur les disques de la colonne vertébrale. - Filtrage de taille 2D (étape 194) Filtrage de taille dans le but d'enlever des voxels isolés, de petits vaisseaux ou des parties d'organes ayant une densité au-dessus de 120 H.U. - Obtention du profil (étape 195) 10 On fait une somme de tous les voxels qui restent sur chaque coupe pour obtenir la surface selon : Surface(cm2) = n • res, (cm) • res, (cm) où n est le nombre de voxels et res, et resy les résolutions en x et en y. Dans le profil, l'axe y représente la surface et l'axe x représente la position 15 (soit en mm, soit en coupes). Néanmoins, la courbe obtenue est assez bruitée et irrégulière parce que l'on garde encore des parties qui ne sont pas des os et surtout à cause du remplissage. En effet, parfois le remplissage n'arrive pas à boucher les trous si les contours ne sont pas fermés, donc la surface va être plus petite dans ce cas. 20 La recherche et la caractérisation de maxima et minima sur une courbe bruitée peuvent être facilitées grâce à un filtrage gaussien 1D. Cela permet de ne garder que les valeurs extrêmes vraiment importantes et d'enlever le bruit. Après le filtrage, la courbe est plus parlante et représentative de l'anatomie. On peut identifier plusieurs zones d'intérêt (figure 11), à savoir : 25 • le fémur ; • le tibia ; • le pelvis ; • les genoux ; • les pieds. 30 D'un point de vue des protocoles de Volume Viewer toutes ces parties appartiennent à la même anatomie : Membres inférieurs. Par ailleurs, la hanche fait aussi partie de l'anatomie Abdomen. Le fémur et le tibia correspondent aux vallées de la courbe et logiquement le maximum entre les deux sont les genoux. Au bout de la courbe et après le tibia se trouvent les pieds (dernier maximum). Les deux maxima les plus proches à gauche du fémur font partie du pelvis. Effectivement, on en identifie un sur la partie supérieure de la hanche et l'autre sur la partie inférieure (juste sur la tête du fémur). - Détection du fémur et du tibia (étape 196) Les premières parties à détecter une fois que l'on a obtenu la courbe sont le fémur et le tibia. L'une des raisons pour faire cela est que l'on peut détecter les vallées de façon indépendante de la position relative sur le volume puisqu'elles n'apparaissent que sur le fémur et le tibia. Cela serait beaucoup plus difficile pour les maxima. D'ailleurs, on utilisera la position du fémur afin d'identifier les maxima correctement. Vis-à-vis de l'algorithme, une vallée est définie comme une région du profil où : • les valeurs sont en dessous d'un seuil (la surface est forcément plus petite) ; • elle a une longueur déterminée ; • les valeurs sont très proches les unes des autres tout au long de la vallée ; • il n'y pas de fortes variations entre des valeurs consécutives (la dérivée est petite). - Détection du pelvis et début des membres inférieurs (étape 197) La détection du fémur sert comme point de repère sur la partie inférieure du corps. Cette information est exploitable pour identifier le pelvis et fixer le début de l'anatomie Membres inférieurs. Effectivement, on va essayer de localiser les deux maxima à gauche des jambes en parcourant le profil dans le sens inverse à partir du début du fémur. Les conditions pour que les maxima soient retenus sont : • la valeur doit être au-dessus d'un seuil (la surface est plus grande que sur les jambes, par exemple) ; • la distance entre le maximum et le début du fémur doit être plus petite qu'une certaine distance maximale (utilisation de l'information anatomique). Il existe trois possibilités : retenir deux maxima, retenir un maximum et ne rien retenir. Dans le premier cas, on va établir le début des membres inférieurs sur le minimum le plus proche à gauche du dernier maximum (juste au-dessus du pelvis). Si l'on ne retient qu'un maximum, on agit analoguement au cas précédent parce que le maximum se trouvant plus en haut sur le pelvis est souvent plus robuste. Finalement, si aucun maximum n'été trouvé on considère que tout fait partie des membres inférieurs. Détection de l'abdomen (étape 198) A partir des donnés obtenues avec le profil des os, on ne fait pas directement un traitement pour détecter l'abdomen mais on prend plutôt une décision sur cette partie. Si l'on a bien identifié les jambes, le pelvis et le début des membres inférieurs, on peut dire sans risque d'erreur que toutes les coupes se trouvant au-dessus des membres inférieurs (les poumons n'y sont pas compris) appartiennent à l'abdomen. D'ailleurs, la limite inférieure de l'abdomen est fixée sur le début du fémur. Des résultats obtenus d'images d'abdomen ou de membres inférieurs sont là aussi illustrés sur les figures. DETECTION DE LA TETE ET DU COU La détection de la tête et du cou est basée sur l'air à l'intérieur de l'un et de l'autre. Etapes - Soustraction de la table d'examen (étape 180) La soustraction de la table est un point important pour le calcul de la boîte englobante puisqu'il est très simplifié si la seule structure 2D qui reste est la tête ou le cou. - Calcul des boîtes englobantes pour la tête et le cou (étape 181) Un seuillage est appliqué afin d'enlever l'air, puis on sélectionne sur chaque coupe la plus grande CC 2D (laquelle doit correspondre à la tête ou au cou). On parcourt ensuite les différentes CC obtenues et on garde pour chacune les couples (xmin,xmax) et (ymin,ymax) dont les combinaisons représentent les coins de la boîte. Les boîtes nous servent d'approximation afin de pouvoir déterminer si un objet se trouve à l'intérieur d'un autre objet (par exemple l'air et la tête). - Enregistrement de l'air à l'intérieur des boîtes (étape 182) Les boîtes ainsi calculées sont comparées à celles de l'air à l'intérieur du corps (tête/cou). On garde l'air avec un seuillage entre -1023 H.U. et -500 H.U. et on calcul pour chaque CC 2D son centre de masse et sa boîte englobante. Une composante est finalement retenue si : • sa boîte englobante est contenue dans la boîte englobante de la tête/cou ; • le centre de masse est à peu près centré sur l'image (si l'on considère que l'acquisition doit normalement contenir la plupart de la tête, l'air doit être centré). Les composantes qui restent après la sélection correspondent à l'air à l'intérieur (par exemple sinus ou trachée). -Obtention et analyse du profil d'air (étape 183) Un profil est obtenu de manière analogue à celui qui avait été calculé pour la densité des os. L'analyse de ce profil (figure 13) permet d'avoir un repère physique des sinus grâce à l'obtention du maximum de la courbe car logiquement ce maximum se trouve toujours au même endroit de la tête. Pour bien le détecter, un filtrage gaussien est appliqué sur la courbe originale. -Limite entre le cou et la tête (étape 184) Une fois le maximum détecté, on garde la valeur en z où la surface tombe jusqu'à 30%. Cette valeur en z est prise comme le début du cou (dont la limite inférieure est le début des poumons). Par ailleurs, la limite inférieure de la tête est fixée 2cm en dessous du début du cou (donc il y aura une petite partie commune aux deux anatomies). REVENDICATIONS 1. Procédé pour le traitement d'images d'anatomie en volume acquises par un système d'imagerie médicale, selon lequel on partitionne les images en sous-volumes correspondant à différentes parties d'anatomie identifiées sur lesdites images, caractérisé en ce que, pour réaliser ce partitionnement, on traite les images acquises pour identifier sur celles-ci les poumons et en ce que, si les poumons sont identifiés (110), on détermine (150a ; 150b) des coupes supérieure et inférieure délimitant ceux-ci selon l'axe de l'anatomie sur laquelle les images ont été acquises et on met en oeuvre un traitement distinct (160, 170, 180 ; 190, 200) sur les images correspondant à des zones d'anatomie respectivement au-dessus et en dessous de ces deux coupes. 2. Procédé selon la revendication 1, caractérisé en ce que, pour identifier les poumons sur les images, on met en oeuvre un traitement (133a à 113e) visant à soustraire l'air dans la trachée des images, selon lequel : • on détermine les différentes composantes connexes d'air qui apparaissent sur les images acquises, • on en détermine une excentricité moyenne 3D, • on compare à un seuil les excentricités moyennes 3D ainsi calculées pour ne conserver que celles de ces composantes connexes dont l'excentricité 3D est inférieure à un premier seuil et dont la longueur selon l'axe z est supérieure à un deuxième seuil, les composantes connexes ainsi isolées étant ensuite soustraites des images. 3. Procédé selon la revendication 2, caractérisé en ce qu'on met en oeuvre préalablement un seuillage (111), ainsi qu'un filtrage (112) sur la taille des composantes connexes. 4. Procédé selon l'une des revendications 1 à 3, caractérisé en ce qu'on sélectionne parmi différentes composantes connexes candidates une composante connexe correspondant aux poumons en fonction d'au moins un des descripteurs suivants : centre de masse, longueur, moyenne et écart-type, ratio volume/longueur. 5. Procédé selon l'une des revendications 1 à 4, caractérisé en ce que, pour déterminer la coupe correspondant à la limite inférieure des poumons, on parcourt les coupes axiales du bas de l'anatomie vers le haut et on détermine la première coupe sur laquelle apparaît un voxel des poumons. 6. Procédé selon l'une des revendications 1 à 4, caractérisé en ce que, pour déterminer la coupe correspondant à la limite supérieure des poumons lorsque la présence de deux poumons est détectée, on parcourt les coupes axiales du haut de l'anatomie vers le bas et on détermine la première coupe sur laquelle apparaissent deux composantes connexes. 7. Procédé selon l'une des revendications 1 à 4, caractérisé en ce que, pour déterminer la coupe correspondant à la limite supérieure des poumons lorsque la présence de deux poumons est détectée, on parcourt les coupes axiales du haut de l'anatomie vers le bas et on détermine la première coupe sur laquelle la surface totale des composantes connexes est supérieure à un seuil donné. 8. Procédé selon la revendication 7, caractérisé en ce que, en l'absence de coupe sur laquelle la surface totale des composantes connexes est supérieure audit seuil, on détermine la première coupe pour laquelle le nombre de composantes connexes est supérieur à 1. 9. Procédé selon l'une des revendications précédentes, caractérisé en ce que, lorsque des poumons sont détectés, on met en oeuvre un traitement (131 à 140) en vue d'identifier et de partitionner les images correspondant au coeur. 10. Procédé selon la revendication 9, caractérisé en ce qu'on compare préalablement la longueur de l'anatomie figurant sur les images à un seuil donné, ledit traitement (131 à 140) en vue d'identifier et de partitionner les images correspondant au coeur n'étant mis en oeuvre que lorsque ladite longueur est supérieure audit seuil. 11. Procédé selon l'une des revendications 9 ou 10, caractérisé en ce que, pour déterminer la coupe correspondant à la limite supérieure du coeur, on détermine au moins une coupe axiale sur laquelle on identifie l'aorte ascendante et l'aorte descendante. 12. Procédé selon la revendication 11, caractérisé en ce que, pour déterminer la coupe correspondant à la limite supérieure du coeur, on détermine au moins trois coupes axiales successives (139) sur lesquelles on identifie l'aorte ascendante et l'aorte descendante. 13. Procédé selon l'une des revendications 9 à 12, caractérisé en ce que, pour identifier une aorte sur une coupe axiale (138), on identifie les composantes connexes sur ladite coupe, on sélectionne parmi celles-ci la composante connexe dont l'excentricité est la plus petite et on compare la variance et la taille de cette composante respectivement à une variance seuil et un rayon maximal attendu pour l'aorte, cette composante n'étant considérée comme une aorte que si sa variance et sa taille sont inférieures à ladite variance et audit rayon maximal. 14. Procédé selon la revendication 13, caractérisé en ce que pour identifier une autre aorte sur une coupe sur laquelle une première composante connexe a déjà été identifiée comme correspondant à une aorte, on parcourt les composantes connexes de la coupe et on calcule leur distance par rapport à la première composante ainsi que la pente de la droite les reliant, et on sélectionne une composante connexe dont l'excentricité est en dessous d'un seuil donné, la distance et la pente déterminées pour cette composante connexe se trouvant dans des intervalles donnés. 15. Procédé selon la revendication 14, caractérisé en ce que, pour identifier une autre aorte, on compare la taille et la statistique des composantes connexes à des valeurs de taille et statistique fonction de celle de la composante connexe identifiée comme correspondant à la première aorte. 16. Procédé selon l'une des revendications 9 à 15, caractérisé en ce qu'on met en oeuvre préalablement sur les images un traitement comportant les étapes suivantes (132, 133) : • filtrage gaussien, • seuillage du volume, • soustraction des contours. 17. Procédé selon la revendication 16, caractérisé en ce que ledit traitement comporte également les étapes suivantes (142, 143): • ouverture, • filtrage 2D. 18. Procédé selon l'une des revendications précédentes, caractérisé en ce que, lorsque les poumons sont détectés, on détermine (195 à 198) un profil axial de surface occupée par les voxels sur les coupes en dessous de la coupe limite inférieure des poumons et on détermine la présence de 20 l'abdomen et de différentes parties des jambes en analysant ce profil. 19. Procédé selon l'une des revendications précédentes, caractérisé en ce que, lorsque les poumons ne sont pas détectés, on compare (210) la longueur de l'anatomie ayant fait l'objet de l'acquisition des images à un seuil et si cette longueur est supérieure audit seuil, on détermine un profil 25 axial de surface occupée par les voxels sur les coupes et on détermine (190, 200) la présence de l'abdomen et de différentes parties des jambes en analysant ce profil. 20. Procédé selon la revendication 19, caractérisé en ce que, lorsque l'analyse du profil ne permet pas de déterminer la présence de l'abdomen et de 30 différentes parties des jambes, on détermine (182, 183) un profil axial de 10 15l'air apparaissant sur les coupes et on analyse ce profil pour en déduire (180) une information sur la zone de la tête ou du cou présente sur les images. 21. Procédé selon l'une des revendications précédentes, caractérisé en ce que, lorsque les poumons sont détectés, on détermine (180 à 184) un profil axial de l'air apparaissant sur les coupes au- dessus de la coupe limite supérieure des poumons et on analyse ce profil pour en déduire une information sur la zone de la tête ou du cou présente sur les images. 22. Procédé selon la revendication 20, caractérisé en ce qu'on compare au préalable à un seuil donné la longueur d'anatomie au-dessus de la coupe supérieure limite des poumons, la détermination du profil axial de l'air n'étant mise en oeuvre que lorsque ladite longueur est supérieure audit seuil. 23. Procédé selon l'une des revendications précédentes, caractérisé en ce que, lorsque les poumons ne sont pas détectés, on compare la longueur de l'anatomie ayant fait l'objet de l'acquisition des images à un seuil et si cette longueur est inférieure audit seuil, on détermine un profil axial de l'air apparaissant sur les coupes et on analyse ce profil pour en déduire une information sur la zone de la tête ou du cou présente sur les images. 24. Procédé selon l'une des revendications précédentes, pour le traitement d'images en volume acquises par un système d'imagerie médicale, caractérisé en ce que (141 à 147): • on parcourt les coupes axiales perpendiculaires à l'axe de l'anatomie pour identifier une coupe axiale correspondant à 25 une limite supérieure du coeur, • on détermine sur cette coupe deux composantes connexes correspondant à l'aorte ascendante et à l'aorte descendante, • on détermine les composantes connexes en volume qui correspondent à ces composantes connexes, 30 • on fusionne lesdites composantes connexes en volume, 2921177 - 29 • on détermine sur la composante connexe en volume ainsi obtenue les centres des composantes connexes qui lui correspondent sur une pluralité de coupes axiales, • on détermine un plan optimisant une distance par rapport au 5 nuage de points constitués par ces différents centres. 25. Système d'imagerie médicale comportant des moyens (10) d'acquisition d'images d'anatomie en volume et des moyens de traitement de ces images aptes à partitionner les images en sous-volumes correspondant à différentes parties d'anatomie identifiées sur lesdites images, caractérisé 10 en ce que lesdits moyens (24, 36, 38) comportent des moyens aptes à mettre en oeuvre le procédé de traitement selon l'une des revendications précédentes. 26. Programme mémorisé sur des moyens aptes à être lus par un ordinateur, caractérisé en ce qu'il comporte des données et/ou instructions pour 15 mettre en oeuvre un procédé selon l'une des revendications 1 à 24.