Skip to main content

BIDs App to retrieve the haemodynamic response function from resting state fMRI data

Project description

Resting state HRF estimation and deconvolution.

PyPI version

Please refer to https://github.com/compneuro-da/rsHRF for MATLAB version

BOLD HRF

The basic idea

This toolbox is aimed to retrieve the onsets of pseudo-events triggering an hemodynamic response from resting state fMRI BOLD voxel-wise signal. It is based on point process theory, and fits a model to retrieve the optimal lag between the events and the HRF onset, as well as the HRF shape, using either the canonical shape with two derivatives, or a (smoothed) Finite Impulse Response.

BOLD HRF

Once that the HRF has been retrieved for each voxel, it can be deconvolved from the time series (for example to improve lag-based connectivity estimates), or one can map the shape parameters everywhere in the brain (including white matter), and use the shape as a pathophysiological indicator.

HRF map

How to use the toolbox

The input is voxelwise BOLD signal, already preprocessed according to your favorite recipe. Important thing are:

  • bandpass filter in the 0.01-0.08 Hz interval (or something like that)
  • z-score the voxel BOLD time series

To be on the safe side, these steps are performed again in the code.

The input can be images (3D or 4D), or directly matrices of [observation x voxels].

It is possible to use a temporal mask to exclude some time points (for example after scrubbing).

The demos allow you to run the analyses on several formats of input data.

Python Package and BIDS-app

A BIDS-App has been made for easy and reproducible analysis. Its documentation can be accessed at:

http://bids-apps.neuroimaging.io/rsHRF/

Collaborators

  • Guorong Wu

  • Nigel Colenbier

  • Sofie Van Den Bossche

  • Daniele Marinazzo

  • Madhur Tandon (Python - BIDS)

  • Asier Erramuzpe (Python - BIDS)

References

  1. Guo-Rong Wu, Wei Liao, Sebastiano Stramaglia, Ju-Rong Ding, Huafu Chen, Daniele Marinazzo*. "A blind deconvolution approach to recover effective connectivity brain networks from resting state fMRI data." Medical Image Analysis, 2013, 17:365-374. PDF

  2. Guo-Rong Wu, Daniele Marinazzo. "Sensitivity of the resting state hemodynamic response function estimation to autonomic nervous system fluctuations." Philosophical Transactions of the Royal Society A, 2016, 374: 20150190. PDF

  3. Guo-Rong Wu, Daniele Marinazzo. "Retrieving the Hemodynamic Response Function in resting state fMRI: methodology and applications." PeerJ PrePrints, 2015. PDF

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for rsHRF, version 1.0.1
Filename, size File type Python version Upload date Hashes
Filename, size rsHRF-1.0.1-py3-none-any.whl (19.5 kB) File type Wheel Python version py3 Upload date Hashes View hashes
Filename, size rsHRF-1.0.1.tar.gz (17.1 kB) File type Source Python version None Upload date Hashes View hashes

Supported by

Elastic Elastic Search Pingdom Pingdom Monitoring Google Google BigQuery Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN SignalFx SignalFx Supporter DigiCert DigiCert EV certificate StatusPage StatusPage Status page