Improved MUSSELS reconstruction for high-resolution multi-shot diffusion weighted imaging

MUSSELS is a one-step iterative reconstruction method for multishot diffusion weighted (msDW) imaging. The current work presents an efficient implementation, termed IRLS MUSSELS, that enables faster reconstruction to enhance its utility for high-resolution diffusion MRI studies.

The recently proposed MUSSELS reconstruction belongs to a new class of parallel imaging-based methods that recover artifact-free DWIs from msDW data without needing phase compensation. The reconstruction is achieved via structured low-rank matrix completion algorithms, which are computationally demanding due to the large size of the Hankel matrices and their associated computations involving singular value decompositions. Because of this, computational demands of the MUSSELS reconstruction scales as the matrix size and the number of shots increases, which hinders its practical utility for high-resolution applications. In this work, we derive a computationally efficient MUSSELS formulation by modifying the iterative reweighted least squares (IRLS) method that were proposed earlier to solve such problems. Using whole-brain in vivo data, we show the utility of the IRLS MUSSELS for routine high-resolution studies with reduced computational burden.

IRLS MUSSELS provides about five times faster reconstruction for matrix sizes 192 × 192 and 256 × 256 compared to the earlier MUSSELS implementation. The widely employed conjugate symmetry priors can also be incorporated into IRLS MUSSELS to reduce blurring of the partial Fourier acquisitions, without incurring much computational burden.

The proposed method is observed to be computationally efficient to enable routine high-resolution studies. The computational complexity matches the traditional msDWI reconstruction methods and provides improved reconstruction results with the additional constraints.