open-discussion
open-discussion > RE: SLOMOCO and ART Toolbox
May 21, 2015 05:05 PM | Erik Beall
RE: SLOMOCO and ART Toolbox
Hi Eduardo,
Yes, that is the correct file for the PESTICA+SLOMOCO regressed EPI. Ignore this if you've already done this, but if not, check the timeseries to make sure cardiac fluctuations experience more reduction than non-cardiac fluctuations. I do this by looking at the timeseries in an artery before and after correction in AFNI with the "Plot" viewer up and switching back and forth between raw and corrected - after correction, everything will have reduced variance and this is true even if you're regressing pure random noise, but if you're regressing cardiac, then the spikier stuff in arteries should experience more reduction than adjacent voxel timeseries (note, there will still be spikey timeseries stuff in arteries even after the best existing physio corrections). After SLOMOCO, you shouldn't have to do any more regression for motion purposes.
If you have to do motion censoring, I do indeed recommend using the slomoco.TDzmetric.txt file instead of a volumetric file like the FSL_motion_outliers. One way you can do it is to take the maximum over each slice group (like volumetric_max_head_motion=squeeze(max(reshape(slomoco,[zdim tdim]))); - I might have the zdim and tdim backwards, I don't have a linux box with matlab available at the moment so I can't check) and then censor by volume with existing thresholds.
Erik
PS: will you be at ISMRM in Toronto in two weeks? I have a poster on SLOMOCO and I'd be happy to talk in more detail in person then.
Yes, that is the correct file for the PESTICA+SLOMOCO regressed EPI. Ignore this if you've already done this, but if not, check the timeseries to make sure cardiac fluctuations experience more reduction than non-cardiac fluctuations. I do this by looking at the timeseries in an artery before and after correction in AFNI with the "Plot" viewer up and switching back and forth between raw and corrected - after correction, everything will have reduced variance and this is true even if you're regressing pure random noise, but if you're regressing cardiac, then the spikier stuff in arteries should experience more reduction than adjacent voxel timeseries (note, there will still be spikey timeseries stuff in arteries even after the best existing physio corrections). After SLOMOCO, you shouldn't have to do any more regression for motion purposes.
If you have to do motion censoring, I do indeed recommend using the slomoco.TDzmetric.txt file instead of a volumetric file like the FSL_motion_outliers. One way you can do it is to take the maximum over each slice group (like volumetric_max_head_motion=squeeze(max(reshape(slomoco,[zdim tdim]))); - I might have the zdim and tdim backwards, I don't have a linux box with matlab available at the moment so I can't check) and then censor by volume with existing thresholds.
Erik
PS: will you be at ISMRM in Toronto in two weeks? I have a poster on SLOMOCO and I'd be happy to talk in more detail in person then.
Threaded View
Title | Author | Date |
---|---|---|
Eduardo A Garza Villarreal | May 4, 2015 | |
Erik Beall | May 7, 2015 | |
Eduardo A Garza Villarreal | May 12, 2015 | |
Erik Beall | May 14, 2015 | |
Eduardo A Garza Villarreal | May 18, 2015 | |
Erik Beall | May 21, 2015 | |
Eduardo A Garza Villarreal | Dec 8, 2015 | |
Eduardo A Garza Villarreal | Dec 9, 2015 | |
Eduardo A Garza Villarreal | May 15, 2015 | |
Eduardo A Garza Villarreal | May 7, 2015 | |