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