Get area index corresponding to atlas map in the given area [ix_list, ix_name] = vb_get_atlas_label(ix_area,atlasfile,atlas_id) --- Input ix_area : Vertex index for given area atlasfile : File name of subject atlas file (*.atlas.mat) atlas_id : Atlas name ('aal', 'brodmann') --- Output ix_list : set of area index corresponding to atlas region ix_name : name of atlas region corresponding to 'ix_list' % 2006/11/14 M.Sato Copyright (C) 2011, ATR All Rights Reserved. License : New BSD License(see VBMEG_LICENSE.txt)
0001 function [ix_list, ix_name] = vb_get_atlas_label(ix_area,atlasfile,atlas_id) 0002 % Get area index corresponding to atlas map in the given area 0003 % 0004 % [ix_list, ix_name] = vb_get_atlas_label(ix_area,atlasfile,atlas_id) 0005 % --- Input 0006 % ix_area : Vertex index for given area 0007 % atlasfile : File name of subject atlas file (*.atlas.mat) 0008 % atlas_id : Atlas name ('aal', 'brodmann') 0009 % --- Output 0010 % ix_list : set of area index corresponding to atlas region 0011 % ix_name : name of atlas region corresponding to 'ix_list' 0012 % 0013 %% 2006/11/14 M.Sato 0014 % 0015 % Copyright (C) 2011, ATR All Rights Reserved. 0016 % License : New BSD License(see VBMEG_LICENSE.txt) 0017 0018 % Get atlas label 0019 Act = vb_get_act(atlasfile,atlas_id); 0020 Vlabel = Act.xxP; 0021 label = Act.label; 0022 label_name = Act.label_name; 0023 dd_label = Act.dd_label; 0024 0025 Npoint = length(Vlabel); 0026 flg = zeros(Npoint,1); 0027 flg(ix_area) = 1; 0028 0029 Vlabel = Vlabel.*flg; 0030 Nlabel = length(label); 0031 ix_list = []; 0032 ix_name = []; 0033 Narea = 0; 0034 0035 for n=1:Nlabel 0036 ix = find( Vlabel == label(n)); 0037 0038 if ~isempty(ix) 0039 Narea = Narea + 1; 0040 ix_list{Narea} = ix; 0041 ix_name{Narea} = label_name{n}; 0042 end 0043 end 0044 0045 %return