FR 2812741 A1 20020208 FR 0010211 A 20000802 Procédé et dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste. L'invention concerne la reconstruction d'une image tridimensionnelle d'un objet parcouru par un produit de contraste. Elle trouve une application particulièrement intéressante dans le domaine médical dans lequel on procède à la reconstruction des structures internes de patients 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. L'invention peut néanmoins trouver des applications dans d'autres domaines, notamment dans le contrôle non destructif industriel dans lequel des examens du même type que les examens médicaux sont pratiqués. Dans le domaine médical, les images projetées bidimensionnelles de l'objet, par exemple la tête d'un patient, sont généralement obtenues par la rotation d'un appareil de prises de vues à rayons X tournant autour de l'objet. Il existe essentiellement deux types d'algorithmes de reconstruction en imagerie par rayons X. Un premier type prévoit un calcul de rétroprojection et de filtrage ou encore une reconstruction par transformée de Fourier à plusieurs dimensions. Un deuxième type concerne les méthodes itératives de reconstruction dites encore algébriques. Le principe d'un tel algorithme algébrique est parfaitement connu de l'homme du métier et a déjà fait l'objet de nombreuses publications. On citera notamment les demandes de brevets français n0 89 03 606, 89 16 906, ou encore 98 07 371 qui décrit une application d'un algorithme itératif de reconstruction algébrique d'images sur un volume multi-résolutions. En bref, après un étalonnage de l'appareil utilisé pour déterminer notamment les paramètres de la projection dans les plans de projection des images acquises, d'un volume observé décomposé en éléments volumiques élémentaires ou voxels (ces paramètres d'étalonnage formant des matrices de projection), l'algorithme de reconstruction algébrique d'image est utilisé pour reconstruire le volume tridimensionnel à partir de ces images projetées bidimensionnelles. Le principe de base de cet algorithme est d'initialiser les voxels du volume à une valeur initiale prédéterminée, par exemple une valeur nulle, et d'itérer un certain nombre de fois les opérations suivantes: projection des voxels dans le plan de chaque image acquise de façon à obtenir une image virtuelle, détermination de la différence entre le volume projeté (image virtuelle) et l'image acquise correspondante puis rétroprojection de cette différence dans le volume. Après un certain nombre d'itérations, on obtient pour chaque voxel une valeur estimée représentative de la densité de produit de contraste injecté dans les vaisseaux radiographiés, ce qui permet de visualiser en trois dimensions la cartographie de ces vaisseaux radiographiés. Ces images tridimensionnelles sont d'une aide précieuse pour le neurologue et le chirurgien, que ce soit pour un diagnostic, la planification des gestes thérapeutiques ou l'évaluation de la forme et de la taille d'objets. En revanche, de telles images volumiques reconstruites présentent un inconvénient majeur. En effet, elles ne permettent pas de visualiser la propagation du produit de contraste injecté dans les artères puisque les images bidimensionnelles acquises, et à partir desquelles on a reconstruit l'image tridimensionnelle, sont acquises pour une quasi stationnarité du produit de contraste. En effet, pour la reconstruction d'une image tridimensionnelle, que l'on désignera ici sous le terme de "statique" en raison du caractère quasi stationnaire du produit de contraste, on cherche plutôt à obtenir un ensemble d'images qui correspondent à un même degré de propagation du produit de contraste. Par conséquent, on ne peut pas, avec une telle reconstruction, analyser le réseau vasculaire du patient en disposant à la fois d'une information anatomique tridimensionnelle et d'une information de flux sanguin (propagation du produit de contraste). L'invention vise à remédier à cette lacune. L'invention a pour but de reconstruire une image tridimensionnelle, que l'on qualifiera ici de "dynamique", par opposition à l'image tridimensionnelle statique évoquée ci-avant, permettant de visualiser la propagation en trois dimensions du produit de contraste dans l'objet. L'invention propose donc un procédé de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste, ce procédé comprenant - une phase de reconstruction (ou de reconstitution) d'une image tridimensionnelle statique de l'objet à partir d'un premier ensemble d'images projetées bidimensionnelles numérisées dudit objet respectivement obtenues pour différents positions d'un appareil de prises de vues autour de l'objet, - une phase d'acquisition d'au moins un deuxième ensemble de n images projetées bidimensionnelles statiques respectivement obtenues pour une même position de l'appareil de prises de vues et à n instants successifs de propagation du produit de contraste, et - une phase de reconstruction de l'image tridimensionnelle dynamique de l'objet, à partir de chaque image bidimensionnelle statique dudit deuxième ensemble et de l'image tridimensionnelle statique reconstruite (ou reconstituée). Selon un mode de mise en oeuvre du procédé, celui-ci comprend un étalonnage de l'appareil dans lequel on élabore un volume virtuel entourant l'objet et décomposé en voxels. L'image tridimensionnelle statique reconstruite est alors formée de valeurs estimées, dites valeurs statiques estimées, respectivement associées aux voxels du volume virtuel, chaque valeur statique estimée étant représentative pour le voxel correspondant de la densité de produit de contraste injecté dans l'objet. L'image tridimensionnelle dynamique se compose en fait de n images tridimensionnelles élémentaires correspondant aux n instants de propagation du produit de contraste. Et, une image tridimensionnelle élémentaire courante correspondant à un instant courant de propagation est reconstruite à partir - desdites valeurs statiques estimées de la densité de produit de contraste injecté dans l'objet, - pour chaque voxel, de l'intensité du pixel de l'image bidimensionnelle statique correspondant audit instant de propagation courant, sur lequel se projette ledit voxel, et - pour chacun de ces pixels, de l'intégrale des valeurs estimées de densité des voxels du volume virtuel situés le long de la ligne de prise de vue associée à ce pixel. Dans certains cas, plusieurs vaisseaux peuvent se projeter sur le même pixel. Ceci correspond généralement à 10% des voxels. Il peut en résulter une ambiguité dans l'estimation de l'image tridimensionnelle dynamique reconstruite. Une telle ambiguité peut être résolue en utilisant une autre série d'images bidimensionnelles statiques acquises depuis un point de vue différent de celui des n images bidimensionnelles du deuxième ensemble. En d'autres termes, selon un autre mode de mise en oeuvre, il est prévu une phase d'acquisition d'un troisième ensemble d'images projetées bidimensionnelles respectivement obtenues pour une même deuxième position de l'appareil de prises de vues, distincte de ladite première position (celle correspondant à l'acquisition des images projetées bidimensionnelles du deuxième ensemble), et aux mêmes intants successifs de propagation du produit de contraste. On reconstruit alors l'image tridimensionnelle dynamique de l'objet, à partir de l'image tridimensionnelle statique reconstituée, de chaque image dudit deuxième ensemble et de chaque image dudit troisième ensemble. L'invention propose également un dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste, comprenant - des premiers moyens de reconstruction aptes à reconstruire une image tridimensionnelle statique de l'objet à partir d'un premier ensemble d'images projetées bidimensionnelles numérisées dudit objet, respectivement obtenues pour différentes positions d'un appareil de prises de vues autour de l'objet, - des moyens d'acquisition aptes à acquérir au moins un deuxième ensemble de n images projetées bidimensionnelles statiques depuis une même première position de l'appareil de prises de vues et correspondant à n instants successifs de propagation du produit de contraste, et - des deuxièmes moyens de reconstruction aptes à reconstruire l'image tridimensionnelle dynamique de l'objet, à partir de chaque image bidimensionnelle statique dudit deuxième ensemble et de l'image tridimensionnelle statique reconstruite. L'invention a également pour objet un dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste, adapté pour mettre en oeuvre le procédé selon l'invention, tel que défini ci-avant. Les différents moyens permettant la mise en oeuvre du procédé sont des moyens logiciels exécutés par un microprocesseur. Ces moyens logiciels peuvent être disposés sur un support, tel qu'une mémoire morte et/ou une disquette, et ils peuvent être aisément réalisés par l'homme du métier à partir de la description des étapes du procédé. L'invention propose donc également un programme d'ordinateur comprenant des moyens de code-programme mettant en oeuvre le procédé selon l'invention tel que défini ci-avant, lorsque ledit programme est exécuté au sein d'un processeur. L'invention propose encore un support, tel qu'une mémoire morte ou une disquette, capable d'être lue par un processeur, et contenant des moyens de code-programme aptes à mettre en oeuvre le procédé selon l'invention lorsque ledit programme est exécuté au sein du processeur. D'autres avantages et caractéristiques de l'invention apparaîtront à l'examen de la description détaillée de modes de mise en oeuvre et de réalisation, nullement limitatifs, et des dessins annexés, sur lesquels: - la figure 1 illustre schématiquement un jeu d'images projetées bidimensionnelles autour d'un objet; - la figure 2 illustre plus en détail l'acquisition de ces images projetées bidimensionnelles; - la figure 3 illustre schématiquement un jeu d'images statiques projetées acquises depuis un même point de vue; - les figures 4 et 5 représentent un organigramme de mode de mise en oeuvre du procédé selon l'invention; et la figure 6 illustre schématiquement une variante de mise en oeuvre de l'invention. Bien que l'invention n'y soit pas limitée, on va maintenant décrire l'application du procédé selon l'invention à la reconstruction d'une image angiographique tridimensionnelle dynamique d'un patient, en particulier sa tête. Si l'on se réfère plus particulièrement aux figures 1 et 2, on voit que le système d'imagerie utilisable pour mettre en oeuvre l'invention permet d'obtenir notamment un ensemble d'images acquises bidimensionnelles IA 1 IAm obtenues en l'espèce par la rotation autour de la tête I d'un patient d'une source 2 à rayons X. En fait, comme il est classique en angiographie, chaque image acquise IAi 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, tout d'abord statique puis dynamique au sens de la présente invention. Plus précisément, deux séquences sont acquises par rotation autour de la tête du patient, avant et après injection du produit de contraste. Pour ce faire, l'appareil de radiographie effectue par exemple deux rotations de 193,5 en 10 secondes environ, permettant l'acquisition de 44 images bidimensionnelles effectuées tous les 4,5 . Cette acquisition des images IAi est destinée à la reconstruction tridimensionnelle d'une image de l'arbre vasculaire, dite statique, c'est-à-dire dans laquelle on ne s'intéresse pas à la propagation du produit de contraste. Dans ces conditions, cette acquisition est effectuée de préférence lorsque le produit de contraste est quasi stationnaire dans les vaisseaux étudiés. Chaque image acquise IAi est obtenue à partir d'un détecteur de rayonnement bidimensionnel, par exemple du type amplificateur de luminescence 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. La normale XOi au plan de projection PPi définit l'axe optique de l'image acquise IAi. Le détecteur est relié à des moyens de traitement 3 comportant notamment des moyens d'échantillonnage reliés à un microprocesseur incorporant de façon logicielle dans sa mémoire de programme associée, l'algorithme de reconstruction d'image 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, tels que les différents moyens de reconstruction d'images. Dans le cas d'un système d'imagerie à rayons X, composé d'une source de rayons X à détecteur bidimensionnel, l'opération géométrique intervenant dans la production de l'image acquise est une projection conique d'un objet analysé, déployée dans un espace à trois dimensions, sur un espace à deux dimensions qui est celui du plan de projection correspondant au plan de détection. Les paramètres géométriques décrivant complètement les diverses projections coniques doivent être connus. Or, il est souvent impossible et trop imprécis d'accéder à ces paramètres de façon directe, c'est-à-dire, par exemple, en mesurant directement sur le système d'acquisition la distance entre la source de rayons X et le détecteur. On appelle "étalonnage" d'un système d'imagerie, l'opération qui aboutit à la connaissance indirecte précise des paramètres géométriques qui interviennent dans la production d'une image. Le principe, classique et connu, est basé sur l'utilisation d'un fantôme géométrique connu dans l'espace tridimensionnel, et dont on acquiert la projection bidimensionnelle. Plus précisément, l'étalonnage comporte les étapes suivantes: - on dispose d'un objet connu, le fantôme d'étalonnage, présentant un certain nombre de points caractéristiques dont la position dans l'espace est connue par des coordonnées mesurées par rapport à un repère propre à cet objet; - on acquiert l'image de ce fantôme dans les conditions géométriques d'un point de vue (ou incidence) que l'on veut étalonner; - on reconnaît les projections des points caractéristiques dans l'image. Pour cela, on associe chaque point caractéristique de l'objet à sa trace dans l'image acquise projetée; - on inverse, au sens mathématique, le système d'équation décrivant la projection; - et on obtient finalement l'ensemble des paramètres de la projection pour le point de vue donné. Une forme de fantôme d'étalonnage géométrique souvent utilisée est celle d'un cube, aux huit coins duquel sont disposées des billes métalliques opaques aux rayons X. L'étalonnage étant une opération connue de l'homme du métier, elle ne sera pas décrite plus en détail. On peut néanmoins citer à cet égard 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 00 804. En bref, on utilise pour un tel étalonnage automatique un fantôme dans lequel les billes sont distribuées, de proche en proche, en une succession telle que des altitudes de billes, mesurées le long de l'axe de rotation du système d'imagerie, et surtout un axe du fantôme, soient monotones croissantes (ou décroissantes) avec un numéro d'ordre des billes dans la succession. L'étalonnage du système d'magerie permet notamment de déterminer l'axe moyen estimé Ax de rotation de l'appareil de prises de vues autour de l'objet, ainsi que la position de la source 2 et les caractéristiques géométriques des axes optiques des différentes images acquises. L'étalonnage permet également de définir un volume virtuel VV (intersection des différents cônes de projection) entourant l'objet 1 et décomposé en éléments volumiques élémentaires Vi ou "voxels". Ce volume VV, et donc chaque voxel Vi, est spatialement repéré dans un référentiel, dénommé ci-après référentiel de calibration, dont l'un des axes, en l'espèce l'axe Z, est confondu avec l'axe de rotation estimé Ax. Il convient de noter ici que les plans de projection PPi dans lesquels se projettent les images acquises IAi ne sont généralement pas parallèles à l'axe Z. L'étalonnage permet également de définir pour chaque image acquise IAi une matrice de projection Pi permettant de déterminer, pour chaque voxel Vi, les coordonnées de sa projection (pixel) dans l'image acquise IAi correspondante. On va maintenant décrire, en se référant plus particulièrement à la figure 4, la phase de reconstruction de l'image tridimensionnelle dite "statique", à partir du premier ensemble d'images bidimensionnelles statiques projetées obtenu IAi. On va appliquer directement sur les images IAi (obtenues après étalonnage 30 et acquisition 31) un algorithme de reconstruction algébrique d'image de façon à délivrer l'image tridimensionnelle statique reconstituée IF, typiquement après trois itérations. On rappelle ici que le principe de base de l'algorithme de reconstruction d'image utilisé ici consiste à initialiser les voxels du volume à une valeur initiale prédéterminée, par exemple la valeur nulle, et d'effectuer ensuite un certain nombre d'itérations, typiquement trois. Chaque itération comporte, pour chaque image acquise, la projection du voxel dans le plan de chaque image acquise de façon à obtenir une image virtuelle, la détermination de la différence entre le volume projeté (image virtuelle) et l'image acquise correspondante, puis la rétroprojection de cette différence dans le volume. Généralement, après trois itérations, on obtient pour chaque voxel de coordonnées (i, j, k) une valeur estimée cijk, désignée ici sous le vocable de valeur estimée statique, représentative de la densité de produit de contraste injecté dans les vaisseaux radiographiés, ce qui permet de visualiser en trois dimensions la cartographie de ces vaisseaux radiographiés (image IF). Bien entendu, dans le cadre de la présente invention, il n'est pas nécessaire de visualiser sur l'écran de visualisation du système d'imagerie, l'image tridimensionnelle statique IF ainsi obtenue puisqu'en fait elle ne consiste ici qu'en une étape intermédiaire pour obtenir, comme on va le voir maintenant plus en détail, l'image tridimensionnelle dynamique reconstruite. L'homme du métier pourra se référer, s'il le souhaite, pour plus de détails concernant la mise en oeuvre d'un tel algorithme algébrique de reconstitution d'image, à la demande de brevet français n 98 07 371 précitée. Outre l'acquisition du premier ensemble d'images IAi, on procède, selon l'invention, à l'acquisition d'un deuxième ensemble d'images projetées bidirectionnelles statiques IBI-IBn (figure 3). Cette série d'images IBi est acquise pour une même position de l'appareil de prises de vues, correspondant à un même point de vue PDV1. Les n images IBt ainsi acquises correspondent à n instants de propagation successifs du produit de contraste. Comme illustré sur la figure 5, c'est à partir de ces images IBi acquises au cours d'une étape d'acquisition 50, et à partir de l'image tridimensionnelle statique reconstituée IF que l'on va, dans une étape 51, reconstituer l'image dynamique tridimensionnelle de l'arbre vasculaire. Il convient de noter ici que l'ordre dans lequel sont effectuées, après étalonnage, les acquisitions 31 (figure 4) et 50 (figure 5), est sans importance pour la mise en oeuvre du procédé selon l'invention. L'image dynamique tridimensionnelle reconstituée IFD se compose en fait de n images dynamiques tridimensionnelles élémentaires sur chacune desquelles on va pouvoir visualiser l'état de propagation à l'instant de propagation t donné, du produit de contraste dans l'arbre vasculaire. En visualisant successivement les n images tridimensionnelles élémentaires, on va pouvoir visualiser de façon dynamique la propagation du produit de contraste dans l'arbre vasculaire. D'une façon analogue à ce qui a été décrit ci-avant pour l'image tridimensionnelle reconstituée statique, la visualisation de l'image tridimensionnelle reconstituée dynamique est obtenue par une valeur estimée, pour chaque voxel de coordonnées (i,j,k), de la densité de produit de contraste à l'instant de propagation t. Cette valeur estimée valable à l'instant de propagation t, est désignée par cijk(t). Et, cette valeur estimée est définie par la formule (1) ci-dessous: cijk. p(t) Ci'j'k(t) = E (1) E dans laquelle cijk désigne la valeur statique estimée de produit de contraste de ce même voxel dans l'image tridimensionnelle reconstituée statique IF, p(t) désigne l'intensité du pixel de l'image IBt acquise à l'instant de propagation t et sur lequel se projette le voxel de coordonnées i,j,k, et il E désigne pour ce pixel p, l'intégrale des valeurs statiques estimées de densité des voxels du volume virtuel VV situés le long de la ligne de prise de vue associée à ce pixel p. E est défini par la formule (2) ci-dessous E = z cij,k ij,k (2) ij,k dans laquelle Xi,j,k vaut 0 si le voxel de coordonnées i,j,k ne se situe pas sur la ligne de prise de vue et vaut I s'il se situe sur la ligne de prise de vue. Il convient de noter ici que la détermination du pixel p est possible en raison de la connaissance de la matrice de projection associée à la série d'images statiques prises selon une ligne de vue prédéterminée. De même, les valeurs X dans la formule (2) sont également connues à partir des valeurs de cette matrice de projection. Il convient de noter ici que si à l'instant t, le voxel de coordonnées i, j,k n'est pas encore atteint par le produit de contraste, alors on doit avoir une valeur faible pour cijk(t). Si la valeur p(t) est faible, c'està-dire si la densité du pixel projeté est faible, la formule (1) va donc bien conduire à une densité de voxel faible à l'instant t. En revanche, si la valeur p(t) est élevée, cela signifie que la ligne de vue concernée traverse d'autres voxels, lumineux, qui viennent parasiter la densité due au voxel concerné de coordonnées i,j,k. La formule (1) rend alors le voxel plus atteint que ce qu'il n'est réellement. Si le voxel de coordonnées i,j,k est atteint par le produit de contraste, la formule (1) peu créer des artefacts comme dans le cas précédent, si le voxel considéré se superpose avec un voxel non encore atteint par le produit de contraste. Cependant, la formule (1) ci-dessus permet une très bonne estimation de l'image reconstituée dynamique, étant donné que les conditions de création d'artefacts mentionnées ci-dessus ne concernent environ que 10% des voxels du volume. Cependant, afin de résoudre cette ambiguïté et d'améliorer par conséquent la reconstruction dynamique de l'image, on peut, comme mentionné sur la figure 6, procéder à l'acquisition d'une autre série de n images bidimensionnelles statiques projetées ici, à partir d'un point de vue différent du point de vue relatif à la première série d'images IBi. Les images IBi ICi sont respectivement acquises aux mêmes instants t. Dans ce cas, comme indiqué sur la figure 5 en tiretés, la reconstitution de l'image dynamique tridimensionnelle dans l'étape 51 s'effectue également à partir des images ICi acquises à l'étape 52. On calculera alors une première valeur cijk(t) pour le voxel i,j,k à l'instant t, en utilisant les images IBi, et une deuxième valeur cik(t) en utilisant cette fois-ci la formule (1) en liaison avec les images ICi. Et, la valeur finale cidk(t) associée au voxel de coordonnées i,j,k sera alors égale à la plus petite des deux valeurs mentionnées ci-avant. REVENDICATIONS 1. Procédé de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste, comprenant une phase de reconstruction (33) d'une image tridimensionnelle statique (IF) de l'objet à partir d'un premier ensemble d'images projetées bidimensionnelles numérisées (IAi) dudit objet respectivement obtenues pour différentes positions d'un appareil de prises de vues autour de l'objet, une phase d'acquisition (50) d'au moins un deuxième ensemble de n images projetées bidimensionnelles statiques (IBi) respectivement obtenues pour une même première position de l'appareil de prises de vues (PDV1) et à n instants successifs de propagation du produit de contraste, et une phase de reconstruction (51) de l'image tridimensionnelle dynamique (IFD) de l'objet, à partir de chaque image bidimensionnelle statique (IBi) dudit deuxième ensemble et de l'image tridimensionnelle statique reconstruite (IF). 2. Procédé selon la revendication 1, caractérisé par le fait qu'il comprend un étalonnage (30) de l'appareil dans lequel on élabore un volume virtuel (VV) entourant l'objet et décomposé en voxels (Vi), par le fait que l'image tridimensionnelle statique reconstruite (IF) est formée de valeurs statiques estimées (cijk) respectivement associées aux voxels du volume virtuel, chaque valeur statique estimée étant représentative pour le voxel correspondant de la densité de produit de contraste injecté dans l'objet, par le fait que l'image tridimensionnelle dynamique (IFD) se compose de n images tridimensionnelles élémentaires correspondant aux n instants de propagation du produit de contraste, et par le fait qu'une image tridimensionnelle élémentaire courante correspondant à un instant courant de propagation est reconstruite à partir - desdites valeurs statiques estimées (cijk) de la densité de produit de contraste injecté dans l'objet, - pour chaque voxel (i,j,k) de l'intensité du pixel (p(t)) de l'image bidimensionnelle statique correspondant audit instant de propagation courant, sur lequel se projette ledit voxel (i,j,k), et - pour chacun de ces pixels, de l'intégrale (E) des valeurs estimées de densité des voxels du volume virtuel situés le long de la ligne de prise de vue associée à ce pixel. 3. Procédé selon la revendication 1 ou 2, caractérisé par le fait qu'il comprend en outre une phase d'acquisition (52) d'un troisième ensemble d'images projetées bidimensionnelles statiques (ICi) respectivement obtenues pour une même deuxième position de l'appareil de prises de vues, distincte de ladite première position, et aux mêmes instants successifs de propagation du produit de contraste, et par le fait qu'on reconstruit l'image tridimensionnelle dynamique de l'objet, à partir de l'image tridimensionnelle statique, de chaque image dudit deuxième ensemble et de chaque image dudit troisième ensemble. 4. Dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste, comprenant des premiers moyens de reconstruction (3) aptes à reconstruire une image tridimensionnelle statique (IF) de l'objet à partir d'un premier ensemble d'images projetées bidimensionnelles numérisées (IAi) dudit objet respectivement obtenues pour différentes positions d'un appareil de prises de vues autour de l'objet, des moyens d'acquisition aptes à acquérir au moins un deuxième ensemble de n images projetées bidimensionnelles statiques (IBi) depuis une même première position de l'appareil de prises de vues et correspondant à n instants successifs de propagation du produit de contraste, et des deuxièmes moyens de reconstruction (3) aptes à reconstruire l'image tridimensionnelle dynamique de l'objet, à partir de chaque image bidimensionnelle statique dudit deuxième ensemble et de l'image tridimensionnelle statique reconstruite. 5. Dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste, adapté pour mettre en oeuvre le procédé selon l'une des revendications 1 à 3.