Brevet FR2799028A1 — publié le 20010330

Procédé de reconstruction d'une image tridimensionnelle d'éléments de fort contraste. La présente invention concerne un procédé de reconstruction d'une image tridimensionnelle d'éléments de fort contraste, en particulier une prothèse présente dans une artère d'un patient, à partir d'un ensemble d'images bidimensionnelles du patient obtenues pour différentes positions d'un appareil de prises de vues. Dans le domaine médical, on procède généralement à la reconstruction des structures internes de patient sous examen, en particulier la reconstruction d'images angiographiques, c'està-dire l'obtention d'images d'arbres vasculaires opacifiés par injection d'un produit de contraste. Les images projetées bidimensionnelles d'un patient sont obtenues par la rotation de l'appareil de prises de vues à rayons X autour du patient. D'une manière classique en angiographie, chaque image acquise est une image soustraite qui est par exemple obtenue par une technique classique de soustraction logarithmique de deux radiographies prises sous la même incidence avant et après une injection d'un produit de contraste dans l'arbre vasculaire dont on souhaite reconstruire l'image tridimensionnelle. Ensuite, on exécute un algorithme de reconstruction d'image tridimensionnelle à partir des images projetées bidimensionnelles de l'objet pour reconstruire le volume tridimensionnel. Il existe plusieurs types d'algorithmes de reconstruction en imagerie par rayons X. On citera notamment l'article de GORDON, BENDER, et HERMAN intitulé "Algebraic reconstruction technic for three dimensional électron microcopy and X ray photography", Journal THEO.BIOL. 29, pages 471 à 781 (1970); l'article "Practical cone-beam algorithm" L.A. Feldkamp, L.C. Davis et J.W. Kress, Journal Optical Society of America A Vol 1 No 6, juin 1984; ou bien l'ouvrage de Anil K. JAIN intitulé "Fondamentals of digital image processing" Prentice Hall information and system sciences series, Thomas Kailath Series Edition; ou encore les demandes de brevets français n 89 03606, n 89 16906 ou n 9610774, auxquels le lecteur est invité à se reporter. Or, une telle reconstruction d'images angiographiques permet de visualiser uniquement les artères dans lesquelles un produit opaque a été introduit, elle ne permet pas la visualisation d'autres éléments intéressants tels que les prothèses. L'invention a pour but de remédier à cet inconvénient en proposant un procédé qui permet la visualisation de prothèse en trois dimensions. Un autre but de l'invention est de permettre la visualisation de prothèse et des artères sur un même écran et selon leurs dispositions réelles. L'invention propose donc un procédé de reconstruction d'une image tridimensionnelle d'éléments de fort contraste à partir d'un ensemble d'images bidimensionnelles d'un objet comprenant lesdits éléments de fort contraste. Pour chaque position différente d'un appareil de prises de vues à rayons X autour de l'objet, une image bidimensionnelle est prise, et la mise en oeuvre d'un algorithme pour la reconstruction de l'image tridimensionnelle est précédée d'une étape de filtrage de l'ensemble d'images bidimensionnelles. L'image tridimensionnelle obtenue ne comporte que les éléments de fort contraste. L'appareil de prises de vues peut ne fournir qu'une seule série d'images bidimensionnelles, dans ce cas aucun produit opaque n'est utilisé, mais il peut également fournir deux séries d'images bidimensionnelles obtenues par deux séquences d'acquisition intercallées d'une injection de produit opaque dans les artères. Les deux séries d'images peuvent servir pour une éventuelle reconstruction angiographique et la première série peut servir pour la reconstruction d'une image tridimensionnelle des éléments de fort contraste selon l'invention. Selon une caractéristique avantageuse de l'invention, le filtrage de chaque image bidimensionnelle acquise est du type passe- bas de façon à obtenir une image bidimensionnelle filtrée ne contenant plus les éléments de fort contraste. On réalise à la suite du filtrage une soustraction logarithmique entre cette image ne contenant plus les éléments de fort contraste et l'image bidimensionnelle acquise de façon à ne conserver principalement que lesdits éléments de fort contraste. Plus précisément, selon l'invention, le filtrage passe-bas sur une image bidimensionnelle acquise consiste à: - prendre l'image bidimensionnelle acquise, - déterminer sa moyenne, - dilater l'image bidimensionnelle acquise de façon à éliminer des objets non désirés de fort contraste et de faible taille par rapport auxdits élément de fort contraste, et - prendre le maximum entre l'image bidimensionnelle acquise dilatée et la moyenne de l'image bidimensionnelle acquise de façon à obtenir une image bidimensionnelle filtrée exempte de tout élément et objet de fort contraste. Selon un mode de mise en oeuvre avantageux de l'invention, lesdits éléments sont des prothèses tels que des spires placées dans une artère d'un patient ou encore des clips vasculaires. De préférence, à la suite de la reconstruction d'image tridimensionnelle contenant les prothèses, on réalise une segmentation par hystérésis de cette image tridimensionnelle de façon à éliminer des voxels inutiles et visualiser uniquement lesdites prothèses dans une image tridimensionnelle primaire. Selon une variante avantageuse de l'invention, pour chaque prothèse: - on réalise une rétro-projection de l'image tridimensionnelle primaire sur les images bidimensionnelles acquises de façon à déterminer les pixels représentant la prothèse sur ces images, - on effectue une approximation en remplaçant l'intensité de chacun de ces pixels par une intensité déterminée à partir des intensités des pixels voisins afin d'obtenir des images bidimensionnelles primaires exemptes de la prothèse, - on réalise une soustraction logarithmique entre les images bidimensionnelles primaires exemptes de la prothèse et les images bidimensionnelles acquises respectives de façon à obtenir de nouvelles images bidimensionnelles comportant uniquement la prothèse, et - on effectue une reconstruction d'image tridimensionnelle suivie d'une segmentation par hystérésis à partir des nouvelles images bidimensionnelles afin d'obtenir une image tridimensionnelle secondaire comportant uniquement la prothèse. Cette variante permet d'obtenir avantageusement une image tridimensionnelle secondaire dans laquelle les spires présentent des niveaux d'intensité ou niveaux de gris plus précis. Selon l'invention, on effectue indépendamment une reconstruction d'image tridimensionnelle angiographique des artères et on réalise une fusion entre l'image tridimensionnelle primaire ou secondaire contenant la prothèse et l'image tridimensionnelle contenant des artères de façon à obtenir une image tridimensionnelle fusionnée montrant la répartition spatiale entre les artères et la prothèse. Selon un mode de réalisation avantageux, au cours de la fusion, lorsque dans l'image tridimensionnelle contenant les artères un voxel comporte une intensité différente de celle que comporte le voxel équivalent dans l'image tridimensionnelle contenant la prothèse, on attribue au voxel résultant de la fusion l'intensité la plus élevée. Avantageusement, avant la fusion de l'image tridimensionnelle primaire ou secondaire contenant la prothèse avec l'image tridimensionnelle contenant les artères, on élimine sur l'image tridimensionnelle contenant les artères des voxels ayant une faible valeur d'atténuation en ne conservant qu'un nombre prédéterminé de voxels de grande valeur d'atténuation. Selon une autre variante avantageuse de l'invention, pour chaque prothèse: - on réalise une rétro-projection de l'image tridimensionnelle primaire sur les images bidimensionnelles acquises de façon à déterminer des pixels représentant la prothèse sur ces images, - on effectue une approximation en remplaçant l'intensité de chacun de ces pixels par une intensité déterminée à partir des intensités des pixels voisins afin d'obtenir des images bidimensionnelles primaires exemptes de la prothèse, - on réalise une soustraction logarithmique entre les images bidimensionnelles primaires exemptes de la prothèse et des images bidimensionnelles contenant des artères et la prothèse acquises après injection de produit opaque aux rayons X de façon à obtenir de nouvelles images bidimensionnelles contenant les artères et la prothèse, et - on effectue une reconstruction d'image tridimensionnelle à partir des nouvelles images bidimensionnelles afin d'obtenir une image tridimensionnelle tertiaire contenant les artères et la prothèse. On évite ainsi d'effectuer l'opération de fusion de deux images tridimensionnelles. D'autres avantages et caractéristiques de l'invention apparaîtront à l'examen de la description détaillée d'un mode de mise en oeuvre nullement limitatif, et des dessins annexés, sur lesquels: - la figure 1 est une image d'une artère comportant un anévrisme dans lequel sont logées des spires, - la figure 2 illustre schématiquement un jeu d'images projetées bidimensionnelles autour d'un objet, - la figure 3 illustre plus en détail l'acquisition de l'une de ces images projetées bidimensionnelles, - les figures 4a, 4b et 4c illustrent une opération de filtrage passe-bas selon l'invention, - la figure 5 est un organigramme d'un mode de mise en oeuvre du procédé selon l'invention, - la figure 6 est une illustration schématique du procédé illustré sur l'organigramme de la figure 5, - la figure 7 est un organigramme d'un autre mode de mise en oeuvre du procédé selon l'invention, et - la figure 8 est un organigramme d'un autre mode avantageux de mise en oeuvre du procédé selon l'invention. La figure 1 illustre une artère 1 comportant un anévrisme 2 dans la tête d'un patient. Le sang passant dans cette artère 1 vient se loger dans l'anévrisme 2 qui forme une poche dans laquelle la pression tend à augmenter. Pour limiter cette augmentation de la pression, on place des spires 3 dans l'anévrisme 2. Les spires 3 sont réalisées à l'aide de matériaux capables de former des caillots de sang, ce qui permet une baisse de la pression dans l'anévrisme 2. On désire visualiser les spires ainsi placées. Si l'on se réfère plus particulièrement aux figures 2 et 3, on voit que le système d'imagerie utilisable pour mettre en oeuvre l'invention permet d'obtenir un ensemble d'images acquises bidimensionnelles Al1-An obtenues en l'espèce par la rotation sur 180 autour de la tête 4 du patient d'une source 5 à rayons X. Chaque image acquise Ai est obtenue à partir d'un détecteur de rayonnement bidimensionnel, par exemple du type amplificateur de luminance utilisé en radiologie, disposé en vis-à-vis du tube à rayons X dans un plan dit plan de projection PPi. Les différents plans de projection sont obtenus par les différentes positions angulaires du détecteur en rotation autour de la tête du patient. Le détecteur est relié à des moyens de traitement 6 comportant notamment des moyens d'échantillonnage reliés à un microprocesseur incorporant de façon logicielle dans sa mémoire de programme associé l'algorithme de reconstruction d'images utilisé dans l'invention et d'une façon générale tous les moyens fonctionnels permettant la mise en oeuvre du procédé selon l'invention. Un étalonnage du système d'imagerie permet notamment de définir un volume virtuel VV entourant l'objet 4 et décomposé en éléments volumiques élémentaires Vi ou "voxels". Il est possible d'utiliser un procédé connu d'étalonnage géométrique automatique d'un système d'imagerie par rayons X, tel que celui décrit dans la demande de brevet français n 93 00804. Les figures 4a, 4b et 4c illustrent schématiquement différentes étapes comprises dans l'opération de filtrage passe-bas selon l'invention. La courbe 7 sur la figure 4a représente par exemple l'intensité de rayons X reçue au niveau du détecteur 6 par une ligne de pixels pour une image bidimensionnelle acquise. La zone 8 correspond à une baisse d'intensité due à la présence des spires 3 dans le faisceau de rayons X. Les spires 3 sont dotées d'un si fort contraste qu'elles absorbent une bonne partie d'intensité de rayons X les traversant. L'image acquise peut comporter quelques objets non désirés de fort contraste, autrement dit du bruit. Un de ces objets non désirés est représenté par un creux 9 de faible dimension. La ligne en pointillé 10 représente la moyenne de l'image en intensité sur cette ligne de pixels. La courbe 11 sur la figure 4b est le résultat d'une opération de dilatation subi par la courbe 7 de la figure 4a. L'opération de dilatation, bien connu de l'homme du métier, consiste à prendre chaque point de la courbe 7 et à remplacer son intensité par un maximum des intensités des points voisins. Le creux 9 disparaît et le creux 8 est légèrement rogné. La courbe 12 sur la figure 4c est le maximum entre la courbe 7 et la moyenne 10. La courbe 12 est le résultat final du filtrage passe- bas. Elle correspond à une image acquise que l'on obtiendrait dans le cas o le patient ne comporte aucun élément de fort contraste. L'image bidimensionnelle obtenue ne contient plus les spires 3. La figure 5 est un organigramme dans lequel des étapes du procédé selon l'invention sont représentées sous forme rectangulaire et les produits de ces étapes, en l'espèce des images, sont représentés sous forme elliptique. Le système d'imagerie permet de réaliser l'acquisition 13 d'un ensemble d'images bidimensionnelles 14 par rotation de la source 5 autour de la tête du patient. Pour chaque position différente de la source 5 autour de la tête du patient, une seule image bidimensionnelle est prise. Les spires 3, dotées d'un fort contraste, sont bien visibles sur les images bidimensionnelles 14. Eventuellement, quelques objets non désirés fortement contrastés sont également présents sur les images bidimensionnelles 14. Chaque image bidimensionnelle 14 acquise subit un filtrage passe-bas 15 selon l'invention. L'image bidimensionnelle filtrée 16 ne contient ni les spires 3 ni les objets non désirés de fort contraste. Ensuite, pour obtenir une image bidimensionnelle ne contenant que des éléments de fort contraste, on réalise une soustraction logarithmique 17 entre l'image bidimensionnelle acquise 14 et l'image bidimensionnelle filtrée 16. L'image soustraite 18 contient les spires 3 et les éventuelles objets non désirés de fort contraste. La reconstruction tridimensionnelle permet d'obtenir une image tridimensionnelle à partir de l'ensemble des images soustraite 18 bidimensionnelles. On utilise une reconstruction de l'image tridimensionnelle à partir des images bidimensionnelles acquises et, par exemple, d'un algorithme itératif de reconstruction algébrique d'image. On effectue une première itération de l'algorithme avec une résolution d'image initiale prédéterminée de façon à obtenir à l'issue de cette première itération des premières valeurs de densité pour les voxels du volume virtuel VV déterminé à l'étalonnage de l'appareil de prises de vues. On subdivise une partie au moins desdits voxels du volume virtuel en plusieurs ensembles correspondant respectivement à des résolutions d'image différentes, multiples ou sous- multiples de la résolution initiale. Et, au sein de chaque itération suivante de l'algorithme, on applique successivement ledit algorithme à partir de chacun desdits ensembles de voxels. L'homme du métier comprendra aisément que d'autres méthodes de reconstruction d'image tridimensionnelle peuvent être utilisées. En d'autres termes, on applique l'algorithme itératif de reconstruction algébrique d'images sur un volume multi-résolutions. Dans ce volume, on sélectionne des voxels qui sont a priori représentatifs des objets d'intérêt à visualiser et on les divise de façon à augmenter la résolution. Les autres voxels, qui sont moins intéressants, car ne concernant pas directement les objets à visualiser, sont soit laissés identiques, soit regroupés au moins une fois de façon à diminuer la valeur de la résolution, mais sont toutefois utilisés pour les calculs de reconstitution d'image, ce qui permet d'obtenir en final des images de très bonne qualité avec un temps de calcul réduit. Au cours de l'étape 19, on réalise également une opération de segmentation par hystérésis à la suite de la reconstruction tridimensionnelle. La segmentation utilise la technique de seuil par hystérésis qui consiste à sélectionner tous les voxels qui ont une intensité plus grande qu'un seuil bas, 0.012 par exemple, et tels que la zone connexe associée inclut au moins un voxel avec une intensité plus grande qu'un seuil haut, 0.025 par exemple. De plus, le nombre total de voxels dans la zone connexe devrait être plus grand qu'un nombre minimum prédéfini, par exemple 100. Une zone connexe concernant une propriété donnée P, qui peut être le fait qu'un voxel possède une intensité supérieure à un seuil donné, est définie comme un ensemble maximum de voxels tel qu'il est possible d'aller de tout voxel de la zone connexe vers tout autre voxel de la zone connexe par un chemin qui inclus seulement des voxels ayant la propriété P. Le nombre minimum prédéfini est un nombre minimum de voxels dans la zone connexe choisi de façon à rejeter un nombre de voxels qui sont dus aux artéfacts aux limites de champ de vision. La segmentation permet d'éliminer des voxels inutiles tels que ceux concernant les objets non désirés de fort contraste. L'image obtenue à l'issue de la segmentation est une image tridimensionnelle primaire 20 représentant les spires 3 en trois dimensions. Par ailleurs, on réalise une reconstruction d'image tridimensionnelle angiographique 21 de façon à obtenir une image tridimensionnelle 22 contenant des artères du patient. La reconstruction d'image tridimensionnelle angiographique 21 est réalisée de manière classique à partir d'images bidimensionnelles obtenues par une double acquisition séparée par une injection de produit opaque dans les artères. A l'étape 23, on applique un seuil sur l'image tridimensionnelle 22 afin d'éliminer des voxels inutiles, c'est-à-dire ceux qui possèdent une faible intensité, et d'obtenir une image tridimensionnelle limitée 24. On peut conserver un nombre compris par exemple dans un intervalle non limitatif de 1 à 5 millions de voxels ayant les plus grandes intensités. On réalise à l'étape 25 une fusion des images tridimensionnelles 20 et 24 pour obtenir l'image fusionnée 26. Lorsque dans l'image tridimensionnelle limitée 24 un voxel comporte une intensité différente de celle que comporte un voxel équivalent, c'est-à- dire de mêmes coordonnées, dans l'image tridimensionnelle primaire , on attribue au voxel équivalent dans l'image fusionnée 26 l'intensité la plus élevée. L'étape 27 permet de visualiser une image tridimensionnelle dans laquelle on distingue les spires 3 et les artères. La figure 6 montre très schématiquement les différentes images obtenues au cours du procédé illustré par l'organigramme de la figure 5. L'image bidimensionnelle acquise 14 contient la tête du patient 34 dans laquelle se trouve les spires 3. Le filtrage passe-bas 15 produit une image bidimensionnelle 16 ne contenant pas de spires 3. La soustraction logarithmique 17 entre ces deux images 14 et 16 conduit à l'image bidimensionnelle soustraite 18 contenant uniquement les spires 3. La reconstruction tridimensionnelle et la segmentation 19 aboutissent à l'image tridimensionnelle primaire 20. On fusionne l'image tridimensionnelle primaire 20 et l'image tridimensionnelle 24 obtenue par reconstruction d'image tridimensionnelle angiographique pour créer l'image fusionnée tridimensionnelle 26 comportant les spires 3, l'artère 1 et l'anévrisme 2. La figure 7 est un organigramme comportant tous les éléments de la figure 5. Elle présente une variante entre les étapes 17 soustraction logarithmique et 25 fusion. En effet les étapes 17 et 19 sont réalisées deux fois. Une première fois identiquement à la figure 5 par le câblage 16a de façon à obtenir l'image tridimensionnelle primaire 20. Mais, contrairement à la figure 5, l'image tridimensionnelle primaire subit une rétroprojection 28 sur les images bidimensionnelles acquises 14 de façon à déterminer l'emplacement des spires sur ces images bidimensionnelles acquises 14. On réalise également au cours de l'étape 28 une approximation des emplacements. On prend en compte des surfaces éventuellement supérieures aux emplacements ainsi déterminés. L'approximation consiste à remplacer l'intensité de ces surfaces par une intensité obtenue en approximant une fonction polynomiale à partir des intensités des pixels environnants. L'image bidimensionnelle 29 obtenue ne contient pas les spires 3. On réalise donc une seconde fois les étapes 17 et 19 au l1 moyen du câblage 16b de façon à obtenir une image tridimensionnelle secondaire 30 montrant les spires 3 dans les mêmes dispositions spatiales que l'image tridimensionnelle primaire 20, mais avec une meilleure précision concernant les niveaux d'intensité des voxels. Ensuite on effectue une fusion à partir des images tridimensionnelles 30 et 24 avec les mêmes critères que celles utilisées pour la fusion de la figure 5. La figure 8 illustre un procédé se basant sur le procédé de la figure 7. Par rapport à la figure 7, les étapes 21 et 23 concernant la reconstruction d'image tridimensionnelle angiographique ne sont pas réalisées, ainsi que l'étape de fusion 25. L'étape 19 est scindée en deux étapes 19a et 19b. Identiquement au procédé de la figure 7, on effectue une première réalisation des étapes 17, 19a et 19b (étape 19) à partir des images bidimensionnelles acquises 14 et des images bidimensionnelles filtrées 16 au moyen des câblages 14a et 16a pour déterminer l'image tridimensionnelle primaire 20. Puis on effectue une deuxième fois les étapes 17 et 19a, sans l'étape 19b, d'une part, à partir de l'image bidimensionnelle 29 obtenue par rétroprojection- approximation 28 au moyen du câblage 16b, et d'autre part, contrairement au procédé de la figure 7, à partir d'une image bidimensionnelle 32 contenant les spires 3 et les artères au moyen du câblage 14b. Les images bidimensionnelles 32 sont obtenues au cours d'une acquisition 31 à la suite d'une injection de produit opaque dans les artères. La deuxième application des étapes 17 et 19 permet d'obtenir une image tridimensionnelle tertiaire 33 directement visualisable au cours de l'étape 27 car cette image tridimensionnelle 33 montre les spires 3 et les artères. La présente invention concerne un procédé aboutissant à la visualisation d'une image tridimensionnelle des spires et des artères. Un opérateur peut ainsi voir la disposition spatiale des spires dans les artères. REVENDICATIONS 1. Procédé de reconstruction d'une image tridimensionnelle (20) d'éléments de fort contraste (3) à partir d'un ensemble d'images bidimensionnelles (14) d'un objet (34) comprenant lesdits éléments de fort contraste (3), dans lequel pour chaque position différente d'un appareil de prises de vues à rayons X autour de l'objet, une image bidimensionnelle (14) est prise, et la mise en oeuvre d'un algorithme pour la reconstruction de l'image tridimensionnelle est précédée d'une étape de filtrage (15) de l'ensemble d'images bidimensionnelles (14). 2. Procédé selon la revendication précédente, caractérisé par le fait que le filtrage (15) de chaque image bidimensionnelle acquise (14) est du type passe-bas de façon à obtenir une image bidimensionnelle filtrée (16) ne contenant plus les éléments de fort contraste (3), et par le fait qu'on réalise à la suite du filtrage une soustraction logarithmique (17) entre cette image ne contenant plus les éléments de fort contraste et l'image bidimensionnelle acquise de façon à ne conserver principalement que lesdits éléments de fort contraste. 3. Procédé selon l'une des revendications précédentes, caractérisé par le fait que le filtrage passe-bas sur une image bidimensionnelle acquise consiste à: - prendre l'image bidimensionnelle acquise (7, 14), - déterminer sa moyenne (10), - dilater l'image bidimensionnelle acquise de façon à éliminer des objets non désirés (9) de fort contraste et de faible taille par rapport auxdits éléments (3) de fort contraste, et - prendre le maximum (12) entre l'image bidimensionnelle acquise dilatée (11) et la moyenne (10) de l'image bidimensionnelle acquise de façon à obtenir une image bidimensionnelle filtrée exempte de tout élément et objet de fort contraste. 4. Procédé selon l'une des revendications précédentes, caractérisé par le fait que lesdits éléments (3) sont des prothèses. 5. Procédé selon la revendication 4, caractérisé par le fait qu'à la suite de la reconstruction d'image tridimensionnelle contenant les prothèses, on réalise une segmentation par hystérésis (19) de cette image tridimensionnelle de façon à éliminer des voxels inutiles et visualiser uniquement lesdites prothèses (3) dans une image tridimensionnelle primaire (20). 6. Procédé selon la revendication 5, caractérisé par le fait que pour chaque prothèse (3): - on réalise une rétro-projection de l'image tridimensionnelle primaire (20) sur les images bidimensionnelles acquises de façon à déterminer les pixels représentant la prothèse sur ces images, - on effectue une approximation en remplaçant l'intensité de chacun de ces pixels par une intensité déterminée à partir des intensités des pixels voisins afin d'obtenir des images bidimensionnelles primaires (29) exemptes de la prothèse (3), - on réalise une soustraction logarithmique (17) entre les images bidimensionnelles primaires (29) exemptes de la prothèse (3) et les images bidimensionnelles acquises (14) respectives de façon à obtenir de nouvelles images bidimensionnelles (18) comportant uniquement la prothèse, et - on effectue une reconstruction d'image tridimensionnelle suivie d'une segmentation par hystérésis (19) à partir des nouvelles images bidimensionnelles (18) afin d'obtenir une image tridimensionnelle secondaire (30) comportant uniquement la prothèse (3). 7. Procédé selon la revendication 5 ou 6, caractérisé par le fait qu'on effectue indépendamment une reconstruction (21) d'image tridimensionnelle angiographique (21, 24) des artères (1) et on réalise une fusion (25) entre l'image tridimensionnelle primaire (20) ou secondaire (30) contenant la prothèse (3) et l'image tridimensionnelle (24) contenant des artères de façon à obtenir une image tridimensionnelle fusionnée (26) montrant la répartition spatiale entre les artères (1) et la prothèse (3). 8. Procédé selon la revendication 7, caractérisé par le fait qu'au cours de la fusion (26), lorsque dans l'image tridimensionnelle (24) contenant les artères (1) un voxel comporte une intensité différente de celle que comporte le voxel équivalent dans l'image tridimensionnelle (20, 30) contenant la prothèse (3), on attribue au voxel résultant de la fusion l'intensité la plus élevée. 9. Procédé selon la revendication 7, caractérisé par le fait qu'avant la fusion (25) de l'image tridimensionnelle primaire (20) ou secondaire (30) contenant la prothèse (3) avec l'image tridimensionnelle (24) contenant les artères, on élimine sur l'image tridimensionnelle (24) contenant les artères des voxels ayant une faible valeur d'atténuation en ne conservant qu'un nombre prédéterminé de voxels de grande valeur d'atténuation. 10. Procédé selon la revendication 5, caractérisé par le fait que pour chaque prothèse (3): - on réalise une rétro-projection de l'image tridimensionnelle primaire (20) sur les images bidimensionnelles acquises de façon à déterminer les pixels représentant la prothèse sur ces images, - on effectue une approximation en remplaçant l'intensité de chacun de ces pixels par une intensité déterminée à partir des intensités des pixels voisins afin d'obtenir des images bidimensionnelles primaires (29) exemptes de la prothèse (3), - on réalise une soustraction logarithmique (17) entre les images bidimensionnelles primaires (29) exemptes de la prothèse (3) et des images bidimensionnelles (32) contenant des artères et la prothèse (3) acquises après injection de produit opaque aux rayons X de façon à obtenir de nouvelles images bidimensionnelles (18) contenant les artères et la prothèse, et - on effectue une reconstruction d'image tridimensionnelle à partir des nouvelles images bidimensionnelles (18) afin d'obtenir une image tridimensionnelle tertiaire (33) contenant les artères (1) et la prothèse (3).