DE102012216652A1 - Angiographisches Untersuchungsverfahren - Google Patents

Angiographisches Untersuchungsverfahren Download PDF

Info

Publication number
DE102012216652A1
DE102012216652A1 DE102012216652.1A DE102012216652A DE102012216652A1 DE 102012216652 A1 DE102012216652 A1 DE 102012216652A1 DE 102012216652 A DE102012216652 A DE 102012216652A DE 102012216652 A1 DE102012216652 A1 DE 102012216652A1
Authority
DE
Germany
Prior art keywords
motion map
volume images
images
corrected
image
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
Application number
DE102012216652.1A
Other languages
English (en)
Other versions
DE102012216652B4 (de
Inventor
Yiannis Kyriakou
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Healthineers Ag De
Original Assignee
Siemens AG
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Siemens AG filed Critical Siemens AG
Priority to DE102012216652.1A priority Critical patent/DE102012216652B4/de
Priority to CN201310371018.8A priority patent/CN103654823B/zh
Priority to US14/012,014 priority patent/US20140081131A1/en
Publication of DE102012216652A1 publication Critical patent/DE102012216652A1/de
Application granted granted Critical
Publication of DE102012216652B4 publication Critical patent/DE102012216652B4/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03—Computed tomography [CT]
    • A61B6/032—Transmission computed tomography [CT]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/44—Constructional features of apparatus for radiation diagnosis
    • A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4435—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure
    • A61B6/4441—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit and the detector unit being coupled by a rigid structure the rigid structure being a C-arm or U-arm
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/44—Constructional features of apparatus for radiation diagnosis
    • A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4458—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit or the detector unit being attached to robotic arms
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/44—Constructional features of apparatus for radiation diagnosis
    • A61B6/4429—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units
    • A61B6/4464—Constructional features of apparatus for radiation diagnosis related to the mounting of source units and detector units the source unit or the detector unit being mounted to ceiling
    • 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/486—Diagnostic techniques involving generating temporal series of image data
    • 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/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • A61B6/5264—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise due to motion
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5288—Devices using data or image processing specially adapted for radiation diagnosis involving retrospective matching to a physiological signal

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Veterinary Medicine (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Pulmonology (AREA)
  • Theoretical Computer Science (AREA)
  • Cardiology (AREA)
  • Automation & Control Theory (AREA)
  • Robotics (AREA)
  • Physiology (AREA)
  • Vascular Medicine (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

Die Erfindung betrifft ein angiographisches Untersuchungsverfahren eines Organs, Gefäßsystems oder anderer Körperregionen als Untersuchungsobjekt eines Patienten mittels einer 4-D-Rotationsangiographie mit folgenden Verfahrensschritten: S1) Akquisition von Projektionsbildern (24) in verschiedenen Herzphasen (c0 bis cN), S2) Rekonstruktion von 3-D-Volumenbildern (26) in den verschiedenen Herzphasen (c0 bis cN), S3) Berechnung einer Motion-Map (28, 38), S4) Bildkombination der 3-D-Volumenbilder (26) mit der Motion-Map (28, 38) zur Erzeugung von resultierenden, korrigierten 3-D-Volumenbildern (40) in den verschiedenen Herzphasen (c0 bis cN) und S5) Darstellung der resultierenden, korrigierten 3-D-Volumenbilder (40).

Description

  • Die Erfindung betrifft ein angiographisches Untersuchungsverfahren eines Organs, Gefäßsystems oder anderer Körperregionen als Untersuchungsobjekt eines Patienten mittels einer 4-D-Rotationsangiographie.
  • Ein derartiges oben genanntes angiographisches Untersuchungsverfahren lässt sich beispielsweise mit einem aus der US 7,500,784 B2 bekanntem Angiographiesystem durchführen, das anhand der 1 nachfolgend erläutert ist.
  • Die Standard-4-D-Rotationsangiographie resultiert in Rekonstruktionen von einzelnen Volumina pro Herz-Phase. Typischerweise sind diese Einzelvolumina sehr stark von Streifen-Artefakten beeinflusst, die durch die geringe Anzahl an vorhandenen Projektionen pro Herzphase entstehen.
  • Die 4-D-Rotationsangiographie, eine sogenannte 4-D-DynaCT®, kann mit mehreren Rotationen durchgeführt werden, kann aber auch mit nur einer Rotation auskommen. Bei Standard-Verfahren spielt die Anzahl von vorhandenen Projektionen pro Phase eine Rolle. In der Regel sind es bei 4-D-DynaCT® mit einer Rotation ca. 30 Projektionen pro Phase. Dadurch entstehen Streifen-Artefakte in den rekonstruierten Schichten, wie dies nachfolgend noch erläutert wird. Je weniger Projektionen genutzt werden, desto mehr Streifen-Artefakte entstehen bei der Rekonstruktion, da diese Rekonstruktionsart keine redundanten Informationen nutzt.
  • Andere aus der Literatur bekannte Verfahren arbeiten mit iterativer Rekonstruktion und Minimierungsverfahren auf Basis der Rohdaten, wie es beispielsweise in "Prior image constrained compressed sensing (PICCS): A method to accurately reconstruct dynamic CT imgages from highly undersampled projection data sets" von Guang-Hong Chen et al., veröffentlicht in Med Phys. 2008 February, Vol. 35, No. 2, Seiten 660 bis 663 beschrieben ist. Dies ist in der Regel sehr aufwändig und benötigt eine neue Rekonstruktionskette.
  • Die 1 zeigt ein als Beispiel dargestelltes biplanes Röntgensystem zur Durchführung einer 4-D-Rotationsangiographie mit zwei von je einem Ständer 1 und 1' in Form eines sechsachsigen Industrie- oder Knickarmroboters gehaltenen C-Bogen 2 und 2', an deren Enden je eine Röntgenstrahlungsquelle, beispielsweise Röntgenstrahler 3 und 3' mit Röntgenröhren und Kollimatoren, und je ein Röntgenbilddetektor 4 und 4' als Bildaufnahmeeinheit angebracht sind. Der Ständer 1 ist dabei auf dem Fußboden 5 montiert, während der zweite Ständer 1' an der Decke 6 befestigt sein kann.
  • Mittels des beispielsweise aus der US 7,500,784 B2 bekannten Knickarmroboters, welcher bevorzugt sechs Drehachsen und damit sechs Freiheitsgrade aufweist, können die C-Bogen 2 und 2' beliebig räumlich verstellt werden, zum Beispiel indem sie um ihre Drehzentren zwischen den Röntgenstrahlern 3 und 3' sowie den Röntgenbilddetektoren 4 und 4' gedreht werden. Das erfindungsgemäße angiographische Röntgensystem 1 bis 4 ist insbesondere um Drehzentren und Drehachsen in der C-Bogen-Ebene der Röntgenbilddetektoren 4 und 4' drehbar, bevorzugt um den Mittelpunkt der Röntgenbilddetektoren 4 und 4' und um den Mittelpunkt der Röntgenbilddetektoren 4 und 4' schneidende Drehachsen.
  • Der bekannte Knickarmroboter weist ein Grundgestell auf, welches beispielsweise auf dem Boden 5 oder an der Decke 6 fest montiert ist. Daran ist drehbar um eine erste Drehachse ein Karussell befestigt. Am Karussell ist schwenkbar um eine zweite Drehachse eine Roboterschwinge angebracht, an der drehbar um eine dritte Drehachse ein Roboterarm befestigt ist. Am Ende des Roboterarms ist drehbar um eine vierte Drehachse eine Roboterhand angebracht. Die Roboterhand weist ein Befestigungselement für den C-Bogen 2 oder 2' auf, welches um eine fünfte Drehachse schwenkbar und um eine senkrecht dazu verlaufende sechste Rotationsachse rotierbar ist.
  • Die Realisierung der Röntgendiagnostikeinrichtung ist nicht auf den Industrieroboter angewiesen. Es können auch übliche C-Bogen-Geräte Verwendung finden.
  • Die Röntgenbilddetektoren 4 und 4' können rechteckige oder quadratische, flache Halbleiterdetektoren sein, die vorzugsweise aus amorphem Silizium (a-Si) erstellt sind. Es können aber auch integrierende und eventuell zählende CMOS-Detektoren Anwendung finden.
  • Im Strahlengang der Röntgenstrahler 3 und 3' befindet sich eine Tischplatte 7 eines Patientenlagerungstisches 8 zur Aufnahme eines zu untersuchenden Patienten als Untersuchungsobjekt. Der Patientenlagerungstisch 8 ist mit einem Bedienpult 9 versehen. An der Röntgendiagnostikeinrichtung ist eine Systemsteuerungseinheit 10 mit einem Bildsystem 11 angeschlossen, das die Bildsignale der Röntgenbilddetektoren 4 und 4' empfängt und verarbeitet (Bedienelemente sind beispielsweise nicht dargestellt.). Die Röntgenbilder können dann auf Displays einer Monitorampel 12 betrachtet werden. Das Bildsystem 11 weist eine Vorrichtung auf, deren Funktion noch genauer beschrieben wird.
  • Anstelle des in 1 beispielsweise dargestellten Röntgensystems mit den Ständern 1 und 1' in Form des sechsachsigen Industrie- oder Knickarmroboters kann, wie in 2 vereinfacht dargestellt, das angiographische Röntgensystem auch eine normale decken- oder bodenmontierte Halterung für den C-Bogen 2 aufweisen.
  • Anstelle der beispielsweise dargestellten C-Bogen 2 und 2' kann das angiographische Röntgensystem auch getrennte decken- und/oder bodenmontierte Halterungen für die Röntgenstrahler 3 und 3' und die Röntgenbilddetektoren 4 und 4' aufweisen, die beispielsweise elektronisch starr gekoppelt sind.
  • Die Erfindung geht von der Aufgabe aus, ein angiographisches Untersuchungsverfahren der eingangs genannten Art derart auszubilden, dass eine Reduktion von Streifen-Artefakten in der herzkorrelierten 4-D-Rotationsangiographie, der sogenannten DynaCT®, unterdrückt wird.
  • Die Aufgabe wird erfindungsgemäß für ein angiographisches Untersuchungsverfahren der eingangs genannten Art durch die im Patentanspruch 1 angegebenen Merkmale gelöst. Vorteilhafte Ausbildungen sind in den abhängigen Patentansprüchen angegeben.
  • Die Aufgabe wird für ein angiographisches Untersuchungsverfahren erfindungsgemäß durch folgende Schritte gelöst:
    • S1) Akquisition von Projektionsbildern in verschiedenen Herzphasen und Positionen,
    • S2) Rekonstruktion von 3-D-Volumenbildern in den verschiedenen Herzphasen aus den Projektionsbildern,
    • S3) Berechnung einer Motion-Map aus den 3-D-Volumenbildern,
    • S4) Bildkombination der 3-D-Volumenbilder mit der Motion-Map zur Erzeugung von resultierenden, korrigierten 3-D-Volumenbildern in den verschiedenen Herzphasen und
    • S5) Darstellung der resultierenden, korrigierten 3-D-Volumenbilder.
  • Dieses erfindungsgemäße Verfahren benutzt redundante Daten, um die Streifen-Artefakte in den Bildern der herzkorrelierten 4-D-Rotationsangiographie, wie beispielsweise bei DynaCT®, zu reduzieren.
  • Die Erfindung ist nachfolgend anhand von in der Zeichnung dargestellten Ausführungsbeispielen näher erläutert. Es zeigen:
  • 1 ein bekanntes biplanes C-Bogen-Angiographiesystem mit je einem Industrieroboter als Tragvorrichtungen,
  • 2 die Verhältnisse bei einer EKG-korrelierten Akquisition während einer Rotation mit einem Rotations-Angiographiesystem gemäß 1,
  • 3 eine Serie von nach einem Standard-Verfahren der Rotationsangiographie gemäß 2 akquirierten Projektionsbildern,
  • 4 Erstellung einer Motion-Map aus rekonstruierten 3-D-Volumenbildern,
  • 5 bis 8 zeichnerische Erläuterungen zur Nachverarbeitung der gemäß 4 erstellten Motion-Map,
  • 9 zeichnerische Darstellung einer linearen Bildkombination mit linearer Interpolation und
  • 10 bis 13 zeichnerische Erläuterungen zum zeitlichen Ablauf der Nachverarbeitung und deren Ergebnisse.
  • In der 2 sind die Verhältnisse bei einer EKG-korrelierten Akquisition mit einem C-Bogen-Gerät gemäß 1 während einer Rotation veranschaulicht, die bei einer Herzrate von 90 bis 131 bpm, einer Dauer von 10s bis 15s sowie mit oder ohne Herzphasensteuerung (Pacing) durchgeführt wird. Erfolgt kein Pacing, dann wird eine bekannte manuelle Sortierung der Phasen aus dem EKG bewirkt.
  • In dieser Figur ist ein erstes EKG 13 dargestellt, das verschiedene Herzphasen c0 bis cN aufweist. Diesen Herzphasen c0 bis cN sind verschiedene Projektionswinkel θ0 bis θ0 + n·Δθ zugeordnet. So ergibt sich für ein erstes Bild 14 einer ersten Herzphase c0 ein Wert P(θ0, c0), ein erstes Bild 15 einer zweiten Herzphase P(θ0 + Δθ, c1), ein erstes Bild 16 einer dritten Herzphase P(θ0 + 2Δθ, c2) und ein erstes Bild 17 einer N-ten Herzphase P(θ0 + NΔθ, cN) P(θ0 + NΔθ, cN).
  • Dies lässt sich fortsetzen, wie dies durch den Pfeil 18 symbolisiert wird, bis man zu einem zweiten EKG 19 gelangt.
  • Diesen Herzphasen c0 bis cN sind wieder verschiedene Projektionswinkel θ0 + n·Δθ bis θ0 + (n + N)·Δθ zugeordnet. So ergibt sich für ein zweites Bild 20 einer ersten Herzphase c0 ein
  • Wert P(θ0 + nΔθ, c0), ein zweites Bild 21 einer zweiten Herzphase P(θ0 + (n + 1)Δθ, c1), ein zweites Bild 22 einer dritten Herzphase P(θ0 + (n + 2)Δθ, c2) und ein zweites Bild 23 einer N-ten Herzphase P(θ0 + (n + N)Δθ, cN).
  • In der 3 ist die nach einem Standard-Verfahren mit ca. 30 Projektionen pro Herzphase bei 120 bpm und 13 s Scanzeit erstellte Serie von Projektionsbildern 24 dargestellt, die die störenden Streifen-Artefakte 25 aufweisen. Die Indizes c0 bis cN bezeichnen die Projektionsbilder 24 der aktuellen Herzphasen.
  • Die 4 zeigt eine Folge von rekonstruierten 3-D-Volumenbildern 26, erstellt mit ca. 30 Projektionen pro Herzphase, aus denen gemäß der Formel
    Figure DE102012216652A1_0002
    eine Berechnung 27
  • einer bildbasierten Bewegungskarte oder -plan, einer sogenannten Bewegungs-Map oder Motion-Map 28 erfolgt. Die Indizes fc0 bis fcN der 3-D-Volumenbilder 26 bezeichnen das rekonstruierte 3-D-Volumen bei der entsprechenden Herzphase (c0 bis cN) und enthalten die Bildinformationen.
  • Da die Motion-Map 28 auch störende Streifen-Artefakte 25 aufweist, wird eine Nachverarbeitung der Motion-Map 28 durchgeführt, die anhand der 5 bis 8 näher erläutert wird.
  • Ein Verfahren ist die Analyse im Frequenzbereich. In 5 werden in einem 3-D-Volumenbild 26 sowie der Motion-Map 28 repräsentativ gewählte zwei Pixel 29 und 30 betrachtet, von denen das erste Pixel 29 eine große Bewegung mit niedriger Frequenz und das zweite Pixel 30 eine geringe Bewegung mit hoher Frequenz aufweisen.
  • Die 6 gibt die Signalverläufe der Pixel 29 und 30 wieder, wobei der Signalverlauf 31 des ersten Pixels 29 eine niedrigere Frequenz als der Signalverlauf 32 des zweiten Pixels 30 aufweist.
  • In der 7 sind nun die Gegebenheiten der Modulation der Herzbewegung und der Streifen-Artefakte 25 über der Ortsfrequenz u aufgetragen, wobei ein modulierter Signalverlauf 33 des ersten Pixels 29 und ein modulierter Signalverlauf 34 des zweiten Pixels 30 dargestellt sind, die eine Modulationsrichtung 35 aufweisen.
  • Die 8 stellt die Gegebenheiten nach einer Demodulation der Herzbewegung und der Streifen-Artefakte 25 aufgetragen über der Ortsfrequenz u mit einem demodulierten Signalverlauf 36 des ersten Pixels 29 und einem demodulierter Signalverlauf 37 des zweiten Pixels 30 dar.
  • Bei dem Prinzip der Modulation und Demodulation geht es im Wesentlichen darum, dass sich die Pixelwerte an einigen Stellen, beispielsweise beim zweiten Pixel 30, nur aufgrund der Streifen-Artefakte 25 quasiperiodisch ändern. Diese quasiperiodischen Änderungen der Streifen-Artefakte 25 beruhen auf dem sogenannten Windmühlen-Effekt. Sie sind Abtastartefakte als Funktion der Zeit. An anderen Stellen, beispielsweise beim ersten Pixel 29, ist die Änderung dieses Pixels 30 als Funktion der Zeit aufgrund des Windmühlen-Effekts und Herzbewegungs-Artefakten zurückzuführen. Diese Art der Änderung gilt es zu erkennen und diese mit Filtern, z.B. Demodulation, selektive Weichzeichnung zu bearbeiten.
  • Die Prinzipien der Modulation und Demodulation sind allgemein aus der Signaltheorie oder Signalverarbeitung bekannt; dabei können Fourier-Analyse oder Bandfilterung eingesetzt werden. Die Modulation ist durch die Aufnahme selber gegeben; die Demodulation wird dazu genutzt, um das "Träger"-Signal von dem "Echten"-Signal zu trennen. Bei dieser hier vorliegenden bestimmten Aufnahmeform geht das relativ einfach, da die Windmühlen-Artefakte eine ziemlich definierte Frequenz aufweisen, die lediglich von der Aufnahmegeometrie abhängt und demnach einfach vorher zu berechnen ist.
  • Als weitere Verfahren zur Nachverarbeitung der Motion-Map 28 sind morphologische Operationen wie beispielsweise Erosion und/oder Dilatation der Motion-Map 28 anwendbar.
  • Auch kann das Verfahren des Sub-Sampling und der Interpolation, beispielsweise bilinear oder spline, der Motion-Map 28 zur Nachverarbeitung angewandt werden.
  • Als Resultat der Nachverarbeitung der Motion-Map 28 durch eines dieser Verfahren erhält man eine korrigierte Motion-Map, die nahezu frei von Streifen-Artefakten 25 ist.
  • Ein in 9 dargestelltes Beispiel einer Bildkombination ist eine lineare Kombination mit linearer Interpolation Es sind jedoch auch andere Kombinationsarten wie beispielsweise polynomische oder quadratische Bildkombinationen möglich. Auch Bildkombinationen mit einem Faltungsoperator sind denkbar.
  • Anhand der 9 wird nun eine der möglichen Bildkombinationen erklärt, die sich allgemein aus folgender Gleichung ergibt: F(x, y, z, cn) = f(x, y, z, cn)·MM(x, y, z) + f(x, y, z)·(1 – MM(x, y, z)) wobei cn die jeweilige Herzphase c0 bis cN darstellt.
  • Die Pixel der rekonstruierten 3-D-Volumenbilder 26 f(x, y, z, cn) werden mit den Pixeln der korrigierten Motion-Map 38 MM(x, y, z) multipliziert. Dazu wird das Produkt aus Eins minus korrigierter Motion-Map 38 MM(x, y, z) und dem Mittelwertbild 39 f(x, y, z) über alle Phasenbilder addiert. Als Ergebnis F(x, y, z, cn) erhält man resultierende, korrigierte 3-D-Volumenbilder 40.
  • Diese Multiplikation stellt den einfachsten Fall einer Bildkombination dar, bei der eine pixel- oder voxelweise Multiplikation (Gewichtung) der zwei Bilder (oder Volumen) immer pro Phase durchgeführt wird, wobei die Motion-Map nach der Nachverarbeitung konstant bleibt.
  • Anders formuliert würde das Ergebnis für das Beispiel der ersten Herzphase c0 folgendermaßen aussehen: Fc0(x, y, z) = fc0(x, y, z)·MM(x, y, z) + f(x, y, z)·(1 – MM(x, y, z))
  • Dies ist beispielhaft für eine lineare Interpolation so dargestellt. Im Falle einer nicht-linearen Kombination müsste eine entsprechende Funktion f(MM(x, y, z)) definiert werden, z.B. polynomisch. Im vorliegenden Falle geht es hauptsächlich darum, die Einzelvolumina entsprechend der Motion-Map zu gewichten.
  • Das Ergebnis der Nachverarbeitung kann auch anhand der 10 bis 13 näher erläutert und symbolisch gezeigt werden, die die zeitliche Reihenfolge der Bildentstehungen wiedergeben. Der Ausgangspunkt ist die Bildserie "vor Motion-Map – Nachverarbeitung" der rekonstruierten 3-D-Volumenbilder 26. Daraus wird die Motion-Map 28 berechnet. Diese Motion-Map 28 wird dann anhand der in den 5 bis 8 beschriebenen Bearbeitung in eine "Motion-Map – Nachverarbeitung", der korrigierten Motion-Map 38 korrigiert. Zum Schluss werden gemäß oben genannter Gleichung die resultierenden, korrigierten 3-D-Volumenbilder 40 "nach Motion-Map – Nachverarbeitung" berechnet.
  • Das oben vorgeschlagene Verfahren arbeitet auf Basis der rekonstruierten Schichten, der 3-D-Volumenbilder 26.
  • Eine Akquisitionsart ist eine Rotation mit guter Winkelabtastung, beispielsweise einer Abtastzeit von 13 s, 0,5° Winkelinkrement und 2 × 2 Binning. Daraus resultieren in etwa 380 Projektionen über alle Phasen. Vorhandene, redundante Informationen werden ausgenutzt, da sich nur einige der Voxel im Bild ändern. Die Änderung der Voxel wird durch die Motion-Map 28 pro Schicht ausgerechnet. Die Motion-Map 28 gibt den Gehalt der Bewegung oder Änderung der Voxel-Werte über die Zeit wieder. Ein Voxel weist eine unterschiedliche Bewegungsfunktion, d.h. Änderungsfunktion oder Gradient, im Herzen auf, als wenn es sich in einem anderen Körperteil befände.
  • Die Motion-Map 28 ist im ersten Schritt auch durch Streifen-Artefakte 25 beeinflusst. Um diese zu reduzieren, werden drei Nachverarbeitungsmethoden vorgeschlagen, um Änderungen durch Streifen-Artefakte 25 und Änderungen durch reine Herzbewegung zu trennen. Daraus resultiert eine Reduzierung der Streifen-Artefakte 25 in der Motion-Map 28.
  • Die Motion-Map 28 wird als Kombinationsgewicht zwischen der Rekonstruktion einer einzelnen Phase (z.B. c0) und dem Mittelwertbild aus allen Phasen genutzt. Es wird dabei angenommen, dass die Voxel-Werte in der Motion-Map 28 mit einem kleinen Wert weniger zur Herzbewegung beitragen.
  • Die Bildkombination kann durch lineare Interpolation gemacht werden, aber auch andere Kombinationsarten sind möglich.
  • Die resultierenden korrigierten 3-D-Volumenbilder 40 weisen deutlich weniger Streifen-Artefakte 25 auf.
  • Das erfindungsgemäße Verfahren ist für Monoplan- und Biplan-Systeme anwendbar. Entgegen vielen anderen bekannten Verfahren ist es ein reines bild-basiertes Verfahren. Es sind weder Rohdaten, noch Geometrie- oder andere Informationen erforderlich.
  • Durch das erfindungsgemäße Verfahren werden Streifen-Artefakte 25 aus 4-D-Rotationsangiographien, sogenannte 4-D-DynaCT®-Bilder, nahezu vollständig mit begrenzten Raum- und Zeit-Auflösungsverlust eliminiert.
  • Die Generierung und Nachverarbeitung der Motion-Map 28 verringert weiterhin störende Streifen-Artefakte 25.
  • Das erfindungsgemäße Verfahren ist auch anwendbar für andere Protokolle mit Änderungen in der Zeitrichtung, beispielsweise Perfusion.
  • Für die Berechnungen wird die vorhandene Rekonstruktionskette effektiv ausgenutzt.
  • ZITATE ENTHALTEN IN DER BESCHREIBUNG
  • Diese Liste der vom Anmelder aufgeführten Dokumente wurde automatisiert erzeugt und ist ausschließlich zur besseren Information des Lesers aufgenommen. Die Liste ist nicht Bestandteil der deutschen Patent- bzw. Gebrauchsmusteranmeldung. Das DPMA übernimmt keinerlei Haftung für etwaige Fehler oder Auslassungen.
  • Zitierte Patentliteratur
    • US 7500784 B2 [0002, 0007]
  • Zitierte Nicht-Patentliteratur
    • "Prior image constrained compressed sensing (PICCS): A method to accurately reconstruct dynamic CT imgages from highly undersampled projection data sets" von Guang-Hong Chen et al., veröffentlicht in Med Phys. 2008 February, Vol. 35, No. 2, Seiten 660 bis 663 [0005]

Claims (5)

  1. Angiographisches Untersuchungsverfahren eines Organs, Gefäßsystems oder anderer Körperregionen als Untersuchungsobjekt eines Patienten mittels einer 4-D-Rotationsangiographie gekennzeichnet durch folgende Schritte: S1) Akquisition von Projektionsbildern (24) in verschiedenen Herzphasen (c0 bis cN) und Positionen, S2) Rekonstruktion von 3-D-Volumenbildern (26) in den verschiedenen Herzphasen (c0 bis cN) aus den Projektionsbildern (24), S3) Berechnung einer Motion-Map (28, 38) aus den 3-D-Volumenbildern (26), S4) Bildkombination der 3-D-Volumenbilder (26) mit der Motion-Map (28, 38) zur Erzeugung von resultierenden, korrigierten 3-D-Volumenbildern (40) in den verschiedenen Herzphasen (c0 bis cN) und S5) Darstellung der resultierenden, korrigierten 3-D-Volumenbilder (40).
  2. Angiographisches Untersuchungsverfahren nach Anspruch 1, dadurch gekennzeichnet, dass aus den 3-D-Volumenbildern (26) ein Mittelwertbild (39) f(x, y, z) über alle Herzphasen gebildet wird, das in der Bildkombination gemäß Verfahrensschritt S4) mit einbezogen wird.
  3. Angiographisches Untersuchungsverfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass sich die resultierenden, korrigierten 3-D-Volumenbilder (40) gemäß folgender Gleichung berechnen: F(x, y, z, cn) = f(x, y, z, cn)·MM(x, y, z) + f(x, y, z)·(1 – MM(x, y, z)), wobei – cn die jeweilige Herzphase c0 bis cN, – f(x, y, z, cn) rekonstruierte 3-D-Volumenbilder (26), – MM(x, y, z) eine Motion-Map (28, 38), – f(x, y, z) ein Mittelwertbild (39) über alle Phasenbilder und – F(x, y, z, cn) resultierende, korrigierte 3-D-Volumenbilder (40) darstellen.
  4. Angiographisches Untersuchungsverfahren nach Anspruch 3, dadurch gekennzeichnet, dass die Motion-Map (28, 38) eine nachverarbeitete, korrigierte Motion-Map (38) ist.
  5. Angiographisches Untersuchungsverfahren nach einem der Ansprüche 1 bis 4, dadurch gekennzeichnet, dass sich die Motion-Map (28) gemäß Verfahrensschritt S3) folgenderweise berechnet:
    Figure DE102012216652A1_0003
DE102012216652.1A 2012-09-18 2012-09-18 Angiographisches Untersuchungsverfahren Active DE102012216652B4 (de)

Priority Applications (3)

Application Number Priority Date Filing Date Title
DE102012216652.1A DE102012216652B4 (de) 2012-09-18 2012-09-18 Angiographisches Untersuchungsverfahren
CN201310371018.8A CN103654823B (zh) 2012-09-18 2013-08-23 血管造影检查方法
US14/012,014 US20140081131A1 (en) 2012-09-18 2013-08-28 Angiographic examination method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102012216652.1A DE102012216652B4 (de) 2012-09-18 2012-09-18 Angiographisches Untersuchungsverfahren

Publications (2)

Publication Number Publication Date
DE102012216652A1 true DE102012216652A1 (de) 2014-03-20
DE102012216652B4 DE102012216652B4 (de) 2023-01-26

Family

ID=50181664

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102012216652.1A Active DE102012216652B4 (de) 2012-09-18 2012-09-18 Angiographisches Untersuchungsverfahren

Country Status (3)

Country Link
US (1) US20140081131A1 (de)
CN (1) CN103654823B (de)
DE (1) DE102012216652B4 (de)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9427200B2 (en) 2014-03-21 2016-08-30 Siemens Aktiengesellschaft Determination of physiological cardiac parameters as a function of the heart rate
US10561383B2 (en) 2014-12-18 2020-02-18 Koninklijke Philips N.V. Imaging system for imaging an elongated region of interest of an object

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016083068A1 (en) * 2014-11-25 2016-06-02 Koninklijke Philips N.V. Digital subtraction angiography
DE102016219709B3 (de) * 2016-10-11 2018-03-01 Siemens Healthcare Gmbh Verfahren zur Ermittlung eines Perfusionsdatensatzes, sowie Röntgenvorrichtung, Computerprogramm und elektronisch lesbarer Datenträger
CN109767468B (zh) * 2019-01-16 2021-04-20 上海长征医院 脏器体积检测方法及装置
US11510642B2 (en) * 2019-02-06 2022-11-29 William E. Butler Spatiotemporal reconstruction in higher dimensions of a moving vascular pulse wave from a plurality of lower dimensional angiographic projections

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102007029731A1 (de) * 2007-06-27 2009-01-02 Siemens Ag Verfahren zur automatischen Bestimmung einer optimalen Herzphase für eine Cardio-CT-Rekonstruktion und Cardio-CT-Rekonstruktion
US7500784B2 (en) 2005-03-18 2009-03-10 Siemens Aktiengesellschaft X-ray device

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5696848A (en) * 1995-03-09 1997-12-09 Eastman Kodak Company System for creating a high resolution image from a sequence of lower resolution motion images
US6539074B1 (en) * 2000-08-25 2003-03-25 General Electric Company Reconstruction of multislice tomographic images from four-dimensional data
DK1494579T3 (da) * 2002-04-02 2011-11-14 Yeda Res & Dev Karakterisering af bevægelige objekter i en stationær baggrund
US6904118B2 (en) * 2002-07-23 2005-06-07 General Electric Company Method and apparatus for generating a density map using dual-energy CT
DE602005012583D1 (de) * 2004-06-18 2009-03-19 Koninkl Philips Electronics Nv Reduktion von artefakten
US7782998B2 (en) 2004-12-21 2010-08-24 General Electric Company Method and apparatus for correcting motion in image reconstruction
EP1869643B1 (de) * 2005-04-07 2014-05-07 Koninklijke Philips N.V. Bildverarbeitungsgerät und -verfahren zur blutfluss-bildgebung
DE102008016892A1 (de) * 2008-04-02 2009-10-15 Siemens Aktiengesellschaft Betriebsverfahren für eine Bildgebungsanlage zur zeitaufgelösten Abbildung eines sich iterativ bewegenden Untersuchungsobjekts
CN103349556B (zh) * 2009-01-21 2015-09-23 皇家飞利浦电子股份有限公司 用于大视场成像及运动伪影的探测和补偿的方法和装置
EP2513864A1 (de) * 2009-12-15 2012-10-24 Koninklijke Philips Electronics N.V. Bewegungskompensation mit gewebedichteretention
US20110148928A1 (en) * 2009-12-17 2011-06-23 General Electric Company System and method to correct motion in gated-pet images using non-rigid registration
US8761467B2 (en) * 2010-10-04 2014-06-24 General Electric Company Method and apparatus for assessing motion correction
US8600132B2 (en) * 2011-05-03 2013-12-03 General Electric Company Method and apparatus for motion correcting medical images
US8897527B2 (en) * 2011-06-07 2014-11-25 Varian Medical Systems, Inc. Motion-blurred imaging enhancement method and system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7500784B2 (en) 2005-03-18 2009-03-10 Siemens Aktiengesellschaft X-ray device
DE102007029731A1 (de) * 2007-06-27 2009-01-02 Siemens Ag Verfahren zur automatischen Bestimmung einer optimalen Herzphase für eine Cardio-CT-Rekonstruktion und Cardio-CT-Rekonstruktion

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Prior image constrained compressed sensing (PICCS): A method to accurately reconstruct dynamic CT imgages from highly undersampled projection data sets" von Guang-Hong Chen et al., veröffentlicht in Med Phys. 2008 February, Vol. 35, No. 2, Seiten 660 bis 663
Schirra C O; u.a: Improvement of cardiac CT reconstruction using local motion vector fields. In: Computerized Medical Imaging and Graphics, 2009, 122 - 130. *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9427200B2 (en) 2014-03-21 2016-08-30 Siemens Aktiengesellschaft Determination of physiological cardiac parameters as a function of the heart rate
US10561383B2 (en) 2014-12-18 2020-02-18 Koninklijke Philips N.V. Imaging system for imaging an elongated region of interest of an object

Also Published As

Publication number Publication date
CN103654823B (zh) 2016-04-20
CN103654823A (zh) 2014-03-26
US20140081131A1 (en) 2014-03-20
DE102012216652B4 (de) 2023-01-26

Similar Documents

Publication Publication Date Title
DE69831742T2 (de) System zur rekonstruktion bei kegelstrahltomographie
DE102010027227B4 (de) Verfahren und Computertomographiegerät zur Durchführung einer angiographischen Untersuchung
DE102007061935A1 (de) Verfahren zur Qualitätssteigerung von computertomographischen Aufnahmeserien durch Bildverarbeitung und CT-System mit Recheneinheit
DE102012204977B4 (de) Verfahren zur iterativen Bildrekonstruktion für Bi-Modale CT-Daten sowie entsprechende Recheneinheit, Computertomographiesystem, Computerprogramm und Datenträger
DE102012216652B4 (de) Angiographisches Untersuchungsverfahren
DE102006041033A1 (de) Verfahren zur Rekonstruktion eines dreidimensionalen Bildvolumens und Röntgengeräte
DE102010022791A1 (de) 3-D-Röntgenbildgebung von Koronargefäßen mit EKG-Gating und Bewegungskorrektur
DE102011006188B4 (de) Verfahren und Computertomographie-System zur Erstellung tomographischer Bilddarstellung mit mindestens zwei Strahler-Detektor-Systemen
DE19647435A1 (de) Systeme, Verfahren und Vorrichtungen zur Rekonstruktion von Bildern in eine Wendelabtastung verwendenden Computer-Tomographie-Systemen
DE102012214351A1 (de) Angiographisches Untersuchungsverfahren eines Gefäßsystems
DE102006006451A1 (de) Verfahren zum Betrieb einer Röntgendiagnostikeinrichtung zur Erzeugung hochaufgelöster Subtraktionsangiographie-Bilder
DE102009047867B4 (de) Verfahren und Vorrichtung zur Korrektur von trunkierten Projektionsdaten
DE102012220028A1 (de) Angiographisches Untersuchungsverfahren
DE102010013360B4 (de) Verfahren zur Rekonstruktion von Bilddaten eines zyklisch sich bewegenden Untersuchungsobjektes
DE102007061934A1 (de) Verfahren zur Qualitätssteigerung von computertomographischen Aufnahmeserien durch Projektionsdatenverarbeitung und CT-System mit Recheneinheit
DE102008016892A1 (de) Betriebsverfahren für eine Bildgebungsanlage zur zeitaufgelösten Abbildung eines sich iterativ bewegenden Untersuchungsobjekts
DE102007045313B4 (de) Verfahren zur getrennten dreidimensionalen Darstellung von Arterien und Venen in einem Untersuchungsobjekt
DE102012205245B4 (de) Angiographisches Untersuchungsverfahren eines Patienten zur 3-D-Rotationsangiographie
DE102010062975B4 (de) Verfahren zur Erzeugung einer vierdimensionalen Darstellung eines einer periodischen Bewegung unterworfenen Zielgebiets eines Körpers
DE102004017478B4 (de) Vorrichtung für die Gewinnung von Strukturdaten eines sich bewegenden Objekts
DE102009007236A1 (de) CT-Bildrekonstruktion eines sich bewegenden Untersuchungsobjektes
DE102019204765B3 (de) Verfahren zur Ermittlung eines dreidimensionalen Tomosynthesedatensatzes, Röntgeneinrichtung, Computerprogramm und elektronisch lesbarer Datenträger
DE102012202648B3 (de) Angiographisches Untersuchungsverfahren zur Abtastung von kleinen beweglichen Objekten
DE102011005554A1 (de) Verfahren zur Reduzierung von Bewegungsartefakten bei der Dual-Energy-Computertomographie
DE102009037242A1 (de) Roadmap-Verfahren zur überlagerten Darstellung von Bildern

Legal Events

Date Code Title Description
R012 Request for examination validly filed
R016 Response to examination communication
R081 Change of applicant/patentee

Owner name: SIEMENS HEALTHCARE GMBH, DE

Free format text: FORMER OWNER: SIEMENS AKTIENGESELLSCHAFT, 80333 MUENCHEN, DE

R016 Response to examination communication
R018 Grant decision by examination section/examining division
R020 Patent grant now final
R081 Change of applicant/patentee

Owner name: SIEMENS HEALTHINEERS AG, DE

Free format text: FORMER OWNER: SIEMENS HEALTHCARE GMBH, MUENCHEN, DE