FMRI: Difference between revisions
No edit summary |
|||
| Line 1: | Line 1: | ||
=====[[MEICA]] and [[Tedana]]===== | =====[[MEICA]] and [[Tedana]]===== | ||
==Motion-related artefacts in fMRI== | |||
Participant in-scanner motion is one of the prominent sources of artefacts in fMRI data. This was recognized as an issue in fMRI long ago (Hajnal et al., 1994; Friston et al., 1996). Motion has complex effects on the signal - including increases, decreases or complex waveforms, depending on factors such as the timing, duration and trajectory of motion (Power et al., 2015, Byrge and Kennedy 2018). Crucially, motion can substantially affect estimates of functional connectivity (Power et al., 2012; Satterthwaite et al., 2012). Initial motion-correction approaches involved regression of motion parameters and their derivatives from voxel-wise BOLD time series (Friston et al., 1996), or regression of the average or "global" signal (Aguirre et al., 1998). Later solutions included censoring of motion-affected frames (Power et al., 2012), removal of non-BOLD (artefactual) signal using independent component analysis applied to multi-echo fMRI data, based on the dependence of the BOLD signal on echo time (Kundu et al., 2012, 2013), or "despiking" of motion-related non-stationary events from a wavelet decomposition of the signal (Patel et al., 2014). The latter method was subsequently extended to provide revised estimates of voxel-wise effective degrees of freedom (df ) of the BOLD time series, which due to denoising are lower than the nominal N(df ) = N(time-points), and which affect estimates of edge probability when incorporated into network analysis (Patel and Bullmore, 2016). | |||
Other prominent artifacts, particularly in fMRI, relate to respiration, vasculature and arousal (e.g.: Murphy et al., 2013; Caballero-Gaudes and Reynolds, 2017). | |||
==Diagnostic measures== | ==Diagnostic measures== | ||
Revision as of 11:35, 10 May 2018
MEICA and Tedana
Participant in-scanner motion is one of the prominent sources of artefacts in fMRI data. This was recognized as an issue in fMRI long ago (Hajnal et al., 1994; Friston et al., 1996). Motion has complex effects on the signal - including increases, decreases or complex waveforms, depending on factors such as the timing, duration and trajectory of motion (Power et al., 2015, Byrge and Kennedy 2018). Crucially, motion can substantially affect estimates of functional connectivity (Power et al., 2012; Satterthwaite et al., 2012). Initial motion-correction approaches involved regression of motion parameters and their derivatives from voxel-wise BOLD time series (Friston et al., 1996), or regression of the average or "global" signal (Aguirre et al., 1998). Later solutions included censoring of motion-affected frames (Power et al., 2012), removal of non-BOLD (artefactual) signal using independent component analysis applied to multi-echo fMRI data, based on the dependence of the BOLD signal on echo time (Kundu et al., 2012, 2013), or "despiking" of motion-related non-stationary events from a wavelet decomposition of the signal (Patel et al., 2014). The latter method was subsequently extended to provide revised estimates of voxel-wise effective degrees of freedom (df ) of the BOLD time series, which due to denoising are lower than the nominal N(df ) = N(time-points), and which affect estimates of edge probability when incorporated into network analysis (Patel and Bullmore, 2016).
Other prominent artifacts, particularly in fMRI, relate to respiration, vasculature and arousal (e.g.: Murphy et al., 2013; Caballero-Gaudes and Reynolds, 2017).
Diagnostic measures
ME-ICA processing for NSPN and BioDep
The following pipeline describes processing steps for multi-echo ICA (ME-ICA) data, as agreed by Manfred, Frantisek and others in May 2018. This is the pipeline applied to NSPN and BioDep projects.
Briefly, the processing pipeline includes the following steps:
...
The specific ME-ICA script to be run is located on the BCNI server at: "/home/rr480/Code/MATLAB/fMRI/MEICA/meica_manfred_2018_05_10/meica.py". This includes the options listed below (some of which are specific to NSPN data). For a detailed description of all processing steps, see MEICA.
-d 'echo[1,2,3].nii'
-a 'brainmask.nii'
-e 13.0,30.55,48.1
--TR=2.43
--tpattern=seq-z
-b 15s
--fres=2.5
--mask_mode='anat'
--no_skullstrip
--OVERWRITE
