DE102004048209B3 - Verfahren und Vorrichtung zur Erzeugung eines dreidimensionalen Bilddatensatzes eines bewegten Objekts mittels Röntgentomographie - Google Patents
Verfahren und Vorrichtung zur Erzeugung eines dreidimensionalen Bilddatensatzes eines bewegten Objekts mittels Röntgentomographie Download PDFInfo
- Publication number
- DE102004048209B3 DE102004048209B3 DE102004048209A DE102004048209A DE102004048209B3 DE 102004048209 B3 DE102004048209 B3 DE 102004048209B3 DE 102004048209 A DE102004048209 A DE 102004048209A DE 102004048209 A DE102004048209 A DE 102004048209A DE 102004048209 B3 DE102004048209 B3 DE 102004048209B3
- Authority
- DE
- Germany
- Prior art keywords
- image data
- data set
- motion
- raw images
- preliminary
- 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.)
- Expired - Fee Related
Links
- 230000033001 locomotion Effects 0.000 title claims abstract description 60
- 239000011159 matrix material Substances 0.000 title claims abstract description 26
- 238000003325 tomography Methods 0.000 title claims description 6
- 238000000034 method Methods 0.000 claims abstract description 54
- 238000011156 evaluation Methods 0.000 claims description 19
- 238000012545 processing Methods 0.000 claims description 8
- 230000000747 cardiac effect Effects 0.000 claims description 6
- 125000004122 cyclic group Chemical group 0.000 claims description 5
- 238000009499 grossing Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 2
- 210000004351 coronary vessel Anatomy 0.000 description 10
- 238000009825 accumulation Methods 0.000 description 7
- 238000002583 angiography Methods 0.000 description 6
- 210000004204 blood vessel Anatomy 0.000 description 6
- 238000001514 detection method Methods 0.000 description 6
- 230000002123 temporal effect Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 238000002591 computed tomography Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000002203 pretreatment Methods 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 210000001367 artery Anatomy 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 208000010201 Exanthema Diseases 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000000746 body region Anatomy 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 201000005884 exanthem Diseases 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 206010037844 rash Diseases 0.000 description 1
- 230000001373 regressive effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/33—Heart-related electrical modalities, e.g. electrocardiography [ECG] specially adapted for cooperation with other devices
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- 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/025—Tomosynthesis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
- A61B5/352—Detecting R peaks, e.g. for synchronising diagnostic apparatus; Estimating R-R interval
-
- 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/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/503—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 the heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/504—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of blood vessels, e.g. by angiography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/54—Control of apparatus or devices for radiation diagnosis
- A61B6/541—Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/412—Dynamic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/421—Filtered back projection [FBP]
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Signal Processing (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Optics & Photonics (AREA)
- High Energy & Nuclear Physics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Cardiology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Es werden ein präzises und vergleichsweise einfach zu realisierendes Verfahren zur Erzeugung eines dreidimensionalen Bilddatensatzes (E, E¶1¶) eines bewegten Objekts (2) mittels Röntgentomographie sowie eine zur Durchführung des Verfahrens besonders geeignete Vorrichtung (1) mit einer drehbar gelagerten Röntgenstrahler-Detektor-Einheit (4, 5) sowie mit einer Auswerteeinheit (16) angegeben. Erfindungsgemäß ist vorgesehen, eine Anzahl von zweidimensionalen Rohbildern (R¶i¶) nach Maßgabe einer zyklischen Relativzeit (t') zu gruppieren, jeweils aus einander gemäß dieser Gruppierung entsprechenden Rohbildern (R¶i¶) mindestens zwei vorläufige 3-D-Bilddatensätze (V¶j¶) zu erzeugen, durch Vergleich jeweils zweier als Quelldatensatz (V¶l¶) bzw. Zieldatensatz (V¶k¶) herangezogener vorläufiger 3-D-Bilddatensätze (V¶j¶) mindestens eine Bewegungsmatrix (B¶lk¶) abzuleiten, durch Anwendung der oder jeder Bewegungsmatrix (B¶lk¶) auf den zugehörigen Zieldatensatz (V¶k¶) einen einer Referenzzeit (T¶l¶) des Quelldatensatzes (V¶l¶) entsprechenden bewegungskompensierten 3-D-Bilddatensatz (V¶k¶'') zu erzeugen und den oder jeden bewegungskompensierten 3-D-Bilddatensatz (V¶k¶'') mit mindestens einem weiteren bewegungskompensierten 3-D-Bilddatensatz (V¶k¶'') oder einem sonstigen, der gleichen Referenzeinheit (T¶l¶) entsprechenden vorläufigen 3-D-Bilddatensatz (V¶l¶) aufzusummieren.
Description
- Die Erfindung bezieht sich auf ein Verfahren zur Erzeugung eines dreidimensionalen Bilddatensatzes eines bewegten Objekts mittels Röntgentomographie. Die Erfindung bezieht sich des Weiteren auf eine Vorrichtung zur Durchführung des genannten Verfahrens. Als Röntgentomographie wird allgemein ein bildgebendes Röntgenverfahren zur Erzeugung eines dreidimensionalen Bilddatensatzes (Tomogramm) eines Objekts verstanden. Der Begriff Röntgentomographie umfasst hierbei neben klassischer Computertomographie insbesondere die so genannte 3D-Rotationsangiographie.
- Als Angiographie wird allgemein die Darstellung der Blutgefäße, insbesondere der Arterien, Venen und Herzkammern, des (im Allgemeinen menschlichen) Körpers mittels bildgebender Röntgenverfahren bezeichnet. Bei der 3D-Rotationsangiographie werden ähnlich wie bei der Computertomographie mehrere zweidimensionale Röntgenbilder (Rohbilder) eines Patienten unter unterschiedlichen Aufnahmeprojektionen aufgenommen. Aus diesen Rohbildern wird anschließend ein dreidimensionaler Bilddatensatz der zu untersuchenden Blutgefäße erzeugt. Insbesondere bei der Untersuchung des Herzens und der dieses umgebenen Koronararterien besteht hierbei das aufnahmetechnische Problem, dass diese Blutgefäße in Folge des Herzschlags einer ständigen Bewegung ausgesetzt sind. Diese Bewegung führt herkömmlicherweise zu starken Bildartefakten bei der Rekonstruktion des 3D-Bilddatensatzes, die die Bildqualität des letzteren mitunter erheblich beeinträchtigen.
- Eine Reduzierung dieser Bewegungsartefakte bei der 3D-Rotationsangiographie des Herzens und herznaher Blutgefäße wurde bisher mittels einer als „EKG-Gating" bezeichneten Methode erzielt. Bei einem derartigen Verfahren, wie es beispielsweise aus O. Wink et al. „Coronary Intervention Planning, Using Hybrid 3D Reconstruction", MICCAI 2002, LNCS 2488 (Springer), S. 604-611, 2002 bekannt ist, werden die Rohbilder mittels eines EKG-Signals mit dem Herzrhythmus des Patienten synchronisiert, wobei nur solche Rohbilder zur Rekonstruktion des 3D-Bilddatensatzes herangezogen werden, die bei einer vorgegebenen bewegungsarmen Phase des Herzzyklus aufgenommen sind. Dies führt jedoch zu einer starken Reduktion der Anzahl verwendbarer Rohbilder. So können bei Anwendung von EKG-Gating von den üblicherweise mehreren hundert Rohbildern, die während einer kompletten Rotation der Röntgenstrahler-Detektoreinheit eines 3D-Rotationsangiographen aufgenommen werden, nur wenige Rohbilder verwendet werden. Diese Reduktion der Bildanzahl hat wiederum eine starke Verminderung der Bildqualität, sowohl in Hinblick auf die räumliche Auflösung als auch auf die Kontrastauflösung, zur Folge.
- Aus Ch. Blondel, et al., „4D-Deformation Field of Coronary Arteries from Monoplane Rotational X-Ray Angiography", Internat. Congress Ser. 1256 (Elsevier), S. 1073-1087, 2003 und Ch. Blondel, et al. „4D-Tomographic Representation of Coronary Arteries from One Rotational X-Ray Sequence", MICCAI 2003, LNCS 2878 (Springer), S. 416-423, 2003 sind des Weiteren Verfahren bekannt, bei welchen den Bildqualität einer tomographischen Darstellung der Koronararterien durch Kompensation der Arterienbewegung weiter verbessert wird. Diese Verfahren beruhen auf einer Segmentierung der Koronararterien aus den bei einem bestimmten Punkt des Herzzyklus aufgenommenen Rohbildern, Modellierung der Koronararterienbewegung anhand eines aus dieser Segmentierung gewonnenen 3D-Punktmodells des Koronararterienbaums und anschließender Erstellung eines deformationskompensierten Tomogramms aus den Rohbildern. Nachteiligerweise ist ein derartiges Verfahren mit vergleichsweise hohem numerischem Aufwand verbunden.
- Aus der
DE 101 29 631 A1 ist ein Verfahren zur Rekonstruktion eines hoch aufgelösten 3D-Bildes bekannt, wonach aus gruppierten Rohbildern vorläufige 3D-Bilddatensätze erzeugt werden, aus denen Bewegungsinformationen im 3D-Raum gewonnen werden. - In der WO 2004/081877 A1 ist ein Verfahren zur Rekonstruktion bewegungskompensierter 3D-Volumenbilder beschrieben, bei dem 2D-Rohbilder bewegungskompensiert werden, aus denen ein 3D-Bild generiert wird.
- Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren anzugeben, das die Erzeugung eines hochqualitativen dreidimensionalen Bilddatensatzes eines bewegten Objekts mittels Röntgentomographie auf vergleichsweise einfache Weise ermöglicht. Der Erfindung liegt weiterhin die Aufgabe zugrunde, eine zur Durchführung des genannten Verfahrens besonders geeignete Vorrichtung anzugeben.
- Bezüglich des Verfahrens wird die Aufgabe erfindungsgemäß gelöst durch die Merkmale des Anspruchs 1. Bezüglich der zugehörigen Vorrichtung wird die Aufgabe erfindungsgemäß gelöst durch die Merkmale des Anspruchs 10.
- Die erfindungsgemäße Vorrichtung umfasst danach eine Röntgenstrahler-Detektoreinheit, die an einem C-Bogen od.dgl. drehbar bevorzugt um eine isozentrische Achse gelagert ist, so dass mittels der Röntgenstrahler-Detektoreinheit Rohbilder eines etwa im Isozentrum gelagerten Objekts unter verschiedenem Projektionswinkel aufgenommen werden können. Bei dem Objekt handelt es sich insbesondere um einen zu untersuchenden Patienten oder ein Körperteil desselben. Die Vorrichtung umfasst weiterhin eine Auswerteeinheit zur Erzeugung eines dreidimensionalen Bilddatensatzes des Objekts anhand der aufgenommenen Rohbilder.
- Die Rohbilder werden hierzu zunächst nach Maßgabe einer zyklischen Relativzeit synchronisiert, d.h. in Gruppen mit vergleichbarer Relativzeit aufgeteilt. Als Relativzeit wird hierbei eine zeitlich zyklisch variierende Referenzgröße bezeichnet, die ggf. bevorzugt aus dem Herzzyklus des Patienten abgeleitet wird.
- Aus mindestens zwei solcher Gruppen von Rohbildern wird jeweils ein zugehöriger vorläufiger 3D-Bilddatensatz erzeugt, der entsprechend jeweils einem (nachfolgend als Referenzzeit bezeichneten) Zeitpunkt der Relativzeit entspricht. Als 3D-Bilddatensatz wird ein dreidimensionales Feld von Bildpunkten (Voxeln) bezeichnet, deren jedem ein Farbwert oder Grauwert zugeordnet ist. Zur Erzeugung der vorläufigen 3D-Bilddatensätze aus den zweidimensionalen Rohbildern wird 3D-Rekonstruktionsverfahren, insbesondere ein gefiltertes Rückprojektionsverfahren angewendet, wie es auch in der herkömmlichen 3D-Rotationsangiographie eingesetzt wird.
- Durch Vergleich zweier solcher vorläufiger 3D-Bilddatensätze wird nun eine zugehörige Bewegungsmatrix abgeleitet. Als Bewegungsmatrix wird ein dreidimensionales Vektorfeld bezeichnet, durch welches Voxel oder Voxelgruppen des einen 3D-Bilddatensatzes mit Voxeln oder Voxelgruppen des anderen 3D-Bilddatensatzes verknüpft werden. Die Bewegungsmatrix kennzeichnet hierbei eine „Bewegung", d.h. räumliche Veränderung einander entsprechender Bildstrukturen zwischen dem als Ausgangspunkt herangezogenen ersten 3D-Bilddatensatz (Quelldatensatz) und dem zeitlich späteren oder früheren zweiten 3D-Bilddatensatz (Zieldatensatz). Die Erstellung der Bewegungsmatrix aus dem Quelldatensatz und dem Zieldatensatz erfolgt bevorzugt mittels einer Korrelationsmethode, z.B. Block-Matching oder einer auf optischem Fluss basierenden Methode.
- Durch anschließende Anwendung der oder jeder Bewegungsmatrix auf den zugehörigen Zieldatensatz wird nun eine zeitliche Rückentwicklung (oder Vorausentwicklung) der Bildinformation des Zieldatensatzes simuliert. Das Resultat dieser Operation ist ein bewegungskompensierter 3D-Bilddatensatz, der der Referenzzeit des zugehörigen Quelldatensatzes entspricht.
- Der oder jeder bewegungskompensierte 3D-Bilddatensatz wird schließlich mit jedem weiteren, der gleichen Referenzzeit entsprechenden bewegungskompensierten 3D-Bilddatensatz aufsummiert. Die Summation erstreckt sich alternativ oder zusätzlich auch über mindestens einen der gleichen Referenzzeit entsprechenden vorläufigen 3D-Bilddatensatz, insbesondere über den zugehörigen Quelldatensatz.
- Das vorstehend beschriebene Verfahren ist mit vergleichsweise geringem Aufwand vollautomatisch durchzuführen und ermöglicht insbesondere eine hochpräzise dreidimensionale Darstellung der herznahen menschlichen Blutgefäße. Die Korrelation zweier vorläufiger 3D-Bilddatensätze zur Erstellung der Bewegungsmatrix ist numerisch vergleichsweise unaufwändig und stabil. Weiterhin ist das Verfahren auch dann noch durchführbar, wenn in den als Quellen- bzw. Zieldatensatz herangezogenen 3D-Bilddatensätzen keine eindeutig segmentierbaren Strukturen vorhanden oder erkennbar sind.
- In bevorzugter Ausführung des Verfahrens und der zugehörigen Vorrichtung ist vorgesehen, dass vor der Ableitung der oder jeder Bewegungsmatrix die hierzu herangezogenen vorläufigen 3D-Bilddatensätze mit Mitteln der elektronischen Bildverarbeitung, insbesondere Kontrastglättung und Grauwert- bzw. Farbwert-Fensterung vorbehandelt werden. Hierdurch werden Rauscheffekte und Rekonstruktionsartefakte, die in den vorläufigen 3D-Bilddatensätzen in vergleichsweise starkem Maße vorhanden sind, unterdrückt.
- Eine weitere Verbesserung der Bildqualität wird optional dadurch erreicht, dass die oder jede Bewegungsmatrix regularisiert wird. Als Regularisierung wird in diesem Sinne jedes Verfahren bezeichnet, das eine Homogenisierung, Rausch- und Artefaktunterdrückung, etc. bei einem Vektorfeld ermöglicht, insbesondere Glättungsverfahren in Orts- und Zeitdomäne, Fensterungsverfahren, etc. Vorzugsweise wird die Regularisierung anhand hinterlegter Information über typische Bewegungsmuster des Objekts parametriert. Insbesondere werden bevorzugt aus Erfahrungswerten der räumlichen Amplitude der Arterienbewegung im menschlichen Körper und der Herzfrequenz Schwellwerte für eine Fensterung der Bewegungsmatrix abgeleitet. Auf diese Weise können insbesondere untypisch langreichweitige oder schnelle Strukturveränderungen in der Bewegungsmatrix, die mit der herzschlagbedingten Arterienbewegung nicht vereinbar sind, als Artefakte erkannt und aus der Bewegungsmatrix entfernt werden.
- Die Relativzeit wird insbesondere aus dem Herzzyklus des zu untersuchenden Patienten abgeleitet, indem als Relativzeit eine mit dem Herzzyklus synchronisierte sägezahnartige Funktion herangezogen wird. Der Herzzyklus wird dabei insbesondere durch eine EKG-Einheit gemessen und als EKG-Signal zur Verfügung gestellt. Alternativ hierzu ist vorgesehen, dass der Herzzyklus aus den aufgenommenen Rohbildern abgeschätzt wird, z.B. indem die Bildinformation der Rohbilder einer zeitlichen Spektralanalyse, einem regressiven Anpassungsverfahren od.dgl. unterworfen wird.
- Bevorzugt ist vorgesehen, dass zur Erstellung der vorläufigen 3D-Bilddatensätze die Relativzeit zunächst in eine Anzahl gleichmäßiger Zeitfenster unterteilt wird, wobei jeweils die dem selben Zeitfenster zugeordneten Rohbilder als Basis eines vorläufigen 3D-Bilddatensatzes herangezogen werden.
- Nachfolgend werden Ausführungsbeispiele der Erfindung anhand einer Zeichnung näher erläutert. Darin zeigen:
-
1 in schematischer Darstellung eine 3D-Rotationsangiographieeinrichtung (Röntgentomograph) mit einer Patientenlagerung, mit einer um diese herum rotierbaren Röntgenstrahler-Detektor-Einheit sowie mit einem Steuer- und Auswertesystem, -
2 in einem schematischen Blockschaltbild eine Auswerteeinheit des Steuer- und Auswertesystems gemäß1 , -
3 in einem schematischen Flussdiagramm ein mittels des Röntgentomographen gemäß1 durchführbares Verfahren zur Erzeugung eines dreidimensionalen Bilddatensatzes, mit den Verfahrensschritten Roh bildaufnahme, Rohbildsynchronisation, Rohbildgruppierung (EKG-Gating), Rückprojektion, Datenvorbehandlung, Bewegungsdetektion, Datennachbehandlung, Bewegungskompensation und Datenkumulation, -
4 in einem schematischen zeitlichen Diagramm ein EKG-Signal und eine für das Verfahren gemäß3 daraus abgeleitete zyklische Relativzeit, -
5 eine schematische Prinzipdarstellung der Rohbildgruppierung und Rückprojektion im Zuge des Verfahrens gemäß3 , -
6 eine schematische Prinzipdarstellung der Bewegungsdetektion, Bewegungskompensation und Datenkumulation im Zuge des Verfahrens gemäß3 und -
7 in Darstellung gemäß3 eine iterativ durchgeführte Variante des Verfahrens. - Einander entsprechende Teile und Größen sind in allen Figuren stets mit den gleichen Bezugszeichen versehen.
- In
1 ist schematisch eine Vorrichtung1 zur Erzeugung eines dreidimensionalen Bilddatensatzes eines Objekts dargestellt. Bei der Vorrichtung1 handelt es sich um einen Röntgentomographen, insbesondere eine 3D-Rotationsangiographieanlage. Bei dem zu untersuchenden Objekt handelt es sich um einen Patienten2 , insbesondere den Brustbereich eines Patienten2 . - Die Vorrichtung
1 umfasst eine Aufnahmeeinheit3 mit einem Röntgenstrahler4 und einem Röntgendetektor5 . Röntgenstrahler4 und Röntgendetektor5 sind in Gegenüberstellung an den Enden eines so genannten C-Bogens6 angebracht. Der C-Bogen6 ist wiederum etwa mittig an einem Stativ7 um eine isozentrische Achse8 drehbar gelagert. Der Röntgenstrahler4 und der Röntgendetektor5 sind dabei durch Verschwenkung des C-Bogens6 gegenüber dem Stativ7 derart verdrehbar, dass ein Zentralstrahl9 der von dem Röntgenstrahler4 in Richtung des Röntgendetektors5 emittierten Röntgenstrahlung innerhalb einer zu der isozentrischen Achse8 senkrechten Aufnahmeebene gegenüber dem umgebenden Raum in einen beliebigen Projektionswinkel θ verschwenkbar ist, wobei der Zentralstrahl9 stets auf die isozentrische Achse8 ausgerichtet ist. Die in1 dargestellte Stellung der Aufnahmeeinheit3 , in welcher der Zentralstrahl9 senkrecht von oben nach unten gerichtet ist, entspricht einem Projektionswinkel von θ = 0°. - Die Vorrichtung
1 umfasst weiterhin einen Patiententisch10 mit einer Tischplatte11 , auf welcher der Patient2 während einer Untersuchung derart gelagert wird, dass seine Körperlängsachse etwa mit der isozentrischen Achse8 der Aufnahmeeinheit3 fluchtet. Die Tischplatte11 ist für die Untersuchung derart in die Öffnung des C-Bogens6 einschiebbar, dass der zu untersuchende Körperbereich des Patienten2 zwischen Röntgenstrahler4 und Röntgendetektor5 zu liegen kommt. - Die Vorrichtung
1 umfasst weiterhin eine EKG-Einheit12 mit einer Anzahl von EKG-Sensoren13 , die zur Aufnahme eines Elektrokardiogramms (EKG), d.h. eines die Herzaktivität des Patienten2 wiedergebenden elektrischen Signals, in bekannter Weise an dem Körper des Patienten2 befestigt werden. - Die Vorrichtung
1 umfasst weiterhin ein Steuer- und Auswertesystem14 . Das Steuer- und Auswertesystem14 umfasst eine Datenverarbeitungsanlage15 , in welcher zusätzlich zu (nicht näher dargestellten) Bedien- und Steuerfunktionen eine Auswerteeinheit16 zur Erzeugung eines dreidimensionalen (3D-) Bilddatensatzes des untersuchten Körperbereichs des Patienten2 implementiert ist. Das Steuer- und Auswertesystem14 umfasst weiterhin Ein-/Ausgabemittel17 , wie z.B. Bildschirm, Tastatur, Maus od.dgl. zur Eingabe von Steueranweisungen so wie zur Anzeige von Zustandsgrößen, Untersuchungsergebnissen, etc. - Im Zuge des von der Vorrichtung
1 durchgeführten Verfahrens werden der Auswerteeinheit16 durch die Aufnahmeeinheit3 digitale Bilddaten von zweidimensionalen Röntgenbildern, im Folgenden als Rohbilder R bezeichnet, zugeführt. Weiterhin wird der Auswerteeinheit16 seitens der EKG-Einheit12 ein EKG-Signal S des Patienten2 zugeführt. - Der Aufbau der Auswerteeinheit
16 ist in2 näher dargestellt. Danach umfasst die Auswerteeinheit16 eine Eingangsschnittstelle20 für das EKG-Signal S sowie eine weitere Eingangsschnittstelle21 für die den Rohbildern R entsprechenden Bilddaten. Die Eingangsschnittstellen20 ,21 sind durch eine Synchronisationsverbindung22 untereinander synchronisiert, so dass jedem Rohbild R ein zeitlich entsprechender Punkt des EKG-Signals S zugeordnet werden kann. - Die Auswerteeinheit
16 umfasst weiterhin ein Analysemodul23 , in welchem Algorithmen zur Bearbeitung und Modifikation der Rohbilder R und des EKG-Signals S gemäß dem nachfolgend näher beschriebenen Verfahren implementiert sind. Zur Speicherung der Rohbilder R des EKG-Signals S und nachfolgend näher beschriebener aufbereiteter Daten umfasst die Auswerteeinheit16 weiterhin ein Speichermodul24 . Die Auswerteeinheit16 weist außerdem eine Ein-/Ausgabeschnittstelle25 zur Datenkommunikation mit den Ein-/Ausgabemitteln17 auf. Die Eingabeschnittstellen20 ,21 , das Analysemodul23 , das Speichermodul24 und die Ein-/Ausgabeschnittstelle25 sind durch Hardwarekomponenten der Datenverarbeitungsanlage15 und zugehöriger Software gebildet und stehen über ein gemeinsames Bussystem26 in bilateralem Datenaustausch. - Das mittels der Vorrichtung
1 durchgeführte Verfahren umfasst gemäß3 die Verfahrensschritte Rohbildaufnahme30 , Rohbildsynchronisation31 , Rohbildgruppierung oder EKG-Gating32 , gefilterte Rückprojektion33 , Datenvorbehandlung34 , Bewegungsdetektion35 , Datennachbehandlung36 , Bewegungskompensation37 und Datenkumulation38 . - Im Zuge der Rohbildaufnahme
30 werden zunächst bei rotierender Aufnahmeeinheit3 eine Anzahl von typischerweise mehreren hundert Rohbildern R während einer Rotation des C-Bogens6 um einen Winkelbereich von typischerweise mehr als 180° aufgenommen. Jedes einzelne Rohbild Ri (i=1, 2, 3,..., N) ist hierbei durch eine Aufnahmezeit ti und einen Projektionswinkel θi gekennzeichnet, was im folgenden durch die funktionale Schreibweise Ri(ti,θi) ausgedrückt ist. Jedes Rohbild Ri umfasst ein zweidimensionales Datenfeld von Bildpunkten (Pixeln), denen jeweils ein Grau- oder Farbwert zugeordnet ist. - Im Zuge der Rohbildsynchronisation
31 wird zunächst anhand des EKG-Signals S eine so genannte zyklische Relativzeit t' ermittelt. Als Relativzeit t' wird eine zwischen 0 und 1 variierende Größe herangezogen, die sägezahnartig, und damit insbesondere periodisch mit der Aufnahmezeit t variiert. Die Periodendauer τ der Relativzeit t' wird dabei an die Periodizität des EKG-Signals S angepasst. Wie anhand von4 verdeutlicht ist, lässt sich die einem beliebigen Aufnahmezeitpunkt ti eines Rohbilds Ri zugeordnete Relativzeit ti' durch die Formel bestimmen, wobei t< und t> die durch den so genannten R-Ausschlag des EKG-Signals angezeigten Zeitpunkte bezeichnen, die die aktuelle Herzphase begrenzen. Im Zuge der Rohbildsynchronisation31 wird somit der Aufnahmezeitpunkt ti auf den entsprechenden Relativzeitpunkt ti' abgebildet. - Bei dem nachfolgenden EKG-Gating
32 werden die Rohbilder Ri(ti',θi) nach Maßgabe der Relativzeit t' gruppiert. Hierbei wird die Relativzeit t' in eine Anzahl von Zeitfenstern Fj (j = 1, 2,...., M) unterteilt, welchen jeweils ein Bemittelter Relativzeitpunkt als Referenzzeit Tj sowie eine Breite Δ zugeordnet ist. Die Zeitfenster Fj sind insbesondere derart gewählt, dass sie zeitlich unmittelbar aneinander angrenzen (5 ). Für die Referenzzeiten Tj gilt dabei mit j=1, 2,..., M. Für die Breite Δ gilt Δ = M–1. Die Breite Δ kann jedoch alternativ auch größer gewählt sein, so dass sich benachbarte Zeitfenster Fj überlappen. Zusätzlich können die Rohbilder Ri in dem jeweiligen Zeitfenster Fj gewichtet werden, so dass z.B Rohbilder Ri an den Grenzen des Zeitfensters Fj glatt zu Null gedämpft werden. - Alle einem gegebenen Zeitfenster Fj zugeordneten Rohbilder Ri(ti'∊Fj, θi) werden gruppiert, d. h. ausgewählt. Aus diesen gruppierten Rohbilder Ri wird anschließend durch gefilterte Rückprojektion
33 ein vorläufiger 3D-Bilddatensatz Vj erzeugt. Die Rückprojektion erfolgt nach einem gängigen Verfahren, insbesondere nach der Feldkamp-Methode. Der resultierende vorläufige 3D-Bilddatensatz Vj entspricht somit der Referenzzeit Tj. - Das in
5 anhand des Zeitfensters F4 beispielhaft verdeutlichte EKG-Gating32 und die anschließende Rückprojektion33 werden für jedes Zeitfenster Fj durchgeführt, so dass nach Abschluss der Rückprojektion33 zu jeder Referenzzeit Tj ein zugehöriger vorläufiger 3D-Bilddatensatz Vj(Tj) existiert. - Wie aus
5 erkennbar ist gehen in die Rekonstruktion eines vorläufigen 3D-Bilddatensatses Vj in Folge des EKG-Gatings32 nur Rohbilder Ri für bestimmte Ausschnitte des Projektionswinkels θ ein. In Folge dieser ausschnitthaften Rückprojektion enthalten die vorläufigen 3D-Bilddatensätze Vj Artefakte, die die räumliche Bildinformation beeinträchtigen. Um diese Artefakte zu reduzieren, werden die vorläufigen 3D-Bilddatensätze Vj im Zuge der Datenvorbehandlung34 mit Mitteln der elektronischen Bildverarbeitung, insbesondere Glättung und Fensterung, aufbereitet. Dabei werden insbesondere Voxel mit geringem Grau- oder Farbwert unterdrückt. - Als Ergebnis der Datenvorbehandlung
34 wird für jede Referenzzeit Tj ein aufbereiteter 3D-Bilddatensatz Vj' ausgegeben, in welchem die Rekonstruktionsartefakte unterdrückt sind, so dass stark kontrastierende Strukturen in den 3D-Bilddatensätzen Vj' stärker hervortreten. Allerdings werden infolge der Bildbearbeitung auch schwach kontrastierende Strukturen, wie sie z.B. durch kleine Nebenarterien hervorgerufen werden, unterdrückt, so dass die aufbereiteten 3D-Bilddatensätze Vj' für diagnostische Zwecke nur unzureichend geeignet wären. - Hingegen sind die aufbereiteten 3D-Bilddatensätze Vj' infolge der Kontrastschärfung zur Erkennung der Strukturbewegung zwischen zwei unterschiedlichen Referenzzeiten Tj besonders gut geeignet. Dies geschieht im Zuge der anschließenden Bewegungsdetektion
35 , deren Prinzip in6 dargestellt ist. Hierbei werden durch ein entsprechendes Detektormodul das Analysemoduls23 zwei aufbereitete 3D-Bilddatensätze Vj' im Hinblick auf die räumliche Veränderung der dreidimensionalen Bildinformation verglichen. Dies geschieht mittels gängiger Korrelationsmethoden oder durch Ermittlung des optischen Flusses. Als Ergebnis der Bewegungsdetektion35 werden eine Anzahl von Bewegungsmatrizen Blk (l,k = 1, 2,...., M) erzeugt, die die räumliche Veränderung der dreidimensionalen Bildinformation zwischen einem Quelldatensatz Vl' und einem Zieldatensatz Vk' kennzeichnen. Unter Berücksichtigung der dem Quelldatensatz Vl' und dem Zieldatensatz Vk' jeweils zugeordneten Referenzzeiten Tl bzw. Tk ist aus dieser räumlichen Strukturänderung das Geschwindigkeitsfeld einer Strukturbewegung ableitbar. - Jede Bewegungsmatrix Blk ist ein dreidimensionales Datenfeld von 3D-Vektoren, deren jeder eine Korrelation zwischen einem Voxel oder einem Voxelbereich des Quelldatensatzes Vl (bzw. Vl') und einem Voxel oder Voxelbereich des Zieldatensatzes Vk (bzw. Vk') herstellt. In
6 sind aus Gründen der Übersichtlichkeit nur Korrelationen des Quelldatensatzes Vj → Vl mit zeitlich benachbarten Zieldatensätzen Vj–1, Vj+1 → Vk dargestellt. Prinzipiell werden auf diese Weise aber bevorzugt Bewegungsmatrizen Blk für beliebige Zieldatensätze Vk mit k = 1,2,...,M und k ≠ l ermittelt. - Die Bewegungsmatrizen Blk werden nun im Zuge der Datennachbehandlung
36 regularisiert, insbesondere in Ort- und Zeit-Domäne geglättet. Bei der Regularisierung werden weiterhin typische Bewegungsmuster der Koronararterien berücksichtigt. Beispielsweise wird aus der Herzfrequenz des Patienten2 und der typischen räumlichen Bewegungsamplitude der Koronararterien im menschlichen Körper Maximalbeträge für die Vektorkomponenten der Bewegungsmatrizen Blk abgeschätzt, anhand derer Artefakte innerhalb der Bewegungsmatrizen Blk erkannt und eliminiert werden. - Als Ergebnis der Datennachbehandlung
36 resultieren aufbereitete Bewegungsmatrizen Blk'. Durch Anwendung einer Bewegungsmatrix Blk' auf den zugehörigen Zieldatensatz Vk kann nun die Zeitentwicklung der Bildinformation des jeweiligen 3D-Bilddatensatzes Vk auf die Referenzzeit Tl des Quelldatensatzes Vl zurück- bzw. vorausentwickelt werden. Zu beachten ist dabei, dass die oder jede Bewegungsmatrix Blk' nicht auf den aufbereiteten Zieldatensatz Vk', aus dem sie berechnet wurde, angewandt wird, sondern den Zieldatensatz Vk in seiner ursprünglichen Form, um die volle Strukturinformation der vorläufigen 3D-Bilddatensätze Vj zu nutzen. - Das Resultat einer solchen Anwendung, die im Zuge der Bewegungskompensation
37 vorgenommen wird, ist ein bewegungskom pensierter 3D-Bilddatensatz Vk'', der der Referenzzeit Tl des jeweiligen Quelldatensatzes Vl entspricht. Die Anwendung der Bewegungsmatrix Blk' auf den Zieldatensatz Vk wird nachfolgend symbolisch durch den Formelausdruck Vk''(Tl) = Blk'· Vk(Tk) ausgedrückt. - Als Resultat der Bewegungskompensation
37 werden somit eine Anzahl von bewegungskompensierten 3D-Bilddatensätzen Vk'' ausgegeben, die alle der selben Referenzzeit Tl, nämlich der des gemeinsamen Quelldatensatzes Vl entsprechen. Diese bewegungskompensierten 3D-Bilddatensätze V''k(Tl) werden im Zuge der Datenkumulation38 zusammen mit dem zugehörigen Quelldatensatz Vl(Tl) voxelweise aufsummiert. Der Begriff Kumulation umfasst in diesem Sinne jede Form der Summenbildung, insbesondere einfache Summation, arithmetische Mittelwertbildung, gewichtete Mittelwertbildung, etc. Als Ergebnis der Datenkumulation38 , und somit als Endresultat des in3 dargestellten Verfahrens wird ein finaler 3D-Bilddatensatz E ausgegeben, der der Referenzzeit Tl entspricht. - Die Wirkung des vorstehend beschriebenen Verfahrens beruht insbesondere darauf, dass alle für die Datenkumulation
38 herangezogenen 3D-Bilddatensätze Vk''(Tl) und Vl(Tl) aufgrund der Bewegungskompensation hinsichtlich der bildlichen Strukturinformation übereinstimmen. Alle Bildstrukturen, die auf Strukturen im Körper des Patienten2 , insbesondere dessen Blutgefäße zurückgehen, erscheinen also in diesen 3D-Bilddatensätzen stets an dem gleichen Ort, während aber die Rekonstruktionsartefakte für alle diese 3D-Bilddatensätze verschieden sind. Bei der Datenkumulation38 werden deshalb die Rekonstruktionsartefakte ausgemittelt, während die strukturelle Bildinformation hervorgehoben wird. In dem finalen 3D-Bilddatensatz E sind deshalb auch vergleichsweise schwach kontrastierende Strukturen erkennbar, die in den vorläufigen Bilddatensätzen Vj durch die Rekonstruktionsartefakte überdeckt werden. -
- Die Summe in GLG
3 läuft aus Gründen der Vereinfachung über alle k = 1,2,...,M und umfasst damit auch formal auch einen Term für k = 1, für welchen aber trivialerweise die Identität Bll' · Vl(Tl) = Vl(Tl) gilt. - Bevorzugt wird zu jeder Referenzzeit Tl mit 1 = 1,2,....,M ein entsprechender finaler 3D-Bilddatensatz El(Tl) berechnet, so dass die komplette vierdimensionale Information über die räumliche Anordnung der Koronararterien im Körper des Patienten
2 sowie die zeitliche Veränderung dieser Anordnung während eines Herzzyklusses zur Verfügung steht. - In
7 ist eine iterativ durchgeführte Variante des vorstehend beschriebenen Verfahrens dargestellt. Dabei werden zu jedem vorläufigen 3D-Bilddatensatz Vj nur benachbarte Bilddatensätze Vj–1 und Vj+1 bewegungskompensiert und aufsummiert, wobei die resultierenden 3D-Bilddatensätze El anstelle der vorläufigen 3D-Bilddatensätze Vj in das Verfahren zurückgespeist werden. Dieser Iterationszyklus wird solange durchlaufen, bis ein vorgegebenes Abbruchkriterium erfüllt ist. Ein Iterationsschritt der in7 dargestellten Verfahrensvariante wird ausgedrückt durch die Formel mit j = 1,2,..., M. Optional sind Mischtypen der Verfahrensvarianten gemäß3 und7 vorgesehen. Insbesondere kann im Laufe aufeinander folgender Iterationsschritte die Bewegungskompensation über eine zunehmende Anzahl benachbarter Zeitfenster erstreckt werden. Des Weiteren ist optional vor gesehen, die Relativzeit t' nicht anhand eines EKG-Signals S zu bestimmen, sondern aus den Rohbildern Ri direkt abzuschätzen. Hierzu wird die zeitliche Änderung der Bildinformation der Rohbilder Ri betrachtet und hieraus durch zeitliche Spektralanalyse, Anpassung mit einer vorgegebenen Modellfunktion oder der in der Computertomographie gebräuchlichen Kymogramm-Methode der Herzzyklus des Patienten2 abgeleitet. -
- 1
- Vorrichtung
- 2
- Patient
- 3
- Aufnahmeeinheit
- 4
- Röntgenstrahler
- 5
- Röntgendetektor
- 6
- C-Bogen
- 7
- Stativ
- 8
- Isozentrische Achse
- 9
- Zentralstrahl
- 10
- Patiententisch
- 11
- Tischplatte
- 12
- EKG-Einheit
- 13
- EKG-Sensor
- 14
- Steuer- und Auswertesystem
- 15
- Datenverarbeitungsanlage
- 16
- Auswerteeinheit
- 17
- Ein-/Ausgabemittel
- 20
- Eingabeschnittstelle
- 21
- Eingabeschnittstelle
- 22
- Synchronisation
- 23
- Analysemodul
- 24
- Speichermodul
- 25
- Ein-/Ausgabeschnittstelle
- 26
- Bussystem
- 30
- Rohbildaufnahme
- 31
- Rohbildsynchronisation
- 32
- EKG-Gating
- 33
- Rückprojektion
- 34
- Datenvorbehandlung
- 35
- Bewegungsdetektion
- 36
- Datennachbehandlung
- 37
- Bewegungskompensation
- 38
- Datenkumulation
Claims (11)
- Verfahren zur Erzeugung eines dreidimensionalen Bilddatensatzes (E, El) eines bewegten Objekts (
2 ) mittels Röntgentomographie, – bei welchem eine Anzahl von zweidimensionalen röntgentomographischen Rohbildern (Ri) nach Maßgabe einer zyklischen Relativzeit (t') gruppiert werden, – bei welchem jeweils aus einander gemäß dieser Gruppierung entsprechenden Rohbildern (Ri) mindestens zwei vorläufige 3D-Bilddatensätze (Vj) erzeugt werden, – bei welchem durch Vergleich jeweils zweier als Quelldatensatz (Vl) bzw. Zieldatensatz (Vk) herangezogener vorläufiger 3D-Bilddatensätze (Vj) mindestens eine Bewegungsmatrix (Blk) abgeleitet wird, – bei welchem durch Anwendung der oder jeder Bewegungsmatrix (Blk) auf den zugehörigen Zieldatensatz (Vk) ein einer Referenzzeit (Tl) des Quelldatensatzes (Vl) entsprechender bewegungskompensierter 3D-Bilddatensatz (Vk'') erzeugt wird und – bei welchem der oder jeder bewegungskompensierte 3D-Bilddatensatz (Vk'') mit mindestens einem weiteren bewegungskompensierten 3D-Bilddatensatz (Vk'') oder einem sonstigen, der gleichen Referenzzeit (Tl) entsprechenden vorläufigen 3D-Bilddatensatz (Vl) aufsummiert wird. - Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass vor Ableitung der oder jeder Bewegungsmatrix (Blk) die hierzu herangezogen vorläufigen 3D-Datensätze (Vl, Vk) mit Mitteln der elektronischen Bildverarbeitung vorbehandelt werden.
- Verfahren nach Anspruch 2, dadurch gekennzeichnet, dass die zur Vorbehandlung der vorläufigen 3D-Bilddatensätze (Vl, Vk) eingesetzten Mittel eine Glättung des Voxelkontrastes und/oder eine Grau- bzw. Farbwertfensterung umfassen.
- Verfahren nach einem der Ansprüche 1 bis 3, dadurch gekennzeichnet, dass die oder jede Bewegungsmatrix (Blk) vor Anwendung auf den zugehörigen Zieldatensatz (Vk) regularisiert wird.
- Verfahren nach Anspruch 4, dadurch gekennzeichnet, dass die oder jede Bewegungsmatrix (Blk) zur Regularisierung in der Ortsdomäne und/ oder in der Zeitdomäne geglättet wird.
- Verfahren nach Anspruch 4 oder 5, dadurch gekennzeichnet, dass die Regularisierung der oder jeder Bewegungsmatrix (Blk) nach Maßgabe hinterlegter Information über typische Bewegungsmuster des Objekts (
2 ) erfolgt. - Verfahren nach einem der Ansprüche 1 bis 6, dadurch gekennzeichnet, dass die Relativzeit (t') aus dem Herzzyklus eines das Objekt bildenden Patienten (
2 ) abgeleitet wird. - Verfahren nach Anspruch 7, dadurch gekennzeichnet, dass der Herzzyklus des Patienten (
2 ) anhand eines EKG-Signals (S) bestimmt wird. - Verfahren nach Anspruch 7, dadurch gekennzeichnet, dass der Herzzyklus anhand der Rohbilder (Ri) abgeschätzt wird.
- Vorrichtung zur Durchführung des Verfahrens nach einem der Ansprüche 1 bis 9, mit einer um eine isozentrische Achse (
8 ) drehbar gelagerten Röntgenstrahler-Detektor-Einheit (4 ,5 ) zur Erzeugung von zweidimensionalen Rohbildern (Ri) eines Objekts (2 ), sowie mit einer Auswerteeinheit (16 ) zur Erzeugung eines dreidimensionalen Bilddatensatzes (E, El) des Objekts (2 ) anhand der Rohbilder (Ri), wobei die Auswerteeinheit (16 ) dazu ausgebildet ist, – die Rohbilder (Ri) nach Maßgabe einer zyklischen Relativzeit (t') zu gruppieren, – jeweils aus einander gemäß dieser Gruppierung entsprechenden Rohbildern (Ri) mindestens zwei vorläufige 3D-Bilddatensätze (Vj) zu erzeugen, – durch Vergleich jeweils zweier als Quelldatensatz (Vl) bzw. Zieldatensatz (Vk) herangezogener vorläufiger 3D-Bilddatensätze (Vj) mindestens eine Bewegungsmatrix (Blk) abzuleiten, – durch Anwendung der oder jeder Bewegungsmatrix (Blk) auf den zugehörigen Zieldatensatz (Vk) einen der Referenzzeit (Tl) des zugehörigen Quelldatensatzes (Vl) entsprechenden bewegungskompensierten 3D-Bilddatensatz (Vk'') zu erzeugen und – den oder jeden bewegungskompensierten 3D-Bilddatensatz (Vk'') mit mindestens einem weiteren bewegungskompensierten 3D-Bilddatensatz (Vk'') oder einem sonstigen, der gleichen Referenzzeit (Tl) entsprechenden vorläufigen 3D-Bilddatensatz (Vl) aufzusummieren. - Vorrichtung nach Anspruch 10, gekennzeichnet durch eine EKG-Einheit (
12 ) zur Messung des Herzzyklus eines das Objekt bildenden Patienten (2 ) als Maß für die Referenzzeit (t').
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102004048209A DE102004048209B3 (de) | 2004-09-30 | 2004-09-30 | Verfahren und Vorrichtung zur Erzeugung eines dreidimensionalen Bilddatensatzes eines bewegten Objekts mittels Röntgentomographie |
US11/238,552 US7899223B2 (en) | 2004-09-30 | 2005-09-29 | Method and device for generating a three-dimensional image data record of a moving object by means of X-ray tomography |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE102004048209A DE102004048209B3 (de) | 2004-09-30 | 2004-09-30 | Verfahren und Vorrichtung zur Erzeugung eines dreidimensionalen Bilddatensatzes eines bewegten Objekts mittels Röntgentomographie |
Publications (1)
Publication Number | Publication Date |
---|---|
DE102004048209B3 true DE102004048209B3 (de) | 2005-09-01 |
Family
ID=34813764
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE102004048209A Expired - Fee Related DE102004048209B3 (de) | 2004-09-30 | 2004-09-30 | Verfahren und Vorrichtung zur Erzeugung eines dreidimensionalen Bilddatensatzes eines bewegten Objekts mittels Röntgentomographie |
Country Status (2)
Country | Link |
---|---|
US (1) | US7899223B2 (de) |
DE (1) | DE102004048209B3 (de) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005045073A1 (de) * | 2005-09-21 | 2007-03-22 | Siemens Ag | Verfahren zum visuellen Unterstützen einer invasiven Untersuchung oder Behandlung des Herzens mit Hilfe eines invasiven Instruments |
WO2007072286A2 (en) * | 2005-12-20 | 2007-06-28 | Philips Intellectual Property & Standards Gmbh | Method for movement compensation of image data |
DE102006051919A1 (de) * | 2006-10-31 | 2008-05-15 | Siemens Ag | Verfahren zum Bereitstellen eines 3D-Röntgenbilddatensatzes des Herzens eines Patienten |
DE102008008611A1 (de) * | 2008-02-12 | 2009-09-03 | Siemens Aktiengesellschaft | Verfahren zur dreidimensionalen Bildrekonstruktion eines dynamisch bewegten Objektes aus Projektionsdaten einer bildgebender Einrichtung sowie zugehörige bildgebende Einrichtung |
DE102008016892A1 (de) | 2008-04-02 | 2009-10-15 | Siemens Aktiengesellschaft | Betriebsverfahren für eine Bildgebungsanlage zur zeitaufgelösten Abbildung eines sich iterativ bewegenden Untersuchungsobjekts |
DE102008033137A1 (de) | 2008-07-15 | 2010-02-04 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Einstellung einer dynamisch anpassbaren Position eines bildgebenden Systems |
DE102009021521A1 (de) | 2009-05-15 | 2010-12-02 | Siemens Aktiengesellschaft | Verfahren zur bewegungskompensierten Rekonstruktion eines dreidimensionalen Bilddatensatzes |
DE102010022791A1 (de) | 2010-06-04 | 2011-12-08 | Siemens Aktiengesellschaft | 3-D-Röntgenbildgebung von Koronargefäßen mit EKG-Gating und Bewegungskorrektur |
DE202008018167U1 (de) | 2008-07-15 | 2011-12-14 | Siemens Aktiengesellschaft | Vorrichtung zur Einstellung einer dynamisch anpassbaren Position eines bildgebenden Systems |
US8103075B2 (en) | 2007-02-23 | 2012-01-24 | Siemens Aktiengesellschaft | Method for providing extended possibilities when imaging a patient's heart |
DE102011075287A1 (de) * | 2011-05-05 | 2012-11-08 | Siemens Aktiengesellschaft | Verfahren zum Gewinnen eines 3D-Röntgenbilddatensatzes zu einem sich periodisch bewegenden Bildobjekt |
DE102011078517A1 (de) * | 2011-07-01 | 2013-01-03 | Siemens Aktiengesellschaft | Verfahren zum Bereitstellen einer Darstellung eines sich nahezu zyklisch bewegenden Objekts |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8175356B2 (en) * | 2006-11-03 | 2012-05-08 | Koninklijke Philips Electronics N.V. | Cardiac phase determination |
WO2008146273A1 (en) * | 2007-05-25 | 2008-12-04 | Yissum Research Development Company Of The Hebrew University Of Jerusalem | Method for imaging during invasive procedures performed on organs and tissues moving in a rhythmic fashion |
DE102009037478B4 (de) | 2009-08-13 | 2017-02-02 | Siemens Healthcare Gmbh | Verfahren zur 3-D-Datenerfassung mit einem Biplan-C-Bogen-System mit Biplan-Akquisitions-Multiplexing |
DE102009037477B4 (de) | 2009-08-13 | 2016-11-17 | Siemens Healthcare Gmbh | Verfahren zur 3-D-Datenerfassung mit einem Biplan-C-Bogen-System mit verbesserter zeitlicher Auflösung der Sekundärrekonstruktion |
DE102011082075B4 (de) * | 2011-09-02 | 2014-06-26 | Siemens Aktiengesellschaft | Röntgeneinrichtung mit an Tragearmen gehalterter Aufnahmeanordnung |
US9129417B2 (en) * | 2012-02-21 | 2015-09-08 | Siemens Aktiengesellschaft | Method and system for coronary artery centerline extraction |
US9526468B2 (en) * | 2014-09-09 | 2016-12-27 | General Electric Company | Multiple frame acquisition for exposure control in X-ray medical imagers |
DE102014223658A1 (de) * | 2014-11-20 | 2016-05-25 | Siemens Aktiengesellschaft | Verfahren zum Bestimmen eines Ergebnisbildes, Computerprogramm, maschinenlesbarer Datenträger sowie bildgebendes Gerät |
EP3236854B1 (de) * | 2015-02-13 | 2019-11-06 | St. Jude Medical International Holding S.à r.l. | Verfolgungsbasierte 3d-modell-erweiterung |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10129631A1 (de) * | 2001-06-20 | 2003-01-02 | Philips Corp Intellectual Pty | Verfahren zur Rekonstruktion eines hoch aufgelösten 3D-Bildes |
WO2004081877A1 (en) * | 2003-03-14 | 2004-09-23 | Koninklijke Philips Electronics N.V. | Motion-corrected three-dimensional volume imaging method |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6486671B1 (en) * | 2001-05-14 | 2002-11-26 | Ge Medical Systems Global Technologies Company Llc | MRI image quality improvement using matrix regularization |
US7372983B2 (en) * | 2001-10-16 | 2008-05-13 | Koninklijke Philips Electronics N.V. | Method for automatic branch labeling |
US6865250B2 (en) * | 2002-12-23 | 2005-03-08 | Ge Medical Systems Global Technology Company Llc | High pitch cardiac helical scan with extended reconstruction windows |
US7587074B2 (en) * | 2003-07-21 | 2009-09-08 | Paieon Inc. | Method and system for identifying optimal image within a series of images that depict a moving organ |
-
2004
- 2004-09-30 DE DE102004048209A patent/DE102004048209B3/de not_active Expired - Fee Related
-
2005
- 2005-09-29 US US11/238,552 patent/US7899223B2/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10129631A1 (de) * | 2001-06-20 | 2003-01-02 | Philips Corp Intellectual Pty | Verfahren zur Rekonstruktion eines hoch aufgelösten 3D-Bildes |
WO2004081877A1 (en) * | 2003-03-14 | 2004-09-23 | Koninklijke Philips Electronics N.V. | Motion-corrected three-dimensional volume imaging method |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102005045073B4 (de) * | 2005-09-21 | 2012-03-22 | Siemens Ag | Verfahren zum visuellen Unterstützen einer invasiven Untersuchung oder Behandlung des Herzens mit Hilfe eines invasiven Instruments |
US8583214B2 (en) | 2005-09-21 | 2013-11-12 | Siemens Aktiengesellschaft | Method for visually supporting an invasive examination or therapy of the heart with the aid of an invasive instrument |
DE102005045073A1 (de) * | 2005-09-21 | 2007-03-22 | Siemens Ag | Verfahren zum visuellen Unterstützen einer invasiven Untersuchung oder Behandlung des Herzens mit Hilfe eines invasiven Instruments |
WO2007072286A2 (en) * | 2005-12-20 | 2007-06-28 | Philips Intellectual Property & Standards Gmbh | Method for movement compensation of image data |
WO2007072286A3 (en) * | 2005-12-20 | 2007-10-18 | Philips Intellectual Property | Method for movement compensation of image data |
DE102006051919A1 (de) * | 2006-10-31 | 2008-05-15 | Siemens Ag | Verfahren zum Bereitstellen eines 3D-Röntgenbilddatensatzes des Herzens eines Patienten |
DE102006051919B4 (de) * | 2006-10-31 | 2016-09-01 | Siemens Healthcare Gmbh | Verfahren zum Bereitstellen eines 3D-Röntgenbilddatensatzes des Herzens eines Patienten |
US8027526B2 (en) | 2006-10-31 | 2011-09-27 | Siemens Aktiengesellschaft | Method for providing a 3D X-ray image dataset of a patient's heart |
US8103075B2 (en) | 2007-02-23 | 2012-01-24 | Siemens Aktiengesellschaft | Method for providing extended possibilities when imaging a patient's heart |
DE102008008611A1 (de) * | 2008-02-12 | 2009-09-03 | Siemens Aktiengesellschaft | Verfahren zur dreidimensionalen Bildrekonstruktion eines dynamisch bewegten Objektes aus Projektionsdaten einer bildgebender Einrichtung sowie zugehörige bildgebende Einrichtung |
DE102008016892A1 (de) | 2008-04-02 | 2009-10-15 | Siemens Aktiengesellschaft | Betriebsverfahren für eine Bildgebungsanlage zur zeitaufgelösten Abbildung eines sich iterativ bewegenden Untersuchungsobjekts |
US8855391B2 (en) | 2008-04-02 | 2014-10-07 | Siemens Aktiengesellschaft | Operating method for an imaging system for the time-resolved mapping of an iteratively moving examination object |
DE202008018167U1 (de) | 2008-07-15 | 2011-12-14 | Siemens Aktiengesellschaft | Vorrichtung zur Einstellung einer dynamisch anpassbaren Position eines bildgebenden Systems |
DE102008033137A1 (de) | 2008-07-15 | 2010-02-04 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Einstellung einer dynamisch anpassbaren Position eines bildgebenden Systems |
DE102009021521A1 (de) | 2009-05-15 | 2010-12-02 | Siemens Aktiengesellschaft | Verfahren zur bewegungskompensierten Rekonstruktion eines dreidimensionalen Bilddatensatzes |
DE102010022791A1 (de) | 2010-06-04 | 2011-12-08 | Siemens Aktiengesellschaft | 3-D-Röntgenbildgebung von Koronargefäßen mit EKG-Gating und Bewegungskorrektur |
US9013471B2 (en) | 2010-06-04 | 2015-04-21 | Siemens Aktiengesellschaft | 3D X-ray imaging of coronary vessels with ECG gating and motion correction |
DE102010022791B4 (de) * | 2010-06-04 | 2015-08-06 | Siemens Aktiengesellschaft | Verfahren zur dreidimensionalen Darstellung einer bewegten Struktur durch ein Rotationsangiographie-Verfahren |
DE102011075287A1 (de) * | 2011-05-05 | 2012-11-08 | Siemens Aktiengesellschaft | Verfahren zum Gewinnen eines 3D-Röntgenbilddatensatzes zu einem sich periodisch bewegenden Bildobjekt |
DE102011075287B4 (de) * | 2011-05-05 | 2016-07-21 | Siemens Healthcare Gmbh | Verfahren zum Gewinnen eines 3D-Röntgenbilddatensatzes zu einem sich periodisch bewegenden Bildobjekt |
DE102011078517A1 (de) * | 2011-07-01 | 2013-01-03 | Siemens Aktiengesellschaft | Verfahren zum Bereitstellen einer Darstellung eines sich nahezu zyklisch bewegenden Objekts |
DE102011078517B4 (de) * | 2011-07-01 | 2013-09-05 | Siemens Aktiengesellschaft | Verfahren zum Bereitstellen einer Darstellung eines sich nahezu zyklisch bewegenden Objekts |
Also Published As
Publication number | Publication date |
---|---|
US7899223B2 (en) | 2011-03-01 |
US20060067459A1 (en) | 2006-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
DE102004048209B3 (de) | Verfahren und Vorrichtung zur Erzeugung eines dreidimensionalen Bilddatensatzes eines bewegten Objekts mittels Röntgentomographie | |
DE102005027963B3 (de) | Verfahren und Vorrichtung zur Rekonstruktion eines 3D-Bilddatensatzes eines bewegten Objektes | |
DE60215964T2 (de) | Verfahren zur rekonstruktion eines 3d bildes mit hoher auflösung | |
DE102013217351B4 (de) | Bildbasierte Bewegungskompensation von Bilddaten | |
EP3185215A1 (de) | Automatisierte ermittlung von konturen auf basis einer iterativen rekonstruktion | |
DE102010019016B4 (de) | Verfahren zur Rekonstruktion von Bilddaten eines bewegten Untersuchungsobjektes aus Messdaten nebst zugehöriger Gegenstände | |
DE102012200715B4 (de) | Verfahren zur Aufnahme und Darstellung wenigstens zweier 3-D-Subtraktionsbilddatensätze sowie C-Bogen-Röntgenvorrichtung hierfür | |
DE102011083643A1 (de) | Verfahren, Rechensystem und CT-System zur Bestimmung eines Bewegungsfeldes und zur bewegungskompensierenden Rekonstruktion mit diesem Bewegungsfeld | |
DE102012205935B4 (de) | Verfahren zur Aufnahme eines vierdimensionalen Angiographie-Datensatzes | |
DE102013201136A1 (de) | Vorhersage eines voraussichtlichen Kontrastmittelverlaufs | |
DE102012205351A1 (de) | Darstellung von Blutgefäßen und Gewebe im Herzen | |
DE102011083647A1 (de) | Verfahren, Rechensystem und CT-System zur Erzeugung eines bewegungskompensierten CT-Bilddatensatzes eines sich teilweise und zyklisch bewegenden Untersuchungsobjektes | |
DE102012209410A1 (de) | Ermittlung einer patientenspezifischen Kontrastmittel-Impulsantwortfunktion | |
DE102015206362B3 (de) | Multizyklische dynamische CT-Bildgebung | |
EP3542721A1 (de) | Verfahren zur bearbeitung von parametern eines maschinenlernverfahrens sowie rekonstruktionsverfahren | |
DE60032569T2 (de) | Verfahren und Geräte zur ungleichförmigen zeitlichen Herzbildgebung | |
DE102010013360B4 (de) | Verfahren zur Rekonstruktion von Bilddaten eines zyklisch sich bewegenden Untersuchungsobjektes | |
EP3457360B1 (de) | Reduzierung von bewegungsartefakten in computertomographischen bilddaten | |
DE102009043633A1 (de) | Verbesserte Abtastung eines zyklisch bewegten Untersuchungsobjektes unter Einsatz eines Kontrastmittels im Rahmen einer Voruntersuchung mittels eines CT-Gerätes | |
DE102008016892A1 (de) | Betriebsverfahren für eine Bildgebungsanlage zur zeitaufgelösten Abbildung eines sich iterativ bewegenden Untersuchungsobjekts | |
DE102011083646A1 (de) | Verfahren, Rechensystem und CT-System zur Bestimmung eines Bewegungsfeldes und zur Erzeugung eines bewegungskompensierten CT-Bilddatensatzes eines sich teilweise bewegenden Objektes | |
DE102007045313B4 (de) | Verfahren zur getrennten dreidimensionalen Darstellung von Arterien und Venen in einem Untersuchungsobjekt | |
DE102011017710B4 (de) | Verfahren zum Bereitstellen eines 3D-Röntgenbilddatensatzes zu einem sich bewegenden Objekt mit stark absorbierendem Material | |
DE102010062975B4 (de) | Verfahren zur Erzeugung einer vierdimensionalen Darstellung eines einer periodischen Bewegung unterworfenen Zielgebiets eines Körpers | |
DE102005005919A1 (de) | Verfahren und CT-Gerät zur Erstellung von Röntgen-CT-Aufnahmen von einem schlagenden Herzen eines Patienten |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
8100 | Publication of patent without earlier publication of application | ||
8364 | No opposition during term of opposition | ||
8320 | Willingness to grant licences declared (paragraph 23) | ||
R081 | Change of applicant/patentee |
Owner name: SIEMENS HEALTHCARE GMBH, DE Free format text: FORMER OWNER: SIEMENS AKTIENGESELLSCHAFT, 80333 MUENCHEN, DE |
|
R119 | Application deemed withdrawn, or ip right lapsed, due to non-payment of renewal fee |