Hi-

This is probably a question for Koen, but perhaps others could help me out.

I want to read into Matlab the ‘posterior_*.mgz’ files generated by the subfield segmentation. I’m not exactly sure the format/content of these files, but in addition to any information (e.g., some coordinate reference), I’m assuming it contains the posterior probability (0 – 255, scaled) for each (sub) voxel for each of the regions/fields.

I tried the obvious Matlab file (i.e., load_mgh.m), but that results in an error (the current Matlab directory contains the files and if I include the full path the same error obtains):

>> [vol M] = load_mgh('hippoposterior_right_CA4-DG.mgz');

gzip: hippoposterior_right_CA4-DG.mgz.gz: No such file or directory

ERROR: problem reading fname

SWITCH expression must be a scalar or string constant.

Error in load_mgh (line 156)

switch type

Given what I know (or think I know) about the procedure by which the subfields are labeled, it likely isn’t surprising the genearl Matlab scripts don’t properly deal with the posterior subfield files.

My objective is to be able to evaluate the distribution of voxels values (posterior probabilities) which exist for each of the subfields/regions.

Thanks,

Wil