EXTRACTING THE HAEMODYNAMIC RESPONSE FUNCTION FROM FMRI TIME SERIES USING FOURIER-WAVELET REGULARISED DECONVOLUTION WITH ORTHOGONAL SPLINE WAVELETS

We describe a method to extract the haemodynamic response func-tion (HRF) from functional magnetic resonance imaging (fMRI) time series based on Fourier-wavelet regularised deconvolution (ForWaRD), and introduce a simple model for the HRF. The HRF extraction algorithm extends the ForWaRD algorithm by introduc-ing a more efficient computation in the case of very long wavelet filters. We compute shift-invariant discrete wavelet transforms (SI-DWT) in the frequency domain, and apply ForWaRD using orthog-onal spline wavelets. Extraction and modelling of subject-specific HRFs is demonstrated, as well as the use of these HRFs in a sub-sequent brain activation analysis. Temporal responses are predicted by using the extracted HRF coefficients. The resulting activation maps show the effectiveness of the proposed method. 1