Hi Meenal,
You can do this from the Matlab command line using "spm_global". Here's the help for that function:
>> help spm_global
returns the global mean for a volume image - a compiled routine
FORMAT GX = spm_global(V)
V - memory mapped volume
GX - mean global activity
__________________________________________________________________________
spm_global returns the mean counts integrated over all the
slices from the volume.
The mean is estimated after discounting voxels outside the object
using a criteria of greater than > (global mean)/8.
__________________________________________________________________________
Copyright (C) 2008 Wellcome Trust Centre for Neuroimaging
The input "V" is an image file that has been loaded using spm_vol. Here a simple script that will let you select multiple images and will return their global means:
function glob=get_global()
vol=spm_vol(spm_select());
for g=1:length(vol)
glob(g,1)=spm_global(vol(g));
end
Hope this helps,
Paul
|