A GPU-enhanced workflow for non-Fourier SENSE reconstruction

TL;DR

This paper presents a GPU-accelerated workflow for non-Fourier SENSE MRI reconstruction, enabling fast and accurate image reconstruction in challenging scenarios with complex signal models.

Contribution

It introduces a GPU-compatible implementation of non-Fourier SENSE reconstruction with practical analysis, available as open-source code.

Findings

GPU implementation significantly speeds up reconstruction times

Proper stopping criteria are essential for optimal image quality

Workflow accurately computes coil sensitivities and off-resonance maps

Abstract

Purpose: Image reconstruction in challenging scenarios requires accurate characterisations of coil sensitivity profiles, local off-resonances (B0) and effective encoding fields. Reconstruction methods utilising all of this information rely on signal models that are not compatible with the classical Fourier/k-space interpretation of the coil data. Hence, the FFT and related techniques are no more applicable, rendering image reconstruction computationally demanding. Methods: This article contains a workflow for accurate sensitivity and B0 mapping as well as other required processing steps. An implementation of non-Fourier SENSE reconstruction is provide that is well suited for execution on a GPU using the FFT. Important practical aspects like stopping criteria and sources of image artifacts are analyzed and documented. Results: Highly performant image reconstruction could be demonstrated on a 2D and 3D spiral dataset. These datasets contain trajectories featuring readout durations up to 71.5ms and undersampling factors up to R = 7. Running the reconstruction on a GPU greatly boosts reconstruction speed. Stopping the reconstruction at the right moment is crucial for image quality. All methods included in this article are available in a public code repository. Conclusion: The provided implementation of non-Fourier SENSE reconstruction is highly performant. When it is executed on GPU, runtimes reach a duration feasible in practice. The presented workflow ensures robust and accurate computation of coil sensitive profiles and off-resonance maps.