mfbox_read_spmfiles.m

toolbox for spm 5 for data, model free analysis

function data=mfbox_read_spmfiles(filename)
% data=mfbox_read_spmfiles(filename)
% read in image data from 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

data = {};

%load file infos
warning off all
V = spm_vol(filename);
warning query all

data.timesteps = 0;
data.dim = [0,0];

while (1)
    no_volumes = numel(V);
    if (no_volumes == 0)
        return
    end

    n = 1;
    while (all(V(1).dim == V(n).dim))
        if (n == no_volumes) 
            break; 
        end, 
        n=n+1; 
    end
    if (n == no_volumes)
        break;
    else
        questdlg(sprintf('Not all volumes share the same dimensions\nCannot import\n'),'Info','Ok','Ok');
        return
    end
end

imagedim = V(1).dim; % the new image dimensions

data.private = V;
data.timesteps = no_volumes;
data.dim = imagedim;

return