500 Naturwissenschaften und Mathematik
Refine
Has Fulltext
- yes (5)
Is part of the Bibliography
- yes (5)
Document Type
- Doctoral Thesis (5)
Keywords
- Kernspintomografie (5) (remove)
Institute
The subject of this work was to develop, implement, optimize and apply methods for quantitative MR imaging of tumors. In the context of functional and physiological characterization, this implied transferring techniques established in tumor model research to human subjects and assessing their feasibility for use in patients. In the context of the morphologic assessment and parameter imaging of tumors, novel concepts and techniques were developed, which facilitated the simultaneous quantification of multiple MR parameters, the generation of “synthetic” MR images with various contrasts, and the fast single-shot acquisition of purely T2-weighted images.
Das Ziel der Promotion war die Entwicklung eines Flusssensors mit dem Fokus auf Flussmessungen an Pflanzen. Dazu musste zunächst die Hardware in Form eines räumlich zugänglichen Magneten und einer Sende- und Empfangseinheit entworfen werden. Um die MR-Konsole ansteuern zu können, musste eine Software entwickelt werden. Die AC-Methode wurde für Flussmessungen mit niedrigen Geschwindigkeiten angepasst und die entsprechende Theorie dazu erweitert. Mit dieser weiterentwickelten AC-Methode wurde die Flussmessung an Pflanzen demonstriert. Dafür wurden im Rahmen einer Kooperation mit der Arbeitsgruppe „Lipid Motobolism“ der IPK-Gatersleben Flussstudien an Weizenpflanzen durchgeführt. Darüber hinaus wurde in dieser Arbeit eine neue Technik zur Wirbelstromvermeidung bei Permanentmagneten entwickelt, um Problemen mit diesen bei Flussmessungen entgegenzuwirken.
Sensorbau
Es wurde ein zugänglicher, mobiler Magnet mit einer Feldstärke von 0,42 T gebaut. Die Feldhomogenität beträgt 0,5 ppm in 1 cm³. Im Vergleich zu dem am Lehrstuhl der EP5 bestehenden, geschlossenen, mobilen Magnetsystem erreicht das in dieser Arbeit gebaute System ein 40fach homogeneres Magnetfeld. Erzielt wurden diese Verbesserungen durch ein spezielles Design, welches durch Computersimulationen sukzessiv optimiert wurde. Durch angepasste Polschuhe konnte darüber hinaus die Induktion von Wirbelströmen im Mittel um einen Faktor 7 reduziert werden, wodurch phasensensitive Flussmessungen ermöglicht wurden.
Um die Zugänglichkeit zu dem Innenraum der HF-Spulen zu gewährleisten, wurde eine Klappspule weiterentwickelt und als Sende- und Empfangseinheit für den Tomographen gebaut. Ferner wurde ein System gebaut, dass direkt um die Pflanze gewickelt werden kann und sich somit für besonders dünne Pflanzenstängel eignet.
Weiterhin wurden zwei Systeme zur Rauschunterdrückung für die Messungen an Pflanzen entwickelt. Dadurch konnte das Rauschen um einen Faktor 92 gesenkt werden. Dies war notwendig, weil die länglichen Pflanzen durch ihre Ausdehnung über das Gehäuse hinweg ein Rauschen in die Empfangsspule induziert haben. Die beiden Rauschunterdrückungssysteme, die elektrische Schirmung und die Gleichtaktunterdrückung, entfernten das Rauschen dabei gleichermaßen.
Flussmessung
Die im Rahmen der Arbeit erfolgte Weiterentwicklung der AC-Methode [102] erlaubte es erstmals mit der Methode quantitative Flussprofile aufzunehmen. In Folge dessen war es
außerdem möglich Geschwindigkeiten unter 200 µm/s zu messen. Die Vorrausetzung dafür war die Implementierung von trapezförmigen Gradienten, welche kürzere Rampzeiten und eine stärkere Kodierung zulassen. Dadurch sind außerdem Intervalle ohne Gradienten realisierbar, die effizientere Refokussierungspulse und die Aufnahme mehrerer Datenpunkte ermöglichen. Die zu erwartenden und simulierten Flussprofile entsprachen den gemessenen Profilen durch die Verwendung einer neuen Auswertungstechnik.
Die neu entwickelte Erweiterung zur Bildgebung ermöglicht die ortsaufgelöste, spektroskopische Flussmessung und so können die Bereiche von Xylem und Phloem voneinander getrennt werden. Dies wurde durch Messungen einer Schwarzerle gezeigt, bei der die im Abschnitt 5.1 beschriebene Struktur dikotyler Pflanzen aufgelöst werden konnte. Zusätzlich können qualitativ genauere Aussagen über die Flussgeschwindigkeit getroffen werden.
Bei Messungen an Pflanzen konnte mit der optimierten AC-Methode die Flussänderungen aufgrund äußerer Einflüsse, wie der Beleuchtung, beobachtet werden. Langzeitmessungen über 9 Tage zeigten einen der Beleuchtung folgenden Flussverlauf - auch bei sehr geringen mittleren Flussänderungen von unter 200 µm/s.
Bloch-Siegert Phasenkodierung
Um eine Phasenkodierung ohne die Induktion von Wirbelströmen zu erhalten, wurde im Rahmen der Arbeit die ortsabhängige Phasenkodierung mittels B1-Gradienten entwickelt. Diese Technik basiert auf HF-Wechselfeldern und benutzt den sogenannten BS-Shift um einen B1-feldabhängigen Frequenzshift zu induzieren. Zwei Rekonstruktionstechniken wurden entwickelt, um die Rekonstruktion von entzerrten Bildern zu ermöglichen. Dies war notwendig, da die Kodierung mittels BS-Shift von B1² abhängt. Infolgedessen wird bei der Verwendung von konstanten HF-Gradienten eine vom Quadrat des Ortes abhängige Phasenkodierung induziert. Als Alternative zu diesem Verfahren wurde ein Gradient entwickelt, der einen wurzelförmigen Feldverlauf hat und somit die lineare Kodierung ohne angepasste Rekonstruktionstechniken ermöglicht.
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.
Das Ziel dieser Promotion ist der Aufbau eines dreMR Setups für einen klinischen 1,5T Scanner, das die Relaxations-Dispersions-Bildgebung ermöglicht, und die anschließende Ergründung von möglichst vielen Anwendungsfeldern von dreMR. Zu der Aufgabe gehört die Bereitstellung der zugrunde liegenden Theorie, der Bau des experimentellen Setups (Offset-Spule und Stromversorgung) sowie die Programmierung der nötigen Software. Mit dem gebauten Setup konnten zwei große Anwendungsfelder — dreMR Messungen mit und ohne Kontrastmitteln — untersucht werden.
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.