Brevet FR2776798A1 — publié le 19991001

"PRODE DE TRAI TEMENT D'I MAGES I NCLUANT DES ETAPES DE SECGMENTATI ON D'UNE I MAGE MULTI Dl MENSI ONNRI FE ET APPARE] L D'I MAGERI E MEDI CALE UTI U SANT CE PROCEDE" Description L'invention concerne un procédé de traitement d'une image multidimensionnelle, formée de points, représentant des objets d'intensités différentes d'un objet à l'autre, limités par leurs bords vis-à-vis d'un fond, ce procédé indcluant des étapes de segmentation d'image. L'invention concerne également un appareil d'imagerie médicale mettant en ceuvre le procédé de segmentation, et en particulier un appareil à rayons X L'invention trouve son application dans les systèmes d'imagerie médicale pour la visualisation en relief du système vasculaire. Par image multidimensionnelle, il faut entendre par exemple une image d'intensité formée de points qui sont des pixels disposés selon une matrice bidimensionnelle dans l'espace avec l'intensité des pixels en troisième dimension ou une image d'intensité formée de points qui sont des voxels disposés selon une matrice tridimensionnelle dans l'espace avec l'intensité des voxels en quatrième dimension. Dans un procédé de traitement d'image, la segmentation détermine des segments qui sont des régions de l'image d'intensités différentes. Les limites des segments sont appelées interfaces. Un procédé de traitement d'images pour réaliser une segmentation d'image est déjà connu de la publication intitulée" A Multiscale Algorithm for Image Segmentation by Variational Method" par G. KOEPFRIER, C. LOPEZ, J. M. MOREL dans" SIAM J. NUMER ANAL, vol.31, No.1, pp.282- 299, February 1994 ". Le procédé décrit dans cette publication comprend un traitement d'une image bidimensionnelle qui comprend au départ des régions formées de pixels délimitées par des interfaces et dont un certain nombre a à peu près la même intensité. A cet effet, ce procédé comprend une étape qui met en oeuvre un algorithme de segmentation. Le but de l'algorithme proposé est de supprimer le plus grand nombre possible d'interfaces afin de fusionner des régions adjacentes dont les intensités sont pratiquement identiques. Cet algorithme est une fonction de coût appelée Energie. La fusion de deux régions adjacentes est possible seulement dans le cas o elle résulte en une minimisation de la fonction Energie. Cette fonction Energie comprend deux termes: un premier terme prenant en compte la variance d'intensité dans chaque région de l'image, et un second terme prenant en compte la longueur totale des frontières dans l'image pondérée par un paramètre, appelé paramètre d'échelle. La mise en ceuvre de l'algorithme consiste à donner d'abord au paramètre d'échelle la valeur 1 et à fusionner deux régions adjacentes, s'il s'en trouve, qui minimisent la fonction Energie. Les régions résultantes sont alors 2 2776798 réorganisées en éliminant l'interface des deux régions fusionnées, les termes de la fonction Energie sont à nouveau calculés et un nouvel essai de fusion utilisant le paramètre d'échelle 2 = 1 est effectué. Cette opération est recommencée jusqu'à ce qu'il ne se trouve plus aucune région ayant une région adjacente pour fusionner quand le paramètre d'échelle = 1. Les régions résultantes sont réorganisées après chaque fusion en éliminant les interfaces. On effectue ensuite les mêmes opérations avec le paramètre d'échelle / = 2, etc., jusqu'à ce qu'on ne puisse plus minimiser la fonction Energie. Pour finir, les interfaces des régions adjacentes ayant des intensités différentes sont concrétisées par des bords et sont reportées sur l'image d'origine. Dans cet algorithme, lorsque le terme lié à la longueur des interfaces a un coût important, cela incite à enlever l'interface des deux régions adjacentes qui produit ce coût important, à condition que la variance d'intensité entre ces régions soit faible, et donc cela incite à fusionner ces deux régions adjacentes afin de minimiser la fonction Energie. Du fait que l'on cherche à fusionner les régions en commençant par utiliser les plus faibles valeurs du paramètre d'échelle,2, cela signifie que l'on réalise d'abord une très faible pénalisation des interfaces, et donc que l'on permet à de nombreuses régions adjacentes d'intensité très voisines de conserver très longtemps leur interface. Puis, au fur et à mesure que la valeur du paramètre d'échelle, augmente, les régions sont fusionnées deux à deux de manière lente et graduelle jusqu'à ce que la fonction Energie ne puisse plus être minimisée. Cet algorithme est valable à plusieurs échelles de l'image, d'o la dénomination de, comme paramètre d'échelle, et celle de l'algorithme comme pyramidal. Un problème important est que le procédé de segmentation connu du document cité est très coûteux en temps de calcul du fait que les régions sont fusionnées de manière très graduelle en utilisant d'abord la plus basse valeur possible du paramètre,2, et en augmentant la valeur de ce paramètre lentement pas à pas, alors qu'une réorganisation des régions résultantes, avec un nouveau calcul des termes de la fonction Energie, doit être effectuée après chaque fusion pour pouvoir essayer à nouveau de minimiser la fonction Energie et, occasionnellement, d'effectuer une autre fusion. Un autre problème encore plus important est que l'on n'est pas assuré de pouvoir empêcher une fusion entre deux régions adjacentes d'intensités substantiellement égales appartenant néanmoins à des objets différents. Un but de l'invention est de fournir un procédé de segmentation d'image multidimensionnelle qui soit beaucoup moins gourmand en temps de calcul que le procédé connu du document cité. Un autre but de l'invention est de fournir un tel procédé de segmentation qui assure de pouvoir d'une part distinguer d'un fond et d'autre part distinguer 3 2776798 entre eux des objets différents quel que soit leur niveau d'intensité vis-à-vis les uns des autres et vis-à-vis du fond. Ces buts sont atteints et les problèmes cités sont résolus au moyen d'un procédé de traitement d'une image multidimensionnelle formée de points représentant des objets d'intensités différentes d'un objet à l'autre, limités par leurs bords vis-à-vis d'un fond, ce procédé incluant des étapes de segmentation de l'image en cellules appartenant au fond et en cellules appartenant aux objets, mettant en ceuvre des critères qui imposent que lesdites cellules ont respectivement des intensités substantiellement uniformes, que ces cellules n'intersectent ni n'englobent aucun bord d'objet et que les bords des objets coïncident avec10 des bords de cellules. Dans une application, ce procédé de traitement d'image comprend, préalablement à l'étape de segmentation en cellules objets et cellules fond, des étapes d'acquisition d'images d'intensité bidimensionnelle (2D) d'un volume, et de reconstruction d'une image tridimensionnelle bruitée (3-D), et comprend, ultérieurement à l'étape de segmentation de l'image tridimensionnelle en cellules objets et cellules fond, des étapes de reconstruction d'une image tridimensionnelle par une méthode mettant en ceuvre des rendus de volumes o chaque objet se voit attribuer respectivement l'intensité individuelle de la cellule correspondante. Dans une autre application, un appareil d'imagerie médicale comprend un système d'acquisition d'une séquence d'images, et un système de traitement d'image qui a accès aux données des images et à un système d'affichage des images et qui comprend un processeur mettant en ceuvre un tel procédé de traitement d'image incluant ces étapes de segmentation d'image. Dans un mode d'application, de ce procédé, l'image multidimensionnelle est une image médicale vasculaire. Dans une application particulière, un appareil à rayons X est associé à un tel appareil d'imagerie médicale. L'invention est décrite ci-après en détail, en référence aux figures schématiques annexées, dont la liste est la suivante: les FIGs.1A et 1 B qui représentent sous forme de diagrammes de blocs fonctionnels les étapes d'un procédé de reconstruction d'une image multidimensionnelle segmentée; la FR G2 qui représente sous forme d'un diagramme de blocs fonctionnels les sous-étapes de pré-segmentation du procédé de reconstruction d'image; la RF G3 qui représente sous forme d'un diagramme de blocs fonctionnels les sous- étapes de segmentation proprement dite de ce procédé; 4 2776798 les RGs 4A, 4B et 4C qui illustrent respectivement la formation des bassins-versants, et l'élimination des lignes de partage des eaux durant l'étape de pré- segmentation; les R Gs.5A et 5B qui représentent respectivement des régions avant et après fusion dans une image bidimensionnelle durant l'étape de fusion; les FI Gs.6A et 6B qui montrent respectivement un histogramme simple et un histogramme cumulé du paramètre de pondération de la fonction de coût; les R Gs.7A, 7B, 7C qui reproduisent respectivement une image médicale 2-D appelée projection représentant une partie du système vasculaire du cerveau, une image 3-D reconstruite à partir de plusieurs projections 2-D obtenue sans traitement de segmentation et une image 3-D reconstruite après segmentation; la RF GL8 qui montre un appareil d'imagerie à rayons X L'invention concerne un procédé de traitement d'image pour la reconstruction d'une image tridimensionnelle, notée 3-D ayant une quatrième dimension d'intensité, à partir d'une image 3-D bruitée, au moyen d'étapes incluant un dé-bruitage de cette image 3-D, une pré- segmentation et une segmentation pour séparer d'un fond rendu uniforme des objets ayant des intensités différentes, mais qui peuvent être substantiellement proches de celle du fond, et pour séparer ces objets entre eux. Les étapes de ce procédé de traitement d'images sont illustrées sous forme de diagrammes de blocs fonctionnels aux RFGs.1 A, 1 B, 2A, 2B, et 3. L'invention peut être généralisée à des images ayant une ou plusieurs dimensions spatiales avec une dimension d'intensité, appelées images multidimensionnelles. Dans un exemple illustré par la FIRG.7B, I'image d'intensité à traiter est une image médicale vasculaire tridimensionnelle donc constituée de voxels ayant trois coordonnées spatiales, avec une quatrième dimension formée par l'intensité des voxels. Dans l'application envisagée, l'image peut contenir 1283 voxels. En référence à la FIG.1A, cette image tridimensionnelle notée 3-D provient d'une reconstruction 10 d'images vasculaires induant des sous- étapes de: acquisition 11 d'un certain nombre d'images bidimensionnelles de rayons X appelées projections 2-D; ces projections 2-D, dont l'une est illustrée par la FIRG.7A, sont les sorties d'un système d'acquisition d'images de rayons X qui tourne sur 180 dans un plan de rotation autour d'un volume constitué par une partie du corps d'un patient, telle que le cerveau, et qui sont des images 2-D du système vasculaire, soustraction 12 du fond appliquée aux projections 2-D comprenant une première acquisition de ces projections à des angles donnés du plan de rotation lorsqu'un produit de contraste a été injecté dans le système vasculaire du patient, une seconde acquisition aux mêmes angles sans produit de contraste et la soustraction des 2776798 intensités des projections correspondantes à l'aide de marqueurs, pour extraire le fond, résultant en des projections 2-D DSA (de l'anglais Digital Substracted Angiography), reconstruction 13 d'image 3-D du système vasculaire du cerveau du patient en combinant les projections 2-D DSA de manière appropriée par des méthodes de rétro- projection connues de l'homme du métier. En référence à la RG.7A, on constate que les projections 2-D permettent difficilement au praticien de distinguer les différents vaisseaux du système vasculaire. En référence à la RF G7B, I'image 3-D, qui est affichée en utilisant une technique connue qui réalise un rendu de surfaces volumiques mieux approprié à réaliser un bon diagnostic. En particulier, en observant comparativement les RGs.7A et 7B, on constate que le vaisseau en forme de S le plus à gauche dans l'image 3-D de la FI G7B est la reconstruction visible du vaisseau pratiquement indistinct situé dans la zone correspondante de l'image 2-D de la RG7A En outre, une telle image volumique permet au praticien de "voyager " dans l'espace, c'est-à-dire dans le système vasculaire. Les projections 2-D DSA, illustrées par la RG7A, sont bruitées; donc l'image 3-D reconstruite, illustrée par la RG7B, est également bruitée. Ceci provient en particulier du fait que le système d'acquisition n'est pas parfait, que durant l'acquisition la distribution du produit contrastant varie, que la mise en coïncidence de deux projections, pour l'extraction du fond, n'est pas parfaite, et que l'algorithme de reconstruction 3-D lui- même n'est pas parfait. Les imperfections de l'image de la FRGl7B gênent le praticien qui veut évaluer la présence ou l'ampleur d'anomalies. En référence à la R G 1 A, I'image 3-D reconstruite doit alors subir de préférence deux phases d'optimisation qui sont une première phase 14 de réduction du bruit, une seconde phase 100 de segmentation. La première phase 14 de réduction du bruit met en ceuvre des méthodes connues de l'homme du métier qui ne sont pas décrites ici. La seconde phase 100 met en oeuvre un procédé de segmentation pour différentier puis séparer des premières régions comprenant exclusivement les informations essentielles, constituées par les vaisseaux dans l'application envisagée, et des secondes régions comprenant tout le reste des éléments de l'image, c'est- à-dire tout ce qui n'est pas des vaisseaux et qui est considéré comme le fond. Les premières et secondes régions de l'image séparées sont appelées sur un fond uniforme. Un problème particulier apparaît durant cette seconde phase 100 de segmentation du fait que, dans l'image 3-D à segmenter, peuvent se trouver d'une part des vaisseaux de forte intensité et d'autre part des vaisseaux d'intensité plus faible peu différente de l'intensité du fond. Il n'est donc pas possible de segmenter l'image en 6 2776798 uniquement une classe "vaisseaux" et une classe des vaisseaux dans l'image 3-D reconstruite. La présente invention propose un procédé de segmentation 100 qui permet d'obtenir une reconstruction avantageuse de l'image 3-D o: tous les vaisseaux sont séparés du fond et sont conservés avec leur propre intensité dans la même image 3-D, quels que soient leurs différents niveaux d'intensité vis-à-vis du fond, dans l'image à segmenter, ce qui est très important pour établir un diagnostic médical d'insuffisance vasculaire, les dimensions, en particulier les diamètres, et les formes des vaisseaux sont reproduites de manière très fidèle, ce qui est très important pour établir un diagnostic médical de présence d'anomalies du diamètre telles que les sténoses et d'anomalies des formes telles que les anévrismes. En outre, le procédé de segmentation 100 proposé selon l'invention nécessite un temps de calcul beaucoup plus court que le procédé déjà connu. Par exemple, le procédé de segmentation avec reconstruction d'image 3-D selon l'invention peut être mis en ceuvre en utilisant des machines SPARC ou ULTRASPARC, connues à ce jour de l'homme du métier, qui réalisent les calculs nécessaires à la reconstruction d'une image 3- D segmentée de 128 3 voxels en environ 1 min. En référence à la FIG 1B, ce procédé de segmentation 100 comprend- deux étapes générales: une opération de pré-segmentation 30 qui consiste à déterminer, dans cette image 3-D, une pluralité de cellules élémentaires, non géométriques, définies par une première condition qui exclut qu'une cellule élémentaire contient ou coupe un vaisseau ou un bord de vaisseau et par une seconde condition qui indclut qu'une cellule élémentaire contient uniquement des voxels d'intensités substantiellement identiques, une opération de segmentation proprement dite 40 dans laquelle des paires de cellules élémentaires adjacentes de l'image pré-segmentée sont combinées par une méthode de fusion (en anglais: merging) pour former des cellules résultantes, en conservant la première condition qu'une cellule résultante ne contient ni ne coupe aucun vaisseau ou bord de vaisseau et la seconde condition qu'une cellule élémentaire ne contient 7 2776798 que des voxels d'intensités substantiellement proches, et l'opération de fusion est poursuivie jusqu'à ce que l'image contienne un nombre de cellules diminué par rapport à l'image pré- segmentée, ce nombre étant choisi nécessaire et suffisant pour séparer les vaisseaux les uns des autres et du fond.5 En référence à la R G.1 B, ces étapes sont suivies d'une troisième étape générale 60 dans laquelle, en référence à la FI G7C, une image 3-D optimisée est reconstruite, par exemple par une méthode de rendu de surfaces volumiques, représentant tous les vaisseaux avec leurs différentes intensités, leurs dimensions et formes fidèlement reproduites, sur un fond uniforme.10 En référence à la RFI G.2, I'étape 30 de pré-segmentation comprend une détermination, dans une image 3-D dé-bruitée issue de l'étape 14, de cellules élémentaires illustrées par la RFl G5A, sans forme prédéterminée, répondant à la première condition précédemment définie qu'aucune cellule élémentaire ne contient ni ne coupe de vaisseaux ou bords de vaisseaux et à la seconde condition d'uniformité d'intensité. Un problème essentiel consiste donc à localiser les positions exactes des parois des vaisseaux. Dans le cas de données parfaites exhaustives et non bruitées, le module du gradient d'intensité, évalué dans une image 2-D ou 3-D, est différent de zéro (0), si ce gradient d'intensité est évalué dans le voisinage d'une paroi d'un vaisseau. Cette propriété peut aider à résoudre le problème de localisation des parois des vaisseaux. Cependant, la solution à ce problème n'est pas simple, car, dans le cas d'une image 3-D, une paroi de vaisseau est définie comme un ensemble de voxels o le module du gradient atteint un maximum local le long de sa propre direction. En référence à la RF G2, une solution appropriée pour résoudre ce problème et pour effectuer la pré-segmentation 30 comprend une méthode particulière de classification des voxels indcluant une transformation de l'image 3-D en image paysage. Lorsque l'on a affaire à une image plane ou 2-D, il est facile de représenter cette image comme un paysage à 3 dimensions spatiales, dont les deux premières dimensions spatiales x, y, sont celles de l'image plane 2-D et dont la troisième dimension spatiale est constituée par l'intensité I (x, y) des pixels formant des reliefs. Par analogie, on peut également construire un tel paysage en partant d'une image volumique 3- D o une fonction de l'intensité devient une quatrième dimension spatiale pour former des reliefs. Dans une image 3-D, ce sont les modules des gradients qui forment des surfaces d'une image paysage que l'on peut se représenter et ce sont les variations des modules des gradients qui forment les reliefs, vallées CB ou crêtes WS de cette image paysage. La R G4A illustre un tel paysage 3-D, en coupe selon l'axe des x. Ce paysage comprend par exemple une montagne M avec une crête WS, et des vallées CB o se recueille l'eau de pluie coulant le long des versants de la montagne M, de part et d'autre de la ligne de crête WS qui est la ligne de partage des eaux. Chaque vallée CB de part et d'autre de la crête WS est un bassin 8 2776798 hydrographique appelé bassin-versant. On détermine des bassins-versants CB, comme des classes indépendantes et on définit une méthode de dclassification des voxels de l'image par laquelle deux voxels appartiennent à une même dclasse, c'est-à-dire à un même bassin versant CB, si des gouttes d'eau fictives qui tombent en ces voxels sur une surface de 5 I'image paysage atteignent le même minimum de ladite surface. Par cette méthode de classification, les voxels qui sont en dehors de tout bassin versant, c'est-à-dire en dehors de toute dclasse, sont sur les lignes de partage des eaux WS, donc sur les frontières des classes définies par les bassinsversants. On décide ensuite que les cellules élémentaires C de l'image pré-segmentée JPS montrée sur la R G.5A10 sont identifiées avec les bassins-versants C8. En référence à la RF G2, en pratique, cette étape 30 de pré- segmentation comprend des sous-étapes de: détermination 31, en chaque voxel, du module du gradient JJJ, d'un vecteur gradient j de ce module de gradient I, i et de son module FI détermination 32 d'une classe des bassins-versants CB effectuée en suivant, en chaque voxel, la direction des vecteurs gradients j des gradients È(l[ jusqu'à un minimum local Min I(G du module du gradient l; I'ensemble des voxels qui convergent vers un même minimum local Min 1 de module du gradient I(1 sont regroupés comme appartenant au même bassin-versant CB, détermination 33 d'une dasse des lignes de partages des eaux WS: les voxels dont les modules du gradient de gradient FIj s'annulent sans correspondre à des minima locaux Min passent entre des couples de voxels notés VBD formés d'un voxel dans chacun des bassins- versants adjacents; comme illustré par les RFl G.4B, 4C, les voxels des lignes de partage des eaux WS sont simplement attribués lors de balayages alternatifs SCN aux bassins-versants situés de part et d'autre l'interface BD, décision 36 que les bassins-versants résultants comportant les voxels 9 2776798 des bassins-versants C8 plus les voxels répartis des lignes de partage des eaux WS sont les cellules élémentaires recherchées, ayant les numéros 1 à n attribués aux bassins-versants. De l'utilisation de cette méthode de classification des voxels et d'identification des cellules élémentaires aux bassins-versants, il résulte que ni les parois des vaisseaux ni les vaisseaux ne peuvent traverser les dclasses des cellules ou bassins-versants CB, puisque les parois des vaisseaux sont aux frontières des dclasses, c'est-à-dire en coinddence avec les crêtes WS et il en résulte que les vaisseaux forment eux-mêmes des classes indépendantes. Il résulte aussi que ces parois de vaisseaux peuvent être détectées avec la précision de détection des lignes de partage des eaux. Cette méthode de classification et d'identification répond donc bien aux deux conditions imposées dans l'étape de pré-segmentation. En référence à la FIG5A, I'image est maintenant formée de cellules élémentaires numérotées C1,...q,...C, ne franchissant et ne contenant aucun vaisseau. Après l'étape de pré-segmentation 30, le nombre 1 à n des cellules élémentaires est trop grand pour reconstruire une image 3-D par rendu de surfaces volumiques avec profit. On applique maintenant, en référence aux FlGs lB, 3A, 3B, une étape de segmentation 40 proprement dite. Cette segmentation 40 est obtenue par l'application sur l'image pré-segmentée d'un algorithme de fusion pour minimiser le nombre d'interfaces en réunissant ensemble des cellules adjacentes ayant des propriétés de similarité. Cette similarité est mesurée par une fonction de coût, notée F, qui comprend une somme de deux termes. Dans une forme générale, la fonction de coût s'écrit F = fi[cr2(C1)]+i;f49G) (1) et en particulier: = 2(Ci) + 1/ + 2 (2a) Une fusion peut être effectuée entre deux cellules adjacentes, si et seulement si cette fusion résulte en ce que la fonction de coût F est minimisée. Les termes de la fonction de coût constituent des critères pour réduire le nombre des cellules par fusiori sous les conditions que l'intensité des voxels à l'intérieur des cellules fusionnées est substantiellement uniforme, et qu'aucune cellule résultant de la fusion ne contient ou traverse des vaisseaux. Le premier terme est un critère qui prend en compte les différences d'intensité des cellules résultant de la pré-segmentation, et qui est constitué par une fonction de sommation d'une fonction f [mo2 (Ci)] des variances notées 2(C;) des intensités dans les différentes cellules C, o les variances sont multipliées par le nombre m1 des voxels dans la cellule C, considérée. La variance o2(C) est calculée dans chacune des cellules en prenant la valeur moyenne du carré de l'écart type d'intensité par rapport à la moyenne d'intensité dans la cellule. Pour évaluer facilement le premier critère, l'intensité interne aux cellules est 2776798 rendue substantiellement uniforme en calculant individuellement la moyenne d'intensité iMOY dans chaque cellule de l'image pré-segmentée et en attribuant cette moyenne d'intensité à chacun des voxels de la cellule considérée. La fonction de coût peut alors être écrite de façon équivalente: F =;[(x,y)-ImOY]2 +iXl/[+ lGlI/lIGotI)2i (2b) o la fonction de la variance f (.2) dans le premier terme est la fonction de sommation des carrés des différences entre l'intensité I (x, y) d'un voxel de coordonnées x, y et l'intensité moyenne iMOY de la cellule C, à laquelle ce voxel appartient. On rappelle que, dans tous les cas, lors d'une fusion, la variance d'intensité de la cellule résultant de la fusion est toujours supérieure à la somme des variances d'intensité des cellules avant fusion (loi statistique). Donc ce critère défavorise toute fusion mais il est contrebalancé par le second terme. Le second terme est un critère qui est pondéré par un facteur de pondération ,, et qui est constitué par une fonction de sommation, sur tous les couples de voxels situés à l'interface des cellules, d'une mesure spécifique notée f2,GI;) prenant en compte le module du gradient IIGII en chacun des voxels de ces couples à l'interface, o IIGO)i est un facteur de normalisation. On effectue, dans le second critère, sur les couples de voxels VBD, situés de part et d'autre d'une interface BD, dans les cellules adjacentes la somme; f2(IGII) prenant en compte des modules de gradients 1G[l déjà calculés. Pour une valeur donnée fixée du facteur de pondération,, ce second critère constitue une pénalité qui est d'autant plus forte que les modules desgradients des voxels d'interface entre deux cellules sont faibles. Ce second critère incite donc à fusionner les cellules dont les interfaces montrent des modules du gradient faibles, car, par cette fusion, le terme de pénalité est minimisé d'o il résulte que la fonction de coût F tend vers une minimisation, même en tenant compte d'une augmentation de la variance d'intensité due à la fusion. En revanche, ce second critère apporte une pénalité faible ou nulle lorsque les voxels d'interface de deux cellules adjacentes montrent des modules du gradient forts, ce qui n'incite pas à fusionner ces cellules, compte tenu du fait que la variance d'intensité qui augmente lors de la fusion empêche de minimiser la fonction de coût F. Ceci permet donc de préserver la séparation des classes vaisseaux et des classes fond, puisque l'on sait que le module du gradient a une valeur forte sur les bords des vaisseaux. Ce type de fonction de coût pour effectuer des fusions n'était pas connu de la publication citée au titre d'état de la technique. Selon l'invention la première étape de pré-segmentation détermine, grâce aux lignes de partages des eaux WS, des interfaces BD de localisations très précises. La seconde étape de segmentation sélectionne, 11 2776798 grâce à la fonction de coût F spécifique, d'une part les parties des interfaces BD qui sont réellement des bords de vaisseaux ayant une localisation précise, et d'autre part les différentes parties de vaisseaux qui constituent chacune une dclasse à part entière bien délimitée des classes formant le fond par les interfaces BD qui coïncident avec les bords des vaisseaux. Cest pour cette raison qu'aucune dclasse ne peut couper un vaisseau et qu'aucun bord de vaisseau ne peut être éliminé lors d'une fusion. Dans l'étape de segmentation 40 proprement dite, on opère des fusions de cellules adjacentes deux à deux lorsque les voxels de ces cellules présentent des caractéristiques prises en compte par les deux critères, et qui sont telles que la fonction de coût F est minimisée par la fusion. Si la fonction de coût F n'est pas minimisée mais au contraire croît par la fusion, cette fusion n'est pas opérée. Il faut noter que le second critère de la fonction de coût F est en outre pondéré par un facteur de pondération variable . Pour mettre en oeuvre l'opération de segmentation 40 au moyen de la fonction de coût F, on fait croître le paramètre, graduellement depuis une valeur de départ, notée 2*0 S par exemple, la valeur,o est prise égal à zéro, alors il ne reste que le premier terme dans la fonction F, et il est tout à fait dclair qu'il n'y a aucun intérêt à fusionner deux cellules, car la variance d'une cellule fusionnée est toujours plus grande que la somme des variances de deux cellules individuelles (loi statistique). Pour /t # 0, les cellules ayant été préalablement numérotées, on visite toutes les cellules C1 à C qui sont adjacentes et qui forment donc des paires. Pour chaque paire, la fonction F opère une compétition entre le premier terme relatif à la variance et le second terme relatif à l'inverse du module du gradient des voxels d'interface, pondéré par le facteur,. Comme on l'a vu précédemment le premier terme relatif à la variance augmente en cas de fusion, alors que le second terme diminue du fait que l'interface des cellules fusionnées est supprimée, ces cellules n'en faisant plus qu'une, et donc que le nombre de voxels sur lequel le second terme doit être calculé diminue. Pour une paire de cellules adjacentes données, la compétition entre le premier et le second termes de la fonction de coût F peut résulter en une minimisation de la fonction de coût F et donc en une fusion des cellules de cette paire, en prenant en compte non seulement la fonction de variance d'intensité et la fonction des modules des gradients mais aussi la valeur variable du paramètre 2 de pondération du second terme. FPour un nombre donné de cellules issues de la pré-segmentation, le premier terme de la fonction de coût F est fixé, et donc plus la valeur de est faible, moins il y a de possibilité de minimiser la fonction de coût F et de fusionner des cellules, alors que plus la valeur de, est grande, plus il y a de possibilité de minimiser la fonction F et de fusionner des cellules. 12 2776798 Après chaque fusion, I'intensité de la cellule résultante est calculée comme la moyenne des intensités des voxels des deux cellules fusionnées et les couples de voxels d'interface des deux cellules fusionnées deviennent des voxels internes ordinaires. Il faut maintenant rappeler qu'un des buts de l'invention est de réaliser les opérations de segmentation en un temps plus court que ce qui est actuellement connu de l'état de la technique. Cest pourquoi, I'étape 40 de segmentation proprement dite qui utilise la fonction de coût F est mise en ceuvre en utilisant des sous-étapes spécifiques décrites plus loin qui indcluent, pour alléger le temps de calcul, la construction d'un premier histogramme simple à partir de l'évaluation des seuils de fusion de couples de deux cellules adjacentes pour fournir le nombre de couples de cellules qui fusionnent pour un même seuil de fusion donné; et la construction d'un second histogramme appelé histogramme cumulé, à partir du premier histogramme simple, pour fournir le nombre de couples de cellules qui fusionnent pour un seuil de fusion inférieur ou égal à un seuil particulier. Purement à titre d'exemple, pour illustrer cette description, on suppose qu'à l'issue de la pré-segmentation 30, on a obtenu 20.00015 dclasses avec 50.000 interfaces dans l'image 3-D pré- segmentée. Théoriquement, la segmentation 40 devrait comporter les sous-étapes suivantes: évaluation de la fonction de coût F avant toute fusion, alors notée Fo: dans le premier terme de Fo, la fonction de sommation est calculée sur les 20 000 cellules issues de la pré-segmentation, et dans le second terme, la fonction de sommation est calculée sur les 50 000 interfaces, sélection pour examen d'un couple de deux cellules C, q, évaluation de la fonction de coût F après la fusion des deux cellules de ce couple, la fonction de coût étant alors notée F1,: dans le premier terme de F1, la somme est effectuée sur les 20 000 cellules moins une, et dans le second terme, la somme est effectuée sur les 50 000 interfaces moins une, détermination d'un seuil de fusion À (j, k) à partir du paramètre de' pondération A qui constitue le seuil auquel l'interface entre les deux cellules du couple examiné disparaît et qui permet la fusion des deux cellules du couple examiné, cette fusion résultant en la fonction de coût F: on définit ' (j, k) comme la valeur qui permet d'annuler une fonction Eo calculée comme la différence entre les fonctions de coût F1 et Fo, dans lesquelles le paramètre de pondération est 2' (j, k) et on écrit: E = F - Fo= 0 avec,2.0 ( j, k) (3) Il en résulte que 2o (j, k) est aisément calculé pour deux cellules particulières C, Ck avec: tO. ..s rx. . (4) 13 2776798 o IMoY' et IMoY sont les valeurs moyennes d'intensité et m mk les nombres de voxels j k j k dans les cellules. Dans la pratique, l'étape 40 de segmentation proprement dite ne comprend pas l'évaluation d'une fonction de coût F telle que celle des formules (1) ou (2) car ce calcul est trop long. L'étape 40 de segmentation est fondée sur l'évaluation du seuil de fusion)]' (j, k) et comprend les sous-étapes suivantes: détermination 41 des couples de cellules adjacentes q, C dans l'image pré-segmentée; il y a autant d'interfaces BDjk que de couples, détermination 42 des intensités moyennes IMOY, 1mOY et des j k nombres m j, mk de voxels des cellules adjacentes q, k des couples, détermination 43 d'autant de seuils de fusion 2o (j, k) qu'il y a de couples de cellules adjacentes C, k, en utilisant la formule (4), dclassification 44 des seuils 20 (j, k), du plus petit au plus grand, que l'on note,il à oMAX, o o construction d'un premier diagramme 51, noté h/ illustré par la FIRG6A, appelé histogramme simple, basé sur le fait que, en pratique, un certain nombre de couples de cellules montrent, dans l'évaluation de la fonction de coût F, le même premier terme et le même. second terme, d'o il résulte que les cellules de ces couples peuvent fusionner pour le même seuil: On porte en abscisse les valeurs des seuils l à,oMAx qui sont des nombres pouvant théoriquement être compris entre 1 et 50.000, s'il y a 50.000 interfaces, et on porte en ordonnée le nombre d'occurrence hp c'est-à-dire le nombre hp de fois o l'on trouve des interfaces qui fournissent E_ = 0 avec A de /4 à A,0x o 1 14 2776798 histogramme simple 51: ainsi les valeurs]0 à 255 sont regroupées dans une tranche notée et, les valeurs,2256 à 2sl sont regroupées dans une tranche notée e2, et ainsi de suite, jusqu'à une dernière tranche notée enax; dans chaque tranche les nombres d'occurrences hp sont additionnées et notées par analogie hl à ho ax; ce sont alors en fait les valeurs p 0 2776798 fusionner un nouveau pourcentage de 10 ou 15 % par exemple de couples de cellules; on opère ainsi jusqu'à ce qu'on ait obtenu un nombre approprié de cellules; on peut choisir de garder par exemple 10 ou 15% des cellules initiales: dans les cellules restantes les voxels d'interfaces éliminés sont des voxels internes ordinaires et les intensités sont les intensités5 moyennes des cellules adjacentes fusionnées. En référence à la RG. 1B et à la FIG7Q l'image 3-D segmentee peut alors subir une opération 60 de reconstruction avec rendu des surfaces volumiques pour être affichée ou enregistrée. Les vaisseaux sont reproduits, grâce à la segmentation selon l'invention, chacun avec un niveau d'intensité spécifique qui permet de les distinguer les uns10 des autres et du fond. En référence à la R G.8, un appareil de radiographie numérisée comprend une source 1 de rayons X, une table 2 pour recevoir le patient, un dispositif 3 intensificateur d'image couplé à un tube video 4, qui fournit des données à un système 5 de traitement numérique d'image comprenant un processeur pour mettre en ceuvre le procédé de15 traitement d'image selon l'invention. Ce processeur comprend plusieurs sorties dont une sortie 6 couplée à un moniteur 7 pour visualiser des images radiographiques, ou à un système enregistreur approprié 8. Le système de traitement d'image peut faire partie intégrante de l'appareil à rayons X ou bien en être éloigné tout en étant couplé. D'une manière générale, le procédé de traitement d'image précédemment décrit peut être exploité partiellement ou entièrement par tout système d'imagerie médicale, et en particulier appliqué à des appareils à rayons X tels que décrits plus haut pour le traitement d'image multidimensionnelle, avec ou sans reconstruction des surfaces volumiques. 16 2776798 REVENDICATI ONS: 1. Procédé de traitement d'une image multidimensionnelle, formée de points, représentant des objets d'intensités différentes d'un objet à l'autre, limités par leurs bords vis-à-vis d'un fond, ce procédé incluant des étapes de segmentation (100; 30,40) de I'image en cellules appartenant au fond et en cellules appartenant aux objets, mettant en ceuvre des critères qui imposent que lesdites cellules ont respectivement des intensités substantiellement uniformes, que ces cellules n'intersectent ni n'englobent aucun bord d'objet et que les bords des objets coïncident avec des bords de cellules. 2. Procédé selon la revendication 1, comprenant une étape de pré- segmentation (30) de l'image en un premier nombre de cellules en effectuant une dclassification des points de l'image dans les différentes cellules selon un critère qui est une fonction des gradients d'intensité (IIGII) de ces points, et une étape de segmentation proprement dite pour réduire le premier nombre des cellules en classant les cellules par couples de cellules adjacentes (q, Ck), en évaluant un seuil de fusion,; (j, k) des cellules des couples comme une fonction de la variance d'intensité ( caractéristiques d'intensité qui conduisent à un seuil de fusion inférieur au seuil particulier. 3. Procédé selon la revendication 2, o dans l'étape de pré-segmentation (30), la détermination des cellules est effectuée en transformant l'image multidimensionnelle en une image paysage par des fonctions des gradients d'intensité des points de l'image, en identifiant (C0) les cellules recherchées à des bassins-versants (C8) et en identifiant les bords délimitant lesdites cellules à des lignes de partage des eaux (W)de l'image paysage. 4. Procédé selon la revendication 3, o les points des lignes de partages des eaux (W) sont versés dans les bassins versants (CB) adjacents situés de part et d'autre de ces lignes, o les cellules reçoivent une intensité uniforme calculée comme la moyenne ( ImoY) d'intensité des points à l'intérieur des bassins versants et forment une classe objets, et o des interfaces (BD) sont définies entre des couples de points de bords (VBD) des couples de cellules adjacentes et coïncident avec des bords d'objets. 5. Procédé selon la revendication 2, ou selon l'une des revendications 3 ou 4, o, dans l'étape de segmentation proprement dite, on définit une fonction de coût (F) comprenant un premier terme qui est une fonction de la variance d'intensité ( C.2) à l'intérieur des cellules présentes dans I'image, et un second terme qui est pondéré par un facteur de pondération variable (,2) et qui 17 2776798 est une fonction des gradients d'intensité aux points de bord des cellules, et on évalue le seuil de fusion (,o (j, k)) pour effectuer la fusion de deux cellules adjacentes (q, C) d'un couple comme la valeur du facteur de pondération variable (,) qui annule la différence (Et) entre la fonction de coût (Fo) calculée avant la fusion des deux cellules et la fonction de coût (F1) calculée après cette fusion. 6. Procédé selon la revendication 5, o la fonction de coût est donnée par la formule: F =[I(x,y)-I2 Y]2i; IIGo)2 dans laquelle le premier terme est la fonction de variance d'intensité exprimée par une première fonction de sommation, sur toutes les cellules (C,,), des carrés des différences entre I'intensité I (x, y) d'un voxel de coordonnées (x, y) et l'intensité moyenne (I MOY) dans les cellules (C) à laquelle ce voxel appartient, et dans laquelle le second terme est pondéré par le facteur de pondération (,), et est constitué par une seconde fonction de sommation, sur tous les couples de voxels situés à l'interface des cellules, d'une mesure spécifique prenant en compte le module du gradient normalisé ( IIGII IIGO) des voxels de ces couples d'interface, et o le seuil de fusion ( 20 (j, k)) qui permet de réaliser la fusion de deux cellules (Cl, C) adjacentes entre elles, est évalué comme une fonction des différences des intensité moyennes (IMOY, IM Ot) des deux cellules à fusionner et du nombre de points j k ( m, mk) dans ces cellules, et de la seconde fonction de sommation, selon une formule ko (j,k)- m)mt (iMoY _ Y t, (m; +mi)vl/[l+ 101i)/110ol|] 7. Procédé selon l'une des revendications 5 ou 6, o on évalue le seuil de fusion particulier (o P) qui permet de réaliser concomitamment plusieurs fusions de chacune O deux cellules adjacentes (q, Ck) entre elles, par une suite d'étapes comprenant une détermination des seuils des fusion ( 2' (j, k)) des cellules de tous les couples de l'image, une classification de ces valeurs (,0 (j, k)) du seuil de fusion par valeurs croissantes (,p), ces valeurs étant en nombre au plus égal au nombre total des interfaces présentes dans l'image pré-segmentée (1 18 2776798 d'occurrences de fusions de cellules adjacentes respectivement dans chaque tranche ( P) des o valeurs de seuil de fusion ( o (j, k)), une construction d'un histogramme cumulé (52) en inscrivant dans chaque tranche (!p) la somme cumulée (HP) du nombre desdites occurrences (h P) depuis la o o o première tranche (Cl) jusqu'à une tranche considérée (ô P), et ce dans toutes les tranches o o (tl - tAux) des valeurs de seuil de fusion ( 2), o o O et une définition du dit seuil de fusion particulier comme la tranche (t) o qui comprend le nombre d'occurrences (H P) correspondant au nombre prédéterminé de o fusions concomitantes recherché. 8. Procédé selon la revendication 7, o, à la suite de la détermination du seuil particulier, on sélectionne un nombre prédéterminé d'interfaces à éliminer de façon concomitante entre des cellules adjacentes de couples, on effectue une visite des couples de cellules adjacentes de l'image, et on procède à une fusion des cellules adjacentes des couples chaque fois que le seuil de fusion relatif à ces cellules ( o (j, k)) est inférieur au dit seuil de fusion particulier (op). o 9. Procédé selon l'une des revendications précédentes comprenant préalablement à l'étape de segmentation des étapes d'acquisition d'images d'intensité bidimensionnelle (2-D) d'un volume, et de reconstruction d'une image tridimensionnelle bruitée (3-D), et comprenant, ultérieurement à l'étape de segmentation de l'image tridimensionnelle en cellules objets et cellules fond, des étapes de reconstruction d'une image tridimensionnelle par une méthode mettant en oeuvre des rendus de volumes o chaque objet se voit attribuer respectivement l'intensité individuelle de la cellule correspondante. 10. Procédé selon l'une des revendications précédentes o l'image multidimensionnelle est une image médicale vasculaire. 11. Appareil d'imagerie médicale comprenant un système d'acquisition d'une d'image, et un système de traitement d'image qui a accès aux données de l'image et à un système d'affichage de l'image et qui comprend un processeur mettant en oeuvre un procédé de traitement d'image selon l'une des revendications 1 à 10. 12. Appareil à rayons X associé à un appareil d'imagerie médicale selon la revendication 11.