Hi Ivan,
Here is a potential plan of attack
- get a fixel mask of the region (for example, CST). This can be computed using
tck2fixel $TractSeg_region_tck $fixel_directory_in $fixel_masks_dir $Region_fixel_mask
(assuming that TractSeg was used to generate the CST bundle). - compute the average FDC metric for each subject in the region. This can be done using
mrstats <FDC_subject#.mif> -mask $Region_fixel_mask
. To print out only the mean, consider adding in the option-output mean
. Loop through the subject list and print (and append) the output to a text file (along with the subject label).
See here and here for more information on this topic
Cheers,
Arkiev