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

Improved l1-SPIRiT using 3D walsh transform-based sparsity basis

l1-SPIRiT is a fast magnetic resonance imaging (MRI) method which combines parallel imaging (PI) with compressed sensing (CS) by performing a joint l1-norm and l2-norm optimization procedure. The original l1-SPIRiT method uses two-dimensional (2D) Wavelet transform to exploit the intra-coil data redundancies and a joint sparsity model to exploit the inter-coil data redundancies. In this work, we propose to stack all the coil images into a three-dimensional (3D) matrix, and then a novel 3D Walsh transform-based sparsity basis is applied to simultaneously reduce the intra-coil and inter-coil data redundancies. Both the 2D Wavelet transform-based and the proposed 3D Walsh transform-based sparsity bases were investigated in the l1-SPIRiT method. The experimental results show that the proposed 3D Walsh transform-based l1-SPIRiT method outperformed the original l1-SPIRiT in terms of image quality and computational efficiency

Accelerating Magnetic Resonance Parametric Mapping Using Simultaneously Spatial Patch-based and Parametric Group-based Low-rank Tensors (SMART)

Quantitative magnetic resonance (MR) parametric mapping is a promising approach for characterizing intrinsic tissue-dependent information. However, long scan time significantly hinders its widespread applications. Recently, low-rank tensor has been employed and demonstrated good performance in accelerating MR parametricmapping. In this study, we propose a novel method that uses spatial patch-based and parametric group-based low rank tensors simultaneously (SMART) to reconstruct images from highly undersampled k-space data. The spatial patch-based low-rank tensor exploits the high local and nonlocal redundancies and similarities between the contrast images in parametric mapping. The parametric group based low-rank tensor, which integrates similar exponential behavior of the image signals, is jointly used to enforce the multidimensional low-rankness in the reconstruction process. In vivo brain datasets were used to demonstrate the validity of the proposed method. Experimental results have demonstrated that the proposed method achieves 11.7-fold and 13.21-fold accelerations in two-dimensional and three-dimensional acquisitions, respectively, with more accurate reconstructed images and maps than several state-of-the-art methods. Prospective reconstruction results further demonstrate the capability of the SMART method in accelerating MR quantitative imaging.Comment: 15 pages, 12 figure

Structured low-rank methods for robust 3D multi-shot EPI

Magnetic resonance imaging (MRI) has inherently slow acquisition speed, and Echo-Planar Imaging (EPI), as an efficient acquisition scheme, has been widely used in functional magnetic resonance imaging (fMRI) where an image series with high temporal resolution is needed to measure neuronal activity. Recently, 3D multi-shot EPI which samples data from an entire 3D volume with repeated shots has been drawing growing interest for fMRI with its high isotropic spatial resolution, particularly at ultra-high fields. However, compared to single-shot EPI, multi-shot EPI is sensitive to any inter-shot instabilities, e.g., subject movement and even physiologically induced field fluctuations. These inter-shot inconsistencies can greatly negate the theoretical benefits of 3D multi-shot EPI over conventional 2D multi-slice acquisitions. Structured low-rank image reconstruction which regularises under-sampled image reconstruction by exploiting the linear dependencies in MRI data has been successfully demonstrated in a variety of applications. In this thesis, a structured low-rank reconstruction method is optimised for 3D multi-shot EPI imaging together with a dedicated sampling pattern termed seg-CAIPI, in order to enhance the robustness to physiological fluctuations and improve the temporal stability of 3D multi-shot EPI for fMRI at 7T. Moreover, a motion compensated structured low-rank reconstruction framework is also presented for robust 3D multi-shot EPI which further takes into account inter-shot instabilities due to bulk motion. Lastly, this thesis also investigates into the improvement of structured low-rank reconstruction from an algorithmic perspective and presents the locally structured low-rank reconstruction scheme

Calibrationless Multi-coil Magnetic Resonance Imaging with Compressed Sensing

We present a method for combining the data retrieved by multiple coils of a Magnetic Resonance Imaging (MRI) system with the a priori assumption of compressed sensing to reconstruct a single image. The final image is the result of an optimization problem that only includes constraints based on fundamental physics (Maxwell's equations and the Biot-Savart law) and accepted phenomena (e.g. sparsity in the Wavelet domain). The problem is solved using an alternating minimization approach: two convex optimization problems are alternately solved, one with the Fast Iterative Shrinkage Threshold Algorithm (FISTA) and the other with the Primal-Dual Hybrid Gradient (PDHG) method. We show results on simulated data as well as data of the knee, brain, and ankle. In all cases studied, results from the new algorithm show higher quality and increased detail when compared to conventional reconstruction algorithms