Fast joint reconstruction of dynamic R2∗R_2^* and field maps in functional MRI.

Blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) is conventionally done by reconstructing T2 * -weighted images. However, since the images are unitless they are nonquantifiable in terms of important physiological parameters. An alternative approach is to reconstruct R2 * maps which are quantifiable and have comparable BOLD contrast as T2* -weighted images. However, conventional R2 * mapping involves long readouts and ignores relaxation during readout. Another problem with fMRI imaging is temporal drift/fluctuations in off-resonance. Conventionally, a field map is collected at the start of the fMRI study to correct for off-resonance, ignoring any temporal changes. Here, we propose a new fast regularized iterative algorithm that jointly reconstructs R2 * and field maps for all time frames in fMRI data. To accelerate the algorithm we linearize the MR signal model, enabling the use of fast regularized iterative reconstruction methods. The regularizer was designed to account for the different resolution properties of both R2 * and field maps and provide uniform spatial resolution. For fMRI data with the same temporal frame rate as data collected for T2 * -weighted imaging the resulting R2 * maps performed comparably to T2 * -weighted images in activation detection while also correcting for spatially global and local temporal changes in off-resonance.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/86002/1/Fessler23.pd

Lung Imaging and Function Assessment using Non-Contrast-Enhanced Magnetic Resonance Imaging

Measurement of pulmonary ventilation and perfusion has significant clinical value for the diagnosis and monitoring of prevalent lung diseases. To this end, non-contrast-enhanced MRI techniques have emerged as a promising alternative to scintigraphical measurements, computed tomography, and contrast-enhanced MRI. Although these techniques allow the acquisition of both structural and functional information in the same scan session, they are prone to robustness issues related to imaging artifacts and post-processing techniques, limiting their clinical utilization. In this work, new acquisition and post-processing techniques were introduced for improving the robustness of non-contrast-enhanced MRI based functional lung imaging. Furthermore, pulmonary functional maps were acquired in 2-year-old congenital diaphragmatic hernia (CDH) patients to demonstrate the feasibility of non-contrast-enhanced MRI methods for functional lung imaging. In the first study, a multi-acquisition framework was developed to improve robustness against field inhomogeneity artifacts. This method was evaluated at 1.5T and 3T field strengths via acquisitions obtained from healthy volunteers. The results demonstrate that the proposed acquisition framework significantly improved ventilation map homogeneity p<0.05. In the second study, a post-processing method based on dynamic mode decomposition (DMD) was developed to accurately identify dominant spatiotemporal patterns in the acquisitions. This method was demonstrated on digital lung phantoms and in vivo acquisitions. The findings indicate that the proposed method led to a significant reduction in dispersion of estimated ventilation and perfusion map amplitudes across different number of measurements when compared with competing methods p<0.05. In the third study, the free-breathing non-contrast-enhanced dynamic acquisitions were obtained from 2-year-old patients after CDH repair, and then processed using the DMD to obtain pulmonary functional maps. Afterwards, functional differences between ipsilateral and contralateral lungs were assessed and compared with results obtained using contrast-enhanced MRI measurements. The results demonstrate that pulmonary ventilation and perfusion maps can be generated from dynamic acquisitions successfully without the need for ionizing radiation or contrast agents. Furthermore, lung perfusion parameters obtained with DMD MRI correlate very strongly with parameters obtained using dynamic contrast-enhanced MRI. In conclusion, the presented work improves the robustness and accuracy of non-contrast-enhanced functional lung imaging using MRI. Overall, the methods introduced in this work may serve as a valuable tool in the clinical adaptation of non-contrast-enhanced imaging methods and may be used for longitudinal assessments of pulmonary functional changes

Low-rank and sparse reconstruction in dynamic magnetic resonance imaging via proximal splitting methods

Dynamic magnetic resonance imaging (MRI) consists of collecting multiple MR images in time, resulting in a spatio-temporal signal. However, MRI intrinsically suffers from long acquisition times due to various constraints. This limits the full potential of dynamic MR imaging, such as obtaining high spatial and temporal resolutions which are crucial to observe dynamic phenomena. This dissertation addresses the problem of the reconstruction of dynamic MR images from a limited amount of samples arising from a nuclear magnetic resonance experiment. The term limited can be explained by the approach taken in this thesis to speed up scan time, which is based on violating the Nyquist criterion by skipping measurements that would be normally acquired in a standard MRI procedure. The resulting problem can be classified in the general framework of linear ill-posed inverse problems. This thesis shows how low-dimensional signal models, specifically lowrank and sparsity, can help in the reconstruction of dynamic images from partial measurements. The use of these models are justified by significant developments in signal recovery techniques from partial data that have emerged in recent years in signal processing. The major contributions of this thesis are the development and characterisation of fast and efficient computational tools using convex low-rank and sparse constraints via proximal gradient methods, the development and characterisation of a novel joint reconstruction–separation method via the low-rank plus sparse matrix decomposition technique, and the development and characterisation of low-rank based recovery methods in the context of dynamic parallel MRI. Finally, an additional contribution of this thesis is to formulate the various MR image reconstruction problems in the context of convex optimisation to develop algorithms based on proximal splitting methods

Dynamic Image and Fieldmap Joint Estimation Methods for MRI Using Single-Shot Trajectories.

In susceptibility-weighted MRI, ignoring the magnetic field inhomogeneity can lead to severe reconstruction artifacts. Correcting for the effects of magnetic field inhomogeneity requires accurate fieldmaps. Especially in functional MRI, dynamic updates are desirable, since the fieldmap may change in time. Also, susceptibility effects that induce field inhomogeneity often have non-zero through-plane gradients, which, if uncorrected, can cause signal loss in the reconstructed images. Most image reconstruction methods that compensate for field inhomogeneity, even using dynamic fieldmap updates, ignore through-plane fieldmap gradients. Furthermore, standard optimization methods, like CG-based algorithms, may be slow to converge and recently proposed algorithms based on the Augmented Lagrangian (AL) framework have shown the potential to lead to more efficient optimization algorithms, especially in MRI reconstruction problems with non-quadratic regularization. In this work, we propose a computationally efficient, model-based iterative method for joint reconstruction of dynamic images and fieldmaps in single coil and parallel MRI, using single-shot trajectories. We first exploit the fieldmap smoothness to perform joint estimation using less than two full data sets and then we exploit the sensitivity encoding from parallel imaging to reduce the acquisition length and perform joint reconstruction using just one full k-space dataset. Subsequently, we extend the proposed method to account for the through-plane gradients of the field inhomogeneity. To improve the efficiency of the reconstruction algorithm we use a linearization technique for fieldmap estimation, which allows the use of the conjugate gradient algorithm. The resulting method allows for efficient reconstruction by applying fast approximations that allow the use of the conjugate gradient algorithm along with FFTs. Our proposed method can be computationally efficient for quadratic regularizers, but the CG-based algorithm is not directly applicable to non-quadratic regularization. To improve the efficiency of our method for non-quadratic regularization we propose an algorithm based on the augmented Lagrangian (AL) framework with variable splitting. This new algorithm can also be used for the non-linear optimization problem of fieldmap estimation without the need for the linearization approximation.PhDElectrical Engineering-SystemsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttp://deepblue.lib.umich.edu/bitstream/2027.42/102449/1/amatakos_1.pd

Magnetic resonance multitasking for motion-resolved quantitative cardiovascular imaging.

Quantitative cardiovascular magnetic resonance (CMR) imaging can be used to characterize fibrosis, oedema, ischaemia, inflammation and other disease conditions. However, the need to reduce artefacts arising from body motion through a combination of electrocardiography (ECG) control, respiration control, and contrast-weighting selection makes CMR exams lengthy. Here, we show that physiological motions and other dynamic processes can be conceptualized as multiple time dimensions that can be resolved via low-rank tensor imaging, allowing for motion-resolved quantitative imaging with up to four time dimensions. This continuous-acquisition approach, which we name cardiovascular MR multitasking, captures - rather than avoids - motion, relaxation and other dynamics to efficiently perform quantitative CMR without the use of ECG triggering or breath holds. We demonstrate that CMR multitasking allows for T1 mapping, T1-T2 mapping and time-resolved T1 mapping of myocardial perfusion without ECG information and/or in free-breathing conditions. CMR multitasking may provide a foundation for the development of setup-free CMR imaging for the quantitative evaluation of cardiovascular health

Three-dimensional echo-shifted EPI with simultaneous blip-up and blip-down acquisitions for correcting geometric distortion

Purpose: Echo-planar imaging (EPI) with blip-up/down acquisition (BUDA) can provide high-quality images with minimal distortions by using two readout trains with opposing phase-encoding gradients. Because of the need for two separate acquisitions, BUDA doubles the scan time and degrades the temporal resolution when compared to single-shot EPI, presenting a major challenge for many applications, particularly functional MRI (fMRI). This study aims at overcoming this challenge by developing an echo-shifted EPI BUDA (esEPI-BUDA) technique to acquire both blip-up and blip-down datasets in a single shot. Methods: A three-dimensional (3D) esEPI-BUDA pulse sequence was designed by using an echo-shifting strategy to produce two EPI readout trains. These readout trains produced a pair of k-space datasets whose k-space trajectories were interleaved with opposite phase-encoding gradient directions. The two k-space datasets were separately reconstructed using a 3D SENSE algorithm, from which time-resolved B0-field maps were derived using TOPUP in FSL and then input into a forward model of joint parallel imaging reconstruction to correct for geometric distortion. In addition, Hankel structured low-rank constraint was incorporated into the reconstruction framework to improve image quality by mitigating the phase errors between the two interleaved k-space datasets. Results: The 3D esEPI-BUDA technique was demonstrated in a phantom and an fMRI study on healthy human subjects. Geometric distortions were effectively corrected in both phantom and human brain images. In the fMRI study, the visual activation volumes and their BOLD responses were comparable to those from conventional 3D echo-planar images. Conclusion: The improved imaging efficiency and dynamic distortion correction capability afforded by 3D esEPI-BUDA are expected to benefit many EPI applications.Comment: 8 figures, peer-reviewed journal pape

Parallel Magnetic Resonance Imaging as Approximation in a Reproducing Kernel Hilbert Space

In Magnetic Resonance Imaging (MRI) data samples are collected in the spatial frequency domain (k-space), typically by time-consuming line-by-line scanning on a Cartesian grid. Scans can be accelerated by simultaneous acquisition of data using multiple receivers (parallel imaging), and by using more efficient non-Cartesian sampling schemes. As shown here, reconstruction from samples at arbitrary locations can be understood as approximation of vector-valued functions from the acquired samples and formulated using a Reproducing Kernel Hilbert Space (RKHS) with a matrix-valued kernel defined by the spatial sensitivities of the receive coils. This establishes a formal connection between approximation theory and parallel imaging. Theoretical tools from approximation theory can then be used to understand reconstruction in k-space and to extend the analysis of the effects of samples selection beyond the traditional g-factor noise analysis to both noise amplification and approximation errors. This is demonstrated with numerical examples.Comment: 28 pages, 7 figure