some kind of a script to do it. A possible example that you could copy and
paste into Matlab would be:
%____________________________________________________________
V = spm_vol(spm_get(1,'*.img'));
sm = 0;
for i=1:V.dim(3),
M = spm_matrix([0 0 i]);
img = spm_slice_vol(V,M,V.dim(1:2),0);
sm = sm + sum(img(:));
end;
sm
fprintf('%g litres\n', sm*det(V.mat(1:3,1:3))/100^3)
%____________________________________________________________
The value of sm is the integral of intensity throughout the image. If this
is an integral over a grey matter image, then you will need to scale by the
volume of each voxel ( det(V.mat(1:3,1:3)) )in order to get an answer in
units of cubic millimetres.
--
Sophie Josee Lafaille, M.Sc.
Research Coordinator
Law and Mental Health Program
Specialty Clinics, 4th Floor, RM 440
Centre for Addiction and Mental Health
250 College Street, Toronto, Ontario, Canada
M5T 1R8
Office: 416-535-8501 Ext. 6278