title: Data-driven haemodynamic response function extraction using Fourier-wavelet regularised deconvolution.
abstract: BACKGROUND
  We present a simple, data-driven method to extract haemodynamic response functions (HRF) from functional magnetic resonance imaging (fMRI) time series, based on the Fourier-wavelet regularised deconvolution (ForWaRD) technique. HRF data are required for many fMRI applications, such as defining region-specific HRFs, effciently representing a general HRF, or comparing subject-specific HRFs.
  
  
  RESULTS
  ForWaRD is applied to fMRI time signals, after removing low-frequency trends by a wavelet-based method, and the output of ForWaRD is a time series of volumes, containing the HRF in each voxel. Compared to more complex methods, this extraction algorithm requires few assumptions (separability of signal and noise in the frequency and wavelet domains and the general linear model) and it is fast (HRF extraction from a single fMRI data set takes about the same time as spatial resampling). The extraction method is tested on simulated event-related activation signals, contaminated with noise from a time series of real MRI images. An application for HRF data is demonstrated in a simple event-related experiment: data are extracted from a region with significant effects of interest in a first time series. A continuous-time HRF is obtained by fitting a nonlinear function to the discrete HRF coeffcients, and is then used to analyse a later time series.
  
  
  CONCLUSIONS
  With the parameters used in this paper, the extraction method presented here is very robust to changes in signal properties. Comparison of analyses with fitted HRFs and with a canonical HRF shows that a subject-specific, regional HRF significantly improves detection power. Sensitivity and specificity increase not only in the region from which the HRFs are extracted, but also in other regions of interest.
authors.affiliation: Robert Steiner MR Unit, Imaging Sciences Department Imperial College, and MRC Clinical Sciences Centre, Hammersmith Hospital, London, UK. amwink@imperial.ac.uk
authors.forename: Alle Meije
authors.initials: AM
authors.lastname: Wink
authors.forename: Hans
authors.initials: H
authors.lastname: Hoogduin
authors.forename: Jos B T M
authors.initials: JB
authors.lastname: Roerdink
mesh_terms: D001783:Blood Flow Velocity / Q000502:physiology*; D001921:Brain / Q000098:blood supply* / Q000502:physiology*; D002560:Cerebrovascular Circulation / Q000502:physiology*; D005583:Fourier Analysis; D006801:Humans; D007090:Image Interpretation, Computer-Assisted / Q000379:methods*; D008279:Magnetic Resonance Imaging / Q000379:methods*; D010100:Oxygen / Q000097:blood*; D010101:Oxygen Consumption / Q000502:physiology
chemical_list: D010100:Oxygen