A Modified Generalized Series Approach: Application to Sparsely Sampled fMRI

In functional MRI, it is often desirable to reduce the readout duration to make the acquired data less prone to T2^* susceptibility artifacts. In addition, a shorter readout length allows for a shorter minimum TE, which is important for optimizing SNR. This can be achieved by undersampling the k-space. However, the conventional Fourier transform-based reconstruction method suffers from under-sampling artifacts such as high-frequency ringing and loss of resolution. To address this problem, we revisit the constrained-model approach using the generalized-series (GS) which has been proposed to address the undersampling problem for dynamic MRI. We propose a modification to the conventional use of the model in order to reflect small hemodynamic signal changes typical in fMRI. Specifically, while realizing that having high model order is necessary to capture missing information, we found that it is not necessary to span all frequencies of GS basis functions uniformly. Instead, having k-space and GS “sampling” trajectories covering low-frequencies uniformly while spanning high-frequencies sparsely, was observed to be an efficient strategy. The ability of the method over the conventional GS approach in improving resolution of functional images and activation maps while reducing undersampling ringing is demonstrated by simulations and experiments at 3T. Reduction in the readout time allowed an increase of statistical signal power as compared to the fully sampled acquisition. Unlike compressed sensing approaches, the proposed method is linear and hence has lower computational complexity. The method could prove useful for other imaging modalities where the signal change is smaller than the baseline component.