Refine
Has Fulltext
- yes (54)
Is part of the Bibliography
- yes (54) (remove)
Year of publication
Document Type
- Doctoral Thesis (54)
Keywords
- Kernspintomografie (54) (remove)
Institute
- Physikalisches Institut (28)
- Graduate School of Life Sciences (12)
- Institut für diagnostische und interventionelle Radiologie (Institut für Röntgendiagnostik) (7)
- Medizinische Klinik und Poliklinik I (5)
- Institut für diagnostische und interventionelle Neuroradiologie (ehem. Abteilung für Neuroradiologie) (2)
- Deutsches Zentrum für Herzinsuffizienz (DZHI) (1)
- Graduate School of Science and Technology (1)
- Institut für Molekulare Infektionsbiologie (1)
- Klinik und Poliklinik für Allgemein-, Viszeral-, Gefäß- und Kinderchirurgie (Chirurgische Klinik I) (1)
- Klinik und Poliklinik für Psychiatrie, Psychosomatik und Psychotherapie (1)
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.
Herzkreislauferkrankungen stellen die häufigsten Todesursachen in den Industrienationen dar. Die Entwicklung nichtinvasiver Bildgebungstechniken mit Hilfe der Magnetresonanz-Tomografie (MRT) ist daher von großer Bedeutung, um diese Erkrankungen frühzeitig zu erkennen und um die Entstehungsmechanismen zu erforschen. In den letzten Jahren erwiesen sich dabei genetisch modifzierte Mausmodelle als sehr wertvoll, da sich durch diese neue Bildgebungsmethoden entwickeln lassen und sich der Krankheitsverlauf im Zeitraffer beobachten lässt.
Ein große Herausforderung der murinen MRT-Bildgebung sind die die hohen Herzraten und die schnelle Atmung. Diese erfordern eine Synchronisation der Messung mit dem Herzschlag und der Atmung des Tieres mit Hilfe von Herz- und Atemsignalen. Konventionelle Bildgebungstechniken verwenden zur Synchronisation mit dem Herzschlag EKG Sonden, diese sind jedoch insbesondere bei hohen Feldstärken (>3 T) sehr störanfällig. In dieser Arbeit wurden daher neue Bildgebungsmethoden entwickelt, die keine externen Herz- und Atemsonden benötigen, sondern das MRT-Signal selbst zur Bewegungssynychronisation verwenden. Mit Hilfe dieser Technik gelang die Entwicklung neuer Methoden zur Flussbildgebung und der 3D-Bildgebung, mit denen sich das arterielle System der Maus qualitativ und quantitativ erfassen lässt, sowie einer neuen Methode zur Quantisierung der longitudinalen Relaxationszeit T1 im murinen Herzen. Die in dieser Arbeit entwickelten Methoden ermöglichen robustere Messungen des Herzkreislaufsystems. Im letzten Kapitel konnte darüber hinaus gezeigt werden dass sich die entwickelten Bildgebungstechniken in der Maus auch auf die humane Bildgebung übertragen lassen.
In dieser Arbeit wurden für spezielle Anwendungen an klinischen MR-Geräten optimierte Phased-Array-Spulen entwickelt. Das Ziel war, durch die Verwendung neuer Spulen entweder neue Anwendungsgebiete für klinische MR-Geräte zu eröffnen oder bei bestehenden Applikationen die Diagnosemöglichkeiten durch eine Kombination von höherem SNR und kleineren g-Faktoren im Vergleich zu bestehenden Spulen zu verbessern. In Kapitel 3 wurde untersucht, ob es durch den Einsatz neu entwickelter, dedizierter Kleintierspulen sinnvoll möglich ist, Untersuchungen an Kleintieren an klinischen MR-Geräten mit einer Feldstärke von 1,5T durchzuführen. Der Einsatz dieser Spulen verspricht dem klinischen Anwender Studien an Kleintieren durchführen zu können, bei denen er den gleichen Kontrast wie bei einer humanen Anwendung erhält und gleichzeitig Kontrastmittel sowie Sequenzen, die klinisch erprobt sind, einzusetzen. Durch die gewählten geometrischen Abmessungen der Spulen ist es möglich, Zubehör von dedizierten Tier-MR-Geräten, wie z. B. Tierliegen oder EKG- bzw. Atemtriggereinheiten, zu verwenden. Durch Vorversuche an für Ratten dimensionierten Spulen wurden grundlegende Zusammenhänge zwischen verwendetem Entkopplungsmechanismus und SNR bzw. Beschleunigungsfähigkeit erarbeitet. Für Ratten wurde gezeigt, dass in akzeptablen Messzeiten von unter fünf Minuten MR-Messungen des Abdomens in sehr guter Bildqualität möglich sind. Ebenfalls gezeigt wurde die Möglichkeit durch den Einsatz von paralleler Bildgebung sowie Kontrastmitteln hochaufgelöste Angiographien durchzuführen. Es stellte sich heraus, dass bei 1,5T dedizierte Mäusespulen bei Raumtemperatur von den SNR-Eigenschaften am Limit des sinnvoll Machbaren sind. Trotzdem war es möglich, auch für Mäuse ein 4-Kanal-Phased-Array zu entwickeln und den Einsatz bei kontrastmittelunterstützten Applikationen zu demonstrieren. Insgesamt wurde gezeigt, dass durch den Einsatz von speziellen, angepassten Kleintierspulen auch Tieruntersuchungen an klinischen MR-Geräten mit niedriger Feldstärke durchführbar sind. Obwohl sich die Bestimmung der Herzfunktion an MR-Geräten im klinischen Alltag zum Goldstandard entwickelt hat, ist die MR-Messung durch lange Atemanhaltezyklen für einen Herzpatienten sehr mühsam. In Kapitel 4 wurde deswegen die Entwicklung einer 32-Kanal-Herzspule beschrieben, welche den Komfort für Patienten deutlich erhöhen kann. Schon mit einem ersten Prototypen für 3T war es möglich, erstmals Echtzeitbildgebung mit leicht reduzierter zeitlicher Auflösung durchzuführen und damit auf das Atemanhalten komplett zu verzichten. Dies ermöglicht den Zugang neuer Patientengruppen, z. B. mit Arrythmien, zu MR-Untersuchungen. Durch eine weitere Optimierung des Designs wurde das SNR sowie das Beschleunigungsvermögen signifikant gesteigert. Bei einem Beschleunigungsfaktor R = 5 in einer Richtung erhält man z. B. gemittelt über das gesamte Herz ein ca. 60 % gesteigertes SNR zu dem Prototypen. Die Kombination dieser Spule zusammen mit neuentwicklelten Methoden wie z. B. Compressed- Sensing stellt es in Aussicht, die Herzfunktion zukünftig in der klinischen Routine in Echtzeit quantifizieren zu können. In Kapitel 5 wurde die Entwicklung einer optimierten Brustspulen für 3T beschrieben. Bei Vorversuchen bei 1,5T wurden Vergleiche zwischen der Standardspule der Firma Siemens Healthcare und einem 16-Kanal-Prototypen durchgeführt. Trotz größerem Spulenvolumen zeigt die Neuentwicklung sowohl hinsichtlich SNR als auch paralleler Bildgebungseigenschaften eine signifikante Verbesserung gegenüber der Standardspule. Durch die Einhaltung aller Kriterien für Medizinprodukte kann diese Spule auch für den klinischen Einsatz verwendet werden. Mit den verbesserten Eigenschaften ist es beispielsweise möglich, bei gleicher Messdauer eine höhere Auflösung zu erreichen. Aufgrund des intrinsischen SNR-Vorteils der 3 T-Spule gegenüber der 1,5 T-Spule ist es dort sogar möglich, bei höheren Beschleunigungsfaktoren klinisch verwertbare Schnittbilder zu erzeugen. Zusammenfassend wurden für alle drei Applikationen NMR-Empfangsspulen entwickelt, die im Vergleich zu den bisher verfügbaren Spulen, hinsichtlich SNR und Beschleunigungsvermögen optimiert sind und dem Anwender neue Möglichkeiten bieten.
Ziel dieser Arbeit war es, das gesamte Lungenvolumen in hoher dreidimensionaler Auflösung mittels der MRT darzustellen. Um trotz der niedrigen Protonendichte der Lunge und der geforderten hohen Auflösung ausreichend Signal für eine verlässliche Diagnostik zu erhalten, sind Aufnahmezeiten von einigen Minuten nötig. Um die Untersuchung für den Patienten angenehmer zu gestalten oder auf Grund der eingeschränkten Fähigkeit eines Atemstopps überhaupt erst zu ermöglichen, war eine Anforderung, die Aufnahmen in freier Atmung durchzuführen. Dadurch entstehen allerdings Bewegungsartefakte, die die Diagnostik stark beeinträchtigen
und daher möglichst vermieden werden müssen. Für eine Bewegungskompensation
der Daten muss die auftretende Atembewegung detektiert werden. Die Bewegungsdetektion
kann durch externe Messgeräte (Atemgurt oder Spirometer) oder durch eine
zusätzliche Anregungen erfolgen (konventionelle Navigatoren) erfolgen. Nachteile
dieser Methoden bestehen darin, dass die Bewegung während der Atmung nicht
direkt verfolgt wird, dass elektronische Messgeräte in die Nähe des Tomographen
gebracht werden und das die Patienten zusätzlich vorbereitet und eingeschränkt
werden. Des Weiteren erfordert eine zusätzliche Anregung extra Messzeit und kann
unter Umständen die Magnetisierung auf unterwünschte Weise beeinflussen.
Um die angesprochenen Schwierigkeiten der Bewegungsdetektion zu umgehen,
wurden in dieser Arbeit innerhalb einer Anregung einer 3d FLASH-Sequenz sowohl
Bilddaten- als auch Navigatordaten aufgenommen. Als Navigator diente dabei das
nach der Rephasierung aller bildgebenden Gradienten entstehende Signal (DC Signal).
Das DC Signal entspricht dabei der Summe aller Signale, die mit einem bestimmten
Spulenelement detektiert werden können. Bewegt sich beispielsweise die Leber
bedingt durch die Atmung in den Sensitivitätsbereich eines Spulenelementes, wird
ein stärkeres DC Signal detektiert werden. Je nach Positionierung auf dem Körper
kann so die Atembewegung mit einzelnen räumlich lokalisierten Spulenelementen
nachverfolgt werden. Am DC Signalverlauf des für die Bewegungskorrektur ausgewählten
Spulenelementes sind dann periodische Signalschwankungen zu erkennen.
Zusätzlich können aus dem Verlauf Expirations- von Inspirationszuständen unterschieden
werden, da sich Endexpirationszustände im Regelfall durch eine längere
Verweildauer auszeichnen.
Grundsätzlich kann das DC Signal vor oder nach der eigentlichen Datenaufnahme
innerhalb einer Anregung aufgenommen werden. Auf Grund der kurzen Relaxationszeit
T∗2 des Lungengewebes fällt das Signal nach der RF Anregung sehr schnell ab. Um
möglichst viel Signal zu erhalten sollten, wie in dieser Arbeit gezeigt wurde, innerhalb einer Anregung zuerst die Bilddaten und danach die Navigatordaten aufgenommen
werden. Dieser Ansatz führt zu einer Verkürzung der Echozeit TE um 0.3 ms und
damit zu einem SNR Gewinn von etwa 20 %. Gleichzeitig ist das verbleibende Signal
nach der Datenakquisition und Rephasierung der bildgebenden Gradienten noch
ausreichend um die Atembewegung zu erfassen und somit eine Bewegungskorrektur
der Daten (Navigation) zu ermöglichen.
Um eine retrospektive Bewegungskorrektur durchführen zu können, müssen Akzeptanzbedingungen
(Schwellenwerte) für die Datenauswahl festgelegt werden. Bei
der Wahl des Schwellenwertes ist darauf zu achten, dass weder zu wenige noch zu
viele Daten akzeptiert werden. Akzeptiert man sehr wenige Daten, zeichnen sich die
Rekonstruktionen durch einen scharfen Übergang zwischen Lunge und Diaphragma
aus, da man sehr wenig Bewegung in den Rekonstruktionen erlaubt. Gleichzeitig
erhöht sich allerdings das Risiko, dass nach der Navigation Linien fehlen. Dies führt
zu Einfaltungsartefakten, die in Form von gestörten Bildintensitäten in den Rekonstruktionen
zu sehen sind und die diagnostische Aussagekraft einschränken. Um
Einfaltungsartefakte zu vermeiden sollte der Schwellenwert so gewählt werden, dass
nach der Datenauswahl keine Linien fehlen. Aus dieser Anforderung lässt sich ein
maximaler Schwellenwert ableiten. Akzeptiert man dagegen sehr viele Daten, zeichnen
sich die Rekonstruktionen durch erhöhtes Signal und das vermehrte Auftreten
von Bewegungsartefakten aus. In diesem Fall müsste der Arzt entscheiden, ob Bewegungsartefakte
die Diagnostik zu stark beeinflussen. Wählt man den Schwellenwert
so, dass weder Linien fehlen noch zu viel Bewegung erlaubt wird, erhält man Rekonstruktionen
die sich durch einen scharfen Diaphragmaübergang auszeichnen und in
denen noch kleinste Gefäße auch in der Nähe des Diaphragmas deutlich zu erkennen
sind. Hierfür haben sich Schwellenwerte, die zu einer Datenakzeptanz von ca. 40 %
führen als günstig erwiesen.
Um Einfaltungsartefakte auf Grund der retrospektiven Datenauswahl zu verhindern,
muss das Bildgebungsvolumen mehrfach abgetastet werden. Dadurch wird
gewährleistet, dass für die letztendliche Rekonstruktion ausreichend Daten zur Verfügung
stehen, wobei mehrfach akzeptierte Daten gemittelt werden. Dies spielt auf
Grund der niedrigen Protonendichte der Lunge eine wesentliche Rolle in der Rekonstruktion
hochaufgelöster Lungendatensätze. Weiterhin führt das Mitteln von
mehrfach akzeptierten Daten zu einer Unterdrückung der sogenannten Ghost Artefakte,
was am Beispiel der Herzbewegung in der Arbeit gezeigt wird.
Da die Messungen unter freier Atmung durchgeführt werden und keine zusätzlichen
externen Messgeräte angeschlossen werden müssen, stellte die Untersuchung
für die Patienten in dieser Arbeit kein Problem dar. Im ersten Teil dieser wurde Arbeit
gezeigt, dass sich mit Hilfe des DC Signales als Navigator und einer retrospektiven
Datenauswahl das gesamte Lungenvolumen in hoher dreidimensionaler Auflösung
von beispielsweise 1.6 x 1.6 x 4 mm3 innerhalb von 13 min. darstellen lässt. Die Anwendbarkeit der vorgestellten Methode zur Bewegungskorrektur wurde neben
Probanden auch an Patienten demonstriert.
Da wie bereits beschrieben das Bildgebungsvolumen mehrfach abgetastet werden
muss, wiederholt sich auch die Abfolge der für die Bildgebung verantwortlichen
Gradienten periodisch. Da sich der Atemzyklus aber auch periodisch wiederholt,
kann es zu Korrelationen zwischen der Atmung und den wiederholten Messungen
kommen. Dies führt dazu, dass auch nach vielen wiederholten Messungen immer
noch größere Bereiche fehlender Linien im k-Raum bleiben, was zu Artefakten in
den Rekonstruktionen führt. Dies konnte im Falle der konventionellen Bewegungskorrektur
in den Gatingmasken, die die Verteilung und Häufigkeit der einzelnen
akzeptierten Phasenkodierschritte im k-Raum zeigen, beobachtet werden.
Da eine vorsätzliche Unterbrechung der Atemperiodizität (der Patient wird dazu
angehalten, seine Atemfrequenz während der Messung absichtlich zu variieren) zur
Vermeidung der angesprochenen Korrelationen nicht in Frage kommt, musste die
Periodizität in der Datenaufnahme unterbrochen werden. In dieser Arbeit wurde
dies durch eine quasizufällige Auswahl von Phasen- und Partitionskodiergradienten
erreicht, da Quasizufallszahlen so generiert werden, dass sie unabhängig von ihrer
Anzahl einen Raum möglichst gleichförmig ausfüllen. Die quasizufällige Datenaufnahme
führt deshalb dazu, das sowohl akzeptierte als auch fehlende Linien nach der
Bewegungskorrektur homogen im k-Raum verteilt auftreten.
Vergleicht man das auftreten von Ghosting zeichnen sich die quasizufälligen Rekonstruktionen
im Vergleich zur konventionellen Datenaufnahme durch eine verbesserte
Reduktion von Ghost Artefakten aus. Dies ist auf die homogene Verteilung mehrfach
akzeptierter Linien im k-Raum zurückzuführen. Die homogenere Verteilung von
fehlenden Linien im k-Raum führt weiterhin zu einer wesentlich stabileren Rekonstruktion
fehlender Linien mit parallelen MRT-Verfahren (z.B. iterativem Grappa).
Dies wird umso deutlicher je höher der Anteil fehlender Linien im k-Raum wird. Im
Falle der konventionellen Datenaufnahme werden die zusammenhängenden Bereiche
fehlender Linien immer größer, was eine erfolgreiche Rekonstruktion mit iterativem
Grappa unmöglich macht. Im Falle der quasizufälligen Datenaufnahme dagegen
können auch Datensätze in denen 40% der Linien fehlen einfaltungsartefaktfrei
rekonstruiert werden.
Im weiteren Verlauf der Arbeit wurde gezeigt, wie die Stabilität der iterativen Grappa
Rekonstruktion im Falle der quasizufälligen Datenaufnahme für eine erhebliche
Reduktion der gesamten Messzeit genutzt werden kann. So ist in einer Messzeit von
nur 74s die Rekonstruktion eines artefaktfreien und bewegungskorrigierten dreidimensionalen
Datensatzes der menschlichen Lunge mit einer Auflösung von 2 x 2 x
5 mm3 möglich. Des Weiteren erlaubt die quasizufällige Datenaufnahme in Kombination
mit iterativem Grappa die Rekonstruktion von Datensätzen unterschiedlicher
Atemphasen von Inspiration bis Expiration (4D Bildgebung). Nach einer Messzeit
von 15min. wurden 19 unterschiedliche Atemzustände rekonstruiert, wobei sich der Anteil der fehlenden Linien zwischen 0 und 20 % lag. Im Falle der konventionellen
Datenaufnahme wäre eine wesentlich längere Messzeit nötig gewesen, um ähnliche
Ergebnisse zu erhalten.
Zum Schluss soll noch ein Ausblick über mögliche Weiterentwicklungen und Anwendungsmöglichkeiten,
die sich aus den Erkenntnissen dieser Arbeit ergeben haben,
gegeben werden. So könnte das quasizufällige Aufnahmeschema um eine Dichtegewichtung
erweitert werden. Hierbei würde der zentrale k-Raum Bereich etwas
häufiger als die peripheren Bereiche akquiriert werden. Dadurch sollte die iterative
Grappa Rekonstruktion noch stabiler funktionieren und Ghost Artefakte besser reduziert
werden. Die Verteilung der Linien sollte allerdings nicht zu inhomogen werden,
um größere Lücken im k-Raum zu vermeiden.
Darüber hinaus könnte die vorgestellte Methode der Bewegungskompensation
auch für die Untersuchung anderer Organe oder Körperteile verwendet werden.
Voraussetzung wäre lediglich das Vorhandensein dezidierter Spulenanordnungen,
mit denen die Bewegung nachverfolgt werden kann. So ist beispielsweise eine dynamische
Bildgebung des frei und aktiv bewegten Knies möglich, wobei zwischen
Beugung und Streckung durch die erste Ableitung des zentralen k-Raum Signales
unterschieden werden kann. Dies kann zusätzliche Diagnoseinformationen liefern
oder für Verlaufskontrollen nach Operationen benutzt werden [15].
Eine Weiterentwicklung mit hohem klinischen Potential könnte die Kombination
der in dieser Arbeit vorgestellten retrospektiven Bewegungskorrektur mit einer Multi-
Gradienten-Echo Sequenz darstellen. Hierzu musste die bestehende Sequenz lediglich
um eine mehrfache Abfolge von Auslesegradienten innerhalb einer Anregung erweitert
werden. Dies ermöglicht eine bewegungskorrigierte voxelweise Bestimmung der
transversalen Relaxationszeit T∗2 in hoher räumlicher Auflösung. Unter zusätzlicher
Sauerstoffgabe kann es zu einer Veränderung von T∗2 kommen, die auf den sogenannten
BOLD Effekt (Blood Oxygen Level Dependent) zurückzuführen ist. Aus dieser
Änderung könnten Rückschlüsse auf hypoxische Tumorareale gezogen werden. Da
diese eine erhöhte Strahlenresistenz aufweisen, könnte auf diese Bereiche innerhalb
des Tumors eine erhöhte Strahlendosis appliziert und so möglicherweise Behandlungsmisserfolge
reduziert werden. Gleichzeitig kann durch die 4D Bildgebung eine
mögliche Tumorbewegung durch die Atmung erfasst und diese Information ebenfalls
in der Bestrahlungsplanung benutzt werden. Die Lungen MRT könnte somit um eine
hochaufgelöste dreidimensionale funktionelle Bildgebung erweitert 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.
Diffusionstensorbildgebung im Vergleich zu anderen Parametermethoden für die Infarktcharakterisierung Ziel dieses Teils der Arbeit war die Klärung der Frage, welches Potential verschiedene MR-Parametersequenzen bei der Charakterisierung eines myokardialen Infarkts sowohl im akuten als auch im chronischen Fall haben. Dazu wurde eine Studie mit akut und chronisch infarzierten Rattenherzen durchgeführt. Untersucht wurden die Parameter T1, T2 und T2* sowie die aus der Diffusionstensorbildgebung berechneten Parameter ADC, FA, cs, cp und cl . Es zeigte sich, dass es kein Analogon zum bei einer cerebralen Ischämie bekannten Mismatch-Konzept gibt. Weder im akuten noch im chronischen war Fall eine ausgewiesene Differenz im diagnostizierten Infarktareal zwischen verschiedenen Sequenzen feststellbar. Alles in allem eignen sich zur detaillierten Charakterisierung der Infarktnarbe am besten eine T2*- oder eine Diffusionstensorsequenz. Die T2*-Sequenz liefert optisch das aufschlussreichere Bild, die aufwendigere Diffusionstensorsequenz dagegen bietet aufgrund der vielfachen Darstellungsmöglichkeiten im Postprocessing ein Mehr an Information und zeigt dazu eine Veränderung der Narbe im Zeitverlauf. Oxygenierungsmessung am Mäuseherz in vivo Die Charakterisierung einer Infarktnarbe kann auch über die Darstellung morphologischer Strukturen hinaus erfolgen. Die Oxygenierung ist ein komplexer Parameter, der funktionelle Auskunft über die Vaskularisierung und Viabilität des Gewebes geben kann. Zugang zu diesem Parameter erhält man über T2*-Messungen, da der Parameter T2* sensitiv auf chemisch gebundenen Sauerstoff reagiert. Hier wurden der Einfluss von reiner Sauerstoffatmung im Gegensatz zu normaler Raumluftatmung auf die Oxygenierung bei gesunden und infarzierten Mäusen untersucht. Die Messungen wurden trotz der Schwierigkeiten, die durch die Bewegung durch Atmung und Herzschlag entstehen, in vivo bei 17,6 Tesla implementiert und durchgeführt. Die Auflösung war ausreichend, um auch nach Infarkt extrem ausgedünnte Myokardwände gut auflösen und charakterisieren zu können. Der Effekt auf das Oxygenierungslevel ist stark unterschiedlich zwischen normalen und infarzierten Herzen, woraus auf eine noch nicht weit fortgeschrittene Revaskularisierung der Narbe eine Woche nach Infarzierung geschlossen werden kann. Die Methode wurde darüber hinaus an einem 7,0 Tesla-Magneten zur Verwendung an Ratten implementiert und auf das im Gegensatz zur Maus veränderte Atmungsverhalten der Ratte angepasst. Zum einen kann dadurch der Einfluss des hohen Magnetfeldes auf die Oxygenierungsmessung untersucht werden, zum anderen ist das Herz als zu untersuchendes Objekt bei der Ratte größer. Diffusionswichtung mittels Hole-Burning Die in dieser Arbeit zur Charakterisierung des Herzens verwendete Diffusionsmethode kann im Grenzfall von kurzen T2-Relaxationszeiten an ihre Grenzen stoßen: Bei den verwendeten starken Magnetfeldern klingt das messbare Signal aufgrund der Relaxationszeit T2 oft sehr schnell ab. Daher wurde eine Methode entwickelt, die einen völlig neuen Ansatz zur diffusionsgewichteten Bildgebung verfolgt, bei dem die Informationen über die Diffusion unabhängig von der limitierenden T2-Zeit gewonnen werden können. Die sog. Hole-Burning-Diffusionssequenz verwendet in einem Vorexperiment lediglich die Longitudinalmagnetisierung zur Diffusionswichtung. Das Signal wird dann mit einer schnellen Auslesesequenz akquiriert. Bei der Präparation werden zunächst auf Subvoxel-Niveau Streifen "gebrannt", d.h. die Magnetisierung wird dort gesättigt. Bis zur nächsten Sättigung ist das Verhalten der Magnetisierung abhängig von der T1-Relaxation in diesem Bereich und vom Diffusionsverhalten. Durch rasches Wiederholen des selektiven Pulszugs wird schließlich eine Gleichgewichtsmagnetisierung erreicht, die von der Diffusionskonstanten D und der T1-Relaxationszeit abhängt. Im Rahmen dieser Arbeit wurden die Abhängigkeiten verschiedener Sequenzparameter untersucht und diese mittels Simulationen optimiert. Außerdem wurde die Sequenz an einem Scanner implementiert und erste Experimente damit durchgeführt. Mit Hilfe von Simulationen konnten dazu Lookup-Tabellen generiert werden, mit denen in bestimmten Bereichen (insbesondere bei nicht zu kurzen T1-Relaxationszeiten) sowohl die Diffusionskonstante D als auch die T1-Relaxationszeit quantifiziert werden konnte.
In dieser Arbeit werden für die Anwendung in der menschlichen Lunge
optimierte Methoden zur Bestimmung von T1- und T2*-Karten diskutiert:
Dc-Gating ermöglicht die Quantifizierung in freier Atmung, wobei für die
T1-Quantifizierung mittels Inversion Recovery eine Korrektur des dc-Signals
entwickelt wurde. Dies hat den Vorteil, dass Parameterkarten aus mehreren
Messungen anhand ihrer dc-Signale passend überlagert werden können. Da T1
und T2* auf unterschiedliche Art und Weise von der Sauerstoffkonzentration
abhängen, verbessert dies die Möglichkeit, ΔT1- und ΔT2*- Differenzkarten aus
Messungen mit unterschiedlichen O2-Konzentrationen im Atemgas zu erstellen.
Die Parameterquantifizierung ist in erster Linie für die Beobachtung von
Krankheitsverläufen interessant, da T1 und T2* absolute, vergleichbare Zahlen
sind. Da T2* deutlich vom Atemzustand abhängt, ist es auch hierfür sinnvoll,
durch Gating identische Atemzustände abzubilden. Um die unterschiedlichen
Einflüsse des Sauerstoffs auf T1 und T2* besser vergleichbar zu machen, wurde
in dieser Arbeit weiterhin eine kombinierte Messung für beide Parameter
implementiert: Da auch diese in freier Atmung stattfindet, profitieren nicht
nur die Differenzkarten von der Überlagerung der Bilder, sondern auch der
Vergleich der ΔT1- und ΔT2*-Karten untereinander.
Messungen mit einer konventionellen kartesischen Methode an COPD-Patienten
unter Raumluft- und 100% Sauerstoffatmung ergaben bei Verwendung identischer
Atemmasken ein deutlich geringeres ΔT1 als in gesunden Probanden. Dass T1 in
der Lunge nicht nur von der Sauerstoffkonzentration sondern auch von der
Gewebezusammensetzung und insbesondere auch dem Blutvolumenanteil abhängt,
zeigte sich hierbei aber auch an den bei COPD im Mittel sehr viel kürzeren
T1-Zeiten bei Raumluft. Die aufgrund emphysematischer Veränderung noch
zusätzlich reduzierte Protonendichte im Parenchym kranker Lungen macht diese
Messungen allerdings besonders schwierig.
Die oben erwähnten Optimierungen der T1-Quantifizierung zielen daher auch
darauf ab, das Signal aus der Lunge zu maximieren, um Patientenmessungen
einfacher zu machen: Messungen in freier Atmung sind für Patienten nicht nur
einfacher, sondern erlauben effektiv auch längere Messzeiten. Insbesondere
wurde aber durch die Entwicklung einer radialen Methode die Echozeit zur
Messung reduziert, um die kurze T2*-Zeit in der Lunge auszugleichen.
Schließlich wurde durch Implementation einer 2D UTE Sequenz die Messung bei
der kürzesten vom Scanner erlaubten Echozeit ermöglicht.
Die Messungen bei ultrakurzen Echozeiten in Probanden zeigten allerdings
deutlich kürzere T1-Zeiten als die zuvor gefundenen oder in der Literatur
dokumentierten. In weiteren Experimenten wurde das sichtbare T1 zu mehreren
Echozeiten mit Hilfe der zur kombinierten Quantifizierung entwickelten
Methode bestimmt. Dabei ergab sich eine Zunahme des gemessenen T1 mit der
Echozeit. Aus diesem Verhalten sowie den gefundenen kürzesten und längsten T1
lässt sich schließen, dass das intra- und extravaskuläre Lungenwasser, also
Blut bzw. das umgebende Gewebe, mit unterschiedlichen T1- und T2*-Zeiten zum
Signal und damit auch dem effektiven T1 beitragen.
Dass das TE der Messung die Gewichtung dieser Kompartimente bestimmt, hat
dabei mehrere Auswirkungen: Einerseits bedeutet dies, dass beim Vergleich von
T1-Messungen in der Lunge stets auch das TE mitbetrachtet werden muss, bei
dem diese durchgeführt wurden. Andererseits lässt sich die Möglichkeit, die
Messung auf die unterschiedlichen Kompartimente abzustimmen, potentiell
ausnutzen, um zusätzliche diagnostische Informationen zu gewinnen: Da T1 vom
Blutvolumenanteil und der Gewebezusammensetzung abhängt, könnte dieser Effekt
helfen, diese beiden Einflüsse zu differenzieren.
Während die in dieser Arbeit beschriebenen Experimente die TE-Abhängigkeit
des sichtbaren T1 in Probanden aufzeigen, liefern sie allerdings noch keine
genaue Erklärung für die möglichen Ursprünge dieses Effekts. Um diese weiter
zu untersuchen, könnten allerdings gezielte Phantom- und in vivo-Experimente
Aufschluss geben: Ein Aufbau, der die Feldverzerrung durch luftgefüllte
Alveolen in Lösungen mit entsprechenden verschiedenen Suszeptibilitäten
nachbildet, reduziert den Unterschied zwischen den Kompartimenten auf T1 und
χ. Eine in vivo-Messung mit möglichst großer Differenz zwischen Ex- und
Inspiration hingegen könnte den Einfluss der Abstände der Kompartimente vom
Gasraum aufzeigen, da die Alveolarwände in tiefer Inspiration am weitesten
gedehnt und daher am dünnsten sind.
Im Rahmen dieser Arbeit sollten die Möglichkeiten der MR Tomographie erkundet werden bakterielle Infektionen im Zeitverlauf darzustellen. Genauer gesagt sollte das Potential der MR Tomographie anhand eines durch eine Infektion induzierten lokalisierten Abszesses unter Verwendung dreier unterschiedlicher MRT Methoden untersucht werden: Mittels nativem \(T_2\) Kontrast; der Verwendung von superparamagnetischen Eisenoxid Partieln (USPIO) als \(T_2^*\) Kontrastmittel; und dem Einsatz von Perfluorkarbonen (PFC) als \(^{19}F\) MRT Marker (siehe Kapitel 3).
Wie erwartet führte die durch die Infektion hervorgerufene Entzündung zu veränderten \(T_2\)-Zeiten, welche auf \(T_2\)-gewichteten MR Bildern eine Lokalisierung des Abszessbereiches erlauben. Jedoch eigneten sich diese Daten aufgrund der graduellen Änderung der \(T_2\)-Zeiten nicht, um eine klare Grenze zwischen Abszess und umliegendem Gewebe zu ziehen.
Superparamagnetische Eisenoxidpartikel andererseit haben als MRT Kontrastmittel bereits in den letzten Jahren ihre Fähigkeit unter Beweis gestellt Entzündungen [53, 58, 64] darzustellen. Die Anreicherung dieser Partikel am Rande des Abszesses [53], wie sie auch in unseren MR Daten zu beobachten war, erlaubte eine relativ scharfe Abgrenzung gegenüber dem umgebenden Gewebe in der chronischen Phase der Infektion (Tag 9 p.i.). Hingegen genügte die nur sehr spärlichen Anreicherung von USPIO Partikeln in der akuten Phase der Infektion (Tag 3 p.i.) nicht für eine entsprechende Abgrenzung [58].
Aufgrund der sehr geringen biologischen Häufigkeit und den sehr kurzen Relaxationszeiten von endogenem Fluor eignen sich Perfluorkarbone als Markersubstanz in der MR Tomographie von biologischen Systemen. Insbesondere da PFC Emulsionen durch phagozytierende Zellen aufgenommen werden und im Bereich von Entzündungen akkumulieren [30, 59]. In dieser Arbeit konnte anhand der erhaltenen MRT Daten eine Akkumulation von Perfluorkarbonen nicht nur in der chronischen Phase, sondern auch in der akuten Phase nachgewiesen werden. Diese Daten erlauben somit zu allen untersuchten Zeitpunkten eine Abgrenzung zwischen Infektion und umliegenden Gewebe.
Aufgrund der besagten Vorteile wurden die Perfluorkarbone gewählt, um die Möglichkeiten der MR Tomographie zu testen, quantitative Informationen über die schwere der Infektion zu liefern. Als Referenz für die Bakterienbelastung wurden die Biolumineszenzbildgebung (BLI) [49, 50] und die Standardmethode zur Bestimmung der Bakterienbelastung cfu (koloniebildenden Einheiten) herangezogen. Eine Gegenüberstellung der zeitlichen Verläufe der durch die Biolumineszenzbildgebung und durch die cfu erhaltenen Daten liefert eine qualitative Übereinstimmung mit den durch die 19F MR Tomographie erhaltenen Daten. Dies trifft hierbei sowohl auf die über den gesamten Infektionsbereich hinweg summierten Signalamplituden, als auch auf das Volumen zu, in dem Fluor am Ort der Infektion akkumuliert wurde. Im Gegensatz zur Methode der cfu Bestimmung sind die MR Tomographie und die Biolumineszenzbildgebung nicht invasiv und erlauben die Verfolgung des Infektionsverlaufes an einem einzelnen Individuum. Hierzu benötigt, im Gegensatz zur MR Tomographie, die Methode der Biolumineszenzbildgebung jedoch einen speziellen Pathogenstamm. Darüber hinaus ist hervorzuheben, dass die MR Tomographie zudem die Möglichkeit bietet auch morphologische Informationen über den Infektionsbereich und seine Umgebung zu akquirieren.
Gerade weil jede dieser Methoden die mit der Infektion einhergehenden Prozesse aus einer leicht anderen Blickrichtung betrachtet, erscheint es sinnvoll diese etablierte Untersuchungsplattform bestehend aus MRT, BLI und cfu über die in dieser Arbeit bearbeitete Fragestellung hinaus näher zu untersuchen. Insbesondere der Aspekt inwieweit die drei Methoden sich gegenseitig ergänzen, könnte einen tieferen Einblick in die Wechselwirkung zwischen Pathogen und Wirt erlauben.
Auch wenn für die betrachtete Fragestellung bereits der hierdurchgeführte semiquanitative Ansatz zur Bestimmung der relativen Fluormengen am Ort der Infektion ausreichte, so ist doch im Allgemeinen wünschenswert probenbezogen die Sensitivität der Spule und damit die Güte der Spulenabstimmung zu bestimmen. Hierzu ist jedoch die Aufnahme von \(B_1\)-Karten unabdingbar und wird entsprechend im Kapitel 4 \(Bloch-Siegert B_1^+-Mapping\) näher addressiert. Der Schwerpunkt liegt hierbei, wie der Kapitelname bereits andeutet, auf der Bloch-Siegert Methode, die insbesondere in der präsentierten Implementierung in einer Turbo/ Multi Spin Echo Sequenz eine effiziente Nutzung der relativ langen \(T_\)2-Zeiten der Perfluorkarbone erlaubt. Da zudem die Bloch-Siegert-Methode eine rein phasenbasierte Methode ist, kann neben der aus den Daten erzeugten \(B_1\)-Karte zugleich ein unverfälschtes Magnitudenbild generiert werden, wodurch eine sehr effiziente Nutzung der vorhandenen Messzeit ermöglicht wird. Diese Eigenschaft ist insbesondere für \(^{19}F\) Bildgebung von besonderem Interesse, da hier für jede Messung, aufgrund der üblicherweise relativ geringen Konzentration an Fluoratomen, lange Messzeiten benötigt werden.
Zusammenfassend konnte anhand des untersuchten Tiermodells sowohl die Fähigkeit der MR Tomographie nachgewiesen werden Infektionen im Zeitverlauf darzustellen, als auch die Fähigkeit der MR Tomographie quantitative Informationen über den Verlauf der Infektion zu liefern. Desweiteren konnte eine Möglichkeit aufgezeigt werden, welche das Potential hat in vertretbarem Zeitrahmen auch in vivo B1+-Karten auf dem Fluorkanal zu erstellen und so einen zentralen Unsicherheitsfaktor, für Relaxometry und absolute Quantifizierung von \(^{19}F\) Daten in vivo, zu beseitigen.
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.
Was bedeutet eine Änderung der Kontrastmittelaufnahme in niedrig malignen Gliomen bei Kindern?
(2014)
Bei niedriggradigen Gliomen WHO I° und II° ist das Kontrastmittelverhalten variabel und nicht mit der Prognose assoziiert. Andererseits wurde bei spontanen Regressionen von niedriggradigen Gliomen bei Patienten sowohl mit als auch ohne eine Neurofibromatose Typ I eine Abnahme der Kontrastmittelaufnahme berichtet. Vielleicht deswegen verleitet eine Zunahme des Enhancements oder ein neu aufgetretenes Enhancement nicht selten zur Diagnose einer Tumorprogression. Es stellt sich also die Frage, ob eine Kontrastmittelaufnahme bei LGGs mit einer Größenzunahme somit einem Tumorprogress assoziiert ist.
Es lässt sich eine überdurchschnittliche Assoziation einer Änderung der Kontrastmittelaufnahme mit dem Wachstumsverhalten bei LGGs nachweisen. Allerdings ist das Maß für eine Progression oder Regression weiterhin die Größenänderung des Tumors. Auch eine neue Kontrastmittelaufnahme innerhalb eines Tumors hat keine Bedeutung für das aktuelle Staging. Sie findet sich zwar häufiger bei Wachstum aber auch bei regredienten Tumoren.
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.
Der Thoraxmagen beschreibt eine zirkuläre Schwachstelle der phrenikoösophagealen Membran mit einer schrittweisen Dislozierung der Magenkardia und des Ösophagus nach mediastinal. Die Therapie des Thoraxmagens kann konservativ im Sinne des „watchful waiting“ oder operativ erfolgen. Aufgrund der möglichen Komplikationen wird die elektive Operation durch die amerikanischen Leitlinien empfohlen. Ein zentrales Problem der Hiatushernienchirurgie stellt die hohe Anzahl an Rezidiven dar. Ob die Gründe hierfür in der Zwerchfellrekonstruktion, Speiseröhrenlänge, Fundoplicatio oder Netzaugmentation liegen, wird nach wie vor kontrovers diskutiert.
In dieser Arbeit wurde die operative Versorgung des Thoraxmagens von 124 Patienten des Universitätsklinikums Würzburg im Zeitraum von September 2008 bis Juni 2015 untersucht. Hierfür war neben den perioperativen Daten auch die Rezidiv- und Letalitätsrate von Relevanz. Das Patientenkollektiv wurde sowohl in Hinblick auf das Lebensalter als auch auf die verschiedenen Versorgungsarten analysiert. Um die postoperative Lebensqualität zu beurteilen, erfolgte die Patientenbefragung mit Hilfe eines Symptomfragebogens und dem Gastrointestinalen Lebensqualitätsindex nach Eypasch (GIQLI). Zusätzlich wurden 17 Patienten postoperativ mittels MRT untersucht, um eine optimierte MRT-Sequenz zur Beurteilung der Hiatusregion zu evaluieren.
Im Vergleich der Altersgruppen zeigte sich trotz einer erhöhten Komorbiditätsrate bei dem Patientenkollektiv ≥ 75 Jahre (p=0,002) kein signifikanter Unterschied bei Betrachtung der intraoperativen Komplikationen. Die Rezidivrate lag unabhängig vom Alter bei 20,2% im Untersuchungszeitraum, jedoch konnte eine verminderte Rezidivrate bei Patienten mit U-Shape Versorgung (p=0,015) festgestellt werden.
In der postoperativen Patientenbefragung zeigten sich 87,0% der Patienten, unabhängig vom Alter und der Versorgungsart, zufrieden mit dem Operationsergebnis und beschrieben ihren Zustand im Vergleich zu präoperativ als gebessert. Die Ergebnisse des GIQLI erbrachten in dem untersuchten Patientenkollektiv ein gegenüber der Allgemeinbevölkerung erniedrigten Wert mit 95,4 Punkten.
Die optimierte MRT-Sequenz zeichnete sich durch eine hohe diagnostische Konfidenz bei guter Bildqualität, kurzer Untersuchungsdauer und gleichzeitig hoher Akzeptanz der Patienten gegenüber dieser Art der Diagnostik aus.
Zusammenfassend stellt die operative Versorgung von Thoraxmägen, unabhängig des Patientenalters, eine sichere Therapieform dar, die zu einer hohen Patientenzufriedenheit führt. Die modifizierte MRT-Untersuchung hat sich als diagnostische Methode bewährt und stellt eine Alternative zu strahlenexponierenden oder von Seiten der Patienten weniger gut tolerierten Untersuchungsmodalitäten dar.
Zusammenfassung
Ziel:
Pathomorphologische Veränderungen der Kiefergelenke treten im Rahmen einer JIA häufig auf. In dieser Arbeit sollten anhand von einer kieferorthopädischen zweidimensionalen Röntgenbildgebung (Orthopantomographie - OPG) pathologische Veränderungen der Kiefergelenke bei JIA-Patienten differenziert und Asymmetrien des Unterkiefers in Abhängigkeit vom Grad der Kondylendestruktion bestimmt werden. Darüber hinaus sollte geprüft werden, wie sich die Krankheitsdauer auf den Befall der Kiefergelenke auswirkt.
Zusätzlich wurden zur differenziellen Analyse der Pathomorphologie die klinische Funktionsanalyse (FAL), die Joint Vibration Analysis (JVA) sowie die 3D-Stereophotogrammetrie (3d-Scan) eingesetzt. Ziel war mittels non-invasiver, kostengünstiger und schnell verfügbarer Untersuchungsmethoden eindeutige Parameter zum Befall der Kiefergelenke zu detektieren.
Patienten und Methodik:
In dieser Arbeit wurden 46 Patientin (28 weiblich; 18 männlich) kaukasischer Herkunft mit diagnostizierter JIA eingeschlossen. Die Kiefergelenke (n = 92) wurden einzeln nach dem Grad ihrer Kondylendestruktion (Grad 0 - 4 nach Billiau et. al. [78]) befundet und in eine leicht betroffene Gruppe 1 (Grad 0, 1 und 2 nach Billiau: röntgenologisch unauffällig, Kondylenerosionen, Kondylenabflachungen) und in eine schwer betroffene Gruppe 2 (Grad 3 und 4 nach Billiau: Kondylusabflachungen mit Erosionen, Komplettverlust des Kondylus) eingeteilt.
Zur Quantifizierung von Unterkieferasymmetrien wurde das Seitenverhältnis aus Kondylus-, Ramus- und Mandibulahöhe bestimmt. Der Vergleich der einzelnen klinischen, instrumentellen und bildgebenden Untersuchungsverfahren (OPG, FAL, JVA, 3d-Scan) erfolgte jeweils durch die Gegenüberstellung der schwer betroffenen und der leicht betroffenen Patientengruppe.
Ergebnisse:
Erkrankungsdauer:
Anhand des Grades der Kondylendestruktion wurden 36 Patienten in die leicht betroffene Gruppe 1 und 10 Patienten in die schwer betroffenen Gruppe 2 eingeteilt. Dabei war die Erkrankungsdauer in der schwer betroffenen Patientengruppe (8,9 ± 5,2 Jahre) signifikant länger, als in der leicht betroffenen Patientengruppe (4,6 ± 4,7; Jahre) (p = 0,031).
FAL:
Die Ergebnisse der FAL zeigten ausgeprägtere funktionelle Abweichungen in der schwer betroffenen Patientengruppe (Gruppe 2). Es wurde jedoch kein signifikanter Unterschied zur Gruppe 1 ermittelt. Die schwer betroffene Patientengruppe zeigten einen höheren Prozentsatz an von Schmerzen bei Palpation der Kiefergelenke (Gruppe 2: 70,0 % vs. Gruppe 1: 61,1 %) oder Schmerzen bei Mundöffnung (Gruppe 2: 10,0% vs. Gruppe 1: 8,3 %), Deflexionen des Unterkiefers (Gruppe 2: 50,0% vs. Gruppe 1: 33,3 %), Gelenkgeräuschen (Gruppe 2: 80,0 % vs. Gruppe 1: 63,9 %) und Mundöffnungseinschränkungen (Gruppe 2: 60,0 % vs. Gruppe 1: 25,0 %). Die durchschnittliche Mundöffnung betrug in Gruppe 2 40,6 mm, während sie in Gruppe 1 43,5 mm betrug. Bei Patienten mit einer Mundöffnung < 40 mm wurde in Gruppe 2 eine durchschnittliche Mundöffnung von 35,3 mm und in Gruppe 1 von 34,1 mm gemessen.
JVA:
Sowohl in der gelenkbezogenen, als auch in der patientengruppenbezogenen Analyse der JVA deuteten die Messparameter in der schwer betroffenen Patientengruppe vermehrt auf chronisch-degenerative oder bestehende Ergüsse im Kiefergelenk hin. In der gelenkbezogenen Auswertung zeigte sich dies insbesondere durch eine reduzierte Signalstärke in der schwer betroffenen Patientengruppe (Total power: p = 0,005; Power < 300 Hz: p = 0,006; Power > 300 Hz: 0,003;) sowie in einer signifikant erhöhten Peak Frequenz (p = 0,036).
OPG:
In der Auswertung der OPGs war die Ratio von Kondylus-, Ramus- und Mandibulahöhe in der schwer Patientengruppe (Ratio 79,6 %, 85,9 %, 86,5 %) signifikant kleiner (Kondylushöhe: p = 0,0005; Ramushöhe: p = 0,0030; Mandibulahöhe: p = 0,0002), als in der leicht betroffenen Patientengruppe (Ratio 93,8 %, Ratio 96,0 %, 95,6 %). Somit ergaben sich in der schwer betroffenen Patientengruppe signifikant stärker ausgeprägte Unterkieferasymmetrien, als in der leicht betroffenen Patientengruppe.
3d-Scan:
Im 3d-Scan kam es bei Patienten mit schwer betroffenen Kiefergelenken signifikant häufiger zu Abweichungen des Weichteilkinns von der Medianebene (Gruppe 2: 3,0 mm vs. Gruppe 1: 1,2 mm; p = 0,041) und zu Asymmetrien des Unterkiefers (Gruppe 2: 62,5 % vs. Gruppe 1: 14,8 %; p = 0,015) als bei Patienten mit leicht betroffenem Kiefergelenk.
Schlussfolgerung:
Es zeigte sich, dass mittels einfacher und schnell verfügbarer Untersuchungsmethoden wie der klinischen Funktionsanalyse, der Joint Vibration Analysis und der OPG-Aufnahme eine Pathologie im Kiefergelenk dargestellt werden kann. Die Methoden können als wichtige Referenz zur Kontrolle der Krankheitsprogression bei Patienten mit JIA dienen. Darüber hinaus ist eine Klassifikation der Kondylen in schwer und leicht betroffene Kiefergelenke mittels pathomorphologischer Analyse möglich. Dabei ist von einem direkten Zusammenhang zwischen Grad der Destruktion, Ausmaß der Unterkieferasymmetrie und Dauer der Erkrankung bei Patienten mit JIA auszugehen. Insgesamt konnte die Wertigkeit klinischer, instrumenteller und bildgebender Untersuchungsverfahren der Kiefergelenksdiagnostik bei Patienten mit juveniler idiopathischer Arthritis dargestellt werden.
Herz-Kreislauf-Erkrankungen stellen weiterhin die Todesursache Nummer eins in Deutschland dar, welche hauptsächlich durch Atherosklerose verursacht werden. Als ein Prädiktor der atherosklerotischen Plaqueentwicklung wird die Wandschubspannung (WSS) diskutiert. Ziel dieser Arbeit war es daher, anhand von Flussphantomen und C57bl/6 Mäusen eine 2D-Gradientenecho-Bildgebungsmethode mit einer 3D-Phasenkontrast-Flusskodierung zu optimieren und anschließend eine longitudinale Kleintierstudie mit ApoE-/- Mäusen durchzuführen, um den Zusammenhang der WSS und der atherosklerotischen Plaqueentwicklung näher zu untersuchen. Zunächst wurden Flussphantome mit einem Schlauchdurchmesser von 4 mm und 1 mm zur Optimierung der Messmethode verwendet. Anschließend wurde die Messmethode weiter angepasst, um in vivo Messungen an C57bl/6 Mäusen durchführen zu können. Nach erfolgter Optimierung wurde eine longitudinale Kleintierstudie mit zwei verschiedenen Diäten, Western Diät und Chow Diät, durchgeführt. Im Rahmen der Studie erfolgten nach einer, acht und zwölf Wochen MR-Messungen sowie histologische Analysen. Es konnte gezeigt werden, dass die Wandschubspannung in Mäusen bei 17,6 Tesla quantifiziert werden kann. Es zeigte sich eine Tendenz, dass Plaqueformationen mit einer höheren Wandschubspannung einhergehen.
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.
In Magnetic Resonance Imaging (MRI), acquisition of dynamic data may be highly complex due to rapid changes occurred in the object to be imaged. For clinical diagnostic, dynamic MR images require both high spatial and temporal resolution. The speed in the acquisition is a crucial factor to capture optimally dynamics of the objects to obtain accurate diagnosis. In the 90’s, partially parallel MRI (pMRI) has been introduced to shorten scan times reducing the amount of acquired data. These approaches use multi-receiver coil arrays to acquire independently and simultaneously the data.
Reduction in the amount of acquired data results in images with aliasing artifacts. Dedicated methods as such Sensitivity Encoding (SENSE) and Generalized Autocalibrating Partially Parallel Acquisition (GRAPPA) were the basis of a series of algorithms in pMRI.
Nevertheless, pMRI methods require extra spatial or temporal information in order to optimally reconstruct the data. This information is typically obtained by an extra scan or embedded in the accelerated acquisition applying a variable density acquisition scheme.
In this work, we were able to reduce or totally eliminate the acquisition of the training data for kt-SENSE and kt-PCA algorithms obtaining accurate reconstructions with high temporal fidelity.
For dynamic data acquired in an interleaved fashion, the temporal average of accelerated data can generate an artifact-free image used to estimate the coil sensitivity maps avoiding the need of extra acquisitions. However, this temporal average contains errors from aliased components, which may lead to signal nulls along the spectra of reconstructions when methods like kt-SENSE are applied. The use of a GRAPPA filter applied to the temporal average reduces these errors and subsequently may reduce the null components in the reconstructed data. In this thesis the effect of using temporal averages from radial data was investigated. Non-periodic artifacts performed by undersampling radial data allow a more accurate estimation of the true temporal average and thereby avoiding undesirable temporal filtering in the reconstructed images. kt-SENSE exploits not only spatial coil sensitivity variations but also makes use of spatio-temporal correlations in order to separate the aliased signals. Spatio-temporal correlations in kt-SENSE are learnt using a training data set, which consists of several central k-space lines acquired in a separate scan. The scan of these extra lines results in longer acquisition times even for low resolution images. It was demonstrate that limited spatial resolution of training data set may lead to temporal filtering effects (or temporal blurring) in the reconstructed data.
In this thesis, the auto-calibration for kt-SENSE was proposed and its feasibility was tested in order to completely eliminate the acquisition of training data. The application of a prior TSENSE reconstruction produces the training data set for the kt-SENSE algorithm. These training data have full spatial resolution. Furthermore, it was demonstrated that the proposed auto-calibrating method reduces significantly temporal filtering in the reconstructed images compared to conventional kt-SENSE reconstructions employing low resolution training images. However, the performance of auto-calibrating kt-SENSE is affected by the Signal-to-Noise Ratio (SNR) of the first pass reconstructions that propagates to the final reconstructions.
Another dedicated method used in dynamic MRI applications is kt-PCA, that was first proposed for the reconstruction of MR cardiac data. In this thesis, kt-PCA was employed for the generation of spatially resolved M0, T1 and T2 maps from a single accelerated IRTrueFISP or IR-Snapshot FLASH measurement. In contrast to cardiac dynamic data, MR relaxometry experiments exhibit signal at all temporal frequencies, which makes their reconstruction more challenging. However, since relaxometry measurements can be represented by only few parameters, the use of few principal components (PC) in the kt-PCA algorithm can significantly simplify the reconstruction. Furthermore, it was found that due to high redundancy in relaxometry data, PCA can efficiently extract the required information from just a single line of training data.
It has been demonstrated in this thesis that auto-calibrating kt-SENSE is able to obtain high temporal fidelity dynamic cardiac reconstructions from moderate accelerated data avoiding the extra acquisition of training data. Additionally, kt-PCA has been proved to be a suitable method for the reconstruction of highly accelerated MR relaxometry data.
Furthermore, a single central training line is necessary to obtain accurate reconstructions. Both reconstruction methods are promising for the optimization of training data acquisition and seem to be feasible for several clinical applications.
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.
Messungen mit Magnetresonanztomographen sind seit jeher mit hohen Lautstärken verbunden. Deshalb wird das Gerät im Volksmund auch als „laute Röhre“ bezeichnet. Bisher wurde das Problem mit Kopfhörern, Ohrenstöpseln und akustischer Dämmung des MRT-Scanners angegangen. Auch in der Fachliteratur wird das Problem als gegeben angesehen und es werden kaum wissenschaftliche Lösungsansätze zur Lautstärkereduktion beschrieben.
Das Ziel der vorliegenden Arbeit war es, Bildgebungs-Sequenzen für schwer‑optimierbare Bildkontraste und sogenannte Standard-Kontraste aus dem klinischen Umfeld hinsichtlich der Lautstärke zu optimieren. Viele dieser Kontraste können bereits mit einfachen Algorithmen wie dem Gradientenglättungsalgorithmus erfolgreich in Hinblick auf die Lautstärke optimiert werden. Allerdings existieren auch Sequenzen beziehungsweise Kontraste, die aufgrund ihrer Eigenschaften nicht von einem solchen Algorithmus profitieren können. Die Optimierungen und Änderungen sollten software-seitig erfolgen, das heißt durch Änderung der Gradientenformen und Datenakquisition.
In der Arbeit wurden die grundlegenden Zusammenhänge zwischen den verwendeten Geräteparametern und der Lautstärke untersucht und zudem die physikalischen Ursachen der Lautstärkeentwicklung hergeleitet. Diese konnten anhand der Lorentz-Kräfte quantitativ beschrieben werden. Somit konnten die Hauptursachen der Lautstärkeentwicklung identifiziert werden. Diese sind abhängig von der Gradienten-Steig-Rate, aber auch von der Amplitude der Gradienten. Es konnte gezeigt werden, dass eine Minimierung dieser Gradientenparameter zu einer geringeren Lautstärkeentwicklung führt. Allerdings führt diese Minimierung in den meisten Fällen auch zu einer systematischen Verlangsamung des Sequenzablaufs, was das Erreichen bestimmter Echozeiten und Bildkontraste unmöglich macht.
Zu den problematischen Kontrasten bezüglich der Lautstärkereduktion zählten der T1- und PD‑Kontrast einer Turbo-Spin-Echo-Sequenz. Durch die Kombination von mehreren Maßnahmen, wie der Adaption der k-Raum-Akquisition, der HF-Pulse-Parameter und den Gradientenformen, war es möglich, die Lautstärke in Beispielmessungen um bis zu 16,8 dB(A) zu reduzieren. Wie bei der kürzlich veröffentlichten Methode zur Reduktion für die T2‑gewichteten Kontraste, wurde dies zulasten einer Messzeitverlängerung von bis zu 50% erreicht. Die Endlautstärke betrug dabei circa 81 dB(A). Mit der Lautstärkeoptimierung der klinisch bedeutsamen T1- und PD‑Kontraste wurde die Palette an leisen, mit der Turbo-Spin-Echo‑Sequenz erzielbaren, Standard-Kontrasten (T1, T2 und PD) nun vervollständigt.
In einem anderen Ansatz wurde die Anwendbarkeit des CAT-Konzepts auf die Lautstärkereduktion untersucht. Beim CAT-Konzept wird die Messung in Einzelmessungen mit verschiedenen Parametern unterteilt. Bisher wurde dieser Ansatz zur SAR-Reduktion verwendet. Das Zentrum des k-Raums wird mit einer SAR-intensiven, kontrastgebenden Messung aufgenommen. Der verbleibende Teil des k-Raums wird mit einer SAR-reduzierten, bildstrukturrelevanten Messung aufgenommen. In dieser Arbeit wurde die Übertragung des CAT-Konzepts auf die Lautstärkereduktion untersucht. Anstelle von SAR-intensiven und SAR‑reduzierten Messungen, wurde hier die Unterteilung in „laute“ und „leise“ Messungen untersucht. Dabei wurden Überlegungen angestellt, die es für eine Vielzahl an Messungen ermöglichen, einen großen Teil der Messung leise zu gestalten ohne die Bildqualität oder den Bildkontrast zu verändern. In einem weiteren Schritt wurden Überlegungen für die Lautstärkereduktion der lauten Messungen vorgestellt. Anschließend wurden für eine GRE- und TSE-Sequenz Optimierungsschritte evaluiert und die Lautstärke gemessen.
Der hinsichtlich der Lautstärkeoptimierung herausforderndste Bildkontrast ist die diffusionsgewichtete Bildgebung. Diese besitzt eine Diffusions-Präparation zur Sichtbarmachung der Diffusivität, bei der die maximal mögliche Gradienten-Amplitude verwendet wird. Ebenso werden nach der Präparation die Daten mit einem EPI‑Akquisitionsmodul mit Blip-Gradienten akquiriert, das mit einem charakteristischem „Pfeifton“ einhergeht. Zum einen wurden die Gradientenformen konsequent angepasst. Zum anderen wurde eine Segmentierung der k-Raum-Akquisition in Auslese-Richtung verwendet, um die Gradienten‑Steig-Raten zu reduzieren. Auch hier konnte eine deutliche Lautstärkereduktion von bis zu 20,0 dB(A) erzielt werden. Dies wurde zulasten einer Messzeitverlängerung von 27% ‑ 34% im Vergleich zur Standard-Sequenz erreicht. Durch eine weitere Messzeitverlängerung um bis zu 23% kann die Lautstärke um weitere 0,9 dB(A) reduziert werden. Dabei hängt die genaue Messzeitverlängerung vom verwendeten GRAPPA-Faktor und der Anzahl der Auslese-Segmente ab. Die entstandene Sequenz wurde in mehreren Kliniken erfolgreich erprobt.
Bisher mussten bei MRT-Messungen stets Kompromisse zwischen „hoher Auflösung“, „hohem SNR“ und „geringer Messzeit“ getroffen werden. Als Anschauung dafür wurde das „Bermuda‑Dreieck der MRT“ eingeführt. Da alle drei Größen sich gegenseitig ausschließen, muss stets ein Mittelweg gefunden werden. Einige der in dieser Arbeit erzielten Erfolge bei der Lautstärkereduktion wurden auf Kosten einer verlängerten Messzeit erreicht. Daher ist es naheliegend, das „Bermuda-Dreieck der MRT“ um die Dimension der „geringen Lautstärke“ zu einer „Bermuda-Pyramide der MRT“ zu erweitern. Damit muss die Lautstärkeentwicklung in die Mittelweg‑Findung miteinbezogen werden.
Die in dieser Arbeit erzielten Lautstärken liegen in der Größenordnung zwischen 80 ‑ 85 dB(A). Somit können Messungen bei Verwendung von Gehörschutz angenehm für den Patienten durchgeführt werden. Durch neue Techniken der Zukunft wird es wahrscheinlich sein, höhere Auflösungen, höheres SNR oder kürzere Aufnahmedauern zu erzielen, beziehungsweise stattdessen diese in eine geringe Lautstärke „umzuwandeln“.
Ebenso werden möglicherweise auf der hardware-technischen Seite Fortschritte erzielt werden, so dass in neueren MRT-Scannergenerationen mehr Wert auf die Lärmdämmung gelegt wird und somit der softwarebasierten Lautstärkereduktion einen Schritt entgegen gekommen wird. Damit könnten zukünftige Patienten-Messungen gänzlich ohne störenden Gehörschutz durchgeführt werden.