mfbox_write_spmfiles.m

toolbox for spm 5 for data, model free analysis

function mfbox_write_spmfiles(filename,X,V)
% mfbox_write_spmfiles(filename)
% write image data to spm ANALYZE files
%
% Copyright by Peter Gruber, Fabian J. Theis and Ingo R. Keck
% Signal Processing & Information Theory group
% Institute of Biophysics, University of Regensburg, Germany
% Homepage: http://research.fabian.theis.name
%           http://www-aglang.uni-regensburg.de
%
% This file is free software, subject to the 
% GNU GENERAL PUBLIC LICENSE, see gpl.txt

s = size(X);
X = reshape(X,[],s(end));
if (prod(s(1:(end-1)))~=prod(V(1).dim))
    questdlg(sprintf('Not all volumes share the same dimensions\nCannot export\n'),'Info','Ok','Ok');
    return
end

sf = max(max(X(:)),-min(X(:)))/32767;
ni = nifti;
%ni.dat = file_array(filename,[V(1).dim,s(end)],'INT32-BE',0,sf,0);
ni.dat = file_array(filename,[V(1).dim,s(end)],'FLOAT32',0,sf,0);
ni.mat = V(1).private.mat;
ni.mat0 = V(1).private.mat;
ni.descrip = '4D image';
create(ni);

for i=1:s(end)
    ni.dat(:,:,:,i) = reshape(X(:,i),V(1).dim);
    spm_get_space([ni.dat.fname ',' num2str(i)], V(1).mat);
end