Detecting spontaneous brain activity in functional magnetic resonance imaging using finite rate of innovation

Several methods have been developed for the sampling and reconstruction of specific classes of signals known as signals with finite rate of innovation (FRI). It is possible to recover the innovations of the signals from very low-rate samples by using adequate exponential reproduction sampling kernels. Recently, the FRI theory has been extended to arbitrary sampling kernels that reproduce approximate exponentials. In this paper, we develop the method for the detection of spontaneous brain activity in functional magnetic resonance imaging (fMRI) data. We model the fMRI timecourse for every voxel as a convolution between the innovation signal - a stream of Diracs- and the hemodynamic response function (HRF). Relaxing the exact exponential reproduction constraint given by Strang-Fix condition, we design an adequate FRI sampling kernel using the canonical HRF model that allows us to retrieve the innovation instants in continuous domain. We illustrate the feasibility of our method by detecting spontaneous brain activity on the simulated and degraded fMRI data using an iterative denoising scheme.

[1]  Karl J. Friston,et al.  Analysis of fMRI Time-Series Revisited—Again , 1995, NeuroImage.

[2]  Bharat B. Biswal,et al.  Resting state fMRI: A personal history , 2012, NeuroImage.

[3]  Natalia Petridou,et al.  Paradigm free mapping with sparse regression automatically detects single‐trial functional magnetic resonance imaging blood oxygenation level dependent responses , 2011, Human brain mapping.

[4]  Thierry Blu,et al.  Extrapolation and Interpolation) , 2022 .

[5]  M. Vetterli,et al.  Sparse Sampling of Signal Innovations , 2008, IEEE Signal Processing Magazine.

[6]  Karl J. Friston,et al.  Statistical parametric mapping , 2013 .

[7]  Thierry Blu,et al.  FRI Sampling With Arbitrary Kernels , 2013, IEEE Transactions on Signal Processing.

[8]  Arie Feuer,et al.  On perfect conditioning of Vandermonde matrices on the unit circle , 2007 .

[9]  Dimitri Van De Ville,et al.  Total activation: fMRI deconvolution through spatio-temporal regularization , 2013, NeuroImage.

[10]  Thierry Blu,et al.  Cardinal exponential splines: part I - theory and filtering algorithms , 2005, IEEE Transactions on Signal Processing.

[11]  Mohamed-Jalal Fadili,et al.  Activelets: Wavelets for sparse representation of hemodynamic responses , 2011, Signal Process..

[12]  Thierry Blu,et al.  Sampling signals with finite rate of innovation , 2002, IEEE Trans. Signal Process..

[13]  Karl J. Friston,et al.  Analysis of fMRI Time-Series Revisited , 1995, NeuroImage.

[14]  Thierry Blu,et al.  Approximate Strang-Fix: sampling infinite streams of Diracs with any kernel , 2013, Optics & Photonics - Optical Engineering + Applications.

[15]  Andrea Bergmann,et al.  Statistical Parametric Mapping The Analysis Of Functional Brain Images , 2016 .

[16]  Pier Luigi Dragotti,et al.  A finite rate of innovation algorithm for fast and accurate spike detection from two-photon calcium imaging , 2013, Journal of neural engineering.