mfbox_getmask.m

toolbox for spm 5 for data, model free analysis

function [tg,mask,smask,xmask]=mfbox_getmask(grid,mask)
% Copyright by Peter Gruber 
% 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

sgr = size(grid);
sm = size(mask);
g = mat2cell(min(repmat(sm',1,sgr(2)),max(ones(sgr), ...
    round(grid))),ones(1,sgr(1)),sgr(2));
tg = sub2ind(sm,g{:});
[m,smask,xmask] = unique(tg);
xmask = smask;
mask = false(sm); mask(m) = true; % to optimize if no resort is done
if (issorted(tg) && length(smask(:))==length(tg(:)))
    smask = [];
end

return