 |
 |
Dear Radu
You could try the Volumes toolbox: http://sourceforge.net/projects/spmtools
Alternatively, the following code should do the trick, once you converted your data into nifti-files (first open these using spm_display and identify the slices which are of interest to you).
Vi = spm_vol(spm_select(inf,
'image','Select volumes'));for
i=1:numel(Vi);cut = spm_input(
'Specify lowest/highest slice','!+1','e',[1 Vi(i).dim(3)],2);dat = spm_read_vols(Vi(i));
dat = dat(:,:,cut(1):cut(2));
Vo = Vi(i);
Vo. fname = [Vi(i).fname(1:end-4)
'_' int2str(cut(1)) '-' int2str(cut(2)) '.nii'];Vo.dim(3) = cut(2)-cut(1)+1;
Vo.mat(3,4) = -Vo.dim(3)/2;
Vo = rmfield(Vo,
'pinfo');Vo = spm_write_vol(Vo,dat);
end
Best
Simon