Refine
Has Fulltext
- yes (20)
Is part of the Bibliography
- yes (20)
Year of publication
Document Type
- Doctoral Thesis (15)
- Journal article (5)
Keywords
- MRI (20) (remove)
Institute
- Physikalisches Institut (20) (remove)
Sonstige beteiligte Institutionen
Within this thesis, three main approaches for the assessment and investigation of altered hemodynamics like wall shear stress, oscillatory shear index and the arterial pulse wave velocity in atherosclerosis development and progression were conducted:
1. The establishment of a fast method for the simultaneous assessment of 3D WSS and PWV in the complete murine aortic arch via high-resolution 4D-flow MRI
2. The utilization of serial in vivo measurements in atherosclerotic mouse models using high-resolution 4D-flow MRI, which were divided into studies describing altered hemodynamics in late and early atherosclerosis
3. The development of tissue-engineered artery models for the controllable application and variation of hemodynamic and biologic parameters, divided in native artery models and biofabricated artery models, aiming for the investigation of the relationship between atherogenesis and hemodynamics
Chapter 2 describes the establishment of a method for the simultaneous measurement of 3D WSS and PWV in the murine aortic arch at, using ultra high-field MRI at 17.6T [16], based on the previously published method for fast, self-navigated wall shear stress measurements in the murine aortic arch using radial 4D-phase contrast MRI at 17.6 T [4]. This work is based on the collective work of Dr. Patrick Winter, who developed the method and the author of this thesis, Kristina Andelovic, who performed the experiments and statistical analyses. As the method described in this chapter is basis for the following in vivo studies and undividable into the sub-parts of the contributors without losing important information, this chapter was not split into the single parts to provide fundamental information about the measurement and analysis methods and therefore better understandability for the following studies. The main challenge in this chapter was to overcome the issue of the need for a high spatial resolution to determine the velocity gradients at the vascular wall for the WSS quantification and a high temporal resolution for the assessment of the PWV without prolonging the acquisition time due to the need for two separate measurements. Moreover, for a full coverage of the hemodynamics in the murine aortic arch, a 3D measurement is needed, which was achieved by utilization of retrospective navigation and radial trajectories, enabling a highly flexible reconstruction framework to either reconstruct images at lower spatial resolution and higher frame rates for the acquisition of the PWV or higher spatial resolution and lower frame rates for the acquisition of the 3D WSS in a reasonable measurement time of only 35 minutes. This enabled the in vivo assessment of all relevant hemodynamic parameters related to atherosclerosis development and progression in one experimental session. This method was validated in healthy wild type and atherosclerotic Apoe-/- mice, indicating no differences in robustness between pathological and healthy mice.
The heterogeneous distribution of plaque development and arterial stiffening in atherosclerosis [10, 12], however, points out the importance of local PWV measurements. Therefore, future studies should focus on the 3D acquisition of the local PWV in the murine aortic arch based on the presented method, in order to enable spatially resolved correlations of local arterial stiffness with other hemodynamic parameters and plaque composition.
In Chapter 3, the previously established methods were used for the investigation of changing aortic hemodynamics during ageing and atherosclerosis in healthy wild type and atherosclerotic Apoe-/- mice using the previously established methods [4, 16] based on high-resolution 4D-flow MRI. In this work, serial measurements of healthy and atherosclerotic mice were conducted to track all changes in hemodynamics in the complete aortic arch over time. Moreover, spatially resolved 2D projection maps of WSS and OSI of the complete aortic arch were generated. This important feature allowed for the pixel-wise statistical analysis of inter- and intragroup hemodynamic changes over time and most importantly – at a glance. The study revealed converse differences of local hemodynamic profiles in healthy WT and atherosclerotic Apoe−/− mice, with decreasing longWSS and increasing OSI, while showing constant PWV in healthy mice and increasing longWSS and decreasing OSI, while showing increased PWV in diseased mice. Moreover, spatially resolved correlations between WSS, PWV, plaque and vessel wall characteristics were enabled, giving detailed insights into coherences between hemodynamics and plaque composition. Here, the circWSS was identified as a potential marker of plaque size and composition in advanced atherosclerosis. Moreover, correlations with PWV values identified the maximum radStrain could serve as a potential marker for vascular elasticity. This study demonstrated the feasibility and utility of high-resolution 4D flow MRI to spatially resolve, visualize and analyze statistical differences in all relevant hemodynamic parameters over time and between healthy and diseased mice, which could significantly improve our understanding of plaque progression towards vulnerability. In future studies the relation of vascular elasticity and radial strain should be further investigated and validated with local PWV measurements and CFD.
Moreover, the 2D histological datasets were not reflecting the 3D properties and regional characteristics of the atherosclerotic plaques. Therefore, future studies will include 3D plaque volume and composition analysis like morphological measurements with MRI or light-sheet microscopy to further improve the analysis of the relationship between hemodynamics and atherosclerosis.
Chapter 4 aimed at the description and investigation of hemodynamics in early stages of atherosclerosis. Moreover, this study included measurements of hemodynamics at baseline levels in healthy WT and atherosclerotic mouse models. Due to the lack of hemodynamic-related studies in Ldlr-/- mice, which are the most used mouse models in atherosclerosis research together with the Apoe-/- mouse model, this model was included in this study to describe changing hemodynamics in the aortic arch at baseline levels and during early atherosclerosis development and progression for the first time. In this study, distinct differences in aortic geometries of these mouse models at baseline levels were described for the first time, which result in significantly different flow- and WSS profiles in the Ldlr-/- mouse model. Further basal characterization of different parameters revealed only characteristic differences in lipid profiles, proving that the geometry is highly influencing the local WSS in these models. Most interestingly, calculation of the atherogenic index of plasma revealed a significantly higher risk in Ldlr-/- mice with ongoing atherosclerosis development, but significantly greater plaque areas in the aortic arch of Apoe-/- mice. Due to the given basal WSS and OSI profile in these two mouse models – two parameters highly influencing plaque development and progression – there is evidence that the regional plaque development differs between these mouse models during very early atherogenesis.
Therefore, future studies should focus on the spatiotemporal evaluation of plaque development and composition in the three defined aortic regions using morphological measurements with MRI or 3D histological analyses like LSFM. Moreover, this study offers an excellent basis for future studies incorporating CFD simulations, analyzing the different measured parameter combinations (e.g., aortic geometry of the Ldlr-/- mouse with the lipid profile of the Apoe-/- mouse), simulating the resulting plaque development and composition. This could help to understand the complex interplay between altered hemodynamics, serum lipids and atherosclerosis and significantly improve our basic understanding of key factors initiating atherosclerosis development.
Chapter 5 describes the establishment of a tissue-engineered artery model, which is based on native, decellularized porcine carotid artery scaffolds, cultured in a MRI-suitable bioreactor-system [23] for the investigation of hemodynamic-related atherosclerosis development in a controllable manner, using the previously established methods for WSS and PWV assessment [4, 16]. This in vitro artery model aimed for the reduction of animal experiments, while simultaneously offering a simplified, but completely controllable physical and biological environment. For this, a very fast and gentle decellularization protocol was established in a first step, which resulted in porcine carotid artery scaffolds showing complete acellularity while maintaining the extracellular matrix composition, overall ultrastructure and mechanical strength of native arteries. Moreover, a good cellular adhesion and proliferation was achieved, which was evaluated with isolated human blood outgrowth endothelial cells. Most importantly, an MRI-suitable artery chamber was designed for the simultaneous cultivation and assessment of high-resolution 4D hemodynamics in the described artery models. Using high-resolution 4D-flow MRI, the bioreactor system was proven to be suitable to quantify the volume flow, the two components of the WSS and the radStrain as well as the PWV in artery models, with obtained values being comparable to values found in literature for in vivo measurements. Moreover, the identification of first atherosclerotic processes like intimal thickening is achievable by three-dimensional assessment of the vessel wall morphology in the in vitro models. However, one limitation is the lack of a medial smooth muscle cell layer due to the dense ECM. Here, the utilization of the laser-cutting technology for the generation of holes and / or pits on a microscale, eventually enabling seeding of the media with SMCs showed promising results in a first try and should be further investigated in future studies. Therefore, the proposed artery model possesses all relevant components for the extension to an atherosclerosis model which may pave the way towards a significant improvement of our understanding of the key mechanisms in atherogenesis.
Chapter 6 describes the development of an easy-to-prepare, low cost and fully customizable artery model based on biomaterials. Here, thermoresponsive sacrificial scaffolds, processed with the technique of MEW were used for the creation of variable, biomimetic shapes to mimic the geometric properties of the aortic arch, consisting of both, bifurcations and curvatures. After embedding the sacrificial scaffold into a gelatin-hydrogel containing SMCs, it was crosslinked with bacterial transglutaminase before dissolution and flushing of the sacrificial scaffold. The hereby generated channel was subsequently seeded with ECs, resulting in an easy-to-prepare, fast and low-cost artery model. In contrast to the native artery model, this model is therefore more variable in size and shape and offers the possibility to include smooth muscle cells from the beginning. Moreover, a custom-built and highly adaptable perfusion chamber was designed specifically for the scaffold structure, which enabled a one-step creation and simultaneously offering the possibility for dynamic cultivation of the artery models, making it an excellent basis for the development of in vitro disease test systems for e.g., flow-related atherosclerosis research. Due to time constraints, the extension to an atherosclerosis model could not be achieved within the scope of this thesis. Therefore, future studies will focus on the development and validation of an in vitro atherosclerosis model based on the proposed bi- and three-layered artery models.
In conclusion, this thesis paved the way for a fast acquisition and detailed analyses of changing hemodynamics during atherosclerosis development and progression, including spatially resolved analyses of all relevant hemodynamic parameters over time and in between different groups. Moreover, to reduce animal experiments, while gaining control over various parameters influencing atherosclerosis development, promising artery models were established, which have the potential to serve as a new platform for basic atherosclerosis research.
Increased aortic stiffness is known to be associated with atherosclerosis and has a predictive value for cardiovascular events. This study aims to investigate the local distribution of early arterial stiffening due to initial atherosclerotic lesions. Therefore, global and local pulse wave velocity (PWV) were measured in ApoE\(^{-/-}\) and wild type (WT) mice using ultrahigh field MRI. For quantification of global aortic stiffness, a new multi-point transit-time (TT) method was implemented and validated to determine the global PWV in the murine aorta. Local aortic stiffness was measured by assessing the local PWV in the upper abdominal aorta, using the flow/area (QA) method. Significant differences between age matched ApoE\(^{-/-}\) and WT mice were determined for global and local PWV measurements (global PWV: ApoE\(^{-/-}\): 2.7 ±0.2m/s vs WT: 2.1±0.2m/s, P<0.03; local PWV: ApoE\(^{-/-}\): 2.9±0.2m/s vs WT: 2.2±0.2m/s, P<0.03). Within the WT mouse group, the global PWV correlated well with the local PWV in the upper abdominal aorta (R\(^2\) = 0.75, P<0.01), implying a widely uniform arterial elasticity.
In ApoE\(^{-/-}\) animals, however, no significant correlation between individual local and global PWV was present (R\(^2\) = 0.07, P = 0.53), implying a heterogeneous distribution of vascular stiffening in early atherosclerosis. The assessment of global PWV using the new multi-point TT measurement technique was validated against a pressure wire measurement in a vessel
phantom and showed excellent agreement. The experimental results demonstrate that vascular stiffening caused by early atherosclerosis is unequally distributed over the length of large vessels. This finding implies that assessing heterogeneity of arterial stiffness by multiple local measurements of PWV might be more sensitive than global PWV to identify early atherosclerotic lesions.
Nuclear Magnetic Resonance (NMR) provides a highly flexible platform for non invasive analysis and imaging biological samples, since the manipulation of nuclear spin allows the tailoring of experiments to maximize the informativeness of the data. MRI is capable of visualizing a holistic picture of the lipid storage in living plant/seed. This review has sought to explain how the technology can be used to acquire functional and physiological data from plant samples, and how to exploit it to characterize lipid deposition in vivo. At the same time, we have referred to the current limitations of NMR technology as applied to plants, and in particular of the difficulty of transferring methodologies optimized for animal/medical subjects to plant ones. A forward look into likely developments in the field is included, anticipating its key future role in the study of living plant.
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.
Low field NMR has been successfully used for the evaluation of seed composition and quality, but largely only in crop species. We show here that 1.5T NMR provides a reliable means for analysing the seed lipid fraction present in a wide range of species, where both the seed size and lipid concentration differed by >10 fold. Little use of high field NMR has been made in seed research to date, even though it potentially offers many opportunities for studying seed development, metabolism and storage. Here we demonstrate how 17.5T and 20T NMR can be applied to image seed structure, and analyse lipid and metabolite distribution. We suggest that further technical developments in NMR/MRI will facilitate significant advances in our understanding of seed biology.
Direct cooling of the catheter tip increases safety for CMR-guided electrophysiological procedures
(2012)
Background: One of the safety concerns when performing electrophysiological (EP) procedures under magnetic resonance (MR) guidance is the risk of passive tissue heating due to the EP catheter being exposed to the radiofrequency (RF) field of the RF transmitting body coil. Ablation procedures that use catheters with irrigated tips are well established therapeutic options for the treatment of cardiac arrhythmias and when used in a modified mode might offer an additional system for suppressing passive catheter heating.
Methods: A two-step approach was chosen. Firstly, tests on passive catheter heating were performed in a 1.5 T Avanto system (Siemens Healthcare Sector, Erlangen, Germany) using a ASTM Phantom in order to determine a possible maximum temperature rise. Secondly, a phantom was designed for simulation of the interface between blood and the vascular wall. The MR-RF induced temperature rise was simulated by catheter tip heating via a standard ablation generator. Power levels from 1 to 6 W were selected. Ablation duration was 120 s with no tip irrigation during the first 60 s and irrigation at rates from 2 ml/min to 35 ml/min for the remaining 60 s (Biotronik Qiona Pump, Berlin, Germany). The temperature was measured with fluoroscopic sensors (Luxtron, Santa Barbara, CA, USA) at a distance of 0 mm, 2 mm, 4 mm, and 6 mm from the catheter tip. Results: A maximum temperature rise of 22.4 degrees C at the catheter tip was documented in the MR scanner. This temperature rise is equivalent to the heating effect of an ablator's power output of 6 W at a contact force of the weight of 90 g (0.883 N). The catheter tip irrigation was able to limit the temperature rise to less than 2 degrees C for the majority of examined power levels, and for all examined power levels the residual temperature rise was less than 8 degrees C.
Conclusion: Up to a maximum of 22.4 degrees C, the temperature rise at the tissue surface can be entirely suppressed by using the catheter's own irrigation system. The irrigated tip system can be used to increase MR safety of EP catheters by suppressing the effects of unwanted passive catheter heating due to RF exposure from the MR scanner.
The focus of the work concerned the development of a series of MRI techniques that were specifically designed and optimized to obtain quantitative and spatially resolved information about characteristic parameters of the lung. Three image acquisition techniques were developed. Each of them allows to quantify a different parameter of relevant diagnostic interest for the lung, as further described below:
1) The blood volume fraction, which represents the amount of lung water in the intravascular compartment expressed as a fraction of the total lung water. This parameter is related to lung perfusion.
2) The magnetization relaxation time T\(_2\) und T*\(_2\)
, which represents the component of T\(_2\) associated with the diffusion of water molecules through the internal magnetic field gradients of the lung. Because the amplitude of these internal gradients is related to the alveolar size, T\(_2\) und T*\(_2\) can be used to obtain information about the microstructure of the lung.
3) The broadening of the NMR spectral line of the lung. This parameter depends on lung inflation and on the concentration of oxygen in the alveoli. For this reason, the spectral line broadening can be regarded as a fingerprint for lung inflation; furthermore, in combination with oxygen enhancement, it provides a measure for lung ventilation.
Proton magnetic resonance imaging (MRI) has recently emerged as a clinical tool to image the lungs. This paper outlines the current technical aspects of MRI pulse sequences, radiofrequency (RF) coils and MRI system requirements needed for imaging the pulmonary parenchyma and vasculature. Lung MRI techniques are presented as a “technical toolkit”, from which MR protocols will be composed in the subsequent papers for comprehensive imaging of lung disease and function (parts 2 and 3). This paper is pitched at MR scientists, technicians and radiologists who are interested in understanding and establishing lung MRI methods. Images from a 1.5 T scanner are used for illustration of the sequences and methods that are highlighted.
Main Messages
• Outline of the hardware and pulse sequence requirements for proton lung MRI
• Overview of pulse sequences for lung parenchyma, vascular and functional imaging with protons
• Demonstration of the pulse-sequence building blocks for clinical lung MRI protocols
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 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.