Refine
Has Fulltext
- yes (5)
Is part of the Bibliography
- yes (5)
Document Type
- Doctoral Thesis (5)
Language
- German (5) (remove)
Keywords
- MRT (5) (remove)
Institute
- Physikalisches Institut (5) (remove)
Sonstige beteiligte Institutionen
Die Aufgabenstellung dieser Arbeit bestand in der Entwicklung und Umsetzung von Verfahren, mit denen die Durchblutung des menschlichen Herzmuskels quantitativ bestimmt werden kann. Im Rahmen dieser Arbeit wurden dazu zwei Ansätze verfolgt, das kontrastmittelfreie Spin-Labeling Verfahren und die kontrastmittelgestützte First-Pass Messung
Die Arbeit befaßt sich mit Methoden der 23Na-NMR-Bildgebung zur Diagnose am ischämischen und infarzierten Herzmuskel. Der erste Teil beschreibt eine Methode zur lokalisierten Messung des intra- und extrazellulären Natriumgehaltes und T1. Die Methode kam in einer Studie zum Einsatz, in der intra- und extrazellulärer Natriumgehalt sowie die T1-Werte an den Tagen 1, 3 und 21 nach Infarkt gemessen wurden.Im zweiten Teil der Arbeit wird die Dynamik des 23Na bei freier Präzession im stationären Zustand (SSFP) sowohl in numerischen Simulationen als auch experimentell untersucht.
Das Ziel dieser Promotion ist der Aufbau eines dreMR Setups für einen klinischen 1,5T Scanner, das die Relaxations-Dispersions-Bildgebung ermöglicht, und die anschließende Ergründung von möglichst vielen Anwendungsfeldern von dreMR. Zu der Aufgabe gehört die Bereitstellung der zugrunde liegenden Theorie, der Bau des experimentellen Setups (Offset-Spule und Stromversorgung) sowie die Programmierung der nötigen Software. Mit dem gebauten Setup konnten zwei große Anwendungsfelder — dreMR Messungen mit und ohne Kontrastmitteln — untersucht werden.
Die MRT hat sich in den letzten Jahren zu einem wichtigen Instrument in der Diagnostik von Herzerkrankungen entwickelt. Da sie ohne ionisierende Strahlung auskommt, stellt sie vor allem auch eine nichtinvasive Alternative zu den nuklearmedizinischen Verfahren und der Computertomographie dar. Im speziellen ermöglicht die kardiale MRT die ortsaufgelöste Darstellung des Herzens mit einer Vielzahl an Kontrasten. Neben der Morphologie können damit auch zahlreiche Funktionsparameter des Herzens, wie die Ejektionsfraktion des linken Ventrikels, oder die Viabilität und Perfusion des Herzmuskels untersucht werden. Atmung und Herzbewegung stellen allerdings große Anforderungen an die MR-Herzbildgebung. Die beiden Störfaktoren limitieren den Zeitraum, der zur Bildakquisition zur Verfügung steht und erzeugen so Konflikte zwischen räumlicher Auflösung, anatomischer Abdeckung, zeitlicher Auflösung und dem Signal-zu-Rausch-Verhältnis (SNR). Ferner ergibt sich für die meisten eingesetzten Verfahren eine erhöhte Komplexität. Die Bildgebungssequenzen müssen mittels EKG an den Herzrhythmus des Patienten angepasst und die Bildakquisitionen im Atemanhaltezustand durchgeführt werden. In manchen Fällen ist sogar eine Aufspaltung der Messung in mehrere Einzelakquisitionen nötig, was wiederum die Dauer der Untersuchungen verlängert und den Patientenkomfort reduziert.
Mit technischen Entwicklungen im Bereich der Gradienten und der Empfangsspulen sowie durch den Einsatz dedizierter Bildgebungstechniken konnten in den letzten Jahren signifikante Verbesserungen erzielt und der Stellenwert der MR-Bildgebung in der Herzdiagnostik erhöht werden. Von großer Bedeutung sind dabei auch Beschleunigungsverfahren wie die Parallele Bildgebung, die eine deutliche Verkürzung der Datenakquisition ermöglichen und so den Einfluss von Atmung und Herzbewegung wirksam reduzieren. Die Beschleunigung wird dabei grundsätzlich durch eine unvollständige Datenakquisition bzw. Unterabtastung des k-Raums erzielt, welche im Zuge der Bildrekonstruktion durch Ausnutzen zusätzlich vorhandener Informationen kompensiert wird. Bei der Parallelen Bildgebung ersetzen beispielsweise mehrere um das Objekt herum angeordnete Empfangsspulen die zum Teil unvollständig durchgeführte Gradientenbasierte Ortskodierung. Die Beschleunigungsverfahren sind allerdings wegen der verringerten Datenaufnahme auch immer mit einer Reduktion des SNR verbunden.
Eine alternative Strategie zur Beschleunigung der 2D-Bildgebung mit mehreren Schichten stellt die simultane Multischichtbildgebung mit Multi-Slice Controlled Aliasing In Parallel Imaging Results In Higher Acceleration(MS-CAIPIRINHA) dar. Anders als bei der konventionellen Parallelen Bildgebung wird die Beschleunigung hier nicht durch eine reduzierte Datenaufnahme erzielt. Vielmehr werden Multiband-RF-Pulse eingesetzt, um die Spins in mehreren Schichten gleichzeitig anzuregen. Durch Anwenden schichtspezifischer RF-Phasenzyklen wird die Phase der Spins individuell in jeder Schicht moduliert, wodurch sich eine gegenseitige Verschiebung der Schichten im FOV ergibt. Die Verschiebung erleichtert die Separation der gleichzeitig angeregten Schichten mit Verfahren der Parallelen Bildgebung. Sie erlaubt außerdem eine Minimierung der bei der Rekonstruktion entstehenden Rauschverstärkung. Die Multischichtbildgebungstechnik zeichnet sich gegenüber der konventionellen Parallelen Bildgebung durch ein wesentlich höheres SNR und durch eine Bildrekonstruktion mit geringeren Rekonstruktionsfehlern aus.
In dieser Dissertation wurden verschiedene Strategien zur Anwendung von MS-CAIPIRINHA in der MRT des Herzens präsentiert sowie ihre Vorund Nachteile gegenübergestellt. Im Allgemeinen ermöglichen die vorgestellten Konzepte eine hinsichtlich des SNR sehr effiziente Erweiterung der
anatomischen Abdeckung. Unter anderem wurde eine Möglichkeit vorgestellt, mit der es uneingeschränkt gelingt, MS-CAIPIRINHA in der Bildgebung mit bSSFP-Sequenzen anzuwenden. Die Steady-State-Sequenz wird aufgrund ihres hohen intrinsischen SNR und vorteilhaften Kontrastverhaltens sehr häufig in der MRT des Herzens bei 1,5T eingesetzt. Wie auch die simultane Multischichtbildgebung erfordert sie zum Halten der Magnetisierung im stationären Zustand die Applikation eines dedizierten RF-Phasenzyklus während der Datenakquisition. Der Phasenzyklus der Sequenz ist allerdings nicht ohne Weiteres mit den Phasenzyklen der Multischichttechnik kompatibel, so dass eine Verknüpfung der beiden Verfahren bisher nur durch Aufspalten der Bildakquisition in mehrere Teilmessungen gelang. Mit dem in Kapitel 5 vorgestellten Konzept ist diese zumeist impraktikable Segmentierung nicht mehr erforderlich. Generalisierte RF-Phasenzyklen, die sowohl die Anforderungen der Sequenz, als auch die der Multischichtbildgebung erfüllen, ermöglichen eine uneingeschränkte Anwendung der Multischichttechnik in der Bildgebung mit bSSFP oder vergleichbaren Steady-State-Sequenzen. Die Multischichttechnik ist damit auch bei Untersuchungen in Echtzeit oder mit Magnetisierungspräparation – Verfahren, die unter anderem in der MR-Herzdiagnostik Verwendung finden – einsetzbar. Anhand von Echtzeit-, Cine- und First-Pass-Herzperfusionsuntersuchungen am menschlichen Herzen konnte die Anwendbarkeit des Konzepts erfolgreich demonstriert werden. Durch die Akquisition zweier Schichten in der Zeit, die normalerweise zur Bildgebung einer einzelnen Schicht benötigt wird, gelang eine Verdoppelung der anatomischen Abdeckung bei unverändert hoher Bildqualität. Bei den Herzperfusionsuntersuchungen konnten je RR-Intervall sechs Schichten akquiriert werden. Bei Echtzeit- und Cine-Messungen erlaubt das Konzept eine signifikante Reduktion der Anzahl der Atemanhaltezustände und dementsprechend eine wirksame Verkürzung der Patientenuntersuchung und eine Verbesserung des Patientenkomforts.
In Kapitel 6 wurde eine effiziente Strategie zur Anwendung der simultanen Multischichtbildgebung in der First-Pass-Herzperfusionsbildgebung bei 3T vorgestellt. Es wurde gezeigt, dass durch den Einsatz von MS-CAIPIRINHA mit Beschleunigungsfaktoren, die größer sind als die Anzahl der simultan angeregten Schichten, neben der anatomischen Abdeckung auch die räumliche Auflösung innerhalb der Bildgebungsschicht erhöht werden kann. Beide Verbesserungen sind für die MR-gestützte Diagnostik der Koronaren Herzerkrankung von Bedeutung. Während mit einer hohen räumlichen Auflösung subendokardiale und transmurale Infarktareale unterschieden werden können, erleichtert eine hohe anatomische Abdeckung die genaue Eingrenzung hypoperfundierter Bereiche. Das grundsätzliche Prinzip der vorgestellten Strategie besteht in der Kombination zweier unterschiedlicher Beschleunigungsansätze: Zur Verbesserung der anatomischen Abdeckung kommt die simultane Multischichtbildgebung zum Einsatz. Zusätzlich zur gleichzeitigen Anregung mehrerer Schichten wird der k-Raum regelmäßig unterabgetastet. Die dabei erzielte Beschleunigung wird zur Verbesserung der räumlichen Auflösung eingesetzt. Die Bildrekonstruktion erfolgt mit Verfahren der Parallelen Bildgebung. Der Vorteil des Konzepts liegt insbesondere im vollständigen Erhalt der Datenakquisitionszeit gegenüber einer unbeschleunigten Messung mit Standardabdeckung und -auflösung. Anders als bei konventionellen Beschleunigungsverfahren wirken sich lediglich die Verkleinerung der Voxelgröße sowie die Rauschverstärkung der Bildrekonstruktion SNR-reduzierend aus. Die Rauschverstärkung wird dabei, durch die gegenseitige Verschiebung der simultan angeregten Schichten im FOV, so gering wie möglich gehalten. Die Anwendbarkeit des Konzepts konnte anhand von Simulationen sowie Untersuchungen an Probanden und Herzinfarktpatienten erfolgreich demonstriert werden. Simultanes Anregen zweier Schichten und 2,5-faches Unterabtasten des k-Raums ermöglichte die Durchführung von Untersuchungen mit einer anatomischen Abdeckung von sechs bis acht Schichten je RR-Intervall und einer räumlichen Auflösung von 2,0×2,0×8,0mm3. Es konnte gezeigt werden, dass die angewandte GRAPPA-Rekonstruktion, trotz der effektiv fünffachen Beschleunigung, robust und im Wesentlichen mit geringer Rauschverstärkung durchführbar ist. Bildqualität und SNR waren für eine sektorweise Absolutquantifizierung der Myokardperfusion ausreichend, während die hohe räumliche Auflösung die Abgrenzung kleiner subendokardialer Perfusionsdefizite ermöglichte. Aufgrund seiner großen Flexibilität und recht einfachen Implementierbarkeit ist das Beschleunigungskonzept vielversprechend hinsichtlich einer Anwendung in der klinischen Routine. Die diesbezügliche Tauglichkeit ist allerdings in weiterführenden Patientenstudien noch zu evaluieren.
Alternativ zu diesem Konzept wurde in Kapitel 7 noch eine weitere, ebenfalls auf MS-CAIPIRINHA basierende Strategie für die First-Pass-Herzperfusionsbildgebung bei 3T mit großer anatomischer Abdeckung und hoher räumlicher Auflösung vorgestellt. Wie zuvor bestand die Grundidee des Konzepts darin, MS-CAIPIRINHA mit Beschleunigungsfaktoren anzuwenden, welche größer sind als die Anzahl der simultan angeregten Schichten und die Vergrößerung der anatomischen Abdeckung durch simultanes Anregen mehrerer Schichten zu realisieren. Um allerdings die bei der Bildrekonstruktion und Schichtseparation entstehende Rauschverstärkung zu minimieren, wurde zur Verbesserung der räumlichen Auflösung innerhalb der Schicht das nichtlineare Beschleunigungsverfahren Compressed Sensing zum Einsatz gebracht. Die erst in den letzten Jahren entwickelte Technik ermöglicht die exakte Rekonstruktion zufällig unterabgetasteter Daten, sofern bekannt ist, dass sich das rekonstruierte Bild in eine wohldefinierte sparse Darstellung überführen lässt. Neben der Erreichbarkeit hoher Beschleunigungsfaktoren bietet Compressed Sensing den Vorteil einer Bildrekonstruktion ohne signifikante Rauscherhöhung. Zur Einbindung des Verfahrens in das Multischichtbildgebungskonzept erfolgt die für die Verbesserung der Auflösung nötige Unterabtastung des k-Raums, zufällig und inkohärent. Zur Bildrekonstruktion sind zwei Teilschritte erforderlich. Im ersten Teilschritt werden die durch die zufällige Unterabtastung entstandenen inkohärenten Artefakte mit Compressed Sensing entfernt, im zweiten die gleichzeitig angeregten Schichten mit Verfahren der Parallelen MRT separiert. Es konnte gezeigt werden, dass die Kombination aus Compressed Sensing und MS-CAIPIRINHA eine Reduktion der inhomogenen Rauschverstärkung ermöglicht und zur Durchführung von qualitativen First-Pass-Herzperfusionsuntersuchungen mit einer Abdeckung von sechs bis acht Schichten je RR-Intervall sowie einer räumlichen Auflösung von 2,0 × 2,0 × 8,0mm3 geeignet ist. Des Weiteren konnte gezeigt werden, dass das angewandte Multischicht-Bildgebungskonzept einer Anwendung des entsprechenden Compressed-Sensing-Konzepts ohne simultane Multischichtanregung überlegen ist. Es stellte sich allerdings auch heraus, dass die rekonstruierten Bilder mit systematischen Fehlern behaftet sind, zu welchen auch ein signifikanter rekonstruktionsbedingter Verlust an zeitlicher Auflösung zählt. Dieser kann zu einer Verzerrung quantitativ bestimmter Perfusionswerte führen und verhindert so robuste quantitative Messungen der Myokardperfusion. Es ist außerdem davon auszugehen, dass auch abrupte Signalveränderungen, die bei Arrhythmien oder Bewegung auftreten, nur sehr ungenau rekonstruiert werden können. Die Systematischen Rekonstruktionsfehler konnten anhand zweier Verfahren, einer Monte-Carlo-Simulation sowie einer Analyse der lokalen Punktantworten präzise Untersucht werden. Die beiden Analysemethoden ermöglichten einerseits die genaue Bestimmung systematischer und statistischer Abweichungen der Signalamplitude und andererseits die Quantifizierung rekonstruktionsbedingter zeitlicher und räumlicher Auflösungsverluste. Dabei konnte ein Mangel an Sparsität als grundlegende Ursache der Rekonstruktionsfehler ermittelt werden. Die bei der Analyse eingesetzten Verfahren erleichtern das Verständnis von Compressed Sensing und können beispielsweise bei der Entwicklung nichtlinearer Beschleunigungskonzepte zur Bildqualitätsanalyse eingesetzt werden.
Quantifizierung myokardialer Mikrostruktur und Perfusion mittels longitudinaler NMR Relaxation
(2018)
Ziel der Arbeit war es die Quantifizierung funktioneller bzw. mikrostruktureller Parameter des Herzmuskels mit Hilfe T1-basierter Methoden zu verbessern. Diese Methoden basieren darauf, die gewünschte Information durch eine geeignete Präparation der Magnetisierung bzw. durch die Gabe von Kontrastmittel in den Zeitverlauf der longitudinalen Relaxation zu kodieren. Aus der Änderung der Relaxationszeit läßt sich dann die gewünschte Information bestimmen. Dafür sollte sowohl der Einfluß der Anatomie als auch derjenige der Meßmethodik auf die Bestimmung der longitudinalen Relaxationszeit und damit auf die Quantifizierung der Funktion bzw. Mikrostrukturparameter untersucht werden.
Speziell der Einfluß der Bildgebungssequenz führt dazu, daß nur eine scheinbare Relaxationszeit gemessen wird. Während dies keinen Einfluß auf die T1-basierte Bestimmung der untersuchten Mikrostrukturparameter hatte, ergab sich für die Perfusionsquantifizierung eine deutliche Abhängigkeit von den Parametern der verwendeten IRLL-Sequenz. Um diesen Einfluß gerecht zu werden, wurden an die Meßmethodik angepaßte Gleichungen zur Bestimmung der Perfusion gefunden mit denen die systematischen Abweichungen korrigiert werden können. Zusätzlich reduzieren die angepaßten Gleichungen die Anforderungen bezüglich der Inversionsqualität im schichtselektiven Experiment. Dies wurde in einem weiteren Projekt bei der Bestimmung der Nierenperfusion im Mausmodell ausgenutzt.
Neben der Untersuchung der Auswirkungen der Meßmethode wurde auch der Einfluß der anatomischen Besonderheiten des Blutkreislaufs am Herzen auf die Parameterquantifizierung mittels T1-basierter Methoden untersucht. Es konnte gezeigt werden, daß auf Grund der Anatomie des Herzens bei typischen Orientierungen der Bildgebungsschicht, auch bei der schichtselektiven Inversionspräparation der Magnetisierung des Herzmuskels ein Anteil des Blutpools invertiert wird. Daraus folgt, daß die vereinfachende Annahme, nach welcher bei schichtselektiver Präparation in Folge von Perfusion nur Blut mit Gleichgewichtsmagnetisierung den Herzmuskel erreicht, nicht erfüllt ist. Es konnte gezeigt werden, daß dies bei Perfusion zu einer deutlichen Unterschätzung der berechneten Perfusionswertes führt. Um mit diesem Problem umgehen zu können, wurde aufbauend auf einem vereinfachten Modell der zeitlichen Entwicklung der Blutmagnetisierung eine Korrektur für die Bestimmung der Perfusionswerte gefunden welche den Einfluß der anatomischen Besonderheiten berücksichtigt.
Das für die Perfusionskorrektur eingeführte Model prognostiziert ebenso, daß auch bei schichtselektiver Inversion die T1-basierte Bestimmung der untersuchten Mikrostrukturparameter von der Perfusion abhängig wird und eine systematische Überschätzung der quantifizierten Werte verursacht. Da die Perfusion im Kleintier deutlich höher ist als im Menschen, ist dieser Einfluß besonders in der präklinischen Forschung zu beachten. So können dort allein durch verminderte Perfusion deutliche Änderungen in den bestimmten Werten der Mikrostrukturparameter erzeugt werden, welche zu einer fehlerhaften Interpretation der Ergebnisse führen und somit ein falsches Bild für die Vorgänge im Herzmuskel suggerieren. Dabei bestätigt der Vergleich mit experimentellen Ergebnissen aus der Literatur die Vorhersagen für das Rattenmodell. Beim Menschen ist der prognostizierte Effekt deutlich kleiner. Der prognostizierte Fehler bspw. im RBV-Wert liegt in diesem Fall bei etwa 10% und wird üblicherweise in der aktuellen Forschung vernachlässigt. Inwieweit dies in er klinischen Forschung gerechtfertigt ist, muß in weiteren Untersuchungen geklärt werden.
Den untersuchten Methoden zur Bestimmung von funktionellen und mikrostrukturellen Parametern ist gemein, daß sie eine exakte Quantifizierung der longitudinalen Relaxationszeit T1 benötigen. Dabei ist im Kleintierbereich die klassische IRLL-Methode als zuverlässige Sequenz zur T1-Quantifizierung etabliert. In der klinischen Bildgebung werden auf Grund der unterschiedlichen Zeitskalen und anderer technischer Voraussetzungen andere Anforderungen an die Datenakquisition gestellt. Dabei hat in den letzten Jahren die MOLLI-Sequenz große Verbreitung gefunden. Sie ist eine Abwandlung der IRLL-Sequenz, bei der mit einer bSSFP-Bildgebungssequenz getriggert ganze Bilder während eines Herzschlages aufgenommen werden. Die MOLLI-Sequenz reagiert dabei empfindlich auf die Wartezeiten zwischen den einzelnen Transienten. Um mit diese Problematik in den Griff zu bekommen und gleichzeitig die Meßzeit verkürzen zu können wurde eine neue Methode zum Fitten der Daten entwickelt, welche die Abhängigkeit der scheinbaren Relaxationszeit von der Wartezeit zwischen den einzelnen Transienten, sowie der mittleren Herzrate fast vollständig eliminiert. Diese Methode liefert für das ganze klinisch Spektrum an erwarteten T1-Zeiten, vor und nach Kontrastmittelgabe, stabile Ergebnisse und erlaubte ein deutliche Verkürzung der Meßzeit, ohne die Anzahl der aufgenommenen Meßzeitpunkte zu reduzieren. Dies wurde in einer initialen klinischen Studie genutzt, um ECV-Werte in Patienten zu bestimmen.
Ein Nachteil der Verwendung der MOLLI-Sequenz ist, daß nur die scheinbare Relaxationszeit aus den Fit der Meßdaten bestimmt wird. Die standardmäßig genutzte Korrektur benutzt aber dem gefitteten Wert der Gleichgewichtsmagnetisierung um den wahren T1-Wert zu bestimmen. Somit ist es für die Bestimmung des T1-Wertes notwendig, die Qualität der Inversionspräparation zu kennen. Auf Basis der neuen Fitmethode wurde eine Anpassung der MOLLI-Sequenz demonstriert, welche die Bestimmung der Gleichgewichtsmagnetisierung unabhängig von der Qualität der Inversionspräparation erlaubt. Dafür verlängert sich die Meßdauer lediglich um einen Herzschlag um in geeigneter Weise ein zusätzliches Bild aufnehmen zu können.
Abschließend wurde in dieser Arbeit der Signal-Zeit-Verlauf der MOLLI-Sequenz eingehend theoretische untersucht um ein besseres Verständnis der getriggerten IRLL-Sequenzen zu entwickeln. In diesem Zusammenhang konnte eine einfache Interpretation der scheinbaren Relaxationszeit gefunden werden. Ebenso konnte erklärt werden, warum die für ungetriggerte IRLL-Sequenzen abgeleitete Korrekturgleichung auch im getriggerten Fall erstaunlich gute Ergebnisse liefert. Weiterhin konnten Fehlerquellen für die verbleibenden Abweichungen identifiziert werden, welche als Ausgangspunkt für die Ableitung verbesserter Korrekturgleichungen genutzt werden können.