Brevet FR2869438A1 — publié le 20051028

Procédé de segmentation automatique des cavités cardiaques La présente invention a pour objet un procédé de segmentation automatique des cavités cardiaques. Le domaine de l'invention de celui de l'imagerie numérique, et plus particulièrement celui de l'imagerie numérique médicale. Les systèmes d'imagerie destinés au diagnostic médical comprennent une variété de types d'imagerie, tels que les systèmes à rayon X, les systèmes de tomodensitométrie (CT), les systèmes à ultrasons, la tomographie à faisceau d'électrons (EBT), les systèmes à résonance magnétique (MR), et similaire. Les systèmes d'imagerie destinés au diagnostic médical génèrent des images d'un objet, tel qu'un patient, par exemple, en l'exposant à une source d'énergie, comme par exemple, les rayons X traversant un patient. Les images générées peuvent être utilisées à plusieurs fins. On peut détecter, par exemple, les défauts internes d'un objet. En outre, on peut déterminer les changements dans une structure ou un alignement internes. On peut représenter également les flux de fluide dans un objet. En outre, l'image peut montrer la présence ou l'absence d'objets dans un objet. L'information obtenue à partir d'une imagerie de diagnostic médical a des applications dans plusieurs domaines comprenant la médecine et la fabrication. Grâce au progrès réalisés dans les techniques de reconstruction d'images synchronisées par électrocardiogramme, l'imagerie par tomodensitométrie fournit maintenant aux radiologistes des vues fiables de l'anatomie cardiaque. Les artères coronaires, ainsi que les cavités droite et gauche du coeur sont rendues visibles grâce à l'injection intraveineuse d'un produit de contraste. Le ventricule gauche et l'oreillette gauche sont séparées par la valve mitrale. Le ventricule gauche est connecté à l'aorte ascendante via la valve aortique. Plusieurs veines pulmonaires (généralement de 3 à 6) sont connectées à l'oreillette gauche. Les cavités gauches sont définies comme l'union de l'oreillette gauche et du ventricule gauche. Elles n'incluent pas l'aorte ascendante, mais incluent le début des veines pulmonaires. Une technique de segmentation 3D automatique des cavités gauches a plusieurs applications importantes: - Facilitation de la visualisation 3D. La compréhension de l'anatomie du ventricule, de l'oreillette, et de I'ostium des veines pulmonaires est beaucoup plus facile en trois dimensions (3D) qu'en deux dimensions (2D). De plus une vue non segmentée des cavités gauches montre beaucoup d'autres structures (côtes, poumons, cavité droite, ...) qui nuisent à la lisibilité de l'anatomie des cavités gauches en masquant leurs surfaces. - La connaissance des contours de l'oreillette gauche rend les opérations d'ablation de veines pulmonaires plus sûres. - La segmentation du ventricule gauche peut être facilement déduite de la segmentation des cavités gauches en laissant l'utilisateur définir à la 10 main la valve mitrale. Dans l'état de la technique on obtient de bonnes images volumiques du muscle cardiaques et de son environnement veineux et artériel grâce à l'injection d'un produit de contraste dans le système sanguin du patient. Cependant l'image résultant d'un tel examen radiologique, dont le but est de visualiser le coeur du patient, fait aussi apparaître de nombreux artefacts qui perturbent la lecture de l'image volumique. Une telle image volumique comporte en effet, en plus du muscle cardiaque, des os (côtes, colonne vertébrale, ...), des veines et artères, et du poumon, pour ne citer que les plus facilement identifiables. Dans la mesure ou seul le coeur intéresse le praticien, une telle abondance d'information est nuisible. L'invention résout ces problèmes en réalisant une segmentation automatique des cavités, gauches ou droites, du muscle cardiaque. Une fois décidées quelles cavités doivent être isolées, l'invention comporte essentiellement trois étapes. Une première étape dans laquelle on détermine, dans l'image volumique issue de l'examen, quel est le volume comportant les cavités cardiaques. Une deuxième étape dans laquelle on effectue une identification des cavités gauches et droites. Et une troisième étape dans laquelle on effectue une reconstruction précise des contours des cavités, gauches ou droites, dont on veut réaliser la segmentation. Si on souhaite réaliser une segmentation des cavités gauches, alors le procédé comporte une quatrième étape optionnelle permettant de retirer l'aorte des cavités gauches, permettant ainsi d'obtenir une image du muscle cardiaque seul. L'invention a donc pour objet un procédé de segmentation d'une partie 35 du coeur dans une première image volumique correspondant à une zone du corps comportant le coeur, caractérisé en ce que le procédé comporte des étapes mises en oeuvre par une unité de traitement: - une première étape d'isolation du coeur dans la première image volumique, cette première étape comportant les étapes suivantes: - une étape de premier seuillage de la première image volumique avec un niveau prédéfini de manière à isoler les éléments les plus clairs de l'image volumique, cette étape produit une deuxième image volumique, - une étape de réalisation d'une première érosion de la deuxième image volumique, cette étape produit une troisième image volumique, - une étape de sélection dans la troisième image volumique de la plus grande composante connexe, cette composante correspondant alors au coeur, cette étape produit une quatrième image volumique, - une deuxième étape de séparation, dans la quatrième image volumique, les cavités gauches des cavités droites du coeur, cette deuxième étape comportant les étapes suivantes: - une étape de deuxième seuillage de la quatrième image volumique avec un niveau prédéfini plus élevé que pour le premier seuillage, - une étape dans laquelle on effectue des érosions successives de la quatrième image volumique seuillée jusqu'à obtenir une image comportant deux composantes connexes principales, en taille, dont le rapport des tailles est compris dans un intervalle prédéfini, cette étape produit une cinquième image volumique, - une étape de d'identification des composantes connexes principales dans la cinquième image volumique, - une troisième étape de reconstruction des contours exacts des composantes connexes principales, cette troisième étape comportant les étapes suivantes: - une étape de calcul du gradient de la quatrième image volumique, cette étape produit une sixième image volumique, - une étape d'application de l'algorithme de ligne de partage des eaux sur la sixième image volumique, cette application se faisant en utilisant comme graines pour l'algorithme les composantes connexes principales identifiées, cette étape produit une septième image volumique, - une étape de sélection, dans la septième image volumique, de l'élément issu du traitement par l'algorithme de ligne de partage des eaux correspondant à la graine ayant été identifiée comme faisant partie des cavités qu'un praticien souhaite isoler, cette sélection permettant d'extraire, de la troisième image volumique, les voxels correspondant aux cavités que le praticien souhaite isoler, cette étape produit une huitième image volumique des cavités isolées, - une étape de dilatation sur la huitième image, cette dilatation ayant une taille égale a celle de la première érosion, cette étape produit une neuvième image volumique. Avantageusement l'invention est aussi caractérisée en ce que les cavités que le praticien souhaite isoler sont les cavités gauches, le procédé comportant alors les étapes suivantes: - une étape de réalisation d'une troisième érosion forte sur la neuvième volumique, cette étape produit une dixième image volumique, - une étape d'identification de 2 points distincts appartenant respectivement a l'aorte et au ventricule gauche dans le nuages de voxels résultant de la troisième érosion, dans ce nuage de voxels, celui situé le plus haut correspondant a l'aorte, tandis que le voxel situé le plus bas correspond au ventricule gauche, - une étape de calcul du gradient de la neuvième image volumique, cette étape produisant une onzième image volumique, - une étape d'application de l'algorithme de ligne de partage des eaux sur la onzième image volumique, cette application se faisant en utilisant comme graine pour l'algorithme les deux voxels identifies à l'étape d'identification, cette étape produisant une douzième image volumique, - une étape de sélection, dans la douzième image volumique, de l'élément issu du traitement par l'algorithme de ligne de partage des eaux correspondant à la graine ayant été identifiée comme faisant partie des cavités gauche, cette sélection permettant d'extraire, de la première image volumique, les voxels correspondant aux cavités gauches, cette étape produisant une treizième image volumique des cavités gauches sans l'aorte. Avantageusement l'invention est aussi caractérisée en ce que le seuil du premier seuillage est de 130 sur une échelle allant de -1024 à 4096. Avantageusement l'invention est aussi caractérisée en ce que, pour les érosions, les éléments structurant sont des croix. Avantageusement l'invention est aussi caractérisée en ce que la première érosion est une érosion de l'ordre de 2 voxels. Avantageusement l'invention est aussi caractérisée en ce que les érosions successives sont effectuées par un élément structurant permettant une érosion de l'ordre de 1 voxel. Avantageusement l'invention est aussi caractérisée en ce que le rapport de la taille de la plus grande des composantes connexes principales sur la taille de la plus petite des composantes connexes principales est compris dans l'intervalle [1, 10]. Avantageusement l'invention est aussi caractérisée en ce que le seuil du deuxième seuillage est de 160 sur une échelle allant de -1024 à 4096. Avantageusement l'invention est aussi caractérisée en ce que l'identification des composantes connexes principales se fait en calculant au moins le centre de gravité de l'une des composantes connexes principales, une droite horizontale, lorsque l'on est face au patient, passant par ce centre de gravité intercepte alors en premier les cavités droites. Avantageusement l'invention est aussi caractérisée en ce que les cavités cardiaques sélectionnées sont les cavités gauches. Avantageusement l'invention est aussi caractérisée en ce que la troisième érosion est une érosion de l'ordre de 10 voxels. L'invention sera mieux comprise à la lecture de la description qui suit et à l'examen des figures qui l'accompagnent. Celles-ci sont présentées à titre indicatif et nullement limitatif de l'invention. Les figures montrent: - Figure 1: une illustration d'étapes du procédé de segmentation selon l'invention; - Figure 2: une illustration d'étapes du procédé de segmentation selon l'invention permettant la séparation des cavités gauches et de l'aorte; - Figure 3: une illustration d'un élément structurant permettant la réalisation des érosions; - Figure 4: une illustration d'une image volumique issue d'un examen après la réalisation d'un premier seuillage selon l'invention; - Figure 5: une illustration d'une image volumique issue de la première étape du procédé selon l'invention; - Figure 6: une illustration des cavités gauches segmentées par le procédé selon l'invention; - Figure 7: une illustration des cavités gauches segmentées et sans l'aorte obtenues par le procédé selon l'invention. Dans la présente invention les étapes de segmentation automatique sont mises en oeuvre par une unité de traitement appartenant, ou étant connecté, à un appareil radiologique de type scanner. Une telle unité de traitement comporte au moins un microprocesseur et des mémoires. Les étapes sont donc mises en oeuvre par ce microprocesseur commandé par des codes instructions enregistrés dans les dites mémoires. Préalablement à la mise en oeuvre de l'invention un praticien procède à un examen radiologique produisant une première image volumique de la zone du corps comportant le coeur d'un patient. Pour une efficacité maximale du procédé selon l'invention cette première image volumique est prise à 70% du cycle cardiaque, ce qui correspond à la diastole. Une image volumique est un ensemble, ou nuage, de points aussi appelés voxels. Chaque voxel est caractérisé par ses coordonnées dans l'espace et une valeur qui lui est attribué par le scanner. Cette valeur est exprimée en HU pour Houndsfield Unit. Le HU est une unité proportionnelle au coefficient d'absorption des rayons X par les tissus. Chaque voxel est donc associé à une valeur exprimée en HU, la dynamique d'une telle valeur étant [-1024, 3072]. Les valeurs les plus faibles correspondant aux tissus les moins absorbant. Les cavités sont artificiellement rendues plus absorbantes par l'utilisation d'un produit de contraste injecté dans le sang. La première image volumique issue de l'examen comporte beaucoup d'information en plus de celles qui intéressent le praticien. Dans le cas de l'invention, seul l'information, donc les voxels, correspondant aux cavités gauches ou droites du coeur intéressent le praticien. En sortie d'examen, la première image volumique comporte des informations concernant non seulement le coeur, mais aussi les veines, artères, os, poumons et autres organes situées à proximité. Dans cette description, toutes les images volumiques ont les même dimensions spatiales et orientation, et sont donc superposables. La figure 1 montre une première étape 101 produisant, à partir de la première image volumique, une quatrième image volumique ne comportant plus que des informations relatives aux cavités cardiaques. C'est à dire que cette quatrième image volumique permet répartir les voxels de la première image volumique en deux catégories. Une première catégorie correspondant aux voxels appartenant aux cavités cardiaques, et une deuxième catégorie comportant tous les voxels de la première image volumique n'étant pas dans la première catégorie. L'étape 101 se subdivise elle-même en 3 étapes. Dans une étape 102 l'unité de traitement effectue un seuillage de la première image volumique en ne conservant que les voxels dont la valeur associée est supérieure à un premier seuil. Dans un exemple préféré, ce premier seuil vaut 130 HU. En d'autres termes on considère ici que tous les voxels ayant une valeur associée supérieure à 130 HU ont une probabilité non nulle d'appartenir au coeur. En procédant ainsi on isole les éléments les plus clairs de l'image volumique. La figure 4 illustre le résultat de cette première étape qui est une deuxième image volumique. Cette deuxième image volumique comporte encore de nombreuse structure en plus du coeur, notamment les os qui absorbent bien les rayons X. Après l'étape 102, l'unité de traitement réalise une étape 103 d'érosion de la deuxième image volumique. Une érosion est une opération de morphologie mathématique connue consistant à traiter chaque voxel par un élément structurant. On centre l'élément structurant sur la voxel à traiter, et si l'élément structurant ainsi placé comporte des voxels n'appartenant pas à la structure que l'on cherche à éroder, alors, le voxel est supprimé de l'image volumique résultat de l'érosion. Par supprimer on entend qu'il est éteint, c'est à dire que la valeur qui lui est associée vaut O. Pour l'invention, l'érosion est isotrope. Ce résultat, d'érosion isotrope, est obtenu, par exemple, par l'utilisation d'un élément structurant en croix. La figure 3 montre un tel élément 300 structurant. Un tel élément structurant comporte 3 branches orthogonales d'une longueur de 3 voxels chacune et se coupant en leur milieu. Les structures à éroder sont celles présentes dans la deuxième image volumique. Le résultat de l'étape 103 est une troisième image volumique comportant un certain nombre de structures dissociées les unes des autres, l'érosion de l'étape 103 ayant contribué à dissocier les structures faiblement liées. Dans une mise en oeuvre préférée l'érosion de l'étape 103 est de 2 voxels. Pour la réaliser on effectue en fait deux érosions successives de 1 voxel avec l'élément structurant précédemment décrit. Les structures de la troisième image volumique sont aussi appelées des composantes connexes. Une composante connexe est un ensemble de voxels, chaque voxel de l'ensemble ayant pour voisin au moins un autre voxel dudit ensemble. Dans l'invention, un tel ensemble est borné par des points supprimés. Dans une étape 104 de sélection, l'unité de traitement associe à chaque composante connexe une taille. Une telle taille est, par exemple, le nombre de voxels que comporte la composante connexe. Une fois la taille des composantes connexes de la 3eme image volumique déterminée, l'unité de traitement sélectionne la plus grande d'entre elle et supprime toutes les autres, produisant ainsi une quatrième image volumique ne comportant que des voxels appartenant aux cavités cardiaques. Une telle image volumique est illustrée par la figure 5. Une telle image volumique montre les ventricules gauche et droit, les oreillettes gauche et droite, les veines pulmonaires et l'aorte. L'étape 104 est la dernière étape de l'étape 101. La figure 1 montre que l'étape 101 est suivie d'une étape 105 de séparation des cavités gauches et droites du coeur. L'étape 105 comporte elle aussi plusieurs étapes. La figure 1 montre une étape 106, suivant l'étape 104, dans laquelle l'unité de traitement effectue un deuxième seuillage dans la quatrième image volumique. Ce deuxième seuillage est effectué avec un seuil supérieur à celui utilisé pour le premier seuillage et valant, dans un exemple préféré, 160. Ce seuillage produit alors une cinquième image volumique initiale. Dans une étape 107 suivant l'étape106, l'unité de traitement effectue des érosions successives en partant de la cinquième image volumique, chaque érosion s'effectuant sur l'image volumique résultant de la précédente érosion. Un des intérêts des érosions successives est que l'on évite ainsi d'endommager plus que nécessaire les structures à séparer. Les érosions successives sont réalisées tant qu'un critère de fin d'itération n'est pas atteint. Dans un exemple préféré ce critère est le suivant pour chaque image volumique issue d'une érosion, on mesure la taille des composantes connexes présentes dans l'image volumique. Puis on calcule le rapport des deux plus grandes composantes connexes. Si ce rapport est compris dans un intervalle prédéterminé, alors on interrompt les érosions et on ne conserve, dans une cinquième image volumique définitive, que les deux plus grandes composantes connexes dont on vient de faire le rapport des tailles. La cinquième image volumique définitive est celle issue des érosions successives et est désignée par la suite comme la cinquième image volumique. Dans un exemple préféré, on calcule le rapport de la taille de la plus grande composante connexe sur la taille de la seconde plus grande composante connexe, et on interrompt les érosions successives lorsque ce rapport est inférieur à 5. Dans une variante, ce rapport doit être inférieur à 10, ou à un nombre compris dans l'intervalle [1, 10]. Une fois le critère de fin d'itération atteint, l'unité de traitement supprime toutes les composantes qui ne sont pas les deux plus grandes composantes connexes. On obtient alors une cinquième image volumique ne comportant plus que deux composantes connexes. De l'étape 107 on passe à une étape 108 d'identification des composantes connexes sélectionnées à l'étape 107. Dans cette étape on considère la cinquième image volumique orientée. Cette orientation est celle issue de l'examen ayant produit la première image volumique. Ainsi le haut d'une image volumique correspond à la partie supérieure du corps à partir duquel l'image volumique a été produite, et ce lorsque la personne à laquelle appartient le corps se tient debout. La gauche d'une image volumique correspond à la partie gauche du corps à partir duquel l'image volumique a été produite, et ce lorsque l'on regarde de face la personne à laquelle appartient ce corps. Ces deux indications permettent de totalement orienter une image volumique. Dans l'étape 107, et dans une variante préférée, l'unité de traitement calcule le centre de gravité d'une des deux composantes connexes restant à la fin de l'étape 107. Puis l'unité de traitement considère une droite horizontale passant par ce centre de gravité et parcourt cette droite en allant de la gauche vers la droite, la gauche et la droite étant celle d'une personne qui fait face au patient. Lors de ce parcours, la première composante connexe rencontrée, parmi, les deux restantes à la fin de l'étape 107, correspond aux cavités gauches du coeur. Dans une autre variante, à l'étape 107, l'unité de traitement calcule les centres de gravité des deux composantes connexes issues de l'étape 107. Les positions relatives de ces deux centres de gravités permettent alors l'identification des composantes connexes, donc des cavités cardiaques. Le centre de gravité situé le plus haut correspond à la composante connexe appartenant aux cavités droites du coeur. A la fin de l'étape 108, les deux composantes connexes issues de 10 l'étape 107 sont chacune identifiée comme correspondant soit aux cavités gauches, soit aux cavités droites. L'étape 105 est suivie d'une étape 109 de reconstruction des contours des cavités cardiaques. L'étape 105 comporte elle-même plusieurs étapes. La figure 1 montre que l'étape 108 est suivie d'une étape 110 de calcul de l'image gradient de la quatrième image volumique. Cette image gradient est une sixième image volumique. La sixième image volumique est obtenue à partir de la quatrième par l'application d'un filtre gradient. L'étape 110 est suivie par une étape 111 d'application de l'algorithme de partage des eaux sur la sixième image volumique. Cet algorithme est appliqué en utilisant comme graines les deux composantes connexes issues de l'étape 108. Une graine correspond donc aux cavités gauches, l'autre graine correspondant aux cavités droites. L'algorithme de partage des eaux à partir de graines peut être vu comme une simulation d'inondation. Le nom "partage des eaux" vient de la manière dont l'algorithme segmente les régions d'un volume en bassins. Ces bassins sont, au départ, les points de faibles intensités dans le volume devant être segmenté. Ces basins partagent des frontières les uns avec les autres. Si on imagine qu'il pleut dans ces basins, le niveau d'eau va monter. Au fur et à mesure que le niveau d'eau monte, les basins se remplissent et finissent par se rejoindre au-dessus de leurs frontières pour former des bassins de plus en plus grand. Dans l'invention l'image qui est segmentée est une image gradient, ce qui a pour effet de renforcer les contours des objets présents dans l'image, tout en creusant l'intérieur de ces objets. Les graines utilisées dans l'invention correspondent à l'intérieur des objets à segmenter, et c'est à partir de ces graines que l'on procède à l'inondation. En d'autres termes ces graines sont les premiers basins, chacun de ces deux premiers basins étant associé à un label. Au cours de l'inondation, tout nouveau bassin entrant en contact avec l'un de ces deux premiers basins acquiert le label de ce basin. Le processus d'inondation, donc l'algorithme de partage des eaux, s'arrête lorsque chaque voxel s'est vu attribuer un label. Dans l'invention l'inondation est bornée par la composante connexe de la quatrième image volumique. C'est à dire que le processus d'inondation n'affecte un label qu'aux voxels de cette composante connexe. A la fin de l'étape 111, tous les voxels de la composante connexe de la quatrième image volumique se sont vu attribué un label, via la segmentation de la sixième image volumique. Donc à la fin de l'étape 111 un voxel non supprimé de la quatrième image volumique appartient soit aux cavités droites, soit aux cavités gauches, du coeur. On passe alors à une étape 112 de sélection dans laquelle dans laquelle on supprime tous les voxels correspondant à l'un uniquement des labels selon que l'on souhaite obtenir une image des cavités gauches ou droites. Les coordonnées des voxels non supprimés sont connues. Cette connaissance permet d'extraire de la troisième image volumique les informations d'intensité correspondant à ces voxels et ainsi de produire une huitième image volumique. A la fin de l'étape 112 on dispose alors de la huitième image volumique ne comportant que des voxels correspondant aux cavités cardiaques que l'on souhaite visualiser. On dit aussi que l'on a isolé les cavités cardiaques, qu'elles soient gauches ou droites. Cette sélection est réalisée en fonction d'un paramétrage effectué par un praticien préalablement à l'étape 101. Ce paramétrage consiste à désigner les cavités, gauches ou droites, que le praticien souhaite isoler. L'unité de traitement utilise alors cette désignation, gauche ou droite, comme critère de sélection à l'étape 112. Dans une variante préférée, l'étape 112 est suivie d'une étape 113 de dilatation. Une dilatation est l'opération inverse d'une érosion. La dilatation appliquée à l'étape 113 est de taille égale à la taille de l'érosion réalisée à l'étape 103. Dans cette variante, l'unité de traitement sauvegarde donc les paramètres de la première érosion de l'étape 103 dans une mémoire. Cette dilatation permet de récupérer l'information qui avait été supprimée à l'étape 103 dans un souci de robustesse. Cette dilatation se fait donc avec le même élément structurant que celui utilisé pour l'étape 103, et comporte le même nombre de passes. L'érosion de l'étape 103 étant dans l'exemple préférée composée de deux érosions de 1 voxel, la dilatation de l'étape 113 est composée de deux dilatations de 1 voxel. Cette dilatation permet donc de récupérer, dans la première image volumique les voxels voisin, selon la dilatation, de la structure isolée dans la huitième image volumique. Dans la variante préférée, la fin de l'étape 113 est aussi la fin de l'étape 109. A la fin de cette étape 109, un praticien dispose donc d'une image volumique, désigné dans cette description comme étant la neuvième, qui ne comporte plus que des voxels correspondant aux cavités gauches du coeur. En effet dans notre exemple on considère que le praticien agissant sur l'unité de traitement souhaite isoler les cavités gauches du coeur et a donc, à l'étape 112, sélectionné le label correspondant aux cavités gauches. La neuvième image volumique est illustrée par la figure 6. La figure 6 montre une image volumique comportant des informations relatives aux oreillette et ventricule gauches, aux veines pulmonaires, et à l'aorte. Si le praticien avait sélectionné le label correspondant aux cavités droites, alors le procédé serait terminé. Cependant dans le cas des cavités gauches l'étape 109 est suivie d'une étape 200 de suppression de l'aorte. L'étape 200 comporte elle-même plusieurs étapes. La figure 2 montre que l'étape 200 comporte une première étape 201 de réalisation d'une troisième érosion forte. Cette troisième érosion est réalisée sur la neuvième image volumique. Cette troisième érosion forte est une érosion de 10 voxels. Comme précédemment cette érosion forte est en fait réalisée par 10 petites érosions de 1 voxel. Le résultat de cette érosion forte est une dixième image volumique. La dixième image volumique comporte un nuages de voxels composé des voxels qui n'ont pas été supprimés par l'érosion forte. L'étape 201 est suivie d'une étape 202 d'identification de 2 voxels distincts appartenant respectivement a l'aorte et au ventricule gauche dans la dixième image volumique. Cette identification se fait grâce à l'orientation de l'image volumique. Le voxel du nuage de voxels de la dixième image volumique situé le plus haut et le plus proche de l'avant du patient, c'est à dire de ses cotes, appartient a l'aorte, tandis que le voxel le plus bas appartient au ventricule gauche. L'étape 201 est aussi suivie d'une étape 203 dans laquelle l'unité de traitement calcul l'image gradient de la neuvième de la neuvième image volumique. Le résultat de l'étape 203 est une onzième image volumique. Une fois les étapes 202 et 203 effectuées, l'unité de traitement passe à une étape 204 de mise en oeuvre de l'algorithme de partage des eaux sur la onzième image volumique en utilisant comme graines les voxels identifiés à l'étape 202. La mise en oeuvre de l'étape 204 est identique à la mise en oeuvre de l'étape 111, si ce n'est que l'on utilise pas la même image de départ, ni les même graines. Le résultat de l'étape 204 est une douzième image volumique correspondant à la neuvième image volumique segmentée en deux parties. Une partie correspondant à l'aorte, l'autre partie correspondant aux ventricule et oreillette gauches. A la fin de l'étape 204, chacun des voxels non supprimés de la onzième image volumique appartient soit à l'aorte, soit aux ventricule et oreillette gauches. L'étape 204 est suivie d'une étape 205 de sélection identique à l'étape 112, si ce n'est que le critère de sélection est cette fois de n'extraire que les voxels correspondant aux ventricule et oreillette gauches. Dans l'étape 205 la partie correspondant à l'aorte est donc supprimée ce qui permet de produire une treizième image volumique, illustrée par la figure 7, ne comportant que les ventricule et oreillette gauches. Cette treizième image volumique est particulièrement utile dans le domaine de l'électrophysiologie. Dans ce domaine, chez certaine personne, il apparaît un signal perturbateur naissant dans les veines pulmonaires. Ce signal entraîne une arythmie qui doit être corrigée par une ablation desdites veines pulmonaires. La treizième image volumique est particulièrement bien adaptée à la préparation de cette ablation car elle ne conserve que les parties du coeur concernées par l'opération. Cette invention présente aussi un intérêt dans la mesure des débits sanguin, et en particulier dans la mesure de la fraction d'éjection des ventricules. 14 REVENDICATIONS 1 - Procédé de segmentation d'une partie du coeur dans une première image volumique correspondant à une zone du corps comportant le coeur, caractérisé en ce que le procédé comporte des étapes mises en oeuvre par une unité de traitement: - une première étape (101) d'isolation du coeur dans la première image volumique, cette première étape comportant les étapes suivantes: - une étape de (102) premier seuillage de la première image volumique avec un niveau prédéfini de manière à isoler les éléments les plus clairs de l'image volumique, cette étape produit une deuxième image volumique, - une étape (103) de réalisation d'une première érosion de la deuxième image volumique, cette étape produit une troisième image 15 volumique, - une étape (104) de sélection dans la troisième image volumique de la plus grande composant connexe, cette composante correspondant alors au coeur, cette étape produit une quatrième image volumique, - une deuxième étape (105) de séparation, dans la quatrième image volumique, les cavités gauches des cavités droites du coeur, cette deuxième étape comportant les étapes suivantes: - une étape (106) de deuxième seuillage de la quatrième image volumique avec un niveau prédéfini plus élevé que pour le premier 25 seuillage, - une étape (107) dans laquelle on effectue des érosions successives de la quatrième image volumique seuillée jusqu'à obtenir une image comportant deux composantes connexes principales, en taille, dont le rapport des tailles est compris dans un intervalle prédéfini, cette étape produit une cinquième image volumique, - une étape (108) de d'identification des composantes connexes principales dans la cinquième image volumique, - une troisième étape (109) de reconstruction des contours exacts des composantes connexes principales, cette troisième étape comportant les 35 étapes suivantes: une étape (110) de calcul du gradient de la quatrième image volumique, cette étape produit une sixième image volumique, - une étape (111) d'application de l'algorithme de ligne de partage des eaux sur la sixième image volumique, cette application se faisant en utilisant comme graines pour l'algorithme les composantes connexes principales identifiées, cette étape produit une septième image volumique, - une étape (112) de sélection, dans la septième image volumique, de l'élément issu du traitement par l'algorithme de ligne de partage des eaux correspondant à la graine ayant été identifiée comme faisant partie des cavités qu'un praticien souhaite isoler, cette sélection permettant d'extraire, de la troisième image volumique, les voxels correspondant aux cavités que le praticien souhaite isoler, cette étape produit une huitième image volumique des cavités isolées, - une étape (113) de dilatation sur la huitième image, cette dilatation ayant une taille égale à celle de la première érosion, cette étape produit une neuvième image volumique. 2 - Procédé selon la revendication 1, caractérisé en ce que les cavités que le praticien souhaite isoler sont les cavités gauches, le procédé 20 comportant alors les étapes suivantes: - une étape (201) de réalisation d'une troisième érosion forte sur la neuvième volumique, cette étape produit une dixième image volumique, - une étape (202) d'identification de 2 voxels distincts appartenant respectivement a l'aorte et au ventricule gauche dans le nuages de voxels résultant de la troisième érosion, dans ce nuage de voxels, celui situé le plus haut correspondant a l'aorte, tandis que le voxel situé le plus bas correspond au ventricule gauche, - une étape (203) de calcul du gradient de la neuvième image volumique, cette étape produisant une onzième image volumique, - une étape (204) d'application de l'algorithme de ligne de partage des eaux sur la onzième image volumique, cette application se faisant en utilisant comme graine pour l'algorithme les 2 voxels identifiés à l'étape d'identification, cette étape produisant une douzième image volumique, - une étape (205) de sélection, dans la douzième image volumique, de l'élément issu du traitement par l'algorithme de ligne de partage des eaux correspondant à la graine ayant été identifiée comme faisant partie des cavités gauche, cette sélection permettant d'extraire, de la première image volumique, les voxels correspondant aux cavités gauches, cette étape produisant une treizième image volumique des cavités gauches sans l'aorte. 3 - Procédé selon l'une des revendications 1 ou 2 caractérisé en ce que le seuil du premier seuillage est de 130 sur une échelle allant de -1024 à 4096. 4 - Procédé selon l'une des revendications 1 à 3, caractérisé en ce que, pour les érosions, les éléments (300) structurant sont des croix. 5 - Procédé selon l'une des revendications 1 à 4, caractérisé en ce que la première érosion est une érosion de l'ordre de 2 voxels. 6 - Procédé selon l'une des revendications 1 à 5, caractérisé en ce que les érosions successives sont effectuées par un élément structurant permettant une érosion de l'ordre de 1 voxel. 7 - Procédé selon l'une des revendications 1 à 6, caractérisé en ce que le rapport de la taille de la plus grande des composantes connexes principales sur la taille de la plus petite des composantes connexes principales est compris dans l'intervalle [1, 10]. 8 - Procédé selon l'une des revendications 1 à 7, caractérisé en ce que le seuil pour le deuxième seuillage est de 160 sur une échelle allant de 1024 à 4096. 9 - Procédé selon l'une des revendications 1 à 8, caractérisé en ce que l'identification des composantes connexes principales se fait en calculant au moins le centre de gravité de l'une des composantes connexes principales, une droite horizontale, lorsque l'on est face au patient, passant par ce centre de gravité intercepte alors en premier les cavités droites. - Procédé selon l'une des revendications 1 à 9, caractérisé en ce que les cavités cardiaques sélectionnées sont les cavités gauches. 11 - Procédé selon l'une des revendications 2 à 10, caractérisé en ce que la troisième érosion est une érosion de l'ordre de 10 voxels.