PEAR: PEriodic And fixed Rank separation for fast fMRI

Purpose: In functional MRI (fMRI), faster acquisition via undersampling of data can improve the spatial‐temporal resolution trade‐off and increase statistical robustness through increased degrees‐of‐freedom. High‐quality reconstruction of fMRI data from undersampled measurements requires proper modeling of the data. We present an fMRI reconstruction approach based on modeling the fMRI signal as a sum of periodic and fixed rank components, for improved reconstruction from undersampled measurements. Methods: The proposed approach decomposes the fMRI signal into a component which has a fixed rank and a component consisting of a sum of periodic signals which is sparse in the temporal Fourier domain. Data reconstruction is performed by solving a constrained problem that enforces a fixed, moderate rank on one of the components, and a limited number of temporal frequencies on the other. Our approach is coined PEAR ‐ PEriodic And fixed Rank separation for fast fMRI. Results: Experimental results include purely synthetic simulation, a simulation with real timecourses and retrospective undersampling of a real fMRI dataset. Evaluation was performed both quantitatively and visually versus ground truth, comparing PEAR to two additional recent methods for fMRI reconstruction from undersampled measurements. Results demonstrate PEAR's improvement in estimating the timecourses and activation maps versus the methods compared against at acceleration ratios of R = 8,10.66 (for simulated data) and R = 6.66,10 (for real data). Conclusions: This paper presents PEAR, an undersampled fMRI reconstruction approach based on decomposing the fMRI signal to periodic and fixed rank components. PEAR results in reconstruction with higher fidelity than when using a fixed‐rank based model or a conventional Low‐rank + Sparse algorithm. We have shown that splitting the functional information between the components leads to better modeling of fMRI, over state‐of‐the‐art methods.

[1]  Jong Chul Ye,et al.  Performance evaluation of accelerated functional MRI acquisition using compressed sensing , 2009, 2009 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[2]  Essa Yacoub,et al.  High-field fMRI unveils orientation columns in humans , 2008, Proceedings of the National Academy of Sciences.

[3]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[4]  Peter Boesiger,et al.  k‐t BLAST and k‐t SENSE: Dynamic MRI with high frame rate exploiting spatiotemporal correlations , 2003, Magnetic resonance in medicine.

[5]  Urs Niesen,et al.  Adaptive Alternating Minimization Algorithms , 2007, IEEE Transactions on Information Theory.

[6]  Mark W. Woolrich,et al.  FSL , 2012, NeuroImage.

[7]  Yonina C. Eldar,et al.  The SPURS Algorithm for Resampling an Irregularly Sampled Signal onto a Cartesian Grid , 2016, IEEE Transactions on Medical Imaging.

[8]  Olaf Dössel,et al.  An Optimal Radial Profile Order Based on the Golden Ratio for Time-Resolved MRI , 2007, IEEE Transactions on Medical Imaging.

[9]  Essa Yacoub,et al.  The WU-Minn Human Connectome Project: An overview , 2013, NeuroImage.

[10]  P. Boesiger,et al.  SENSE: Sensitivity encoding for fast MRI , 1999, Magnetic resonance in medicine.

[11]  Dieter Jaeger,et al.  Quasi-periodic patterns (QPP): Large-scale dynamics in resting state fMRI that correlate with local infraslow electrical activity , 2014, NeuroImage.

[12]  Peter Boesiger,et al.  Compressed sensing in dynamic MRI , 2008, Magnetic resonance in medicine.

[13]  Justin P. Haldar,et al.  Low rank matrix recovery for real-time cardiac MRI , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[14]  A J Sederman,et al.  Compressed sensing reconstruction improves sensitivity of variable density spiral fMRI , 2013, Magnetic resonance in medicine.

[15]  Stephen M. Smith,et al.  Temporally-independent functional modes of spontaneous brain activity , 2012, Proceedings of the National Academy of Sciences.

[16]  Jong Chul Ye,et al.  k‐t FOCUSS: A general compressed sensing framework for high resolution dynamic MRI , 2009, Magnetic resonance in medicine.

[17]  T. Blumensath,et al.  Theory and Applications , 2011 .

[18]  Daniel K Sodickson,et al.  Low‐rank plus sparse matrix decomposition for accelerated dynamic MRI with separation of background and dynamic components , 2015, Magnetic resonance in medicine.

[19]  Lewis Saks,et al.  Sense as in Sensitivity. , 1969 .

[20]  Michael Herbst,et al.  Improving temporal resolution in fMRI using a 3D spiral acquisition and low rank plus sparse (L+S) reconstruction , 2017, NeuroImage.

[21]  Stephen M Smith,et al.  k-t FASTER: Acceleration of functional MRI data acquisition using low rank constraints , 2014, Magnetic resonance in medicine.

[22]  Stephen M. Smith,et al.  Investigations into resting-state connectivity using independent component analysis , 2005, Philosophical Transactions of the Royal Society B: Biological Sciences.

[23]  Yonina C. Eldar,et al.  Fast reference based MRI , 2015, 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[24]  Priya Aggarwal,et al.  Optshrink LR + S: accelerated fMRI reconstruction using non-convex optimal singular value shrinkage , 2017, Brain Informatics.

[25]  Peter J. Koopmans,et al.  Multi-echo fMRI of the cortical laminae in humans at 7T , 2011, NeuroImage.

[26]  D. Donoho,et al.  The Optimal Hard Threshold for Singular Values is 4 / √ 3 , 2013 .

[27]  I. Daubechies,et al.  An iterative thresholding algorithm for linear inverse problems with a sparsity constraint , 2003, math/0307152.

[28]  M Desco,et al.  Exploitation of temporal redundancy in compressed sensing reconstruction of fMRI studies with a prior-based algorithm (PICCS). , 2015, Medical physics.

[29]  ManKin Choy,et al.  High spatial resolution compressed sensing (HSPARSE) functional MRI , 2016, Magnetic resonance in medicine.

[30]  David L. Donoho,et al.  The Optimal Hard Threshold for Singular Values is 4/sqrt(3) , 2013, 1305.5870.

[31]  Stephen M. Smith,et al.  Multiplexed Echo Planar Imaging for Sub-Second Whole Brain FMRI and Fast Diffusion Imaging , 2010, PloS one.

[32]  Yonina C. Eldar,et al.  Compressed sensing for longitudinal MRI: An adaptive-weighted approach. , 2014, Medical physics.

[33]  E.J. Candes Compressive Sampling , 2022 .

[34]  Jeffrey A. Fessler,et al.  Nonuniform fast Fourier transforms using min-max interpolation , 2003, IEEE Trans. Signal Process..

[35]  N J Pelc,et al.  Unaliasing by Fourier‐encoding the overlaps using the temporal dimension (UNFOLD), applied to cardiac imaging and fMRI , 1999, Magnetic resonance in medicine.

[36]  D. Donoho,et al.  Sparse MRI: The application of compressed sensing for rapid MR imaging , 2007, Magnetic resonance in medicine.

[37]  Yonina C. Eldar,et al.  Reference-based MRI. , 2015, Medical physics.

[38]  Stephen M. Smith,et al.  Accelerating functional MRI using fixed‐rank approximations and radial‐cartesian sampling , 2016, Magnetic resonance in medicine.

[39]  N. Filippini,et al.  Group comparison of resting-state FMRI data using multi-subject ICA and dual regression , 2009, NeuroImage.

[40]  Emmanuel J. Candès,et al.  A Singular Value Thresholding Algorithm for Matrix Completion , 2008, SIAM J. Optim..

[41]  Jong Chul Ye,et al.  Compressed sensing fMRI using gradient-recalled echo and EPI sequences , 2014, NeuroImage.

[42]  Priya Aggarwal,et al.  Accelerated fMRI reconstruction using Matrix Completion with Sparse Recovery via Split Bregman , 2016, Neurocomputing.

[43]  K. Schittkowski,et al.  NONLINEAR PROGRAMMING , 2022 .

[44]  S. Frick,et al.  Compressed Sensing , 2014, Computer Vision, A Reference Guide.

[45]  Yi Ma,et al.  Robust principal component analysis? , 2009, JACM.

[46]  Tülay Adali,et al.  Comparison of multi‐subject ICA methods for analysis of fMRI data , 2010, Human brain mapping.

[47]  Vince D Calhoun,et al.  Optimal compressed sensing reconstructions of fMRI using 2D deterministic and stochastic sampling geometries , 2012, BioMedical Engineering OnLine.

[48]  Patrick L. Combettes,et al.  Signal Recovery by Proximal Forward-Backward Splitting , 2005, Multiscale Model. Simul..

[49]  Shiqian Ma,et al.  Convergence of Fixed-Point Continuation Algorithms for Matrix Rank Minimization , 2009, Found. Comput. Math..

[50]  Ahmed H. Tewfik,et al.  Under-sampled functional MRI using low-rank plus sparse matrix decomposition , 2015, 2015 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[51]  Mark Chiew,et al.  Recovering task fMRI signals from highly under-sampled data with low-rank and temporal subspace constraints , 2018, NeuroImage.

[52]  Jennifer A McNab,et al.  Motion correction for functional MRI with three‐dimensional hybrid radial‐Cartesian EPI , 2016, Magnetic resonance in medicine.

[53]  N. Schuff,et al.  Accelerated fMRI using Low-Rank Model and Sparsity Constraints , 2012 .