rest_loadmask.m

While resting-state fMRI is drawing more and more attention, there has not been a software for its d

function mask=rest_loadmask(AX, AY, AZ, AMaskFilename)
%Load mask for REST by Xiao-Wei Song
% AX, AY, AZ should be the same size as the Volume, this size should be set automatically, but I use this for compatibility. And I specially give the code  for many files' use already. And another reason is the null mask need it
% AMaskFilename should be the mask filename
%------------------------------------------------------------------------------------------------------------------------------
%	Copyright(c) 2007~2010
%	State Key Laboratory of Cognitive Neuroscience and Learning in Beijing Normal University
%	Written by Xiao-Wei Song 
%	http://resting-fmri.sourceforge.net
% dawnsong, 20070509
%------------------------------------------------------------------------------------------------------------------------------
[e-mail] found Ynng He's bug in loading a null-mask
% 	<a href="[e-mail]">Mail to Author</a>: Xiaowei Song
%	Version=1.0;
%	Release=20070903;	
	
	
%Load mask, copy from reho.m revised by Dawnwei.Song, 20070504
[pathstr, name, ext, versn] = fileparts(mfilename('fullpath'));

if ( strcmp(AMaskFilename, '')|| (isnumeric(AMaskFilename) && AMaskFilename==0) )% like the old parameter, back-compatible  Xiaowei Song, 20070421
	mask=ones(AX, AY, AZ);
elseif( strcmpi(AMaskFilename, 'Default')||( isnumeric(AMaskFilename) && AMaskFilename==1) ) % like the old parameter	, back-compatible  Xiaowei Song, 20070421
	switch int2str([AX, AY, AZ])
	    case '79  95  69'   % 'default''[2 2 2]'
	        load([pathstr '/mask/reho_mask02.mat']);
	    case '53  63  46'   % 'default''[3 3 3]'	        
			load([pathstr '/mask/reho_mask03.mat']);
	    case '91  109  91'  % 'template''[2 2 2]'	        
			load([pathstr '/mask/reho_mask12.mat']);
	    case '61  73  61'   % 'template' '[3 3 3]'	        
			load([pathstr '/mask/reho_mask13.mat']);
	    otherwise
	        error(sprintf('There are no appropriate default mask file:\n\tVolume size=79*95*69 ,Voxel size=2*2*2;\n\tVolume size=53*63*46, Voxel size=3*3*3;\n\tVolume size=91*109*91, Voxel size=2*2*2;\n\tVolume size=61*73*61, Voxel size=3*3*3;\n Please set bMask = 0.'));
	end %end switch
else		% new,  Xiaowei Song, 20070421
	if (ischar(AMaskFilename))
		%if is img file, read and load
		if strcmpi(AMaskFilename(end-3:end), '.img')
			[mask, vsizeTmp, origin]=rest_readfile(AMaskFilename);            
		else%if is mat file, direct load			
			load(AMaskFilename);
		end
	else
		error('There are no appropriate mask file. Please set bMask = 0.');
	end
end %mask select end

%Check whether mask is 3d
if ndims(mask)~=3
	error('mask error, mask is not 3d');
end	

% Brutely check the mask's size to make sure whether the mask's size same to the required size
if ~all(size(mask)==[AX, AY, AZ]),	
	%warning(sprintf('\n\tMask does not match. Brutely use "No Mask".\n\tMask size is %dx%dx%d, not equal to required size %dx%dx%d',size(mask), [AX, AY, AZ]));
	
	%20070820 Zang's advice
	error(sprintf('\n\tMask does not match.\n\tMask size is %dx%dx%d, not same with required size %dx%dx%d',size(mask), [AX, AY, AZ]));
	
	mask = ones(AX, AY, AZ);	
end