Save annatomical connectivity matrix w.r.t. connection strength and mean fiber length Input dir_name: Directory name where connection parameters are saved Nlabel: Number of label files cfile_base: basename of MAT files containing connection parameters cmatfile: MAT file containning connectivity matrices to be saved [absolute path] 2012/06/19 M.Fukushima 2012/06/29 M.Fukushima 2012/06/30 M.Fukushima Copyright (C) 2011, ATR All Rights Reserved. License : New BSD License(see VBMEG_LICENSE.txt)
0001 function save_anat_matrix(dir_name, Nlabel, cfile_base, cmatfile) 0002 % Save annatomical connectivity matrix w.r.t. connection strength and 0003 % mean fiber length 0004 % Input 0005 % dir_name: Directory name where connection parameters are saved 0006 % Nlabel: Number of label files 0007 % cfile_base: basename of MAT files containing connection parameters 0008 % cmatfile: MAT file containning connectivity matrices to be saved 0009 % [absolute path] 0010 % 0011 % 2012/06/19 M.Fukushima 0012 % 2012/06/29 M.Fukushima 0013 % 2012/06/30 M.Fukushima 0014 % 0015 % Copyright (C) 2011, ATR All Rights Reserved. 0016 % License : New BSD License(see VBMEG_LICENSE.txt) 0017 0018 c_count = zeros(Nlabel); 0019 f_leng = zeros(Nlabel); 0020 Ntrack = zeros(Nlabel,1); 0021 t_count = zeros(Nlabel,1); 0022 count = []; 0023 0024 for ii = 1:Nlabel 0025 cfile = fullfile(dir_name, filesep, [cfile_base num2str(ii) '.mat']); 0026 load(cfile) 0027 count = count(1:Nlabel); 0028 mfiber_length = mfiber_length(1:Nlabel); 0029 0030 c_count(:,ii) = count; 0031 f_leng(:,ii) = mfiber_length; 0032 Ntrack(ii) = length(Noverlap); % 2012/09/29 M.Fukushima 0033 if isfield(track_parm,'total_count')==0 0034 t_count(ii) = 0; 0035 else 0036 t_count(ii) = str2num(track_parm.total_count); 0037 end 0038 0039 if mod(ii,10)==0 0040 fprintf('%d/%d done.\n',ii,Nlabel); 0041 elseif ii == Nlabel 0042 fprintf('%d/%d done.\n',ii,Nlabel); 0043 end 0044 end 0045 0046 save(cmatfile, 'c_count', 'f_leng', 'Ntrack', 't_count')