Dual-TRACER: High resolution fMRI with constrained evolution reconstruction

Subspace-constrained approaches to low-rank fMRI acceleration

Acceleration methods in fMRI aim to reconstruct high fidelity images from under-sampled k-space, allowing fMRI datasets to achieve higher temporal resolution, reduced physiological noise aliasing, and increased statistical degrees of freedom. While low levels of acceleration are typically part of standard fMRI protocols through parallel imaging, there exists the potential for approaches that allow much greater acceleration. One such existing approach is k-t FASTER, which exploits the inherent low-rank nature of fMRI. In this paper, we present a reformulated version of k-t FASTER which includes additional L2 constraints within a low-rank framework.

We evaluated the effect of three different constraints against existing low-rank approaches to fMRI reconstruction: Tikhonov constraints, low-resolution priors, and temporal subspace smoothness. The different approaches are separately tested for robustness to under-sampling and thermal noise levels, in both retrospectively and prospectively-undersampled finger-tapping task fMRI data. Reconstruction quality is evaluated by accurate reconstruction of low-rank subspaces and activation maps.

The use of L2 constraints was found to achieve consistently improved results, producing high fidelity reconstructions of statistical parameter maps at higher acceleration factors and lower SNR values than existing methods, but at a cost of longer computation time. In particular, the Tikhonov constraint proved very robust across all tested datasets, and the temporal subspace smoothness constraint provided the best reconstruction scores in the prospectively-undersampled dataset. These results demonstrate that regularized low-rank reconstruction of fMRI data can recover functional information at high acceleration factors without the use of any model-based spatial constraints.

Improved statistical efficiency of simultaneous multi-slice fMRI by reconstruction with spatially adaptive temporal smoothing

We introduce an approach to reconstruction of simultaneous multi-slice (SMS)-fMRI data that improves statistical efficiency. The method incorporates regularization to adjust temporal smoothness in a spatially varying, encoding-dependent manner, reducing the g-factor noise amplification per temporal degree of freedom. This results in a net improvement in tSNR and GLM efficiency, where the efficiency gain can be derived analytically as a function of the encoding and reconstruction parameters. Residual slice leakage and aliasing is limited when fMRI signal energy is dominated by low frequencies. Analytical predictions, simulated and experimental results demonstrate a marked improvement in statistical efficiency in the temporally regularized reconstructions compared to conventional slice-GRAPPA reconstructions, particularly in central brain regions. Furthermore, experimental results confirm that residual slice leakage and aliasing errors are not noticeably increased compared to slice-GRAPPA reconstruction. This approach to temporally regularized image reconstruction in SMS-fMRI improves statistical power, and allows for explicit choice of reconstruction parameters by directly assessing their impact on noise variance per degree of freedom.