脳全体量算出のスクリプト

Vmask = spm_vol(spm_select(1,'image','Select the binary mask'));
mask=spm_load_float(Vmask);
Vdata = spm_vol(spm_select(Inf,'image','Select modulated warped GM'));
volume = zeros(numel(Vdata),1);
for ii=1:length(Vdata)
dat = spm_load_float(Vdata(ii,:));
dat = dat.*mask;
dat(~isfinite(dat)) = 0;
voxvol = abs(det(Vdata(ii).mat))/100^3;
volume(ii) = sum(dat(:))*voxvol;
end;