TY - THES A1 - Kleineisel, Jonas T1 - Variational networks in magnetic resonance imaging - Application to spiral cardiac MRI and investigations on image quality T1 - Variational Networks in der Magnetresonanztomographie - Anwendung auf spirale Herzbildgebung und Untersuchungen zur Bildqualität N2 - 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. N2 - Eine Beschleunigung des Bildgebungsprozesses ist heute ein wichtiges Ziel von klinischer und technischer Forschung in der Magnetresonanztomographie (MRT). Dadurch könnten Robustheit, Verfügbarkeit und Patientenkomfort erhöht, Kosten gesenkt und ganz neue Arten von Untersuchungen möglich gemacht werden. Da sich die meisten modernen Ansätze hierfür auf eine fortgeschrittene Signal- und Bildverarbeitung stützen, ist die Bildrekonstruktion ein zentraler Baustein. In diesem Bereich haben Deep Learning (DL)-basierte Methoden in der jüngeren Vergangenheit bemerkenswertes Potenzial gezeigt und eine Vielzahl an Publikationen konnte deren Nutzen in der MRT-Rekonstruktion feststellen. Allerdings besteht dabei das Risiko von subtilen und doch kritischen Fehlern. Daher ist das Ziel dieser Arbeit, die DL-basierte MRT-Rekonstruktion weiterzuentwickeln, während gleichzeitig hohe Bildqualität und Treue der erzeugten Bilder mit den gemessenen Daten gewährleistet wird. Eine Netzwerkarchitektur, die dafür besonders geeignet ist, ist das Variational Network (VN). Um den Nutzen dieser Netzwerke für nicht-kartesische Herzbildgebung zu untersuchen, beschreibt der erste Teil dieser Arbeit eine Anwendung von VNs, welche spezifisch für die Rekonstruktion von beschleunigten Akquisitionen mit spiralen Auslesetrajektorien angepasst wurden. Die vorgeschlagene Methode wird mit einer segmentierten Rekonstruktion, einem U-Net, und einem Compressed Sensing (CS)-Modell anhand von qualitativen und quantitativen Metriken verglichen. Während das U-Net schlecht abschneidet, zeigen die VN- und CS-Methoden eine gute Bildqualität. In der funktionalen Herzbildgebung beschleunigt die vorgeschlagene Echtzeit-Methode mit VN-Rekonstruktion die Aufnahme gegenüber dem Goldstandard wesentlich, von etwa zehn zu nur einer Minute. Klinische Parameter stimmen im Mittel überein. Die Bewertung von Bildqualität in der MRT-Rekonstruktion ist im Allgemeinen komplex, vor allem für moderne, nichtlineare Methoden. Daher wurden anschließend forgeschrittene Techniken zur präsizen Analyse von Bildqualität demonstriert. Mit zwei separaten Methoden wurde einerseits die Auflösung und andererseits die Verstärkung oder Unterdrückung von Rauschen in jedem Pixel eines untersuchten Bildes lokal quantifiziert. Damit wurden lokale Karten von Auflösung und Rauschen in Rekonstruktionen durch Parallele Bildgebung (GRAPPA), CS, U-Net und VN für MR-Aufnahmen des Gehirns berechnet. In den untersuchten Bildern zeigte GRAPPA gleichmäßig eine ideale Auflösung, aber merkliche Rauschverstärkung. Die anderen Methoden verhalten sich lokal unterschiedlich je nach Struktur des untersuchten Bildes. Die gemessene lokale Unschärfe unterschied sich an den Kanten gegenüber homogenen Bildbereichen, und Rauschen wurde überall außer an Kanten unterdrückt. Insgesamt wurde für VNs eine Kombination von verschiedenen günstigen Eigenschaften festgestellt, unter anderem ein guter Kompromiss zwischen Auflösung und Rauschen, schnelle Laufzeit, und hohe Qualität und Datentreue der erzeugten Bilder. Daher erscheint diese Netzwerkarchitektur als ein äußerst vielversprechender Ansatz für MRT-Rekonstruktion. KW - Kernspintomografie KW - Convolutional Neural Network KW - Maschinelles Lernen KW - Bildgebendes Verfahren KW - magnetic resonance imaging KW - convolutional neural network KW - variational network KW - cardiac imaging KW - machine learning KW - local point-spread function KW - resolution KW - g-factor Y1 - 2024 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-347370 ER - TY - THES A1 - Portmann, Johannes T1 - Accelerated inversion recovery MRI of the myocardium using spiral acquisition T1 - Beschleunigte Inversion-Recovery MR-Bildgebung des Myokards mit spiralen Auslesezügen N2 - 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. N2 - Diese Arbeit behandelt die Beschleunigung der kardiovaskulären MRT zum Erfassen funktioneller Information bei Steady-State-Kontrast und zur Unter- suchung der Vitalität bei Wiederherstellung der Magnetisierung nach ihrer Inversion. Zwei Ansätze werden eingeführt und im Detail diskutiert: MOCO- MAP, welches ein exponentielles Modell nutzt, um dynamische Daten zu rekonstruieren, und IR-CRISPI, welches mit seinem “low-rank plus sparse"- Algorithmus mit Compressed Sensing verwandt ist. MOCO-MAP ist der Nachfolger der modellbasierten Beschleunigung des Parameter-Mappings (MAP) für die Anwendung im Bereich des Myokards. Hierzu wurde es mit einer Bewegungskorrektur (MOCO) versehen, um expo- nentielles Fitten eines unbewegten Objects in Zeitrichtung zu ermöglichen. Das Einbringen dieses physikalischen Vorwissens zusammen mit dem Erzwin- gen von Konsistenz mit den Messdaten wird dazu genutzt, iterativ eine Bildfolge aus Daten einer deutlich kürzeren Messung als herkömmlich zu rekonstruieren (< 3 s gegenüber ca. 10 s). Die Ergebnisse zeigen die Umsetz- barkeit der Methode sowie die Nachweisbarkeit von Delayed Enhancements im Myokard, aber deutliche funktionelle Abweichungen und Artefakte bereits aufgrund von kleinen Ungenauigkeiten der Bewegungskorrektur. IR-CRISPI geht aus CRISPI hervor, welches zur Auswertung von funk- tionellen Echtzeitdaten bei konstantem Kontrast dient. Mit der Rekon- struktion durch getrennte Berechnung von niedrigrangigem und dünnbe- setztem Matrixanteil wird hier bei der Datenrekonstruktion weniger stark eingeschränkt als bei einem strikten exponentiellen Modell. Die pirale Auf- nahmeweise erlaubt hierzu ausreichend effiziente k-Raumabdeckung. Die “low-rank plus sparse"-Rekonstruktion ist kompatibel mit dynamischen und mit Inversion-Recovery-Daten. Eine Bewegungskorrektur ist folglich nicht nötig. IR-CRISPI wurde mit einer Rauschunterdrückung durch räumliche Wavelet- Filterung versehen. Eine Studie, die 10 Patienten einschließt, zeigt die Eignung für die medizinische Anwendung. Der Vergleich mit herkömm- lichen Aufnahmetechniken lässt auf den gewonnenen diagnostischen Nutzen schließen. Besonders für Patienten, die Schwierigkeiten mit dem Luftanhal- ten haben, eröffnet diese Echtzeitaufnahmemethode eine wertvolle Alterna- tive und erhöhte Stabilität. Am Ende konnte gerade mittels IR-CRISPI eine bemerkenswerte Beschleu- nigung der kardiovaskulären MR-Untersuchung verwirklicht werden. Trotz der kurzen Aufnahmezeit von weniger als 100 s für den kompletten linken Ven- trikel schließt es nicht nur die Information zweier herkömmlicher Protokolle mit ein (Cine und LGE), die zusammen mehr als 9,6 min dauern, sondern es erlaubt zusätzlich auch das Einstellen der TI-Zeit im Nachhinein und liefert Ergebnisse mit geringerem Artefaktlevel bei ähnlicher Bildqualität KW - Kernspintomografie KW - Herzfunktion KW - Herzmuskel KW - Bildgebendes Verfahren KW - Echtzeit KW - cine loop KW - late enhancement KW - late gadolinium-enhancement KW - magnetic resonance imaging KW - real-time imaging KW - spiral trajectory Y1 - 2023 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-302822 ER - TY - THES A1 - Wech, Tobias T1 - Compressed Sensing in der funktionellen kardialen Magnetresonanztomographie T1 - Compressed sensing in functional cardiac magnetic resonance imaging N2 - 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. N2 - Because of its high reproducibility and its low variability, magnetic resonance imaging (MRI) of the heart is considered the gold-standard for assessing the cardiac function. In preclinical research, magnetic resonance imaging equally provides an accurate characterization of the cardiac function and enables an excellent analysis of modeled diseases. However, there is still a need for improvement in both applications. Clinical cardiac MRI represents a sophisticated procedure featuring long scan times. This renders the examination comparatively expensive. In preclinical imaging of small animals, long scan times have to be accepted to obtain the required high spatial and temporal resolution. Fast imaging is thus essential for an effective application of cardiac MRI in large collectives. Besides the improvement of the scanner hardware and the optimization of imaging sequences, research in the last decade concentrated on procedures to accelerate MR-data acquisition by exploiting information theory. While numerous publications were associated with parallel imaging (PI) at the beginning of this millennium, the compressed sensing theory (CS) recently gained more and more interest. The latter technique enables the reconstruction of signals from undersampled linear measurements (e.g. the Fourier basis) by exploiting the sparsity of the signal in any known transform domain. As MRI is perfectly qualified for an application of CS, a lot of publications already report on dedicated research. However, the algorithms needed for the reconstruction of undersampled data significantly alter the imaging process of MRI. Classical MRI could be assumed to be linear and stationary in a sufficiently good approximation. The introduction of CS into MRI means a change towards a non-linear and non-stationary transformation. Object information is no longer transferred into an image independently from its location and proportional to its intensity. The image is rather the result of an optimization process maximizing both the fidelity to measured data as well as the sparsity of the signal. The first chapter of this thesis describes a method to objectively evaluate the image quality of MR images reconstructed by CS algorithms. The acceleration with CS promises a reduction of scan time while preserving the image quality. The latter, however, has only been assessed qualitatively or in an insufficient quantitative manner. While classical (linear and stationary) MRI could be validated robustly and effectively by determining a point spread function (PSF), CS algorithms prohibit a corresponding analysis in an analogous manner due to their non-linear and non-stationary behavior. Therefore, the application of a local point spread function (LPSF) was investigated for the iterative soft thresholding (IST) algorithm used in this thesis, to enable a comparative evaluation for imaging systems including CS. The LPSF considers the local dependency of the CS algorithm and thus has to be determined in every location (pixel) of an image. In addition, the LPSF was defined in the linear part of the CS transformation. Small local perturbation on the image to be evaluated were reconstructed for this purpose. The width of the main lobe of the LPSF was used to perform spatially resolved studies on the resolution. The influence of typical undersampling schemes for CS as well as the usage of a discrete gradient transform for a further sparsification were investigated. Subsequently, the procedure was used to assess the spatial and temporal resolution in cardiac MRI. For all CS reconstructions performed in this work, the method allowed a solid and objective analysis of the image resolution. While up to now, comparisons to a fully sampled reference are widely used for a quality assessment, the proposed resolution evaluation represents a step towards a standardized analysis of images obtained by exploiting CS acceleration. The study on sampling schemes revealed, that also for CS accelerated acquisitions, the highest frequencies of the desired k-space have to be included. Former publications proposed undersampling patterns which partly featured a strong weight towards the center of k-space. The results of this thesis put these findings into perspective, as a loss in resolution has been observed for according approaches, at least for the simulations performed in this work. The dynamic acquisitions which were reconstructed exploiting x-f-sparsity proved to be prone to temporal blurring. This is substantiated by the suppression of high temporal frequencies and was analyzed by means of spatially resolved maps. Besides the resolution, an investigation of potential aliasing artifacts as well as the signal-to-noise-ratio (SNR) is essential for a comprehensive quality evaluation. While aliasing may also be investigated by means of the entries of the LPSF outside the main lobe, a modification of the multi-replica method proposed by Robson et al. was presented in chapter 5 to analyze the noise in CS reconstructed images. Taking into account all quality parameters, a robust evaluation of image quality is possible, even when CS is included in the imaging process. This allows a more objective comparison between new developments and present standard procedures and thus may aid the introduction of CS in clinical imaging. After the theoretical analysis on image quality, the next part of this thesis reports on the first application of CS to accelerate functional cardiac MRI of small animals. The studies were performed in cooperation with the British Heart Foundation Experimental Magnetic Resonance Unit (BMRU) of the University of Oxford. The algorithms needed for the CS acceleration were developed and optimized by means of the data acquired by the BMRU at their 9,4 T scanner.\\ An acceleration solely based on CS was investigated first. For this purpose, an ECG- and respiratory gated Cartesian cine acquisition was undersampled in phase encoding direction and reconstructed using CS. The dynamic time series was sparsified by determining temporal difference images for every time frame. A fully sampled temporal average image was obtained by varying the sampling pattern in the temporal dimension. Subsequently, this average image was subtracted from the images of individual heart phases, yielding the sparse temporal difference images. In the validation stage of the study, fully sampled cine acquisitions of mouse hearts were retrospectively undersampled in order to figure out the maximum possible CS acceleration. Cardiac functional parameters were determined for each group of a certain undersampling factor and compared by a statistical analysis. It was shown that a three-fold acceleration is possible without any degradation in the accuracy of the method. This undersampling factor was then validated in an accelerated measurement with a subsequent CS reconstruction. For this purpose, the sampling patterns were further optimized using the transform point spread function. In the subsequent chapter, the CS theory was combined with PI to further increase the acceleration. Again, the phase encoding direction of a Cartesian trajectory was undersampled. The acquisitions were performed using a 9,4 T scanner equipped with an 8 channel mouse coil. In order to exploit the potential of both techniques, CS and PI were combined in a serial manner. First, the k-space was equidistantly undersampled to enable the application of PI. An additional undersampling according to pseudo random numbers was then performed on the resulting sub-grid to allow an acceleration by CS. In consequence, the reconstruction was performed in a serial manner, too. CS was first applied to reconstruct the equidistantly undersampled sub-grid. GRAPPA was used subsequently to compute the still missing data. The equidistantly undersampled sub-grid was shifted from heart phase to heart phase in order to obtain a fully sampled low temporal resolution k-space for a calibration of the GRAPPA-weights. This procedure spares the acquisition time of a separate calibration scan. The following combinations were investigated by retrospectively undersampling a fully sampled cine dataset: R_CS x R_PI = 2 x 2, 2 x 3, 3 x 2 and 3 x 3. The analysis of the noise behavior, the systematic error and the resolution leads to the conclusion that a six-fold acceleration is possible using the proposed hybrid technique. While an increasing factor of the CS acceleration resulted in a slightly larger systematic error, a higher PI acceleration factor led to a slight noise enhancement. However, noise was suppressed for increasing CS acceleration at the same time. In summary, the deviations were at a level which allowed accelerations of up to R_CS x R_PI = 3 x 3. The determination of LPSFs showed a loss in spatial resolution of approximately 50% for a six-fold and up to 64% for nine-fold acceleration. Obviously, the observed suppression of noise was paid by a reduced image resolution for the comparatively noisy acquisitions in small animals. The increased blurring at the endocardial border impedes the segmentation and represents a possible source of error. Taking into consideration all results, a six-fold acceleration (R_CS x R_PI = 2 x 3, 3 x 2) seems reasonable. The additional usage of PI thus enables a further acceleration by a factor of 2 in comparison to an exclusive application of CS. In summary, CS enables a distinct reduction of scan time in preclinical functional cardiac MRI of small animals. Even if no phased-array-coils are available, the necessary amount of data can be reduced to one third without impairing the accuracy of left-ventricular volumes and mass measurements. For acquisitions performed with phased-array-coils, the three-fold acceleration by PI can be extended by an additional two-fold CS acceleration to a joint factor of R=6. Therefore, CS may contribute to an effective application of cardiac MRI in small animals for large collectives. In the last part of the thesis, a modality for clinical functional MRI of the heart was developed. CS was used to enable the acquisition of the whole heart in a single breath-hold of the patient. The current method of choice usually acquires 10-15 2D-slices of the heart, while each measurement requires a separate breath-hold. An undersampled 3D-trajectory was used to reach the necessary acceleration. Phase-encoding in one direction and radial projections in the two remaining ones allowed for an effective acquisition below the Nyquist-criterion. The sparsification of the image series was achieved by subtracting a temporal average image as performed for the preclinical studies. The functionality of the reconstruction technique at an acceleration factor of R ~ 10 was validated in a simulation based on a retrospectively undersampled dataset. The differences between the CS reconstructed and the fully sampled dataset were negligible and thus, the proposed trajectory was implemented at the scanner. An image series depicting the cardiac function with coverage of the full heart was acquired in a single breath-hold of a healthy volunteer using this sequence. Data not covered by the trajectory were reconstructed by the algorithm developed in the validation stage. Due to the slice profile of the fast excitation pulses, additional aliasing artifacts were present in the outer partitions with respect to the images obtained in the simulation. Streaking artifacts of a low intensity were still visible. Apart from that, the dynamics of the heart were excellently captured. The high contrast between the blood pool and the myocardium perfectly qualifies the images for the assessment of cardiac functional parameters. Therefore, the method allows for a higher patient comfort and throughput compared to the gold standard by drastically reducing the amount of necessary breath-holds to a single one. KW - Kernspintomografie KW - Magnetresonanztomographie KW - Kardiale MR-Bildgebung KW - Compressed Sensing KW - magnetic resonance imaging KW - cardiac magnetic resonance imaging KW - compressed sensing KW - Signalregenerierung KW - NMR-Tomographie Y1 - 2012 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-77179 ER - TY - THES A1 - Ye, Yuxiang T1 - Molecular and Cellular Magnetic Resonance Imaging of Myocardial Infarct Healing T1 - Molekulare und zelluläre Magnetresonanztomographie der Wundheilung nach Myokardinfarkt N2 - Myokardinfarkte (MI) sind eine der häufigsten Todesursachen weltweit. Eine rechtzeitige Wiederherstellung des koronaren Blutflusses im ischämischen Myokard reduziert signifikant die Sterblichkeitsrate akuter Infarkte und vermindert das ventrikuläre Remodeling. Überlebende MI-Patienten entwickeln jedoch häufig eine Herzinsuffizienz, die mit einer reduzierten Lebensqualität, hohen Sterblichkeitsrate (10% jährlich), sowie hohen Kosten für das Gesundheitssystem einhergeht. Die Entwicklung der Herzinsuffizienz nach einem MI ist auf den hohen Verlust kontraktiler Kardiomyozyten, während der Ischämie-Reperfusion zurückzuführen. Anschließende komplexe strukturelle und funktionelle Veränderungen resultieren aus Modifikationen des infarzierten und nicht infarzierten Myokards auf molekularer und zellulärer Ebene. Die verbesserte Überlebensrate von Patienten mit akutem MI und das Fehlen effizienter Therapien, die die Entwicklung und das Fortschreiten des ventrikulären Remodelings verhindern, führen zu einer hohen Prävalenz der Herzinsuffizienz. Die kardiale Magnetresonanztomographie (MRT) ist eine wichtige Methode zur Diagnose und Beurteilung des Myokardinfarktes. Mit dem technologischen Fortschritt wurden die Grenzen der MRT erweitert, so dass es heute möglich ist, auch molekulare und zelluläre Ereignisse in vivo und nicht-invasiv zu untersuchen. In Kombination mit kardialer Morphologie und Funktion könnte die Visualisierung essentieller molekularer und zellulärer Marker in vivo weitreichende Einblicke in den Heilungsprozess infarzierter Herzen liefern, was zu neuen Erkenntnissen für ein besseres Verständnis und bessere Therapien des akuten MI führen könnte. In dieser Arbeit wurden Methoden für die molekulare und zelluläre kardiale MRT-Bildgebung der Inflammation und des Kalziumstroms im Heilungsprozess des akuten Myokardinfarktes in vivo in einem Rattenmodel mit klinischer Relevanz etabliert. N2 - 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. KW - Kernspintomografie KW - Herzinfarkt KW - magnetic resonance imaging KW - myocardial infarction KW - Wundheilung KW - NMR-Tomographie Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:bvb:20-opus-72514 ER -