Motion robust reconstruction of multi-shot diffusion-weighted images without phase estimation via locally low-rank regularization
Multi-shot diffusion-weighted magnetic resonance imaging acquires multiple k-space segments of diffusion-weighted MRI data, estimates reconstructed multi-shot diffusion weighted images, and combines the estimated images to obtain a final reconstructed MRI image. The estimation of images iteratively calculates updated multi-shot images from the multiple k-space segments and current multi-shot images using a convex model without estimating motion-induced phase, constructs multiple locally low-rank spatial-shot matrices from the updated multi-shot images, and calculates current multi-shot images from spatial-shot matrices.
1. A method for multi-shot diffusion-weighted magnetic resonance imaging, the method comprising:
a) acquiring with a magnetic resonance imaging system multiple k-space segments of diffusion-weighted MRI data;
b) estimating reconstructed multi-shot diffusion weighted images by iteratively performing steps comprising:
i) calculating updated multi-shot images from the multiple k-space segments and current multi-shot images using a convex model without estimating motion-induced phase;
ii) constructing multiple locally low-rank spatial-shot matrices from the updated multi-shot images;
iii) calculating the current multi-shot images from the spatial-shot matrices;
c) combining the reconstructed multi-shot diffusion weighted images to obtain a reconstructed MRI image.
2. The method of claim 1 wherein the diffusion-weighted MRI data is acquired in response to an echo planar imaging sequence after applying diffusion encoding gradients.
3. The method of claim 2 wherein the echo planar imaging sequence is a 2D single-refocused Stejskal-Tanner spin-echo planar imaging sequence.
4. The method of claim 1 wherein constructing the spatial-shot matrices from the updated multi-shot image comprises calculating a low-rank approximation by doing singular value decomposition (SVD) and soft-thresholding.
5. The method of claim 1 where each of the spatial-shot matrices corresponds to a spatial block in an imaging plane, where each column of the matrix corresponds to a different shot of the multiple k-space segments.
6. The method of claim 1 where calculating updated multi-shot images uses a sampling operator, a Fourier transform, and an encoding operator constructed from sensitivity maps, and their conjugate operators, and the acquired data.
7. The method of claim 6 where the sensitivity maps are estimated from T2 images acquired without applying diffusion encoding gradients.
8. The method of claim 1 further comprising flipping and conjugating the multiple k-space segments of diffusion-weighted MRI data and treating them as virtual shots.
9. The method of claim 8 further comprising estimating sensitivity maps for the virtual shots from conjugated T2 images.