Refine
Has Fulltext
- yes (12)
Is part of the Bibliography
- yes (12)
Year of publication
Document Type
- Doctoral Thesis (12)
Keywords
- Kernspintomografie (12) (remove)
Institute
- Graduate School of Life Sciences (12) (remove)
Sonstige beteiligte Institutionen
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.
Myocardial infarction (MI) is a leading cause of death worldwide. Timely restoration of coronary blood flow to ischemic myocardium significantly reduces acute infarct mortality and attenuates ventricular remodeling. However, surviving MI patients frequently develop heart failure, which is associated with reduced quality of life, high mortality rate (10% annually), as well as high healthcare expenditures. The main processes involved in the evolution of heart failure post-MI are the great loss of contractile cardiomyocytes during ischemia-reperfusion and the subsequent complex structural and functional alterations, which are rooted in modifications at molecular and cellular levels in both the infarcted and non-infarcted myocardium. However, we still lack efficient treatments to prevent the development and progression of left ventricular remodeling. The improved survival rate of acute MI patients combined with the lack of effective therapy for post-MI remodeling contributes to the high prevalence of heart failure. Cardiac Magnetic Resonance Imaging (MRI) is an important tool for diagnosis and assessment of MI. With the advancement of this technology, the frontier of MRI has been extended to probing molecular and cellular events in vivo and non-invasively. In combination with assessment of morphology and function, the visualization of essential molecular and cellular markers in vivo could provide comprehensive, multifaceted views of the healing process in infarcted hearts, which might give new insight for the treatment of acute MI. In this thesis, molecular and cellular cardiac MRI methods were established to visualize and investigate inflammation and calcium flux in the healing process of acute MI in vivo, in a clinically relevant rat model.
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.
This work considered the frequency-modulated balanced steady-state free precession (fm-bSSFP) sequence as a tool to provide banding free bSSFP MR images. The sequence was implemented and successfully applied to suppress bandings in various in vitro and in vivo examples. In combination with a radial trajectory it is a promising alternative for standard bSSFP applications. First, two specialized applications were shown to establish the benefits of the acquisition strategy in itself. In real time cardiac imaging, it was shown that the continuous shift in frequency causes a movement of the bandings across the FOV. Thus, no anatomical region is constantly impaired, and a suitable timeframe can be found to examine all important structures. Furthermore, a combination of images with different artifact positions, similar to phase-cycled acquisitions is possible. In this way, fast, banding-free imaging of the moving heart was realized. Second, acquisitions with long TR were shown. While standard bSSFP suffers from increasing incidence of bandings with higher TR values, the frequency-modulated approach provided banding free images, regardless of the TR.
A huge disadvantage of fm-bSSFP, in combination with the radial trajectory, is the decrease in signal intensity. In this work a specialized reconstruction method, the multifrequency reconstruction for frequency-modulated bSSFP (Muffm), was established, which successfully compensated that phenomena. The application of Muffm to several anatomical sites, such as inner ear, legs and cardiac acquisitions, proofed the advantageous SNR of the reconstruction.
Furthermore, fm-bSSFP was applied to the clinically highly relevant task of water-fat separation. Former approaches of a phase-sensitive separation procedure in combination with standard bSSFP showed promising results but failed in cases of high inhomogeneity or high field strengths where banding artifacts become a major issue. The novel approach of using the fm-bSSFP acquisition strategy with the separation approach provided robust, reliable images of high quality. Again, losses in signal intensity could be regained by Muffm, as both approaches are completely compatible.
Opposed to conventional banding suppression techniques, like frequency-scouts or phase-cycling, all reconstruction methods established in this work rely on a single radial acquisition, with scan times similar to standard bSSFP scans. No prolonged measurement times occur and patient time in the scanner is kept as short as possible, improving patient comfort, susceptibility to motion or physiological noise and cost of one scan.
All in all, the frequency-modulated acquisition in combination with specializes reconstruction methods, leads to a completely new quality of images with short acquisition times.
Epidemiologische Studien schätzen die Inzidenz chronischer Obstipation auf bis zu 27% der Gesamtbevölkerung. Betroffenen Patienten ist die Stuhlentleerung nicht oder nur unter großer Anstrengung und nicht selten nur unter Zuhilfenahme der Hand möglich. Häufig sind funktionelle Pathologien, welche sich nur während der Defäkation ausbilden, hierfür verantwortlich. Daher ist für die Diagnose und Evaluation dieser Pathologien ein bildgebendes Verfahren notwendig, welches die dynamische Darstellung der Defäkation ermöglicht. Der Goldstandard zur Untersuchung von Patienten mit funktionellen Beckenbodenstörungen ist die Entero-Colpo-Cysto-Defäkographie (ECCD). Diese Durchleuchtungsmethode erfordert die Applikation ionisierender Strahlung im Bereich des Beckens. Außerdem müssen für die Untersuchung Rektum und Vagina mit bariumhaltigem Kontrastmittel, der Dünndarm mit barium- und iodhaltigem Kontrastmittel und zusätzlich die Blase mit iodhaltigem Kontrastmittel gefüllt werden. Bei der MR-Defäkographie hingegen ist keine ionisierende Strahlung notwendig und nur eine rektale Füllung mit Ultraschallgel als Kontrastmittel erforderlich. Zudem ermöglichen statische Aufnahmen aufgrund des hohen Weichteilkontrasts der MR-Bildgebung eine detaillierte Darstellung des gesamten Beckenbodens. Die MR-Bildgebung ist jedoch im Vergleich zu anderen Bildgebungsmodalitäten, wie beispielsweise der radiographischen Durchleuchtung, langsam. Besonders zur Darstellung dynamischer Prozesse ist daher eine starke Beschleunigung des Akquisitionsprozesses notwendig. Bei der Standard 2D MR-Defäkographie wird für die Beschleunigung der Datenakquisition eine regelmäßige zweifache Unterabtastung des k-Raums vorgenommen. Hierdurch lassen sich aber nur drei räumlich voneinander getrennte
68
2D Schichten mit einer zeitlichen Aktualisierungsrate der drei Schichten von ca. 1s akquirieren. Dadurch ist aber besonders die Diagnose lateral lokalisierter Pathologien eingeschränkt oder gar nicht möglich. Daher wurde in dieser Arbeit eine 3D MR-Defäkographie zur dynamischen Darstellung der Defäkation innerhalb eines vollständigen 3D Volumens entwickelt, implementiert und anhand von 9 Patientenmessungen optimiert. Die letzten 4 Patienten wurden mit den optimierten Sequenzparametern untersucht. Ausgehend von der kartesischen Datenakquisition der bestehenden 2D MRDefäkographie wurden zunächst dreidimensionale kartesische Trajektorien zur Datenakquisition und dafür geeignete Algorithmen zur Datenrekonstruktion untersucht. In diesem Zusammenhang wurde ein GRAPPA Centric-Out Akquisitionsschema in Kombination mit einer GRAPPA Datenrekonstruktion vorgestellt. Es zeigte sich jedoch, dass eine Stack-of-Stars Trajektorie in Bezug auf die stabile, rauscharme, dynamische Darstellung der Defäkation, vorteilhaft gegenüber der untersuchten kartesischen GRAPPA Centric-Out Trajektorie ist. Zur weiteren Optimierung der Messsequenz wurden daher drei radiale Stackof-Stars Akquisitionsschemata untersucht: Das Standard Stack-of-Stars Schema sowie zwei mit View-Sharing und zwei unterschiedlichen Dichtegewichtungen modifizierte Stack-of-Stars Schemata (DW-Sampling 1 und DW-Sampling 2). Das View-Sharing ermöglicht durch die Umstellung der Reihenfolge der akquirierten Partitionen nahezu eine Verdopplung der rekonstruierten Zeitpunkte der dynamisch gemessenen Zeitserie. Die Dichtegewichtung bewirkt, dass in den zentralen Partitionen mehr radiale Speichen gemessen werden und damit das k-Raum Zentrum dichter abgetastet wird als in den äußeren Partitionen. Beim Dichtegewichtungsschema DW-Sampling 2 ist der Abfall der Anzahl der innerhalb einer Partition gemessenen Speichen stärker als beim DW-Sampling 1. Trotzdem führte das mit View-Sharing und DW-Sampling 2 modifizierte Stackof-Stars Akquisitionsschema in Verbindung mit der FISTA Compressed Sensing Datenrekonstruktion zum besten Kompromiss zwischen erreichbarer räumlicher
69
und zeitlicher Auflösung. Dieses optimierte Setup ermöglicht die dynamische Darstellung der Defäkation in 7 Schichten eines vollständigen 3D Volumens mit einer Volumenaktualisierungsrate von 1,3s. Im Vergleich zur standardmäßig durchgeführten 2D MR-Defäkographie ist daher eine mehr als doppelt so große Abdeckung mit einer vergleichbaren zeitlichen Aktualisierungsrate und einer etwas geringeren räumlichen Auflösung gewährleistet. Hierdurch lassen sich zusätzlich zu den gewöhnlichen zentral gelegenen Pathologien auch lateral ausgeprägte Pathologien besser abdecken und diagnostizieren.
This work deals with the acceleration of cardiovascular MRI for the assessment
of functional information in steady-state contrast and for viability assessment
during the inversion recovery of the magnetization. Two approaches
are introduced and discussed in detail. MOCO-MAP uses an exponential
model to recover dynamic image data, IR-CRISPI, with its low-rank plus
sparse reconstruction, is related to compressed sensing.
MOCO-MAP is a successor to model-based acceleration of parametermapping
(MAP) for the application in the myocardial region. To this end, it
was augmented with a motion correction (MOCO) step to allow exponential
fitting the signal of a still object in temporal direction. Iteratively, this
introduction of prior physical knowledge together with the enforcement of
consistency with the measured data can be used to reconstruct an image
series from distinctly shorter sampling time than the standard exam (< 3 s
opposed to about 10 s). Results show feasibility of the method as well as
detectability of delayed enhancement in the myocardium, but also significant
discrepancies when imaging cardiac function and artifacts caused already by
minor inaccuracy of the motion correction.
IR-CRISPI was developed from CRISPI, which is a real-time protocol
specifically designed for functional evaluation of image data in steady-state
contrast. With a reconstruction based on the separate calculation of low-rank
and sparse part, it employs a softer constraint than the strict exponential
model, which was possible due to sufficient temporal sampling density via
spiral acquisition. The low-rank plus sparse reconstruction is fit for the use on
dynamic and on inversion recovery data. Thus, motion correction is rendered
unnecessary with it.
IR-CRISPI was equipped with noise suppression via spatial wavelet filtering.
A study comprising 10 patients with cardiac disease show medical
applicability. A comparison with performed traditional reference exams offer
insight into diagnostic benefits. Especially regarding patients with difficulty
to hold their breath, the real-time manner of the IR-CRISPI acquisition provides
a valuable alternative and an increase in robustness.
In conclusion, especially with IR-CRISPI in free breathing, a major acceleration
of the cardiovascular MR exam could be realized. In an acquisition
of less than 100 s, it not only includes the information of two traditional
protocols (cine and LGE), which take up more than 9.6 min, but also allows
adjustment of TI in retrospect and yields lower artifact level with similar
image quality.
Fulminant myocarditis is rare but a potentially life-threatening disease. Acute or mild myocarditis following acute ischemia is generally associated with a profound activation of the host’s immune system. On one hand this is mandatory to protect the host’s heart by fighting the invading agents (i.e., bacteria, viruses or other microbial agents) and/or to induce healing and repair processes in the damaged myocardium. On other hand, uncontrolled activation of the immune system may result in the generation of auto-reactive (not always beneficial) immune cells.
Myocarditis or inflammatory cardiomyopathy is characterized by focal or diffuse infiltrates, myocyte necrosis and/or apoptosis and subsequent fibrotic replacement of the heart muscle. In humans, about 30% of the myocarditis-patients develop dilated cardiomyopathy. As the clinical picture of myocarditis is multifaceted, it is difficult to diagnose the disease. Therefore, the main goal of the present work was to test and further develop novel non-invasive methods for the detection of myocardial inflammation by employing both contrast enhanced MRI techniques as well as novel nuclear tracers that are suitable for in vivo PET/ SPECT imaging.
As a part of this thesis, a pre-clinical animal model was successfully established by immunizing female Lewis rats with whole-porcine cardiac myosin (CM). Induction of Experimental Autoimmune Myocarditis (EAM) is considered successful when anti-myosin antibody titers are increased more than 100-fold over control animals and pericardial effusion develops. In addition, cardiac tissues from EAM-rats versus controls were analyzed for the expression of various pro-inflammatory and fibrosis markers. To further exploit non-invasive MRI techniques for the detection of myocarditis, our EAM-rats were injected either with (1) ultra-small Paramagnetic iron oxide particles (USPIO’s; Feraheme®), allowing for in vivo imaging , (2) micron sized paramagnetic iron oxide particles (MPIO) for ex vivo inflammatory cell-tracking by cMRI, or (3) with different radioactive nuclear tracers (67gallium citrate, 68gallium-labeled somatostatin analogue, and 68gallium-labeled cyclic RGD-peptide) which in the present work have been employed for autoradiographic imaging, but in principle are also suitable for in vivo nuclear imaging (PET/SPECT). In order to compare imaging results with histology, consecutive heart sections were stained with hematoxylin & eosin (HE, for cell infiltrates) and Masson Goldner trichrome (MGT, for fibrosis); in addition, immuno-stainings were performed with anti-CD68 (macrophages), anti-SSRT2A (somatostatin receptor type 2A), anti-CD61 (β3-integrins) and anti-CD31 (platelet endothelial cell adhesion molecule 1).
Sera from immunized rats strongly reacted with cardiac myosin. In immunized rats, echocardiography and subsequent MRI revealed huge amounts of pericardial effusion (days 18-21). Analysis of the kinetics of myocardial infiltrates revealed maximal macrophage invasion between days 14 and 28. Disappearance of macrophages resulted in replacement-fibrosis in formerly cell-infiltrated myocardial areas. This finding was confirmed by the time-dependent differential expression of corresponding cytokines in the myocardium. Immunized animals reacted either with an early or a late pattern of post-inflammation fibrosis. Areas with massive cellular infiltrates were easily detectible in autoradiograms showing a high focal uptake of 67gallium-citrate and 68gallium labeled somatostatin analogues (68Ga DOTA-TATE). Myocardium with a loss of cardiomyocytes presented a high uptake of 68gallium labeled cyclic RGD-peptide (68Ga NOTA-RGD). MRI cell tracking experiments with Feraheme® as the contrast-agent were inconclusive; however, strikingly better results were obtained when MPIOs were used as a contrast-agent: histological findings correlated well with in vivo and ex vivo MPIO-enhanced MRI images.
Imaging of myocardial inflammatory processes including the kinetics of macrophage invasion after microbial or ischemic damage is still a major challenge in, both animal models and in human patients. By applying a broad panel of biochemical, histological, molecular and imaging methods, we show here that different patterns of reactivity may occur upon induction of myocarditis using one and the same rat strain. In particular, immunized Lewis rats may react either with an early or a late pattern of macrophage invasion and subsequent post-inflammation fibrosis. Imaging results achieved in the acute inflammatory phase of the myocarditis with MPIOs, 67gallium citrate and 68gallium linked to somatostatin will stimulate further development of contrast agents and radioactive-nuclear tracers for the non-invasive detection of acute myocarditis and in the near future perhaps even in human patients.
Morphological and Functional Ultrashort Echo Time (UTE) Magnetic Resonance Imaging of the Human Lung
(2019)
In this thesis, a 3D Ultrashort echo time (3D-UTE) sequence was introduced in the Self-gated Non-Contrast-Enhanced Functional Lung Imaging (SENCEFUL) framework. The sequence was developed and implemented on a 3 Tesla MR scanner. The 3D-UTE technique consisted of a nonselective RF pulse followed by a koosh ball quasi-random sampling order of the k-space. Measurements in free-breathing and without contrast agent were performed in healthy subjects and a patient with lung cancer.
A gating technique, using a combination of different coils with high signal correlation, was evaluated in-vivo and compared with a manual approach of coil selection. The gating signal offered an estimation of the breathing motion during measurement and was used as a reference to segment the acquired data into different breathing phases.
Gradient delays and trajectory errors were corrected during post-processing using the Gradient Impulse Response Function. Iterative SENSE was then applied to determine the fully sampled data.
In order to eliminate signal changes caused by motion, a 3D image registration was employed, and the results were compared to a 2D image registration method.
Ventilation was assessed in 3D and regionally quantified by monitoring the signal changes in the lung parenchyma. Finally, image quality and quantitative ventilation values were compared to the standard 2D-SENCEFUL technique.
3D-UTE, combined with an automatic gating technique and SENCEFUL MRI, offered ventilation maps with high spatial resolution and SNR. Compared to the 2D method, UTE-SENCEFUL greatly improved the clinical quality of the structural images and the visualization of the lung parenchyma.
Through‐plane motion, partial volume effects and ventilation artifacts were also reduced with a three-dimensional method for image registration.
UTE-SENCEFUL was also able to quantify regional ventilation and presented similar results to previous studies.
Acceleration is a central aim of clinical and technical research in magnetic resonance imaging (MRI) today, with the potential to increase robustness, accessibility and patient comfort, reduce cost, and enable entirely new kinds of examinations. A key component in this endeavor is image reconstruction, as most modern approaches build on advanced signal and image processing. Here, deep learning (DL)-based methods have recently shown considerable potential, with numerous publications demonstrating benefits for MRI reconstruction. However, these methods often come at the cost of an increased risk for subtle yet critical errors. Therefore, the aim of this thesis is to advance DL-based MRI reconstruction, while ensuring high quality and fidelity with measured data. A network architecture specifically suited for this purpose is the variational network (VN). To investigate the benefits these can bring to non-Cartesian cardiac imaging, the first part presents an application of VNs, which were specifically adapted to the reconstruction of accelerated spiral acquisitions. The proposed method is compared to a segmented exam, a U-Net and a compressed sensing (CS) model using qualitative and quantitative measures. While the U-Net performed poorly, the VN as well as the CS reconstruction showed good output quality. In functional cardiac imaging, the proposed real-time method with VN reconstruction substantially accelerates examinations over the gold-standard, from over 10 to just 1 minute. Clinical parameters agreed on average.
Generally in MRI reconstruction, the assessment of image quality is complex, in particular for modern non-linear methods. Therefore, advanced techniques for precise evaluation of quality were subsequently demonstrated.
With two distinct methods, resolution and amplification or suppression of noise are quantified locally in each pixel of a reconstruction. Using these, local maps of resolution and noise in parallel imaging (GRAPPA), CS, U-Net and VN reconstructions were determined for MR images of the brain. In the tested images, GRAPPA delivers uniform and ideal resolution, but amplifies noise noticeably. The other methods adapt their behavior to image structure, where different levels of local blurring were observed at edges compared to homogeneous areas, and noise was suppressed except at edges. Overall, VNs were found to combine a number of advantageous properties, including a good trade-off between resolution and noise, fast reconstruction times, and high overall image quality and fidelity of the produced output. Therefore, this network architecture seems highly promising for MRI reconstruction.
Cardiovascular disease is one of the leading causes of death worldwide and, so far, echocardiography, nuclear cardiology, and catheterization are the gold standard techniques used for its detection. Cardiac magnetic resonance (CMR) can replace the invasive imaging modalities and provide a "one-stop shop" characterization of the cardiovascular system by measuring myocardial tissue structure, function and perfusion of the heart, as well as anatomy of and flow in the coronary arteries. In contrast to standard clinical magnetic resonance imaging (MRI) scanners, which are often operated at a field strength of 1.5 or 3 Tesla (T), a higher resolution and subsequent cardiac parameter quantification could potentially be achieved at ultra-high field, i.e., 7 T and above.
Unique insights into the pathophysiology of the heart are expected from ultra-high field MRI, which offers enhanced image quality in combination with novel contrast mechanisms, but suffers from spatio-temporal B0 magnetic field variations. Due to the resulting spatial misregistration and intra-voxel dephasing, these B0-field inhomogeneities generate a variety of undesired image artifacts, e.g., artificial image deformation. The resulting macroscopic field gradients lead to signal loss, because the effective transverse relaxation time T2* is shortened. This affects the accuracy of T2* measurements, which are essential for myocardial tissue characterization. When steady state free precession-based pulse sequences are employed for image acquisition, certain off-resonance frequencies cause signal voids. These banding artifacts complicate the proper marking of the myocardium and, subsequently, systematic errors in cardiac function measurements are inevitable. Clinical MR scanners are equipped with basic shim systems to correct for occurring B0-field inhomogeneities and resulting image artifacts, however, these are not sufficient for the advanced measurement techniques employed for ultra-high field MRI of the heart.
Therefore, this work focused on the development of advanced B0 shimming strategies for CMR imaging applications to correct the spatio-temporal B0 field variations present in the human heart at 7 T. A novel cardiac phase-specific shimming (CPSS) technique was set up, which featured a triggered B0 map acquisition, anatomy-matched selection of the shim-region-of-interest (SROI), and calibration-based B0 field modeling. The influence of technical limitations on the overall spherical harmonics (SH) shim was analyzed. Moreover, benefits as well as pitfalls of dynamic shimming were debated in this study. An advanced B0 shimming strategy was set up and applied in vivo, which was the first implementation of a heart-specific shimming approach in human UHF MRI at the time.
The spatial B0-field patterns which were measured in the heart throughout this study contained localized spots of strong inhomogeneities. They fluctuated over the cardiac cycle in both size and strength, and were ideally addressed using anatomy-matched SROIs. Creating a correcting magnetic field with one shim coil, however, generated eddy currents in the surrounding conducting structures and a resulting additional, unintended magnetic field. Taking these shim-to-shim interactions into account via calibration, it was demonstrated for the first time that the non-standard 3rd-order SH terms enhanced B0-field homogeneity in the human heart. However, they were attended by challenges for the shim system hardware employed in the presented work, which was indicated by the currents required to generate the optimal 3rd-order SH terms exceeding the dynamic range of the corresponding shim coils. To facilitate dynamic shimming updated over the cardiac cycle for cine imaging, the benefit of adjusting the oscillating CPSS currents was found to be vital. The first in vivo application of the novel advanced B0 shimming strategy mostly matched the simulations.
The presented technical developments are a basic requirement to quantitative and functional CMR imaging of the human heart at 7 T. They pave the way for numerous clinical studies about cardiac diseases, and continuative research on dedicated cardiac B0 shimming, e.g., adapted passive shimming and multi-coil technologies.
Für die Verwendung von zellbasierten Therapeutika ist vor allem die korrekt Identifikation
sowohl vom Ausgangsmaterial wie auch dem produziertem Material von
zentraler Wichtigkeit. In dieser Arbeit wurde eine Methodik entwickelt, welche eine
nicht-invasive Klassifizierung von Zellen und zellulärer Entwicklung aufgrund ihrer
zweidimensionalen Magnetresonanz-Korrelationsspektren ermöglichte.
Hierzu wurde ein mobiler MR-Scanner mit einer Feldstärke von 0.5T und einem Isozentrum
von 1 cm3 verwendet. Aufgrund der kompakten und leichten Bauweise war
es möglich, das System in normalen Zellkulturlaboren zu verwenden. Von den Proben
wurde ein zweidimensionales T1/T2 -Korrelationsspektrum aufgenommen, anhand
dessen die Zellen klassifiziert werden sollten. Mithilfe von Agarose-Dotagraf® -Zell-
Phantomen konnte die Stabilität und Reproduzierbarkeit des Messsystems und der
verwendeten Sequenz validiert werden.
Aufgrund der unter Umständen recht langen Messzeiten der MR-Technologie war
auch die Handhabung und Kultur der Zellproben während des Messprozesses von
großer Bedeutung. Um hierfür den Durchsatz an Proben zu erhöhen, wurde eine kostengünstige
und ebenfalls mobile Robotikanlage entwickelt. Diese basierte auf dem
kommerziell erhältlichen Roboterarm Braccio, welcher durch einen Arduino Mega
Mikrocontroller gesteuert wurde. Mit bis zu 24 Proben pro Tag konnte durch die
Automatisierung der Durchsatz an Proben um den Faktor 3 – 4 gesteigert werden.
Durch den entwickelten Prozess war es möglich, eine umfangreiche Datenbank –
bestehend aus 362 unabhängigen Messungen (biologische Replikate) – aufzubauen.
Die Datenbank enthielt Messungen von zehn unterschiedlichen Zelllinien. Zusätzlich
wurden T1/T2 -Korrelationsspektren von mesenchymalen Stromazellen (MSCs)
vor und nach deren Differenzierung zu Adipocyten aufgenommen, um ihre zelluläre
Entwicklung nicht-invasiv charakterisieren zu können.
Die aufgenommenen Daten wurden mithilfe einer geeigneten Support Vector Machine
wie auch angepassten künstlichen neuronalen Netzwerken klassifiziert. Mithilfe
dieser Methoden konnten die Zelllinien und MSCs anhand ihrer aufgenommenen
Korrelationsspektren mit einer Genauigkeit von bis zu 98% klassifiziert werden.
Diese hohe Treffsicherheit legte den Schluss nahe, dass die Kombination aus nichtinvasiver,
zweidimensionaler T1/T2 -MR-Relaxometrie und der Verwendung von geeigneten
Methoden des machine learning und der künstlichen Intelligenz eine effiziente
Methodik für die nicht-invasive Klassifizierung von Zellen sowie zellulärer
Entwicklung darstellt.
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.