The HCP has received requests for slice-timing parameters for the HCP fMRI data and, therefore, we are providing this information with some background as to why the HCP functional preprocessing pipeline does not include slice-timing correction (see Smith et al. "Resting-state fMRI in the Human Connectome Project" NeuroImage, vol. 80, pp. 144-168, Oct. 2013 and  Glasser et al. “The minimal preprocessing pipelines for the Human Connectome Project” Neuroimage, vol. 80, pp. 105–124, Oct. 2013).

The slice timing of the HCP fMRI is non-trivial due to the multi-band nature of the acquisition (MB=8 was used on rfMRI and tfMRI). Examples of the timing for both an even or odd number of total slices, and differing MB factors are provided in the following file: CMRR_MB_Slice_Order.pdf   

The current version of the MB sequence is available here:

The TR for HCP fMRI is quite short at 0.72 sec. Because of the short TR, HCP consortium members felt the value of slice timing correction (STC) for HCP fMRI data is limited.  Therefore we have opted not to slice-time correct the task or resting state data in our HCP processing.

If you still feel it is necessary to do slice timing correction, you will need to account for the non-trivial nature of slice acquisition in the MB acquisitions. Although we aren't distributing DICOMs (due to their size and the presence of privileged information in them), we do have slice time order/acquisition information available in the original DICOM headers.

In the following text file, we have extracted the slice timing order info for an exemplar subject. Note that these values are only accurate to 2.5 ms:

BOLD acquisition times extracted from {Subject_ID}/REST1_RL, 2nd DICOM: MosaicRefAcqTimes.txt

The issue of when to apply STC in a processing stream is the subject of some debate.

AFNI and SPM apply STC before motion correction, whereas FSL applies it after.  The FSL rationale is that applying STC after motion correction (and relatedly, gradient distortion correction) is preferable because you don't want to "bake in" errors related to spatial misalignment as part of STC, since spatial intensity differences between voxels can be much larger (e.g., CSF vs. gray matter) than temporal intensity changes over a TR.

This would be even more true for the short TR data of the HCP.  So, if STC is applied to HCP fMRI data, it is likely best to apply it to the minimally preprocessed data (i.e., after gradient nonlinearity correct, motion correction, and B0 distortion correction).



  • No labels