FR2862790A1 - Procede et dispositif pour l'observation d'un phenomene transitoire a l'aide d'un arceau vasculaire - Google Patents
Procede et dispositif pour l'observation d'un phenomene transitoire a l'aide d'un arceau vasculaire Download PDFInfo
- Publication number
- FR2862790A1 FR2862790A1 FR0313787A FR0313787A FR2862790A1 FR 2862790 A1 FR2862790 A1 FR 2862790A1 FR 0313787 A FR0313787 A FR 0313787A FR 0313787 A FR0313787 A FR 0313787A FR 2862790 A1 FR2862790 A1 FR 2862790A1
- Authority
- FR
- France
- Prior art keywords
- voxel
- function
- projections
- acquired
- given
- 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 40
- 238000003384 imaging method Methods 0.000 title claims abstract description 19
- 210000003462 vein Anatomy 0.000 title description 2
- 230000002123 temporal effect Effects 0.000 claims abstract description 19
- 238000005259 measurement Methods 0.000 claims description 21
- 239000002872 contrast media Substances 0.000 claims description 19
- 230000002792 vascular Effects 0.000 claims description 18
- 238000012905 input function Methods 0.000 claims description 13
- 238000002347 injection Methods 0.000 claims description 12
- 239000007924 injection Substances 0.000 claims description 12
- 230000003247 decreasing effect Effects 0.000 claims description 8
- 238000009792 diffusion process Methods 0.000 claims description 8
- 238000005316 response function Methods 0.000 claims description 8
- 238000010521 absorption reaction Methods 0.000 claims description 4
- 238000011161 development Methods 0.000 claims description 3
- 230000001052 transient effect Effects 0.000 description 17
- 238000002591 computed tomography Methods 0.000 description 15
- 238000001802 infusion Methods 0.000 description 5
- 238000002059 diagnostic imaging Methods 0.000 description 4
- 230000010412 perfusion Effects 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000302 ischemic effect Effects 0.000 description 2
- 238000002603 single-photon emission computed tomography Methods 0.000 description 2
- 230000001732 thrombotic effect Effects 0.000 description 2
- ZCYVEMRRCGMTRW-UHFFFAOYSA-N 7553-56-2 Chemical compound [I] ZCYVEMRRCGMTRW-UHFFFAOYSA-N 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 230000003042 antagnostic effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 210000004204 blood vessel Anatomy 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000001361 intraarterial administration Methods 0.000 description 1
- 238000010253 intravenous injection Methods 0.000 description 1
- 229910052740 iodine Inorganic materials 0.000 description 1
- 239000011630 iodine Substances 0.000 description 1
- 238000012633 nuclear imaging Methods 0.000 description 1
- 230000035479 physiological effects, processes and functions Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001225 therapeutic effect Effects 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- 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/44—Constructional features of apparatus for radiation diagnosis
- A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
- A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
- A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
-
- 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/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- 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/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
- A61B6/5229—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image
- A61B6/5235—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data combining image data of a patient, e.g. combining a functional image with an anatomical image combining images from the same or different ionising radiation imaging techniques, e.g. PET and CT
-
- 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- 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/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
-
- 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
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- General Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pulmonology (AREA)
- Algebra (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
L'invention propose un procédé et un dispositif pour l'obtention d'une représentation d'un phénomène physique qui varie dans le temps et dans l'espace à travers un objet organique, ledit objet étant représenté par un ensemble de voxels et une fonction paramétrique étant associée à chaque voxel pour modéliser l'évolution temporelle dudit phénomène en ce voxel, ledit procédé étant adapté à un dispositif d'imagerie par rayons X apte à élaborer des images 3D dudit objet à partir de projections 2D acquises, lors de la rotation dudit dispositif autour de l'objet, à des temps de mesure et selon des directions de projection donnés.On accède pour cela aux projections 2D acquises par ledit dispositif, auxdits temps de mesure et selon lesdites orientations donnés, généralement uniquement en vue de l'élaboration des images 3D et on détermine de manière exacte la fonction paramétrique associée à chaque voxel de l'objet à partir desdites projections 2D.
Description
En d'autres termes, l'obtention d'une représentation de phénomènes
physiologiques transitoires nécessite l'obtention de données quadridimensionnelles représentant à la fois les distributions spatiales et temporelles d'une substance à travers un objet organique.
On décrit généralement un phénomène dynamique à l'aide d'un modèle paramétrique.
On peut à ce titre considérer une fonction F(x, y, z, t) dans laquelle (x, y, z) sont les trois coordonnées spatiales et t est le temps. Cette fonction F décrit par exemple la densité d'un produit donné dans l'espace io tridimensionnel, en fonction du temps.
Comme mentionné ci-dessus, un exemple typique est la quantité d'agent de contraste injecté dans le flux artériel ou veineux d'un patient.
Cet agent de contraste est alors naturellement diffusé dans l'ensemble du corps du patient. L'observation, au cours du temps, de la concentration de cet agent de contraste dans les différents tissus du patient peut permettre d'établir un diagnostic relatif à l'état de ces tissus.
La concentration au point de coordonnées (x, y, z) du produit injecté peut typiquement être représentée par une exponentielle décroissante égale à 0 lorsque t est inférieur au moment to auquel l'injection est opérée, et de la forme A(t-to)exp (t-t0" lorsque t est plus grand que to. Cette fonction exponentielle a ainsi pour paramètres les quantités A, p et to.
Ce type de fonction exponentielle décroissante est observé pour une injection impulsionnelle de l'agent de contraste.
De manière plus générale, on observe la convolution entre la fonction de l'arrivée de l'agent de contraste (ci-après dénommée fonction d'entrée de contraste) audit point de coordonnées (x, y, z) et la fonction de réponse impulsionnelle (c'est-à-dire la fonction exponentielle décroissante dans le cas mentionné précédemment).
II n'existe aucune méthode directe permettant d'accéder à la fonction F(x, y, z, t).
Une approche indirecte est basée sur l'utilisation d'un système d'imagerie médicale par rayons X. Un tel système d'imagerie fournit effectivement les intégrales de la fonction F(x, y, z, t) selon un certain nombre de directions dans l'espace et l'utilisation d'algorithmes tomographiques bien connus de l'homme du métier permet de reconstruire indirectement la fonction F(x, y, z, t).
Une telle reconstruction est cependant possible seulement si les intégrales sont mesurées le long d'un ensemble de directions qui représente une révolution de 180 autour du corps imagé.
io Cet ensemble de mesures doit en outre être réalisé dans un intervalle temporel suffisamment restreint pour que les variations de la fonction F(x, y, z, t) au cours de l'acquisition des images médicales puissent être négligées.
La tomodensitométrie est une technique d'imagerie médicale assistée par ordinateur qui permet d'examiner un objet organique, plan par plan, en mesurant les différences de densité d'absorption des rayons X par les tissus biologiques.
Un tomodensitomètre (ou scanner CT selon l'expression anglo-saxonne Computed Tomography) réalise un balayage permettant l'acquisition d'une image 3D en une seconde environ (voire même en 0,5 seconde sur les équipements les plus récents).
La figure 1 illustre de manière schématique l'élaboration d'une image dite de perfusion à l'aide d'un tel scanner CT.
Le corps à examiner est radiographié selon différents angles (on parle de balayage tomographique autour du patient) de manière à réaliser un ensemble de projections appelé sinogramme (c'est-à-dire une coupe vue sous différentes incidences). Chaque colonne du sinogramme correspond ainsi à une projection selon un angle donné.
Et comme cela est mentionné ci-dessus, le temps nécessaire pour acquérir, à l'aide d'un scanner CT, toutes les projections d'un même balayage (réalisation d'un sinogramme par une rotation de l'appareil autour du patient) est faible devant la dynamique du phénomène transitoire observé.
Il est de la sorte possible de considérer que les projections d'un même balayage sont acquises à un même instant, sans qu'il n'y ait d'évolutions conséquentes du phénomène physique lors dudit balayage (le produit diffusant peu pendant ce laps de temps).
En réalisant une série temporelle de balayages, on acquiert un ensemble de sinogrammes successifs. Cet ensemble de sinogrammes réalisés à des temps de mesure successifs est référencé par la lettre a sur io la figure 1.
A partir de chaque sinogramme, la reconstruction tomographique calcule l'image d'une coupe (ou tranche) axiale de 1 à 10 millimètres d'épaisseur au travers de la structure interne du corps imagé.
On obtient donc (cf. référence b sur la figure 1) un ensemble de coupes CT en succession dans le temps, pour un angle d'incidence donné.
La fonction paramétrique est ensuite déterminée à partir de l'ensemble des reconstructions F(x, y, z, t;), où t; représente les temps successifs de balayage.
Une courbe de densité temporelle (cf. référence c) permet ainsi de déterminer les paramètres du modèle considéré afin de représenter l'évolution du phénomène transitoire dans le temps et dans l'espace.
Une image de perfusion est finalement obtenue (cf. référence d) qui représente l'évolution du phénomène transitoire pour la tranche considérée.
Finalement, lorsque le dispositif d'imagerie est un scanner CT, l'acquisition des données est très rapide, ce qui permet de négliger la dynamique du phénomène transitoire. Les aspects spatiaux et temporels peuvent alors être découplés.
D'un point de vue pratique, il existe un besoin pour réaliser ce type de d'observation d'un phénomène transitoire (détermination d'une fonction F(x, y, z, t)) sur un dispositif d'imagerie autre qu'un scanner CT (par exemple sur un arceau vasculaire).
L'arceau vasculaire est en effet le dispositif d'imagerie de choix pour connaître l'état de perfusion de certaines zones du cerveau et suivre ainsi un patient en situation d'urgence ischémique.
Si cela s'avère nécessaire, l'arceau vasculaire permet alors de réaliser une intervention invasive par introduction d'un cathéter à travers le réseau artériel jusqu'à la zone ischémique, et par injection in situ d'agent thrombotique.
L'efficacité de ces agents thrombotiques doit être contrôlée par un examen de type perfusion. Il est donc souhaitable de pouvoir réaliser cet io examen sur le dispositif même qui a permis de réaliser l'intervention, ce que permet l'arceau vasculaire.
Un arceau vasculaire peut, de manière similaire à un scanner CT, être entraîné en rotation autour du patient et acquérir un ensemble d'images en vue de permettre la construction d'images tridimensionnelles 3D.
Cependant, les performances des arceaux vasculaires aujourd'hui disponibles, ainsi que certaines raisons de sécurité, n'autorisent la rotation de l'arceau vasculaire qu'à une vitesse limitée. L'arceau tourne ainsi généralement autour de l'objet à une vitesse de l'ordre de 40 degrés par seconde.
En outre, ce type d'équipement peut seulement exécuter une rotation selon une couverture angulaire limitée, qui est typiquement d'environ 200 degrés.
On peut finalement estimer à environ 5 secondes le temps nécessaire pour une rotation de 180 degrés de l'arceau autour du patient.
La dynamique du phénomène temporel ne peut plus alors être négligée, du fait de sa non stationnarité lors d'une rotation de 180 degrés de l'arceau, et les algorithmes classiques de reconstruction tomographique ne peuvent plus être appliqués.
On notera en outre qu'il est nécessaire de disposer d'un nombre suffisant de mesures pour permettre l'estimation des paramètres de la fonction paramétrique F(x, y, z, t). Il faut pour cela réaliser un certain nombre de rotations de l'arceau. II a pour cela été proposé de procéder à des rotations alternatives dans le sens et dans le sens inverse des aiguilles d'une montre (on parle alors d'acquisition de type essui-glace ).
La multiplicité nécessaire des acquisitions essuie-glace ne fait alors qu'augmenter la durée d'acquisition, ce qui est antagonique avec la problématique de la dynamique du phénomène.
On note ici qu'il pourrait être envisagé de réduire la durée d'acquisition dé'l'ensemble des données en diminuant la couverture angulaire. Une telle solution entraînerait cependant une réduction importante de la qualité de la io reconstruction (données statistiquement non fiables, mauvaise résolution spatiale).
On précise également ici qu'il a été proposé un mode de rotation continu d'un arceau vasculaire autour d'un patient. Ce mode de rotation continu permet d'optimiser le temps d'acquisition en éliminant la perte de temps relative à l'inversion du sens de rotation lorsqu'un balayage essuieglace est réalisé. Malgré tout, le temps nécessaire pour réaliser une acquisition selon une couverture angulaire de 180 degrés est toujours relativement long devant la dynamique du phénomène. Et cette dynamique ne peut, y compris dans ce mode continu, être négligée.
Un objet de l'invention est de passer outre les limitations susmentionnées et de permettre l'observation d'un phénomène transitoire à l'aide d'un appareil d'imagerie médicale qui ne peut acquérir des images 3D d'un objet en un temps suffisamment court vis-à-vis de la dynamique du phénomène observé.
A cet effet, l'invention propose, selon un premier aspect, un procédé pour l'obtention d'une représentation d'un phénomène physique qui varie dans le temps et dans l'espace à travers un objet organique, ledit objet étant représenté par un ensemble de voxels et une fonction paramétrique étant associée à chaque voxel pour modéliser l'évolution temporelle dudit phénomène en ce voxel, ledit procédé étant adapté à un dispositif d'imagerie par rayons X apte à élaborer des images 3D dudit objet à partir de projections 2D acquises, lors de la rotation dudit dispositif autour de l'objet, à des temps de mesure et selon des directions de projection donnés, ledit procédé étant caractérisé en ce qu'il comporte les étapes consistant à : É accéder aux projections 2D acquises par ledit dispositif, auxdits temps de mesure et selon lesdites orientations donnés, en vue de l'élaboration des images 3D; É déterminer de manière exacte la fonction paramétrique associée à chaque voxel de l'objet à partir desdites projections 2D.
Des aspects préférés, mais non limitatifs, du procédé selon l'invention io sont les suivants: l'étape de détermination comporte une opération associant, pour chaque image représentant une projection 2D, l'intensité de chaque pixel composant ladite image à l'absorption des rayons X, au temps de mesure donné pour l'acquisition de ladite projection, par les voxels de l'objet situés sur une ligne de projection passant le pixel considéré et correspondant à l'orientation donnée de la mesure; - l'étape de détermination comporte la résolution de l'ensemble des équations du type Ip(x, y) = J J J a;lk(x,y)h;lk(pi,p2,...pn,tp), dans i j k lesquelles Ip(x, y) représente l'intensité d'un pixel d'une image représentant une projection 2D acquise à un temps de mesure donné tp, selon une direction de projection donnée, (i, j, k) représentent les coordonnées spatiales, hi J,k(p1,p2,...pn,t) représente la fonction de paramètres p1,p2,...pn associée au voxel de coordonnées (i, j, k), et a;, j,k est un coefficient dont la valeur est 0 ou 1 selon que le voxel de coordonnées (i, j, k) est sur la ligne de projection ou non; - l'étape de détermination peut être réalisée selon la méthode du gradient; - le procédé peut en outre comporter une étape classique de reconstruction tomographique des projections 2D acquises, de manière à identifier un point de départ pour l'application de l'étape de détermination.
le procédé peut en outre comporter, préalablement à l'étape de détermination, une étape de soustraction visant à éliminer des informations communes à des images; le phénomène observé peut être celui de la diffusion d'un agent de contraste injecté de manière intraartérielle ou intraveineuse dans ledit objet; la fonction paramétrique associée à chaque voxel de l'objet est la convolution entre la fonction de réponse impulsionnelle du tissu et la fonction d'entrée de contraste audit voxel; lo la fonction de réponse impulsionnelle est une fonction exponentielle décroissante; le procédé peut en outre comporter les étapes d'élaboration d'un modèle local du vaisseau dans lequel diffuse l'agent de contraste sous forme d'un cylindre, de modélisation du flux dudit agent de contraste à travers ledit cylindre et d'estimation de la fonction d'entrée de contraste à partir de ladite modélisation et des différentes projections disponibles; l'injection de l'agent de contraste étant réalisée de manière intraartérielle, la fonction d'entrée de contraste est une fonction de Dirac.
L'invention propose également, selon un second aspect, un dispositif d'obtention d'une représentation d'un phénomène physique qui varie dans le temps et dans l'espace à travers un objet organique, ledit objet étant représenté par un ensemble de voxels et une fonction paramétrique étant associée à chaque voxel pour modéliser l'évolution temporelle dudit phénomène en ce voxel, ledit dispositif comprenant un dispositif d'imagerie par rayons X aptes à élaborer des images 3D dudit objet à partir de projections 2D acquises lors de la rotation dudit dispositif autour de l'objet, à des temps de mesure et selon des directions de projection données, ledit dispositif étant caractérisé en ce qu'il comporte en outre: É des moyens d'accès auxdites projections 2D; et des moyens de détermination de manière exacte de la fonction paramétrique associée à chaque voxel de l'objet à partir desdites projections 2D.
Des aspects préférés, mais non limitatifs, du dispositif selon l'invention sont les suivants: - le dispositif peut comporter un module d'affichage pour afficher une représentation de l'évolution temporelle du phénomène physique à travers l'objet; - le dispositif d'imagerie est apte à observer la diffusion d'un agent de io contraste injecté dans ledit objet de manière intraveineuse ou intraartérielle; - le dispositif d'imagerie peut être un arceau vasculaire.
D'autres caractéristiques, buts et avantages de l'invention apparaîtront à la lecture de la description détaillée qui va suivre, et au regard des dessins annexés, donnés à titre d'exemples non limitatifs et sur lesquels: - la figure 1, déjà commentée, représente de manière schématique le procédé classique d'élaboration d'une image de perfusion à l'aide d'un scanner CT; - la figure 2 représente de manière schématique le procédé selon l'invention pour l'élaboration d'un volume de perfusion à l'aide d'un dispositif d'imagerie à rayons X, du type arceau vasculaire.
Comme on l'a vu précédemment, la vitesse d'acquisition d'un scanner CT est telle qu'il est possible de négliger l'évolution du phénomène transitoire observé lors d'un balayage (acquisition d'un sinogramme). Les aspects temporels et spatiaux du problème de l'observation dudit phénomène peuvent alors être découplés.
Dans le cadre de l'invention, la dynamique du phénomène ne peut être négligée. Les aspects temporels et spatiaux doivent alors être traités conjointement et le problème doit être considéré comme un problème quadridimensionnel. io
On considère ici l'injection (par exemple de manière intraveineuse ou intra-artérielle) d'un produit de contraste (par exemple à base d'iode) dans un objet organique (par exemple le corps d'un patient). Le phénomène transitoire à observer est alors celui de la diffusion dudit produit de contraste dans une zone d'intérêt de l'objet organique.
On considère plus particulièrement l'angiographie, c'est-à-dire l'injection d'un produit de contraste dans les vaisseaux sanguins (veines, artères) afin de les opacifier et de les visualiser dans un but diagnostique ou thérapeutique.
io Bien entendu, l'invention n'est pas limitée à l'observation de la diffusion d'un produit de contraste, mais s'étend à tout phénomène transitoire pouvant être observé à l'aide d'un dispositif d'imagerie apte à réaliser les fonctions ici décrites.
On considère une représentation discrète de l'espace selon laquelle la zone d'intérêt est constituée d'une pluralité de régions cubiques élémentaires. Une région cubique élémentaire est également appelée voxel (pixel tri-dimensionnel).
Dans le cadre de l'invention, on considère que l'activité de chaque voxel peut être modélisée par une fonction paramétrique.
Dans la suite de la description, les coordonnées spatiales i, j et k identifient un voxel particulier de la zone d'intérêt.
Une fonction paramétrique h;e(pl,p2,...pn,t) est ainsi associée au voxel de coordonnées spatiales i, j, k; pi, p2,....pn étant les paramètres de la fonction paramétrique.
De manière non limitative, cette fonction paramétrique peut être une fonction parabolique, une fonction similaire à l'exponentielle décroissante dont il a été fait mention précédemment, ou encore une fonction plus complexe, par exemple une bi-exponentielle décroissante.
L'invention met en oeuvre un dispositif d'imagerie médicale par rayons 30 X, par exemple un arceau vasculaire. Ledit dispositif d'imagerie est apte à tourner autour du patient en décrivant un arc, préférablement d'au moins 180 degrés.
Comme on l'a vu précédemment, un arceau vasculaire balaie un arc de 180 degrés autour du patient en 5 secondes environ.
Ledit arceau vasculaire est apte à élaborer une image 3D de la zone d'intérêt à partir d'une pluralité de projections 2D acquises à des positions données lors d'une rotation autour de l'objet le long de l'arc.
Lesdites projections 2D sont acquises à des temps de mesure et selon des directions de projection données.
Environ 10 projections 2D sont typiquement acquises en une seconde, et un balayage de 180 degrés autour du patient permet ainsi d'acquérir environ 50 projections 2D.
Chacune de ces images 2D est ainsi acquise en un temps suffisamment court devant la dynamique du phénomène transitoire observé et 5 secondes environ sont nécessaires pour acquérir un nombre suffisant de projections pour l'élaboration d'une image 3D de qualité.
Bien entendu, afin de disposer d'un nombre suffisant de mesures, l'invention peut mettre en oeuvre une pluralité de rotations de l'arceau autour de l'objet.
Une trentaine de secondes sont ainsi nécessaires pour réaliser 6 balayages autour du patient au cours desquels environ 300 images 2D du patient sont acquises.
Ce nombre important d'images 2D permet, comme cela sera plus particulièrement détaillé par la suite, d'obtenir un volume de perfusion 25 décrivant précisément l'évolution du phénomène transitoire.
Ladite pluralité de rotations de l'arceau autour du patient est préférablement réalisée en effectuant successivement des acquisitions de type essuie-glace .
Bien entendu, l'invention s'applique également lorsque l'arceau 30 vasculaire dispose d'un mode de rotation continu. L'échantillonnage temporelle des données versus leur répartition angulaire est dans ce cas avantageusement amélioré.
L'invention propose d'utiliser directement lesdites projections 2D qui sont généralement acquises uniquement en vue de l'élaboration des images 3D.
Comme on l'a vu précédemment, ces projections 2D sont chacune associées à un temps de mesure et à une direction de projection donnés.
En référence à la figure 2, on a représenté sous la référence e l'ensemble de ces projections, chacune étant associée à un couple (temps io t, direction de projection 8).
L'accès aux différentes mesures permet ainsi d'obtenir une succession d'images lp, chacune de ses images étant acquise à un temps de mesure donné tp, selon une orientation donnée.
L'intensité de chaque pixel composant l'image lp d'une des projections 2D correspond ainsi à l'absorption des rayons X, au temps de mesure tp, par les voxels de la zone d'intérêt situés sur une ligne de projection passant par le pixel considéré et correspondant à la direction de projection donnée de la mesure.
En d'autres termes, l'intensité de chaque pixel de ces images correspond à l'intégrale de la fonction paramétrique h;i,k le long d'une ligne dite de projection définie par le centre de projection et le pixel considéré.
On a ainsi la relation: Ip(x, y) = f f f a;j,k(x,y)h;i,k(p1,p2,...pn,tp), dans i j k laquelle a;i,k(x,y) est un coefficient dont la valeur est 0 ou 1 selon que le voxel de coordonnées (i,j,k) est sur la ligne de projection ou non.
On notera que le document CA 2 288 584 (ci-après CA) propose la description d'un phénomène dynamique par la mise en oeuvre d'un traitement d'images réalisées lors d'un balayage tomographique.
Plus précisément, il est proposé dans ce document CA un procédé et un appareil permettant de produire une image représentant le changement de radioactivité d'un objet en utilisant des équipements d'imagerie nucléaire standards (du type gammatomographe ou SPECT selon l'expression anglo-saxonne Single-Photon Emission Computed Tomography).
La dynamique du phénomène observé dans le cadre du document CA (diffusion d'un traceur radioactif) est particulièrement lente et autorise ainsi des temps d'acquisition particulièrement longs (une dizaine de secondes pour obtenir une image 2D, une vingtaine de minutes pour réaliser un balayage de 180 degrés autour du patient).
Comme on l'a vu précédemment, la présente invention s'intéresse quant à elle à des phénomènes transitoires dont la dynamique est nettement supérieure.
On rappelle que dans le cadre de l'invention, une projection 2D est acquise en un dixième de seconde, soit une différence d'un facteur 100 dans le temps d'acquisition, ce qui permet typiquement de réaliser 6 balayages et d'acquérir plus de 300 projections 2D en seulement une trentaine de secondes.
On notera également que la solution proposée dans le document CA décrit les contraintes temporelles sur les intensités reconstruites au fil du temps par des inégalités et n'impose aucune forme prédéterminée de fonction pour décrire l'évolution du phénomène dynamique en chaque voxel.
Et si l'utilisation de fonctions paramétriques est mentionnée dans ce document CA dans le cadre de l'imagerie 3D, elle n'est aucunement transposée dans le cadre de l'imagerie projective; le document CA la rejetant même en invoquant notamment son inefficacité en terme de calcul.
L'utilisation des projections 2D n'était donc ni enseignée, ni suggérée dans ce document CA afin de résoudre la problématique des phénomènes transitoires dans le cas d'une technique d'acquisition pour laquelle la solution proposée dans le document CA n'est a priori pas transposable.
Revenant à la description de l'invention, on détermine par la suite, afin de caractériser le phénomène transitoire, les paramètres de la fonction paramétrique associée à chacun des voxels de la zone d'intérêt.
Cette détermination est réalisée en résolvant les différentes relations intégrales du type Ip(x, y) = f f f a;e(x,y)h;lk(p1,p2,...pn,tp). i j k
En fonction du nombre de paramètres pris en compte et de la complexité de l'expression mise en jeux (linéaire, quadratique), différentes méthodes de résolution peuvent être envisagées.
Selon un mode de réalisation préféré de l'invention, la méthode du gradient est utilisée afin de déterminer les paramètres des fonctions paramétriques associées à chaque voxel.
Bien entendu, d'autres types de méthode d'identification des paramètres des fonctions paramétriques peuvent être envisagés, l'invention n'étant en aucun cas limitée au choix d'une méthode de résolution.
Finalement, l'invention ne nécessite pas la mise en oeuvre d'une étape intermédiaire visant à reconstruire des coupes CT à partir de sinogrammes, comme cela est classiquement réalisé avec les scanners CT.
Selon un mode de réalisation avantageux, des algorithmes tomographiques classiques, du type de ceux mentionnés précédemment notamment à l'examen de la figure 1, peuvent toutefois être mis en oeuvre.
Une telle résolution du problème tomographique par des méthodes régulières permet de procéder à une première reconstruction, sans inclure la dimension temporelle, et d'ajuster ainsi le modèle paramétrique choisi pour chaque voxel.
La méthode de résolution choisie (par exemple la méthode du gradient) peut ensuite être utilisée afin d'identifier l'optimum des paramètres du modèle préalablement ajusté.
Dans le cadre de ce mode de réalisation avantageux, les reconstructions tomographiques classiques sont ainsi utilisées lors d'une étape d'initialisation afin d'identifier un point de départ pour l'algorithme de la méthode de résolution choisie.
On notera que des méthodes de régularisation peuvent également être employées, et cela en utilisant en particulier le fait que le comportement de voxels voisins est quasi-similaire.
Finalement, et en référence à la figure 2, les aspects spatiaux et 5 temporels du problème de l'observation d'un phénomène transitoire sont couplés et traités conjointement dans le cadre de l'invention.
Les paramètres des fonctions paramétriques étant déterminés, on obtient finalement un volume dit de perfusion (cf. référence f) décrivant l'évolution spatiale et temporelle du phénomène transitoire à travers la zone d'intérêt de l'objet organique.
Bien entendu, le dispositif selon l'invention peut également comprendre un module d'affichage permettant d'afficher une représentation de l'évolution temporelle du phénomène physique à travers l'objet.
L'invention trouve en particulier application pour l'observation de la diffusion d'un produit de contraste, que ce produit ait été injecté de manière intraveineuse ou intra-artérielle dans le corps d'un patient.
Lorsque l'injection de l'agent de contraste est réalisée de manière intraveineuse, la fonction paramétrique h;i,k(p1,p2,...pn,t) est relativement complexe.
Il s'agit en effet de la convolution entre la fonction de réponse impulsionnelle du tissu et la fonction d'entrée artérielle.
De manière similaire, lorsque l'injection de l'agent de contraste est réalisée de manière intra-artérielle, la fonction paramétrique hiJ,k(p1, p2,É..Pn,t) correspond également à la convolution de la fonction de réponse impulsionnelle du tissu et la fonction d'entrée de contraste.
La détermination des paramètres est alors envisageable seulement si la fonction d'entrée artérielle peut être directement estimée.
Il est à cet effet possible d'élaborer un modèle local du vaisseau sous forme d'un cylindre et de modéliser le flux de produit de contraste à travers ce cylindre. La fonction d'entrée de contraste peut alors être estimée en utilisant cette modélisation et les différentes projections disponibles, le modèle cylindrique élaboré permettant effectivement de s'affranchir des variations dues au changement d'angles de projection.
Selon un mode de réalisation particulier de l'invention, lorsque l'injection de l'agent de contraste est réalisée de manière intraartérielle, la fonction d'entrée de contraste peut être considérée comme étant proche d'une fonction de Dirac. De telle sorte sa convolution avec la fonction de réponse impulsionnelle est égale à la fonction de réponse impulsionnelle elle-même. L'influence de la fonction d'entrée de contraste peut alors être négligée
10 dans les calculs. L'injection intra-artérielle permet ainsi un meilleur contraste et un bruit plus réduit que l'injection intraveineuse.
L'invention s'applique également en modes d'acquisition en soustraction, modes selon lesquels on élimine les informations communes à des images de manière à permettre notamment une lecture plus aisée des informations radiographiques.
Le mode soustrait est particulièrement utile en imagerie vasculaire puisqu'il permet la suppression des éléments anatomiques qui ne sont pas liés au phénomène étudié, ce qui rend la tache de reconstruction moins complexe.
On peut à titre d'exemple réaliser une soustraction d'une image de la zone d'intérêt réalisée, selon une orientation donnée, lorsque par exemple aucun produit de contraste ne diffuse dans ladite zone à une image acquise, selon la même orientation donnée, lors de la diffusion d'un produit de contraste dans ladite zone.
Dans le cadre de l'invention, la soustraction peut être réalisée sur les projections 2D elles-mêmes, avant le démarrage de la création du volume de perfusion.
L'invention procure ainsi un avantage supplémentaire vis-à-vis des techniques antérieures, pour lesquelles, du fait de limitations techniques, la reconstruction tomographique est réalisée sur des acquisitions non soustraites.
Claims (16)
1. Procédé pour l'obtention d'une représentation d'un phénomène physique qui varie dans le temps et dans l'espace à travers un objet organique, ledit objet étant représenté par un ensemble de voxels et une fonction paramétrique étant associée à chaque voxel pour modéliser l'évolution temporelle dudit phénomène en ce voxel, ledit procédé étant adapté à un dispositif d'imagerie par rayons X apte à élaborer des images 3D dudit objet à partir de projections 2D acquises, lors de la io rotation dudit dispositif autour de l'objet, à des temps de mesure et selon des directions de projection donnés, ledit procédé étant caractérisé en ce qu'il comporte les étapes consistant à: É accéder aux projections 2D acquises par ledit dispositif, auxdits 15 temps de mesure et selon lesdites orientations donnés, en vue de l'élaboration des images 3D; É déterminer de manière exacte la fonction paramétrique associée à chaque voxel de l'objet à partir desdites projections 2D.
2. Procédé selon la revendication 1, caractérisé en ce que l'étape de détermination comporte une opération visant à associer, pour chaque image représentant une projection 2D, l'intensité de chaque pixel composant ladite image à l'absorption des rayons X, au temps de mesure donné pour l'acquisition de ladite projection, par les voxels de l'objet situés sur une ligne de projection passant le pixel considéré et correspondant à l'orientation donnée de la mesure.
3. Procédé selon la revendication 2, caractérisé en ce que l'étape de détermination de la fonction paramétrique associée à un voxel comporte 30 la résolution de l'ensemble des équations du type Ip(x, y) = J J J ÀIj, k(x,Y)hi,j,k(p1,p2,...pn,tp), dans lesquelles Ip(x, y) représente i j k l'intensité d'un pixel d'une image représentant une projection 2D acquise à un temps de mesure donné tp, selon une direction de projection donnée, (i, j, k) représentent les coordonnées spatiales, h;,j,k(p1,p2,...pn,t) représente la fonction de paramètres p1,p2,...pn associée au voxel de coordonnées (i, j, k), et À;,j,k est un coefficient dont la valeur est 0 ou 1 selon que le voxel de coordonnées (i, j, k) est sur la ligne de projection ou non.
io
4. Procédé selon l'une des revendications précédentes, caractérisé en ce que l'étape de détermination est réalisée selon la méthode du gradient.
5. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'il comporte en outre une étape classique de reconstruction tomographique des projections 2D acquises, de manière à identifier un point de départ pour l'application de l'étape de détermination.
6. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'il comporte en outre, préalablement à l'étape de détermination, une étape de soustraction visant à soustraire de chaque image représentant une projection 2D acquise en présence du phénomène physique selon une orientation donnée, une image de l'objet acquise, en l'absence dudit phénomène, selon la même orientation donnée, de manière à éliminer les informations communes auxdites images.
7. Procédé selon l'une des revendications précédentes, dans lequel le phénomène observé est celui de la diffusion d'un agent de contraste, caractérisé en ce que la fonction paramétrique associée à chaque voxel de l'objet est la convolution entre la fonction de réponse impulsionnelle du tissu et la fonction d'entrée de contraste audit voxel.
8. Procédé selon la revendication précédente, caractérisé en ce que la fonction de réponse impulsionnelle est une fonction exponentielle décroissante, de paramètres A, p et t , égale à 0 avant le moment t auquel l'injection est opérée, et de la forme A(t-t0)exp"ct-t0)i' lorsque le temps t est plus grand que to.
9. Procédé selon l'une des deux revendications précédentes, caractérisé en ce qu'il comporte en outre les étapes d'élaboration d'un modèle local du vaisseau dans lequel diffuse l'agent de contraste sous forme d'un cylindre, de modélisation du flux dudit agent de contraste à travers ledit cylindre et d'estimation de la fonction d'entrée de contraste à partir de ladite modélisation et des différentes projections disponibles.
10. Procédé selon l'une des revendications 7 ou 8, dans lequel l'injection de l'agent de contraste est réalisée de manière intra- artérielle, ledit procédé étant caractérisé en ce que la fonction d'entrée de contraste est une fonction de Dirac.
11. Dispositif d'obtention d'une représentation d'un phénomène physique qui varie dans le temps et dans l'espace à travers un objet organique, ledit objet étant représenté par un ensemble de voxels et une fonction paramétrique étant associée à chaque voxel pour modéliser l'évolution temporelle dudit phénomène en ce voxel, ledit dispositif comprenant un dispositif d'imagerie par rayons X aptes à élaborer des images 3D dudit objet à partir de projections 2D acquises lors de la rotation dudit dispositif autour de l'objet, à des temps de mesure et selon des directions de projection données, ledit dispositif étant caractérisé en ce qu'il comporte en outre: É des moyens d'accès auxdites projections 2D; É et des moyens de détermination de manière exacte de la fonction paramétrique associée à chaque voxel de l'objet à partir desdites projections 2D.
12. Dispositif selon la revendication 11, caractérisé en ce qu'il comporte en outre un module d'affichage pour afficher une représentation de l'évolution temporelle du phénomène physique à travers l'objet.
io
13. Dispositif selon l'une des revendications 11 ou 12, caractérisé en ce que ledit dispositif d'imagerie est apte à observer la présence d'un agent de contraste injecté dans ledit objet de manière intraveineuse ou infra-artérielle.
14. Dispositif selon l'une des revendications 11 à 13, caractérisé en ce que le dispositif d'imagerie est un arceau vasculaire.
15. Dispositif selon la revendication 14, caractérisé en ce que l'arceau vasculaire est entraîné en rotation autour dudit objet le long d'un arc d'au moins 180 degrés alternativement dans un sens puis dans l'autre, de manière à réaliser une série de balayages de type essuie-glace dudit l'objet.
16. Dispositif selon la revendication 14, caractérisé en ce que l'arceau vasculaire est entraîné autour dudit objet selon un mode de rotation continu, de manière à réaliser une série de balayages complets dudit objet.
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0313787A FR2862790B1 (fr) | 2003-11-25 | 2003-11-25 | Procede et dispositif pour l'observation d'un phenomene transitoire a l'aide d'un arceau vasculaire |
DE102004055770A DE102004055770A1 (de) | 2003-11-25 | 2004-11-18 | Verfahren und Vorrichtung zur medizinischen Bildgebung |
JP2004340297A JP2005160102A (ja) | 2003-11-25 | 2004-11-25 | 医用イメージングの方法及び装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0313787A FR2862790B1 (fr) | 2003-11-25 | 2003-11-25 | Procede et dispositif pour l'observation d'un phenomene transitoire a l'aide d'un arceau vasculaire |
Publications (2)
Publication Number | Publication Date |
---|---|
FR2862790A1 true FR2862790A1 (fr) | 2005-05-27 |
FR2862790B1 FR2862790B1 (fr) | 2006-05-05 |
Family
ID=34531243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR0313787A Expired - Fee Related FR2862790B1 (fr) | 2003-11-25 | 2003-11-25 | Procede et dispositif pour l'observation d'un phenomene transitoire a l'aide d'un arceau vasculaire |
Country Status (3)
Country | Link |
---|---|
JP (1) | JP2005160102A (fr) |
DE (1) | DE102004055770A1 (fr) |
FR (1) | FR2862790B1 (fr) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101126725B (zh) * | 2007-09-24 | 2010-12-15 | 舒嘉 | 采用x射线容积摄影实现图像重建的方法 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102007045527A1 (de) | 2007-09-24 | 2009-04-16 | Siemens Ag | Verfahren und Vorrichtung zum Bestimmen eines Durchblutungsparameters von einem Gewebe sowie dessen Verwendung |
DE102009004184B3 (de) * | 2009-01-09 | 2010-07-29 | Siemens Aktiengesellschaft | Verfahren zur Bestimmung der arteriellen Inputfunktion für Perfusionsmessungen und C-Bogen-Röntgengerät |
-
2003
- 2003-11-25 FR FR0313787A patent/FR2862790B1/fr not_active Expired - Fee Related
-
2004
- 2004-11-18 DE DE102004055770A patent/DE102004055770A1/de not_active Withdrawn
- 2004-11-25 JP JP2004340297A patent/JP2005160102A/ja not_active Withdrawn
Non-Patent Citations (7)
Title |
---|
AKBUDAK E ET AL: "Contrast-agent phase effects: an experimental system for analysis of susceptibility, concentration, and bolus input function kinetics", MAGN. RESON. MED. (USA), MAGNETIC RESONANCE IN MEDICINE, DEC. 1997, WILLIAMS & WILKINS, USA, vol. 38, no. 6, December 1997 (1997-12-01), pages 990 - 1001, XP001173049, ISSN: 0740-3194 * |
FARNCOMBE T ET AL: "An evaluation of dynamic SPECT imaging methods", 1998, PISCATAWAY, NJ, USA, IEEE, USA, 1998, pages 1959 - 1963 vol., XP010338152, ISBN: 0-7803-5021-9 * |
FENG D ET AL: "A TECHNIQUE FOR EXTRACTING PHYSIOLOGICAL PARAMETERS AND THE REQUIRED INPUT FUNCTION SIMULTANEOUSLY FROM PET IMAGE MEASUREMENTS:THEORY AND SIMULATION STUDY", IEEE TRANSACTIONS ON INFORMATION TECHNOLOGY IN BIOMEDICINE, IEEE SERVICE CENTER, PISCATAWAY, NJ, US, vol. 1, no. 4, 1 December 1997 (1997-12-01), pages 243 - 254, XP000740698, ISSN: 1089-7771 * |
HEBBER E ET AL: "Direct estimation of dynamic parameters in SPECT tomography", IEEE TRANS. NUCL. SCI. (USA), IEEE TRANSACTIONS ON NUCLEAR SCIENCE, DEC. 1997, IEEE, USA, vol. 44, no. 6, pt.2, December 1977 (1977-12-01), pages 2425 - 2430, XP002287474, ISSN: 0018-9499 * |
LIMBER M A ET AL: "Direct reconstruction of functional parameters for dynamic SPECT", IEEE TRANS. NUCL. SCI. (USA), IEEE TRANSACTIONS ON NUCLEAR SCIENCE, AUG. 1995, USA, vol. 42, no. 4, pt.1, August 1995 (1995-08-01), pages 1249 - 1256, XP002287475, ISSN: 0018-9499 * |
OLLINGER J M ET AL: "A preliminary evaluation of the use of the EM algorithm for estimating parameters in dynamic tracer-studies", IEEE TRANS. NUCL. SCI. (USA), IEEE TRANSACTIONS ON NUCLEAR SCIENCE, FEB. 1985, USA, vol. NS-32, no. 1, February 1985 (1985-02-01), pages 848 - 854, XP001194453, ISSN: 0018-9499 * |
SPINELLI A E ET AL: "Arterial input function measurements using radiotracer and microbubbles: preliminary results", 2002, PISCATAWAY, NJ, USA, IEEE, USA, 2002, pages 2018 - 2022 vol., XP002287476, ISBN: 0-7803-7324-3 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101126725B (zh) * | 2007-09-24 | 2010-12-15 | 舒嘉 | 采用x射线容积摄影实现图像重建的方法 |
Also Published As
Publication number | Publication date |
---|---|
DE102004055770A1 (de) | 2005-12-15 |
FR2862790B1 (fr) | 2006-05-05 |
JP2005160102A (ja) | 2005-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
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 | |
FR2812741A1 (fr) | Procede et dispositif de reconstruction d'une image tridimensionnelle dynamique d'un objet parcouru par un produit de contraste | |
US6983182B2 (en) | Time resolved computed tomography angiography | |
EP0668570B1 (fr) | Procédé de traitement d'images pour déterminer localement le centre et la demi-largeur d'objets contrastant sur un fond et arrangement pour mettre en oeuvre ce procédé | |
WO1998010378A2 (fr) | Procede de reconstruction d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique | |
FR2862786A1 (fr) | Procede de recalage elastique d'une image applique a l'angiographie soustraite numerisee | |
FR2779853A1 (fr) | Procede de reconstruction d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique | |
JPH09253079A (ja) | X線断層撮影装置 | |
FR2641099A1 (fr) | ||
FR2800189A1 (fr) | Procede de reconstruction multi-resolutions d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique | |
FR2848007A1 (fr) | Procede et systeme pour accentuer une image de tomosynthese en utilisant un filtrage transversal | |
FR2884340A1 (fr) | Procede et dispositif de traitement d'images en angiographie soustraite | |
EP0927953A1 (fr) | Procédé de traitement d'une image multidimensionelle bruitée et appareil d'imagerie médicale mettant en oeuvre ce procédé | |
EP0752684A1 (fr) | Procédé de reconstruction d'images tridimensionnelles sur un objet mobile ou déformable | |
EP1639551B1 (fr) | Procede de reconstruction d'une image tomographique par une methode analytique comprenant une modelisation amelioree du mouvement de l'objet | |
FR2863086A1 (fr) | Procede de production d'une sequence d'images volumiques d'une zone d'un organe d'un etre vivant. | |
WO2016059329A1 (fr) | Système et procédé pour estimer une quantité d'intérêt d'un système dynamique artère/tissu/veine | |
Michel et al. | Weighted schemes applied to 3D-OSEM reconstruction in PET | |
FR2862790A1 (fr) | Procede et dispositif pour l'observation d'un phenomene transitoire a l'aide d'un arceau vasculaire | |
FR2965651A1 (fr) | Reconstruction tomographique d'un objet en mouvement | |
EP0814430A1 (fr) | Procédé de production d'une séquence restaurée d'images d'un objet en mouvement à partir de mesures bruitées | |
EP3509483B1 (fr) | Systeme et procede pour reconstruire un signal physiologique d'un systeme dynamique artere/tissu/veine d'un organe dans un espace surfacique | |
FR2733841A1 (fr) | Procede de realisation de la cartographie d'emission d'un corps corrigee de l'attenuation par ce corps | |
FR3023155A1 (fr) | Systeme et procede d'imagerie tridimensionnelle en profondeur | |
EP1521214A2 (fr) | Procédé de reconstruction d'image avec acquisition de donnees selon plusieurs modes d'imagerie et synchronisation des series d'images |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
ST | Notification of lapse |
Effective date: 20090731 |