Refine
Has Fulltext
- yes (38)
Is part of the Bibliography
- yes (38)
Year of publication
Document Type
- Doctoral Thesis (38)
Keywords
- Kernspintomografie (23)
- MRI (5)
- Kernspintomographie (4)
- Lunge (4)
- MRT (4)
- Biophysik (3)
- Compressed Sensing (3)
- Magnetresonanztomographie (3)
- Relaxationszeit (3)
- 3D-Bildgebung (2)
Institute
- Physikalisches Institut (33)
- Graduate School of Life Sciences (3)
- Graduate School of Science and Technology (3)
- Institut für Molekulare Infektionsbiologie (1)
- Institut für Theoretische Physik und Astrophysik (1)
- Lehrstuhl für Tissue Engineering und Regenerative Medizin (1)
- Medizinische Klinik und Poliklinik I (1)
Sonstige beteiligte Institutionen
- Universitätsklinikum Würzburg (3)
- Biomedizinische NMR Forschungs GmbH am Max-Planck-Institut fuer biophysikalische Chemie (1)
- Deutsches Zentrum für Herzinsuffizienz (1)
- Fraunhofer Institut für Integrierte Schaltungen (IIS) (1)
- Fraunhofer Institute for Integrierte Schaltungen (IIS) (1)
- Institut für Medizintechnik Schweinfurt (IMeS) (1)
- MRB Research Center for Magnetic-Resonance-Bavaria (1)
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.
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.
Das Ziel dieser Arbeit war die Entwicklung und die Anfertigung eines 3D Erdfeld-NMR Tomographen, um damit die benötigte Technik der MR eines MR-MPI-Tomographen am Lehrstuhl zu etablieren. Daraufhin wurden alle nötigen Komponenten für ein komplettes 3D Erdfeld-NMR-System entwickelt, gebaut und getestet. Mit diesem Wissen wurde in enger Zusammenarbeit mit der MPI-Arbeitsgruppe am Lehrstuhl ein multimodaler MR-MPI-Tomograph angefertigt und die prinzipielle Machbarkeit der technischen Kombination dieser zwei Modalitäten (MRT/MPI) in einer einzigen Apparatur gezeigt.
Auf diesem Entwicklungsweg sind zusätzlich innovative Systemkomponenten entstanden, wie der Bau eines neuen Präpolarisationssystems, mit dem das Präpolarisationsfeld kontrolliert und optimiert abgeschaltet werden kann. Des Weiteren wurde ein neuartiges 3D Gradientensystem entwickelt, das parallel und senkrecht zum Erdmagnetfeld ausgerichtet werden kann, ohne die Bildgebungseigenschaften zu verlieren. Hierfür wurde ein 3D Standard-Gradientensystem mit nur einer weiteren Spule, auf insgesamt vier Gradientenspulen erweitert. Diese wurden entworfen, gefertigt und anhand von Magnetfeldkarten ausgemessen. Anschließend konnten diese Ergebnisse mit der hier präsentierten Theorie und den Simulationsergebnissen übereinstimmend verglichen werden.
MPI (Magnetic Particle Imaging) ist eine neue Bildgebungstechnik mit der nur Kontrastmittel detektiert werden können. Das hat den Vorteil der direkten und eindeutigen Detektion von Kontrastmitteln, jedoch fehlt die Hintergrundinformation der Probe. Wissenschaftliche Arbeiten prognostizieren großes Potential, die Hintergrundinformationen der MRT mit den hochauflösenden Kontrastmittelinformationen mittels MPI zu kombinieren. Jedoch war es bis jetzt nicht möglich, diese beiden Techniken in einer einzigen Apparatur zu etablieren. Mit diesem Prototyp konnte erstmalig eine MR-MPI-Messung ohne Probentransfer durchgeführt und die empfindliche Lokalisation von Kontrastmittel mit der Überlagerung der notwendigen Hintergrundinformation der Probe gezeigt werden. Dies ist ein Meilenstein in der Entwicklung der Kombination von MRT und MPI und bringt die Vision eines zukünftigen, klinischen, multimodalen MR-MPI-Tomographen ein großes Stück näher.
Die MRT hat sich in den letzten Jahren zu einem wichtigen Instrument in der Diagnostik von Herzerkrankungen entwickelt. Da sie ohne ionisierende Strahlung auskommt, stellt sie vor allem auch eine nichtinvasive Alternative zu den nuklearmedizinischen Verfahren und der Computertomographie dar. Im speziellen ermöglicht die kardiale MRT die ortsaufgelöste Darstellung des Herzens mit einer Vielzahl an Kontrasten. Neben der Morphologie können damit auch zahlreiche Funktionsparameter des Herzens, wie die Ejektionsfraktion des linken Ventrikels, oder die Viabilität und Perfusion des Herzmuskels untersucht werden. Atmung und Herzbewegung stellen allerdings große Anforderungen an die MR-Herzbildgebung. Die beiden Störfaktoren limitieren den Zeitraum, der zur Bildakquisition zur Verfügung steht und erzeugen so Konflikte zwischen räumlicher Auflösung, anatomischer Abdeckung, zeitlicher Auflösung und dem Signal-zu-Rausch-Verhältnis (SNR). Ferner ergibt sich für die meisten eingesetzten Verfahren eine erhöhte Komplexität. Die Bildgebungssequenzen müssen mittels EKG an den Herzrhythmus des Patienten angepasst und die Bildakquisitionen im Atemanhaltezustand durchgeführt werden. In manchen Fällen ist sogar eine Aufspaltung der Messung in mehrere Einzelakquisitionen nötig, was wiederum die Dauer der Untersuchungen verlängert und den Patientenkomfort reduziert.
Mit technischen Entwicklungen im Bereich der Gradienten und der Empfangsspulen sowie durch den Einsatz dedizierter Bildgebungstechniken konnten in den letzten Jahren signifikante Verbesserungen erzielt und der Stellenwert der MR-Bildgebung in der Herzdiagnostik erhöht werden. Von großer Bedeutung sind dabei auch Beschleunigungsverfahren wie die Parallele Bildgebung, die eine deutliche Verkürzung der Datenakquisition ermöglichen und so den Einfluss von Atmung und Herzbewegung wirksam reduzieren. Die Beschleunigung wird dabei grundsätzlich durch eine unvollständige Datenakquisition bzw. Unterabtastung des k-Raums erzielt, welche im Zuge der Bildrekonstruktion durch Ausnutzen zusätzlich vorhandener Informationen kompensiert wird. Bei der Parallelen Bildgebung ersetzen beispielsweise mehrere um das Objekt herum angeordnete Empfangsspulen die zum Teil unvollständig durchgeführte Gradientenbasierte Ortskodierung. Die Beschleunigungsverfahren sind allerdings wegen der verringerten Datenaufnahme auch immer mit einer Reduktion des SNR verbunden.
Eine alternative Strategie zur Beschleunigung der 2D-Bildgebung mit mehreren Schichten stellt die simultane Multischichtbildgebung mit Multi-Slice Controlled Aliasing In Parallel Imaging Results In Higher Acceleration(MS-CAIPIRINHA) dar. Anders als bei der konventionellen Parallelen Bildgebung wird die Beschleunigung hier nicht durch eine reduzierte Datenaufnahme erzielt. Vielmehr werden Multiband-RF-Pulse eingesetzt, um die Spins in mehreren Schichten gleichzeitig anzuregen. Durch Anwenden schichtspezifischer RF-Phasenzyklen wird die Phase der Spins individuell in jeder Schicht moduliert, wodurch sich eine gegenseitige Verschiebung der Schichten im FOV ergibt. Die Verschiebung erleichtert die Separation der gleichzeitig angeregten Schichten mit Verfahren der Parallelen Bildgebung. Sie erlaubt außerdem eine Minimierung der bei der Rekonstruktion entstehenden Rauschverstärkung. Die Multischichtbildgebungstechnik zeichnet sich gegenüber der konventionellen Parallelen Bildgebung durch ein wesentlich höheres SNR und durch eine Bildrekonstruktion mit geringeren Rekonstruktionsfehlern aus.
In dieser Dissertation wurden verschiedene Strategien zur Anwendung von MS-CAIPIRINHA in der MRT des Herzens präsentiert sowie ihre Vorund Nachteile gegenübergestellt. Im Allgemeinen ermöglichen die vorgestellten Konzepte eine hinsichtlich des SNR sehr effiziente Erweiterung der
anatomischen Abdeckung. Unter anderem wurde eine Möglichkeit vorgestellt, mit der es uneingeschränkt gelingt, MS-CAIPIRINHA in der Bildgebung mit bSSFP-Sequenzen anzuwenden. Die Steady-State-Sequenz wird aufgrund ihres hohen intrinsischen SNR und vorteilhaften Kontrastverhaltens sehr häufig in der MRT des Herzens bei 1,5T eingesetzt. Wie auch die simultane Multischichtbildgebung erfordert sie zum Halten der Magnetisierung im stationären Zustand die Applikation eines dedizierten RF-Phasenzyklus während der Datenakquisition. Der Phasenzyklus der Sequenz ist allerdings nicht ohne Weiteres mit den Phasenzyklen der Multischichttechnik kompatibel, so dass eine Verknüpfung der beiden Verfahren bisher nur durch Aufspalten der Bildakquisition in mehrere Teilmessungen gelang. Mit dem in Kapitel 5 vorgestellten Konzept ist diese zumeist impraktikable Segmentierung nicht mehr erforderlich. Generalisierte RF-Phasenzyklen, die sowohl die Anforderungen der Sequenz, als auch die der Multischichtbildgebung erfüllen, ermöglichen eine uneingeschränkte Anwendung der Multischichttechnik in der Bildgebung mit bSSFP oder vergleichbaren Steady-State-Sequenzen. Die Multischichttechnik ist damit auch bei Untersuchungen in Echtzeit oder mit Magnetisierungspräparation – Verfahren, die unter anderem in der MR-Herzdiagnostik Verwendung finden – einsetzbar. Anhand von Echtzeit-, Cine- und First-Pass-Herzperfusionsuntersuchungen am menschlichen Herzen konnte die Anwendbarkeit des Konzepts erfolgreich demonstriert werden. Durch die Akquisition zweier Schichten in der Zeit, die normalerweise zur Bildgebung einer einzelnen Schicht benötigt wird, gelang eine Verdoppelung der anatomischen Abdeckung bei unverändert hoher Bildqualität. Bei den Herzperfusionsuntersuchungen konnten je RR-Intervall sechs Schichten akquiriert werden. Bei Echtzeit- und Cine-Messungen erlaubt das Konzept eine signifikante Reduktion der Anzahl der Atemanhaltezustände und dementsprechend eine wirksame Verkürzung der Patientenuntersuchung und eine Verbesserung des Patientenkomforts.
In Kapitel 6 wurde eine effiziente Strategie zur Anwendung der simultanen Multischichtbildgebung in der First-Pass-Herzperfusionsbildgebung bei 3T vorgestellt. Es wurde gezeigt, dass durch den Einsatz von MS-CAIPIRINHA mit Beschleunigungsfaktoren, die größer sind als die Anzahl der simultan angeregten Schichten, neben der anatomischen Abdeckung auch die räumliche Auflösung innerhalb der Bildgebungsschicht erhöht werden kann. Beide Verbesserungen sind für die MR-gestützte Diagnostik der Koronaren Herzerkrankung von Bedeutung. Während mit einer hohen räumlichen Auflösung subendokardiale und transmurale Infarktareale unterschieden werden können, erleichtert eine hohe anatomische Abdeckung die genaue Eingrenzung hypoperfundierter Bereiche. Das grundsätzliche Prinzip der vorgestellten Strategie besteht in der Kombination zweier unterschiedlicher Beschleunigungsansätze: Zur Verbesserung der anatomischen Abdeckung kommt die simultane Multischichtbildgebung zum Einsatz. Zusätzlich zur gleichzeitigen Anregung mehrerer Schichten wird der k-Raum regelmäßig unterabgetastet. Die dabei erzielte Beschleunigung wird zur Verbesserung der räumlichen Auflösung eingesetzt. Die Bildrekonstruktion erfolgt mit Verfahren der Parallelen Bildgebung. Der Vorteil des Konzepts liegt insbesondere im vollständigen Erhalt der Datenakquisitionszeit gegenüber einer unbeschleunigten Messung mit Standardabdeckung und -auflösung. Anders als bei konventionellen Beschleunigungsverfahren wirken sich lediglich die Verkleinerung der Voxelgröße sowie die Rauschverstärkung der Bildrekonstruktion SNR-reduzierend aus. Die Rauschverstärkung wird dabei, durch die gegenseitige Verschiebung der simultan angeregten Schichten im FOV, so gering wie möglich gehalten. Die Anwendbarkeit des Konzepts konnte anhand von Simulationen sowie Untersuchungen an Probanden und Herzinfarktpatienten erfolgreich demonstriert werden. Simultanes Anregen zweier Schichten und 2,5-faches Unterabtasten des k-Raums ermöglichte die Durchführung von Untersuchungen mit einer anatomischen Abdeckung von sechs bis acht Schichten je RR-Intervall und einer räumlichen Auflösung von 2,0×2,0×8,0mm3. Es konnte gezeigt werden, dass die angewandte GRAPPA-Rekonstruktion, trotz der effektiv fünffachen Beschleunigung, robust und im Wesentlichen mit geringer Rauschverstärkung durchführbar ist. Bildqualität und SNR waren für eine sektorweise Absolutquantifizierung der Myokardperfusion ausreichend, während die hohe räumliche Auflösung die Abgrenzung kleiner subendokardialer Perfusionsdefizite ermöglichte. Aufgrund seiner großen Flexibilität und recht einfachen Implementierbarkeit ist das Beschleunigungskonzept vielversprechend hinsichtlich einer Anwendung in der klinischen Routine. Die diesbezügliche Tauglichkeit ist allerdings in weiterführenden Patientenstudien noch zu evaluieren.
Alternativ zu diesem Konzept wurde in Kapitel 7 noch eine weitere, ebenfalls auf MS-CAIPIRINHA basierende Strategie für die First-Pass-Herzperfusionsbildgebung bei 3T mit großer anatomischer Abdeckung und hoher räumlicher Auflösung vorgestellt. Wie zuvor bestand die Grundidee des Konzepts darin, MS-CAIPIRINHA mit Beschleunigungsfaktoren anzuwenden, welche größer sind als die Anzahl der simultan angeregten Schichten und die Vergrößerung der anatomischen Abdeckung durch simultanes Anregen mehrerer Schichten zu realisieren. Um allerdings die bei der Bildrekonstruktion und Schichtseparation entstehende Rauschverstärkung zu minimieren, wurde zur Verbesserung der räumlichen Auflösung innerhalb der Schicht das nichtlineare Beschleunigungsverfahren Compressed Sensing zum Einsatz gebracht. Die erst in den letzten Jahren entwickelte Technik ermöglicht die exakte Rekonstruktion zufällig unterabgetasteter Daten, sofern bekannt ist, dass sich das rekonstruierte Bild in eine wohldefinierte sparse Darstellung überführen lässt. Neben der Erreichbarkeit hoher Beschleunigungsfaktoren bietet Compressed Sensing den Vorteil einer Bildrekonstruktion ohne signifikante Rauscherhöhung. Zur Einbindung des Verfahrens in das Multischichtbildgebungskonzept erfolgt die für die Verbesserung der Auflösung nötige Unterabtastung des k-Raums, zufällig und inkohärent. Zur Bildrekonstruktion sind zwei Teilschritte erforderlich. Im ersten Teilschritt werden die durch die zufällige Unterabtastung entstandenen inkohärenten Artefakte mit Compressed Sensing entfernt, im zweiten die gleichzeitig angeregten Schichten mit Verfahren der Parallelen MRT separiert. Es konnte gezeigt werden, dass die Kombination aus Compressed Sensing und MS-CAIPIRINHA eine Reduktion der inhomogenen Rauschverstärkung ermöglicht und zur Durchführung von qualitativen First-Pass-Herzperfusionsuntersuchungen mit einer Abdeckung von sechs bis acht Schichten je RR-Intervall sowie einer räumlichen Auflösung von 2,0 × 2,0 × 8,0mm3 geeignet ist. Des Weiteren konnte gezeigt werden, dass das angewandte Multischicht-Bildgebungskonzept einer Anwendung des entsprechenden Compressed-Sensing-Konzepts ohne simultane Multischichtanregung überlegen ist. Es stellte sich allerdings auch heraus, dass die rekonstruierten Bilder mit systematischen Fehlern behaftet sind, zu welchen auch ein signifikanter rekonstruktionsbedingter Verlust an zeitlicher Auflösung zählt. Dieser kann zu einer Verzerrung quantitativ bestimmter Perfusionswerte führen und verhindert so robuste quantitative Messungen der Myokardperfusion. Es ist außerdem davon auszugehen, dass auch abrupte Signalveränderungen, die bei Arrhythmien oder Bewegung auftreten, nur sehr ungenau rekonstruiert werden können. Die Systematischen Rekonstruktionsfehler konnten anhand zweier Verfahren, einer Monte-Carlo-Simulation sowie einer Analyse der lokalen Punktantworten präzise Untersucht werden. Die beiden Analysemethoden ermöglichten einerseits die genaue Bestimmung systematischer und statistischer Abweichungen der Signalamplitude und andererseits die Quantifizierung rekonstruktionsbedingter zeitlicher und räumlicher Auflösungsverluste. Dabei konnte ein Mangel an Sparsität als grundlegende Ursache der Rekonstruktionsfehler ermittelt werden. Die bei der Analyse eingesetzten Verfahren erleichtern das Verständnis von Compressed Sensing und können beispielsweise bei der Entwicklung nichtlinearer Beschleunigungskonzepte zur Bildqualitätsanalyse eingesetzt werden.
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.
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.
Verschiedene Konzepte der Röntgenmikroskopie haben sich mittlerweile im Labor etabliert und ermöglichen heute aufschlussreiche Einblicke in eine Vielzahl von Probensystemen. Der „Labormaßstab“ bezieht sich dabei auf Analysemethoden, die in Form von einem eigenständigen Gerät betrieben werden können. Insbesondere sind sie unabhängig von der Strahlerzeugung an einer Synchrotron-Großforschungseinrichtung und einem sonst kilometergroßen Elektronen-speicherring. Viele der technischen Innovationen im Labor sind dabei ein Transfer der am Synchrotron entwickelten Techniken. Andere wiederum basieren auf der konsequenten Weiterentwicklung etablierter Konzepte. Die Auflösung allein ist dabei nicht entscheidend für die spezifische Eignung eines Mikroskopiesystems im Ganzen. Ebenfalls sollte das zur Abbildung eingesetzte Energiespektrum auf das Probensystem abgestimmt sein. Zudem muss eine Tomographieanalage zusätzlich in der Lage sein, die Abbildungsleistung bei 3D-Aufnahmen zu konservieren.
Nach einem Überblick über verschiedene Techniken der Röntgenmikroskopie konzentriert sich die vorliegende Arbeit auf quellbasierte Nano-CT in Projektionsvergrößerung als vielversprechende Technologie zur Materialanalyse. Hier können höhere Photonenenergien als bei konkurrierenden Ansätzen genutzt werden, wie sie von stärker absorbierenden Proben, z. B. mit einem hohen Anteil von Metallen, zur Untersuchung benötigt werden. Das bei einem ansonsten idealen CT-Gerät auflösungs- und leistungsbegrenzende Bauteil ist die verwendete Röntgen-quelle. Durch konstruktive Innovationen sind hier die größten Leistungssprünge zu erwarten. In diesem Zuge wird erörtert, ob die Brillanz ein geeignetes Maß ist, um die Leistungsfähigkeit von Röntgenquellen zu evaluieren, welchen Schwierigkeiten die praktische Messung unterliegt und wie das die Vergleichbarkeit der Werte beeinflusst. Anhand von Monte-Carlo-Simulationen wird gezeigt, wie die Brillanz verschiedener Konstruktionen an Röntgenquellen theoretisch bestimmt und miteinander verglichen werden kann. Dies wird am Beispiel von drei modernen Konzepten von Röntgenquellen demonstriert, welche zur Mikroskopie eingesetzt werden können. Im Weiteren beschäftigt sich diese Arbeit mit den Grenzen der Leistungsfähigkeit von Transmissionsröntgenquellen. Anhand der verzahnten Simulation einer Nanofokus-Röntgenquelle auf Basis von Monte-Carlo und FEM-Methoden wird untersucht, ob etablierte Literatur¬modelle auf die modernen Quell-konstruktionen noch anwendbar sind. Aus den Simulationen wird dann ein neuer Weg abgeleitet, wie die Leistungsgrenzen für Nanofokus-Röntgenquellen bestimmt werden können und welchen Vorteil moderne strukturierte Targets dabei bieten.
Schließlich wird die Konstruktion eines neuen Nano-CT-Gerätes im Labor-maßstab auf Basis der zuvor theoretisch besprochenen Nanofokus-Röntgenquelle und Projektionsvergrößerung gezeigt, sowie auf ihre Leistungsfähigkeit validiert. Es ist spezifisch darauf konzipiert, hochauflösende Messungen an Materialsystemen in 3D zu ermöglichen, welche mit bisherigen Methoden limitiert durch mangelnde Auflösung oder Energie nicht umsetzbar waren. Daher wird die praktische Leistung des Gerätes an realen Proben und Fragestellungen aus der Material¬wissenschaft und Halbleiterprüfung validiert. Speziell die gezeigten Messungen von Fehlern in Mikrochips aus dem Automobilbereich waren in dieser Art zuvor nicht möglich.
X-ray dark-field imaging allows to resolve the conflict between the demand for centimeter scaled fields of view and the spatial resolution required for the characterization of fibrous materials structured on the micrometer scale. It draws on the ability of X-ray Talbot interferometers to provide full field images of a sample's ultra small angle scattering properties, bridging a gap of multiple orders of magnitude between the imaging resolution and the contrasted structure scale. The correspondence between shape anisotropy and oriented scattering thereby allows to infer orientations within a sample's microstructure below the imaging resolution. First demonstrations have shown the general feasibility of doing so in a tomographic fashion, based on various heuristic signal models and reconstruction approaches. Here, both a verified model of the signal anisotropy and a reconstruction technique practicable for general imaging geometries and large tensor valued volumes is developed based on in-depth reviews of dark-field imaging and tomographic reconstruction techniques.
To this end, a wide interdisciplinary field of imaging and reconstruction methodologies is revisited. To begin with, a novel introduction to the mathematical description of perspective projections provides essential insights into the relations between the tangible real space properties of cone beam imaging geometries and their technically relevant description in terms of homogeneous coordinates and projection matrices. Based on these fundamentals, a novel auto-calibration approach is developed, facilitating the practical determination of perspective imaging geometries with minimal experimental constraints. A corresponding generalized formulation of the widely employed Feldkamp algorithm is given, allowing fast and flexible volume reconstructions from arbitrary tomographic imaging geometries. Iterative reconstruction techniques are likewise introduced for general projection geometries, with a particular focus on the efficient evaluation of the forward problem associated with tomographic imaging. A highly performant 3D generalization of Joseph's classic linearly interpolating ray casting algorithm is developed to this end and compared to typical alternatives. With regard to the anisotropic imaging modality required for tensor tomography, X-ray dark-field contrast is extensively reviewed. Previous literature is brought into a joint context and nomenclature and supplemented by original work completing a consistent picture of the theory of dark-field origination. Key results are explicitly validated by experimental data with a special focus on tomography as well as the properties of anisotropic fibrous scatterers. In order to address the pronounced susceptibility of interferometric images to subtle mechanical imprecisions, an efficient optimization based evaluation strategy for the raw data provided by Talbot interferometers is developed. Finally, the fitness of linear tensor models with respect to the derived anisotropy properties of dark-field contrast is evaluated, and an iterative scheme for the reconstruction of tensor valued volumes from projection images is proposed. The derived methods are efficiently implemented and applied to fiber reinforced plastic samples, imaged at the ID19 imaging beamline of the European Synchrotron Radiation Facility. The results represent unprecedented demonstrations of X-ray dark-field tensor tomography at a field of view of 3-4cm, revealing local fiber orientations of both complex shaped and low-contrast samples at a spatial resolution of 0.1mm in 3D. The results are confirmed by an independent micro CT based fiber analysis.
Die vorliegende Arbeit beschäftigt sich mit der volumenbildgebenden Untersuchung von mechanischen Ermüdungsprozessen in Titan-Zahnimplantaten. Im Vordergrund steht die Entwicklung einer neuen Messmethode der In-situ-Mikrotomografie am Synchrotron.
Zahnimplantate werden beim Gebrauch mechanisch wiederholt belastet (Wechsellast). Nach vielen zyklischen Belastungen können aufgrund von mikroplastische Verformungen Ermüdungsschäden auftreten. Diese können im Extremfall zum Versagen und Verlust eines Implantats führen. Die Computertomographie ist eine sehr geeignete zerstörungsfrei Prüfmethode, um Zahnimplantate zu untersuchen. Diese Arbeit erweitert die bisherige CT-Methode insofern, dass In-situ-Beobachtungen bei mechanischer Belastung möglich sind.
Die in dieser Arbeit untersuchten Zahnimplantate weisen an der Implantat-Abutment-Grenzfläche bei eintretender Ermüdung einen Mikrospalt auf. Dieser wird als Indikator für einsetzende Fatigue- Prozesse benutzt. Der in der Synchrotron CT verfügbare Inlinephasenkontrast ermöglicht eine verbesserte Bestimmung der Mikrospaltgröße. Da die schnellen Bewegungen der Ermüdungsprüfung mittels Standard-CT-Verfahren schwer zu erfassen sind, war die stroboskopische Aufnahmemethode das zielführende Messverfahren, um in-situ-Prüfung zu ermöglichen.
Die 4 kommerziellen Zahnimplantattypen werden neben der In-situ-Fatigue Prüfung auch mittels klassischer Ermüdungsprüfung untersucht und mit der Neuen Messmethode verglichen. Die hier entwickelte In-situ-Fatigue-Prüfstation kann Proben bis zu 345 N tomographisch untersuchen. Neben den experimentellen Untersuchungen wird eine statische FEM-Betrachtung durchgeführt und mit experimentellen Messdaten verglichen. Zuletzt wird mit der entwickelten Messtation Knochenrisse in der Implantat Umgebung untersucht.
Quantifizierung myokardialer Mikrostruktur und Perfusion mittels longitudinaler NMR Relaxation
(2018)
Ziel der Arbeit war es die Quantifizierung funktioneller bzw. mikrostruktureller Parameter des Herzmuskels mit Hilfe T1-basierter Methoden zu verbessern. Diese Methoden basieren darauf, die gewünschte Information durch eine geeignete Präparation der Magnetisierung bzw. durch die Gabe von Kontrastmittel in den Zeitverlauf der longitudinalen Relaxation zu kodieren. Aus der Änderung der Relaxationszeit läßt sich dann die gewünschte Information bestimmen. Dafür sollte sowohl der Einfluß der Anatomie als auch derjenige der Meßmethodik auf die Bestimmung der longitudinalen Relaxationszeit und damit auf die Quantifizierung der Funktion bzw. Mikrostrukturparameter untersucht werden.
Speziell der Einfluß der Bildgebungssequenz führt dazu, daß nur eine scheinbare Relaxationszeit gemessen wird. Während dies keinen Einfluß auf die T1-basierte Bestimmung der untersuchten Mikrostrukturparameter hatte, ergab sich für die Perfusionsquantifizierung eine deutliche Abhängigkeit von den Parametern der verwendeten IRLL-Sequenz. Um diesen Einfluß gerecht zu werden, wurden an die Meßmethodik angepaßte Gleichungen zur Bestimmung der Perfusion gefunden mit denen die systematischen Abweichungen korrigiert werden können. Zusätzlich reduzieren die angepaßten Gleichungen die Anforderungen bezüglich der Inversionsqualität im schichtselektiven Experiment. Dies wurde in einem weiteren Projekt bei der Bestimmung der Nierenperfusion im Mausmodell ausgenutzt.
Neben der Untersuchung der Auswirkungen der Meßmethode wurde auch der Einfluß der anatomischen Besonderheiten des Blutkreislaufs am Herzen auf die Parameterquantifizierung mittels T1-basierter Methoden untersucht. Es konnte gezeigt werden, daß auf Grund der Anatomie des Herzens bei typischen Orientierungen der Bildgebungsschicht, auch bei der schichtselektiven Inversionspräparation der Magnetisierung des Herzmuskels ein Anteil des Blutpools invertiert wird. Daraus folgt, daß die vereinfachende Annahme, nach welcher bei schichtselektiver Präparation in Folge von Perfusion nur Blut mit Gleichgewichtsmagnetisierung den Herzmuskel erreicht, nicht erfüllt ist. Es konnte gezeigt werden, daß dies bei Perfusion zu einer deutlichen Unterschätzung der berechneten Perfusionswertes führt. Um mit diesem Problem umgehen zu können, wurde aufbauend auf einem vereinfachten Modell der zeitlichen Entwicklung der Blutmagnetisierung eine Korrektur für die Bestimmung der Perfusionswerte gefunden welche den Einfluß der anatomischen Besonderheiten berücksichtigt.
Das für die Perfusionskorrektur eingeführte Model prognostiziert ebenso, daß auch bei schichtselektiver Inversion die T1-basierte Bestimmung der untersuchten Mikrostrukturparameter von der Perfusion abhängig wird und eine systematische Überschätzung der quantifizierten Werte verursacht. Da die Perfusion im Kleintier deutlich höher ist als im Menschen, ist dieser Einfluß besonders in der präklinischen Forschung zu beachten. So können dort allein durch verminderte Perfusion deutliche Änderungen in den bestimmten Werten der Mikrostrukturparameter erzeugt werden, welche zu einer fehlerhaften Interpretation der Ergebnisse führen und somit ein falsches Bild für die Vorgänge im Herzmuskel suggerieren. Dabei bestätigt der Vergleich mit experimentellen Ergebnissen aus der Literatur die Vorhersagen für das Rattenmodell. Beim Menschen ist der prognostizierte Effekt deutlich kleiner. Der prognostizierte Fehler bspw. im RBV-Wert liegt in diesem Fall bei etwa 10% und wird üblicherweise in der aktuellen Forschung vernachlässigt. Inwieweit dies in er klinischen Forschung gerechtfertigt ist, muß in weiteren Untersuchungen geklärt werden.
Den untersuchten Methoden zur Bestimmung von funktionellen und mikrostrukturellen Parametern ist gemein, daß sie eine exakte Quantifizierung der longitudinalen Relaxationszeit T1 benötigen. Dabei ist im Kleintierbereich die klassische IRLL-Methode als zuverlässige Sequenz zur T1-Quantifizierung etabliert. In der klinischen Bildgebung werden auf Grund der unterschiedlichen Zeitskalen und anderer technischer Voraussetzungen andere Anforderungen an die Datenakquisition gestellt. Dabei hat in den letzten Jahren die MOLLI-Sequenz große Verbreitung gefunden. Sie ist eine Abwandlung der IRLL-Sequenz, bei der mit einer bSSFP-Bildgebungssequenz getriggert ganze Bilder während eines Herzschlages aufgenommen werden. Die MOLLI-Sequenz reagiert dabei empfindlich auf die Wartezeiten zwischen den einzelnen Transienten. Um mit diese Problematik in den Griff zu bekommen und gleichzeitig die Meßzeit verkürzen zu können wurde eine neue Methode zum Fitten der Daten entwickelt, welche die Abhängigkeit der scheinbaren Relaxationszeit von der Wartezeit zwischen den einzelnen Transienten, sowie der mittleren Herzrate fast vollständig eliminiert. Diese Methode liefert für das ganze klinisch Spektrum an erwarteten T1-Zeiten, vor und nach Kontrastmittelgabe, stabile Ergebnisse und erlaubte ein deutliche Verkürzung der Meßzeit, ohne die Anzahl der aufgenommenen Meßzeitpunkte zu reduzieren. Dies wurde in einer initialen klinischen Studie genutzt, um ECV-Werte in Patienten zu bestimmen.
Ein Nachteil der Verwendung der MOLLI-Sequenz ist, daß nur die scheinbare Relaxationszeit aus den Fit der Meßdaten bestimmt wird. Die standardmäßig genutzte Korrektur benutzt aber dem gefitteten Wert der Gleichgewichtsmagnetisierung um den wahren T1-Wert zu bestimmen. Somit ist es für die Bestimmung des T1-Wertes notwendig, die Qualität der Inversionspräparation zu kennen. Auf Basis der neuen Fitmethode wurde eine Anpassung der MOLLI-Sequenz demonstriert, welche die Bestimmung der Gleichgewichtsmagnetisierung unabhängig von der Qualität der Inversionspräparation erlaubt. Dafür verlängert sich die Meßdauer lediglich um einen Herzschlag um in geeigneter Weise ein zusätzliches Bild aufnehmen zu können.
Abschließend wurde in dieser Arbeit der Signal-Zeit-Verlauf der MOLLI-Sequenz eingehend theoretische untersucht um ein besseres Verständnis der getriggerten IRLL-Sequenzen zu entwickeln. In diesem Zusammenhang konnte eine einfache Interpretation der scheinbaren Relaxationszeit gefunden werden. Ebenso konnte erklärt werden, warum die für ungetriggerte IRLL-Sequenzen abgeleitete Korrekturgleichung auch im getriggerten Fall erstaunlich gute Ergebnisse liefert. Weiterhin konnten Fehlerquellen für die verbleibenden Abweichungen identifiziert werden, welche als Ausgangspunkt für die Ableitung verbesserter Korrekturgleichungen genutzt werden können.
Ziel dieser Arbeit ist es, die quantitative MRT in den Fokus zu rücken. In den letzten Jahren hat sich auf diesem Forschungsgebiet viel weiterentwickelt und es wurden verschiedenste Sequenzen und Methoden vorgestellt, um insbesondere Relaxationszeitparameter quantitativ in kurzer Zeit zu messen. Steady-State-Sequenzen eignen sich besonders für diese Thematik, da sie kurze Messzeiten benötigen und darüber hinaus ein relativ hohes SNR besitzen. Speziell die IR TrueFISP-Sequenz bietet für die Parameterquantifizierung viel Potential. Ursprünglich wurde diese Sequenz an der Universität Würzburg zur simultanen Messung von T1- und T2-Relaxationszeiten vorgestellt und hinsichtlich der Zeiteffizienz weiterentwickelt. In dieser Arbeit wurde ein neuartiger iterativer Rekonstruktionsansatz für die IR TrueFISP-Sequenz entwickelt, der auf einer Hauptkomponentenanalyse (PCA) basiert und sich die glatten Signalverläufe zu Nutze macht. Aufgrund der hohen Zeitauflösung dieser Rekonstruktionstechnik werden dabei auch Gewebekomponenten mit kurzen Relaxationszeiten detektierbar. Weiterhin bewahrt der Rekonstruktionsansatz Informationen mehrerer Gewebekomponenten innerhalb eines Voxels und ermöglicht damit eine relaxographische Untersuchung. Insbesondere beim Menschen führen der Partialvolumeneffekt und die Mikrostruktur des Gewebes zu Signalverläufen, die ein multi-exponentielles Signal liefern. Die MR-Relaxographie, also die Darstellung von Relaxationszeitverteilungen innerhalb eines Voxels, stellt eine Möglichkeit dar, um die beteiligten Gewebekomponenten aus dem überlagerten Signalverlauf zu extrahieren. Insgesamt bilden die optimierte Relaxometrie mit der Möglichkeit der analytischen Korrektur von Magnetfeldinhomogenitäten und die beschleunigte Relaxographie die Hauptteile dieser Dissertation. Die Hauptkapitel werden im Folgenden noch einmal gesondert zusammengefasst.
Die simultane Aufnahme der quantitativen T1- und T2-Parameter-Karten kann mit einem
Goldenen-Winkel-basiertem radialen IR TrueFISP-Readout in ungefähr 7 Sekunden pro
Schicht erreicht werden. Die bisherige Rekonstruktionstechnik mit dem KWIC-Filter ist
durch dessen breite Filter-Bandbreite und somit in der zeitlichen Auflösung limitiert. Besonders bei hohen räumlichen Frequenzen wird eine sehr große Anzahl an Projektionen
zusammengefasst um ein Bild zu generieren. Dies sorgt dafür, dass Gewebekomponenten mit kurzer T1*-Relaxationszeit (z.B. Fett oder Myelin) nicht akkurat aufgelöst werden können. Um dieses Problem zu umgehen, wurde die T1* shuffling-Rekonstruktion entwickelt, die auf dem T2 Shuffling-Ansatz basiert. Diese Rekonstruktionstechnik macht sich die glatten Signalverläufe der IR TrueFISP-Sequenz zu Nutze und ermöglicht die Anwendung einer PCA. Die iterative Rekonstruktion sorgt dafür, dass mit nur acht kombinierten Projektionen pro generiertem Bild eine merklich verbesserte temporäre Auflösung erzielt werden kann. Ein Nachteil ist jedoch das stärkere Rauschen in den ersten Bildern der Zeitserie bedingt durch die angewandte PCA. Dieses verstärkte Rauschen äußert sich in den leicht erhöhten Standardabweichungen in den berechneten Parameter-Karten. Jedoch ist der Mittelwert näher an den Referenzwerten im Vergleich zu den Ergebnissen mit dem KWIC-Filter. Letztendlich kann man sagen, dass die Ergebnisse leicht verrauschter, aber exakter sind.
Mittels zusätzlichen Regularisierungstechniken oder Vorwissen bezüglich des Rauschlevels
wäre es zudem noch möglich, das SNR der ersten Bilder zu verbessern, um dadurch den
beschriebenen Effekt zu verringern.
Grundsätzlich hängt die Genauigkeit von IR TrueFISP vom T1/T2-Verhältnis des betreffenden
Gewebes und dem gewählten Flipwinkel ab. In dieser Arbeit wurde der Flipwinkel besonders für weiße und graue Masse im menschlichen Gehirn optimiert. Mit den verwendeten 35° wurde er außerdem etwas kleiner gewählt, um zudem Magnetisierungstransfereffekte zu minimieren. Mit diesen Einstellungen ist die Präzision vor allem für hohe T1- und niedrige T2-Werte sehr gut, wird jedoch insbesondere für höhere T2-Werte schlechter. Dies ist aber ein generelles Problem der IR TrueFISP-Sequenz und hängt nicht mit der entwickelten Rekonstruktionsmethode zusammen. Außerdem wurde im fünften Kapitel eine
Akquisitionstechnik vorgestellt, die eine 3D-Abdeckung der quantitativen Messungen des Gehirns in klinisch akzeptabler Zeit von unter 10 Minuten erzielt. Dies wird durch Einsatz der parallelen Bildgebung erreicht, da eine Kombination aus radialer Abtastung in der Schicht und kartesischer Aufnahme in Schichtrichtung (Stack-of-Stars) vorliegt.
Ein großes Problem in der Steady-State-Sequenz (und somit auch bei IR TrueFISP) sind
Magnetfeldinhomogenitäten, die durch Suszeptibilitätsunterschiede verschiedener Gewebe und/oder Inhomogenitäten des Hauptmagnetfeldes hervorgerufen werden. Diese führen zu Signalauslöschungen und damit verbunden zu den beschriebenen Banding-Artefakten. Mithilfe der analytisch ermittelten Korrekturformeln ist es nun möglich, die berechneten (T1,T2)-Wertepaare unter Berücksichtigung der tatsächlich auftretenden Off- Resonanzfrequenz für einen großen Bereich zu korrigieren. An den kritischen Stellen, an denen die Bandings auftreten, liefert jedoch auch diese Korrektur keine brauchbaren Ergebnisse. Grundsätzlich ist es für die Genauigkeit der Ergebnisse stets zu empfehlen, die Flipwinkel- und B0-Karte zusätzlich mit aufzunehmen, um diese Parameter für die quantitative Auswertung exakt zu kennen. Mit den beschriebenen Methoden aus Kapitel 6 könnte es prinzipiell auch möglich sein, die Off-Resonanzfrequenz aus dem Signalverlauf zu ermitteln und auf die zusätzliche Messung der B0-Karte zu verzichten. B0-Änderungen während der Messung, die von der Erwärmung der passiven Shim-Elemente im MR-System hervorgerufen werden, sind kaum zu korrigieren. Ein stabiler Scanner ohne B0-Drift ist deshalb für quantitative Auswertungen erforderlich.
Die erwähnte Messzeit von 7 Sekunden pro Schicht garantiert, dass auch Gewebe mit längeren Relaxationskomponenten annähernd im Steady-State sind, was wiederum für das Umkehren des Signals in den abklingenden Verlauf gegen Null und die anschließende
Multikomponentenanalyse (vgl. Kapitel 7) notwendig ist. Mit der inversen Laplace-
Transformation ist es innerhalb eines Voxels möglich, Signalverläufe auf mehrere Komponenten hin zu untersuchen. Der ursprünglich angenommene mono-exponentielle Verlauf wird durch ein multi-exponentielles Verhalten abgelöst, was vor allem in biologischem Gewebe eher der Wahrheit entspricht. Gewebe mit kurzen Relaxationskomponenten (T1* < 200 ms) sind klinisch relevant und mit T1* shuffling detektierbar. Vor allem Myelin innerhalb des Gehirns ist bei neurologischen Fragestellungen ein Indikator zur Diagnose im Frühstadium (z.B. für neurodegenerative Erkrankungen) und
deshalb von besonderem Interesse. Die Integration über verschiedene T1*-Zeitbereiche im
T1*-Spektrum ermöglicht dazu die Erstellung von Gewebekomponenten-Karten, mithilfe
derer klinische Auswertungen sinnvoll wären. Die Erstellung dieser Karten ist prinzipiell
möglich und funktioniert für mittlere und lange Gewebekomponenten recht gut. Die
klinisch relevanten kurzen Gewebekomponenten sind dagegen bei der radialen Aufnahme
mit nur einem Schuss noch nicht befriedigend. Deshalb wurde die Aufnahmetechnik in
eine quasi-zufällige kartesische Akquisition mit mehreren Schüssen weiterentwickelt. Die Ergebnisse wurden in Kapitel 7 vorgestellt und sind vielversprechend. Einzig die Messzeit sollte mit zusätzlichen Beschleunigungen noch weiter verkürzt und auf eine kartesische 3D-Akquisition erweitert werden.
Die Beschränkung auf T1*-Spektren bei der Multikomponentenanalyse und die Tatsache, dass deren Amplitude von einer Kombination von S0 und Sstst abhängen, führen dazu, dass es nicht ohne Weiteres möglich ist für einen einzelnen Gewebetyp an die T1- und T2-Information zu gelangen. In Kapitel 8 wurde gezeigt, dass dies mit einer zusätzlichen
Messung gelingen kann. Das finale Ergebnis dieser Messungen ohne und mit Inversion sind zweidimensionale Spektren, bei der für jede Gewebekomponente innerhalb eines Voxels der T1- und T2-Wert abgelesen werden kann. Wichtig hierbei ist die Tatsache, dass der verwendete Ansatz kein Vorwissen über die Anzahl der zu erwartenden Gewebekomponenten (Peaks) im Voxel voraussetzt. Auch bei dieser Methodik ist die Kenntnis über den tatsächlichen Flipwinkel von Bedeutung, da dieser in den Formeln zur Berechnung von T1 und T2 verwendet wird. Die Stabilität des B0-Feldes ist hier ebenso von enormer Bedeutung, da Änderungen zwischen den beiden Messungen zu einem unterschiedlichen Steady-State und somit zu Abweichungen bei den nachfolgenden Berechnungen führen, die auf den selben Steady-State-Wert ausgelegt sind.
Zusammenfassend lässt sich sagen, dass mit dieser Arbeit die Grundlagen für genauere
und robustere quantitative Messungen mittels Steady-State-Sequenzen gelegt wurden. Es
wurde gezeigt, dass sich Relaxationszeitspektren für jedes einzelne Voxel generieren lassen.
Dadurch ist eine verbesserte Auswertung möglich, um genauere Aussagen über die Zusammensetzung einer Probe (vor allem beim menschlichen Gewebe) treffen zu können. Zudem wurde die Theorie für ultraschnelle 2D-Relaxographie-Messungen vorgestellt. Erste”Proof of Principle“-Experimente zeigen, dass es möglich ist, 2D-Relaxationszeitspektren in sehr kurzer Zeit zu messen und graphisch darzustellen. Diese Aufnahme- und Datenverarbeitungstechnik ist in dieser Form einmalig und in der Literatur kann bis dato keine schnellere Methode gefunden werden.
Neue Therapieansätze durch Tissue Engineering erfordern gleichzeitig angepasste Diagnosemöglichkeiten und nicht-invasive Erfolgskontrollen. Speziell die 3D-MR-Bildgebung ist ein vielversprechendes Instrument, um Parameter mit hoher räumlicher Präzision zu quantifizieren. Vor diesem Hintergrund wurden im Rahmen dieser Arbeit neue Ansätze für die hochauflösende 3D-MRT in vivo entwickelt und deren Eignung im Bereich des Tissue Engineerings gezeigt.
Welchen Vorteil die Quantifizierung von Parametern bietet, konnte im Rahmen einer prä-klinischen Studie an einem Modell der Hüftkopfnekrose gezeigt werden. Der Therapieverlauf wurde zu verschiedenen Zeitpunkten kontrolliert. Trotz der niedrigen räumlichen Auflösung, konnten durch eine systematische Auswertung der Signalintensitäten von T1- und T2-FS-gewichteten Aufnahmen Rückschlüsse über Veränderungen in der Mikrostruktur gezogen werden, die darüber hinaus in guter Übereinstimmung mit Ergebnissen von ex vivo µCT-Aufnahmen waren. Dort konnte eine Verdickung der Trabekelstruktur nachgewiesen werden, welche sehr gut mit einer Signalabnahme in den T1-gewichteten Aufnahmen korrelierte. Die radiale Auswertung der Daten erlaubte dabei eine komprimierte Darstellung der Ergebnisse. Dadurch wurde eine effiziente Auswertung der umfangreichen Daten (verschiedene Tiere an mehreren Zeitpunkten mit einer Vielzahl an Einzelaufnahmen) ermöglicht und eine unabhängige Bewertung erreicht.
Um die Limitationen der begrenzten Auflösung von 2D-Multi-Schichtaufnahmen aufzuheben, wurden neue Ansätze für eine hochaufgelöste 3D-Aufnahme entwickelt. Hierfür wurden Spin-Echo-basierte Sequenzen gewählt, da diese eine genauere Abbildung der Knochenmikrostruktur erlauben als Gradienten-Echo-basierte Methoden. Zum einen wurde eine eigene 3D-FLASE-Sequenz entwickelt und zum anderen eine modifizierte 3D-TSE-Sequenz. Damit an Patienten Aufnahmen bei klinischer Feldstärke von 1,5 T mit einer hohen räumlichen Auflösung innerhalb einer vertretbaren Zeit erzielt werden können, muss eine schnelle und signalstarke Sequenz verwendet werden. Eine theoretische Betrachtung bescheinigte der TSE-Sequenz eine um 25 % höhere Signaleffizienz verglichen mit einer FLASE-Sequenz mit identischer Messzeit. Dieser Unterschied konnte auch im Experiment nachgewiesen werden. Ein in vivo Vergleich der beiden Sequenzen am Schienbein zeigte eine vergleichbare Darstellung der Spongiosa mit einer Auflösung von 160 × 160 × 400 µm.
Für die Bildgebung des Hüftkopfs mit der neuen Sequenz waren jedoch aufgrund der unterschiedlichen Anatomie weitere Modifikationen notwendig. Um längere Messzeiten durch ein unnötig großes Field-of-View zu vermeiden, mussten Einfaltungsartefakte unterdrückt werden. Dies wurde durch die orthogonale Anwendung der Anregungs- und Refokussierungspulse in der TSE-Sequenz effizient gelöst. Technisch bedingt konnte jedoch nicht eine vergleichbare Auflösung wie am Schienbein realisiert werden.
Der Vorteil der 3D-Bildgebung, dass Schichtdicken von deutlich weniger als 1 mm erreicht werden können, konnte jedoch erfolgreich auf den Unterkiefer übertragen werden. Der dort verlaufende Nervus Mandibularis ist dabei eine wichtige Struktur, deren Verlauf im Vorfeld von verschiedenen operativen Eingriffen bekannt sein muss. Er ist durch eine dünne knöcherne Wand vom umgebenden Gewebe getrennt. Im Vergleich mit einer 3D-VIBE-Sequenz zeigte die entwickelte 3D-TSE-Sequenz mit integrierter Unterdrückung von Einfaltungsartefakten eine ähnlich gute Lokalisierung des Nervenkanals über die gesamte Länge der Struktur. Dies konnte in einer Studie an gesunden Probanden mit verschiedenen Beobachtern nachgewiesen werden. Durch die neue Aufnahmetechnik konnte darüber hinaus die Auflösung im Vergleich zu bisherigen Studien deutlich erhöht werden, was insgesamt eine präzisere Lokalisierung des Nervenkanals erlaubt.
Ein Baustein des Tissue Engineerings sind bio-resorbierbare Materialien, deren Abbau- und Einwachsverhalten noch untersucht werden muss, bevor diese für die klinische Anwendung zugelassen werden. Die durchgeführten in vitro µMR-Untersuchungen an Polymerscaffolds zeigten die reproduzierbare Quantifizierung der Porengröße und Wandstärke. Darüber hinaus wurde eine inhomogene Verteilung der Strukturparameter beobachtet. Die Ergebnisse waren in guter Übereinstimmung mit µCT-Aufnahmen als Goldstandard. Unterschiedliche Varianten der Scaffolds konnten identifiziert werden. Dabei bewies sich die MR-Bildgebung als zuverlässige Alternative.
Insgesamt zeigen die Ergebnisse dieser Arbeit, welche Vorteile und Anwendungsmöglichkeiten die 3D-MRT-Bildgebung bietet, und dass auch mit klinischer Feldstärke in vivo Voxelgrößen im Submillimeterbereich für alle Raumrichtungen erreichbar sind. Die erzielten Verbesserungen in der räumlichen Auflösung erhöhen die Genauigkeit der verschiedenen Anwendungen und ermöglichen eine bessere Identifikation von kleinen Abweichungen, was eine frühere und zuverlässigere Diagnose für Patienten verspricht.
In the present work, the energetic structure and coherence properties of the silicon vacancy point defect in the technologically important material silicon carbide are extensively studied by the optically detected magnetic resonance (ODMR) technique in order to verify its high potential for various quantum applications. In the spin vacancy, unique attributes are arising from the C3v symmetry and the spin-3/2 state, which are not fully described by the standard Hamiltonian of the uniaxial model. Therefore, an advanced Hamiltonian, describing well the appearing phenomena is established and the relevant parameters are experimentally determined. Utilizing these new accomplishments, several quantum metrology techniques are proposed.
First, a vector magnetometry scheme, utilizing the appearance of four ODMR lines, allows for simultaneous detection of the magnetic field strength and the tilting angle of the magnetic field from the symmetry axis of the crystal.
The second magnetometry protocol utilizes the appearance of energetic level anticrossings (LAC) in the ground state (GS) energy levels. Relying only on the change in photoluminescence in the vicinity of this GSLACs, this all-optical method does not require any radio waves and hence provides a much easier operation with less error sources as for the common magnetometry schemes utilizing quantum points.
A similar all-optical method is applied for temperature sensing, utilizing the thermal shift of the zero field splitting and consequently the anticrossing in the excited state (ES). Since the GSLACs show no dependence on temperature, the all-optical magnetometry and thermometry (utilizing the ESLACs) can be conducted subsequently on the same defect.
In order to quantify the achievable sensitivity of quantum metrology, as well as to prove the potential of the Si-vacancy in SiC for quantum processing, the coherence properties are investigated by the pulsed ODMR technique. The spin-lattice relaxation time T1 and the spin-spin relaxation time T2 are thoroughly analyzed for their dependence on the external magnetic field and temperature.
For actual sensing implementations, it is crucial to obtain the best signal-to-noise ratio without loss in coherence time. Therefore, the irradiation process, by which the defects are created in the crystal, plays a decisive role in the device performance. In the present work, samples irradiated with electrons or neutrons with different fluences and energies, producing different defect densities, are analyzed in regard to their T1 and T2 times at room temperature.
Last but not least, a scheme to substantially prolong the T2 coherence time by locking the spin polarization with the dynamic decoupling Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence is applied.
Bone marrow dosimetry is a topic of high interest in molecular radiotherapy. Predicting the level of hematological toxicity is one of the most important goals of nuclear medicine radiation dosimetry. To achieve this, it is necessary to quantify the absorbed dose to the active bone marrow, thus aiming at administering the most efficient therapy with a minimum level of adverse effects in the patient. The anatomical complexity of trabecular bone and bone marrow leads to the need of applying non-nuclear medicine imaging methods for determining the spatial distribution of soft tissue, adipose tissue, and bone in spongiosa.
Therefore, the two objectives of this dissertation are: i) to apply magnetic resonance imaging (MRI) for quantification of the fat volume fraction, and ii) to validate a method based on dual-energy quantitative computed tomography (DEQCT) for quantification of the trabecular bone volume fraction.
In a first step, an MRI sequence (two-point Dixon) for fat-water separation was validated in a 3 Tesla system by quantifying the fat volume fraction in a phantom and the lumbar vertebrae of volunteers and comparing with magnetic resonance spectroscopy (MRS). After successful validation, the fat volume fraction was retrospectively measured in the five lumbar vertebrae of 44 patient images acquired in the clinical routine. The two-point Dixon showed a good quantification of the fat volume fraction in the phantom experiment (-9.8% maximum relative error with respect to the nominal values). In the volunteers, a non-significant difference between MRI and MRS was found for the quantification of the fat volume fraction in volumes-of-interest with similar dimensions and position in both quantification methodologies (MRI and MRS). In the study with patient data, the marrow conversion (red → yellow marrow) was found to be age-dependent, and slower in males (0.3% per year) than in females (0.5% per year). Also, considerable variability of the fat volume fraction in patients of similar ages and the same gender was observed.
These results enable the use of two-point Dixon MRI in the quantification of the fat volume fraction in the bone marrow. Additionally, the constant marrow conversion during adulthood suggests that a patient-specific approach should replace the assumption of a constant cellularity volume fraction of 0.7 (reference man) (1,2) as proposed by the International Commission on Radiological Protection (ICRP).
In a second step, a quantification method based on DEQCT was validated in two CT systems: i) a clinical CT integrated into a SPECT/CT and ii) a dual-source computed tomography (DSCT) system. The method was applied in two phantoms: the first was used to validate the DEQCT method by the quantification of the hydroxyapatite volume fraction in three vials of 50 ml each and three different hydroxyapatite concentrations (100 mg/cm3, 200 mg/cm3, 300 mg/cm3). The second phantom was the European spine phantom (ESP), an anthropomorphic spine phantom. It was used to quantify the bone mineral content (BMC) on the whole vertebra and the hydroxyapatite volume fraction (VFHA) in the spongiosa region of each vertebra of the phantom. Lastly, the BMC of lumbar vertebrae 1 (LV1) and 2 (LV2) was measured in a patient using DEQCT and dual-energy X-ray absorptiometry (DEXA). Furthermore, the hydroxyapatite volume fraction (VFHA) and the bone volume fraction (VFB) was calculated for both the whole vertebrae and the spongiosa region of LV1 and LV2.
The measured and nominal hydroxyapatite volume fraction in the vial phantom showed a good correlation (maximum relative error: 14.2%). The quantification of the BMC on the whole vertebra and the VFHA on the spongiosa region showed larger relative errors than in the validation phantom. The quantification of BMC on LV1 and LV2 showed relative errors between DEXA and DSCT equal to 7.6% (LV1) and -8.4% (LV2). Also, the values of the VFHA (mineral bone) were smaller than the VFB. This result is consistent with the bone composition (mineral bone plus organic material).
The DEQCT method enables the quantification of hydroxyapatite (mineral bone) and bone (mineral bone plus organic material) in a clinical setting. However, the method showed an overestimation of the quantified mineral bone volume fraction. This overestimation might be related to the lack of detailed information on the CT X-ray spectra and detector sensitivity. Also, the DEQCT method showed a dependency on the CT reconstruction kernel and the chemical description of the materials to be quantified.
Based on the results of this work, the feasibility for quantifying the fat volume fraction and the bone volume fraction in the spongiosa in a clinical setting has been demonstrated/proven. Furthermore, the differences in fat volume fraction in females and males, as well as the variability of the fat volume fraction in subjects of similar ages, questions the approximation of the cellularity volume fraction by only a single ICRP reference value in bone marrow dosimetry for molecular radiotherapy. Lastly, this study presents the first approach for non-invasive quantification of the bone volume fraction (mineral bone plus organic material) for improved bone marrow dosimetry.
Magnetic Resonance Imaging at field strengths up to 3 T, has become a default diagnostic modality for a variety of disorders and injuries, due to multiple reasons ranging from its non-invasive nature to the possibility of obtaining high resolution images of internal organs and soft tissues. Despite tremendous advances, MR imaging of certain anatomical regions and applications present specific challenges to be overcome. One such application is MR Musculo-Skeletal Imaging. This work addresses a few difficult areas within MSK imaging from the hardware perspective, with coil solutions for dynamic imaging of knee and high field imaging of hand.
Starting with a brief introduction to MR physics, different types of RF coils are introduced in chapter 1, followed by sections on design of birdcage coils, phased arrays and their characterization in chapter 2. Measurements, calculations and simulations, done during the course of this work, have been added to this chapter to give a quantitative feel of the concepts explained.
Chapter 3 deals with the construction of a phased array receiver for dynamic imaging of knee of a large animal model, i.e. minipig, at 1.5 T. Starting with details on the various aspects of an application that need to be considered when an MR RF array is designed, the chapter details the complex geometry of the region of interest in a minipig and reasons that necessitate a high density array. The sizes of the individual elements that constitute the array have been arrived at by studying the ratio of unloaded to loaded Q factors and choosing a size that provides the best ratio but still maintains a uniform SNR throughout the movement of the knee. To have a minimum weight and to allow mechanical movement of the knee, the Preamplifiers were located in a separate box. A movement device was constructed to achieve adjustable periodic movement of the knee of the anesthetized animal. The constructed array has been characterized for its SNR and compared with an existing product coil to show the improvement. The movement device was also characterized for its reproducibility. High resolution static images with anatomical details marked have been presented. The 1/g maps show the accelerations possible with the array. Snapshots of obtained dynamic images trace the cruciate ligaments through a cycle of movement of the animal's knee.
The hardware combination of a high density phased array and a movement device designed for a minipig's knee was used as a 'reference' and extended in chapter 4 for a human knee. In principle the challenges are similar for dynamic imaging of a human knee with regards to optimization of the elements, the associated electronics and the construction of the movement device. The size of the elements were optimized considering the field penetration / sensitivity required for the internal tissues. They were distributed around the curvature of the knee keeping in mind the acceleration required for dynamic imaging and the direction of the movement. The constructed movement device allows a periodic motion of the lower half of the leg, with the knee placed within the coil, enabling visualization of the tissues inside, while the leg is in motion. Imaging has been performed using dynamic interleaved acquisition sequence where higher effective TR and flip angles are achieved due to a combination of interleaving and segmentation of the sequence. The movement device has been characterized for its reproducibility while the SNR distribution of the constructed RF array has been compared with that of a commercially available standard 8 channel array. The results show the improvement in SNR and acceleration with the constructed geometry. High resolution static images, dynamic snapshots and the 3D segmentation of the obtained images prove the usefulness of the complete package provided in the design, for performing dynamic imaging at a clinically relevant field strength.
A simple study is performed in chapter 5 to understand the effects of changes in overlap for coil configurations with different loads and at different frequencies. The noise levels of individual channels and the correlation between them are plotted against subtle changes in overlap, at 64 and 123 MHz. SNR for every overlap setup is also measured and plotted. Results show that achieving critical overlap is crucial to obtain the best possible SNR in those coil setups where the load offered by the sample is low.
Chapter 6 of the thesis work deals with coil design for high field imaging of hand and wrists at 7 T, with an aim to achieve ultra high resolution imaging. At this field strength due to the increase in dielectric effects and the resulting decrease in homogeneity, whole body transmit coils are impractical and this has led engineers to design local transmit coils, for specific anatomies. While transmit or transceive arrays are usually preferred, to mitigate SAR effects, the spatial resolution obtained is limited. It is shown that a solution to this, with regards to hand imaging, can be a single volume transmit coil, along with high density receive arrays optimized for different regions of the hand. The use of a phased array for reception provides an increased SNR / penetration under high resolution. A volume transmit coil could pose issues in homogeneity at 7 T, but the specific anatomy of hand and wrist, with comparatively less water content, limits dielectric effects to have homogeneous B_1+ profile over the hand. To this effect, a bandpass birdcage and a 12 channel receive array are designed and characterized. Images of very high spatial resolution (0.16 x 0.16 x 0.16 mm3) with internal tissues marked are presented. In vivo 1/g maps show that an acceleration of up to 3 is possible and the EM simulation results presented show the uniform field along with SAR hotspots in the hand. To reduce the stress created due to the 'superman' position of imaging, provisions in the form of a holder and a hand rest have been designed and presented. Factors that contributed to the stability of the presented design are also listed, which would help future designs of receive arrays at high field strengths.
In conclusion, the coils and related hardware presented in this thesis address the following two aspects of MSK imaging: Dynamic imaging of knee and High resolution imaging of hand / wrist. The presented hardware addresses specific challenges and provides solutions. It is hoped that these designs are steps in the direction of improving the existing coils to get a better knowledge and understanding of MSK diseases such as Rheumatoid Arthritis and Osteoarthritis. The hardware can aid our study of ligament reconstruction and development. The high density array and transmit coil design for hand / wrist also demonstrates the benefits of the obtained SNR at 7 T while maintaining SAR within limits. This design is a contribution towards optimizing hardware at high field strength, to make it clinically acceptable and approved by regulatory bodies.
Although the contribution to the Isotropic Gamma-Ray Background (IGRB) from unresolved extragalactic objects has been studied for many years, its exact composition and origin are as of yet unknown. It is suspected that diffuse processes such as dark matter annihilation contribute to the total IGRB, as well as unresolved gamma-ray emission from Active Galactic Nuclei (AGN), including radio galaxies. Radio galaxies are a source class that emit strongly at radio wavelengths, some of which have also been detected at gamma-ray wavelengths by the Fermi Large Area Telescope (Fermi-LAT), and by very high energy gamma-ray Cherenkov telescopes. It is thought that due to the orientation of their jets, radio galaxies are detected less numerously at gamma-ray energies than blazars. Furthermore, only a small number of radio galaxies have been detected at gamma-ray energies though it is considered that others do as well. It is for these reasons that gamma-ray emitting radio galaxies, an interesting and elusive class of objects, are selected for investigation in this work.
In order to reach the goal of better understanding diffuse processes, it is necessary to model the radio galaxy spectral energy distributions (SEDs). As AGN emission is variable with respect to time, it is critical to use simultaneously collected observations. Calculation of the SED based on simultaneous, multiwavelength data across the electromagnetic spectrum produces a reasonably accurate representation of the state of an object in a given time range. The gamma-ray emitting radio galaxies M 87, NGC 1275, Pictor A, and Centaurus A are selected here based on having been detected in very high energy gamma-rays by Cherenkov telescopes, as well as in other wavelengths. A uniquely consistent analysis approach is applied, in which each radio galaxy is analyzed the same way using simultaneously collected data. This approach sets it apart from other studies.
Fermi-LAT raw data for each source in the sample is analyzed in time ranges which directly overlap the very high energy gamma-ray Cherenkov observations, as well as several other wavelength ranges. A synchrotron self-Compton (SSC) model is applied, which provides accurate treatment of synchrotron and inverse-Compton processes occurring in the jets of AGN, while estimating physical characteristics of the source. It is found that the spectra of M 87, NGC 1275, Pictor A, and Centaurus A can be well described by the same SSC model, producing values for the physical characteristics such as the doppler factor and magnetic field, which are relatively consistent with each other.
In order to characterize the diffuse emission from dark matter self-annihilation, the radio galaxy SEDs are also fit with a dark matter model, resulting in an estimated dark matter particle mass of around 4.7 TeV which lies within predicted ranges.
The highly dense regions near the black holes of AGN provide the optimal conditions for detecting these signatures. It is also found here that discrepancies between the expected emission and the observed emission in the spectra of some radio galaxies can be explained using the combined SSC and dark matter model. As emission from dark matter annihilation is expected to remain steady with respect to time, a key feature of this work is the novelty of the combined SSC and dark matter model, and the finding that dark matter characteristics may be revealed through similar multiwavelength analyses during future low emission states of the AGN.
The radio galaxy sample is then extended to include all gamma-ray emitting radio galaxies detected by the Fermi-LAT, and a calculation of the core radio, total radio, and gamma-ray luminosities is followed through. A future step in extending this work would be to estimate the gamma-ray luminosity function of radio galaxies and their percent contribution to the total IGRB, based on the widely agreed upon assumption that a reasonable estimate of the gamma-ray luminosity function of a population can be attained by appropriately scaling its radio luminosity function, as gamma-ray luminosities and radio luminosities are strongly linearly correlated. This work has also provided the basis for such a calculation by outlining the theory and initial steps.
It is the hope that the vast scope of the gathered data, its simultaneity, and the use of consistent analysis methods across the sample, will provide an improved foundation for a future calculation of the contribution of this population to the IGRB, as well as encourage stricter requirements for multiwavelength studies.
As a non-destructive testing method, X-ray imaging has proved to be suitable for the examination of a variety of objects. The measurement principle is based on the attenuation of X-rays caused by these objects. This attenuation can be recorded as shades of intensity using X-ray detectors and thus contains information about the inner structure of the investigated object. Since X-rays are electromagnetic waves, they also experience a change of phase in addition to their attenuation while penetrating an object. In general, imaging methods based on this effect are referred to as phase contrast imaging techniques. In the laboratory, the two mainly used methods are the propagation based phase contrast or in-line phase contrast and the grating interferometry.
While in-line phase contrast - under certain conditions - shows edge enhancement at interfaces due to interference, phase contrast in the grating interferometry is only indirectly measurable by the use of several gratings. In addition to phase contrast, grating interferometry provides access to the so-called dark-field imaging contrast, which measures the scattering of X-rays caused by an object.
These two imaging techniques, together with a novel concept of laboratory X-ray sources, the liquid-metal-jet, form the main part of this work. Compared to conventional X-ray sources, the liquid-metal-jet source offers higher brightness. The term brightness is defined by the number of X-ray photons per second, emitting area (area of the X-ray spot) and solid angle at which they are emitted.
On the basis of this source, a high resolution in-line phase contrast setup was partially developed in the scope of this work. Several computed tomographies show the feasibility of in-line phase contrast and the improvement of image quality by applying phase retrieval algorithms.
Moreover, the determination of optimized sample positions for in-line phase contrast imaging is treated at which the edge enhancement is maximized. Based on primitive fiber objects, this optimization has proven to be a good approximation.
With its high brightness in combination with a high spatial coherence, the liquid-metal-jet source is also interesting for grating interferometry. The development of such a setup is also part of this work. The overall concept and the characterization of the setup is presented as well as the applicability and its limits for the investigation of various objects.
Due to the very unique concept of this grating interferometer it was possible to realize a modified interferometer system by using a single grating only. Its concept and results are also presented in this work.
Furthermore, a grating interferometer based on a microfocus X-ray tube was tested regarding its performance. Thereby, parameters like the anode material, acquisition geometry and gratings were altered in order to find the advantages and disadvantages of each configuration.
Für die Verwendung von zellbasierten Therapeutika ist vor allem die korrekt Identifikation
sowohl vom Ausgangsmaterial wie auch dem produziertem Material von
zentraler Wichtigkeit. In dieser Arbeit wurde eine Methodik entwickelt, welche eine
nicht-invasive Klassifizierung von Zellen und zellulärer Entwicklung aufgrund ihrer
zweidimensionalen Magnetresonanz-Korrelationsspektren ermöglichte.
Hierzu wurde ein mobiler MR-Scanner mit einer Feldstärke von 0.5T und einem Isozentrum
von 1 cm3 verwendet. Aufgrund der kompakten und leichten Bauweise war
es möglich, das System in normalen Zellkulturlaboren zu verwenden. Von den Proben
wurde ein zweidimensionales T1/T2 -Korrelationsspektrum aufgenommen, anhand
dessen die Zellen klassifiziert werden sollten. Mithilfe von Agarose-Dotagraf® -Zell-
Phantomen konnte die Stabilität und Reproduzierbarkeit des Messsystems und der
verwendeten Sequenz validiert werden.
Aufgrund der unter Umständen recht langen Messzeiten der MR-Technologie war
auch die Handhabung und Kultur der Zellproben während des Messprozesses von
großer Bedeutung. Um hierfür den Durchsatz an Proben zu erhöhen, wurde eine kostengünstige
und ebenfalls mobile Robotikanlage entwickelt. Diese basierte auf dem
kommerziell erhältlichen Roboterarm Braccio, welcher durch einen Arduino Mega
Mikrocontroller gesteuert wurde. Mit bis zu 24 Proben pro Tag konnte durch die
Automatisierung der Durchsatz an Proben um den Faktor 3 – 4 gesteigert werden.
Durch den entwickelten Prozess war es möglich, eine umfangreiche Datenbank –
bestehend aus 362 unabhängigen Messungen (biologische Replikate) – aufzubauen.
Die Datenbank enthielt Messungen von zehn unterschiedlichen Zelllinien. Zusätzlich
wurden T1/T2 -Korrelationsspektren von mesenchymalen Stromazellen (MSCs)
vor und nach deren Differenzierung zu Adipocyten aufgenommen, um ihre zelluläre
Entwicklung nicht-invasiv charakterisieren zu können.
Die aufgenommenen Daten wurden mithilfe einer geeigneten Support Vector Machine
wie auch angepassten künstlichen neuronalen Netzwerken klassifiziert. Mithilfe
dieser Methoden konnten die Zelllinien und MSCs anhand ihrer aufgenommenen
Korrelationsspektren mit einer Genauigkeit von bis zu 98% klassifiziert werden.
Diese hohe Treffsicherheit legte den Schluss nahe, dass die Kombination aus nichtinvasiver,
zweidimensionaler T1/T2 -MR-Relaxometrie und der Verwendung von geeigneten
Methoden des machine learning und der künstlichen Intelligenz eine effiziente
Methodik für die nicht-invasive Klassifizierung von Zellen sowie zellulärer
Entwicklung darstellt.