Dear Eugenio,
"assign to each voxel the label with the highest posterior
probability" implies that each voxel must belong to one
and only one label, right?
I've found that this is not the case with the posterior_*.mgz files outputted by freesurfer,
where each voxel (especially the ones with the lower prob values)
can belong to multiple labels.
The passages I've done to come up at this conclusion was:
1) binarize all the posterior maps
2) summing them up
3) search for values greater than 1
4) if there are values >1, then there are voxels belonging to
multiple labels.
I've found voxels belonging to 7 different labels at the same time!
I think that the posterior files outputted by freesurfer are missing
the final step of required to each voxel to the single label
with the highest posterior prob.
I think that a clarification from the developers is needed here.
Best regards,
Luigi.