Hello FSL Experts,
I was hoping to get some feed back on my attempt to denoise a set of 4D files using melodic and fsl_regfilt. I am trying to remove any melodic component that has a high correlation with motion. So I followed the following procedure:
1. Ran Single session-ICA melodic on all 4D images
2. Ran a script to run mcflirt on all 4D images and the nextract time courses for each of the 6 motion components (xyz rot/trans) from the mcflirt par file.
3. Used a MATLAB script to run the “regress” matlab function where the motion parts are run against each melodic component. This is repeated for all images/all melodic components.
4. Collected any melodic component with a p-value under .01 for each image
5. Ran fsl_regfilt removing the components found in 4
I am hoping for any feedback as to whether I am doing this properly or if there is a better method to achieve the same effect. I have read that melodic denoising is the preferred method for removing motion components rather than including the 6 motion components from step 2 into the first level GLM.
Thanks,
Michael
|