Refine
Has Fulltext
- yes (60)
Is part of the Bibliography
- yes (60)
Year of publication
Document Type
- Doctoral Thesis (59)
- Journal article (1)
Keywords
- NMR-Tomographie (60) (remove)
Institute
- Institut für diagnostische und interventionelle Radiologie (Institut für Röntgendiagnostik) (19)
- Physikalisches Institut (19)
- Institut für Anorganische Chemie (5)
- Medizinische Klinik und Poliklinik I (5)
- Graduate School of Life Sciences (3)
- Neurochirurgische Klinik und Poliklinik (2)
- Neurologische Klinik und Poliklinik (2)
- Theodor-Boveri-Institut für Biowissenschaften (2)
- Institut für Hygiene und Mikrobiologie (1)
- Klinik und Poliklinik für Hals-, Nasen- und Ohrenkrankheiten, plastische und ästhetische Operationen (1)
Sonstige beteiligte Institutionen
In der Diagnostik und Therapie der KHK sind das frühzeitige Erkennen und die Beurteilung funktioneller Folgen atherosklerotischer Veränderungen von großer Bedeutung. Die First-Pass MR-Bildgebung ermöglicht Aussagen über die myokardiale Perfusion und damit die hämodynamische Relevanz einer Koronarstenose. In der vorliegenden Arbeit wurden quantitative Werte für die myokardiale Durchblutung gesunder Probanden unter Adenosin-induziertem Stress und in Ruhe unter Einsatz der Präbolustechnik bestimmt. Eine exakte Darstellung der arteriellen Inputfunktion wurde durch einen Kontrastmittelbolus in niedriger Dosierung erreicht, die Verwendung höherer Kontrastmitteldosen führte dagegen zu einem verbesserten Signal-zu-Rausch-Verhältnis im Myokard. Die Absolutwerte der myokardialen Perfusion unter Stressbedingungen und in Ruhe wie auch die myokardiale Perfusionsreserve zeigten vergleichbare Mittelwerte, wiesen aber eine geringere Streubreite im Vergleich zu früheren MR Studien auf und waren vergleichbar mit in PET-Studien erzielten Ergebnissen. Weiterhin wurden unter Verwendung dieser Methode Werte für das myokardiale Verteilungsvolumen des Kontrastmittels als wichtiger Parameter in der Differenzierung von gesundem und infarziertem Herzmuskelgewebe ermittelt und die Laufzeit der Boluspassage nach Injektion in Ruhe und unter Stress bestimmt, die zur Unterscheidung von antegrad perfundiertem und von über Kollateralen versorgtem Myokard dienen kann. Mit Hilfe der MRT war es auch möglich, Unterschiede zwischen subendo- und subepimyokardialer Perfusion zu quantifizieren. Die erzielten Ergebnisse entsprechen bisher publizierten Werten, die mit anderen Modalitäten gewonnen wurden. Der Vergleich der absoluten Perfusion bei verminderter zeitlicher Auflösung mit den bei hoher zeitlicher Auflösung gemessenen Werten ergab nur geringfügige Abweichungen der Ergebnisse voneinander. Dadurch eröffnet sich die Möglichkeit, durch die Zeitersparnis mehrere Schichten abwechselnd bei verschiedenen Herzschlägen zu messen und damit eine erweiterte Abdeckung des linksventrikulären Myokards zu erreichen. Durch die quantitative Auswertung der First-Pass MR-Perfusionsmessung stellt die beschriebene Methode eine vielversprechende Option im Bereich der nichtinvasiven Diagnostik verschiedener myokardialer Erkrankungen dar.
Magnetic Resonance Imaging (MRI) is a non-invasive medical imaging technique, that is rou- tinely used in clinical practice for detection and diagnosis of a wide range of different diseases. In MRI, no ionizing radiation is used, making even repeated application unproblematic. This is an important advantage over other common imaging methods such as X-rays and Computer To- mography. One major drawback of MRI, however, are long acquisition times and associated high costs of experiments. Since the introduction of MRI, several important technical developments have been made to successfully reduce acquisition times. In this work, novel approaches were developed to increase the efficiency of MRI acquisitions.
In Chapter 4, an improved radial turbo spin-echo (TSE) combined acquisition and reconstruction strategy was introduced. Cartesian turbo spin-echo sequences [3] are widely used especially for the detection and diagnosis of neurological pathologies, as they provide high SNR images with both clinically important proton density and T2 contrasts. TSE acquisitions combined with radial sampling are very efficient, since it is possible to obtain a number of ETL images with different contrasts from a single radial TSE measurement [56–58]. Conventionally, images with a particular contrast are obtained from both radial and Cartesian TSE acquisitions by combining data from different echo times into a single image. In the radial case, this can be achieved by employing k-space weighted image contrast (KWIC) reconstruction. In KWIC, the center region of k-space is filled exclusively with data belonging to the desired contrast while outer regions also are assembled with data acquired at other echo times. However, this data sharing leads to mixed contrast contributions to both Cartesian and radial TSE images. This is true especially for proton density weighted images and therefore may reduce their diagnostic value.
In the proposed method, an adapted golden angle reordering scheme is introduced for radial TSE acquisitions, that allows a free choice of the echo train length and provides high flexibility in image reconstruction. Unwanted contrast contaminations are greatly reduced by employing a narrow-band KWIC filter, that restricts data sharing to a small temporal window around the de- sired echo time. This corresponds to using fewer data than required for fully sampled images and consequently leads to images exhibiting aliasing artifacts. In a second step, aliasing-free images are obtained using parallel imaging. In the neurological examples presented, the CG-SENSE algorithm [42] was chosen due to its stable convergence properties and its ability to reconstruct arbitrarily sampled data. In simulations as well as in different in vivo neurological applications, no unwanted contrast contributions could be observed in radial TSE images reconstructed with the proposed method. Since this novel approach is easy to implement on today’s scanners and requires low computational power, it might be valuable for the clinical breakthrough of radial TSE acquisitions.
In Chapter 5, an auto-calibrating method was introduced to correct for stimulated echo contribu- tions to T2 estimates from a mono-exponential fit of multi spin-echo (MSE) data. Quantification of T2 is a useful tool in clinical routine for the detection and diagnosis of diseases as well as for tis- sue characterization. Due to technical imperfections, refocusing flip angles in a MSE acquisition deviate from the ideal value of 180○. This gives rise to significant stimulated echo contributions to the overall signal evolution. Therefore, T2 estimates obtained from MSE acquisitions typically are notably higher than the reference. To obtain accurate T2 estimates from MSE acquisitions, MSE signal amplitudes can be predicted using the extended phase graph (EPG, [23, 24]) algo- rithm. Subsequently, a correction factor can be obtained from the simulated EPG T2 value and applied to the MSE T2 estimates. However, EPG calculations require knowledge about refocus- ing pulse amplitudes, T2 and T1 values and the temporal spacing of subsequent echoes. While the echo spacing is known and, as shown in simulations, an approximate T1 value can be assumed for high ratios of T1/T2 without compromising accuracy of the results, the remaining two parameters are estimated from the data themselves. An estimate for the refocusing flip angle can be obtained from the signal intensity ratio of the second to the first echo using EPG. A conventional mono- exponential fit of the MSE data yields a first estimate for T2. The T2 correction is then obtained iteratively by updating the T2 value used for EPG calculations in each step. For all examples pre- sented, two iterations proved to be sufficient for convergence. In the proposed method, a mean flip angle is extracted across the slice. As shown in simulations, this assumption leads to greatly reduced deviations even for more inhomogeneous slice profiles. The accuracy of corrected T2 values was shown in experiments using a phantom consisting of bottles filled with liquids with a wide range of different T2 values. While T2 MSE estimates were shown to deviate significantly from the spin-echo reference values, this is not the case for corrected T2 values. Furthermore, applicability was demonstrated for in vivo neurological experiments.
In Chapter 6, a new auto-calibrating parallel imaging method called iterative GROG was pre- sented for the reconstruction of non-Cartesian data. A wide range of different non-Cartesian schemes have been proposed for data acquisition in MRI, that present various advantages over conventional Cartesian sampling such as faster acquisitions, improved dynamic imaging and in- trinsic motion correction. However, one drawback of non-Cartesian data is the more complicated reconstruction, which is ever more problematic for non-Cartesian parallel imaging techniques. Iterative GROG uses Calibrationless Parallel Imaging by Structured Low-Rank Matrix Completion (CPI) for data reconstruction. Since CPI requires points on a Cartesian grid, it cannot be used to directly reconstruct non-Cartesian data. Instead, Grappa Operator Gridding (GROG) is employed in a first step to move the non-Cartesian points to the nearest Cartesian grid locations. However, GROG requires a fully sampled center region of k-space for calibration. Combining both methods in an iterative scheme, accurate GROG weights can be obtained even from highly undersampled non-Cartesian data. Subsequently, CPI can be used to reconstruct either full k- space or a calibration area of arbitrary size, which can then be employed for data reconstruction with conventional parallel imaging methods.
In Chapter 7, a new 2D sampling scheme was introduced consisting of multiple oscillating effi- cient trajectories (MOET), that is optimized for Compressed Sensing (CS) reconstructions. For successful CS reconstruction of a particular data set, some requirements have to be met. First, ev- ery data sample has to carry information about the whole object, which is automatically fulfilled for the Fourier sampling employed in MRI. Additionally, the image to be reconstructed has to be sparse in an arbitrary domain, which is true for a number of different applications. Last, data sam- pling has to be performed in an incoherent fashion. For 2D imaging, this important requirement of CS is difficult to achieve with conventional Cartesian and non-Cartesian sampling schemes. Ra- dial sampling is often used for CS reconstructions of dynamic data despite the streaking present in undersampled images. To obtain incoherent aliasing artifacts in undersampled images while at the same time preserving the advantages of radial sampling for dynamic imaging, MOET com- bines radial spokes with oscillating gradients of varying amplitude and alternating orientation orthogonal to the readout direction. The advantage of MOET over radial sampling in CS re- constructions was demonstrated in simulations and in in vivo cardiac imaging. MOET provides superior results especially when used in CS reconstructions with a sparsity constraint directly in image space. Here, accurate results could be obtained even from few MOET projections, while the coherent streaking artifacts present in the case of radial sampling prevent image recovery even for smaller acceleration factors. For CS reconstructions of dynamic data with sparsity constraint in xf-space, the advantage of MOET is smaller since the temporal reordering is responsible for an important part of incoherency. However, as was shown in simulations of a moving phantom and in the reconstruction of ungated cardiac data, the additional spatial incoherency provided by MOET still leads to improved results with higher accuracy and may allow reconstructions with higher acceleration factors.
Magnetic Resonance Imaging (MRI) is an imaging modality which provides anatomical or functional images of the human body with variable contrasts in an arbitrarily positioned slice without the need for ionizing radiation. In MRI, data are not acquired directly, but in the reciprocal image space (otherwise known as k-space) through the application of spatially variable magnetic field gradients. The k-space is made up of a grid of data points which are generally acquired in a line-by-line fashion (Cartesian imaging). After the acquisition, the k-space data are transformed into the image domain using the Fast Fourier Transformation (FFT). However, the acquisition of data is not limited to the rectilinear Cartesian sampling scheme described above. Non-Cartesian acquisitions, where the data are collected along exotic trajectories, such as radial and spiral, have been shown to be beneficial in a number of applications. However, despite their additional properties and potential advantages, working with non-Cartesian data can be complicated. The primary difficulty is that non-Cartesian trajectories are made up of points which do not fall on a Cartesian grid, and a simple and fast FFT algorithm cannot be employed to reconstruct images from non-Cartesian data. In order to create an image, the non-Cartesian data are generally resampled on a Cartesian grid, an operation known as gridding, before the FFT is performed. Another challenge for non-Cartesian imaging is the combination of unusual trajectories with parallel imaging. This thesis has presented several new non-Cartesian parallel imaging methods which simplify both gridding and the reconstruction of images from undersampled data. In Chapter 4, a novel approach which uses the concepts of parallel imaging to grid data sampled along a non-Cartesian trajectory called GRAPPA Operator Gridding (GROG) is described. GROG shifts any acquired k-space data point to its nearest Cartesian location, thereby converting non-Cartesian to Cartesian data. The only requirements for GROG are a multi-channel acquisition and a calibration dataset for the determination of the GROG weights. Chapter 5 discusses an extension of GRAPPA Operator Gridding, namely Self-Calibrating GRAPPA Operator Gridding (SC-GROG). SC-GROG is a method by which non-Cartesian data can be gridded using spatial information from a multi-channel coil array without the need for an additional calibration dataset, as required in standard GROG. Although GROG can be used to grid undersampled datasets, it is important to note that this method uses parallel imaging only for gridding, and not to reconstruct artifact-free images from undersampled data. Chapter 6 introduces a simple, novel method for performing modified Cartesian GRAPPA reconstructions on undersampled non-Cartesian k-space data gridded using GROG to arrive at a non-aliased image. Because the undersampled non-Cartesian data cannot be reconstructed using a single GRAPPA kernel, several Cartesian patterns are selected for the reconstruction. Finally, Chapter 7 discusses a novel method of using GROG to mimic the bunched phase encoding acquisition (BPE) scheme. In MRI, it is generally assumed that an artifact-free image can be reconstructed only from sampled points which fulfill the Nyquist criterion. However, the BPE reconstruction is based on the Generalized Sampling Theorem of Papoulis, which states that a continuous signal can be reconstructed from sampled points as long as the points are on average sampled at the Nyquist frequency. A novel method of generating the “bunched” data using GRAPPA Operator Gridding (GROG), which shifts datapoints by small distances in k-space using the GRAPPA Operator instead of employing zig-zag shaped gradients, is presented in this chapter. With the conjugate gradient reconstruction method, these additional “bunched” points can then be used to reconstruct an artifact-free image from undersampled data. This method is referred to as GROG-facilitated Bunched Phase Encoding, or GROG-BPE.
Zur Beurteilung der Lungenanatomie wurde das MT-STIR-Verfahren vorgestellt. Es wurde gezeigt, dass das MT-STIR-Verfahren das störende Signal des umgebenden Muskelgewebes effektiv unterdrückt und damit die Visualisierung des Lungenparenchyms verbessert. Im Vergleich zu konventionellen anatomischen 1H-MR-Verfahren wie IR- und MIR-Verfahren erhöht das MT-STIR-Verfahren das Signal-zu-Rausch-Verhältnis (SNR) des Lungenparenchyms signifikant und vermeidet den Signalausfall des Lungenparenchyms aufgrund der pathologischen Verkürzung der Lungen-T1-Relaxationszeit auf ca. 900 ms wie bei Patienten mit Mukoviszidose (CF), so dass sowohl große Lungenperfusionsdefekte in Patienten mit CF als auch kleine ungefährliche Lungenentzündungen in „gesunden“ Probanden durch das MT-STIR-Verfahren gut dargestellt werden können. Für die indirekte, aber quantitative Beurteilung der Lungenventilation wurde die oben genannte schnelle quantitative Lungen-T1-Mapping-Technik während der Inhalation eines Atemgasgemisches mit verschiedenen O2-Konzentrationen (21%, 40%, 60%, 80% und 100%) eingesetzt. Dabei ist im Blut physikalisch gelöster Sauerstoff leicht paramagnetisch und dient als Blut-T1-verkürzendes MR-Kontrastmittel (KM). In der Lunge ist das Blut die Hauptquelle des freien Wassers, so dass Lungen-T1-Werte nach dem Zwei-Kompartimente-Schnellaustausch-Modell der Lungen-T1-Relaxationszeit durch den Blut-T1-Wert beeinflusst werden. Die zugehörige Theorie, ein O2-gestütztes Lungen-T1-Modell, wurde aus der Lungenphysiologie und den T1-Relaxationsmechanismen hergeleitet und zeigt, dass bei Probanden die Lungen-T1-Verkürzung von 21% O2 zu 100% O2 ca. 11% beträgt und die Beziehung zwischen dem Lungen-R1 (= 1/T1)-Wert und der inhalierten O2-Konzentration linear mit einer Steigung von 0,12 1/s und einem R1-Achsenabschnitt von 0,70 1/s ist. Die Steigung wurde im Rahmen dieser Doktorarbeit als oxygen transfer function (OTF) definiert und ist vom gasaustauschbestimmenden Ventilations-Perfusions- und Diffusions-Perfusions-Verhältnis abhängig, so dass sie praktisch ein Maß für den pulmonalen Gasaustausch darstellt. Experimentell wurde gezeigt, dass Lungen-T1-Werte bei 100% O2 um 10% kürzer als bei 21% O2 sind, was gut mit dem O2-gestützten Lungen-T1-Modell übereinstimmt. Weiterhin wurde die OTF dadurch bestimmt, dass die gemessenen Lungen-R1-Werte gegen die inhalierte O2-Konzentration aufgetragen wurden und eine Gerade an die Messpunkte angepasst wurde. Gesundes Lungenparenchym von Probanden und gut perfundiertes Lungenparenchym von Patienten mit CF zeigten OTF-Werte zwischen 0,10 und 0,14 1/s, R1-Achsenabschnitte zwischen 0,70 und 0,80 1/s und ausgezeichnete Korrelationskoeffizienten von annähernd 1,00, was mit dem O2-gestützten Lungen-T1-Modell übereinstimmt. Schlecht perfundiertes Lungenparenchym von Patienten mit CF zeigte eindeutig erniedrigte OTF-Werte, erhöhte R1-Achsenabschnitte und schlechte Korrelationskoeffizienten. Das O2-gestützte Lungen-T1-Mapping-Verfahren zeigt eine hohe Reproduzierbarkeit. Zur Beurteilung der Lungenperfusion wurde eine quantitative Perfusionsmapping-Technik mittels Protonen-Spin-Labeling, ohne Verwendung eines intravenösen Kontrastmittels wie Gadolinium (Gd)-DTPA, vorgestellt. Aus einer nicht-schichtselektiven (globalen) T1-Map und einer schichtselektiven T1-Map derselben Lungenschicht, die jeweils mit der oben genannten schnellen quantitativen Lungen-T1-Mapping-Technik akquiriert wurde, wurde eine Perfusionamap berechnet, wobei jedes Pixel in der Perfusionsmap eine Perfusionsrate in Einheiten von m/100g/min hat. Es wurde demonstriert, dass die hintere coronale Lungenschicht eine höhere Perfusionsrate als die vordere coronale Lungenschicht desselben Probanden hatte, als er in Rückenlage gemessen wurde, was den Gravitationseffekt auf die Lungenperfusion bestätigt. Die berechneten Perfusionsraten des gut perfundierten Lungenparenchyms von Probanden und von Patienten mit CF lagen zwischen 400 und 600 m/100g/min, die gut mit dem Literaturwert übereinstimmen. Die berechneten Perfusionsraten des schlecht perfundierten Lungenparenchyms von Patienten mit CF waren niedriger als 200 m/100g/min. Die Spin-Labeling-Technik zeigte eine hohe Reproduzierbarkeit und niedrige relative Fehler der berechneten Perfusionsraten.
This thesis focuses on various aspects and techniques of 19F magnetic resonance (MR). The first chapters provide an overview of the basic physical properties, 19F MR and MR sequences related to this work. Chapter 5 focuses on the application of 19F MR to visualize biological processes in vivo using two different animal models. The dissimilar models underlined the wide applicability of 19F MR in preclinical research. A subsection of Chapter 6 shows the application of compressed sensing (CS) to 19F turbo-spin-echo chemical shift imaging (TSE-CSI), which leads to reduced measurement time. CS, however, can only be successfully applied when a sufficient signal-to-noise ratio (SNR) is available. When the SNR is low, so-called spike artifacts occur with the CS algorithm used in the present work. However, it was shown in an additional subsection that these artifacts can be reduced using a CS-based post processing algorithm. Thus, CS might help overcome limitations with time consuming 19F CSI experiments. Chapter 7 deals with a novel technique to quantify the B+1 profile of an MR coil. It was shown that, using a specific application scheme of off resonant pulses, Bloch-Siegert (BS)-based B+1 mapping can be enabled using a Carr Purcell Meiboom Gill (CPMG)-based TSE sequence. A fast acquisition of the data necessary for B+1 mapping was thus enabled. In the future, the application of BS-CPMG-TSE B+1 mapping to improve quantification using 19F MR could therefore be possible.
Die vorliegende Dissertation befasst sich mit der Entwicklung und Anwendung von Schnellbild-Verfahren der Magnetresonanztomografie (MRT), die auf der Messung stimulierter Echosignale beruhen. Die schnellen STEAM (STEAM = stimulated echo acquisition mode) MRT-Technik leidet im Gegensatz zu alternativen Messtechniken auch bei hohen Magnetfeldstärken nicht unter Bildverzerrungen und Begrenzungen durch eine hohe Hochfrequenzbelastung. Da sie jedoch Bilder geringerer Signalintensität liefert, bestand das Primärziel dieser Arbeit in einer Verbesserung des Signal-zu-Rausch-Verhältnisses, um die Vorteile für die Herzbildgebung und die diffusionsgewichtete MRT des Gehirns besser nutzen zu können (Anwendungen am Menschen bei einer Feldstärke von 3 Tesla). Durch eine Kombination verschiedener physikalischer (Messtechnik, Ortskodierung, Unterabtastung) und mathematischer Maßnahmen (Bildrekonstruktion) konnte in dieser Arbeit eine erhebliche Signalsteigerung bei gleichzeitiger Verkürzung des Messzeit von Einzelbildern erreicht werden. Die Ergebnisse der STEAM-MRT am Herzen zeigen Schnittbilder des Herzmuskels ohne den störenden Einfluss der Signale des Blutes. Die diffusionsgewichtete MRT des Gehirns lieferte ohne Suszeptibilitätsartefakte eine verzerrungsfreie Kartierung des Diffusionstensors und – daraus abgeleitet – eine anatomisch korrekte, dreidimensionale Rekonstruktion von Nervenfaserbahnen etwa des Cingulums.
In der vorliegenden Arbeit wurde die myokardiale Perfusion gesunder Probanden in Ruhe und unter Adenosin–induziertem, sowie Kälte-induziertem Stress quantitativ mittels First-Pass-Perfusions MR-Bildgebung beim gleichen Probanden untersucht und die Perfusionswerte alters- und geschlechtsabhängig verglichen. Hierbei wurden zwanzig Probanden in altersabhängige Gruppen eingeteilt, wobei der cut-off vierzig Jahre war. Diese Arbeit zeigt, dass der Einsatz des CPT zur Beurteilung der endothelvermittelten Vasodilatation im MRT möglich ist. Im Vergleich zum Einsatz von Adenosin zeigt die Stress-Untersuchung mittels CPT weniger bis keine unerwünschten Nebenwirkungen und die Kontraindikationen sind geringer. Somit ist der Einsatz des CPT zur Beurteilung der myokardialen Stressperfusion als eine nebenwirkungsarme Untersuchungsmodalität anzusehen, die es ermöglicht, die frühesten Veränderungen des Endothels bei beginnender koronarer Arteriosklerose zu detektieren. Als Ergebnis ist festzuhalten, dass in dieser Arbeit alle Probanden eine statistisch signifikante Perfusionserhöhung beim CPT aufweisen. Es zeigt sich kein statistisch signifikanter Alters- oder Geschlechtsunterschied beim CPT. Auf die Verabreichung von Adenosin reagiert die junge Probandengruppe mit einer statistisch signifikant höheren Myokardperfusion. Dies könnte in erster Linie durch den höheren Anteil an Frauen in der jungen Gruppe zu erklären sein.
In der vorliegenden prospektiven experimentellen Studie wurden jeweils 18 Probanden sowie Patienten, die unter einer genetisch gesicherten Friedrich-Ataxie leiden, mit der Tissue-Phase‐Mapping (TPM) Sequenz im MRT untersucht. Mit der erwähnten Sequenz ist es möglich, die Geschwindigkeit der Herzwandbewegung über einen Herzzyklus zeitlich hoch aufgelöst (13,8 ms) darzustellen. Es wurde in der Vergangenheit gezeigt, dass die Daten der TPM‐Messung reproduzierbar sowie mit denen aus einer Ultraschalluntersuchung gewonnenen Daten vergleichbar sind. Die Aufnahme erfolgt unter freier Atmung in Navigatortechnik. Dadurch ist diese Untersuchung auch bei Patienten möglich, die aufgrund Ihrer Erkrankung sonst nicht ausreichend lange die Luft anhalten könnten. Die Friedreich‐Ataxie ist die häufigste aller Ataxien in der Adoleszenzphase und wird autosomal-rezessiv vererbt. Neben neurologischen Ausfällen kann es auch zu einer kardialen Beteiligung kommen. Dabei zeigt sich bei einem großen Teil der Patienten eine Kardiomyopathie mit asymmetrischer septaler Hypertrophie sowie einer dynamischen linksventrikulären Ausflussobstruktion. Die American Heart Association klassifiziert dieses Krankeitsbild als sekundäre Kardiomyopathie neuromuskulären Ursprungs. Es wurde in der Vergangenheit bereits gezeigt, dass es unter einer Therapie mit dem Medikament Idebenone zu einer Verbesserung der kardialen Funktion kommen kann. Im Zusammenhang mit einer groß angelegten, multizentrischen Phase III Studie (MICONOS-Studie) kamen diese Patienten nach Würzburg für eine MRT-Untersuchung. Im Rahmen dessen erfolgte dann auf freiwilliger Basis eine Untersuchung mit der TPM-Sequenz. Diese Substudie wurde von der lokalen Ethikkommission genehmigt und die Patienten wie auch Probanden haben selbst oder durch Erziehungsberechtigte der Untersuchung zugestimmt. Es konnte durch diese Arbeit gezeigt werden, dass zum einen die Methode des Tissue-Phase‐Mappings erfolgreich am Institut für Röntgendiagnostik des Universitätsklinikums Würzburg eingeführt werden konnte. Zum anderen, dass die Ergebnisse der herzgesunden Probanden mit anderen, in der Vergangenheit durchgeführten Studien, ähnlich und vergleichbar sind. Zudem konnten teilweise signifikante Unterschiede der systolischen wie auch der diastolischen Geschwindigkeiten zwischen Probanden und Patienten der MICONOS Studie in einigen ROIs bzw. global nachgewiesen werden. Dies bedeutet, dass das hier vorgestellte Verfahren in der Lage ist, Unterschiede in den Geschwindigkeiten der regionalen wie auch globalen Herzwandbewegung in allen drei Bewegungsachsen zwischen verschiedenen Kollektiven zu detektieren. Diese Differenzen können ein Frühzeichen für eine pathologische Herzmuskelerkrankung sein und somit helfen, dass frühzeitig eine entsprechende Therapie begonnen wird. Aufgrund der Ergebnisse dieser Dissertation werden weitere Studien folgen, die helfen werden, dass das Tissue-Phase‐Mapping Verfahren Einzug in die klinische Routine erhalten wird.
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.