Refine
Has Fulltext
- yes (5)
Is part of the Bibliography
- yes (5)
Year of publication
- 2018 (5) (remove)
Document Type
- Doctoral Thesis (5)
Keywords
- Kernspintomografie (4)
- 3D-Bildgebung (2)
- Aorta (1)
- Bildgebendes Verfahren (1)
- DNMR-Spektroskopie (1)
- Dynamic magnetic resonance imaging (1)
- Dynamische MR Bildgebung (1)
- Dynamische Messung (1)
- Flussbildgebung (1)
- Herzbildgebung (1)
Institute
Sonstige beteiligte Institutionen
In Magnetic Resonance Imaging (MRI), acquisition of dynamic data may be highly complex due to rapid changes occurred in the object to be imaged. For clinical diagnostic, dynamic MR images require both high spatial and temporal resolution. The speed in the acquisition is a crucial factor to capture optimally dynamics of the objects to obtain accurate diagnosis. In the 90’s, partially parallel MRI (pMRI) has been introduced to shorten scan times reducing the amount of acquired data. These approaches use multi-receiver coil arrays to acquire independently and simultaneously the data.
Reduction in the amount of acquired data results in images with aliasing artifacts. Dedicated methods as such Sensitivity Encoding (SENSE) and Generalized Autocalibrating Partially Parallel Acquisition (GRAPPA) were the basis of a series of algorithms in pMRI.
Nevertheless, pMRI methods require extra spatial or temporal information in order to optimally reconstruct the data. This information is typically obtained by an extra scan or embedded in the accelerated acquisition applying a variable density acquisition scheme.
In this work, we were able to reduce or totally eliminate the acquisition of the training data for kt-SENSE and kt-PCA algorithms obtaining accurate reconstructions with high temporal fidelity.
For dynamic data acquired in an interleaved fashion, the temporal average of accelerated data can generate an artifact-free image used to estimate the coil sensitivity maps avoiding the need of extra acquisitions. However, this temporal average contains errors from aliased components, which may lead to signal nulls along the spectra of reconstructions when methods like kt-SENSE are applied. The use of a GRAPPA filter applied to the temporal average reduces these errors and subsequently may reduce the null components in the reconstructed data. In this thesis the effect of using temporal averages from radial data was investigated. Non-periodic artifacts performed by undersampling radial data allow a more accurate estimation of the true temporal average and thereby avoiding undesirable temporal filtering in the reconstructed images. kt-SENSE exploits not only spatial coil sensitivity variations but also makes use of spatio-temporal correlations in order to separate the aliased signals. Spatio-temporal correlations in kt-SENSE are learnt using a training data set, which consists of several central k-space lines acquired in a separate scan. The scan of these extra lines results in longer acquisition times even for low resolution images. It was demonstrate that limited spatial resolution of training data set may lead to temporal filtering effects (or temporal blurring) in the reconstructed data.
In this thesis, the auto-calibration for kt-SENSE was proposed and its feasibility was tested in order to completely eliminate the acquisition of training data. The application of a prior TSENSE reconstruction produces the training data set for the kt-SENSE algorithm. These training data have full spatial resolution. Furthermore, it was demonstrated that the proposed auto-calibrating method reduces significantly temporal filtering in the reconstructed images compared to conventional kt-SENSE reconstructions employing low resolution training images. However, the performance of auto-calibrating kt-SENSE is affected by the Signal-to-Noise Ratio (SNR) of the first pass reconstructions that propagates to the final reconstructions.
Another dedicated method used in dynamic MRI applications is kt-PCA, that was first proposed for the reconstruction of MR cardiac data. In this thesis, kt-PCA was employed for the generation of spatially resolved M0, T1 and T2 maps from a single accelerated IRTrueFISP or IR-Snapshot FLASH measurement. In contrast to cardiac dynamic data, MR relaxometry experiments exhibit signal at all temporal frequencies, which makes their reconstruction more challenging. However, since relaxometry measurements can be represented by only few parameters, the use of few principal components (PC) in the kt-PCA algorithm can significantly simplify the reconstruction. Furthermore, it was found that due to high redundancy in relaxometry data, PCA can efficiently extract the required information from just a single line of training data.
It has been demonstrated in this thesis that auto-calibrating kt-SENSE is able to obtain high temporal fidelity dynamic cardiac reconstructions from moderate accelerated data avoiding the extra acquisition of training data. Additionally, kt-PCA has been proved to be a suitable method for the reconstruction of highly accelerated MR relaxometry data.
Furthermore, a single central training line is necessary to obtain accurate reconstructions. Both reconstruction methods are promising for the optimization of training data acquisition and seem to be feasible for several clinical applications.
The noninvasive magnetic resonance imaging technique allows for the investigation of functional processes in the living plant. For this purpose during this work, different NMR imaging methods were further developed and applied.
For the localisation of the intrusion of water into the germinating rape seed with the simultaneous depiction of the lipid-rich tissue via a 3D rendering, in Chap. 5 the technique of interleaved chemical selective acquisition of water and lipid was used in the germinating seed. The utilization of high-resolution MR images of germinated seeds enabled the localization of a predetermined water gap in the lipid-rich aleurone layer, which resides directly under the seed coat. The for a long time in biology prevalent discussion, whether such a gap exists or the seed soaks up the water from all sides, rather like a sponge, could hereby, at least for the rapeseed seed, be answered clearly. Furthermore, the segmentation and 3D visualization of the vascular tissue in the rapeseed seeds was enabled by the high-resolution datasets, a multiply branched structure preconstructed in the seed could be shown. The water is directed by the vascular tissue and thus awakens the seed gradually to life. This re-awakening could as well be tracked by means of invasive imaging via an oxygen sensor. In the re-awakened seeds, the lipid degradation starts, other than expected, not in the lipid-rich cotyledons but in the residual endosperm remaining from seed development and in the aleurone layer which previously protected the embryo. Within this layer, the degradation could be verified in the high-resolution MR datasets.
The method presented in Chap. 6 provides a further characteristic trait for phenotyping of seeds and lipid containing plants in general. The visualization of the compounds of fatty acids in plant seeds and fruits could be achieved by the distinct utilization of chemical shift-selective imaging techniques. Via the application of a CSI sequence the fatty acid compounds in an olive were localized in a 2D slice. In conjunction with an individually adjusted CHESS presaturation module Haa85 the high-resolution 3D visualization of saturated and unsaturated fatty acid compounds in different seeds was achieved. The ratio maps calculated from these datasets allow to draw conclusions from the developmental stage or the type of seed. Furthermore, it could be shown that the storage condition of two soybean seeds with different storage time durations lead to no degradation of the fatty acid content.
Additional structural information from inside of dry seeds are now accessible via MRI. In this work the imaging of cereal seeds could be significantly improved by the application of the UTE sequence. The hitherto existing depictions of the lipid distribution, acquired with the spin echo sequence, were always sufficient for examinations of the lipid content, yet defects in the starchy endosperm or differences in the starch concentration within the seed remained constantly unseen with this technique. In a direct comparison of the datasets acquired with the previous imaging technique (spin echo) and with UTE imaging, the advantage of data acquisition with UTE could be shown. By investigating the potential seed compounds (starch, proteins, sugar) in pure form, the constituent parts contributing to the signal could be identified as bound water (residual moisture) and starch. The application of a bi-exponential fit on the datasets of the barley seed enabled the separate mapping of magnetization and of relaxation time of two components contributing to the NMR signal. The direct comparison with histological stainings verified the previous results, thus this technique can be used for the selective imaging of starch in dry seeds.
Conclusions on the translocation characteristics in plants can be drawn by the technique proposed in Chap. 8. The associated translocation velocities can now, even in the range of several um/h, be determined in the living plant. Based on calculated concentrations of an MR contrast agent, which was taken up by the plant, these translocation velocities were estimated both in longitudinal direction, thus along the vascular bundle, and in horizontal direction, thus out of the bundle. The latter velocity is located below the contrast agent's velocity value of free diffusion. By adjusting a dynamic contrast-enhancing imaging technique (DCE-Imaging, Tof91) the acquisition duration of a T1-map was significantly reduced. By means of these maps, local concentrations of the contrast agent in plant stems and the siliques of the rapeseed plant could be determined.
Numerous questions in plant science can only be answered by non-invasive techniques such as MRI. For this reason, besides the experimental results achieved in this work, further NMR methods were tested and provided for the investigation of plants.
As an example, the study on the imaging of magnetic exchange processes are mentioned, which provided the groundwork for a possible transfer of CEST experiments (Chemical Exchange Saturation Transfer) to the plant. The results are presented in the bachelor thesis of A. Jäger Jae17, which was performed under my supervision, they find great interest under biologists.
The development of new technologies, which extend the possibilities for the investigation of living organisms, is of great importance. For this reason, I have contributed to the development of the currently unpublished method RACETE (Refocused Acquisition of Chemical Exchange Transferred Excitations [Jak17, Reu17, Gut18a]). By rephasing the transferred magnetization the utilization of properties which have not been available in chemical "`exchange"' experiments is enabled. With this method a positive contrast is generated, thus a reference experiment is not mandatory. Furthermore, the image phase, which in classical experiments contains no information about the exchanged protons, can be used for the distinct identification of multiple substances which have been excited simultaneously.
This recently at the Department of Experimental Physics V developed method can be used in particular for the identification of lipids and for the localization of sugars and amino acids, thus it can serve the enhancement and improvement of non-invasive analytical methods.
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.
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.