When we compute the aseg.stats, we include partial volume correction, we don't just count up the voxels in hippo.
doug
Guang Zeng wrote:
Hi, there,
I use the following matlab code to find right and left hippocampus in aseg.mgz
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% filename = 'case_1/mri/aseg.mgz'; [vol, M, mr_parms, volsz] = load_mgh(filename); hippo_left = find( (vol > 16.99) & (vol < 17.01) ); hippo_right = find( (vol > 52.99) & (vol < 53.01) ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
I found size(hippo_left, 1) = 3075; size(hippo_right, 1) = 3237; However, the value I found in aseg.mgz is
3162 3162.0 Left-Hippocampus 3390 3390.0 Right-Hippocampus
Why these values are different (3075 < 3162, 3237 < 3390)?
Thanks a lot!
Guang
Hotmail® has ever-growing storage! Don’t worry about storage limits. Check it out. http://windowslive.com/Tutorial/Hotmail/Storage?ocid=TXT_TAGLM_WL_HM_Tutorial_Storage_062009
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer