87.61.-c Magnetic resonance imaging
Refine
Has Fulltext
- yes (11)
Is part of the Bibliography
- yes (11)
Document Type
- Doctoral Thesis (11)
Keywords
- Kernspintomografie (6)
- MRI (5)
- NMR-Tomographie (5)
- NMR-Bildgebung (3)
- Herzbildgebung (2)
- Kernspintomographie (2)
- MRT (2)
- Relaxometrie (2)
- Relaxometry (2)
- Steady-State Sequences (2)
Institute
Sonstige beteiligte Institutionen
In this work, accelerated non-Cartesian Magnetic Resonance Imaging (MRI) methods were established and applied to cardiovascular imaging (CMR) at different magnetic field strengths (3T and 7T).
To enable rapid data acquisition, highly efficient spiral k-space trajectories were created. In addition, hybrid sampling patterns such as the twisting radial lines (TWIRL) k-space trajectory were studied.
Imperfections of the dynamic gradient system of a MR scanner result in k-space sampling errors. Ultimately, these errors can lead to image artifacts in non-Cartesian acquisitions.
Among other reasons such as an increased reconstruction complexity, they cause the lack of spiral sequences in clinical routine compared to standard Cartesian imaging.
Therefore, the Gradient System Transfer Functions (GSTFs) of both scanners were determined and used for k-space trajectory correction in post-correction as well as in terms of a pre-emphasis.
The GSTF pre-emphasis was implemented as a fully automatic procedure, which enabled a precise correction of arbitrary gradient waveforms for double-oblique slice orientations.
Consequently, artifacts due to trajectory errors could be mitigated, which resulted in high image quality in non-Cartesian MRI.
Additionally, the GSTF correction was validated by measuring pre-emphasized spiral gradient outputs, which showed high agreement with the theoretical gradient waveforms.
Furthermore, it could be demonstrated that the performance of the GSTF correction is superior to a simple delay compensation approach.
The developed pulse sequences were applied to gated as well as real-time CMR. Special focus lied on the implementation of a spiral imaging protocol to resolve the beating heart of animals and humans in real time and free breathing.
In order to achieve real-time CMR with high spatiotemporal resolution, k-space undersampling was performed. For this reason, efficient sampling strategies were developed with the aim to facilitate compressed sensing (CS) during image reconstruction.
The applied CS approach successfully removed aliasing artifacts and yielded high-resolution cardiac image series. Image reconstruction was performed offline in all cases such that the images were not available immediately after acquisition at the scanner.
Spiral real-time CMR could be performed in free breathing, which led to an acquisition time of less than 1 minute for a whole short-axis stack.
At 3T, the results were compared to the gold standard of electrocardiogram-gated Cartesian CMR in breath hold, which revealed similar values for important cardiovascular functional and volumetric parameters.
This paves the way to an application of the developed framework in clinical routine of CMR.
In addition, the spiral real-time protocol was transferred to swallowing and speech imaging at 3T, and first images were presented.
The results were of high quality and confirm the straightforward utilization of the spiral sequence in other fields of MRI.
In general, the GSTF correction yielded high-quality images at both field strengths, 3T and 7T.
Off-resonance related blurring was mitigated by applying non-Cartesian readout gradients of short duration. At 7T, however, B1-inhomogeneity led to image artifacts in some cases.
All in all, this work demonstrated great advances in accelerating the MRI process by combining efficient, undersampled non-Cartesian k-space coverage with CS reconstruction.
Trajectory correction using the GSTF can be implemented at any scanner model and enables non-Cartesian imaging with high image quality.
Especially MRI of dynamic processes greatly benefits from the presented rapid imaging approaches.
Ziel dieser Arbeit ist es, die quantitative MRT in den Fokus zu rücken. In den letzten Jahren hat sich auf diesem Forschungsgebiet viel weiterentwickelt und es wurden verschiedenste Sequenzen und Methoden vorgestellt, um insbesondere Relaxationszeitparameter quantitativ in kurzer Zeit zu messen. Steady-State-Sequenzen eignen sich besonders für diese Thematik, da sie kurze Messzeiten benötigen und darüber hinaus ein relativ hohes SNR besitzen. Speziell die IR TrueFISP-Sequenz bietet für die Parameterquantifizierung viel Potential. Ursprünglich wurde diese Sequenz an der Universität Würzburg zur simultanen Messung von T1- und T2-Relaxationszeiten vorgestellt und hinsichtlich der Zeiteffizienz weiterentwickelt. In dieser Arbeit wurde ein neuartiger iterativer Rekonstruktionsansatz für die IR TrueFISP-Sequenz entwickelt, der auf einer Hauptkomponentenanalyse (PCA) basiert und sich die glatten Signalverläufe zu Nutze macht. Aufgrund der hohen Zeitauflösung dieser Rekonstruktionstechnik werden dabei auch Gewebekomponenten mit kurzen Relaxationszeiten detektierbar. Weiterhin bewahrt der Rekonstruktionsansatz Informationen mehrerer Gewebekomponenten innerhalb eines Voxels und ermöglicht damit eine relaxographische Untersuchung. Insbesondere beim Menschen führen der Partialvolumeneffekt und die Mikrostruktur des Gewebes zu Signalverläufen, die ein multi-exponentielles Signal liefern. Die MR-Relaxographie, also die Darstellung von Relaxationszeitverteilungen innerhalb eines Voxels, stellt eine Möglichkeit dar, um die beteiligten Gewebekomponenten aus dem überlagerten Signalverlauf zu extrahieren. Insgesamt bilden die optimierte Relaxometrie mit der Möglichkeit der analytischen Korrektur von Magnetfeldinhomogenitäten und die beschleunigte Relaxographie die Hauptteile dieser Dissertation. Die Hauptkapitel werden im Folgenden noch einmal gesondert zusammengefasst.
Die simultane Aufnahme der quantitativen T1- und T2-Parameter-Karten kann mit einem
Goldenen-Winkel-basiertem radialen IR TrueFISP-Readout in ungefähr 7 Sekunden pro
Schicht erreicht werden. Die bisherige Rekonstruktionstechnik mit dem KWIC-Filter ist
durch dessen breite Filter-Bandbreite und somit in der zeitlichen Auflösung limitiert. Besonders bei hohen räumlichen Frequenzen wird eine sehr große Anzahl an Projektionen
zusammengefasst um ein Bild zu generieren. Dies sorgt dafür, dass Gewebekomponenten mit kurzer T1*-Relaxationszeit (z.B. Fett oder Myelin) nicht akkurat aufgelöst werden können. Um dieses Problem zu umgehen, wurde die T1* shuffling-Rekonstruktion entwickelt, die auf dem T2 Shuffling-Ansatz basiert. Diese Rekonstruktionstechnik macht sich die glatten Signalverläufe der IR TrueFISP-Sequenz zu Nutze und ermöglicht die Anwendung einer PCA. Die iterative Rekonstruktion sorgt dafür, dass mit nur acht kombinierten Projektionen pro generiertem Bild eine merklich verbesserte temporäre Auflösung erzielt werden kann. Ein Nachteil ist jedoch das stärkere Rauschen in den ersten Bildern der Zeitserie bedingt durch die angewandte PCA. Dieses verstärkte Rauschen äußert sich in den leicht erhöhten Standardabweichungen in den berechneten Parameter-Karten. Jedoch ist der Mittelwert näher an den Referenzwerten im Vergleich zu den Ergebnissen mit dem KWIC-Filter. Letztendlich kann man sagen, dass die Ergebnisse leicht verrauschter, aber exakter sind.
Mittels zusätzlichen Regularisierungstechniken oder Vorwissen bezüglich des Rauschlevels
wäre es zudem noch möglich, das SNR der ersten Bilder zu verbessern, um dadurch den
beschriebenen Effekt zu verringern.
Grundsätzlich hängt die Genauigkeit von IR TrueFISP vom T1/T2-Verhältnis des betreffenden
Gewebes und dem gewählten Flipwinkel ab. In dieser Arbeit wurde der Flipwinkel besonders für weiße und graue Masse im menschlichen Gehirn optimiert. Mit den verwendeten 35° wurde er außerdem etwas kleiner gewählt, um zudem Magnetisierungstransfereffekte zu minimieren. Mit diesen Einstellungen ist die Präzision vor allem für hohe T1- und niedrige T2-Werte sehr gut, wird jedoch insbesondere für höhere T2-Werte schlechter. Dies ist aber ein generelles Problem der IR TrueFISP-Sequenz und hängt nicht mit der entwickelten Rekonstruktionsmethode zusammen. Außerdem wurde im fünften Kapitel eine
Akquisitionstechnik vorgestellt, die eine 3D-Abdeckung der quantitativen Messungen des Gehirns in klinisch akzeptabler Zeit von unter 10 Minuten erzielt. Dies wird durch Einsatz der parallelen Bildgebung erreicht, da eine Kombination aus radialer Abtastung in der Schicht und kartesischer Aufnahme in Schichtrichtung (Stack-of-Stars) vorliegt.
Ein großes Problem in der Steady-State-Sequenz (und somit auch bei IR TrueFISP) sind
Magnetfeldinhomogenitäten, die durch Suszeptibilitätsunterschiede verschiedener Gewebe und/oder Inhomogenitäten des Hauptmagnetfeldes hervorgerufen werden. Diese führen zu Signalauslöschungen und damit verbunden zu den beschriebenen Banding-Artefakten. Mithilfe der analytisch ermittelten Korrekturformeln ist es nun möglich, die berechneten (T1,T2)-Wertepaare unter Berücksichtigung der tatsächlich auftretenden Off- Resonanzfrequenz für einen großen Bereich zu korrigieren. An den kritischen Stellen, an denen die Bandings auftreten, liefert jedoch auch diese Korrektur keine brauchbaren Ergebnisse. Grundsätzlich ist es für die Genauigkeit der Ergebnisse stets zu empfehlen, die Flipwinkel- und B0-Karte zusätzlich mit aufzunehmen, um diese Parameter für die quantitative Auswertung exakt zu kennen. Mit den beschriebenen Methoden aus Kapitel 6 könnte es prinzipiell auch möglich sein, die Off-Resonanzfrequenz aus dem Signalverlauf zu ermitteln und auf die zusätzliche Messung der B0-Karte zu verzichten. B0-Änderungen während der Messung, die von der Erwärmung der passiven Shim-Elemente im MR-System hervorgerufen werden, sind kaum zu korrigieren. Ein stabiler Scanner ohne B0-Drift ist deshalb für quantitative Auswertungen erforderlich.
Die erwähnte Messzeit von 7 Sekunden pro Schicht garantiert, dass auch Gewebe mit längeren Relaxationskomponenten annähernd im Steady-State sind, was wiederum für das Umkehren des Signals in den abklingenden Verlauf gegen Null und die anschließende
Multikomponentenanalyse (vgl. Kapitel 7) notwendig ist. Mit der inversen Laplace-
Transformation ist es innerhalb eines Voxels möglich, Signalverläufe auf mehrere Komponenten hin zu untersuchen. Der ursprünglich angenommene mono-exponentielle Verlauf wird durch ein multi-exponentielles Verhalten abgelöst, was vor allem in biologischem Gewebe eher der Wahrheit entspricht. Gewebe mit kurzen Relaxationskomponenten (T1* < 200 ms) sind klinisch relevant und mit T1* shuffling detektierbar. Vor allem Myelin innerhalb des Gehirns ist bei neurologischen Fragestellungen ein Indikator zur Diagnose im Frühstadium (z.B. für neurodegenerative Erkrankungen) und
deshalb von besonderem Interesse. Die Integration über verschiedene T1*-Zeitbereiche im
T1*-Spektrum ermöglicht dazu die Erstellung von Gewebekomponenten-Karten, mithilfe
derer klinische Auswertungen sinnvoll wären. Die Erstellung dieser Karten ist prinzipiell
möglich und funktioniert für mittlere und lange Gewebekomponenten recht gut. Die
klinisch relevanten kurzen Gewebekomponenten sind dagegen bei der radialen Aufnahme
mit nur einem Schuss noch nicht befriedigend. Deshalb wurde die Aufnahmetechnik in
eine quasi-zufällige kartesische Akquisition mit mehreren Schüssen weiterentwickelt. Die Ergebnisse wurden in Kapitel 7 vorgestellt und sind vielversprechend. Einzig die Messzeit sollte mit zusätzlichen Beschleunigungen noch weiter verkürzt und auf eine kartesische 3D-Akquisition erweitert werden.
Die Beschränkung auf T1*-Spektren bei der Multikomponentenanalyse und die Tatsache, dass deren Amplitude von einer Kombination von S0 und Sstst abhängen, führen dazu, dass es nicht ohne Weiteres möglich ist für einen einzelnen Gewebetyp an die T1- und T2-Information zu gelangen. In Kapitel 8 wurde gezeigt, dass dies mit einer zusätzlichen
Messung gelingen kann. Das finale Ergebnis dieser Messungen ohne und mit Inversion sind zweidimensionale Spektren, bei der für jede Gewebekomponente innerhalb eines Voxels der T1- und T2-Wert abgelesen werden kann. Wichtig hierbei ist die Tatsache, dass der verwendete Ansatz kein Vorwissen über die Anzahl der zu erwartenden Gewebekomponenten (Peaks) im Voxel voraussetzt. Auch bei dieser Methodik ist die Kenntnis über den tatsächlichen Flipwinkel von Bedeutung, da dieser in den Formeln zur Berechnung von T1 und T2 verwendet wird. Die Stabilität des B0-Feldes ist hier ebenso von enormer Bedeutung, da Änderungen zwischen den beiden Messungen zu einem unterschiedlichen Steady-State und somit zu Abweichungen bei den nachfolgenden Berechnungen führen, die auf den selben Steady-State-Wert ausgelegt sind.
Zusammenfassend lässt sich sagen, dass mit dieser Arbeit die Grundlagen für genauere
und robustere quantitative Messungen mittels Steady-State-Sequenzen gelegt wurden. Es
wurde gezeigt, dass sich Relaxationszeitspektren für jedes einzelne Voxel generieren lassen.
Dadurch ist eine verbesserte Auswertung möglich, um genauere Aussagen über die Zusammensetzung einer Probe (vor allem beim menschlichen Gewebe) treffen zu können. Zudem wurde die Theorie für ultraschnelle 2D-Relaxographie-Messungen vorgestellt. Erste”Proof of Principle“-Experimente zeigen, dass es möglich ist, 2D-Relaxationszeitspektren in sehr kurzer Zeit zu messen und graphisch darzustellen. Diese Aufnahme- und Datenverarbeitungstechnik ist in dieser Form einmalig und in der Literatur kann bis dato keine schnellere Methode gefunden 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.
Neue Therapieansätze durch Tissue Engineering erfordern gleichzeitig angepasste Diagnosemöglichkeiten und nicht-invasive Erfolgskontrollen. Speziell die 3D-MR-Bildgebung ist ein vielversprechendes Instrument, um Parameter mit hoher räumlicher Präzision zu quantifizieren. Vor diesem Hintergrund wurden im Rahmen dieser Arbeit neue Ansätze für die hochauflösende 3D-MRT in vivo entwickelt und deren Eignung im Bereich des Tissue Engineerings gezeigt.
Welchen Vorteil die Quantifizierung von Parametern bietet, konnte im Rahmen einer prä-klinischen Studie an einem Modell der Hüftkopfnekrose gezeigt werden. Der Therapieverlauf wurde zu verschiedenen Zeitpunkten kontrolliert. Trotz der niedrigen räumlichen Auflösung, konnten durch eine systematische Auswertung der Signalintensitäten von T1- und T2-FS-gewichteten Aufnahmen Rückschlüsse über Veränderungen in der Mikrostruktur gezogen werden, die darüber hinaus in guter Übereinstimmung mit Ergebnissen von ex vivo µCT-Aufnahmen waren. Dort konnte eine Verdickung der Trabekelstruktur nachgewiesen werden, welche sehr gut mit einer Signalabnahme in den T1-gewichteten Aufnahmen korrelierte. Die radiale Auswertung der Daten erlaubte dabei eine komprimierte Darstellung der Ergebnisse. Dadurch wurde eine effiziente Auswertung der umfangreichen Daten (verschiedene Tiere an mehreren Zeitpunkten mit einer Vielzahl an Einzelaufnahmen) ermöglicht und eine unabhängige Bewertung erreicht.
Um die Limitationen der begrenzten Auflösung von 2D-Multi-Schichtaufnahmen aufzuheben, wurden neue Ansätze für eine hochaufgelöste 3D-Aufnahme entwickelt. Hierfür wurden Spin-Echo-basierte Sequenzen gewählt, da diese eine genauere Abbildung der Knochenmikrostruktur erlauben als Gradienten-Echo-basierte Methoden. Zum einen wurde eine eigene 3D-FLASE-Sequenz entwickelt und zum anderen eine modifizierte 3D-TSE-Sequenz. Damit an Patienten Aufnahmen bei klinischer Feldstärke von 1,5 T mit einer hohen räumlichen Auflösung innerhalb einer vertretbaren Zeit erzielt werden können, muss eine schnelle und signalstarke Sequenz verwendet werden. Eine theoretische Betrachtung bescheinigte der TSE-Sequenz eine um 25 % höhere Signaleffizienz verglichen mit einer FLASE-Sequenz mit identischer Messzeit. Dieser Unterschied konnte auch im Experiment nachgewiesen werden. Ein in vivo Vergleich der beiden Sequenzen am Schienbein zeigte eine vergleichbare Darstellung der Spongiosa mit einer Auflösung von 160 × 160 × 400 µm.
Für die Bildgebung des Hüftkopfs mit der neuen Sequenz waren jedoch aufgrund der unterschiedlichen Anatomie weitere Modifikationen notwendig. Um längere Messzeiten durch ein unnötig großes Field-of-View zu vermeiden, mussten Einfaltungsartefakte unterdrückt werden. Dies wurde durch die orthogonale Anwendung der Anregungs- und Refokussierungspulse in der TSE-Sequenz effizient gelöst. Technisch bedingt konnte jedoch nicht eine vergleichbare Auflösung wie am Schienbein realisiert werden.
Der Vorteil der 3D-Bildgebung, dass Schichtdicken von deutlich weniger als 1 mm erreicht werden können, konnte jedoch erfolgreich auf den Unterkiefer übertragen werden. Der dort verlaufende Nervus Mandibularis ist dabei eine wichtige Struktur, deren Verlauf im Vorfeld von verschiedenen operativen Eingriffen bekannt sein muss. Er ist durch eine dünne knöcherne Wand vom umgebenden Gewebe getrennt. Im Vergleich mit einer 3D-VIBE-Sequenz zeigte die entwickelte 3D-TSE-Sequenz mit integrierter Unterdrückung von Einfaltungsartefakten eine ähnlich gute Lokalisierung des Nervenkanals über die gesamte Länge der Struktur. Dies konnte in einer Studie an gesunden Probanden mit verschiedenen Beobachtern nachgewiesen werden. Durch die neue Aufnahmetechnik konnte darüber hinaus die Auflösung im Vergleich zu bisherigen Studien deutlich erhöht werden, was insgesamt eine präzisere Lokalisierung des Nervenkanals erlaubt.
Ein Baustein des Tissue Engineerings sind bio-resorbierbare Materialien, deren Abbau- und Einwachsverhalten noch untersucht werden muss, bevor diese für die klinische Anwendung zugelassen werden. Die durchgeführten in vitro µMR-Untersuchungen an Polymerscaffolds zeigten die reproduzierbare Quantifizierung der Porengröße und Wandstärke. Darüber hinaus wurde eine inhomogene Verteilung der Strukturparameter beobachtet. Die Ergebnisse waren in guter Übereinstimmung mit µCT-Aufnahmen als Goldstandard. Unterschiedliche Varianten der Scaffolds konnten identifiziert werden. Dabei bewies sich die MR-Bildgebung als zuverlässige Alternative.
Insgesamt zeigen die Ergebnisse dieser Arbeit, welche Vorteile und Anwendungsmöglichkeiten die 3D-MRT-Bildgebung bietet, und dass auch mit klinischer Feldstärke in vivo Voxelgrößen im Submillimeterbereich für alle Raumrichtungen erreichbar sind. Die erzielten Verbesserungen in der räumlichen Auflösung erhöhen die Genauigkeit der verschiedenen Anwendungen und ermöglichen eine bessere Identifikation von kleinen Abweichungen, was eine frühere und zuverlässigere Diagnose für Patienten verspricht.
Die bSSFP-Sequenz kombiniert kurze Akquisitionszeiten mit einem hohen Signal-zu-Rausch-Verhältnis, was sie zu einer vielversprechenden Bildgebungsmethode macht. Im klinischen Alltag ist diese Technik jedoch bisher - abgesehen von vereinzelten Anwendungen - kaum etabliert. Die Hauptgründe hierfür sind Signalauslöschungen in Form von Bandingartefakten sowie der erzielte T2/T1-gewichtete Mischkontrast. Das Ziel dieser Dissertation war die Entwicklung von Methoden zur Lösung der beiden genannten Limitationen, um so eine umfassendere Verwendung von bSSFP für die MR-Diagnostik zu ermöglichen.
Magnetfeldinhomogenitäten, die im Wesentlichen durch Suszeptibilitätsunterschiede oder Imperfektionen seitens der Scannerhardware hervorgerufen werden, äußern sich bei der bSSFP-Bildgebung in Form von Bandingartefakten. Mit DYPR-SSFP (DYnamically Phase-cycled Radial bSSFP) wurde ein Verfahren vorgestellt, um diese Signalauslöschungen effizient zu entfernen. Während für bereits existierende Methoden mehrere separate bSSFP-Bilder akquiriert und anschließend kombiniert werden müssen, ist für die Bandingentfernung mittels DYPR-SSFP lediglich die Aufnahme eines einzelnen Bildes notwendig. Dies wird durch die neuartige Kombination eines dynamischen Phasenzyklus mit einer radialen Trajektorie mit quasizufälligem Abtastschema ermöglicht. Die notwendigen Bestandteile können mit geringem Aufwand implementiert werden. Des Weiteren ist kein spezielles Rekonstruktionsschema notwendig, was die breite Anwendbarkeit des entwickelten Ansatzes ermöglicht. Konventionelle Methoden zur Entfernung von Bandingartefakten werden sowohl bezüglich ihrer Robustheit als auch bezüglich der notwendigen Messzeit übertroffen.
Um die Anwendbarkeit von DYPR-SSFP auch jenseits der gewöhnlichen Bildgebung zu demonstrieren, wurde die Methode mit der Fett-Wasser-Separation kombiniert. Basierend auf der Dixon-Technik konnten so hochaufgelöste Fett- sowie Wasserbilder erzeugt werden. Aufgrund der Bewegungsinsensitivät der zugrunde liegenden radialen Trajektorie konnten die Messungen unter freier Atmung durchgeführt werden, ohne dass nennenswerte Beeinträchtigungen der Bildqualität auftraten. Die erzielten Ergebnisse am Abdomen zeigten weder Fehlzuordnungen von Fett- und Wasserpixeln noch verbleibende Bandingartefakte.
Ein Nachteil der gewöhnlichen Dixon-basierten Fett-Wasser-Separation ist es, dass mehrere separate Bilder zu verschiedenen Echozeiten benötigt werden. Dies führt zu einer entsprechenden Verlängerung der zugehörigen Messzeit. Abhilfe schafft hier die Verwendung einer Multiecho-Sequenz. Wie gezeigt werden konnte, ermöglicht eine derartige Kombination die robuste, bandingfreie Fett-Wasser-Separation in klinisch akzeptablen Messzeiten.
DYPR-SSFP erlaubt die Entfernung von Bandingartefakten selbst bei starken Magnetfeldinhomogenitäten. Dennoch ist es möglich, dass Signalauslöschungen aufgrund des Effekts der Intravoxeldephasierung verbleiben. Dieses Problem tritt primär bei der Bildgebung von Implantaten oder am Ultrahochfeld auf. Als Abhilfe hierfür wurde die Kombination von DYPR-SSFP mit der sogenannten z-Shim-Technik untersucht, was die Entfernung dieser Artefakte auf Kosten einer erhöhten Messzeit ermöglichte.
Die mit DYPR-SSFP akquirierten radialen Projektionen weisen aufgrund des angewendeten dynamischen Phasenzyklus leicht verschiedene Signallevel und Phasen auf. Diese Tatsache zeigt sich durch inkohärente Bildartefakte, die sich jedoch durch eine Erhöhung der Projektionsanzahl effektiv reduzieren lassen. Folglich bietet es sich in diesem Kontext an, Anwendungen zu wählen, bei denen bereits intrinsisch eine verhältnismäßig hohe Anzahl von Projektionen benötigt wird. Hierbei hat sich gezeigt, dass neben der hochaufgelösten Bildgebung die Wahl einer 3D radialen Trajektorie eine aussichtsreiche Kombination darstellt. Die in der vorliegenden Arbeit vorgestellte 3D DYPR-SSFP-Technik erlaubte so die isotrope bandingfreie bSSFP-Bildgebung, wobei die Messzeit im Vergleich zu einer gewöhnlichen bSSFP-Akquisition konstant gehalten werden konnte. Verbleibende, durch den dynamischen Phasenzyklus hervorgerufene Artefakte konnten effektiv mit einem Rauschunterdrückungsalgorithmus reduziert werden. Anhand Probandenmessungen wurde gezeigt, dass 3D DYPR-SSFP einen aussichtsreichen Kandidaten für die Bildgebung von Hirnnerven sowie des Bewegungsapparats darstellt.
Während die DYPR-SSFP-Methode sowie die darauf beruhenden Weiterentwicklungen effiziente Lösungen für das Problem der Bandingartefakte bei der bSSFP-Bildgebung darstellen, adressiert die vorgestellte RA-TOSSI-Technik (RAdial T-One sensitive and insensitive Steady-State Imaging) das Problem des bSSFP-Mischkontrasts. Die Möglichkeit der Generierung von T2-Kontrasten basierend auf der bSSFP-Sequenz konnte bereits in vorausgehenden Arbeiten gezeigt werden. Hierbei wurde die Tatsache ausgenutzt, dass der T1-Anteil des Signalverlaufs nach Beginn einer bSSFP-Akquisition durch das Einfügen von Inversionspulsen in ungleichmäßigen Abständen aufgehoben werden kann. Ein so akquiriertes Bild weist folglich einen reinen, klinisch relevanten T2-Kontrast auf. Die im Rahmen dieser Arbeit vorgestellte Methode basiert auf dem gleichen Prinzip, jedoch wurde anstelle einer gewöhnlichen kartesischen Trajektorie eine radiale Trajektorie in Kombination mit einer KWIC-Filter-Rekonstruktion verwendet. Somit können bei gleichbleibender oder sogar verbesserter Bildqualität aus einem einzelnen, mit RA-TOSSI akquirierten Datensatz verschiedene T2-Wichtungen als auch gewöhnliche T2/T1-Wichtungen generiert werden. Mittels Variation der Anzahl der eingefügten Inversionspulse konnte ferner gezeigt werden, dass es neben den besagten Wichtungen möglich ist, zusätzliche Kontraste zu generieren, bei denen verschiedene Substanzen im Bild ausgelöscht sind. Diese Substanzen können am Beispiel der Gehirnbildgebung Fett, graue Masse, weiße Masse oder CSF umfassen und zeichnen sich neben den reinen T2-Kontrasten durch eine ähnlich hohe klinische Relevanz aus. Die mögliche Bedeutung der vorgestellten Methode für die klinische Verwendung wurde durch Messungen an einer Gehirntumorpatientin demonstriert.
Zusammenfassend lässt sich sagen, dass die im Rahmen dieser Dissertation entwickelten Techniken einen wertvollen Beitrag zur Lösung der eingangs beschriebenen Probleme der bSSFP-Bildgebung darstellen. Mit DYPR-SSFP akquirierte Bilder sind bereits mit bestehender, kommerzieller Rekonstruktionssoftware direkt am Scanner rekonstruierbar. Die Software für die Rekonstruktion von RA-TOSSI-Datensätzen wurde für Siemens Scanner implementiert. Folglich sind beide Methoden für klinische Studien einsetzbar, was gleichzeitig den Ausblick dieser Arbeit darstellt.
This thesis focuses on various aspects and techniques of 19F magnetic resonance (MR). The first chapters provide an overview of the basic physical properties, 19F MR and MR sequences related to this work. Chapter 5 focuses on the application of 19F MR to visualize biological processes in vivo using two different animal models. The dissimilar models underlined the wide applicability of 19F MR in preclinical research. A subsection of Chapter 6 shows the application of compressed sensing (CS) to 19F turbo-spin-echo chemical shift imaging (TSE-CSI), which leads to reduced measurement time. CS, however, can only be successfully applied when a sufficient signal-to-noise ratio (SNR) is available. When the SNR is low, so-called spike artifacts occur with the CS algorithm used in the present work. However, it was shown in an additional subsection that these artifacts can be reduced using a CS-based post processing algorithm. Thus, CS might help overcome limitations with time consuming 19F CSI experiments. Chapter 7 deals with a novel technique to quantify the B+1 profile of an MR coil. It was shown that, using a specific application scheme of off resonant pulses, Bloch-Siegert (BS)-based B+1 mapping can be enabled using a Carr Purcell Meiboom Gill (CPMG)-based TSE sequence. A fast acquisition of the data necessary for B+1 mapping was thus enabled. In the future, the application of BS-CPMG-TSE B+1 mapping to improve quantification using 19F MR could therefore be possible.
Magnetic resonance imaging (MRI) is a medical imaging method that involves no ionizing radiation and can be used non-invasively. Another important - if not the most important - reason for the widespread and increasing use of MRI in clinical practice is its interesting and highly flexible image contrast, especially of biological tissue. The main disadvantages of MRI, compared to other widespread imaging modalities like computed tomography (CT), are long measurement times and the directly resulting high costs. In the first part of this work, a new technique for accelerated MRI parameter mapping using a radial IR TrueFISP sequence is presented. IR TrueFISP is a very fast method for the simultaneous quantification of proton density, the longitudinal relaxation time T1, and the transverse relaxation time T2. Chapter 2 presents speed improvements to the original IR TrueFISP method. Using a radial view-sharing technique, it was possible to obtain a full set of relaxometry data in under 6 s per slice. Furthermore, chapter 3 presents the investigation and correction of two major sources of error of the IR TrueFISP method, namely magnetization transfer and imperfect slice profiles. In the second part of this work, a new MRI thermometry method is presented that can be used in MRI-safety investigations of medical implants, e.g. cardiac pacemakers and implantable cardioverter-defibrillators (ICDs). One of the major safety risks associated with MRI examinations of pacemaker and ICD patients is RF induced heating of the pacing electrodes. The design of MRI-safe (or MRI-conditional) pacing electrodes requires elaborate testing. In a first step, many different electrode shapes, electrode positions and sequence parameters are tested in a gel phantom with its geometry and conductivity matched to a human body. The resulting temperature increase is typically observed using temperature probes that are placed at various positions in the gel phantom. An alternative to this local thermometry approach is to use MRI for the temperature measurement. Chapter 5 describes a new approach for MRI thermometry that allows MRI thermometry during RF heating caused by the MRI sequence itself. Specifically, a proton resonance frequency (PRF) shift MRI thermometry method was combined with an MR heating sequence. The method was validated in a gel phantom, with a copper wire serving as a simple model for a medical implant.
Diffusionstensorbildgebung im Vergleich zu anderen Parametermethoden für die Infarktcharakterisierung Ziel dieses Teils der Arbeit war die Klärung der Frage, welches Potential verschiedene MR-Parametersequenzen bei der Charakterisierung eines myokardialen Infarkts sowohl im akuten als auch im chronischen Fall haben. Dazu wurde eine Studie mit akut und chronisch infarzierten Rattenherzen durchgeführt. Untersucht wurden die Parameter T1, T2 und T2* sowie die aus der Diffusionstensorbildgebung berechneten Parameter ADC, FA, cs, cp und cl . Es zeigte sich, dass es kein Analogon zum bei einer cerebralen Ischämie bekannten Mismatch-Konzept gibt. Weder im akuten noch im chronischen war Fall eine ausgewiesene Differenz im diagnostizierten Infarktareal zwischen verschiedenen Sequenzen feststellbar. Alles in allem eignen sich zur detaillierten Charakterisierung der Infarktnarbe am besten eine T2*- oder eine Diffusionstensorsequenz. Die T2*-Sequenz liefert optisch das aufschlussreichere Bild, die aufwendigere Diffusionstensorsequenz dagegen bietet aufgrund der vielfachen Darstellungsmöglichkeiten im Postprocessing ein Mehr an Information und zeigt dazu eine Veränderung der Narbe im Zeitverlauf. Oxygenierungsmessung am Mäuseherz in vivo Die Charakterisierung einer Infarktnarbe kann auch über die Darstellung morphologischer Strukturen hinaus erfolgen. Die Oxygenierung ist ein komplexer Parameter, der funktionelle Auskunft über die Vaskularisierung und Viabilität des Gewebes geben kann. Zugang zu diesem Parameter erhält man über T2*-Messungen, da der Parameter T2* sensitiv auf chemisch gebundenen Sauerstoff reagiert. Hier wurden der Einfluss von reiner Sauerstoffatmung im Gegensatz zu normaler Raumluftatmung auf die Oxygenierung bei gesunden und infarzierten Mäusen untersucht. Die Messungen wurden trotz der Schwierigkeiten, die durch die Bewegung durch Atmung und Herzschlag entstehen, in vivo bei 17,6 Tesla implementiert und durchgeführt. Die Auflösung war ausreichend, um auch nach Infarkt extrem ausgedünnte Myokardwände gut auflösen und charakterisieren zu können. Der Effekt auf das Oxygenierungslevel ist stark unterschiedlich zwischen normalen und infarzierten Herzen, woraus auf eine noch nicht weit fortgeschrittene Revaskularisierung der Narbe eine Woche nach Infarzierung geschlossen werden kann. Die Methode wurde darüber hinaus an einem 7,0 Tesla-Magneten zur Verwendung an Ratten implementiert und auf das im Gegensatz zur Maus veränderte Atmungsverhalten der Ratte angepasst. Zum einen kann dadurch der Einfluss des hohen Magnetfeldes auf die Oxygenierungsmessung untersucht werden, zum anderen ist das Herz als zu untersuchendes Objekt bei der Ratte größer. Diffusionswichtung mittels Hole-Burning Die in dieser Arbeit zur Charakterisierung des Herzens verwendete Diffusionsmethode kann im Grenzfall von kurzen T2-Relaxationszeiten an ihre Grenzen stoßen: Bei den verwendeten starken Magnetfeldern klingt das messbare Signal aufgrund der Relaxationszeit T2 oft sehr schnell ab. Daher wurde eine Methode entwickelt, die einen völlig neuen Ansatz zur diffusionsgewichteten Bildgebung verfolgt, bei dem die Informationen über die Diffusion unabhängig von der limitierenden T2-Zeit gewonnen werden können. Die sog. Hole-Burning-Diffusionssequenz verwendet in einem Vorexperiment lediglich die Longitudinalmagnetisierung zur Diffusionswichtung. Das Signal wird dann mit einer schnellen Auslesesequenz akquiriert. Bei der Präparation werden zunächst auf Subvoxel-Niveau Streifen "gebrannt", d.h. die Magnetisierung wird dort gesättigt. Bis zur nächsten Sättigung ist das Verhalten der Magnetisierung abhängig von der T1-Relaxation in diesem Bereich und vom Diffusionsverhalten. Durch rasches Wiederholen des selektiven Pulszugs wird schließlich eine Gleichgewichtsmagnetisierung erreicht, die von der Diffusionskonstanten D und der T1-Relaxationszeit abhängt. Im Rahmen dieser Arbeit wurden die Abhängigkeiten verschiedener Sequenzparameter untersucht und diese mittels Simulationen optimiert. Außerdem wurde die Sequenz an einem Scanner implementiert und erste Experimente damit durchgeführt. Mit Hilfe von Simulationen konnten dazu Lookup-Tabellen generiert werden, mit denen in bestimmten Bereichen (insbesondere bei nicht zu kurzen T1-Relaxationszeiten) sowohl die Diffusionskonstante D als auch die T1-Relaxationszeit quantifiziert werden konnte.
Der Hauptteil der vorliegenden Arbeit befasste sich mit der Anwendung und der Entwicklung von neuen Methoden der spektroskopischen NMR-Bildgebung zur nicht-invasiven metabolischen Charakterisierung von Xenograft-Tumormodellen bei 17,6 T. In einem weiteren Abschnitt wurden verschiedene etablierte Methoden der lokalisierten NMR-Spektroskopie und der spektroskopischen Bildgebung genutzt, um den Metabolismus von Hülsenfrüchten (Pisum sativum) am Hochfeld zu untersuchen. Im experimentellen Teil der Arbeit wurde der selektive Mehrquantenfilter Sel-MQC zur Laktatbestimmung in neun verschiedenen Xenograft-Tumormodellen verwendet. Diese Werte wurden mit Ergebnissen aus der Biolumineszenz und mit der Tumorkontrolldosis 50 (TCD50) der Tumorlinien korreliert. Der Sel-MQC-Editierungsfilter stellte sich als äußerst robuste Methode heraus das Laktat im NMR-Spektrum eindeutig von koresonanten Lipiden des Unterhautfettgewebes bzw. von tumoreigenen Lipiden zu trennen. Der Vergleich mit dem durch die Biolumineszenz bestimmten Laktat zeigte durchweg niedrigere Werte in den NMR-Messungen. Der Hauptgrund für diesen Unterschied besteht wahrscheinlich darin, dass mit der NMR-Methode nur das freie Laktat bestimmt werden kann, wohingegen die Biolumineszenz das gesamte Laktat erfasst. Das mit der NMR detektierbare freie Laktat zeigte allerdings eine mäßige Korrelation zur TCD50 (R = 0,46), wodurch dieser Parameter nur als bedingt prognostisch wertvoll für die Strahlentherapie von Tumoren angesehen werden kann. Der Informationsgehalt pro Messzeit und damit die Effizienz der Standard-Sel-MQC-Editierungssequenz konnte durch verschiedene methodische Erweiterungen gesteigert werden. Eine zusätzliche spektral selektive Wasserunterdrückung und ein weiteres Aufnahmefenster ermöglichte neben der Messung des Laktatsignals die Akquisition sämtlicher Resonanzen des 1H-Spektrums mit einer kurzen Echozeit. Somit konnten zusätzlich das Gesamtcholin und die Methyl- und Metylengruppen der Lipide aufgenommen werden. Neben dem Laktat erwies sich das Verhältnis von Lipid-Methylensignal zu Gesamtcholin (L1/tCho) als aussagekräftigster Parameter, um zwei untersuchte Xenograft-Tumormodelle zu unter-scheiden. Die spektroskopische Sel-MQC-Bildgebungssequenz, deren k-Raumantastung in der Regel mit reiner Phasenkodierung durchgeführt wird, konnte durch eine Verwendung eines Lesegradienten beschleunigt werden. Die bei dem Sel-MQC-Filter auftretenden typischen Artefakte im Bereich der Wasserresonanz sind durch zwei Aufnahmen nach dem Dixon-Prinzip und einem anschließenden Additionsverfahren unterdrückbar. Bei einer ausreichenden Aufnahmezeit, die abhängig vom T2* der zu editierenden Resonanz ist, kann mit der Methode eine nahezu ähnlich hohe Sensitivität wie mit dem rein phasenkodierten Experiment erreicht werden. Eine in die Sequenz eingefügte frequenzselektive Refokussierung der Laktat-CH3-Gruppe ermöglichte die Aufnahme mehrerer Laktatechos ohne eine Phasenmodulation durch die J-Kopplung im Signal zu erhalten. Die nach einer Anregung erhaltenen Echos können zur weiteren Beschleunigung der Sequenz oder zur Bestimmung der apparenten transversalen Relaxationszeit des editieren Metaboliten verwendet werden. Das Grundprinzip des Sel-MQC-Filters konnte in einem umgekehrten Verfahren dazu verwendet werden mobile Lipide im Tumor ohne das koresonante Laktatsignal zu detektieren, um damit die Lipiddetektion zu spezifizieren. Da zur Unterdrückung des Metabolitensignals nur die J-Kopplung ausgenutzt wird, müssen weder Relaxationszeiten noch Diffusionskoeffizienten für die Editierung bekannt sein. Die Aufnahme des Lipidsignals wird dabei in einer Präparation erreicht, was die Sequenz robust gegenüber Bewegungsartefakten macht. Die Methode kann beispielsweise mit Diffusionsgradienten kombiniert werden, um den apparenten Diffusionskoeffizienten mobiler Lipide im Tumorgewebe zu bestimmen. Das hohe Magnetfeld von 17,6 T und damit die vergrößerte chemische Verschiebung eigneten sich insbesonders dazu spektroskopische Messungen an Pflanzensystemen durchzuführen. Im letzten Teil der Arbeit wurden unterschiedliche lokalisierte 1D-, 2D-NMR-Methoden und die spektroskopische Bildgebung verwendet, um den Wildtyp und eine Mutantenform des Pisum sativum nicht-invasiv metabolisch zu untersuchen. Die mit der NMR bestimmten Metabolitenkonzentrationen im Endosperm des Pisum sativum korrelierten mit Resultaten aus biochemischen Auswertungen. Weiterhin konnten mit den NMR-Methoden auch Ergebnisse gewonnen werden, die mit biochemischen und histologischen Verfahren nicht erreicht werden können. Die Untersuchung von Pflanzen – oder wie hier von Pflanzensamen – mit spektroskopischen NMR-Methoden bieten zusätzliche und für bestimmte Fragestellungen auch einzigartige Ansätze deren Metabolismus in vivo zu untersuchen.
The goal of the work presented in this thesis was to explore the possibilities and limitations of MRI / MRS using an ultra high field of 17.6 tesla. A broad range of specific applications and MR methods, from MRI to MRSI and MRS were investigated. The main foci were on sodium magnetic resonance spectroscopic imaging of rodents, magnetic resonance spectroscopy of the mouse brain, and the detection of small amounts of iron labeled stem cells in the rat brain using MRI Sodium spectroscopic imaging was explored since it benefits tremendously from the high magnetic field. Due to the intrinsically low signal in vivo, originating from the low concentrations and short transverse relaxation times, only limited results have been achieved by other researchers until now. Results in the literature include studies conducted on large animals such as dogs to animals as small as rats. No studies performed on mice have been reported, despite the fact that the mouse is the most important laboratory animal due to the ready availability of transgenic strains. Hence, this study concentrated on sodium MRSI of small rodents, mostly mice (brain, heart, and kidney), and in the case of the brain on young rats. The second part of this work concentrated on proton magnetic resonance spectroscopy of the rodent brain. Due to the high magnetic field strength not only the increasing signal but also the extended spectral resolution was advantageous for such kind of studies. The difficulties/limitations of ultra high field MRS were also investigated. In the last part of the presented work detection limits of iron labeled stem cells in vivo using magnetic resonance imaging were explored. The studies provided very useful benchmarks for future researchers in terms of the number of labeled stem cells that are required for high-field MRI studies. Overall this work has shown many of the benefits and the areas that need special attention of ultra high fields in MR. Three topics in MRI, MRS and MRSI were presented in detail. Although there are significant additional difficulties that have to be overcome compared to lower frequencies, none of the work presented here would have been possible at lower field strengths.