Refine
Has Fulltext
- yes (36) (remove)
Is part of the Bibliography
- yes (36)
Year of publication
Document Type
- Doctoral Thesis (36)
Language
- German (36) (remove)
Keywords
- Kernspintomografie (36) (remove)
Institute
- Physikalisches Institut (19)
- Institut für diagnostische und interventionelle Radiologie (Institut für Röntgendiagnostik) (5)
- Medizinische Klinik und Poliklinik I (5)
- Graduate School of Life Sciences (4)
- Institut für diagnostische und interventionelle Neuroradiologie (ehem. Abteilung für Neuroradiologie) (2)
- Institut für Molekulare Infektionsbiologie (1)
- Klinik und Poliklinik für Allgemein-, Viszeral-, Gefäß- und Kinderchirurgie (Chirurgische Klinik I) (1)
- Klinik und Poliklinik für Anästhesiologie (ab 2004) (1)
- Klinik und Polikliniken für Zahn-, Mund- und Kieferkrankheiten (1)
- Lehrstuhl für Orthopädie (1)
Sonstige beteiligte Institutionen
Die MRT des Herzens wird aufgrund hoher Reproduzierbarkeit und geringer Variabilität als Referenzstandard für die Bestimmung der kardialen Funktion betrachtet. Auch in der präklinischen Forschung bietet die MRT eine ausgezeichnete Charakterisierung der kardialen Funktion und ermöglicht eine exzellente Analyse modellierter Krankheitsbilder. In beiden Fällen besteht jedoch weiterhin Optimierungsbedarf. Die klinische Herz-MRT stellt ein aufwendiges Verfahren mit relativ langer Messzeit dar und ist dadurch mit hohen Untersuchungskosten verbunden. In der präklinischen Kleintierbildgebung müssen zum Erreichen der notwendigen höheren Orts- und Zeitauflösung ebenfalls lange Aufnahmezeiten in Kauf genommen werden. Um die kardiale MRT dort routinemäßig in großen Studienkollektiven anwenden zu können, ist eine schnellere Bildgebung essentiell. Neben einer Verbesserung der Tomographen-Hardware und der Optimierung von Bildgebungssequenzen standen im letzten Jahrzehnt vermehrt informationstheoretische Ansätze zur Beschleunigung der MR-Datenakquisition im Fokus der Entwicklung. Während zu Beginn des Jahrtausends die Parallele Bildgebung (PI) einen Forschungsschwerpunkt repräsentierte, spielte sich in den letzten fünf Jahren vermehrt die von Donoho und Candès eingeführte Compressed Sensing (CS) Theorie in den Vordergrund. Diese ermöglicht eine Signalrekonstruktion aus unvollständig gemessenen Koeffizienten einer linearen Messung (z.B. Fouriermessung) unter Ausnutzung der Sparsität des Signals in einer beliebigen Transformationsbasis. Da sich die MRT hervorragend für den Einsatz von CS eignet, wurde die Technik in der Forschung bereits vielfach angewendet. Die zur Rekonstruktion unterabgetasteter Aufnahmen nötigen CS-Algorithmen haben jedoch eine signifikante Veränderung des Bildgebungsprozesses der MRT zur Folge. Konnte dieser zuvor in guter Näherung als linear und stationär betrachtet werden, so repräsentiert die CS-Rekonstruktion eine nichtlineare und nichtstationäre Transformation. Objektinformation wird nicht mehr ortsunabhängig und proportional zur Intensität in die Abbildung transportiert. Das Bild ist viel mehr das Ergebnis eines Optimierungsprozesses, der sowohl die Konsistenz gegenüber der unterabgetasteten Messung als auch die Sparsität des Signals maximiert. Der erste Teil dieser Dissertation beschreibt eine Methode, die eine objektive Einschätzung der Bildqualität CS-rekonstruierter MR-Bilder ermöglicht. Die CS-Beschleunigung verspricht eine Verkürzung der Messzeit ohne Verlust an Bildqualität, wobei letztere bisher größtenteils qualitativ bzw. quantitativ nur unzureichend beurteilt wurde. Konnte der Bildgebungsprozess der klassischen MRT (linear und stationär) durch die Bestimmung einer Punktspreizfunktion (PSF) robust und effektiv validiert und optimiert werden, erlauben die CS-Algorithmen aufgrund ihres nichtlinearen und nichtstationären Verhaltens ohne Weiteres keine äquivalente Analyse. Um dennoch eine entsprechende Evaluierung des CS-Bildgebungsprozesses zu ermöglichen, wurde die Anwendung einer lokalen Punktspreizfunktion (LPSF) für den in der Folge verwendeten Iterative Soft Thresholding Algorithmus untersucht. Die LPSF berücksichtigt die Ortsabhängigkeit der CS-Rekonstruktion und muss daher für jeden Ort (Pixel) eines Bildes bestimmt werden. Darüber hinaus wurde die LPSF im linearen Bereich der CS-Transformation ermittelt. Dazu wurde das zu bewertende Bild nach Anwenden einer kleinen lokalen Störung rekonstruiert. Die Breite des Hauptmaximums der LPSF wurde schließlich verwendet, um ortsaufgelöste Auflösungsstudien durchzuführen. Es wurde sowohl der Einfluss typischer Unterabtastschemata für CS als auch der Einsatz diskreter Gradienten zur Sparsifizierung eines Phantombildes untersucht. Anschließend wurde die Prozedur zur Bestimmung der räumlichen und zeitlichen Auflösung in der Herzbildgebung getestet. In allen Beispielen ermöglichte das vorgeschlagene Verfahren eine solide und objektive Analyse der Bildauflösung CS-rekonstruierter Aufnahmen. Wurde zuvor meist ausschließlich auf Vergleiche mit einer vollständig abgetasteten Referenz zur Qualitätsbeurteilung zurückgegriffen, so stellt die vorgestellte Auflösungsbestimmung einen Schritt in Richtung einer standardisierten Bildanalyse bei der Verwendung der Beschleunigung mittels CS dar. Die Analyse der Abtastmuster zeigte, dass auch bei der Anwendung von CS die Berücksichtigung der nominell höchsten Frequenzen k_max unerlässlich ist. Frühere Publikationen schlagen Abtastfolgen mit einer teils starken Gewichtung der Messpunkte zum k-Raum-Zentrum hin vor. Die Ergebnisse der vorliegenden Arbeit relativieren ein derartiges Vorgehen, da zumindest bei den durchgeführten Untersuchungen ein Auflösungsverlust bei analoger Vorgehensweise zu verzeichnen war. Ebenso zeigten sich dynamische Aufnahmen, die unter Verwendung des x-f-Raums als sparse Basis rekonstruiert wurden, durchaus anfällig für zeitliches Blurring. Dieses resultiert aus der Unterdrückung hoher zeitlicher Frequenzen und konnte durch die ortsaufgelösten Auflösungskarten sichtbar gemacht werden. Neben der Auflösung ist für eine umfassende Analyse der Bildqualität auch die Untersuchung potentieller Aliasing-Artefakte sowie des Signal-zu-Rausch-Verhältnisses (SNR) notwendig. Während Aliasing mit Hilfe der Einträge der LPSF außerhalb des Hauptmaximums untersucht werden kann, wurde in Kap. 5 eine Modifikation der Multi-Replika-Methode von Robson et al. zur Rauschanalyse bei Verwendung nichtlinearer Algorithmen vorgestellt. Unter Einbeziehung aller genannten Qualitätsparameter ist eine robuste Bewertung der Bildqualität auch bei einer Verwendung von CS möglich. Die differenzierte Evaluierung ebnet den Weg hin zu einem objektiven Vergleich neuer Entwicklungen mit bisherigen Standard-Techniken und kann dadurch den Einzug von CS in die klinische Anwendung vorantreiben. Nach den theoretischen Betrachtungen der Bildqualität behandelt die Dissertation die erstmalige Anwendung von CS zur Beschleunigung der funktionellen Herzdiagnostik in der präklinischen MR-Kleintierbildgebung. Diese Studien wurden in Zusammenarbeit mit der British Heart Foundation Experimental Magnetic Resonance Unit (BMRU) der University of Oxford durchgeführt. Die Algorithmen für eine Beschleunigung mittels der CS-Theorie wurden anhand der dort am 9,4T Tomographen gemessenen (unterabgetasteten) Datensätze entwickelt und optimiert. Zunächst wurde eine Beschleunigung ausschließlich mittels CS untersucht. Dazu wurde die segmentierte, EKG- und Atemgetriggerte kartesische Cine-Aufnahme in Phasenkodierrichtung unterabgetastet und mittels CS rekonstruiert. Die sparse Darstellung wurde durch Ermitteln zeitlicher Differenzbilder für jede Herzphase erhalten. Durch Variation der Abtastmuster in der zeitlichen Dimension konnte ein vollständig abgetastetes zeitliches Mittelbild bestimmt werden, das anschließend von jedem einzelnen Herzphasenbild subtrahiert wurde. In einer Validierungsphase wurden an der Maus vollständig aufgenommene Cine-Akquisitionen retrospektiv unterabgetastet, um die maximal mögliche Beschleunigung mittels CS zu ermitteln. Es wurden u.a. funktionelle Herz-Parameter für jede Gruppe des jeweiligen Beschleunigungsfaktors bestimmt und mittels einer statistischen Analyse verglichen. Die Gesamtheit aller Ergebnisse zeigte die Möglichkeit einer dreifachen Beschleunigung ohne eine Degradierung der Genauigkeit der Methode auf. Die ermittelte Maximalbeschleunigung wurde in einer unterabgetastet gemessenen Bilderserie mit anschließender CS-Rekonstruktion validiert. Die Abtastschemata wurden dazu mit Hilfe der Transformations-Punktspreizfunktion weiter optimiert. In einer Erweiterung der Studie wurde zum Zweck einer noch höheren Beschleunigung die CS-Technik mit der PI kombiniert. Erneut fand eine Unterabtastung der Phasenkodierrichtung einer kartesischen Trajektorie statt. Die Messungen erfolgten mit einer 8-Kanal-Mäusespule an einem 9,4T Tomographen. Um das Potential beider Beschleunigungstechniken auszunutzen, wurden die Methoden CS und PI in serieller Weise implementiert. Für die PI-Beschleunigung wurde der vollständig abgetastete k-Raum zunächst gleichmäßig unterabgetastet. Auf dem resultierenden Untergitter wurde zusätzlich eine Unterabtastung nach Pseudo-Zufallszahlen durchgeführt, um eine Beschleunigung mittels CS zu ermöglichen. Die entwickelte Rekonstruktion erfolgte ebenfalls seriell. Zunächst wurde mittels CS das äquidistante Untergitter rekonstruiert, um anschließend mittels GRAPPA die noch fehlenden Daten zu berechnen. Um eine zusätzliche Messung zur Kalibrierung der GRAPPA-Faktoren zu umgehen, wurde das äquidistant unterabgetastete Untergitter von Herzphase zu Herzphase um je einen Phasenkodierschritt weitergeschoben. Dieses Vorgehen erlaubt die Ermittlung eines vollständig abgetasteten k-Raums mit einer geringeren zeitlichen Auflösung, der die notwendige Bestimmung der Wichtungsfaktoren ermöglicht. Folgende Kombinationen von Beschleunigungsfaktoren wurden mittels retrospektiver Unterabtastung eines vollständig aufgenommenen Datensatzes untersucht: R_CS x R_PI = 2 x 2, 2 x 3, 3 x 2 und 3 x 3. Die Analyse des Bildrauschens, des systematischen Fehlers und der Auflösung führte zu dem Schluss, dass eine sechsfache Beschleunigung mit Hilfe der hybriden Rekonstruktionstechnik möglich ist. Während mit steigender CS-Beschleunigung der systematische Fehler leicht anstieg, führte ein höherer PI-Beschleunigungsfaktor zu einer leichten Verstärkung des statistischen Fehlers. Der statistische Fehler zeigte jedoch ebenfalls eine Verringerung bei steigender Beschleunigung mittels CS. Die Fehler waren allerdings stets auf einem Niveau, das durchaus auch Beschleunigungen bis R_CS x R_PI =3 x 3 zulässt. Die LPSF-Analyse zeigte einen Verlust der räumlichen Auflösung von ca. 50 % bei R=6 sowie einen mittleren Verlust von 64 % bei R=9. Offensichtlich ging die ebenfalls beobachtete Minimierung des Bildrauschens durch den CS-Algorithmus im Falle der relativ stark verrauschten Kleintieraufnahmen zu Lasten der Bildauflösung. Die mit zunehmender Beschleunigung stärker geblurrten Grenzen zwischen Blutpool und Myokardgewebe erschweren die Segmentierung und stellen eine mögliche Fehlerquelle dar. Unter Beachtung aller Ergebnisse ist eine sechsfache Beschleunigung (R_CS x R_PI = 2 x 3, 3 x 2) vertretbar. Die Hinzunahme der PI ermöglicht somit im Vergleich zur alleinigen Verwendung von CS eine weitere Beschleunigung um einen Faktor von zwei. Zusammenfassend ermöglicht der Einsatz von CS in der präklinischen funktionellen Herzbildgebung am Kleintier eine deutliche Reduktion der Messzeit. Bereits ohne Vorhandensein von Mehrkanalspulen kann die notwendige Datenmenge ohne signifikante Beeinflussung der Messergebnisse auf ein Drittel reduziert werden. Ist der Einsatz von Spulenarrays möglich, kann die mit PI mögliche dreifache Beschleunigung um einen weiteren Faktor zwei mittels CS auf R=6 erweitert werden. Dementsprechend kann CS einen wesentlichen Beitrag dazu leisten, dass das Potential Herz-MRT am Kleintier in großen Studienkollektiven effektiver abgerufen werden kann. Im letzten Teil der Arbeit wurde eine Technik für die funktionelle klinische MR-Herzbildgebung entwickelt. Hier wurde eine Beschleunigung mittels CS verwendet, um die Aufnahme des gesamten Herzens innerhalb eines Atemstillstandes des Patienten zu ermöglichen. Bei der derzeitigen Standardmethode werden üblicherweise 10-15 2D-Schichten des Herzens akquiriert, wobei jede einzelne Aufnahme einen Atemstillstand des Patienten erfordert. Für die notwendige Beschleunigung wurde eine unterabgetastete 3D-Trajektorie verwendet. Durch Phasenkodierung einer Richtung sowie radiale Projektionen in den beiden anderen Dimensionen konnte eine effiziente Aufnahme unterhalb des Nyquist-Kriteriums erreicht werden. Die Sparsifizierung erfolgte, wie bereits in der beschriebenen präklinischen Anwendung, durch die Subtraktion eines zeitlichen Mittelbildes. In einer Simulation anhand eines retrospektiv unterabgetasteten Datensatzes konnte die theoretische Funktionalität der Rekonstruktionstechnik bei einer Beschleunigung bezüglich der Nyquist-Abtastung von R ~ 10 validiert werden. Die Unterschiede zum vollständig abgetasteten Datensatz waren vernachlässigbar klein, so dass die vorgeschlagene Abtastfolge am Tomographen implementiert wurde. Mit dieser Sequenz wurde anschließend eine funktionelle Bilderserie an einem gesunden Probanden mit vollständiger Herzabdeckung innerhalb eines Atemstopps aufgenommen. Fehlende Daten wurden analog zur Simulation mit Hilfe des vorgeschlagenen Algorithmus rekonstruiert. Im Vergleich zur Simulation ergaben sich aufgrund des Schichtprofils der 3D-Slab-Anregung zusätzliche Aliasing-Artefakte in den äußeren Partitionen. Die für radiale Aufnahmen typischen Streifenartefakte waren im rekonstruierten Bild, wenn auch mit sehr geringer Amplitude, noch erkennbar. Davon abgesehen wurde die Dynamik jedoch über das gesamte Herz hinweg gut dargestellt. Der hohe Kontrast zwischen Myokard und Blutpool bescheinigt den Bildern eine hervorragende Eignung für die Bestimmung funktioneller Herzparameter mittels einer Segmentierung. Zusammengefasst erlaubt die entwickelte Methode aufgrund der drastischen Reduktion der notwendigen Atemstopps des Patienten einen deutlich erhöhten Patientenkomfort sowie einen schnelleren Durchsatz aufgrund der verkürzten Messzeit.
Ziel dieser Arbeit war es, das gesamte Lungenvolumen in hoher dreidimensionaler Auflösung mittels der MRT darzustellen. Um trotz der niedrigen Protonendichte der Lunge und der geforderten hohen Auflösung ausreichend Signal für eine verlässliche Diagnostik zu erhalten, sind Aufnahmezeiten von einigen Minuten nötig. Um die Untersuchung für den Patienten angenehmer zu gestalten oder auf Grund der eingeschränkten Fähigkeit eines Atemstopps überhaupt erst zu ermöglichen, war eine Anforderung, die Aufnahmen in freier Atmung durchzuführen. Dadurch entstehen allerdings Bewegungsartefakte, die die Diagnostik stark beeinträchtigen
und daher möglichst vermieden werden müssen. Für eine Bewegungskompensation
der Daten muss die auftretende Atembewegung detektiert werden. Die Bewegungsdetektion
kann durch externe Messgeräte (Atemgurt oder Spirometer) oder durch eine
zusätzliche Anregungen erfolgen (konventionelle Navigatoren) erfolgen. Nachteile
dieser Methoden bestehen darin, dass die Bewegung während der Atmung nicht
direkt verfolgt wird, dass elektronische Messgeräte in die Nähe des Tomographen
gebracht werden und das die Patienten zusätzlich vorbereitet und eingeschränkt
werden. Des Weiteren erfordert eine zusätzliche Anregung extra Messzeit und kann
unter Umständen die Magnetisierung auf unterwünschte Weise beeinflussen.
Um die angesprochenen Schwierigkeiten der Bewegungsdetektion zu umgehen,
wurden in dieser Arbeit innerhalb einer Anregung einer 3d FLASH-Sequenz sowohl
Bilddaten- als auch Navigatordaten aufgenommen. Als Navigator diente dabei das
nach der Rephasierung aller bildgebenden Gradienten entstehende Signal (DC Signal).
Das DC Signal entspricht dabei der Summe aller Signale, die mit einem bestimmten
Spulenelement detektiert werden können. Bewegt sich beispielsweise die Leber
bedingt durch die Atmung in den Sensitivitätsbereich eines Spulenelementes, wird
ein stärkeres DC Signal detektiert werden. Je nach Positionierung auf dem Körper
kann so die Atembewegung mit einzelnen räumlich lokalisierten Spulenelementen
nachverfolgt werden. Am DC Signalverlauf des für die Bewegungskorrektur ausgewählten
Spulenelementes sind dann periodische Signalschwankungen zu erkennen.
Zusätzlich können aus dem Verlauf Expirations- von Inspirationszuständen unterschieden
werden, da sich Endexpirationszustände im Regelfall durch eine längere
Verweildauer auszeichnen.
Grundsätzlich kann das DC Signal vor oder nach der eigentlichen Datenaufnahme
innerhalb einer Anregung aufgenommen werden. Auf Grund der kurzen Relaxationszeit
T∗2 des Lungengewebes fällt das Signal nach der RF Anregung sehr schnell ab. Um
möglichst viel Signal zu erhalten sollten, wie in dieser Arbeit gezeigt wurde, innerhalb einer Anregung zuerst die Bilddaten und danach die Navigatordaten aufgenommen
werden. Dieser Ansatz führt zu einer Verkürzung der Echozeit TE um 0.3 ms und
damit zu einem SNR Gewinn von etwa 20 %. Gleichzeitig ist das verbleibende Signal
nach der Datenakquisition und Rephasierung der bildgebenden Gradienten noch
ausreichend um die Atembewegung zu erfassen und somit eine Bewegungskorrektur
der Daten (Navigation) zu ermöglichen.
Um eine retrospektive Bewegungskorrektur durchführen zu können, müssen Akzeptanzbedingungen
(Schwellenwerte) für die Datenauswahl festgelegt werden. Bei
der Wahl des Schwellenwertes ist darauf zu achten, dass weder zu wenige noch zu
viele Daten akzeptiert werden. Akzeptiert man sehr wenige Daten, zeichnen sich die
Rekonstruktionen durch einen scharfen Übergang zwischen Lunge und Diaphragma
aus, da man sehr wenig Bewegung in den Rekonstruktionen erlaubt. Gleichzeitig
erhöht sich allerdings das Risiko, dass nach der Navigation Linien fehlen. Dies führt
zu Einfaltungsartefakten, die in Form von gestörten Bildintensitäten in den Rekonstruktionen
zu sehen sind und die diagnostische Aussagekraft einschränken. Um
Einfaltungsartefakte zu vermeiden sollte der Schwellenwert so gewählt werden, dass
nach der Datenauswahl keine Linien fehlen. Aus dieser Anforderung lässt sich ein
maximaler Schwellenwert ableiten. Akzeptiert man dagegen sehr viele Daten, zeichnen
sich die Rekonstruktionen durch erhöhtes Signal und das vermehrte Auftreten
von Bewegungsartefakten aus. In diesem Fall müsste der Arzt entscheiden, ob Bewegungsartefakte
die Diagnostik zu stark beeinflussen. Wählt man den Schwellenwert
so, dass weder Linien fehlen noch zu viel Bewegung erlaubt wird, erhält man Rekonstruktionen
die sich durch einen scharfen Diaphragmaübergang auszeichnen und in
denen noch kleinste Gefäße auch in der Nähe des Diaphragmas deutlich zu erkennen
sind. Hierfür haben sich Schwellenwerte, die zu einer Datenakzeptanz von ca. 40 %
führen als günstig erwiesen.
Um Einfaltungsartefakte auf Grund der retrospektiven Datenauswahl zu verhindern,
muss das Bildgebungsvolumen mehrfach abgetastet werden. Dadurch wird
gewährleistet, dass für die letztendliche Rekonstruktion ausreichend Daten zur Verfügung
stehen, wobei mehrfach akzeptierte Daten gemittelt werden. Dies spielt auf
Grund der niedrigen Protonendichte der Lunge eine wesentliche Rolle in der Rekonstruktion
hochaufgelöster Lungendatensätze. Weiterhin führt das Mitteln von
mehrfach akzeptierten Daten zu einer Unterdrückung der sogenannten Ghost Artefakte,
was am Beispiel der Herzbewegung in der Arbeit gezeigt wird.
Da die Messungen unter freier Atmung durchgeführt werden und keine zusätzlichen
externen Messgeräte angeschlossen werden müssen, stellte die Untersuchung
für die Patienten in dieser Arbeit kein Problem dar. Im ersten Teil dieser wurde Arbeit
gezeigt, dass sich mit Hilfe des DC Signales als Navigator und einer retrospektiven
Datenauswahl das gesamte Lungenvolumen in hoher dreidimensionaler Auflösung
von beispielsweise 1.6 x 1.6 x 4 mm3 innerhalb von 13 min. darstellen lässt. Die Anwendbarkeit der vorgestellten Methode zur Bewegungskorrektur wurde neben
Probanden auch an Patienten demonstriert.
Da wie bereits beschrieben das Bildgebungsvolumen mehrfach abgetastet werden
muss, wiederholt sich auch die Abfolge der für die Bildgebung verantwortlichen
Gradienten periodisch. Da sich der Atemzyklus aber auch periodisch wiederholt,
kann es zu Korrelationen zwischen der Atmung und den wiederholten Messungen
kommen. Dies führt dazu, dass auch nach vielen wiederholten Messungen immer
noch größere Bereiche fehlender Linien im k-Raum bleiben, was zu Artefakten in
den Rekonstruktionen führt. Dies konnte im Falle der konventionellen Bewegungskorrektur
in den Gatingmasken, die die Verteilung und Häufigkeit der einzelnen
akzeptierten Phasenkodierschritte im k-Raum zeigen, beobachtet werden.
Da eine vorsätzliche Unterbrechung der Atemperiodizität (der Patient wird dazu
angehalten, seine Atemfrequenz während der Messung absichtlich zu variieren) zur
Vermeidung der angesprochenen Korrelationen nicht in Frage kommt, musste die
Periodizität in der Datenaufnahme unterbrochen werden. In dieser Arbeit wurde
dies durch eine quasizufällige Auswahl von Phasen- und Partitionskodiergradienten
erreicht, da Quasizufallszahlen so generiert werden, dass sie unabhängig von ihrer
Anzahl einen Raum möglichst gleichförmig ausfüllen. Die quasizufällige Datenaufnahme
führt deshalb dazu, das sowohl akzeptierte als auch fehlende Linien nach der
Bewegungskorrektur homogen im k-Raum verteilt auftreten.
Vergleicht man das auftreten von Ghosting zeichnen sich die quasizufälligen Rekonstruktionen
im Vergleich zur konventionellen Datenaufnahme durch eine verbesserte
Reduktion von Ghost Artefakten aus. Dies ist auf die homogene Verteilung mehrfach
akzeptierter Linien im k-Raum zurückzuführen. Die homogenere Verteilung von
fehlenden Linien im k-Raum führt weiterhin zu einer wesentlich stabileren Rekonstruktion
fehlender Linien mit parallelen MRT-Verfahren (z.B. iterativem Grappa).
Dies wird umso deutlicher je höher der Anteil fehlender Linien im k-Raum wird. Im
Falle der konventionellen Datenaufnahme werden die zusammenhängenden Bereiche
fehlender Linien immer größer, was eine erfolgreiche Rekonstruktion mit iterativem
Grappa unmöglich macht. Im Falle der quasizufälligen Datenaufnahme dagegen
können auch Datensätze in denen 40% der Linien fehlen einfaltungsartefaktfrei
rekonstruiert werden.
Im weiteren Verlauf der Arbeit wurde gezeigt, wie die Stabilität der iterativen Grappa
Rekonstruktion im Falle der quasizufälligen Datenaufnahme für eine erhebliche
Reduktion der gesamten Messzeit genutzt werden kann. So ist in einer Messzeit von
nur 74s die Rekonstruktion eines artefaktfreien und bewegungskorrigierten dreidimensionalen
Datensatzes der menschlichen Lunge mit einer Auflösung von 2 x 2 x
5 mm3 möglich. Des Weiteren erlaubt die quasizufällige Datenaufnahme in Kombination
mit iterativem Grappa die Rekonstruktion von Datensätzen unterschiedlicher
Atemphasen von Inspiration bis Expiration (4D Bildgebung). Nach einer Messzeit
von 15min. wurden 19 unterschiedliche Atemzustände rekonstruiert, wobei sich der Anteil der fehlenden Linien zwischen 0 und 20 % lag. Im Falle der konventionellen
Datenaufnahme wäre eine wesentlich längere Messzeit nötig gewesen, um ähnliche
Ergebnisse zu erhalten.
Zum Schluss soll noch ein Ausblick über mögliche Weiterentwicklungen und Anwendungsmöglichkeiten,
die sich aus den Erkenntnissen dieser Arbeit ergeben haben,
gegeben werden. So könnte das quasizufällige Aufnahmeschema um eine Dichtegewichtung
erweitert werden. Hierbei würde der zentrale k-Raum Bereich etwas
häufiger als die peripheren Bereiche akquiriert werden. Dadurch sollte die iterative
Grappa Rekonstruktion noch stabiler funktionieren und Ghost Artefakte besser reduziert
werden. Die Verteilung der Linien sollte allerdings nicht zu inhomogen werden,
um größere Lücken im k-Raum zu vermeiden.
Darüber hinaus könnte die vorgestellte Methode der Bewegungskompensation
auch für die Untersuchung anderer Organe oder Körperteile verwendet werden.
Voraussetzung wäre lediglich das Vorhandensein dezidierter Spulenanordnungen,
mit denen die Bewegung nachverfolgt werden kann. So ist beispielsweise eine dynamische
Bildgebung des frei und aktiv bewegten Knies möglich, wobei zwischen
Beugung und Streckung durch die erste Ableitung des zentralen k-Raum Signales
unterschieden werden kann. Dies kann zusätzliche Diagnoseinformationen liefern
oder für Verlaufskontrollen nach Operationen benutzt werden [15].
Eine Weiterentwicklung mit hohem klinischen Potential könnte die Kombination
der in dieser Arbeit vorgestellten retrospektiven Bewegungskorrektur mit einer Multi-
Gradienten-Echo Sequenz darstellen. Hierzu musste die bestehende Sequenz lediglich
um eine mehrfache Abfolge von Auslesegradienten innerhalb einer Anregung erweitert
werden. Dies ermöglicht eine bewegungskorrigierte voxelweise Bestimmung der
transversalen Relaxationszeit T∗2 in hoher räumlicher Auflösung. Unter zusätzlicher
Sauerstoffgabe kann es zu einer Veränderung von T∗2 kommen, die auf den sogenannten
BOLD Effekt (Blood Oxygen Level Dependent) zurückzuführen ist. Aus dieser
Änderung könnten Rückschlüsse auf hypoxische Tumorareale gezogen werden. Da
diese eine erhöhte Strahlenresistenz aufweisen, könnte auf diese Bereiche innerhalb
des Tumors eine erhöhte Strahlendosis appliziert und so möglicherweise Behandlungsmisserfolge
reduziert werden. Gleichzeitig kann durch die 4D Bildgebung eine
mögliche Tumorbewegung durch die Atmung erfasst und diese Information ebenfalls
in der Bestrahlungsplanung benutzt werden. Die Lungen MRT könnte somit um eine
hochaufgelöste dreidimensionale funktionelle Bildgebung erweitert werden.
In dieser Arbeit wurden für spezielle Anwendungen an klinischen MR-Geräten optimierte Phased-Array-Spulen entwickelt. Das Ziel war, durch die Verwendung neuer Spulen entweder neue Anwendungsgebiete für klinische MR-Geräte zu eröffnen oder bei bestehenden Applikationen die Diagnosemöglichkeiten durch eine Kombination von höherem SNR und kleineren g-Faktoren im Vergleich zu bestehenden Spulen zu verbessern. In Kapitel 3 wurde untersucht, ob es durch den Einsatz neu entwickelter, dedizierter Kleintierspulen sinnvoll möglich ist, Untersuchungen an Kleintieren an klinischen MR-Geräten mit einer Feldstärke von 1,5T durchzuführen. Der Einsatz dieser Spulen verspricht dem klinischen Anwender Studien an Kleintieren durchführen zu können, bei denen er den gleichen Kontrast wie bei einer humanen Anwendung erhält und gleichzeitig Kontrastmittel sowie Sequenzen, die klinisch erprobt sind, einzusetzen. Durch die gewählten geometrischen Abmessungen der Spulen ist es möglich, Zubehör von dedizierten Tier-MR-Geräten, wie z. B. Tierliegen oder EKG- bzw. Atemtriggereinheiten, zu verwenden. Durch Vorversuche an für Ratten dimensionierten Spulen wurden grundlegende Zusammenhänge zwischen verwendetem Entkopplungsmechanismus und SNR bzw. Beschleunigungsfähigkeit erarbeitet. Für Ratten wurde gezeigt, dass in akzeptablen Messzeiten von unter fünf Minuten MR-Messungen des Abdomens in sehr guter Bildqualität möglich sind. Ebenfalls gezeigt wurde die Möglichkeit durch den Einsatz von paralleler Bildgebung sowie Kontrastmitteln hochaufgelöste Angiographien durchzuführen. Es stellte sich heraus, dass bei 1,5T dedizierte Mäusespulen bei Raumtemperatur von den SNR-Eigenschaften am Limit des sinnvoll Machbaren sind. Trotzdem war es möglich, auch für Mäuse ein 4-Kanal-Phased-Array zu entwickeln und den Einsatz bei kontrastmittelunterstützten Applikationen zu demonstrieren. Insgesamt wurde gezeigt, dass durch den Einsatz von speziellen, angepassten Kleintierspulen auch Tieruntersuchungen an klinischen MR-Geräten mit niedriger Feldstärke durchführbar sind. Obwohl sich die Bestimmung der Herzfunktion an MR-Geräten im klinischen Alltag zum Goldstandard entwickelt hat, ist die MR-Messung durch lange Atemanhaltezyklen für einen Herzpatienten sehr mühsam. In Kapitel 4 wurde deswegen die Entwicklung einer 32-Kanal-Herzspule beschrieben, welche den Komfort für Patienten deutlich erhöhen kann. Schon mit einem ersten Prototypen für 3T war es möglich, erstmals Echtzeitbildgebung mit leicht reduzierter zeitlicher Auflösung durchzuführen und damit auf das Atemanhalten komplett zu verzichten. Dies ermöglicht den Zugang neuer Patientengruppen, z. B. mit Arrythmien, zu MR-Untersuchungen. Durch eine weitere Optimierung des Designs wurde das SNR sowie das Beschleunigungsvermögen signifikant gesteigert. Bei einem Beschleunigungsfaktor R = 5 in einer Richtung erhält man z. B. gemittelt über das gesamte Herz ein ca. 60 % gesteigertes SNR zu dem Prototypen. Die Kombination dieser Spule zusammen mit neuentwicklelten Methoden wie z. B. Compressed- Sensing stellt es in Aussicht, die Herzfunktion zukünftig in der klinischen Routine in Echtzeit quantifizieren zu können. In Kapitel 5 wurde die Entwicklung einer optimierten Brustspulen für 3T beschrieben. Bei Vorversuchen bei 1,5T wurden Vergleiche zwischen der Standardspule der Firma Siemens Healthcare und einem 16-Kanal-Prototypen durchgeführt. Trotz größerem Spulenvolumen zeigt die Neuentwicklung sowohl hinsichtlich SNR als auch paralleler Bildgebungseigenschaften eine signifikante Verbesserung gegenüber der Standardspule. Durch die Einhaltung aller Kriterien für Medizinprodukte kann diese Spule auch für den klinischen Einsatz verwendet werden. Mit den verbesserten Eigenschaften ist es beispielsweise möglich, bei gleicher Messdauer eine höhere Auflösung zu erreichen. Aufgrund des intrinsischen SNR-Vorteils der 3 T-Spule gegenüber der 1,5 T-Spule ist es dort sogar möglich, bei höheren Beschleunigungsfaktoren klinisch verwertbare Schnittbilder zu erzeugen. Zusammenfassend wurden für alle drei Applikationen NMR-Empfangsspulen entwickelt, die im Vergleich zu den bisher verfügbaren Spulen, hinsichtlich SNR und Beschleunigungsvermögen optimiert sind und dem Anwender neue Möglichkeiten bieten.
Herzkreislauferkrankungen stellen die häufigsten Todesursachen in den Industrienationen dar. Die Entwicklung nichtinvasiver Bildgebungstechniken mit Hilfe der Magnetresonanz-Tomografie (MRT) ist daher von großer Bedeutung, um diese Erkrankungen frühzeitig zu erkennen und um die Entstehungsmechanismen zu erforschen. In den letzten Jahren erwiesen sich dabei genetisch modifzierte Mausmodelle als sehr wertvoll, da sich durch diese neue Bildgebungsmethoden entwickeln lassen und sich der Krankheitsverlauf im Zeitraffer beobachten lässt.
Ein große Herausforderung der murinen MRT-Bildgebung sind die die hohen Herzraten und die schnelle Atmung. Diese erfordern eine Synchronisation der Messung mit dem Herzschlag und der Atmung des Tieres mit Hilfe von Herz- und Atemsignalen. Konventionelle Bildgebungstechniken verwenden zur Synchronisation mit dem Herzschlag EKG Sonden, diese sind jedoch insbesondere bei hohen Feldstärken (>3 T) sehr störanfällig. In dieser Arbeit wurden daher neue Bildgebungsmethoden entwickelt, die keine externen Herz- und Atemsonden benötigen, sondern das MRT-Signal selbst zur Bewegungssynychronisation verwenden. Mit Hilfe dieser Technik gelang die Entwicklung neuer Methoden zur Flussbildgebung und der 3D-Bildgebung, mit denen sich das arterielle System der Maus qualitativ und quantitativ erfassen lässt, sowie einer neuen Methode zur Quantisierung der longitudinalen Relaxationszeit T1 im murinen Herzen. Die in dieser Arbeit entwickelten Methoden ermöglichen robustere Messungen des Herzkreislaufsystems. Im letzten Kapitel konnte darüber hinaus gezeigt werden dass sich die entwickelten Bildgebungstechniken in der Maus auch auf die humane Bildgebung übertragen lassen.
Das Signal-zu-Rausch-Verhältnis (SNR) stellt bei modernen Bildgebungstechniken in der Magnetresonanz-Tomographie heutzutage oftmals die entscheidende Limitation dar. Eine Verbesserung durch Modifikation der Hardware ist kostspielig und führt meistens zu einer Verstärkung anderer Probleme, wie zum Beispiel erhöhte Energiedeposition ins Gewebe. Im Gegensatz dazu ist Dichtegewichtung eine Methode, die eine SNR-Erhöhung durch Modifikation der Aufnahmetechnik ermöglicht. In der MR-Bildgebung erfolgt oftmals eine retrospektive Filterung des aufgenommenen Signalverlaufs, beispielsweise zur Artefaktreduktion. Damit einhergehend findet eine Veränderung der Modulationstransferfunktion (MTF) bzw. ihrer Fouriertransformierten, der räumlichen Antwortfunktion (SRF), statt. Optimales SNR wird nach dem Matched Filter-Theorem erzielt, wenn die nachträgliche Filterung dem aufgenommenen Signalverlauf proportional ist. Dies steht dem Ziel der Artefaktreduktion entgegen. Bei Dichtegewichtung steht durch nicht-kartesische Abtastung des k-Raums mit der k-Raum-Dichte ein zusätzlicher Freiheitsgrad zur Verfügung. Dieser ermöglicht es, im Falle eines konstanten Signalverlaufs eine gewünschte MTF ohne Filterung zu erreichen. Bei veränderlichem Signalverlauf kann ein SNR Matched Filter angewendet werden, dessen negative Einflüsse auf die MTF durch Dichtegewichtung kompensiert werden. Somit ermöglicht Dichtegewichtung eine vorgegebene MTF und gleichzeitig ein optimales SNR. In der vorliegenden Arbeit wurde Dichtegewichtung erstmals bei den schnellen Multi-Echo-Sequenzen Turbo-Spin-Echo und Echoplanar-Bildgebung (EPI) angewendet. Im Gegensatz zu bisherigen Implementierungen muss hier der Signalabfall durch T2- bzw. T2*-Relaxation berücksichtigt werden. Dies führt dazu, dass eine prospektiv berechnete dichtegewichtete Verteilung nur bei einer Relaxationszeit optimal ist. Bei Geweben mit abweichenden Relaxationszeiten können sich wie auch bei den kartesischen Varianten dieser Sequenzen Änderungen an SRF und SNR ergeben. Bei dichtegewichteter Turbo-Spin-Echo-Bildgebung des Gehirns konnte mit den gewählten Sequenzparametern ein SNR-Vorteil von 43 % gegenüber der kartesischen Variante erzielt werden. Die Akquisition wurde dabei auf die T2-Relaxationszeit von weißer Substanz optimiert. Da die meisten Gewebe im Gehirn eine ähnliche Relaxationszeit aufweisen, blieb der visuelle Gesamteindruck identisch zur kartesischen Bildgebung. Der SNR-Gewinn konnte in der dichtegewichteten Implementierung zur Messzeithalbierung genutzt werden. Dichtegewichtete EPI weist eine hohe Anfälligkeit für geometrische Verzerrungen, welche durch Inhomogenitäten des Hauptmagnetfeldes verursacht werden, auf. Die Verzerrungen konnten erfolgreich mit einer Conjugate Phase-Methode korrigiert werden. Dazu muss die räumliche Verteilung der Feldinhomogenitäten bekannt sein. Dazu ist zusätzlich zur eigentlichen EPI-Aufnahme die zeitaufwendige Aufnahme einer sogenannten Fieldmap erforderlich. Im Rahmen dieser Arbeit konnte eine Methode entwickelt werden, welche die zur Erlangung einer Fieldmap notwendige Aufnahmedauer auf wenige Sekunden reduziert. Bei dieser Art der Fieldmap-Aufnahme müssen jedoch durch Atmung hervorgerufene Effekte auf die Bildphase berücksichtigt werden. Die Fieldmap-Genauigkeit kann durch Aufnahme unter Atempause, Mittelung oder retrospektiver Phasenkorrektur erhöht werden. Für die gewählten EPI-Sequenzparameter wurde mit Dichtegewichtung gegenüber der kartesischen Variante ein SNR-Gewinn von 14 % erzielt. Anhand einer funktionellen MRT (fMRI)-Fingertapping-Studie konnte demonstriert werden, dass die SNR-Steigerung auch zu einer signifikant erhöhten Aktivierungsdetektion in Teilen der Hirnareale führt, die bei der Fingerbewegung involviert sind. Die Verwendung von zusätzlicher EPI-Phasenkorrektur und iterativer Optimierung der dichtegewichteten k-Raum-Abtastung führt zu weiteren Verbesserungen der dichtegewichteten Bildgebung mit Multi-Echo-Sequenzen.
Vorliegende Untersuchung am Universitätsklinikum Würzburg sowie die Befragung von Anästhesisten/Anästhesistinnen im Raum der 3 DACH-Länder zeigen, dass bildgebende Verfahren bei Säuglingen mit einer niedrigen Rate an Komplikationen, zumeist in medikamentöser Sedierung mit Propofol, durchgeführt werden. Wie international üblich ist im Säuglingsalter die Magnetresonanztomographie das bildgebende Verfahren der Wahl und wird, mit überzeugender Häufigkeit, erfolgreich durchgeführt.
Die Untersuchung am Universitätsklinikum Würzburg legt nahe, dass männliche Säuglinge häufiger eine Bildgebung benötigen und häufiger höheren ASA-Kategorie zugeschrieben werden. Dabei scheinen sie auch häufiger Komplikationen zu erleben und bedürfen daher besonderer Aufmerksamkeit. Eine eventuelle Alternative zur Sedierung kann dabei die „feed-and-sleep“ Methode darstellen. In unserer Umfrage konnten wir erheben, dass diese Methode bisher wenig verbreitet ist, obwohl in diesem Zusammenhang eventuell Abläufe und Prozesszeiten strukturiert und optimiert werden können, da beispielsweise die Nachüberwachung entfällt. Vorstellbar wäre beispielsweise, mehrere Säuglinge zum gleichen Zeitpunkt ins MRT zu bestellen, um gegebenenfalls den am frühesten eingeschlafenen Säugling vorzuziehen. Diese Methode sollte zukünftig Einzug in die wissenschaftliche Untersuchung von bildgebenden Verfahren bei Säuglingen finden.
Die Umfrage im deutschsprachigen Raum zeigt eine Leitlinien-gerechte Betreuung von Säuglingen für bildgebende Verfahren, die mit einer hohen Qualität, und zumeist erfolgreich von erfahrenen Anästhesisten/Anästhesistinnen durchgeführt wird. Eventuelle Verbesserungen können im Bereich der Ausbildung nachfolgender Ärztinnen/Ärzte und in der häufigeren Verwendung der „feed-and-sleep“ Methode liegen, die vielen Kollegen/Kolleginnen bekannt ist, aber nur selten durchgeführt wird.
Ziel ist eine qualitativ hochwertige, schnellstmöglich durchgeführte Bildgebung, die ohne oder mit der niedrigst möglichen Dosierung eines sedierenden Medikamentes zu erreichen ist.