FR2812741A1 - Procede et dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste - Google Patents
Procede et dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste Download PDFInfo
- Publication number
- FR2812741A1 FR2812741A1 FR0010211A FR0010211A FR2812741A1 FR 2812741 A1 FR2812741 A1 FR 2812741A1 FR 0010211 A FR0010211 A FR 0010211A FR 0010211 A FR0010211 A FR 0010211A FR 2812741 A1 FR2812741 A1 FR 2812741A1
- Authority
- FR
- France
- Prior art keywords
- static
- dimensional
- dimensional image
- image
- dynamic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 230000003068 static effect Effects 0.000 claims abstract description 52
- 101100407422 Arabidopsis thaliana PDV1 gene Proteins 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 description 9
- 238000003384 imaging method Methods 0.000 description 8
- 239000002872 contrast media Substances 0.000 description 7
- 230000002792 vascular Effects 0.000 description 6
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- QXXSGHXNHITGDL-UHFFFAOYSA-N CC(C)CC(C)Cc1ccccc1 Chemical compound CC(C)CC(C)Cc1ccccc1 QXXSGHXNHITGDL-UHFFFAOYSA-N 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000002601 radiography Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/006—Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/583—Calibration using calibration phantoms
- A61B6/584—Calibration using calibration phantoms determining position of components of the apparatus or device using images of the phantom
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/404—Angiography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Optics & Photonics (AREA)
- Algebra (AREA)
- High Energy & Nuclear Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Biophysics (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Generation (AREA)
- Image Analysis (AREA)
Abstract
Le procédé comprend une phase de reconstruction d'une image tridimensionnelle statique IF 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, 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 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.
Description
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.
Claims (5)
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.
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0010211A FR2812741B1 (fr) | 2000-08-02 | 2000-08-02 | Procede et dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste |
| JP2001232982A JP4745551B2 (ja) | 2000-08-02 | 2001-08-01 | 造影剤に覆われた被検体のダイナミック３次元画像を再構成するための方法及びデバイス |
| DE10137521A DE10137521A1 (de) | 2000-08-02 | 2001-08-01 | Verfahren und Vorrichtung zur Rekonstruktion eines dynamischen dreidimensionalen Bilds eines von einem Kontrastmittel durchlaufenen Objekts |
| US09/920,586 US7113631B2 (en) | 2000-08-02 | 2001-08-01 | Method and device for reconstruction of a dynamic three-dimensional image of an object covered by a contrast medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0010211A FR2812741B1 (fr) | 2000-08-02 | 2000-08-02 | Procede et dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2812741A1 true FR2812741A1 (fr) | 2002-02-08 |
| FR2812741B1 FR2812741B1 (fr) | 2003-01-17 |
Family
ID=8853236
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR0010211A Expired - Lifetime FR2812741B1 (fr) | 2000-08-02 | 2000-08-02 | Procede et dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US7113631B2 (fr) |
| JP (1) | JP4745551B2 (fr) |
| DE (1) | DE10137521A1 (fr) |
| FR (1) | FR2812741B1 (fr) |
Families Citing this family (35)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4889903B2 (ja) * | 2000-03-30 | 2012-03-07 | ウイスコンシン アラムナイ リサーチ ファウンデーシヨン | 患者の磁気共鳴血管造影図を作成するためのｍｒｉシステムの作動方法 |
| DE10215808B4 (de) * | 2002-04-10 | 2005-02-24 | Siemens Ag | Verfahren zur Registrierung für navigationsgeführte Eingriffe |
| CA2550842A1 (fr) * | 2003-12-30 | 2005-07-21 | The Trustees Of The Stevens Institute Of Technology | Systeme d'imagerie tridimensionnel utilisant des impulsions optiques, des melangeurs optiques non lineaires et l'etalonnage holographique |
| US7693318B1 (en) | 2004-01-12 | 2010-04-06 | Pme Ip Australia Pty Ltd | Method and apparatus for reconstruction of 3D image volumes from projection images |
| CN1664861A (zh) * | 2004-09-22 | 2005-09-07 | 冯彦 | 一种动态标识的制作及显示方法 |
| US8189002B1 (en) | 2004-10-29 | 2012-05-29 | PME IP Australia Pty, Ltd. | Method and apparatus for visualizing three-dimensional and higher-dimensional image data sets |
| US7778392B1 (en) | 2004-11-02 | 2010-08-17 | Pme Ip Australia Pty Ltd | Method of reconstructing computed tomography (CT) volumes suitable for execution on commodity central processing units (CPUs) and graphics processors, and apparatus operating in accord with those methods (rotational X-ray on GPUs) |
| US7609884B1 (en) | 2004-12-23 | 2009-10-27 | Pme Ip Australia Pty Ltd | Mutual information based registration of 3D-image volumes on GPU using novel accelerated methods of histogram computation |
| EP1846753B1 (fr) * | 2005-02-04 | 2012-11-07 | Koninklijke Philips Electronics N.V. | Systeme de determination de la geometrie vasculaire et des caracteristiques de circulation |
| US7623732B1 (en) | 2005-04-26 | 2009-11-24 | Mercury Computer Systems, Inc. | Method and apparatus for digital image filtering with discrete filter kernels using graphics hardware |
| US8411914B1 (en) * | 2006-11-28 | 2013-04-02 | The Charles Stark Draper Laboratory, Inc. | Systems and methods for spatio-temporal analysis |
| DE602008003941D1 (de) * | 2007-04-13 | 2011-01-27 | Brainlab Ag | Verfahren und Vorrichtung für die Planung lokaler Verabreichung eines Arzneimittels |
| US7853061B2 (en) * | 2007-04-26 | 2010-12-14 | General Electric Company | System and method to improve visibility of an object in an imaged subject |
| DE102007024451B4 (de) * | 2007-05-25 | 2010-07-08 | Siemens Ag | Ermittlungsverfahren für örtlich dreidimensional aufgelöste Anwesenheitsverteilungen einer Substanz in einem Gefäßsystem und hiermit korrespondierende Einrichtungen |
| US8019151B2 (en) * | 2007-06-11 | 2011-09-13 | Visualization Sciences Group, Inc. | Methods and apparatus for image compression and decompression using graphics processing unit (GPU) |
| US8392529B2 (en) | 2007-08-27 | 2013-03-05 | Pme Ip Australia Pty Ltd | Fast file server methods and systems |
| WO2009067680A1 (fr) | 2007-11-23 | 2009-05-28 | Mercury Computer Systems, Inc. | Procédés et appareil de segmentation automatique d'image |
| US9904969B1 (en) | 2007-11-23 | 2018-02-27 | PME IP Pty Ltd | Multi-user multi-GPU render server apparatus and methods |
| US10311541B2 (en) | 2007-11-23 | 2019-06-04 | PME IP Pty Ltd | Multi-user multi-GPU render server apparatus and methods |
| US8319781B2 (en) | 2007-11-23 | 2012-11-27 | Pme Ip Australia Pty Ltd | Multi-user multi-GPU render server apparatus and methods |
| WO2009067675A1 (fr) | 2007-11-23 | 2009-05-28 | Mercury Computer Systems, Inc. | Système de visualisation client-serveur à traitement de données hybride |
| US8693745B2 (en) * | 2009-05-04 | 2014-04-08 | Duke University | Methods and computer program products for quantitative three-dimensional image correction and clinical parameter computation in optical coherence tomography |
| US8538106B2 (en) * | 2009-10-20 | 2013-09-17 | Siemens Aktiengesellschaft | Three-dimensional esophageal reconstruction |
| US8787521B2 (en) * | 2009-12-23 | 2014-07-22 | General Electric Company | System and method of iterative image reconstruction for computed tomography |
| ES2392229B1 (es) | 2010-08-27 | 2013-10-16 | Telefónica, S.A. | Método de generación de un modelo de un objeto plano a partir de vistas del objeto. |
| EP2737852B1 (fr) * | 2012-11-30 | 2015-08-19 | GE Sensing & Inspection Technologies GmbH | Procédé de détection des propriétés d'imagerie géométriques d'un détecteur à panneau plat, systèm de test adapté suivant et corps de calibrage |
| US11244495B2 (en) | 2013-03-15 | 2022-02-08 | PME IP Pty Ltd | Method and system for rule based display of sets of images using image content derived parameters |
| US8976190B1 (en) | 2013-03-15 | 2015-03-10 | Pme Ip Australia Pty Ltd | Method and system for rule based display of sets of images |
| US9509802B1 (en) | 2013-03-15 | 2016-11-29 | PME IP Pty Ltd | Method and system FPOR transferring data to improve responsiveness when sending large data sets |
| US10540803B2 (en) | 2013-03-15 | 2020-01-21 | PME IP Pty Ltd | Method and system for rule-based display of sets of images |
| US11183292B2 (en) | 2013-03-15 | 2021-11-23 | PME IP Pty Ltd | Method and system for rule-based anonymized display and data export |
| US10070839B2 (en) | 2013-03-15 | 2018-09-11 | PME IP Pty Ltd | Apparatus and system for rule based visualization of digital breast tomosynthesis and other volumetric images |
| US9984478B2 (en) | 2015-07-28 | 2018-05-29 | PME IP Pty Ltd | Apparatus and method for visualizing digital breast tomosynthesis and other volumetric images |
| US11599672B2 (en) | 2015-07-31 | 2023-03-07 | PME IP Pty Ltd | Method and apparatus for anonymized display and data export |
| US10909679B2 (en) | 2017-09-24 | 2021-02-02 | PME IP Pty Ltd | Method and system for rule based display of sets of images using image content derived parameters |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5287273A (en) * | 1990-03-15 | 1994-02-15 | Mount Sinai School Of Medicine | Functional organ images |
| US5827187A (en) * | 1996-04-23 | 1998-10-27 | Mayo Foundation For Medical Education And Research | Dynamic MR digital subtraction angiography with complex subtraction |
| WO1999027382A1 (fr) * | 1997-11-21 | 1999-06-03 | Cornell Research Foundation, Inc. | Procede d'angiographie par resonance magnetique avec echantillonnage dynamique de l'espace de moment |
| WO2000037957A1 (fr) * | 1998-12-18 | 2000-06-29 | Wisconsin Alumni Research Foundation | Detection de contrastes et reconstitution guidee en arm a contrastes ameliores |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2644590B1 (fr) * | 1989-03-20 | 1994-08-19 | General Electric Cgr Sa | Procede d'acquisition de donnees radiologiques et de reconstruction de structures correspondant a ce corps |
| FR2656129B1 (fr) * | 1989-12-20 | 1992-03-13 | Gen Electric Cgr | Procede de reconstruction multi-echelle de l'image de la structure d'un corps. |
| JP2622064B2 (ja) * | 1992-04-02 | 1997-06-18 | 株式会社東芝 | Ｃｔ装置 |
| EP0587334B1 (fr) * | 1992-09-09 | 1999-06-30 | Picker International, Inc. | Procédé d'imagerie et appareil associé |
| FR2700909B1 (fr) * | 1993-01-27 | 1995-03-17 | Gen Electric Cgr | Dispositif et procédé automatique de calibration géométrique d'un système d'imagerie par rayons X. |
| US6073042A (en) * | 1997-09-25 | 2000-06-06 | Siemens Medical Systems, Inc. | Display of three-dimensional MRA images in which arteries can be distinguished from veins |
| FR2779853B1 (fr) * | 1998-06-11 | 2000-08-11 | Ge Medical Syst Sa | Procede de reconstruction d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique |
| EP1419484A2 (fr) * | 2001-08-10 | 2004-05-19 | Koninklijke Philips Electronics N.V. | Appareil d'examen aux rayons x de reconstruction d'un ensemble tridimensionnel de donnees a partir d'images de projection |
-
2000
- 2000-08-02 FR FR0010211A patent/FR2812741B1/fr not_active Expired - Lifetime
-
2001
- 2001-08-01 US US09/920,586 patent/US7113631B2/en active Active
- 2001-08-01 DE DE10137521A patent/DE10137521A1/de not_active Ceased
- 2001-08-01 JP JP2001232982A patent/JP4745551B2/ja not_active Expired - Lifetime
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5287273A (en) * | 1990-03-15 | 1994-02-15 | Mount Sinai School Of Medicine | Functional organ images |
| US5827187A (en) * | 1996-04-23 | 1998-10-27 | Mayo Foundation For Medical Education And Research | Dynamic MR digital subtraction angiography with complex subtraction |
| WO1999027382A1 (fr) * | 1997-11-21 | 1999-06-03 | Cornell Research Foundation, Inc. | Procede d'angiographie par resonance magnetique avec echantillonnage dynamique de l'espace de moment |
| WO2000037957A1 (fr) * | 1998-12-18 | 2000-06-29 | Wisconsin Alumni Research Foundation | Detection de contrastes et reconstitution guidee en arm a contrastes ameliores |
Non-Patent Citations (1)
| Title |
|---|
| CHENEVERT T L ET AL: "DYNAMIC THREE-DIMENSIONAL IMAGING WITH PARTIAL K-SPACE SAMPLING: INITIAL APPLICATION FOR GADOLINIUM-ENHANCED RATE CHARACTERIZATION OF BREAST LESIONS", RADIOLOGY,US,OAK BROOK,IL, vol. 196, 1 July 1995 (1995-07-01), pages 135 - 142, XP002035246, ISSN: 0033-8419 * |
Also Published As
| Publication number | Publication date |
|---|---|
| US7113631B2 (en) | 2006-09-26 |
| JP4745551B2 (ja) | 2011-08-10 |
| JP2002177254A (ja) | 2002-06-25 |
| DE10137521A1 (de) | 2002-02-14 |
| FR2812741B1 (fr) | 2003-01-17 |
| US20020123680A1 (en) | 2002-09-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| FR2812741A1 (fr) | Procede et dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste | |
| EP0925556B1 (fr) | Procede de reconstruction d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique | |
| EP0379399B1 (fr) | Procédé de calcul et d'exploitation de l'image en projection conique, par exemple au sens des rayons x, d'un objet tridimensionnel echantillonné, et procédé de reconstruction tridimensionnelle d'un objet étudié utilisant ce procédé de calcul | |
| EP0389333B1 (fr) | Procédé d'acquisition de données radiologiques relatives à un corps irradié et de reconstruction de structures correspondant à ce corps | |
| FR2779853A1 (fr) | Procede de reconstruction d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique | |
| EP2453408B1 (fr) | Procédé de traitement d'images radiologiques pour la détection d'une sténose | |
| FR2799028A1 (fr) | Procede de reconstitution d'une image tridimentionnelle d'elements de fort contraste | |
| FR2802002A1 (fr) | Procede de recalage automatique d'images tridimensionnelles | |
| FR2700909A1 (fr) | Dispositif et procédé automatique de calibration géométrique d'un système d'imagerie par rayons X. | |
| FR2953964A1 (fr) | Procede de traitement d'images obtenues par tomographie ou tomosynthese a faible nombre de projections | |
| FR2799031A1 (fr) | Procede de reconstruction d'une section, par exemple transversale, d'un element d'interet contenu dans un objet, en particulier un vaisseau du coeur humain | |
| FR2954556A1 (fr) | Procede de traitement d'acquisitions de tomosynthese pour obtenir une representation du contenu d'un organe | |
| EP0840252A1 (fr) | Procédé de traitement d'image numérique pour l'extraction automatique d'objets en forme de rubans | |
| FR2963976A1 (fr) | Procede de traitement d'images pour la determination de zones suspectes dans une matrice tissulaire, et son utilisation pour la navigation 3d a travers la matrice tissulaire | |
| FR2797978A1 (fr) | Procede de recalage automatique d'images | |
| FR2848007A1 (fr) | Procede et systeme pour accentuer une image de tomosynthese en utilisant un filtrage transversal | |
| FR2967520A1 (fr) | Procede de traitements d'images radiologiques d'un patient | |
| EP0752684A1 (fr) | Procédé de reconstruction d'images tridimensionnelles sur un objet mobile ou déformable | |
| EP1417646B1 (fr) | Procede de reconstruction d'une image d'un objet evoluant | |
| FR2965651A1 (fr) | Reconstruction tomographique d'un objet en mouvement | |
| FR2799029A1 (fr) | Procede de reconstruction d'une image tridimensionnelle d'un objet en mouvement, en particulier une image tridimensionnelle de vaisseaux du coeur humain | |
| FR2944203A1 (fr) | Systeme et procede pour determiner la position d'un instrument medical | |
| FR2823345A1 (fr) | Procede d'amelioration de la qualite d'une image radiographique tridimensionnelle d'un objet et dispositif radiographique correspondant | |
| FR2927719A1 (fr) | Procede de traitement d'images obtenues par tomosynthese et dispositif associe | |
| WO2006131646A1 (fr) | Procédé et dispositif de reconstruction 3d d'un objet à partir de plusieurs images 2d |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PLFP | Fee payment |
Year of fee payment: 17 |
|
| PLFP | Fee payment |
Year of fee payment: 18 |
|
| PLFP | Fee payment |
Year of fee payment: 19 |
|
| PLFP | Fee payment |
Year of fee payment: 20 |