DE102006045423A1 - Verfahren zur Nachbearbeitung eines dreidimensionalen Bilddatensatzes einer Gefäßstruktur - Google Patents
Verfahren zur Nachbearbeitung eines dreidimensionalen Bilddatensatzes einer Gefäßstruktur Download PDFInfo
- Publication number
- DE102006045423A1 DE102006045423A1 DE102006045423A DE102006045423A DE102006045423A1 DE 102006045423 A1 DE102006045423 A1 DE 102006045423A1 DE 102006045423 A DE102006045423 A DE 102006045423A DE 102006045423 A DE102006045423 A DE 102006045423A DE 102006045423 A1 DE102006045423 A1 DE 102006045423A1
- Authority
- DE
- Germany
- Prior art keywords
- data set
- dsa
- projection
- vessel
- image data
- 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
- 230000002792 vascular Effects 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 title claims description 30
- 238000002583 angiography Methods 0.000 title claims description 11
- 238000012805 post-processing Methods 0.000 title claims description 4
- 238000004590 computer program Methods 0.000 claims abstract description 5
- 238000002408 directed self-assembly Methods 0.000 claims description 19
- 230000011218 segmentation Effects 0.000 claims description 15
- 238000002591 computed tomography Methods 0.000 claims description 5
- 238000002059 diagnostic imaging Methods 0.000 claims description 3
- 241001465754 Metazoa Species 0.000 claims description 2
- 230000002490 cerebral effect Effects 0.000 claims description 2
- 238000003325 tomography Methods 0.000 claims 1
- 230000006978 adaptation Effects 0.000 abstract description 2
- 206010002329 Aneurysm Diseases 0.000 description 11
- 238000002595 magnetic resonance imaging Methods 0.000 description 5
- 238000003384 imaging method Methods 0.000 description 3
- 238000002560 therapeutic procedure Methods 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 2
- 210000000988 bone and bone Anatomy 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 208000019553 vascular disease Diseases 0.000 description 2
- 241000974840 Ellipes Species 0.000 description 1
- 208000031481 Pathologic Constriction Diseases 0.000 description 1
- TZCXTZWJZNENPQ-UHFFFAOYSA-L barium sulfate Chemical compound [Ba+2].[O-]S([O-])(=O)=O TZCXTZWJZNENPQ-UHFFFAOYSA-L 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000001795 light effect Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004393 prognosis Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- 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/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/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
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/36—Image-producing devices or illumination devices not otherwise provided for
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B17/00—Surgical instruments, devices or methods, e.g. tourniquets
- A61B2017/00743—Type of operation; Specification of treatment sites
- A61B2017/00778—Operations on blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/36—Image-producing devices or illumination devices not otherwise provided for
- A61B2090/364—Correlation of different images or relation of image positions in respect to the body
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/36—Image-producing devices or illumination devices not otherwise provided for
- A61B90/37—Surgical systems with images on a monitor during operation
- A61B2090/376—Surgical systems with images on a monitor during operation using X-rays, e.g. fluoroscopy
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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/10072—Tomographic images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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; CALCULATING OR 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/30101—Blood vessel; Artery; Vein; Vascular
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10—TECHNICAL SUBJECTS COVERED BY FORMER USPC
- Y10S—TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y10S378/00—X-ray or gamma ray systems or devices
- Y10S378/901—Computer tomography program or processor
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Theoretical Computer Science (AREA)
- Biophysics (AREA)
- High Energy & Nuclear Physics (AREA)
- Radiology & Medical Imaging (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Vascular Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Physiology (AREA)
- Dentistry (AREA)
- Cardiology (AREA)
- Architecture (AREA)
- Computer Graphics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Die Erfindung betrifft ein Verfahren zur Nachbearbeitung eines 3-D-Bilddatensatzes einer Gefäßstruktur eines menschlichen oder tierischen Körpers, bei welchem eine 2-D-DSA (digitale Subtraktionsangiographie) der Gefäßstruktur aufgenommen und mit dem 3-D-Bilddatensatz registriert wird. Die 2-D-DSA wird mit einem entsprechenden aus dem 3-D-Datensatz berechneten Projektionsbild verglichen und dieser wird verändert, z. B. durch Veränderung der Segmentierungsparameter, um ihn an die 2-D-DSA anzupassen. Damit kann die hervorragende Ortsauflösung der 2-D-DSA zur Verbesserung des 3-D-Bilddatensatzes verwendet werden.
Description
- Die Erfindung betrifft ein Verfahren zur Nachbearbeitung eines dreidimensionalen (3D) Bilddatensatzes einer Gefäßstruktur eines menschlichen oder tierischen Körpers, wobei der dreidimensionale Bilddatensatz mit einer medizinischen Bildgebungsmodalität aufgenommen wurde.
- Zur Darstellung von Blutgefäßen stehen heutzutage Aufnahmemöglichkeiten zur Verfügung, mit denen dreidimensionale Bilddatensätze erzeugt werden können. Hierfür können verschiedene Bildgebungsmodalitäten verwendet werden, z.B. Computertomographie (CT), Magnetresonanztomographie (MR) oder 3D Rotationsangiographie.
- Wichtige Anwendungen sind dabei die Diagnose von Gefäßerkrankungen wie Aneurysmen und Stenosen und die Therapieplanung. Bei der Planung kommt es auf die Prognose von Rupturwahrscheinlichkeiten von Aneurysmen und die davon abhängige Therapiewahl oder die Auswahl von geeigneten therapeutischen Systemen (z.B. Stents) und deren Abmessungen (z.B. Durchmesser und Länge) an. Speziell bei der Einschätzung der Rupturwahrscheinlichkeit von Aneurysmen kann ein dreidimensionaler Bilddatensatz der Gefäßstruktur als Basis für eine Computersimulation (mittels „Computational Fluid Dynamics") dienen, mit der die Rupturwahrscheinlichkeit berechnet werden kann. Eine wichtige Größe ist hier z.B. der Durchmesser des Halses eines Aneurysma. Die exakte Einschätzung des Aneurysma-Halses kann unter Umständen die Entscheidung beeinflussen, ob dieses durch einen operativen Eingriff entfernt wird, oder stattdessen ein so genanntes intravaskuläres Coiling durchgeführt wird.
- Es ist daher sehr wichtig, als Basis für die Therapieplanung ein möglichst exaktes dreidimensionales Abbild des Gefäßes, besonders im Umfeld der Gefäßerkrankung, zu haben.
- Die mit gegenwärtig verwendeten Bildgebungsmodalitäten erzeugten dreidimensionalen Bilddatensätze besitzen jedoch keine hohe Ortsauflösung. Beim MR ist das Signal/Rausch-Verhältnis der limitierende Faktor, der die Auflösung des MR-Bildes auf ca. 1 mm3 begrenzt. Ein Röntgenbild hat zwar im Prinzip eine hohe Ortsauflösung; durch die Rekonstruktion vieler Röntgenbilder zu einem dreidimensionalen Bilddatensatz, z.B. beim CT oder bei der 3D-Rotationsangiographie, geht jedoch ein großer Teil der Ortsauflösung verloren. Um die Röntgendosis für den Patienten zu minimieren, können daher auch mit diesem Verfahren keine Auflösungen von wesentlich mehr als 0,2 bis 1 mm3 erreicht werden.
- Ferner werden dreidimensionale Bilddatensätze nach der Aufnahme und Rekonstruktion in der Regel segmentiert, wobei der Datensatz in Segmente, d.h. Volumenbereiche, aufgeteilt wird, die jeweils der Gefäßstruktur oder dem Hintergrund zugeordnet sind. Die Bildintensität des Hintergrunds wird auf Null gesetzt. Dies dient dazu, die Gefäße ohne das umgebende Gewebe und Knochen darzustellen.
- Speziell im Bereich des Halses eines Aneurysma ist es sehr schwierig, eine Segmentierung korrekt durchzuführen. Hierfür verantwortlich ist wiederum die mangelnde Ortsauflösung im rekonstruierten 3D-Bilddatensatz sowie Rekonstruktionsartefakte. Falls die Segmentierung daher auf ein spezielles Gebiet (z.B. der Aneurysma-Hals) optimiert ist, führt dies unter Umständen zu einer suboptimalen Segmentierung anderer Bereiche (z.B. Zuführende Gefäße, Aneurysma-Bauch).
- Zur Darstellung von Gefäßen ist ferner die so genannte digitale Subtraktions-Angiographie (DSA) bekannt. Dabei werden von der Gefäßstruktur zwei zeitlich aufeinander folgende Röntgenbilder gemacht, üblicherweise mit einem C-Bogen Rönt gengerät. Zwischen den Röntgenbildern wird ein Kontrastmittel in die Blutbahn injiziert. Die beiden Röntgenbilder unterscheiden sich also lediglich in der Abbildung der Gefäße, die in dem ersten Bild (Maskenbild) kaum zu sehen sind, in dem zweiten Bild (Füllungsbild) jedoch stark kontrastiert sind. Die digitalisierten Bilder werden voneinander abgezogen. Somit sind in dem Differenzbild, der DSA, lediglich die kontrastierten Blutgefäße zu sehen. Die DSA liefert somit zweidimensionale (2D) Bilder mit einer hohen Ortsauflösung, jedoch ohne Tiefeninformation. Daher wird die DSA im Folgenden auch mit „2D DSA" bezeichnet.
- Die Erfindung hat sich die Aufgabe gesetzt, ein Verfahren zur Nachbearbeitung eines 3D-Bilddatensatzes einer Gefäßstruktur zur Verfügung zu stellen, mit dem die Qualität der Darstellung der Gefäßstruktur verbessert werden kann.
- Hierzu stellt die Erfindung ein Verfahren gemäß Anspruch 1 bereit, welches die folgenden Schritte aufweist:
- (a) Bereitstellen eines 3D-Bilddatensatzes der Gefäßstruktur, der mit einem medizinischen Bildgebungsmodalität aufgenommen wurde;
- (b) Bereitstellen einer 2D DSA der Gefäßstruktur in einer ersten Projektionsrichtung;
- (c) Segmentieren des 3D-Bilddatensatzes und/oder der 2D DSA, wobei diese in Segmente aufgeteilt werden, die jeweils der Gefäßstruktur oder dem Hintergrund zugeordnet sind;
- (d) Registrieren der 2D DSA mit dem 3D-Bilddatensatz;
- (e) Berechnen einer Projektion des 3D-Bilddatensatzes auf eine Bildebene in der ersten Projektionsrichtung, zur Erzeugung eines berechneten Projektionsbildes; und
- (f) Vergleichen des berechneten Projektionsbildes mit der 2D DSA und automatisches Verändern des 3D-Bilddatensatzes zum Anpassen des berechneten Projektionsbildes an die 2D DSA.
- Der 3D-Bilddatensatz kann von einem CT oder MR stammen, kann aber auch durch 3D Rotationsangiographie in einem Rotations lauf auf der gleichen Angiographie-Anlage gewonnen werden, mit der auch die 2D DSA angefertigt werden. Das Verfahren kann entweder durchgeführt werden, während der Patient in der Angiographie-Anlage, z.B. einem C-Bogen Röntgengerät, gelagert ist, oder als Nachbearbeitung später geschehen. Der 3D-Bilddatensatz kann z.B. gewonnen werden durch zwei voneinander subtrahierte Rotationsläufe eines C-Bogen Röntgengeräts, bei welchem der zweite Lauf (Füllungslauf) kontrastiert ist. Von diesen Bildern wird der erste Lauf ohne Kontrastmittel (Maskenlauf) abgezogen, und die Bildserie zu einer dreidimensionalen Darstellung der Gefäßstruktur rekonstruiert.
- Vorzugsweise wird dann auf dem 3D-Bilddatensatz eine möglichst gute erste Arbeitsprojektionsrichtung gewählt, in der z.B. das Aneurysma möglichst überlagerungsfrei dargestellt ist. Es können auch gleich mehrere Projektionsrichtungen ausgewählt werden, wie unten noch genauer erläutert. Vorzugsweise werden dabei Projektionsrichtungen, die aufgrund der Geometrie des C-Bogens nicht eingestellt werden können, automatisch verboten.
- Daraufhin wird eine 2D DSA der Gefäßstruktur in der ausgewählten ersten Projektionsrichtung bereitgestellt, z.B. mit einem C-Bogen-System, eventuell auch mit einem Biplan-System, aufgenommen. In dem Spezialfall, dass der 3D-Bilddatensatz aus zwei voneinander subtrahierten Rotationsläufen eines C-Bogen Röntgensystems rekonstruiert wurde, sind diese 2D DSAs bereits in den Rotationsläufen enthalten. Vorzugsweise sind nur diese bereits vorhandenen Angulationen einstellbar oder anzeigbar.
- Daraufhin werden der 3D-Bilddatensatz und/oder die 2D DSA segmentiert. Das bedeutet, sie werden in Volumenbereiche bzw. Bildbereich aufgeteilt, die entweder der Gefäßstruktur oder dem Hintergrund zugeordnet sind. Solche Volumen- bzw. Bildbereiche können beliebige Formen aufweisen und werden im Folgenden „Segmente" genannt. Die Segmentierung geschieht vorzugsweise durch ein Schwellwert-Verfahren, d.h. alle Pixel bzw. Voxel, die über dem Schwellwert liegen, entweder auf einen gemeinsamen, hohen Wert gesetzt oder unverändert gelassen, während alle Pixel oder Voxel, die unter dem Schwellwert liegen, auf Null gesetzt werden. Auch dies dient dazu, die Gefäße ohne das umgebende Gewebe und Knochen darzustellen. Der Schwellwert wird hierzu automatisch oder manuell angepasst und ist jeweils für einen Bilddatensatz konstant (globaler Schwellwert). Alternativ kann die Segmentierung auch mittels einer Transferfunktion durchgeführt werden.
- Daraufhin wird die 2D DSA mit dem 3D-Bilddatensatz registriert. Dies bedeutet, dass die jeweiligen Koordinatensysteme der beiden Bildsätze z.B. durch eine Transformationsmatrix zueinander in Beziehung gesetzt werden, sodass man aus der Position eines Voxels im 3D-Bilddatensatz unmittelbar die Position des diesem Raumpunkt entsprechenden Pixels auf der 2D DSA ermitteln kann. Umgekehrt ist die Funktion nicht eindeutig, da zu jedem Pixel der 2D DSA ein Strahl von Voxeln durch den 3D-Bilddatensatz gehört.
- Die Registrierung kann entweder durch so genannte 2D-3D-Registrierungsverfahren erreicht werden, bei denen beispielsweise verschiedene Projektionen aus dem 3D-Bilddatensatz berechnet und mit der 2D DSA verglichen werden. In dem Sonderfall, in dem der 3D-Bilddatensatz mit der gleichen Angiographie-Anlage gewonnen wurde wie die 2D DSA, kann die Registrierung leichter aus der bekannten Anlagengeometrie ermittelt werden.
- Daraufhin wird in Schritt (e) eine Projektion des 3D-Bilddatensatzes auf eine Bildebene in der ausgewählten ersten Projektionsrichtung berechnet und somit ein berechnetes Projektionsbild erzeugt. Dieses wird in Schritt (f) mit der 2D DSA verglichen und an diese angepasst, indem der 3D-Bilddatensatz verändert wird. Somit wird also der 3D-Bilddatensatz nachbearbeitet, was im Folgenden anhand von Ausführungsbeispielen noch genauer erläutert wird. Salopp ausgedrückt, wird eine 2D DSA mit dem entsprechend berechneten Projekti onsbild „überlagert" und Abweichungen zu der 2D DSA werden erkannt und ggf. korrigiert.
- Optional kann der 3D-Bilddatensatz vor Schritt (f) auf eine höhere Ortsauflösung hochskaliert werden, damit die genaueren Ortsinformationen der 2D DSA nicht in der Pixelunschärfe des 3D-Bilddatensatzes untergehen.
- Vorzugsweise werden die Schritte (c) bis (f) mit einer zweiten 2D DSA der Gefäßstruktur wiederholt, die aus einer zweiten Projektionsrichtung aufgenommen wird. Diese zweite Projektionsrichtung liegt vorzugsweise in einem Winkel von ungefähr 30° bis 110°, besonders bevorzugt 80° bis 100°, zur ersten Projektionsrichtung. Am meisten bevorzugt ist ein Winkel von ungefähr 90°. Dies erlaubt eine Anpassung des 3D-Bilddatensatzes an hoch aufgelöste 2D DSAs aus verschiedenen Blickrichtungen.
- Gemäß einer besonders bevorzugten Ausführungsform werden mehrere 2D DSAs der Gefäßstruktur, z.B. 3 bis 5, in mehreren Projektionsrichtungen bereitgestellt und die Schritte (c) bis (f) mit diesen mehreren 2D DSAs und den entsprechenden Projektionsrichtungen durch den 3D-Bilddatensatz wiederholt. Der 3D-Bilddatensatz wird dadurch iterativ verändert und verbessert.
- Der Schritt (f) des Vergleichens und automatischen Veränderns wird vorzugsweise folgendermaßen durchgeführt:
- (f1) die der Gefäßstruktur zugeordneten Segmente auf dem berechneten Projektionsbild und auf der 2D DSA werden miteinander verglichen und Abweichungen erkannt;
- (f2) die abweichenden Segmentbereiche werden durch Rückprojektion von dem berechneten Projektionsbild in das Volumen des 3D-Bilddatensatzes einem Gefäßabschnitt in dem 3D-Bilddatensatz zugeordnet; und
- (f3) der 3D-Bilddatensatz wird durch Verschieben der Segmentgrenzen zwischen dem Gefäßabschnitt und dem Hintergrund so verändert, dass das dem Gefäßabschnitt zugeordnete Segment auf einem neu berechneten Projektionsbild stärker mit dem diesem Gefäßabschnitt zugeordneten Segment auf der 2D DSA übereinstimmt.
- Schritt (f2) ist deshalb notwendig, weil aus der 2D-Projektion auf dem berechneten Projektionsbild auf die 3D-Position des anzupassenden Gefäßes rückgeschlossen werden muss. Daher ist es sinnvoll, eine Projektionsrichtung auszuwählen, in der die Gefäßstruktur möglichst überlagerungsfrei dargestellt wird. Außerdem muss die Anlagengeometrie bekannt sein. Wird dann die Projektion eines Gefäßabschnitts anhand der Anlagengeometrie zum Röntgenfokus rückverfolgt, kann auf diesem Weg nur ein Gefäß liegen, dessen 3D Position dann bekannt ist. Vorzugsweise wird diese Methode zur Bestimmung der 3D-Position der abweichenden Segmentbereiche in dem Volumen des 3D-Bilddatensatzes benutzt.
- Ist dann die 3D-Position des abweichenden Gefäßabschnitts im 3D-Bilddatensatz bekannt, kann der 3D-Bilddatensatz so verändert werden, dass die Abweichungen verringert werden. Hierzu werden im Folgenden zwei Ausführungsbeispiele genauer erläutert. Beim ersten Ausführungsbeispiel werden die Segmentgrenzen zwischen dem Gefäßabschnitt und dem Hintergrund durch pixelweises Verändern eines Segmentierungs-Schwellwertes des 3D-Bilddatensatzes verschoben. Das heißt, die Segmentierungsparameter des 3D-Bilddatensatzes werden lokal (nicht global) angepasst, um eine optimale 3D-Darstellung zu gewährleisten. Optional können die pixelweisen Segmentierungs-Schwellwerte untereinander geglättet werden, um abrupte Übergänge zu vermeiden.
- Bei der zweiten Ausführungsform wird aus dem 3D-Bilddatensatz ein Gefäßbaum-Modell der Gefäßstruktur berechnet, welches die Centrelines der die Gefäßstruktur bildenden Gefäße sowie Werte für deren Durchmesser beinhaltet. Die Segmentgrenzen zwischen dem Gefäßabschnitt und dem Hintergrund werden dann durch Verändern von Durchmesserangaben im Gefäßbaumodell verschoben. Die Durchmesserangaben können sich auch auf Ellip sendurchmesser beziehen, sodass pro Gefäßabschnitt also zwei Durchmesserangaben und eine Winkelangabe, die die Ausrichtung der langen Achse der Ellipse angibt, vorhanden sind und gegebenenfalls verändert werden.
- Besonders bevorzugt wird das Verfahren bei einer zerebralen Gefäßstruktur verwendet.
- Die Erfindung ist auch auf ein Computerprogrammprodukt gerichtet, welches auf einem computerlesbaren Medium gespeicherte Softwarecodeabschnitte enthält. Diese können eine Recheneinheit dazu veranlassen, das erfindungsgemäße Verfahren auszuführen, wenn das Computerprogrammprodukt auf der Recheneinheit installiert ist. Bei der Recheneinheit handelt es sich vorzugsweise um einen in der medizinischen Bildgebung verwendeten Bildverarbeitungscomputer.
- Die Erfindung wird nun anhand von Ausführungsbeispielen in Bezug auf die beiliegenden Zeichnungen näher erläutert. In den Zeichnungen zeigen:
-
1 eine beispielhafte Darstellung einer 2D DSA einer Gefäßstruktur; -
2 eine Darstellung eines 3D-Bilddatensatzes derselben Gefäßstruktur in „volume rendering" Darstellung; -
3 eine Überlagerung der 2D DSA der1 mit einem aus dem 3D-Bilddatensatz der2 berechneten Projektionsbild; -
4 eine schematische, pixelweise Darstellung der Überlagerung der3 ; -
5 eine schematische Darstellung des 3D-Bilddatensatzes mit einer berechneten Projektion und 2D DSA in einer ersten Projektionsrichtung; -
6 die schematische Darstellung des 3D-Bilddatensatzes der5 mit einer berechneten Projektion und einer 2D DSA in einer zweiten Projektionsrichtung; -
7 eine schematische Darstellung eines Gefäßbaum-Modells mit einer berechneten Projektion und einer 2D DSA in einer ersten Projektionsrichtung; und -
8 eine schematische Darstellung des Gefäßbaum-Modells der7 mit einer berechneten Projektion und einer 2D DSA in einer zweiten Projektionsrichtung. -
1 zeigt beispielhaft eine 2D DSA1 einer Gefäßstruktur2 mit einem Aneurysma2a . Wie man erahnen kann, weist die 2D DSA hervorragende Ortsauflösung auf, es fehlt jedoch die Tiefeninformation. -
2 zeigt dagegen eine dreidimensionale Darstellung3 desselben Gefäßbaums, die durch einen 3D-Rotationslauf mit einem C-Bogen Röntgengerät gewonnen wurde. Der 3D-Bilddatensatz wurde mit einem globalen Segmentierungs-Schwellwert segmentiert und nur die über dem Schwellwert liegenden Werte werden dargestellt. Gewählt wurde eine so genannte „volume rendering" Darstellung, in der die Gefäßstruktur im 3D-Volumen mit berechneten Schatten- und Lichteffekten versehen wurde, um einen visuellen dreidimensionalen Eindruck zu erzeugen. - Gerade bei 3D-Bilddatensätzen ist der gewählte globale Segmentierungs-Schwellwert sehr wichtig, da insbesondere ein Aneurysma mit verschiedenen Schwellwerten deutlich unterschiedlich dargestellt wird.
- Erfindungsgemäß wird aus einem 3D-Bilddatensatz wie bei der
3 in der gleichen Projektionsrichtung, in der die 2D DSA aufgenommen wurde, eine Projektion aus dem 3D-Bilddatensatz berechnet. Ein solches berechnetes Projektionsbild4 ist in3 in hellgrau dargestellt und mit der DSA1 der1 (in schwarz dargestellt) überlagert. Wie aus3 erkennbar ist, decken sich die Projektion der Gefäßstruktur11 auf der 2D DSA1 und die Darstellung der Gefäßstruktur14 auf dem berechneten Projektionsbild4 nicht exakt, vielmehr sind Abweichungen vorhanden. Diese Abweichungen werden mit den im Folgenden beschriebenen Verfahren zur Nachbearbeitung und Verbesserung des 3D-Bilddatensatzes verwendet, da man davon ausgeht, dass die 2D DSA im Zweifelsfall genauere Ergebnisse liefert. -
4 stellt schematisch einen fiktiven Ausschnitt aus der3 dar. Der mit11 bezeichnete Kasten soll die Pixel eines Gefäßabschnitts auf der 2D DSA darstellen, während der Kasten14 einen Gefäßabschnitt in dem berechneten Projektionsbild darstellt. Die fetten Linien zeigen jeweils die Segmentgrenzen zwischen Hintergrund9 und den Gefäßabschnitten11 bzw.14 an. Wie man sieht, überlappen sich die Gefäßsegmente11 und14 nur teilweise. In den Abschnitten A und B herrscht keinerlei Abweichung zwischen den beiden Darstellungen, da entweder Hintergrund auf Hintergrund oder Vordergrund auf Vordergrund liegt. In den Bereichen C ist das Gefäß auf der DSA11 dagegen breiter als die berechnete Projektion14 des 3D-Bilddatensatzes. Die entsprechende Segmentgrenze des Gefäßabschnitts14 sollte also verschoben werden. Dies kann beispielsweise dadurch geschehen, dass der lokale Schwellwert für die Pixel im Bereich C verringert wird, sodass die Darstellung des Gefäßabschnitts14 sich in diesem Bereich mit dem Gefäßabschnitt11 zumindest annähernd deckt. - In den Abschnitten D ist der Gefäßabschnitt
14 in dem 3D-Bilddatensatz gegenüber der DSA zu breit, der lokale Schwellwert in den Bereichen D ist also zu klein. Wird der Schwellwert hier erhöht, verschieben sich die Segmentgrenzen rechts und links im Gefäßabschnitt14 in Richtung der Segmentgrenzen des Gefäßabschnitts11 . Auf diese Weise kann die Gefäßdarstellung14 im 3D-Bilddatensatz mit der Gefäßdarstellung11 auf der 2D DSA1 in stärkere Übereinstimmung gebracht werden. - Anhand der
5 und6 wird nun das erste Ausführungsbeispiel genauer erläutert. Bei diesem Ausführungsbeispiel werden die Segmentgrenzen durch die pixelweise Anpassung von Segmentierungs-Schwellwerten verschoben. Es wird also nicht ein globaler Schwellwert zur Segmentierung des 3D-Bilddatensatzes verwendet, sondern ein gesonderter Schwellwert für jedes Voxel, welcher iterativ angepasst wird, bis der 3D-Bilddatensatz der oder den 2D DSAs angepasst ist. -
5 zeigt einen schematischen Schnitt durch einen solchen 3D-Bilddatensatz6 . Die einzelnen Voxel sind wiederum als Kästchen dargestellt. - Die von links unten nach rechts oben schraffierten Voxel
10 stellen die Voxel dar, die nach der Segmentierung einem bestimmten Gefäßabschnitt8 innerhalb des 3D-Bilddatensatzes zugeordnet wurden. Unter8' und8'' sind weitere, kleinere Gefäße dargestellt. Mit7 ist der Röntgenfokus bezeichnet, von dem die Projektionskegel21 und24 ausgehen. Der Projektionskegel24 ist die Projektion des Gefäßes8 im 3D-Bilddatensatz. Auf die Bildebene14 projiziert, ergibt dies eine zweidimensionale Darstellung des Gefäßes8 auf dem berechneten Projektionsbild14 . Dieses berechnete Projektionsbild wird überlagert mit einer 2D DSA11 . Auf dieser ist der Gefäßabschnitt in dieser Projektionsrichtung deutlich schmaler dargestellt. Verfolgt man den Gefäßabschnitt11 somit anhand der dazugehörigen Projektionsstrahlen21 zum Röntgenfokus7 zurück, wird nur ein kleinerer Teil des Gefäßabschnitts8 hiervon abgedeckt. Die quer schraffierten Voxel12 sind somit zwar in dem 3D-Bilddatensatz dem Gefäßabschnitt zugeordnet, nicht jedoch gemäß der 2D DSA. In dem Voxel12 sollte somit der Schwellwert lokal erhöht werden, um die beiden in stärkere Übereinstimmung zu bringen. - Dieses Verfahren wird vorzugsweise aus einer zweiten Projektionsrichtung wiederholt, wie in
6 dargestellt. Hier ist das gleiche Gefäß8 im 3D-Bilddatensatz gezeigt. Eine Projektion der hiervon umfassten Voxel10 in der zweiten Projekti onsrichtung vom Röntgenfokus7 auf die Bildebene14 ergibt somit ein weiteres berechnetes Projektionsbild14' . Der hier zugehörige Projektionskegel (der Pfad der Röntgenstrahlen) ist mit24' bezeichnet. - Auch in dieser Projektionsrichtung wird eine 2D DSA
11' aufgenommen. Wie in der Zeichnung zu erkennen ist, ist auf dieser 2D DSA das Gefäß breiter als im berechneten Projektionsbild14' . Verfolgt man die Röntgenprojektion entlang der Strahlen21' zum Röntgenfokus7 zurück, trifft man auf die senkrecht schraffierten Voxel15 , die gemäß der 2D DSA zum Gefäß gehören sollten. Für diese Voxel15 muss also der lokale Schwellwert vermindert werden. - Genauer kann das Verfahren wie folgt iterativ ablaufen:
Eine DSA11 und eine berechnete Projektion14 werden miteinander verglichen und die Voxel12 bzw.15 ermittelt, für die der Schwellwert lokal erhöht oder vermindert werden soll. Der Schwellwert wird dann um einen Wert erhöht oder vermindert, der z.B. nach der Anzahl der abweichenden Voxel festgelegt wird, oder aber einer vorbestimmten Tabelle entnommen wird. - Optional werden daraufhin die Schwellwerte räumlich geglättet, um abrupte Übergänge zu vermeiden.
- Daraufhin wird eine neue Projektion des mit den veränderten Schwellwerten segmentierten 3D-Bilddatensatzes berechnet. Dieses neue berechnete Projektionsbild wird wieder mit der gleichen oder mit einer zweiten 2D DSA überlagert und verglichen. Daraufhin wird wieder der lokale Schwellwert in den abweichenden Pixeln erhöht oder vermindert. Optional werden die lokalen Schwellwerte räumlich geglättet.
- Diese Schritte werden wiederholt, bis eine maximale Übereinstimmung zwischen berechnetem Projektionsbild und DSA erreicht ist, oder aber bis ein oder mehrere lokale Schwellwer te so weit vom ursprünglichen Schwellwert abweichen, dass eine weitere Verbesserung nicht mehr zu erwarten ist.
- Anhand der
7 und8 wird nun die zweite Ausführungsform näher erläutert.7 zeigt einen Schnitt durch ein Gefäßbaummodell, welches aus dem 3D-Bilddatensatz6 der5 und6 berechnet worden ist. In dem dargestellten Bereich umfasst das Gefäßbaummodell die Centreline16 eines Gefäßes, sowie den ellipsenförmigen Durchmesser dieses Gefäßes34 . Somit ist der Querschnitt des Gefäßes in dem Gefäßbaummodell hierdurch definiert und durch die von links unten nach rechts oben schraffierte Fläche dargstellt. - Aus dem Gefäßbaummodell wird, ähnlich wie bei dem ersten Ausführungsbeispiel, anhand der bekannten Anlagengeometrie ein Projektionsbild
14 berechnet, indem die Umrisse des Gefäßbaummodells34 ausgehend vom Röntgenfokus7 entlang der Strahlen24 auf die Bildebene14 projiziert werden. Dieses Bild wird mit einer 2D DSA11 verglichen. Die Umrisse des Gefäßes auf der 2D DSA werden umgekehrt entlang der Strahlen21 zum Röntgenfokus7 rückverfolgt. Dabei wird deutlich, dass der Durchmesser des Gefäßes in der Richtung quer zu den Strahlen21 auf der 2D DSA schmaler ist als das Gefäßbaummodell34 . Dementsprechend wird in dieser Richtung ein neuer Ellipsendurchmesser berechnet und das Gefäßbaummodell entsprechend angepasst. Das angepasste Modell31 ist senkrecht schraffiert dargestellt. - Diese Anpassung kann gemäß
8 in einer zweiten Projektionsrichtung, welche ungefähr senkrecht zur ersten Projektionsrichtung steht, wiederholt werden. Hier wird nun entlang der Strahlen24' eine Projektion des durch die Centrelines16 und die Ellipse34 definierten Gefäßes auf die Bildebene14' berechnet. Dieses berechnete Projektionsbild14' wird mit einer 2D DSA11' , die in der gleichen Projektionsrichtung aufgenommen wurde, verglichen. Anhand der Strahlen21' zum Röntgenfokus7 wird das Gefäß34 ermittelt, welches in dem Gefäßbaummodell dem Gefäßabschnitt11' der 2D DSA entspricht. Da bei wird festgestellt, dass dieser Gefäßabschnitt auf der 2D DSA11' breiter ist als auf dem berechneten Projektionsbild14' . Somit wird der Ellipsendurchmesser in diese Richtung verbreitert, um das neue Gefäßbaummodell31' zu erhalten. - Obgleich dies auf den Figuren nicht dargestellt wird, könnte zur Anpassung des Gefäßbaummodells auch die Centreline
16 verschoben werden, falls dies notwendig ist.
Claims (14)
- Verfahren zur Nachbearbeitung eines 3D Bilddatensatzes (
6 ) einer Gefäßstruktur (8 ) eines menschlichen oder tierischen Körpers, mit den folgenden Schritten: (a) Bereitstellen eines 3D Bilddatensatzes (6 ) der Gefäßstruktur (8 ), der mit einer medizinischen Bildgebungsmodalität aufgenommen wurde; (b) Bereitstellen einer 2D DSA (11 ) (digitale Subtraktionsangiographie) der Gefäßstruktur (8 ) in einer ersten Projektionsrichtung (21 ); (c) Segmentieren des 3D Bilddatensatzes (6 ) und/oder der 2D DSA (11 ), wobei diese in Segmente aufgeteilt werden, die jeweils der Gefäßstruktur (8 ) oder dem Hintergrund (9 ) zugeordnet sind; (d) Registrieren der 2D DSA (11 ) mit dem 3D Bilddatensatz (6 ); (e) Berechnen einer Projektion des 3D Bilddatensatzes (6 ) auf eine Bildebene in der ersten Projektionsrichtung (24 ), zur Erzeugung eines berechneten Projektionsbildes (14 ); und (f) Vergleichen des berechneten Projektionsbildes (14 ) mit der 2D DSA (11 ) und automatisches Verändern des 3D-Bilddatensatzes zum Anpassen des berechneten Projektionsbildes an die 2D DSA. - Verfahren nach Anspruch 1, bei welchem die Segmentierung des 3D Bilddatensatzes (
6 ) und der 2D DSA (11 ) durch ein Schwellwert-Verfahren erfolgt. - Verfahren nach einem der vorhergehenden Ansprüche, bei welchem eine zweite 2D DSA (
11 ) der Gefäßstruktur in einer zweiten Projektionsrichtung (21' ) aufgenommen wird und die Schritte (c) bis (f) mit dieser zweiten 2D DSA (11' ) und der zweiten Projektionsrichtung (21' ) wiederholt werden. - Verfahren nach Anspruch 3, bei welchem die zweite Projektionsrichtung (
21' ) in einem Winkel von ungefähr 30° bis 110° zur ersten Projektionsrichtung (21 ) steht. - Verfahren nach einem der Ansprüche 3 oder 4, bei welchem mehrere 2D DSAs (
11 ,11' ) der Gefäßstruktur in mehreren Projektionsrichtungen (21 ,21' ) bereitgestellt werden und die Schritte (c) bis (f) mit diesen mehreren 2D DSAs und den entsprechenden mehreren Projektionsrichtungen durch den 3D Bilddatensatz wiederholt werden. - Verfahren nach einem der vorhergehenden Ansprüche, bei welchem Schritt (f) folgendermaßen durchgeführt wird: (f1) die der Gefäßstruktur (
8 ) zugeordneten Segmente auf dem berechneten Projektionsbild (14 ,14' ) und auf der 2D DSA (11 ,11' ) werden miteinander verglichen und Abweichungen erkannt; (f2) die abweichenden Segmentbereiche (12 ,15 ) werden durch Rückprojektion von dem berechneten Projektionsbild (14 ) in das Volumen des 3D Bilddatensatzes (6 ) einem Gefäßabschnitt (8 ) im 3D Bilddatensatz zugeordnet; und (f3) der 3D Bilddatensatz (6 ) wird durch Verschieben der Segmentgrenzen zwischen dem Gefäßabschnitt (8 ) und dem Hintergrund (9 ) so verändert, dass das dem Gefäßabschnitt (8 ) zugeordnete Segment auf einem neu berechneten Projektionsbild stärker mit dem diesem Gefäßabschnitt zugeordneten Segment auf der 2D DSA (11 ,11' ) übereinstimmt. - Verfahren nach Anspruch 6, bei welchem die Segmentgrenzen zwischen dem Gefäßabschnitt (
8 ) und dem Hintergrund (9 ) in Schritt (f3) durch pixelweises Verändern eines Segmentierungs-Schwellwertes des 3D Bilddatensatzes (6 ) verschoben werden. - Verfahren nach Anspruch 7, bei welchem die pixelweisen Segmentierungs-Schwellwerte vor dem Berechnen eines neuen Projektionsbildes in Schritt (f3) räumlich geglättet werden.
- Verfahren nach Anspruch 6, bei welchem aus dem 3D Bilddatensatz (
6 ) ein Gefäßbaummodell (16 ,34 ) der Gefäßstruktur berechnet wird, wobei das Gefäßbaummodell (16 ,34 ) die Centrelines (16 ) der die Gefäßstruktur bildenden Gefäße sowie Werte für deren Durchmesser (34 ) beinhaltet, und wobei die Segmentgrenzen zwischen dem Gefäßabschnitt und dem Hintergrund (9 ) durch Verändern von Durchmesserangaben im Gefäßbaummodell (16 ,34 ) verschoben werden. - Verfahren nach einem der vorhergehenden Ansprüche, bei welchem die Gefäßstruktur (
8 ) eine zerebrale Gefäßstruktur ist. - Verfahren nach einem der vorhergehenden Ansprüche, bei welchem die 2D DSAs (
11 ,11' ) während des Verfahrens mit einem C-Bogen Röntgengerät aufgenommen werden. - Verfahren nach einem der vorhergehenden Ansprüche, bei welchem der 3D Bilddatensatz (
6 ) durch Magnetresonanztomographie oder Computertomographie aufgenommen wurde. - Verfahren nach einem der vorhergehenden Ansprüche, bei welchem der 3D Bilddatensatz (
6 ) durch 3D Rotationsangiographie aufgenommen wurde, wobei bei der Aufnahme der 3D Rotationsangiographie auch die 2D DSAs zumindest teilweise aufgenommen wurden. - Computerprogrammprodukt, enthaltend auf einem Computerlesbaren Medium gespeicherte Softwarecodeabschnitte, welche eine Recheneinheit dazu veranlassen, das Verfahren nach einem der vorhergehenden Ansprüche auszuführen, wenn das Computerprogrammprodukt auf der Recheneinheit installiert ist.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102006045423.5A DE102006045423B4 (de) | 2006-09-26 | 2006-09-26 | 07.09.07Verfahren zur Nachbearbeitung eines dreidimensionalen Bilddatensatzes einer Gefäßstruktur |
US11/901,674 US7903856B2 (en) | 2006-09-26 | 2007-09-18 | Method for post-processing a three-dimensional image data set of vessel structure |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102006045423.5A DE102006045423B4 (de) | 2006-09-26 | 2006-09-26 | 07.09.07Verfahren zur Nachbearbeitung eines dreidimensionalen Bilddatensatzes einer Gefäßstruktur |
Publications (2)
Publication Number | Publication Date |
---|---|
DE102006045423A1 true DE102006045423A1 (de) | 2008-04-10 |
DE102006045423B4 DE102006045423B4 (de) | 2016-07-14 |
Family
ID=39154469
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE102006045423.5A Expired - Fee Related DE102006045423B4 (de) | 2006-09-26 | 2006-09-26 | 07.09.07Verfahren zur Nachbearbeitung eines dreidimensionalen Bilddatensatzes einer Gefäßstruktur |
Country Status (2)
Country | Link |
---|---|
US (1) | US7903856B2 (de) |
DE (1) | DE102006045423B4 (de) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102009011725A1 (de) * | 2009-03-04 | 2010-10-28 | Siemens Aktiengesellschaft | Verfahren zur Bildunterstützung bei der Navigation eines medizinischen Instruments und Vorrichtung zur Durchführung eines minimalinvasiven Eingriffs zur Therapie eines Tumors |
DE102011078220A1 (de) | 2011-06-28 | 2012-08-30 | Siemens Aktiengesellschaft | Verfahren zum Bereitstellen von zumindest zwei Abwandlungen eines 3-D-Bilddatensatzes zu einer Gefäßstruktur eines Patienten |
DE102011005634A1 (de) | 2011-03-16 | 2012-09-20 | Friedrich-Alexander-Universität Erlangen-Nürnberg | Verfahren zum Bereitstellen eines nachbearbeiteten 3D-Bilddatensatzes, insbesondere zu einer Gefäßstruktur |
DE102012205935A1 (de) * | 2012-04-12 | 2013-10-17 | Siemens Aktiengesellschaft | Verfahren zur Aufnahme eines vierdimensionalen Angiographie-Datensatzes |
US10319091B2 (en) | 2015-12-04 | 2019-06-11 | Siemens Healthcare Gmbh | Providing image support to a practitioner |
DE102011083703B4 (de) | 2011-09-29 | 2021-12-30 | Siemens Healthcare Gmbh | Aufnahme von Bildern für einen Bilddatensatz der digitalen Subtraktionsangiographie und einen Perfusionsdatensatz |
Families Citing this family (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9427201B2 (en) | 2007-06-30 | 2016-08-30 | Accuray Incorporated | Non-invasive method for using 2D angiographic images for radiosurgical target definition |
US8200466B2 (en) | 2008-07-21 | 2012-06-12 | The Board Of Trustees Of The Leland Stanford Junior University | Method for tuning patient-specific cardiovascular simulations |
US20100026783A1 (en) * | 2008-08-01 | 2010-02-04 | Real D | Method and apparatus to encode and decode stereoscopic video data |
US9405886B2 (en) | 2009-03-17 | 2016-08-02 | The Board Of Trustees Of The Leland Stanford Junior University | Method for determining cardiovascular information |
WO2011008906A1 (en) * | 2009-07-15 | 2011-01-20 | Mayo Foundation For Medical Education And Research | Computer-aided detection (cad) of intracranial aneurysms |
US8643642B2 (en) | 2009-08-17 | 2014-02-04 | Mistretta Medical, Llc | System and method of time-resolved, three-dimensional angiography |
US8654119B2 (en) * | 2009-08-17 | 2014-02-18 | Mistretta Medical, Llc | System and method for four dimensional angiography and fluoroscopy |
WO2011091300A2 (en) | 2010-01-24 | 2011-07-28 | Mistretta Medical, Llc | System and method for implementation of 4d time-energy subtraction computed tomography |
DE102010013221B4 (de) * | 2010-03-29 | 2012-06-28 | Siemens Aktiengesellschaft | Verfahren zur hervorgehobenen Darstellung von Objekten bei interventionellen angiographischen Untersuchungen sowie Vorrichtung zur Durchführung des Verfahrens |
WO2011128797A1 (en) * | 2010-04-15 | 2011-10-20 | Koninklijke Philips Electronics N.V. | Instrument-based image registration for fusing images with tubular structures |
BR112013001487B1 (pt) * | 2010-07-21 | 2022-05-03 | Armin E. Moehrle | Método e aparelho de criação de relatório de imagem |
US20120022366A1 (en) * | 2010-07-21 | 2012-01-26 | Marcus Pfister | Registration of aorta to patient via two 2d images for placement of a stent |
US8315812B2 (en) | 2010-08-12 | 2012-11-20 | Heartflow, Inc. | Method and system for patient-specific modeling of blood flow |
US8768031B2 (en) | 2010-10-01 | 2014-07-01 | Mistretta Medical, Llc | Time resolved digital subtraction angiography perfusion measurement method, apparatus and system |
US20120128265A1 (en) * | 2010-11-23 | 2012-05-24 | Toshiba Medical Systems Corporation | Method and system utilizing iterative reconstruction with adaptive parameters for computer tomography (ct) images |
EP2468207A1 (de) | 2010-12-21 | 2012-06-27 | Renishaw (Ireland) Limited | Verfahren und Vorrichtung zur Bildanalyse |
DE102011005777B4 (de) * | 2011-03-18 | 2018-04-26 | Siemens Healthcare Gmbh | Angiographiesystem zur angiographischen Untersuchung eines Patienten und angiographisches Untersuchungsverfahren |
WO2012174263A2 (en) | 2011-06-15 | 2012-12-20 | Mistretta Medical, Llc | System and method for four dimensional angiography and fluoroscopy |
US8667430B1 (en) | 2012-10-24 | 2014-03-04 | GlobalFoundries, Inc. | Methods for directed self-assembly process/proximity correction |
US8667428B1 (en) * | 2012-10-24 | 2014-03-04 | GlobalFoundries, Inc. | Methods for directed self-assembly process/proximity correction |
CN104837525B (zh) * | 2012-12-11 | 2019-04-05 | 皇家飞利浦有限公司 | 用于辅助执行近距离放射治疗的辅助装置 |
US9091628B2 (en) | 2012-12-21 | 2015-07-28 | L-3 Communications Security And Detection Systems, Inc. | 3D mapping with two orthogonal imaging views |
US9009634B2 (en) | 2013-07-08 | 2015-04-14 | GlobalFoundries, Inc. | Methods for fabricating integrated circuits including generating photomasks for directed self-assembly |
US9170501B2 (en) | 2013-07-08 | 2015-10-27 | GlobalFoundries, Inc. | Methods for fabricating integrated circuits including generating photomasks for directed self-assembly |
US9023730B1 (en) | 2013-11-05 | 2015-05-05 | GlobalFoundries, Inc. | Methods for fabricating integrated circuits including generating e-beam patterns for directed self-assembly |
DE102014201559A1 (de) * | 2014-01-29 | 2015-07-30 | Siemens Aktiengesellschaft | Angiographisches Untersuchungsverfahren eines Gefäßsystems in einer interessierenden Körperregion eines Patienten |
US9286434B2 (en) | 2014-05-23 | 2016-03-15 | GlobalFoundries, Inc. | Methods for fabricating integrated circuits including generating photomasks for directed self-assembly (DSA) using DSA target patterns |
US9888968B2 (en) * | 2014-07-22 | 2018-02-13 | Siemens Healthcare Gmbh | Method and system for automated therapy planning for arterial stenosis |
US9517042B2 (en) * | 2014-09-12 | 2016-12-13 | General Electric Company | Systems and methods for imaging phase selection for computed tomography imaging |
DE102016215971A1 (de) * | 2016-08-25 | 2018-03-01 | Siemens Healthcare Gmbh | Segmentierung einer Angiographie unter Verwendung einer bestehenden dreidimensionalen Rekonstruktion |
US10770175B2 (en) | 2017-09-15 | 2020-09-08 | Multus Medical Llc | System and method for segmentation and visualization of medical image data |
EP3581111A1 (de) * | 2018-06-13 | 2019-12-18 | Siemens Healthcare GmbH | Verfahren und darstellungsvorrichtung zur nachbearbeitung und anzeige eines dreidimensionalen angiographie-bilddatensatzes, computerprogramm und elektronisch lesbares speichermedium |
US11103142B2 (en) * | 2019-04-02 | 2021-08-31 | Tencent America LLC | System and method for predicting vertebral artery dissection |
CN114723683B (zh) * | 2022-03-22 | 2023-02-17 | 推想医疗科技股份有限公司 | 头颈动脉血管分割方法及其装置、电子设备和存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004072903A2 (en) * | 2003-02-12 | 2004-08-26 | Philips Intellectual Property & Standards Gmbh | Method for the 3d modeling of a tubular structure |
US20060188139A1 (en) * | 2004-12-20 | 2006-08-24 | Ali Khamene | Tree structure based 2D to 3D registration |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7602970B2 (en) * | 2005-03-21 | 2009-10-13 | Siemens Medical Solutions Usa, Inc. | System and method for Kalman filtering in vascular segmentation |
US7715626B2 (en) * | 2005-03-23 | 2010-05-11 | Siemens Medical Solutions Usa, Inc. | System and method for vascular segmentation by Monte-Carlo sampling |
US7822254B2 (en) * | 2006-04-21 | 2010-10-26 | Siemens Medical Solutions Usa, Inc. | Automatic positioning of matching multi-planar image reformatting (MPR) views of multiple 3D medical images |
-
2006
- 2006-09-26 DE DE102006045423.5A patent/DE102006045423B4/de not_active Expired - Fee Related
-
2007
- 2007-09-18 US US11/901,674 patent/US7903856B2/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004072903A2 (en) * | 2003-02-12 | 2004-08-26 | Philips Intellectual Property & Standards Gmbh | Method for the 3d modeling of a tubular structure |
US20060188139A1 (en) * | 2004-12-20 | 2006-08-24 | Ali Khamene | Tree structure based 2D to 3D registration |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102009011725A1 (de) * | 2009-03-04 | 2010-10-28 | Siemens Aktiengesellschaft | Verfahren zur Bildunterstützung bei der Navigation eines medizinischen Instruments und Vorrichtung zur Durchführung eines minimalinvasiven Eingriffs zur Therapie eines Tumors |
DE102011005634A1 (de) | 2011-03-16 | 2012-09-20 | Friedrich-Alexander-Universität Erlangen-Nürnberg | Verfahren zum Bereitstellen eines nachbearbeiteten 3D-Bilddatensatzes, insbesondere zu einer Gefäßstruktur |
DE102011078220A1 (de) | 2011-06-28 | 2012-08-30 | Siemens Aktiengesellschaft | Verfahren zum Bereitstellen von zumindest zwei Abwandlungen eines 3-D-Bilddatensatzes zu einer Gefäßstruktur eines Patienten |
DE102011083703B4 (de) | 2011-09-29 | 2021-12-30 | Siemens Healthcare Gmbh | Aufnahme von Bildern für einen Bilddatensatz der digitalen Subtraktionsangiographie und einen Perfusionsdatensatz |
DE102012205935A1 (de) * | 2012-04-12 | 2013-10-17 | Siemens Aktiengesellschaft | Verfahren zur Aufnahme eines vierdimensionalen Angiographie-Datensatzes |
US9036780B2 (en) | 2012-04-12 | 2015-05-19 | Siemens Aktiengesellschaft | Method for recording a four-dimensional angiography data record |
DE102012205935B4 (de) * | 2012-04-12 | 2018-11-15 | Siemens Healthcare Gmbh | Verfahren zur Aufnahme eines vierdimensionalen Angiographie-Datensatzes |
US10319091B2 (en) | 2015-12-04 | 2019-06-11 | Siemens Healthcare Gmbh | Providing image support to a practitioner |
Also Published As
Publication number | Publication date |
---|---|
DE102006045423B4 (de) | 2016-07-14 |
US20080212857A1 (en) | 2008-09-04 |
US7903856B2 (en) | 2011-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE102006045423B4 (de) | 07.09.07Verfahren zur Nachbearbeitung eines dreidimensionalen Bilddatensatzes einer Gefäßstruktur | |
EP1114615B1 (de) | Vorrichtung zur Darstellung des zeitlichen Verlaufs des Blutflusses | |
DE602004002939T2 (de) | Methode zur dreidimensionalen modellierung von rohrförmigen strukturen | |
DE102013218819B3 (de) | Verfahren zur Reduzierung von Artefakten in einem Bilddatensatz und Röntgeneinrichtung | |
DE102006049865B4 (de) | Auswertungsverfahren für Daten eines Untersuchungsobjekts unter Korrelierung eines Volumendatensatzes mit einer Anzahl von Durchleuchtungsbildern | |
DE102009032059A1 (de) | Sinogrammbearbeitung für die Metallartefaktreduktion in der Computertomographie | |
DE102014201134B4 (de) | Verfahren und Vorrichtung zur Erzeugung eines 2-D-Projektionsbildes eines Gefäßsystems nebst korrespondierenden Gegenständen | |
DE10144020A1 (de) | Dreidimensionales Bilderzeugungsverfahren und -Vorrichtung und zugehöriges Röntgengerät | |
DE102005044653A1 (de) | Verfahren und Vorrichtung zur Rekonstruktion eines dreidimensionalen Bildvolumens aus zweidimensionalen Projektionsbildern | |
DE102014201559A1 (de) | Angiographisches Untersuchungsverfahren eines Gefäßsystems in einer interessierenden Körperregion eines Patienten | |
DE102012212774A1 (de) | Verfahren zur Korrektur von Metallartefakten und Röntgeneinrichtung | |
DE102013202313A1 (de) | Verfahren und Vorrichtung zur Korrektur von Bewegungsartefakten bei einem computertomographischen Bild | |
EP2242023A1 (de) | Verfahren zur bewegungskompensierten Rekonstruktion eines dreidimensionalen endgültigen Bilddatensatzes | |
DE112017001477T5 (de) | Verbesserte bildqualität durch zweistufige temporale rauschunterdrückung | |
DE102008038357B3 (de) | Verfahren zur Erzeugung von 2D-Schnittbildern aus 3D-Projektionsdaten, die mittels eines CT-Systems von einem metallische Anteile enthaltenden Untersuchungsobjekt erfasst wurden | |
DE102008045449B4 (de) | Verfahren zur Erstellung von computertomographischen Aufnahmen eines Patienten mit metallischen Anteilen und Computersystem zur Durchführung des Verfahrens | |
EP3616621B1 (de) | Verfahren zur 3d dsa und vorrichtung | |
DE102010040041B3 (de) | Verfahren zur Korrektur von durch zeitliche Veränderungen von Schwächungswerten auftretenden Artefakten | |
EP1174086A1 (de) | Computertomographie-Verfahren mit kegelförmiger Durchstrahlung eines Objekts | |
DE102008050570A1 (de) | Verfahren zum Erzeugen eines 3D-Bilddatensatzes | |
DE102010022526A1 (de) | Verfahren und Vorrichtung zur Verbesserung einer CT-Bild-Darstellung von zumindest einem medizinischen Instrument in einem nahezu nicht kontrastierbaren Ziel-Bereich | |
DE102011075912A1 (de) | Verfahren zum Bereitstellen eines 3D-Bilddatensatzes mit unterdrückten Bildartefakten und Computertomograph | |
DE102009037251A1 (de) | Verfahren zum Erzeugen von 3D-Bilddaten eines Körpers | |
DE102022201003B3 (de) | Verfahren zum Festlegen einer Aufnahmetrajektorie | |
EP2187353B1 (de) | Verfahren zum Visualisieren von Strukturen in einem Körper, insbesondere im Körper eines Patienten |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
OP8 | Request for examination as to paragraph 44 patent law | ||
8120 | Willingness to grant licences paragraph 23 | ||
R079 | Amendment of ipc main class |
Free format text: PREVIOUS MAIN CLASS: G06T0017400000 Ipc: G06T0019000000 Effective date: 20130314 |
|
R018 | Grant decision by examination section/examining division | ||
R081 | Change of applicant/patentee |
Owner name: SIEMENS HEALTHCARE GMBH, DE Free format text: FORMER OWNER: SIEMENS AKTIENGESELLSCHAFT, 80333 MUENCHEN, DE |
|
R020 | Patent grant now final | ||
R119 | Application deemed withdrawn, or ip right lapsed, due to non-payment of renewal fee |