Hi Pual,
Thank you. I didn't know 3dcalc can be done with roi label directly until I saw your answer.
And this is an adaptation from your answer to meet my need:
# get a list of all roi names( get labeltable|rm summary lines|rm "|print 2ed column names)
set rois = `3dinfo -labeltable mri_aparc+aseg_REN_gm.nii.gz | tail -n +4 | sed 's/\"//g' | awk '{print $2}'`
# extract each roi into single-roi mask
foreach roi ($rois)
3dcalc -a "mri_aparc+aseg_REN_gm.nii.gz<$roi>" -expr 'step(a)' -prefix $roi.nii.gz
end
2086
Edited 1 time(s). Last edit at 06/19/2018 02:17AM by EXP.2086.