FR2949590A1 - Procede de reconstruction tridimensionnelle d'un objet a partir d'une seule vue - Google Patents
Procede de reconstruction tridimensionnelle d'un objet a partir d'une seule vue Download PDFInfo
- Publication number
- FR2949590A1 FR2949590A1 FR0955984A FR0955984A FR2949590A1 FR 2949590 A1 FR2949590 A1 FR 2949590A1 FR 0955984 A FR0955984 A FR 0955984A FR 0955984 A FR0955984 A FR 0955984A FR 2949590 A1 FR2949590 A1 FR 2949590A1
- Authority
- FR
- France
- Prior art keywords
- elementary
- elementary object
- source
- projection image
- recording means
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; 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
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- 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/12—Arrangements for detecting or locating foreign bodies
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30021—Catheter; Guide wire
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/436—Limited angle
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Le procédé évalue l'orientation dans l'espace d'un objet élémentaire (10) de forme sensiblement cylindrique de révolution dans un dispositif de radiographie comportant une source (30) de rayons X et des moyens d'enregistrement (20) situés en regard de la source, avec des étapes de : a- détermination sur l'image de projection d'une quantité A de rayons X absorbés par l'objet élémentaire; et b- évaluation d'un sinus d'un angle θ entre un axe longitudinal de l'objet élémentaire et un plan de l'image de projection en fonction de la quantité A de rayons X absorbés. Le procédé reconstruit alors l'objet formé d'objets élémentaires avec des étapes de : a- détermination d'une longueur de l'objet élémentaire S ; et b- détermination des coordonnées relatives dans l'espace d'une P des extrémités (P , P ) de l'objet élémentaire S par rapport à l'autre P extrémité.
Description
DOMAINE TECHNIQUE GENERAL L'invention concerne l'évaluation de l'orientation et d'une dimension d'un objet élémentaire de forme sensiblement cylindrique de révolution ainsi que la reconstruction tridimensionnelle d'un objet formé d'une série de ces objets élémentaires ou assimilable par approximation à une série de ces objets élémentaires, le tout à partir d'une image de projection bidimensionnel dans un dispositif de radiographie.
ETAT DE LA TECHNIQUE Dans certaines procédures interventionnelles sous rayons X, le praticien est amené à insérer des cathéters et autres fils-guide à l'intérieur du corps du patient qu'il est en train de traiter. Les cathéters peuvent servir pour l'injection d'un liquide de contraste, comme de l'iode, et les fil-guides peuvent être utilisés pour la mise en place d'instruments thérapeutiques comme des ballons, des endoprothèses vasculaires, etc... Ces cathéters et fil-guides ont en commun d'être des structures curvilinéaires avec une section sensiblement circulaire.
Dans le cadre de ces procédures interventionnelles, et particulièrement en neurochirurgie, le réseau vasculaire du patient traité est très compliqué et une projection bidimensionnelle fournit très souvent insuffisamment d'informations. De ce fait, il est alors nécessaire de construire des modèles tridimensionnels de ce réseau vasculaire et de fournir une aide pour la localisation des instruments dans ces modèles tridimensionnels. Ce besoin émerge aussi en cardiologie interventionnelle où le rôle des modèles tridimensionnels est croissant sous l'impulsion du développement des scanners cardiaques.
Dans ce contexte, il est admis qu'une seule vue bidimensionnelle de 30 l'instrument interventionnel n'est pas suffisante pour en construire une représentation tridimensionnelle. Certaines techniques ont été développées à partir d'acquisitions bi-planaires, qui permettent d'utiliser simultanément deux projections bidimensionnelles sous des angles différents. Toutefois, de telles projections ne sont pas suffisantes pour reconstruire en trois dimensions sans aucune ambigüité une structure curvilinéaire, comme un fil-guide ou un cathéter. De plus, cela nécessite une dose de rayonnement supérieur à l'acquisition d'une seule prise de vue, ce qui n'est jamais souhaitable. De telles techniques sont décrites dans les articles de Merle, A.B.et al., intitulé 3D reconstruction of the deformable coronary tree skeleton from twoX-ray angiographic views dans Computers in Cardiology , du 13-16 Sep 1998, pages 757-760, et de Baert, S.A.M. et al., intitulé Three-dimensional guide-wire reconstruction from biplane image sequences for integrated display in 3-D vasculature , dans Medical Imaging, IEEE Transactions on, Volume 22, Issue 10, Oct. 2003, pages 1252-1258. Dans le cadre de la cardiologie interventionnelle, la reconstruction tridimensionnelle avec les méthodes traditionnelles demandant d'acquérir des images sous différentes angulations est généralement limitée par le bougé entre les prises de vues induit par le mouvement cardiaque et respiratoire du patient.
PRESENTATION DE L'INVENTION Un but de l'invention est de fournir un procédé d'évaluation de l'orientation dans l'espace d'un objet élémentaire de forme sensiblement cylindrique de révolution qui soit précis tout en réduisant l'irradiation du patient.
Selon un premier aspect, l'invention concerne un procédé d'évaluation de l'orientation dans l'espace d'un objet élémentaire de forme sensiblement cylindrique de révolution dans un dispositif de radiographie comportant une source de rayons X et des moyens d'enregistrement situés en regard de la source et agencés de sorte à enregistrer une image de projection de l'objet élémentaire alors situé entre la source et les moyens d'enregistrement, le procédé comportant des étapes de : a- détermination sur l'image de projection d'une quantité A de rayons X absorbés par l'objet élémentaire; et b- évaluation du sinus d'un angle e entre un axe longitudinal de l'objet élémentaire et un plan de l'image de projection en fonction de la quantité A de rayons X absorbés.
Ainsi, l'utilisation des propriétés d'absorption en rayons X de l'objet élémentaire permet de déterminer précisément l'orientation dans l'espace de l'objet à partir d'une seule image de projection bidimensionnel de l'objet ce qui évite, de plus, d'irradier inutilement le patient.
Avantageusement, mais facultativement, ce procédé présente au moins l'une des caractéristiques supplémentaires suivantes : - dans l'étape b, sin(9) _ ln A où r est un rayon d'une section de l'objet élémentaire et p une absorbance linéaire radiologique de l'objet élémentaire ; et, - dans l'étape b, sin(9) = ln(A) , où K est une constante d'absorbance 20 dépendante de l'objet élémentaire.
Selon un deuxième aspect, l'invention concerne un procédé d'évaluation d'une dimension dréel dans l'espace d'un objet élémentaire de forme sensiblement cylindrique dans un dispositif de radiographie comportant 25 une source de rayons X et des moyens d'enregistrement situés en regard de la source et agencés de sorte à enregistrer une image de projection de l'objet élémentaire alors situé entre la source et les moyens d'enregistrement, le procédé comportant des étapes de : a- évaluation d'une orientation dans l'espace de l'objet élémentaire 30 par le procédé selon le premier aspect; b- détermination sur l'image de projection d'une dimension projetée dobs; et c- détermination de la dimension dréel à partir de l'orientation évaluée et de la dimension projetée dobs déterminée. Ainsi, un tel procédé permet de déterminer la longueur d'une portion de l'objet en s'affranchissant des effets projectifs inhérents aux rayons X utilisés dans le dispositif de radiographie.
10 Avantageusement, mais facultativement, ce procédé présente la caractéristique supplémentaire suivante : - dans l'étape c, la dimension dréel est évaluée par : dréel=dobs/sin(e).
15 Selon un troisième aspect, l'invention concerne un procédé de reconstruction tridimensionnel d'un objet formé d'une série d'objets élémentaires de forme sensiblement cylindrique de révolution dans un dispositif de radiographie comportant une source de rayons X et des moyens d'enregistrement situés en regard de la source et agencés de sorte 20 à enregistrer une image de projection de l'objet alors situé entre la source et les moyens d'enregistrement, le procédé de reconstruction tridimensionnel comportant, pour chaque objet élémentaire S;=[P;_1,P;] de l'objet, des étapes de : a- détermination d'une orientation dans l'espace de l'objet 25 élémentaire Si par le procédé selon le premier aspect; et b- détermination des coordonnées relatives dans l'espace d'une P; des extrémités (P;_l, P;) de l'objet élémentaire Si par rapport à l'autre P;_l des extrémités (P;_l, P;) de l'objet élémentaire Si.
30 Avantageusement, mais facultativement, ce procédé présente l'une des caractéristiques supplémentaires suivantes :5 - entre l'étape a et l'étape b, le procédé comporte une étape de : al- détermination d'une longueur de l'objet élémentaire Si par le procédé selon le deuxième aspect ; - lors de l'étape a, une valeur de l'angle e1 est choisie parmi les valeurs e1 et u-e1 de sorte à ce que e1 soit le plus proche possible de 6;_1 ; - une fois traité l'ensemble des objet élémentaire S;=[P;_1,P;] de l'objet, le procédé détermine la longueur totale de l'objet reconstruit ; et, - la longueur totale de l'objet reconstruit est la somme des longueurs des objets élémentaires Si déterminées à l'étape b. Selon un quatrième aspect, l'invention concerne un dispositif de radiographie comportant une source de rayons x et des moyens d'enregistrement situés en regard de la source et agencés de sorte à enregistrer une image de projection d'un objet alors situé entre la source et les moyens d'enregistrement, et des moyens de mise en oeuvre d'au moins un des procédés précédents.
PRESENTATION DES FIGURES D'autres caractéristiques et avantages de l'invention apparaîtront dans la description ci-après d'un mode de réalisation de l'invention. Aux dessins annexés : la figure 1 est un schéma tridimensionnel de principe montrant l'absorption de rayons X d'un objet élémentaire, La figure Ibis est un schéma tridimensionnel de principe du calcul, la figure 2 et 2bis sont des schémas de principe du procédé d'évaluation de l'orientation, d'une part, et, d'autre part, d'une dimension d'un objet élémentaire, les figures 3 et 4 sont des schémas tridimensionnels de principe du procédé de reconstruction tridimensionnelle d'un fil-guide, et - la figure 5 illustre schématiquement un dispositif radiographique médical pour la mise en oeuvre du procédé de l'invention
DESCRIPTION DETAILLEE DE L'INVENTION Sur la figure 5, est illustré, schématiquement, un dispositif de radiographie ou d'imagerie médicale 100 pour l'acquisition d'images bidimensionnelles de projection pour l'évaluation de l'orientation d'un objet élémentaire de forme sensiblement cylindrique de révolution et pour la reconstruction tridimensionnelle d'un objet formé d'une série de ces objets élémentaires. Un tel dispositif peut être un dispositif de radiographie pour la détection et le repérage d'objets introduits au sein d'un patient lors d'interventions chirurgicales. Le dispositif de radiographie 100 comprend un système d'acquisition 15 d'images 103, un système de traitement d'images 105 et un système d'affichage 104. Le système d'acquisition 103 permet d'acquérir une pluralité de projections bidimensionnelles d'une région d'intérêt - d'un organe - d'un patient. Le système d'acquisition 103 est notamment constitué par un 20 détecteur 20 situé en regard d'une source de rayons X 30. Le détecteur est par exemple une caméra numérique. Le système d'acquisition est par exemple un système d'acquisition par rayons X, ce dernier comprenant tout moyen connu permettant l'émission de rayons X sur l'objet 102 et l'acquisition d'images résultantes. 25 Le système d'affichage 104 peut être intégré dans le système d'acquisition d'images 103 ou le système de traitement d'images 105, ou être séparé du système d'acquisition 103 et du système de traitement 105. Le système d'affichage 104 est par exemple un écran d'ordinateur, un moniteur, un écran plat, un écran plasma ou tout type de dispositif d'affichage connu 30 du commerce. Le système d'affichage 104 permet à un praticien de contrôler la reconstruction et/ou l'affichage des images bidimensionnelles acquises. Le système de traitement 105 est adapté à la mise en oeuvre de procédé de traitement (par exemple de reconstruction d'une image tridimensionnelle à partir d'images bidimensionnelles). Le système de traitement 105 peut être intégré dans le système d'acquisition 103 d'images ou être séparé du système d'acquisition 103 d'images. Le système de traitement 105 est par exemple un/des ordinateur(s), un/des processeur(s), un/des microcontrôleur(s), un/des micro-ordinateur(s), un/des automate(s) programmable(s), un/des circuit(s) intégré(s) spécifique(s) d'application, d'autres circuits programmables, ou d'autres dispositifs qui incluent un ordinateur tel qu'une station de travail. Le système de traitement 105 est couplé à des moyens mémoires 106 qui peuvent être intégrés ou séparés du système de traitement 105. Ces moyens peuvent être formés par un disque dur ou tout autre moyen de stockage amovible (un cd-rom, une disquette etc.). Ces moyens mémoires peuvent servir à stocker une image tridimensionnelle de la zone de l'organe visualisée comme une image bidimensionnelle acquise ou traitée. Il peut s'agir d'une mémoire ROM/RAM du système de traitement 105, un CD-ROM, une clé USB, une mémoire d'un serveur central. Le système de traitement 105 peut comprendre un dispositif de lecture (non représenté) par exemple un lecteur de disquettes ou un lecteur de CD-ROM, pour lire les instructions du procédé de traitement (qui va être décrit dans la suite) d'un support d'instructions (non montré), comme une disquette ou un CD-ROM. En variante, le système de traitement 105 exécute les instructions du procédé de traitement (qui va être décrit dans la suite) stockées dans des micrologiciels (non représentés).
En référence aux figures 1, 1 bis, 2 et 2bis, nous allons décrire maintenant un procédé d'évaluation de l'orientation d'un objet élémentaire Si (10) qui présente une forme homogène sensiblement cylindrique de révolution de longueur dréel, d'axe (11) et dont une section circulaire présente un rayon r. Pour l'illustration de notre propos, cet objet élémentaire est, ici, une partie d'un fil-guide. Il est bien entendu que ce qui suit est applicable à un cathéter ainsi qu'à tout autre instrument interventionnel pouvant être modélisé en une série d'objets élémentaires Si (10) du type précité.
De manière générale, ce qui va suivre est applicable à tout objet dont l'absorbance est fonction de l'orientation.
Le procédé est destiné à être mis en oeuvre dans un dispositif de radiographie 100 comportant la source 30 de rayons X ainsi que les moyens d'enregistrement 20 situés en regard de cette source de rayons X. Les moyens d'enregistrements sont agencés de sorte à pouvoir prendre une image de projection bidimensionnel d'un objet qui est alors situé entre ladite source de rayons X et lesdits moyens d'enregistrement. Ici, les moyens d'enregistrement 20 sont sensiblement plans et comprennent un détecteur numérique.
Il est connu de la personne de l'Art des procédés qui permettent de déterminer, sur une image de projection bidimensionnelle, l'axe ou la ligne médiane d'un fil-guide. Pour chaque point de projection du fil-guide situé sur cet axe ou ligne médiane, une absorbance radiologique est calculable. Il est entendu qu'une absorbance radiologique représente l'absorption en rayons X du fil-guide au niveau du point considéré. Il est possible d'estimer cette absorbance radiologique en formant le ratio entre une valeur de quantité de photons X au point (Pi)considéré du fil-guide et une valeur de quantité de photons X au niveau d'un arrière plan (AP) proche du point considéré du fil-guide, comme cela est illustré en figure 1 bis.
L'objet élémentaire Si (10) de rayon r présente une absorbance linéaire radiologique p. Cela veut dire que si N photons X (1) traversent une épaisseur l de matière formant un tel cylindre de révolution, seuls
ressortent photons X (2). L'absorbance du cylindre de révolution N.e-P.l est alors A= N = e-'' Toutefois, pour l'objet élémentaire Si, l'épaisseur l de matière traversée par un rayonnement X incident (1) dépend seulement de l'angle e entre l'axe (11) de l'objet élémentaire Si et le rayonnement X incident (1). La relation est donnée par l'équation : sin(e)=2.r/l. De ce fait, l'absorbance A de l'objet élémentaire dépend seulement de l'angle e et inversement. D'où : ln(A) La détermination de la valeur de l'angle e ne requiert que de connaître précisément une valeur du rayon r de la section de l'objet élémentaire Si ainsi qu'une valeur de l'absorbance radiologique linéaire p dudit objet élémentaire Si, valeur dépendant du matériau formant ledit objet élémentaire et du type de rayons X émis par la source 30 vers les moyens d'enregistrement 20. Ces paramètres sont facilement déterminés lors d'une acquisition et à partir de données constructeur aussi bien du fil-guide que du dispositif de radiographie dans lequel le procédé est mis en oeuvre. Il est à noter que, dans cette équation, si l'équivalence orientation - absorbance est connue en un point, elle est utilisable pour estimer le produit pr, de là une constante d'absorbance K=2pr, et ainsi en déduire la relation pour n'importe quelle orientation pour le même objet. Ceci permet de calibrer la méthode à partir d'une image simple acquise dans des conditions particulières telle que l'orientation de l'objet soit connue en au moins un point.
Toutefois, l'équation ci-dessus ne donnant qu'une valeur du sinus de l'angle 0, la valeur de l'angle e en lui-même peut être égale à e ou à u-e. Cette ambiguïté ne permet pas de savoir si l'objet élémentaire Si (10) à 2.r sin(9) _ ù,u. modéliser vient vers les moyens d'enregistrements (20)(voir figure 2) ou s'éloigne desdits moyens d'enregistrement (20) (voir figure 2bis). La limite de cette méthode est alors atteinte.
Une fois la valeur de l'angle e déterminée, et en particulier la valeur de sin(e), selon le procédé précédemment décrit, il est alors possible de déterminer la longueur dréel de l'objet élémentaire Si (10) à partir d'une dimension dobs mesurée sur l'image de projection dudit objet élémentaire. En effet : dréel=dobs/sin(e)
En référence aux figures 3 et 4, nous allons maintenant décrire un procédé de reconstruction tridimensionnel d'un objet présentant une structure curvilinéaire tel un fil-guide ou un cathéter. Une solution pour réduire l'ambiguïté sur e sera abordée et décrite.
Connaissant la forme de la projection du fil-guide par le repérage d'une ligne médiane (22) de ladite projection dans l'image de projection, cette forme est découpée en une série de tronçons élémentaires sensiblement rectilignes. Chacun de ces tronçons élémentaires correspond à la projection sur les moyens d'enregistrement (20) d'un axe médian (11) d'un objet élémentaire Si (10) du fil-guide à reconstruire en trois dimensions. En pratique, ici, la taille de chacun des tronçons élémentaires peut correspondre à la taille d'un pixel de l'image de projection appartenant à la ligne médiane (22) de la projection du fil-guide. Bien entendu, tout autre découpage en tronçons est possible et utilisable. Pour simplifier la description, il est supposé que l'image de projection est obtenue lors d'une projection dite parallèle : les rayons X (1) sont parallèles les uns par rapport aux autres et perpendiculaire à un plan (X,Y) contenant l'image de projection. Ainsi, les rayons X (1) sont selon un axe Z.
En projection parallèle, des coordonnées en X et Y sont invariables, et, selon notre hypothèse, toute translation selon l'axe Z produit la même image. De ce fait, il est impossible de connaître la position du point considéré selon l'axe Z, c'est-à-dire long d'une ligne de projection. De manière analogue, en projection conique, il n'est pas possible de connaître la position d'un point le long de la ligne de projection.
Soient p;(x;,y;)un pixel de la projection du fil-guide sur l'image de projection et P;(X;,Y;,Z;) un point de la reconstruction tridimensionnelle du fil-guide se projetant sur l'image de projection au pixel pi. Le procédé de reconstruction commence en un point po de la projection du fil-guide qui sert alors de pixel de départ. Soit alors po(xo=O,yo=O) l'origine du plan (X,Y). Soient Xo=xo=O, Yo=yo=O et Zo=O les coordonnées du point Po de l'extrémité de la reconstruction tridimensionnelle du fil-guide.
Le procédé passe au pixel suivant de la projection (22) du fil-guide. Dans le cas présent, le pixel suivant est le pixel pi et correspond à la projection du point PI que le procédé va déterminer. Le procédé évalue alors l'absorbance de l'objet élémentaire Si=[Po,PI] à partir du tronçon [po,pi] qui est sa projection sur l'image de projection. Selon l'image et le modèle de projection, le point PI est sur une ligne droite de l'espace qui est parallèle à l'axe Z et qui passe par le pixel pl de l'image. De ce fait, Xi=xi et Yi=y,. Le procédé alors détermine, d'une part la valeur de sin(e) (ainsi que la valeur de la longueur de l'objet élémentaire SI, dréell, en utilisant les procédés précédemment décrits s'il est nécessaire de calculer ultérieurement une longueur de l'objet reconstruit). Comme nous l'avons vu, il existe alors deux possibilités pour le positionnement de PI: l'angle et a soit pour valeur et soit u-e1. Il n'est pas possible de lever cette ambiguïté la méthode de calcul du sinus de l'angle et précédemment décrite. Le procédé va alors avoir dès le début deux reconstructions possibles pour le fil-guide selon la valeur choisie pour el. Le procédé calcule alors de la même manière le point P2 qui se projette alors sur le pixel P2 suivant de la projection (22) du fil-guide. Le procédé détermine l'absorbance de l'objet élémentaire S2=[P1,P2], les valeurs de sin(02) et éventuellement dréel2• Maintenant, le procédé va lever l'ambiguïté sur la valeur de 02. Pour cela, le procédé va considérer que les instruments interventionnels tels que des fil-guides ou des cathéters sont lisses ou réguliers, c'est-à-dire qu'ils ne peuvent pas changer brusquement de direction. Mathématiquement parlant, un angle entre deux objets élémentaires Si_l et Si doit être minimal. Donc, le procédé choisit la valeur de 02 de telle sorte que 02 soit le plus proche possible de e1. Connaissant alors l'angle 02i le procédé détermine alors les coordonnées de P2. Ensuite de manière itérative, le procédé va déterminer chaque point Pi de la reconstruction tridimensionnelle du fil-guide se projetant au pixel pi de la projection (22) sur l'image de projection : le procédé détermine l'absorbance de l'objet élémentaire Si=[Pi_1,Pi], les valeurs de sin(ei) et éventuellement de dréeu. Le choix de la valeur de 0i est tel que 0i et 0i_1 soient le plus proches possibles. De là, le procédé détermine les coordonnées du point Pi.
Une ambigüité peut survenir pour le calcul de 0i lorsque l'objet élémentaire précédent Si_1 est parallèle à l'image de projection (donc des moyens d'enregistrement). Dans ce cas il n'est pas possible de choisir entre u-0i et 0i par la méthode précédente. A chacune de ces ambigüités particulières, il y a donc deux possibilités de reconstruction tridimensionnelle qui existe. Dans la pratique, ces ambigüités sont peu nombreuses et il y a donc un nombre fini et restreint de reconstructions tridimensionnelles possibles à partir d'une seule vue. Soit npara le nombre de points où l'objet est parallèle aux moyens d'acquisition, il y a 2'1ù'1 reconstructions tridimensionnelles possibles. De ce fait, le procédé de reconstruction fournit à l'utilisateur un nombre réduit de reconstructions tridimensionnelles. L'utilisateur sélectionne la bonne reconstruction selon son expérience et ses connaissances de la situation. Par exemple, si l'objectif de la reconstruction est de situer le fil-guide dans un volume tridimensionnel représentant les vaisseaux, certaines reconstructions peuvent être éliminées en se basant sur la forme des vaisseaux.
Il est aussi possible de déterminer la longueur totale du fil-guide, et donc de connaitre la longueur réelle de fil-guide introduit dans le patient lors de l'intervention chirurgicale. Ceci peut être réalisé sans faire de reconstruction tridimensionnelle, en calculant seulement les dréeli.• Pour cela, le procédé intègre l'équation dréel=dobs/sin(e) le long du fil-guide (ou s représente l'abscisse curviligne le long du fil-guide dans l'image de projection): d Longueur/ii . s guide û Pro;eetiousm(O (s En pratique, le procédé peut additionner l'ensemble des valeurs calculées : Longueurfil-guide û ~juréel, Cette information permet un suivi plus précis de l'instrument interventionnel, par la présence d'un invariant. Bien entendu, il est possible d'apporter à l'invention de nombreuses modifications sans sortir du cadre de celle-ci. dréeli20
Claims (10)
- REVENDICATIONS1. Procédé de reconstruction tridimensionnelle d'un objet élémentaire (10) de forme sensiblement cylindrique de révolution dans un dispositif de radiographie comportant une source de rayons X (30) et des moyens d'enregistrement (20) situés en regard de la source et agencés de sorte à enregistrer une image de projection de l'objet élémentaire alors situé entre la source et les moyens d'enregistrement, le procédé comportant une étape d'évaluation de l'orientation dans l'espace de l'objet élémentaire comprenant des sous-étapes de : a- détermination sur l'image de projection d'une quantité A de rayons X absorbés par l'objet élémentaire; et b- évaluation du sinus d'un angle 0 entre un axe longitudinal de l'objet élémentaire et un plan de l'image de projection en fonction de la 15 quantité A de rayons X absorbés.
- 2. Procédé selon la revendication 1, caractérisé en ce que, dans la sous-étape b, sin(0) _ -2,u ln(A) ou r est un rayon d'une section de l'objet élémentaire et p une absorbance linéaire radiologique de l'objet 20 élémentaire.
- 3. Procédé selon la revendication 1, caractérisé en ce que, dans la sous-étape b, sin(0) = K ln(A)' où K est une constante d'absorbance dépendante de l'objet élémentaire. 25
- 4. Procédé selon l'une des revendications 1 à 3, caractérisé en ce qu'il comporte une étape supplémentaire d'évaluation d'une dimension dréel dans l'espace de l'objet élémentaire (10) comprenant des sous-étapes de : c- détermination sur l'image de projection d'une dimension projetée 30 dobs; et d- détermination de la dimension dréet à partir de l'orientation évaluée et de la dimension projetée dobs déterminée.
- 5. Procédé selon la revendication 4, caractérisé en ce que, dans la sous- étape d, la dimension dréet est évaluée par : dréel=lobs/sin(8).
- 6. Procédé de reconstruction tridimensionnelle d'un objet formé d'une série d'objets élémentaires (10) de forme sensiblement cylindrique de révolution dans un dispositif de radiographie comportant une source de rayons X (30) et des moyens d'enregistrement (20) situés en regard de la source et agencés de sorte à enregistrer une image de projection de l'objet alors situé entre la source et tes moyens d'enregistrement, le procédé de reconstruction tridimensionnel comportant, pour chaque objet élémentaire S;=[P;_1,P] de l'objet, des étapes de : a- détermination d'une reconstruction tridimensionnelle de l'objet élémentaire Si par un procédé selon l'une des revendications 1 à 5 ; et b- détermination des coordonnées relatives dans l'espace d'une (P;) des extrémités (P;_1, P;) de l'objet élémentaire Si par rapport à l'autre (Pm) des extrémités (Pm, P;) de l'objet élémentaire Si.
- 7. Procédé selon la revendication 6, caractérisé en ce que, lors de l'étape a, une valeur de l'angle 0; est choisie parmi les valeurs a; et u-01 de sorte à ce que 6; soit le plus proche possible de 6;_i.
- 8. Procédé selon la revendication 7, caractérisé en ce que, une fois traité l'ensemble des objet élémentaire S;=[P;.1,P;] de l'objet, le procédé détermine la longueur totale de l'objet reconstruit.
- 9. Procédé selon la revendication 8, caractérisé en ce que la longueur 30 totale de l'objet reconstruit est la somme des longueurs des objets élémentaires Si déterminées à l'étape b.
- 10. Dispositif de radiographie comportant une source de rayons x (30) et des moyens d'enregistrement (20) situés en regard de la source et agencés de sorte à enregistrer une image de projection d'un objet (10) alors situé entre la source et les moyens d'enregistrement, caractérisé en ce qu'il comporte en outre des moyens de mise en oeuvre d'un procédé selon l'une des revendications 1 à 9.10
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0955984A FR2949590B1 (fr) | 2009-09-02 | 2009-09-02 | Procede de reconstruction tridimensionnelle d'un objet a partir d'une seule vue |
| GB1014014.3A GB2473326B (en) | 2009-09-02 | 2010-08-23 | Process for three-dimentional reconstruction of an object from a single view |
| JP2010188946A JP5702971B2 (ja) | 2009-09-02 | 2010-08-26 | 放射線撮像デバイス |
| US12/873,797 US8588500B2 (en) | 2009-09-02 | 2010-09-01 | Process for three-dimensional reconstruction of an object from a single view |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR0955984A FR2949590B1 (fr) | 2009-09-02 | 2009-09-02 | Procede de reconstruction tridimensionnelle d'un objet a partir d'une seule vue |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| FR2949590A1 true FR2949590A1 (fr) | 2011-03-04 |
| FR2949590B1 FR2949590B1 (fr) | 2012-08-03 |
Family
ID=42261711
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| FR0955984A Expired - Fee Related FR2949590B1 (fr) | 2009-09-02 | 2009-09-02 | Procede de reconstruction tridimensionnelle d'un objet a partir d'une seule vue |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US8588500B2 (fr) |
| JP (1) | JP5702971B2 (fr) |
| FR (1) | FR2949590B1 (fr) |
| GB (1) | GB2473326B (fr) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB201113683D0 (en) * | 2011-08-09 | 2011-09-21 | Imorphics Ltd | Image processing method |
| US9986931B2 (en) | 2011-09-08 | 2018-06-05 | Apn Health, Llc | Automatically determining 3D catheter location and orientation using 2D fluoroscopy only |
| DE102012207261A1 (de) * | 2012-05-02 | 2013-11-07 | Siemens Aktiengesellschaft | Chirurgisches Hilfsmittel |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3589505B2 (ja) * | 1995-06-09 | 2004-11-17 | 株式会社日立メディコ | 3次元画像処理表示装置 |
| US6351513B1 (en) * | 2000-06-30 | 2002-02-26 | Siemens Corporate Research, Inc. | Fluoroscopy based 3-D neural navigation based on co-registration of other modalities with 3-D angiography reconstruction data |
| DE10210647A1 (de) * | 2002-03-11 | 2003-10-02 | Siemens Ag | Verfahren zur Bilddarstellung eines in einen Untersuchungsbereich eines Patienten eingebrachten medizinischen Instruments |
| WO2004055734A1 (fr) * | 2002-12-18 | 2004-07-01 | Koninklijke Philips Electronics N.V. | Procede et appareil de determination de la pose d'un implant |
| JP2007275346A (ja) * | 2006-04-07 | 2007-10-25 | Toshiba Corp | 放射線画像処理装置、放射線画像処理方法及び放射線画像処理プログラム |
-
2009
- 2009-09-02 FR FR0955984A patent/FR2949590B1/fr not_active Expired - Fee Related
-
2010
- 2010-08-23 GB GB1014014.3A patent/GB2473326B/en not_active Expired - Fee Related
- 2010-08-26 JP JP2010188946A patent/JP5702971B2/ja not_active Expired - Fee Related
- 2010-09-01 US US12/873,797 patent/US8588500B2/en not_active Expired - Fee Related
Non-Patent Citations (2)
| Title |
|---|
| BAERT S A M ET AL: "Guide Wire Reconstruction and Visualization in 3DRA Using Monoplane Fluoroscopic Imaging", IEEE TRANSACTIONS ON MEDICAL IMAGING, IEEE SERVICE CENTER, PISCATAWAY, NJ, US LNKD- DOI:10.1109/TMI.2005.844073, vol. 24, no. 5, 1 May 2005 (2005-05-01), pages 612 - 623, XP011131262, ISSN: 0278-0062 * |
| ESTHAPPAN J ET AL: "3D ORIENTATIONS OF CATHETERS FROM SINGLE PROJECTIONS", COMPUTER-AIDED DIAGNOSIS IN MEDICAL IMAGING, ELSEVIER, AMSTERDAM, NL, 20 September 1998 (1998-09-20), pages 357 - 361, XP008056669 * |
Also Published As
| Publication number | Publication date |
|---|---|
| GB2473326A (en) | 2011-03-09 |
| JP5702971B2 (ja) | 2015-04-15 |
| GB201014014D0 (en) | 2010-10-06 |
| GB2473326B (en) | 2015-01-07 |
| FR2949590B1 (fr) | 2012-08-03 |
| US8588500B2 (en) | 2013-11-19 |
| JP2011050738A (ja) | 2011-03-17 |
| US20110135183A1 (en) | 2011-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1222636B1 (fr) | Reconstitution statistique de surfaces en trois dimensions | |
| EP0925556B1 (fr) | Procede de reconstruction d'une image tridimensionnelle d'un objet, en particulier une image tridimensionnelle angiographique | |
| EP0389333B1 (fr) | Procédé d'acquisition de données radiologiques relatives à un corps irradié et de reconstruction de structures correspondant à ce corps | |
| 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 | |
| FR2700909A1 (fr) | Dispositif et procédé automatique de calibration géométrique d'un système d'imagerie par rayons X. | |
| FR2802002A1 (fr) | Procede de recalage automatique d'images tridimensionnelles | |
| EP0485265B1 (fr) | Dispositif et procédé de contrÔle non destructif à acquisition simultanée de données radiographiques et de données tomographiques | |
| EP2453408A1 (fr) | Procédé de traitement d'images radiologiques pour la détection d'une sténose | |
| FR2641099A1 (fr) | ||
| FR2953964A1 (fr) | Procede de traitement d'images obtenues par tomographie ou tomosynthese a faible nombre de projections | |
| FR2797978A1 (fr) | Procede de recalage automatique d'images | |
| EP2446820A1 (fr) | Dispositif d'imagerie par rayons X avec un bras en C et une unité anti-collisions | |
| EP0573364B1 (fr) | Procédé de reconstruction d'images tridimensionnelles d'un objet par des mesures utilisant un rayonnement conique et un réseau bidimensionnel de détecteurs | |
| FR2949590A1 (fr) | Procede de reconstruction tridimensionnelle d'un objet a partir d'une seule vue | |
| FR2954086A1 (fr) | Dispositif et procede d'imagerie medicale | |
| EP3368919B1 (fr) | Collimateur tournant pour determiner la position d'un element muni de capteurs dans un systeme d'imagerie par rayons x | |
| FR2944203A1 (fr) | Systeme et procede pour determiner la position d'un instrument medical | |
| EP4010739B1 (fr) | Collimateur tournant pour un système de détection de rayons x | |
| FR3042106A1 (fr) | Procede de determination d'une dose de rayonnement appliquee a un patient | |
| EP4154035A1 (fr) | Procede de traitement des donnees relatives a un examen radiologique d'un patient | |
| FR2927719A1 (fr) | Procede de traitement d'images obtenues par tomosynthese et dispositif associe | |
| FR2733841A1 (fr) | Procede de realisation de la cartographie d'emission d'un corps corrigee de l'attenuation par ce corps | |
| EP3987312A1 (fr) | Camera compton et procede d'imagerie 3d campton | |
| FR2879327A1 (fr) | Procede et ensemble d'imagerie pour l'assistance d'un utilisateur lors d'une procedure de radiologie interventionnelle | |
| WO1998029762A1 (fr) | Dispositif et procede de localisation nucleaire par calcul de barycentre pondere de detecteurs fonctionnant en parallele, et application aux gamma-cameras |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PLFP | Fee payment |
Year of fee payment: 8 |
|
| PLFP | Fee payment |
Year of fee payment: 9 |
|
| PLFP | Fee payment |
Year of fee payment: 10 |
|
| PLFP | Fee payment |
Year of fee payment: 11 |
|
| PLFP | Fee payment |
Year of fee payment: 12 |
|
| PLFP | Fee payment |
Year of fee payment: 13 |
|
| ST | Notification of lapse |
Effective date: 20230505 |