Spatio-temporal wavelet regularization for parallel MRI reconstruction: application to functional MRI

Abstract

Background Parallel magnetic resonance imaging (MRI) is a fast imaging technique that helps acquiring highly resolved images in space/time. Its performance depends on the reconstruction algorithm, which can proceed either in the k-space or in the image domain. Objective and methods To improve the performance of the widely used SENSE algorithm, 2D regularization in the wavelet domain has been investigated. In this paper, we first extend this approach to 3D-wavelet representations and the 3D sparsity-promoting regularization term, in order to address reconstruction artifacts that propagate across adjacent slices. The resulting optimality criterion is convex but nonsmooth, and we resort to the parallel proximal algorithm to minimize it. Second, to account for temporal correlation between successive scans in functional MRI (fMRI), we extend our first contribution to 3D ? t acquisition schemes by incorporating a prior along the time axis into the objective function. Results Our first method (3D-UWR-SENSE) is validated on T1-MRI anatomical data for gray/white matter segmentation. The second method (4D-UWR-SENSE) is validated for detecting evoked activity during a fast eventrelated functional MRI protocol. Conclusion We show that our algorithm outperforms the SENSE reconstruction at the subject and group levels (15 subjects) for different contrasts of interest (motor or computation tasks) and two parallel acceleration factors (R ¼ 2 and R ¼ 4) on 2 2 3 mm3 echo planar imaging (EPI) images.