Hi Lijie (cc to mailing list),

I wrote a script to do just this in July. Here is a link to a tar file. It requires knowledge of matlab and bash, but is pretty straight forward. It

a) transforms the subfield probabilities to nu.mgz space (bash, using mri_label2vol)

b) applies a threshold to the maps (Koen suggested 50% of max probability, hence 128 of 255) to avoid over-estimation of the subfield volumes (bash)

c) comparing the probabilities of individual subfields in IJK space and assigning a unique label to each subfield

We can continue this discussion off-line if you like