Hello all,

I have an (NHP) cluster map aligned to NMT, which I would like to cross check against the ARM atlas, to produce a report like:
cluster 1: 90% belongs to atlas area 75, 5% to atlas area 133, 5% to nowhere
cluster 2: 14% belongs to atlas area 22, 86 % to nowhere
and then
atlas area 55: 45% in cluster 1, 55% in cluster 3
atlas area 72: 12% in cluster 1, 88% in no cluster
not necessarily exactly in this form, of course. Atlas area names might be included as well.

Is there an AFNI program that would do that for me, ideally with a lot of control over output format?

Alternatively (preferably?) I thought I could just code it in Matlab. But is there a better way to load HEAD/BRIK into Matlab than converting to nifti with 3dAFNItoNIFTI, followed by niftiread in Matlab? I googled up this: [sscc.nimh.nih.gov] but it seems to be old (possibly outdated?) and the link does not work anyway.

Or, if you have other suggestions how to achieve my goal, I'd be happy to hear them.

Two caveats:
If own coding is needed, I would rather stay in Matlab: I know very little of Python and none of R, and coding anything slightly complex in tcsh is not my favorite way of spending time
My Matlab and my AFNI live on the same computer but in different OSes, so I would prefer to avoid mixing up these two, like calling AFNI programs from Matlab script.