Dear SPM experts, Even after reading the related posts on this list I am wondering if I am correctly computing an scalar grey matter volume in a particular ROI. First of all, is there any difference between GM volume and GM density? Here is my pseudocode approach to compute the GM volume for each processed scan (smwc1) in MNI space: gmfile = read ('smwc1XXX.nii'); roifile = read('roi.nii'); %binary roi in MNI space and with the same dimensions that gmfile prod = gmfile .* roifile; GM_volume = sum(prod(:)); Is that correct? Thanks a lot! Fred