Spelling suggestions: "subject:"kernspintomografie"" "subject:"kernspintomographie""
1 |
Wertigkeit der Absoluten Quantitativen Perfusionsauswertung nach Myokardinfarkt im Akutstadium und im Langzeitverlauf mittels Magnetresonanztomographie / Validity of perfusion measurement using MRI in the acute stage of myocardial infarction and upon one year follow upKunz, Julian January 2021 (has links) (PDF)
Es wurden Perfusionsmessungen mittels MRT an Infarktpatienten im Akutstadium und im Langzeitverlauf durchgeführt und quantitativ mittels einem Sektormodell ausgewertet. Hierbei zeigte sich, dass sich die Perfusionswerte im Infarktareal und gesunden Myokard nicht signifikant unterschieden und dass sich diese auch im Jahresverlauf nicht signifikant änderten. Es ergab sich auch keine signifikante Korrelation zwischen der Größe des Infarkareales und den gemessenen Perfusionswerten. / Quantitative cardiac perfusion was measured in the acute stage of myocardial infarction by means of MRI, as well as upon one year follow-up. Using a sector based model, measured perfusion values did not differ significantly between infarcted and non-infarcted tissue. Also, a significant change between the perfusion value upon one year follow-up could not be shown. Furthermore, no significant correlation between perfusion values and the size of the infarcted area could be shown.
|
2 |
Morphologische und funktionelle MRT-Infarktcharakterisierung und Entwicklung einer diffusionsgewichteten MRT-Methode / Morphological and functional MRI infarct characterization and development of a diffusion-weighted MRI methodWeber, Daniel January 2011 (has links) (PDF)
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. / Diffusion tensor imaging for the characterization of myocardial infarction in comparison to other methods The aim of this part of this work was to evaluate the potential of different MR sequences for the characterization of myocardial infarction in both the acute and chronic case. Therefore a study of acute as well as chronic infarcted rat hearts was performed, and the parameters T1, T2, T2* and the parameters ADC, FA, cs, cp and cl calculated from the diffusion tensor images were investigated. It turned out that there is no equivalent to the ischemia. Neither in the acute nor in the chronic case, a notably difference inside the affected area was detectable between different sequences. All in all, for detailed characterization of the infarct scar a T2* or a diffusion tensor sequence are most suitable. The T2* sequence provides a more informative visual image, whereas the more time-consuming diffusion tensor sequence provides a surplus of information due to the multiple display options in post-processing and shows the remodelling of the scar tissue over time. Oxygen level measurements in mouse hearts in vivo The characterization of an infarct scar can also go beyond the representation of morphological structure. The oxygenation is a complex parameter that can provide functional information of the vascularization and viability of the tissue. Access to this parameter is obtained by T2*-measurements, as the parameter T2* is sensitive to chemically bound oxygen. The influence of pure oxygen breathing in contrast to normal room air breathing on the oxygenation level in healthy and infarcted mice have been explored. Despite the difficulties caused by the movement due to respiration and heartbeat the measurements were implemented and carried out at 17.6 Tesla in vivo. The resolution was sufficient to resolve and investigate extremely thinned heart walls after infarction. The effect on the oxygenation level varies considerably between normal and infarcted hearts; that may be caused by a not yet advanced revascularization of the scar. In addition, the method was implemented to a 7.0 Tesla magnet for use in rats and adapted to the respiration of rats, which is different to the respiration of mice. The first reason was that the influence of the higher magnetic field on the measurement of the oxygenation level could be examined. Second, the heart as the examined object is larger in rats. Diffusion weighting using hole burning The MR diffusion method used in this work for the characterization of myocardial infarctions could be limited by extremely short T2 relaxation times. With the strong magnetic fields used here the measurable signal decays very fast due to the relaxation time T2. Therefore, a method for a completely new approach to diffusion-weighted imaging was developed, where the diffusion weighting can be obainted without being limited by the time constant T2. The so-called hole-burning diffusion sequence uses only the longitudinal magnetization for the diffusion weighting in a preliminary experiment. The signal is then acquired with a fast read-out sequence. During the preparation stripes will be "burned" into the magnetization on a subvoxel level, i.e. the magnetization is saturated there. Until the next saturation pulse the behavior of the magnetization depends first on the T1 relaxation time in this area and second on the diffusion. By rapidly repeating the selective pulse train a steady state magnetization dependend on the diffusion constant D and the T1 relaxation time is reached. In this work the dependencies between different sequence parameters were investigated and optimized using simulations. In addition, the sequence was implemented on a MR scanner and first experiments were carried out. With simulated lookup-tables we were able to quantify both the diffusion coefficient D and the T1 relaxation time in the case of not too short relaxation times T1.
|
3 |
Compressed Sensing in der funktionellen kardialen Magnetresonanztomographie / Compressed sensing in functional cardiac magnetic resonance imagingWech, Tobias January 2012 (has links) (PDF)
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. / Because of its high reproducibility and its low variability, magnetic resonance imaging (MRI) of the heart is considered the gold-standard for assessing the cardiac function. In preclinical research, magnetic resonance imaging equally provides an accurate characterization of the cardiac function and enables an excellent analysis of modeled diseases. However, there is still a need for improvement in both applications. Clinical cardiac MRI represents a sophisticated procedure featuring long scan times. This renders the examination comparatively expensive. In preclinical imaging of small animals, long scan times have to be accepted to obtain the required high spatial and temporal resolution. Fast imaging is thus essential for an effective application of cardiac MRI in large collectives. Besides the improvement of the scanner hardware and the optimization of imaging sequences, research in the last decade concentrated on procedures to accelerate MR-data acquisition by exploiting information theory. While numerous publications were associated with parallel imaging (PI) at the beginning of this millennium, the compressed sensing theory (CS) recently gained more and more interest. The latter technique enables the reconstruction of signals from undersampled linear measurements (e.g. the Fourier basis) by exploiting the sparsity of the signal in any known transform domain. As MRI is perfectly qualified for an application of CS, a lot of publications already report on dedicated research. However, the algorithms needed for the reconstruction of undersampled data significantly alter the imaging process of MRI. Classical MRI could be assumed to be linear and stationary in a sufficiently good approximation. The introduction of CS into MRI means a change towards a non-linear and non-stationary transformation. Object information is no longer transferred into an image independently from its location and proportional to its intensity. The image is rather the result of an optimization process maximizing both the fidelity to measured data as well as the sparsity of the signal. The first chapter of this thesis describes a method to objectively evaluate the image quality of MR images reconstructed by CS algorithms. The acceleration with CS promises a reduction of scan time while preserving the image quality. The latter, however, has only been assessed qualitatively or in an insufficient quantitative manner. While classical (linear and stationary) MRI could be validated robustly and effectively by determining a point spread function (PSF), CS algorithms prohibit a corresponding analysis in an analogous manner due to their non-linear and non-stationary behavior. Therefore, the application of a local point spread function (LPSF) was investigated for the iterative soft thresholding (IST) algorithm used in this thesis, to enable a comparative evaluation for imaging systems including CS. The LPSF considers the local dependency of the CS algorithm and thus has to be determined in every location (pixel) of an image. In addition, the LPSF was defined in the linear part of the CS transformation. Small local perturbation on the image to be evaluated were reconstructed for this purpose. The width of the main lobe of the LPSF was used to perform spatially resolved studies on the resolution. The influence of typical undersampling schemes for CS as well as the usage of a discrete gradient transform for a further sparsification were investigated. Subsequently, the procedure was used to assess the spatial and temporal resolution in cardiac MRI. For all CS reconstructions performed in this work, the method allowed a solid and objective analysis of the image resolution. While up to now, comparisons to a fully sampled reference are widely used for a quality assessment, the proposed resolution evaluation represents a step towards a standardized analysis of images obtained by exploiting CS acceleration. The study on sampling schemes revealed, that also for CS accelerated acquisitions, the highest frequencies of the desired k-space have to be included. Former publications proposed undersampling patterns which partly featured a strong weight towards the center of k-space. The results of this thesis put these findings into perspective, as a loss in resolution has been observed for according approaches, at least for the simulations performed in this work. The dynamic acquisitions which were reconstructed exploiting x-f-sparsity proved to be prone to temporal blurring. This is substantiated by the suppression of high temporal frequencies and was analyzed by means of spatially resolved maps. Besides the resolution, an investigation of potential aliasing artifacts as well as the signal-to-noise-ratio (SNR) is essential for a comprehensive quality evaluation. While aliasing may also be investigated by means of the entries of the LPSF outside the main lobe, a modification of the multi-replica method proposed by Robson et al. was presented in chapter 5 to analyze the noise in CS reconstructed images. Taking into account all quality parameters, a robust evaluation of image quality is possible, even when CS is included in the imaging process. This allows a more objective comparison between new developments and present standard procedures and thus may aid the introduction of CS in clinical imaging. After the theoretical analysis on image quality, the next part of this thesis reports on the first application of CS to accelerate functional cardiac MRI of small animals. The studies were performed in cooperation with the British Heart Foundation Experimental Magnetic Resonance Unit (BMRU) of the University of Oxford. The algorithms needed for the CS acceleration were developed and optimized by means of the data acquired by the BMRU at their 9,4 T scanner.\\ An acceleration solely based on CS was investigated first. For this purpose, an ECG- and respiratory gated Cartesian cine acquisition was undersampled in phase encoding direction and reconstructed using CS. The dynamic time series was sparsified by determining temporal difference images for every time frame. A fully sampled temporal average image was obtained by varying the sampling pattern in the temporal dimension. Subsequently, this average image was subtracted from the images of individual heart phases, yielding the sparse temporal difference images. In the validation stage of the study, fully sampled cine acquisitions of mouse hearts were retrospectively undersampled in order to figure out the maximum possible CS acceleration. Cardiac functional parameters were determined for each group of a certain undersampling factor and compared by a statistical analysis. It was shown that a three-fold acceleration is possible without any degradation in the accuracy of the method. This undersampling factor was then validated in an accelerated measurement with a subsequent CS reconstruction. For this purpose, the sampling patterns were further optimized using the transform point spread function. In the subsequent chapter, the CS theory was combined with PI to further increase the acceleration. Again, the phase encoding direction of a Cartesian trajectory was undersampled. The acquisitions were performed using a 9,4 T scanner equipped with an 8 channel mouse coil. In order to exploit the potential of both techniques, CS and PI were combined in a serial manner. First, the k-space was equidistantly undersampled to enable the application of PI. An additional undersampling according to pseudo random numbers was then performed on the resulting sub-grid to allow an acceleration by CS. In consequence, the reconstruction was performed in a serial manner, too. CS was first applied to reconstruct the equidistantly undersampled sub-grid. GRAPPA was used subsequently to compute the still missing data. The equidistantly undersampled sub-grid was shifted from heart phase to heart phase in order to obtain a fully sampled low temporal resolution k-space for a calibration of the GRAPPA-weights. This procedure spares the acquisition time of a separate calibration scan. The following combinations were investigated by retrospectively undersampling a fully sampled cine dataset: R_CS x R_PI = 2 x 2, 2 x 3, 3 x 2 and 3 x 3. The analysis of the noise behavior, the systematic error and the resolution leads to the conclusion that a six-fold acceleration is possible using the proposed hybrid technique. While an increasing factor of the CS acceleration resulted in a slightly larger systematic error, a higher PI acceleration factor led to a slight noise enhancement. However, noise was suppressed for increasing CS acceleration at the same time. In summary, the deviations were at a level which allowed accelerations of up to R_CS x R_PI = 3 x 3. The determination of LPSFs showed a loss in spatial resolution of approximately 50% for a six-fold and up to 64% for nine-fold acceleration. Obviously, the observed suppression of noise was paid by a reduced image resolution for the comparatively noisy acquisitions in small animals. The increased blurring at the endocardial border impedes the segmentation and represents a possible source of error. Taking into consideration all results, a six-fold acceleration (R_CS x R_PI = 2 x 3, 3 x 2) seems reasonable. The additional usage of PI thus enables a further acceleration by a factor of 2 in comparison to an exclusive application of CS. In summary, CS enables a distinct reduction of scan time in preclinical functional cardiac MRI of small animals. Even if no phased-array-coils are available, the necessary amount of data can be reduced to one third without impairing the accuracy of left-ventricular volumes and mass measurements. For acquisitions performed with phased-array-coils, the three-fold acceleration by PI can be extended by an additional two-fold CS acceleration to a joint factor of R=6. Therefore, CS may contribute to an effective application of cardiac MRI in small animals for large collectives. In the last part of the thesis, a modality for clinical functional MRI of the heart was developed. CS was used to enable the acquisition of the whole heart in a single breath-hold of the patient. The current method of choice usually acquires 10-15 2D-slices of the heart, while each measurement requires a separate breath-hold. An undersampled 3D-trajectory was used to reach the necessary acceleration. Phase-encoding in one direction and radial projections in the two remaining ones allowed for an effective acquisition below the Nyquist-criterion. The sparsification of the image series was achieved by subtracting a temporal average image as performed for the preclinical studies. The functionality of the reconstruction technique at an acceleration factor of R ~ 10 was validated in a simulation based on a retrospectively undersampled dataset. The differences between the CS reconstructed and the fully sampled dataset were negligible and thus, the proposed trajectory was implemented at the scanner. An image series depicting the cardiac function with coverage of the full heart was acquired in a single breath-hold of a healthy volunteer using this sequence. Data not covered by the trajectory were reconstructed by the algorithm developed in the validation stage. Due to the slice profile of the fast excitation pulses, additional aliasing artifacts were present in the outer partitions with respect to the images obtained in the simulation. Streaking artifacts of a low intensity were still visible. Apart from that, the dynamics of the heart were excellently captured. The high contrast between the blood pool and the myocardium perfectly qualifies the images for the assessment of cardiac functional parameters. Therefore, the method allows for a higher patient comfort and throughput compared to the gold standard by drastically reducing the amount of necessary breath-holds to a single one.
|
4 |
Relaxations-Dispersions-Bildgebung in der Magnetresonanztomographie / Relaxation Dispersion Magnetic Resonance ImagingHölscher, Uvo Christoph January 2012 (has links) (PDF)
Das Ziel dieser Promotion ist der Aufbau eines dreMR Setups für einen klinischen 1,5T Scanner, das die Relaxations-Dispersions-Bildgebung ermöglicht, und die anschließende Ergründung von möglichst vielen Anwendungsfeldern von dreMR. Zu der Aufgabe gehört die Bereitstellung der zugrunde liegenden Theorie, der Bau des experimentellen Setups (Offset-Spule und Stromversorgung) sowie die Programmierung der nötigen Software. Mit dem gebauten Setup konnten zwei große Anwendungsfelder — dreMR Messungen mit und ohne Kontrastmitteln — untersucht werden. / The goal of this dissertation is the design of a dreMR setup for a clinical 1.5T whole body scanner and the subsequent exploration of possible application fields for the dreMR method. This task includes the investigation of the underlying theory, the design and construction of the dreMR setup (offset-coil and current driver) and the preparation of required software. Two major application fields have been demonstrated: dreMR with and without contrast agents.
|
5 |
Real-time MRI of Moving Spins Using Undersampled Radial FLASH / Echtzeit MRI von bewegten Spins mithilfe der unterabgetasteten radialen FLASH sequenzJoseph, Arun Antony January 2013 (has links) (PDF)
Nuclear spins in motion is an intrinsic component of any dynamic process when studied using magnetic resonance imaging (MRI). Moving spins define many functional characteristics of the human body such as diffusion, perfusion and blood flow. Quantitative MRI of moving spins
can provide valuable information about the human physiology or of a technical system. In particular, phase-contrast MRI, which is based on two images with and without a flow-encoding gradient, has emerged as an important diagnostic tool in medicine to quantify human blood flow. Unfortunately, however, its clinical usage is hampered by long acquisition times which only provide mean data averaged across multiple cardiac cycles and therefore preclude Monitoring the immediate physiological responses to stress or exercise. These limitations are expected to be overcome by real-time imaging which constitutes a primary aim of this thesis.
Short image acquisition times, as the core for real-time phase-contrast MRI, can be mainly realized through undersampling of the acquired data. Therefore the development focused on related technical aspects such as pulse sequence design, k-space encoding schemes and image
reconstruction. A radial encoding scheme was experimentally found to be robust to motion and
less sensitive to undersampling than Cartesian encoding. Radial encoding was combined with a FLASH acquisition technique for building an efficient real-time phase-contrast MRI sequence.
The sequence was further optimized through overlapping of gradients to achieve the shortest possible echo time. Regularized nonlinear inverse reconstruction (NLINV), a technique which jointly estimates the image content and its corresponding coil sensitivities, was used for image
reconstruction. NLINV was adapted specifically for phase-contrast MRI to produce both Magnitude images and phase-contrast maps. Real-time phase-contrast MRI therefore combined two highly undersampled (up to a factor of 30) radial gradient-echo acquisitions with and without a
flow-encoding gradient with modified NLINV reconstructions. The developed method achieved
real-time phase-contrast MRI at both high spatial (1.3 mm) and temporal resolution (40 ms).
Applications to healthy human subjects as well as preliminary studies of patients demonstrated
real-time phase-contrast MRI to offer improved patient compliance (e.g., free breathing) and immediate access to physiological variations of flow parameters (e.g., response to enhanced intrathoracic pressure). In most cases, quantitative blood flow was measured in the ascending aorta as an important blood vessel of the cardiovascular circulation system commonly studied
in the clinic. The performance of real-time phase-contrast MRI was validated in comparison to standard Cine phase-contrast MRI using studies of flow phantoms as well as under in vivo conditions. The evaluations confirmed good agreement for comparable results.
As a further extension to real-time phase-contrast MRI, this thesis implemented and explored a dual-echo phase-contrast MRI method which employs two sequential gradient echoes with and without flow encoding. The introduction of a flow-encoding gradient in between the two echoes
aids in the further reduction of acquisition time. Although this technique was efficient under in vitro conditions, in vivo studies showed the influence of additional motion-induced Phase contributions. Due to these additional temporal phase information, the approach showed Little promise for quantitative flow MRI.
As a further method three-dimensional real-time phase-contrast MRI was developed in this thesis to visualize and quantify multi-directional flow at about twice the measuring time of the standard real-time MRI method, i.e. at about 100 ms temporal resolution. This was achieved
through velocity mapping along all three physical gradient directions. Although the method is still too slow to adequately cover cardiovascular blood flow, the preliminary results were found to be promising for future applications in tissues and organ systems outside the heart. Finally, future developments are expected to benefit from the adaptation of model-based reconstruction
techniques to real-time phase-contrast MRI. / Die Bewegung der Kernspins ist eine wesentliche Eigenschaft von dynamischen Vorgängen, die mit Hilfe der Magnetresonanztomographie (MRT) untersucht werden. Bewegte oder fließende Spins charakterisieren viele Funktionen des menschlichen Körpers, wie z.B. die Gewebeperfusion und den Blutfluss in den Gefäßen. Die quantitative MRT von bewegten Spins kann daher
wertvolle Informationen über die menschliche Physiologie oder auch über ein technisches System geben. Insbesondere die Phasenkontrast-MRT, die auf der Aufnahme von zwei Bildern mit und ohne flusskodierenden Gradienten basiert, hat sich als ein wichtiges diagnostisches Werkzeug in
der Medizin entwickelt, um den Blutfluss funktionell zu quantifizieren. Die klinische Nutzung ist jedoch durch die langen Messzeiten eingeschränkt, da die Daten über mehrere Herzzyklen gemittelt werden müssen und damit die Untersuchung unmittelbarer physiologischer Reaktionen
auf Stress und/oder Muskelbelastung ausgeschlossen ist. Ein primäres Ziel dieser Arbeit war es, diese Einschränkungen durch die Entwicklung einer MRT-Flussmessung in Echtzeit zu überwinden.
Entscheidende Grundlage jeder Echtzeit-MRT sind kurze Aufnahmezeiten, die vor allem durch eine Reduktion der aufgenommenen Daten (Unterabtastung) realisiert werden. Daher konzentrierte sich die hier vorgestellte Entwicklung auf die damit verbundenen technischen Aspekte wie die MRT-Sequenz zur Datenaufnahme, das räumliche Kodierungsschema, und die Bildrekonstruktion.
Experimentell erwies sich ein radiales Kodierungsschema als robust gegenüber Bewegungen und relativ unempfindlich gegenüber milder Unterabtastung. Dieses Kodierungsschema wurde mit der FLASH Aufnahmetechnik für eine effiziente Phasenkontrast-Sequenz in Echtzeit kombiniert. Zusätzlich wurde die Sequenz durch Überlappung von Gradienten hinsichtlich einer kurzen Echozeit optimiert. Für die Bildrekonstruktion wurde die regularisierte nichtlineare inverse Rekonstruktion (NLINV) verwendet, bei der die Bildinformation und die entsprechenden pulensensitivitäten gleichzeitig geschätzt werden. NLINV wurde speziell für die Phasenkontrast-MRT angepasst, um sowohl Betragsbilder als auch robuste Phasenkontrast-Karten mit hoher raumzeitlicher Genauigkeit zu berechnen. Das erarbeitete Verfahren der Phasenkontrast-MRT in Echtzeit kombiniert daher zwei stark unterabgetastete (bis zu einem Faktor von 30) und unterschiedlich flusskodierte, radiale Gradientenecho-Aufnahmen mit einer modifizierten NLINV Rekonstruktion. Mit dieser Methode wurde sowohl eine gute räumliche Auflösung (1.3 mm), als auch eine hohe zeitliche Auflösung (40 ms) erreicht.
Bei Anwendungen an gesunden Probanden sowie vorläufigen Untersuchungen von Patienten konnte nachgewiesen werden, dass die Phasenkontrast-MRT in Echtzeit einen verbesserten
Komfort für die Patienten (z.B. freie Atmung) und unmittelbaren Zugang zu physiologischen Veränderungen der Flussparameter bietet (z.B. Reaktion auf erhöhten Druck im Brustraum). In den meisten Fällen wurden quantitative Blutflussmessungen in der aufsteigenden Aorta, einem
klinisch wichtigen Gefäß des Herz-Kreislauf-Systems, vorgenommen. Die Messungen mit der Phasenkontrast-MRT in Echtzeit wurden mit der EKG-getriggerten Cine Phasenkontrast-MRT (klinischer Standard) an einem Flussphantom und unter in vivo Bedingungen verglichen. Die
Ergebnisse zeigten unter vergleichbaren Bedingungen gute Übereinstimmung.
Im Rahmen dieser Arbeit wurde zusätzlich eine Doppelecho-Variante der Phasenkontrast-MRT in Echtzeit implementiert. Das Einfügen eines flusskodierenden Gradienten zwischen den beiden Echos führte zu einer weiteren Reduzierung der Messzeit. Obwohl sich diese Technik unter in vitro Bedingungen als tauglich erwies, zeigten sich bei in vivo Studien störende Einflüsse
durch bewegungsinduzierte Phasenbeiträge, die wenig Erfolg für quantitative Flussmessungen versprechen.
Als weitere Methode wurde in dieser Arbeit eine dreifach kodierte Sequenz zur Phasenkontrast-MRT entwickelt, um multidirektionalen Fluss zu untersuchen. Die Geschwindigkeitskodierung entlang aller drei physikalischen Gradientenrichtungen führte zu einer verlängerten Messzeit (zeitliche Auflösung � 100 ms) gegenüber der Echtzeit-Flussmessung in nur einer Richtung.
Obwohl das Verfahren noch zu langsam ist, um den kardiovaskulären Blutfluss adäquat zu beschreiben, waren vorläufige Ergebnisse in Körperregionen außerhalb des Herzens für zukünftige klinische Anwendungen sehr vielversprechend. Es ist zu erwarten, dass entsprechende
Weiterentwicklungen von modellbasierten ekonstruktionsverfahren profitieren werden.
|
6 |
Hybridmethoden zur Reduzierung der spezifischen Absorptionsrate für neuroradiologische MRT-Untersuchungen an Hochfeldsystemen / Hybrid methods for reducing specific Absorption rate in neuroradiologic MRI-examinations at high-field MR-systemsCholi, Morwan January 2013 (has links) (PDF)
Die klinische Magnetresonanztomografie (MRT) operiert meist bei einer Magnetfeldstärke von 1,5 Tesla (T). Es halten jedoch immer mehr 3T MRT-Systeme Einzug
im klinischen Alltag und seit kurzem auch 7T Ganzkörper-MRT-Systeme in die Grundlagenforschung. Höhere Magnetfeldstärken führen grundsätzlich zum einem
verbesserten Signal-zu-Rausch- Verhältnis, welches sich gewinnbringend in eine erhöhte Ortsauflösung oder schnellere Bildaufnahme äußert. Ein Nachteil ist aber die dabei im Patienten deponierte Hochfrequenz-Energie (HF-Energie), welche quadratisch mit ansteigender Feldstärke zusammenhängt. Charakterisiert wird diese durch die
spezifische Absorptionsrate (SAR) und ist durch vorgegebene gesetzliche Grenzwerte beschränkt. Moderne, SAR-intensive MRT-Techniken (z.B. Multispinecho-Verfahren)
sind bereits bei 1,5T nahe den zulässigen SAR-Grenzwerten und somit nicht unverändert auf Hochfeld-Systeme übertragbar. In dieser Arbeit soll das Potential modularer
Hybrid-MRT-Techniken genutzt werden, um das SAR bei besonders SAR-intensiven MRT-Verfahren ohne signifikante Einbußen in der Bildqualität erheblich zu verringern. Die Hybrid-Techniken sollen in Verbindung mit zusätzlichen Methoden der SAR-Reduzierung den breiteren Einsatz SAR-intensiver MRT-Techniken an hohen Magnetfeldern ermöglichen. Ziel dieser Arbeit ist es, routinefähige und SAR-reduzierte MRT-Standard-Protokolle für neuroanatomische Humanuntersuchungen mit räumlicher Höchstauflösung bei Magnetfeldern von 3T und 7T zu etablieren. / Spin echo based MRI sequences builds one of the main priorities in the medical/-morphological MRI imaging. Especially T2-weighted spin-echo sequences and the
multi-spin-echo RARE imaging sequence represents one of the basic methods, which allows to minimize measurement times down to minutes or seconds. Modern MRI systems usually have a magnetic field strength beyond 1.5T. In
order to acquire high resolution images with acceptable acquisition times, high RF power is needed. In many cases the intrinsic safety limits of the radiated power
are already exeeded at field strengths of 3T. This leads often to restrictions for the full performance. In research systems of 7T and more T2-weighted images are only feasible with significant amount of time. Especially in RARE imaging methods, which require a large number of refocusing pulses, this is a significant restriction factor.
Therefore, in recent years there were further development of hybrid sequences which combines different acquisition methods together into one, to exploit their fully
advantages of the basic methods and to minimize ...
|
7 |
Advances in Fast MRI Experiments / Neue Methoden in der MR-BildgebungNeumann, Daniel January 2014 (has links) (PDF)
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. / Die Magnetresonanztomographie (MRT) ist ein wichtiges nicht-invasives medizinisches Bildge- bungsverfahren, das im klinischen Alltag zur Entdeckung und Diagnose einer Vielzahl von Krank- heiten verwendet wird. Im Gegensatz zu anderen Methoden wie Röntgen und Computertomo- graphie kommt die MRT ohne den Einsatz ionisierender Strahlung aus, was selbst häufige An- wendungen ohne gesundheitliche Risiken erlaubt. Einer der größten Nachteile der MRT sind lange Messzeiten, die in Kombination mit der teuren Technik hohe Untersuchungskosten bedin- gen. Obwohl in der Vergangenheit durch die Entwicklung von sowohl verbesserter Hardware als auch neuen Rekonstruktionsverfahren bereits bedeutende Fortschritte in Bezug auf die Akquisi- tionsdauer erzielt werden konnten, ist eine weitere Beschleunigung nach wie vor ein wichtiges Forschungsgebiet im Bereich der MRT. Ziel dieser Arbeit war daher die Entwicklung neuer An- sätze zur Steigerung der Effizienz von MRT Experimenten.
In Kapitel 4 wurde eine kombinierte Akquisitions- und Rekonstruktionsstrategie für radiale Turbo Spin-Echo (TSE) Experimente vorgestellt. Im klinischen Alltag sind kartesische TSE Sequenzen zur Untersuchung diverser Krankheitsbilder weit verbreitet, da sie ein hohes SNR aufweisen und die Aufnahme der klinisch wichtigen Bilder mit Protonendichte- und T2-Kontrast erlauben. Im Gegensatz zu kartesischem Abtasten, wo aus einem Datensatz lediglich ein Bild mit bes- timmtem Kontrast erzeugt wird, sind radiale TSE Akquisitionen hocheffizient, da hier aus einem Datensatz mehrere Bilder mit verschiedenem Kontrast gewonnen werden können. In beiden Fällen wird in konventionellen Rekonstruktionsmethoden jedes Bild eines definierten Kontrasts durch das Zusammensetzen eines vollständig abgetasteten k-Raums mit Daten von verschiedenen Echozeiten erzeugt. Im radialen Fall geschieht dies durch die sogenannte "k-space weighted im age contrast" (KWIC) Rekonstruktion. Hierbei wird das Zentrum des k-Raums ausschließlich mit zum gewünschten Kontrast gehörigen Daten gefüllt, während die äußeren Bereiche des k-Raums auch Daten von anderen Echozeiten enthalten. Obwohl der Kontrast von MRT Bildern haupt- sächlich von den Daten im k-Raum Zentrum dominiert wird, führt die Kombination von Daten verschiedener Echozeiten in sowohl radialen als auch kartesischen TSE Bildern zu einem uner- wünschten Mischkontrast. Dieser Effekt wird vor allem in protonendichtegewichteten Bildern sichtbar und kann somit deren diagnostischen Wert deutlich verringern.
Ein unerwünschter Mischkontrast kann verhindert werden, indem die Bandbreite des KWIC- Filters auf ein kleines zeitliches Fenster um die angestrebte Echozeit herum eingeschränkt wird. Um eine freie Wahl der Echozuglänge und hohe Flexibilität in der Bildrekonstruktion zu er- möglichen, wurde für die radiale TSE Akquisition ein angepasstes Abtastschema unter Verwen- dung des goldenen Winkels vorgestellt. Da bei einem KWIC-Filter mit reduzierter Bandbre- ite für jedes Bild weniger Daten zur Verfügung stehen als für einen voll abgetasteten k-Raum benötigt, weisen rekonstruierte Bildern Einfaltungsartefakte auf. Diese werden in einem zweiten Schritt durch die Anwendung paralleler Bildgebung beseitigt. In den gezeigten Beispielen wurde dazu der CG-SENSE Algorithmus verwendet, da er stabile Konvergenz aufweist und für die Rekonstruktion von Daten mit irregulären Abtastschema angewandt werden kann. Anschließend werden bestehende Korrelationen der Bilderserie zur Reduktion verbleibender Artefakte und zu einer Verbesserung des SNR ausgenutzt. Wie mittels Simulationen gezeigt und für neurologische Daten bestätigt, weisen radiale TSE Bilder, die mit dieser Methode rekonstruiert wurden, keinen sichtbaren Mischkontrast mehr auf. Die erreichte Bildqualität ist hierbar vergleichbar mit kon- ventionellen Rekonstruktionsmethoden. Da die vorgestellte Rekonstruktion einfach auf heutigen Scannern implementiert werden kann und lediglich niedrige Rechenkapazitäten benötigt, könnte sie einen wichtigen Beitrag für den klinischen Durchbruch radialer TSE Akquisitionen darstellen.
In Kapitel 5 wurde eine selbstkalibrierende Methode zur Korrektur von aus Multi Spin-Echo (MSE) Bildern gewonnenen T2 Karten vorgestellt. In der klinischen Anwendung spielt die Quan- tifizierung von T2 unter anderem bei der Diagnose von Krankheiten sowie bei der Klassifizierung von Gewebe eine wichtige Rolle. Eine MSE Sequenz verwendet mehrere RF-Pulse, um ein einzelnes Spin-Echo wiederholt zu refokussieren. Idealerweise betragen die Flipwinkel der Re- fokussierungspulse hierbei exakt 180○, um einen exponentiellen Signalabfall zu erhalten. Auf- grund technischer Ungenauigkeiten weichen die Werte der Flipwinkel von Refokussierungspulsen jedoch grundsätzlich von 180○ ab. Niedrigere Flipwinkel führen zu stimulierten Echos, die wesentlich zu den einzelnen Echoamplituden beitragen und den Signalabfall entlang des Echozugs deutlich verlängern können. Somit weisen auch T2 Werte, die aus solchen Bilderserien berech- net werden, eine teilweise deutliche Erhöhung auf. Um exakte Werte zu erhalten, kann der MSE Signalverlauf mittels des "extended phase graph" (EPG) Algorithmus abgeschätzt und so ein Kor- rekturfaktor ermittelt werden. Hierzu müssen die Flipwinkel der Refokussierungspulse, T1 und T2 Werte sowie der zeitliche Abstand der Echos (ESP) bekannt sein. Wie in Simulationen gezeigt wurde, kann T1 für hohe Werte des Quotienten T1/T2 abgeschätzt werden, ohne an Genauigkeit der T2 Ergebnisse einzubüßen. Abschätzungen der verbleibenden benötigten Parameter können direkt aus den Daten selbst gewonnen werden. Während der Flipwinkel aus den Intensitäten der ersten beiden Echos berechnet wird, liefert ein mono-exponentieller Fit der MSE Bilderserie eine erste Näherung für T2. Die Korrektur für die T2 Werte kann anschließend aus den EPG Sig- nalverläufen berechnet werden. Durch Aktualisierung von T2 und erneuter Ausführung des EPG-Algorithmus wird die Genauigkeit der Korrektur iterativ erhöht, wobei schon eine sehr geringe Zahl von Iterationen zu Konvergenz führt. Wie in Simulationen und in Phantomexperimenten für verschiedenste T2-Werte gezeigt, weisen korrigierte T2 Werte eine hohe Genauigkeit auf. Dies gilt auch für niedrigere nominelle Flipwinkel als 180○ und ist somit von speziellem Interesse bei höheren Feldstärken B0, wo Grenzwerte der spezifischen Absorptionsrate die Einstrahlung einer Vielzahl von RF-Pulsen hoher Amplitude verbietet.
In Kapitel 6 wurde iteratives GROG, eine neue selbstkalibrierende iterative parallele Bildge- bungsmethode für die Rekonstruktion von nichtkartesischen Daten vorgestellt. Es sind eine Vielzahl nichtkartesischer Trajektorien für MRT Messungen bekannt, die zahlreiche Vorteile gegenüber kartesischer Bildgebung bieten. Dazu gehören unter anderem eine schnellere Akquisi- tion, verbesserte dynamische Bildgebung sowie die Möglichkeit zur intrinischen Bewegungskor- rektur. Ein Nachteil nichtkartesischer Daten jedoch ist eine aufwendigere Rekonstruktion, sowohl bei voll abgetasteten Datensätzen als insbesondere auch in der parallelen Bildgebung. Itera- tives GROG verwendet Calibrationless Parallel Imaging by Structured Low-Rank Matrix Com- pletion (CPI) zur Rekonstruktion fehlender Daten. Diese Methode benötigt Daten auf karte- sischen Gitterpunkten und kann nicht direkt für nichtkartesische Experimente angewandt wer- den. Stattdessen werden die nichtkartesischen Daten zunächst mittels Grappa Operator Gridding (GROG) in einem ersten Schritt auf ein kartesisches Gitter verschoben. GROG basiert auf paral- leler Bildgebung und benötigt einen voll abgetasteten Teil des k-Raums zur Kalibrierung. Erste Kalibrationsdaten können gewonnen werden, indem die nichtkartesischen Punkte ohne Änderung auf die nächsten kartesischen Gitterpunkte verschoben werden und eine CPI-Rekonstruktion eines zentralen k-Raum Bereichs durchgeführt wird. Anschließend wird GROG angewandt um exakte Werte der kartesischen Gitterpunkte zu erhalten und der Prozess wird iteriert. Nach Kon- vergenz können entweder Kalibrationsdaten gewünschter Größe für eine konventionelle parallele Bildgebungsmethode erzeugt oder artefaktfreie Bilder mit CPI rekonstruiert werden.
In Kapitel 7 wurde ein neues Abtastungsschema für die 2D Bildgebung vorgestellt, das aus Multiplen Oszillierenden Effizienten Trajektorien (MOET) besteht und optimierte Compressed Sensing (CS) Rekonstruktionen ermöglicht. Für eine erfolgreiche Anwendung von Compressed Sensing müssen einige Voraussetzungen erfüllt sein. Erstens muss jeder Datenpunkt Informa- tionen über das ganze Objekt enthalten, was bei der MRT aufgrund der Datenakquisition im Fourier-Raum automatisch erfüllt ist. Weiterhin muss das gemessene Objekt in einer beliebigen Basis sparse sein. Dies ist für viele verschiedene Anwendungen in der MRT der Fall. Drittens muss für CS Rekonstruktionen die Datenakquisition im k-Raum einem inkohärenten Muster fol- gen. Diese wichtige Voraussetzung ist in der zweidimensionalen Bildgebung mit konventionellen kartesischen und nicht-kartesischen Abtastschemata nur schwer zu erreichen. Deshalb wird für CS Rekonstruktionen häufig eine radiale Trajektorie eingesetzt, trotz der kohärenten streaking- Artefakte in unterabgetasteten Bildern. MOET verwendet daher eine Kombination von radialen Projektionen zusammen mit oszillierenden Gradienten auf der zur Ausleserichtung orthogonalen Achse. So erhält man inkohärente Aliasing-Artefakte und bewahrt gleichzeitig die Vorteile der radialen Bildgebung für die dynamische MRT. Die Überlegenheit von MOET gegenüber radi- aler Bildgebung für CS Rekonstruktionen konnte in Simulationen sowie in der Herzbildgebung aufgezeigt werden. Dies gilt insbesondere für CS Rekonstruktionen direkt im Bildraum, wo MOET gute Resultate liefert während die kohärenten Artefakte bei radialer Bildgebung eine genaue Bildwiederherstellung verhindert. Bei Rekonstruktionen dynamischer Daten, wo Sparsität im xf-Raum ausgenutzt wird, ist der Vorteil von MOET weniger ausgeprägt, da hier bere- its die zeitliche Anordnung der Projektionen einen wesentlichen Beitrag zur Charakteristik der Aliasingartefakte liefert. Wie in Simulationen und für die in vivo Herzbildgebung gezeigt werden konnte, erlaubt die zusätzliche räumliche Inkohärenz von MOET jedoch auch in diesem Fall eine höhere Genauigkeit sowie Rekonstruktionen von Daten höherer Beschleunigung.
|
8 |
Neue Steady-State-Techniken in der Magnetresonanztomographie / Novel Steady-State Techniques for Magnetic Resonance ImagingBenkert, Thomas January 2015 (has links) (PDF)
Die bSSFP-Sequenz kombiniert kurze Akquisitionszeiten mit einem hohen Signal-zu-Rausch-Verhältnis, was sie zu einer vielversprechenden Bildgebungsmethode macht. Im klinischen Alltag ist diese Technik jedoch bisher - abgesehen von vereinzelten Anwendungen - kaum etabliert. Die Hauptgründe hierfür sind Signalauslöschungen in Form von Bandingartefakten sowie der erzielte T2/T1-gewichtete Mischkontrast. Das Ziel dieser Dissertation war die Entwicklung von Methoden zur Lösung der beiden genannten Limitationen, um so eine umfassendere Verwendung von bSSFP für die MR-Diagnostik zu ermöglichen.
Magnetfeldinhomogenitäten, die im Wesentlichen durch Suszeptibilitätsunterschiede oder Imperfektionen seitens der Scannerhardware hervorgerufen werden, äußern sich bei der bSSFP-Bildgebung in Form von Bandingartefakten. Mit DYPR-SSFP (DYnamically Phase-cycled Radial bSSFP) wurde ein Verfahren vorgestellt, um diese Signalauslöschungen effizient zu entfernen. Während für bereits existierende Methoden mehrere separate bSSFP-Bilder akquiriert und anschließend kombiniert werden müssen, ist für die Bandingentfernung mittels DYPR-SSFP lediglich die Aufnahme eines einzelnen Bildes notwendig. Dies wird durch die neuartige Kombination eines dynamischen Phasenzyklus mit einer radialen Trajektorie mit quasizufälligem Abtastschema ermöglicht. Die notwendigen Bestandteile können mit geringem Aufwand implementiert werden. Des Weiteren ist kein spezielles Rekonstruktionsschema notwendig, was die breite Anwendbarkeit des entwickelten Ansatzes ermöglicht. Konventionelle Methoden zur Entfernung von Bandingartefakten werden sowohl bezüglich ihrer Robustheit als auch bezüglich der notwendigen Messzeit übertroffen.
Um die Anwendbarkeit von DYPR-SSFP auch jenseits der gewöhnlichen Bildgebung zu demonstrieren, wurde die Methode mit der Fett-Wasser-Separation kombiniert. Basierend auf der Dixon-Technik konnten so hochaufgelöste Fett- sowie Wasserbilder erzeugt werden. Aufgrund der Bewegungsinsensitivät der zugrunde liegenden radialen Trajektorie konnten die Messungen unter freier Atmung durchgeführt werden, ohne dass nennenswerte Beeinträchtigungen der Bildqualität auftraten. Die erzielten Ergebnisse am Abdomen zeigten weder Fehlzuordnungen von Fett- und Wasserpixeln noch verbleibende Bandingartefakte.
Ein Nachteil der gewöhnlichen Dixon-basierten Fett-Wasser-Separation ist es, dass mehrere separate Bilder zu verschiedenen Echozeiten benötigt werden. Dies führt zu einer entsprechenden Verlängerung der zugehörigen Messzeit. Abhilfe schafft hier die Verwendung einer Multiecho-Sequenz. Wie gezeigt werden konnte, ermöglicht eine derartige Kombination die robuste, bandingfreie Fett-Wasser-Separation in klinisch akzeptablen Messzeiten.
DYPR-SSFP erlaubt die Entfernung von Bandingartefakten selbst bei starken Magnetfeldinhomogenitäten. Dennoch ist es möglich, dass Signalauslöschungen aufgrund des Effekts der Intravoxeldephasierung verbleiben. Dieses Problem tritt primär bei der Bildgebung von Implantaten oder am Ultrahochfeld auf. Als Abhilfe hierfür wurde die Kombination von DYPR-SSFP mit der sogenannten z-Shim-Technik untersucht, was die Entfernung dieser Artefakte auf Kosten einer erhöhten Messzeit ermöglichte.
Die mit DYPR-SSFP akquirierten radialen Projektionen weisen aufgrund des angewendeten dynamischen Phasenzyklus leicht verschiedene Signallevel und Phasen auf. Diese Tatsache zeigt sich durch inkohärente Bildartefakte, die sich jedoch durch eine Erhöhung der Projektionsanzahl effektiv reduzieren lassen. Folglich bietet es sich in diesem Kontext an, Anwendungen zu wählen, bei denen bereits intrinsisch eine verhältnismäßig hohe Anzahl von Projektionen benötigt wird. Hierbei hat sich gezeigt, dass neben der hochaufgelösten Bildgebung die Wahl einer 3D radialen Trajektorie eine aussichtsreiche Kombination darstellt. Die in der vorliegenden Arbeit vorgestellte 3D DYPR-SSFP-Technik erlaubte so die isotrope bandingfreie bSSFP-Bildgebung, wobei die Messzeit im Vergleich zu einer gewöhnlichen bSSFP-Akquisition konstant gehalten werden konnte. Verbleibende, durch den dynamischen Phasenzyklus hervorgerufene Artefakte konnten effektiv mit einem Rauschunterdrückungsalgorithmus reduziert werden. Anhand Probandenmessungen wurde gezeigt, dass 3D DYPR-SSFP einen aussichtsreichen Kandidaten für die Bildgebung von Hirnnerven sowie des Bewegungsapparats darstellt.
Während die DYPR-SSFP-Methode sowie die darauf beruhenden Weiterentwicklungen effiziente Lösungen für das Problem der Bandingartefakte bei der bSSFP-Bildgebung darstellen, adressiert die vorgestellte RA-TOSSI-Technik (RAdial T-One sensitive and insensitive Steady-State Imaging) das Problem des bSSFP-Mischkontrasts. Die Möglichkeit der Generierung von T2-Kontrasten basierend auf der bSSFP-Sequenz konnte bereits in vorausgehenden Arbeiten gezeigt werden. Hierbei wurde die Tatsache ausgenutzt, dass der T1-Anteil des Signalverlaufs nach Beginn einer bSSFP-Akquisition durch das Einfügen von Inversionspulsen in ungleichmäßigen Abständen aufgehoben werden kann. Ein so akquiriertes Bild weist folglich einen reinen, klinisch relevanten T2-Kontrast auf. Die im Rahmen dieser Arbeit vorgestellte Methode basiert auf dem gleichen Prinzip, jedoch wurde anstelle einer gewöhnlichen kartesischen Trajektorie eine radiale Trajektorie in Kombination mit einer KWIC-Filter-Rekonstruktion verwendet. Somit können bei gleichbleibender oder sogar verbesserter Bildqualität aus einem einzelnen, mit RA-TOSSI akquirierten Datensatz verschiedene T2-Wichtungen als auch gewöhnliche T2/T1-Wichtungen generiert werden. Mittels Variation der Anzahl der eingefügten Inversionspulse konnte ferner gezeigt werden, dass es neben den besagten Wichtungen möglich ist, zusätzliche Kontraste zu generieren, bei denen verschiedene Substanzen im Bild ausgelöscht sind. Diese Substanzen können am Beispiel der Gehirnbildgebung Fett, graue Masse, weiße Masse oder CSF umfassen und zeichnen sich neben den reinen T2-Kontrasten durch eine ähnlich hohe klinische Relevanz aus. Die mögliche Bedeutung der vorgestellten Methode für die klinische Verwendung wurde durch Messungen an einer Gehirntumorpatientin demonstriert.
Zusammenfassend lässt sich sagen, dass die im Rahmen dieser Dissertation entwickelten Techniken einen wertvollen Beitrag zur Lösung der eingangs beschriebenen Probleme der bSSFP-Bildgebung darstellen. Mit DYPR-SSFP akquirierte Bilder sind bereits mit bestehender, kommerzieller Rekonstruktionssoftware direkt am Scanner rekonstruierbar. Die Software für die Rekonstruktion von RA-TOSSI-Datensätzen wurde für Siemens Scanner implementiert. Folglich sind beide Methoden für klinische Studien einsetzbar, was gleichzeitig den Ausblick dieser Arbeit darstellt. / The bSSFP sequence combines short acquisition times with a high signal-to-noise ratio, and is therefore a promising imaging technique. However, aside from a few applications, this method is hardly established in the clinical routine. The main reasons for this are signal voids that arise as banding artifacts and the obtained T2/T1-weighted mixed contrast. The goal of this dissertation was to develop strategies to overcome these limitations and allow for a more widespread use of bSSFP for MR diagnostics.
In bSSFP imaging, magnetic field inhomogeneities, which are mainly caused by susceptibility differences and imperfections of the scanner hardware, manifest as banding artifacts. In order to efficiently remove these artifacts from the image, DYnamically Phase-cycled Radial bSSFP (DYPR-SSFP) was proposed. While existing methods rely on the acquisition and subsequent combination of several separate bSSFP images, banding removal with DYPR-SSFP requires the acquisition of only a single data set. This is achieved by combining a dynamic phase-cycle with a radial trajectory and a quasi-random acquisition scheme. The individual components of this method can be implemented with little effort. Furthermore, no specific reconstruction scheme is necessary, guaranteeing the broad applicability of the developed approach. DYPR-SSFP outperformed conventional methods for banding removal both in robustness and scan time.
In order to demonstrate the applicability of DYPR-SSFP beyond conventional imaging, the method was also applied to fat-water separation. Based on the Dixon technique, fat and water images were generated with high resolution. Due to the motion robustness of the underlying radial trajectory, measurements could be performed during free-breathing, without notable degradation of image quality. Abdominal images showed neither regional fat-water flipping nor residual banding artifacts.
A drawback of standard Dixon-based fat-water separation is the fact that several separate images with different echo times have to be acquired, therefore prolonging the respective scan time. This can be overcome by using a multiecho sequence. It was demonstrated that the combination of such multiecho sequence and Dixon DYPR-SSFP allows for robust, banding-free fat-water separation in clinically acceptable scan times.
DYPR-SSFP guarantees removal of banding artifacts even for strong magnetic field inhomogeneities. However, signal voids may remain due to intravoxel dephasing. This problem primarily arises when imaging metallic implants or when moving to ultra-high field strengths. To address this issue, the combination of DYPR-SSFP with the so-called z-shim technique was investigated, allowing the removal of these artifacts at the expense of an increased measurement time.
Due to the applied dynamic phase-increment, radial projections which are acquired with DYPR-SSFP exhibit slightly different signal levels and phases. This results in incoherent artifacts, that can be effectively reduced by increasing the number of acquired projections. Therefore, DYPR-SSFP should be preferably applied when many projections are intrinsically necessary. It has been demonstrated that, besides high resolution imaging, the choice of a 3D radial trajectory is a promising combination. The proposed 3D DYPR-SSFP technique allowed isotropic banding-free bSSFP imaging without any expense of additional scan time compared to a conventional bSSFP acquisition. Residual artifacts caused by the dynamic phase-cycle could be effectively mitigated by applying a denoising algorithm. Volunteer measurements showed that 3D DYPR-SSFP is a promising candidate for imaging of the cranial nerves and the musculoskeletal system.
While DYPR-SSFP and all presented resulting methods constitute an efficient solution for banding artifacts in bSSFP imaging, the proposed RAdial T-One sensitive and insensitive Steady-State Imaging (RA-TOSSI) method addresses the problem of the mixed contrast in bSSFP imaging. The possibility to generate T2-contrast with bSSFP has been shown before. The previous approach is based on the fact that T1-relaxation during the transient phase of a bSSFP acquisition can be suppressed by inserting unequally spaced inversion pulses. Thus, the resulting image shows a clinically relevant T2-contrast. The method which was presented as part of this dissertation relies on the same principle. However, instead of the originally proposed Cartesian trajectory, a radial trajectory in combination with a KWIC-filter reconstruction was applied. This allows the generation of several T2-weighted images as well as T2/T1-weighted images from a single RA-TOSSI dataset, while image quality remains comparable or even improves compared with existing methods. It could further be shown that varying the number of inversion pulses allows the generation of additional contrasts, where different tissue types are attenuated in the image. In the case of brain imaging for instance, these tissues comprise fat, gray matter, white matter, and CSF and offer high clinical relevance similar to T2-weighted images. Measurements of a patient with a brain tumor demonstrate the possible impact of the proposed method.
In conclusion, the techniques developed as part of this dissertation present a valuable contribution to the solution of various problems which are associated with bSSFP imaging. Images acquired with DYPR-SSFP can be reconstructed directly at the scanner using existing, commercial reconstruction software. The software for the reconstruction of RA-TOSSI data was implemented for Siemens scanners. Therefore, both methods can be directly employed for clinical studies which remain as future work.
|
9 |
Differential imaging characteristics and dissemination potential of pilomyxoid astrocytomas versus pilocytic astrocytomas / Pilomyxoides Astrozytom und pilozytisches Astrozytom: Vergleich anhand der BildgebungAlkonyi, Balint January 2014 (has links) (PDF)
Background and Aims: PMA is a recently described rare tumor entity occuring most often in young children. Due the worse outcome of PMA-patients as compared to children with pilocytic astrocytoma (PA), it has received a grade II assignment in the latest WHO classification. Nevertheless, increasing evidence suggests that the two tumor types are indeed pathologically and genetically related. The radiological differentiation of PMAs from PAs is challenging and the limited available data could not yet provide unequivocal distinguishing imaging features. Furthermore, it is not completely clarified whether PMA cases are associated with a higher rate of CSF dissemination compared to similarly young patients with PA. The aim of our study was firstly to compare MR/CT imaging features of these tumors, and secondly, to evaluate the occurrence of CSF dissemination.
Material and Methods: The study population included 15 children with PMA and 32 children with PA. A third group consisted of eight children with PAs with focal pilomyxoid features. All cases had been registered in the German multicenter SIOP/HIT-LGG trials. The initial MRIs (and CT scans, if available) at establishing the diagnosis were retrospectively analyzed according to standardized criteria and the findings compared between PMAs and PAs. Furthermore, we compared the occurrence of imaging evidences of CSF tumor dissemination between children with PMA and PA, respectively.
Results: The imaging appearance of PMAs and PAs was very similar. However, PAs tended to show more frequently cystic components (p=0.03). As opposed to PAs, PMAs did not have large tumor cysts. We did not find differences with respect to tumor size and tumor margin. Gadolinium enhancement of PMAs was significantly more frequently homogeneous (p=0.006). PMAs appeared to show more often intratumoral hemorrhages (p=0.047). Furthermore, suprasellar PMAs tended to have a more homogeneus texture on T2-weighted MR images (p=0.026). Within the subgroup < 6 years of age the PMA histology tended to have a larger effect on the occurrence of CSF dissemination than the age (p=0.05 vs.0.12).
Conclusions: Although the radiological appearance of PMAs and PAs is similar, some imaging features, like enhancement pattern or presence of cysts or hemorrhage may help differentiating these low-grade gliomas. Our results corroborate previous scarce data suggesting higher rate of CSF dissemination in PMAs, even in the youngest patient population. Thus, in young children with a chiasmatic-hypothalamic tumor suggestive of a PMA, an intensive search for CSF dissemination along the entire neuraxis should be performed. / Hintergrund und Ziele: Das PMA ist eine seltene benigne Tumorentität (WHO Grad II nach der Klassifikation von 2007), welche vorwiegend bei Kleinkindern vorkommt. Dem PMA wird eine schlechtere Prognose als dem PA zugeschrieben. Die radiologische Unterscheidung zwischen PMA und PA erscheint schwierig. Zudem ist unklar inwieweit das PMA mit einer höheren Meningeoserate assoziiert sein könnte, verglichen zum PA bei Patienten vergleichbaren Alters. Das Ziel unserer Studie war es, bildgebende Kriterien sowie das Auftreten einer Liquordissemination zwischen beiden Tumorarten zu vergleichen.
Material und Methoden: Das untersuchte Kollektiv umfasste 15 Kinder mit PMA sowie 32 Kinder mit PA. Zudem eine dritte Gruppe umfasste 8 Kinder mit PA mit fokalen pilomyxoiden Merkmale. Die initialen MRT-Untersuchungen bei Diagnosestellung wurden anhand standardisierter bildgebender Kriterien analysiert. Anhand der MRT-Bildgebung wurde des Weiteren das Auftreten einer Liquordissemination bei Diagnose und im Verlauf verglichen.
Ergebnisse: PMAs und PAs stellten sich bildmorphologisch ähnlich dar. PAs zeigten öfters zystische Tumorkomponente (p=0.03). Im Gegensatz zu PAs, wiesen PMAs keine großen Tumorzysten auf. Auch die Kontrastmittelanreicherung war häufiger homogen in PMAs (p=0.006). Intratumorale Hämorrhagien konnten häufiger in PMAs beobachtet werden (p=0.047). Desweiteren haben sich supraselläre PMAs in der T2-Wichtung häufiger homogener dargestellt (p=0.026). Die histologische Diagnose PMA schien einen größeren Effekt (p=0.05) auf das Auftreten einer Liquordissemination zu haben als das Alter selbst (p=0.12).
Schlussfolgerungen: Obwohl sich PMA und PA in der Bildgebung sehr ähnlich darstellen, identifizierten wir einige Kriterien welche beitragen könnten diese Tumorentitäten zu unterscheiden. Bei radiomorphologischen Kriterien, die für ein PMA sprechen, sollte bei den jungen Kindern gezielt nach einer Liquordisseminierung gesucht werden.
|
10 |
Magnetresonanztomographie in der Zahnheilkunde - hochauflösende zahnmedizinische Anwendungen in der MRT mit einer Entwicklung zur Bewegungskorrektur / Magnetic Resonance Imaging in Dentistry – high-resolution dental applications in MRI with development of a method for motion correctionHopfgartner, Andreas January 2015 (has links) (PDF)
Die zahnmedizinische Behandlung von Erkrankungen der Zähne oder im Bereich der Mundhöhle erfolgt bei Weitem nicht immer aus optischen Gründen. Diese Erkrankungen werden auch mit ernsthaften Erkrankungen in Zusammenhang gebracht. Studien haben gezeigt, dass einige Erkrankungen im Mund- und Zahnbereich zu Herz- und Lungenkrankheiten oder Diabetes führen können. Oftmals erstreckt sich die Pathologie oder Symptomatik von Mund- und Zahnerkrankungen über einen weiten Bereich. In der zahnmedizinischen Klinik kommen daher viele verschiedene diagnostische Apparate zum Einsatz. Allerdings zählt die Magnetresonanztomographie, die sich in anderen Bereichen bereits zum wichtigsten bildgebenden Diagnosetool entwickelt hat, dort noch nicht zu den Standardverfahren. Dabei liegen ihre Vorteile auf der Hand: sie ist bekannt für sehr gute Bildkontraste vor allem zwischen verschiedenen Weichgewebsarten und kommt ohne gefährliche ionisierende Strahlung aus. Wahrscheinlich ist ersteres der Grund, warum die MRT in der Zahnmedizin noch nicht sonderlich vertreten ist, kommt es dort oft auf die kontrastreiche Darstellung von Hartgeweben an. Neueste Entwicklungen und Studien belegen jedoch die vielseitigen Vorteile der MRT auch in diesem Bereich. Ziel dieser Arbeit von der applikativen Seite betrachtet, war es, das enorme Potential der MRT in den vielseitigen Bereichen der Zahnmedizin weiterhin aufzuzeigen. Viele dieser Anwendungen stellen jedoch sehr hohe Anforderungen an die Systeme. Meist sind die darzustellenden Strukturen sehr klein und erfordern eine hohe Auflösung. Während man beim Röntgenverfahren beispielsweise die Energie des Strahles (Dosis) steigern kann, bedeutet dies in der MRT (ohne das Gerät zu wechseln) eine Verlängerung der Messzeit. Gerade im Bereich des Kopfes kommt es oft zu ungewollten Bewegungen, die das Ergebnis und die Reproduzierbarkeit der gewonnenen diagnostischen Informationen verschlechtern oder gänzlich unbrauchbar machen. Die grösste Herausforderung dabei ist die dreidimensionale Abformung von Zahnoberflächen in der Prothetik. Dieses Verfahren kann eine aufwändige und unangenehme manuelle Abformung der Zähne und die Herstellung eines Zwischengipsmodells ersetzen und ein direktes dreidimensionales Modell der Zahnoberflächen produzieren. Durch die moderne CAD-/CAM-Technik kann daraus vom Zahntechniker direkt eine Zahnrestauration erstellt werden. Daher war ein wichtiger Bestandteil des Projekts dentale MRT die Entwicklung einer Methode zur Erkennung und gleichzeitiger Korrektur von Bewegungen.
Verschiedenste Anforderungen waren an die Methode gestellt. Zum einen muss die Methode bereits Bewegungen im Bereich von ~100 µm erkennen, um die Anforderungen an die finale Bildauflösung zu unterschreiten. Bei der dentalen Abformung wird eine 1-Kanal-Empfängerspule verwendet und je nach Messung kann der Patient dabei auf dem Bauch oder Rücken liegen. Weiterhin muss die Bewegungserkennung ohne zusätzliche externe Geräte wie Kameras, deren Sicht z.B. durch den Patienten verdeckt ist, durchführbar sein.
Die vorliegende Arbeit deckt also zwei größere Themenblöcke ab. Zum einen wurden in der Arbeit neue Applikationen entwickelt oder weiterentwickelt, um verschiedenen Bereichen der Zahnmedizin den Zugang zu MRTUntersuchungen zu eröffnen. Kapitel 4 beschreibt die Möglichkeit, die Bewegung des Kiefergelenks dynamisch zu erfassen. Es stellte sich in der Arbeit heraus, dass sowohl die Bewegung von Weichgewebeanteilen darstellbar waren, als auch der intraartikuläre Abstand im Kiefergelenk unter Kaubelastung in Echtzeit vermessen werden konnte. Dabei wurde die Bildgebungssequenz und der zugehörige Rekonstruktionsalgorithmus so entwickelt, dass die Daten flexibel und ohne Vorwissen akquiriert und aufbereitet werden können. Hierbei konnten verschiedenen Pathologien anhand der dynamischen Bilder sichtbar gemacht werden und die dynamische MRT konnte Erkrankungen erkennen, die mit anderen Mitteln nicht sichtbar waren. Die vielen diagnostischen Möglichkeiten, die dadurch entstehen sind bisher noch nicht untersucht und sollten durch großangelegte Studien untersucht und belegt werden.
Kapitel 5 beschreibt die Ergebnisse einer großangelegten Studie im Bereich der dentomaxillären Bildgebung . Die diagnostischen Möglichkeiten der MRT für die kieferorthopädische Anwendung liegen klar auf der Hand. Die typischen Patienten in der Kieferorthopädie sind Kinder und Jugendliche. Die Abwesenheit von gewebsschädigender Strahlung ist hier ein besonderer Vorteil der MRT. Eine Messung dauert zudem nach diversen Weiterentwicklungen der Methode nur noch 2 (bzw. 4) Minuten. Die Auflösung in den gerenderten Bildern beträgt 0.25x0.25x0.5 mm. Mit der Methode konnte unter anderem die Geminisierung einer Zahnwurzel und der Abstand des Zahnmarks zur Zahnoberfläche (Zahnschmelz) dargestellt und vermessen werden.
Kapitel 6 stellt Neuentwicklungen im Bereich der dentalen Abformung von Zahnoberflächen dar. Hier wurde eine neue Methode entwickelt um den Patientenkomfort bei der Messung zu steigern und so Bewegungen im Vorhinein zu unterbinden. Bei der alten Methode liegt der Patient auf dem Bauch und ein großer Teil der Mundhöhle ist mit Kontrastmittel befüllt. Durch die Verwendung einer präparierten Tiefziehschiene kann das Kontrastmittel nun lokal appliziert
werden und eine Messung in Rückenlage das Patienten ist somit problemlos möglich. Die damit verbundene Reproduzierbarkeit der Abformungsergebnisse wäre durch eine großangelegte Studie zu zeigen. Die Hauptaufgabe der vorliegenden Dissertation war es, eine Methode zur Bewegungskorrektur zu entwickeln, die es ohne eine große Anzahl an Zusatzgeräten ermöglicht, die Bewegung eines Subjekts während der Messung zu erfassen und dementsprechend zu korrigieren. Diese neue Methode, gestützt auf einer Messung eines MRT-aktiven Markers der am Subjekt angebracht wird, beruht außer der Verwendung des Markers nur auf MRT-Hardware. Die Methode wird in Kapitel 8 vorgestellt. Da es sich bei der Methode um eine Neuentwicklung handelt, war es in erster Linie wichtig, die Einflüsse der verschiedenen Parameter, die sich auf die Positionierungsgenauigkeit auswirken, abzuschätzen und letzten Endes festzulegen. Dies wurde in mehreren Vorstudien, Experimenten und Computersimulationen abgehandelt. In der Arbeit konnte durch Validierungsexperimente gezeigt werden dass sich mit dem bildbasierten Navigator Bewegungen im Genauigkeitsbereich von ~50 µm (Translation) und ~0.13◦(Rotation) detektieren lassen. Mit den Positionsinformationen lassen sich MRT-Daten retrospektiv korrigieren oder idealerweise das Bildgebungsvolumen in Echtzeit anpassen um Inkonsistenzen in den Daten im Vorhinein vorzubeugen. Durch Bewegung beeinträchtigte in-vivo Daten konnten so mit der Methode korrigiert werden und anhand eines geeigneten Phantoms konnte die Verbesserung der Erkennung von Kanten, wie sie beispielsweise bei der dentalen Abformung angewandt wird, gezeigt werden. Die kontinuierlichen Entwicklungen in den Bereichen Hard-, Software und Algorithmik ermöglichen weitere hochauflösende Anwendungen. In Kapitel 9 sind die Ergebnisse einer Studie gezeigt, die sich mit der Analyse der Handbewegungen während einer Messung beschäftigt. Für eine hochauflösenden Darstellung der Handanatomie bei 7 T ist eine Unterbindung der Handbewegung sehr wichtig. Um ein geeignetes Design für eine Empfängerspule zu entwerfen, die Bewegungen der Hand unterbindet, wurde eine qualitative Bewegungsanalyse der Hand in mehreren verschiedenen Positionen durchgeführt. Durch Vergleich der Ergebnisse konnte so auf geeignete Designs zurückgeschlossen werden. / The treatment of the teeth or diseased of the oral cavity is by far not only administered for aesthetic reasons. These diseases are sometimes also associated with other serious diseases. Studies have shown that some diseases of the mouth, the gingiva or the surrounding area can lead to heart and lung disease or diabetes.
Oftentimes the pathology or symptomatology of dental or oral diseases extends to a wide area. In the dental clinic many different diagnostic devices are used. However, magnetic resonance imaging, which has developed in other areas as the most important diagnostic imaging tool, is not frequently used in dentistry to the present day, although their advantages are obvious: it is known for excellent image contrast, mainly between different soft tissues and comes without hazardous ionizing radiation. The former is probably the reason why the MRI is not yet a standard method in dentistry: here in most cases the contrast of hard tissues is of relevance. However, recent developments and studies demonstrated the versatile advantages of MRI in this area.
The aim of this work as seen from the perspective of application, was to continuously show the enormous potential of MRI in the diverse areas of dentistry. However, many of these applications put very high requirements on the systems. Usually structures to display are very small and require very high resolution. To improve the resolution while using the X-ray method, e.g., one can increase the beam energy (dose). In MRI (without changing the MRT scanner) this results in an extension of measurement time. Especially in the area of the head this oftentimes leads to unwanted movements during the measurement time that worsen the outcome and reproducibility of the obtained diagnostic information or making it completely useless.
The biggest challenge is the measurement of a three-dimensional impression of the tooth surfaces in prosthetics. This process can replace a complex and unpleasant manual impression of the teeth and avoid the production of an intermediate plaster model. Using MRT techniques, a direct three- dimensional model of the tooth surfaces can be produced. By modern CAD/CAM technology, a dental restoration can be directly manufactured by the dental technician using the digital 3D model. Therefore, an important task of the project was the development of a dental MRT method for the detection and correction of movements. Various requirements were imposed on the method. Firstly, the method must be able to detect movements in the range of ~100 µm to fall below the requirements of the final image resolution. For the acquisition of the contrast agent’s signal, a 1-channel receiver coil is used and depending on the measurement, the patient can lie prone or supine. Furthermore, the motion detection system must work without extensive external devices such as cameras, whose direct vision may be obscured by the patient, e.g.
This thesis covers two major subject areas. Firstly, new applications and methods have been developed and further developed in order to provide the various fields of dentistry access to MRT techniques.
Chapter 4 describes the possibility to image the motion of the temporomandibular joint dynamically in real-time. In this work it turned out that both, the movement of the soft tissue components were represented, as well as the intra-articular distance in the TMJ could be measured during mastication (under load) in real-time. Here, the imaging sequence and the corresponding reconstruction algorithm were designed such that the data can be acquired without a prioiri knowledge and processed flexibly. MRT showed different pathologies in the images and dynamic MRT could detect some diseases that could not be diagnosed by other means. The emerging diagnostic possibilities should be investigated and the results verified by large-scale studies.
Chapter 5 describes the results of dento-maxillary MRT imaging, supported by a large-scale study. The diagnostic capabilities of MRI for orthodontic applications are obvious. The typical patient in orthodontics are children and adolescents. The absence of tissue-damaging radiation is a particular advantage of MRI here. After various developments, the acquisition time of a measurement lasted depending on the method only 2 (4) minutes. The resolution in the rendered images was 0.25x0.25x0.50 mm3. Using the proposed method, among other things a geminisation of a tooth root could be shown and the distance of the dental pulp to the tooth surface (enamel) measured.
Chapter 6 presents new developments in the field of digital impressions of tooth surfaces. Here, a new method was developed in order to increase patient comfort during the measurement. This approach helps to prevent movements of the subject in advance. With the old method, the patient lies prone and a large part of the oral cavity is filled with contrast agent. By using a prepared dental cast, the contrast agent can be applied locally and hence the patient may lay supine during the measurement. The associated reproducibility of dental impressions should be shown through a large-scale study.
The main task of this thesis was to develop a method for motion correction that allows to detect the movement of a subject during the measurement without a large number of additional devices and correct the acquired data accordingly. This new navigator method, based on the measurement of a MRT-active marker attached to the subject, makes use of MRT hardware only, except for the additional marker. The method is described in chapter 8. Since this is a new development, it was important to primarily estimate the effects of the various parameters and their impact on the positioning accuracy. This has been evaluated in several preliminary studies, experiments and computer simulations. By validation experiments it was shown in the studies that the image-based navigator detects movements with an accuracy of ~50 µm(translation) and ~0.13◦ (rotation). With the position information obtained from the navigator, the MRT data can be corrected retrospectively or the volume of interest can be adjusted in real-time during the imaging process to prevent inconsistencies in the data in advance. In-vivo MRT data impaired by motion of a subject during the measurement could be corrected using the MoCoLoCo method. By using an appropriate phantom and simulation a movement, it could be shown that using the proposed method, the quality of edge detection (as used in dental impressions, e.g.) could be restored.
Various new high-resolution applications emerged due to the continuous development in hardware, software and algorithms. In chapter 9, the results of a study are presented, which deals with the analysis of shivering movements of the hand during a measurement. For a high-resolution depiction of hand anatomy at 7 T, a suppression of the hand movement is very important. In order to develop an optimal design for a hand receiver coil, a qualitative analysis of the hand movement in several different positions was performed. By comparison of the results, a suitable coil design could be developed.
|
Page generated in 0.0832 seconds