Hi all, I’m hoping to preprocess a dataset using a combination of fmriprep, tedana, and xcp_d. It is my understanding that there’s plans for these software packages to become better integrated, but for now things are a bit disjointed. I think I’ve figured out what I need to do, so I’m mostly hoping that someone can confirm if this is right (or let me know if there’s a much easier way and/or if someone has a script I can adopt).
Preprocessing steps:
- Run fmriprep using the --me-output-echos argument
- Run tedana on the fmriprep outputs
- Take the optimally combined (but not denoised) tedana output, and warp it to a standard space using an fmriprep transformation matrix (perhaps using the script on this page, under “Warping scanner-space fMRIPrep outputs to standard space”).
- Take the tedana ICA mixing tsv file, and rename the components that are classified as signal (from the metrics tsv file) such that they start with
signal__
- Run xcp_d with the -c argument, specifying the edited ICA mixing file (to regress out the rejected ICs, but make sure they’re orthogonalized with respect to the signal regressors)
- Rejoice
Obviously some quality checking can occur in between steps to make sure the outputs are all sensible.
Hopefully I haven’t completely misunderstood everything! Let me know if anyone has any thoughts.
Thanks!
-Kirk