%ERP_SHOWPLOT_UI display ERP waveforms on a figure % % USAGE: erp_showplot_ui(fig) % % ERP_SHOWPLOT_UI(fig) will display ERP waveforms on 'fig'. % % Called by ERP_PLOT_UI, ERP_PLOT_OPTION_UI % % I (fig) - figure where ERP waveform will be displayed on % % Created on 25-NOV-2002 by Jimmy Shen % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function erp_showplot_ui(fig) if ~exist('fig','var') fig = gcf; end system = getappdata(fig,'system'); time_info = getappdata(fig,'time_info'); timepoint = time_info.timepoint; start_timepoint = time_info.start_timepoint; start_time = time_info.start_time; end_time = time_info.end_time; digit_interval = time_info.digit_interval; init_option = getappdata(fig,'init_option'); if isempty(init_option) init_flag = 1; else init_flag = 0; end ha = getappdata(fig,'ha'); % main axes hm_chan_name = getappdata(fig,'hm_chan_name'); need_chan_name = get(hm_chan_name,'userdata'); hm_chan_axes = getappdata(fig,'hm_chan_axes'); need_chan_axes = get(hm_chan_axes,'userdata'); hm_chan_tick = getappdata(fig,'hm_chan_tick'); need_chan_tick = get(hm_chan_tick,'userdata'); font_size_list = [6 8 10 11 12 14 16 18]; if(init_flag) %% if first time run showplot % delete possible old_legend % old_legend_hdl = getappdata(fig,'legend_hdl'); if ~isempty(old_legend_hdl) try delete(old_legend_hdl{1}); hm_legend = getappdata(fig,'hm_legend'); setappdata(fig,'legend_hdl',[]); set(hm_legend,'Userdata',1,'Label','&Legend on'); catch end end old_txtbox_hdl = getappdata(fig,'txtbox_hdl'); if ~isempty(old_txtbox_hdl) try delete(old_txtbox_hdl); setappdata(fig,'txtbox_hdl',[]); catch end end view_option = getappdata(fig,'view_option'); color_code = getappdata(fig,'color_code'); % color & linestyle bs_color_code = getappdata(fig,'bs_color_code'); wave_selection = getappdata(fig,'wave_selection'); avg_selection = getappdata(fig,'avg_selection'); bs_selection = getappdata(fig,'bs_selection'); wave_selection(find(wave_selection==0)) = []; avg_selection(find(avg_selection==0)) = []; bs_selection(find(bs_selection==0)) = []; color_selection = 1:(length(wave_selection)+length(avg_selection)); bs_color_selection = 1:length(bs_selection); if length(color_selection) > length(color_code) % need more color tmp = []; for i=1:ceil( length(color_selection)/length(color_code) ) tmp = [tmp; color_code]; end color_code = tmp; end if length(bs_color_selection)>length(bs_color_code) % need more color tmp = []; for i=1:ceil( length(bs_color_selection)/length(bs_color_code) ) tmp = [tmp; bs_color_code]; end bs_color_code = tmp; end x_interval_selection = getappdata(fig,'x_interval_selection'); y_interval_selection = getappdata(fig,'y_interval_selection'); % display wave size scale, from 0.05 to 1, 1 means full screen % eta = getappdata(fig,'eta'); % font size % font_size = font_size_list(getappdata(fig,'font_size_selection')); color_code1 = color_code(1:length(avg_selection),:); color_code2 = color_code((length(avg_selection)+1): ... length(color_selection),:); color_code3 = bs_color_code(1:length(bs_selection),:); wave_amplitude = getappdata(fig,'wave_amplitude'); % input wave wave_name = getappdata(fig,'wave_name'); % wave legend selected_wave_info = getappdata(fig,'selected_wave_info'); % [condition(selected), subjects(selected)] selected_wave = getappdata(fig,'selected_wave'); % based on % selected_cond & selected_subj selected_wave_idx = selected_wave(wave_selection); selected_channels = getappdata(fig,'selected_channels'); selected_chan_idx = find(selected_channels); wave_disp = wave_amplitude(:,selected_chan_idx,selected_wave_idx); selected_wave_name = wave_name(selected_wave_idx); selected_wave_name = char(selected_wave_name); % legend string num_wave = size(selected_wave_name,1); % # of wave to display avg_amplitude = getappdata(fig,'avg_amplitude'); % input avg avg_name = getappdata(fig,'avg_name'); % avg legend if ~isempty(avg_amplitude) selected_avg = find(getappdata(fig,'selected_conditions')); selected_avg_idx = selected_avg(avg_selection); avg_disp = avg_amplitude(:,selected_chan_idx,selected_avg_idx); selected_avg_name = avg_name(selected_avg_idx); else selected_avg = []; selected_avg_idx = []; avg_disp = []; selected_avg_name = []; end selected_avg_name = char(selected_avg_name); num_avg = size(selected_avg_name,1); bs_amplitude = getappdata(fig,'bs_amplitude'); % input bootstrap bs_name = getappdata(fig,'bs_name'); % bs legend if ~isempty(bs_amplitude) bs_x_amplitude = bs_amplitude; selected_bs = getappdata(fig,'selected_bs'); % 1 to lv num_lv = length(selected_bs); % re-order the LV sequence, so that LV1 will be put on top, % LV2 will be put on bottom, LV3 will be put next top, ... % bs_y_amplitude = []; for i = 1:ceil(num_lv / 2) bs_y_amplitude = [bs_y_amplitude, num_lv-i+1]; if (num_lv - i + 1) ~= i bs_y_amplitude = [bs_y_amplitude, i]; end end % normalize the new order, and reshape it to the required % y direction vector for bootstrap % if num_lv > 1 bs_y_amplitude = (bs_y_amplitude - 1) / (num_lv - 1); end; bs_y_amplitude = reshape(bs_y_amplitude, [1,1,num_lv]); bs_y_amplitude = repmat(bs_y_amplitude, ... [timepoint, size(bs_amplitude, 2)]); selected_bs_idx = selected_bs(bs_selection); selected_bs_name = bs_name(selected_bs_idx); else bs_y_amplitude = []; selected_bs_idx = []; bs_x_amplitude = []; selected_bs_name = []; end selected_bs_name = char(selected_bs_name); num_bs = size(selected_bs_name,1); if ~isempty(wave_disp) wave_max = max(wave_disp(:)); wave_min = min(wave_disp(:)); else wave_max = 0; wave_min = 0; end if ~isempty(avg_disp) avg_max = max(avg_disp(:)); avg_min = min(avg_disp(:)); else avg_max = 0; avg_min = 0; end wave_top = max(wave_max, avg_max); % top end of y axis if wave_max < 0 wave_top = 0; end wave_bottom = min(wave_min, avg_min); % bottom end of y axis if wave_min > 0 wave_bottom = 0; end if wave_top + wave_bottom == 0 msg = 'ERROR: No wave to display.'; uiwait(msgbox(msg,'ERROR','modal')); return; end if view_option == 3 north = [sprintf('%0.2f',wave_top)]; south = [sprintf('%0.2f',wave_bottom)]; else north = [sprintf('%0.2f',wave_top) ' \muV']; south = [sprintf('%0.2f',wave_bottom) ' \muV']; end west = sprintf('%d ms', start_timepoint * digit_interval); % east = sprintf('%d ms', (start_timepoint+timepoint-1) * digit_interval); east = sprintf('%d ms', (start_timepoint+timepoint) * digit_interval); % belows are 4 lines represent the channel axis % hor_x = [0 1]; hor_y = [0 0]; if start_timepoint < 0 ver_x = ((-start_timepoint)/(timepoint-1)) * [1 1]; else ver_x = [0 0]; end ver_y = [wave_bottom wave_top] / (wave_top - wave_bottom); xtick_length = 1/40; if start_timepoint < 0 xtick_x = ((-start_timepoint)/(timepoint-1)) * [1 1]; else xtick_x = [0 0]; end xtick_y = [-xtick_length xtick_length]; ytick_length = 1/80; if start_timepoint < 0 ytick_x = ((-start_timepoint)/(timepoint-1))+ ... [-ytick_length ytick_length]; else ytick_x = [-ytick_length ytick_length]; end ytick_y = [0 0]; % ------------------- below: set waveform origin ------------- switch system.class case 1 type_str = 'BESAThetaPhi|EGI'; switch system.type case 1 load('erp_chan148'); case 2 load('erp_egi_chan129'); end case 2 type_str = 'CTF-150'; switch system.type case 1 load('meg_ctf_chan150'); end end % convert all channel location from [theta,phi] to [x,y] % for i=1:size(chan_loc,1) [x(i) y(i)] = erp_getxy(chan_loc(i,2),chan_loc(i,1)); end min_x = min(x); max_x = max(x); min_y = min(y); max_y = max(y); % apply channel mask % chan_mask = getappdata(fig,'chan_mask'); chan_nam = chan_nam(chan_mask,:); chan_loc = chan_loc(chan_mask,:); % convert selected locations from [theta,phi] to [x,y] % for i=selected_chan_idx [x(i) y(i)] = erp_getxy(chan_loc(i,2),chan_loc(i,1)); end % normalize & shift x & y array % x = (x+abs(min_x)) / (abs(max_x)+abs(min_x)); y = (y+abs(min_y)) / (abs(max_y)+abs(min_y)); % ------------------- above: set waveform origin ------------- % ------------------- below: start to draw axes & wave ------------- if ~isempty(wave_disp) | ~isempty(avg_disp) % there is wave axes(ha); cla; hold on; % create template for sub axis % axis_template = plot(hor_x,hor_y,'k-',ver_x,ver_y,'k-', ... 'visible','off'); xtick_template_hdl = plot(xtick_x,xtick_y,'k-', ... 'visible','off'); ytick_template_hdl = plot(ytick_x,ytick_y,'k-', ... 'visible','off'); % xtick % x_interval_pattern = [5 2 1]; % nice $ pattern x_interval = x_interval_pattern; while(x_interval(1) <= timepoint) x_interval = [x_interval_pattern.*10, x_interval]; x_interval_pattern = x_interval_pattern.*10; end x_interval(find(x_interval > timepoint))=[]; % take out too large x_interval(1) = []; % make large one close to half timepoint x_interval = x_interval .* digit_interval; % become (ms) if isempty(x_interval_selection) % default x_interval2 = x_interval(1); else x_interval2 = x_interval(x_interval_selection); end xtick_start = ceil(start_time/x_interval2) * x_interval2; xtick_end = floor(end_time/x_interval2) * x_interval2; xtick_pos = [xtick_start:x_interval2:xtick_end]; % (ms) if xtick_start > 0 xtick_pos = xtick_pos - start_time; % shift y axis end xtick_pos = xtick_pos ./ digit_interval; % point again xtick_template = []; % traverse x axis and put tick on % k = 1; for j = xtick_pos xtick_template(k,:) = copyobj(xtick_template_hdl, ha); set(xtick_template(k,:),'xdata', ... get(xtick_template(k,:),'xdata') ... + j/timepoint); if j==0 set(xtick_template(k,:),'user','origin'); end k = k + 1; end if isempty(j) j = 0; end % ytick % if abs(wave_top) > abs(wave_bottom) wave_long = abs(wave_top); else wave_long = abs(wave_bottom); end % choose [2 4 8] as divisor to determin ytick interval % y_interval(1) = rri_ceil(wave_long, 2); y_interval(2) = rri_ceil(wave_long, 4); y_interval(3) = rri_ceil(wave_long, 8); if isempty(y_interval_selection) % interval not selected y_interval2 = y_interval(1); % select 1st one else y_interval2 = y_interval(y_interval_selection); end ytick_start = ceil(wave_bottom/y_interval2) * y_interval2; ytick_end = floor(wave_top/y_interval2) * y_interval2; ytick_pos = [ytick_start:y_interval2:ytick_end]; ytick_template = []; % traverse upper half of y axis and put tick on % k = 1; for j = ytick_pos ytick_template(k,:) = copyobj(ytick_template_hdl, ha); set(ytick_template(k,:),'ydata', ... get(ytick_template(k,:),'ydata') ... +j/(wave_top - wave_bottom)); if j==0 set(ytick_template(k,:),'user','origin'); end k = k + 1; end % create legend axis % offset = 0.08; % legend axis offset to left offset2 = 0; % legend axis offset to bottom offset3 = 0.01; % axis label offset to axis laxis_hdl = copyobj(axis_template, ha); for j=1:2 % just for 2 axis line set(laxis_hdl(j), 'visible', 'on', ... 'xdata',eta*get(laxis_hdl(j),'xdata') ... +offset, ... 'ydata',eta*(get(laxis_hdl(j),'ydata') ... +abs(wave_bottom)/(abs(wave_bottom)+abs(wave_top))) ... +offset2); end % create legend tickmark % if ~isempty(xtick_template) lxtick_hdl = copyobj(xtick_template, ha); for j=1:size(lxtick_hdl,1) % for all the components on xtick set(lxtick_hdl(j), 'visible', 'on', ... 'xdata',eta*get(lxtick_hdl(j),'xdata') ... +offset, ... 'ydata',eta*(get(lxtick_hdl(j),'ydata') ... +abs(wave_bottom)/(abs(wave_bottom)+abs(wave_top))) ... +offset2); if strcmp(get(lxtick_hdl(j),'user'),'origin') set(lxtick_hdl(j), 'visible', 'off'); end % create legend tickmark label % % lxticklabel_x = get(lxtick_hdl(j),'xdata'); % lxticklabel_x = lxticklabel_x(1); % lxticklabel_y = get(lxtick_hdl(j),'ydata'); % lxticklabel_y = lxticklabel_y(1); % lxticklabel_hdl(j) = text(lxticklabel_x, lxticklabel_y, ... % num2str(x_interval*j), ... % 'horizon', 'center', 'vertical', 'top', ... % 'fontsize',12,'clipping','on'); end else lxtick_hdl = []; end if ~isempty(ytick_template) lytick_hdl = copyobj(ytick_template, ha); for j=1:size(lytick_hdl,1) % for all the components on ytick set(lytick_hdl(j), 'visible', 'on', ... 'xdata',eta*get(lytick_hdl(j),'xdata') ... +offset, ... 'ydata',eta*(get(lytick_hdl(j),'ydata') ... +abs(wave_bottom)/(abs(wave_bottom)+abs(wave_top))) ... +offset2); if strcmp(get(lytick_hdl(j),'user'),'origin') set(lytick_hdl(j), 'visible', 'off'); end % create legend tickmark label % % lyticklabel_x = get(lytick_hdl(j),'xdata'); % lyticklabel_x = lyticklabel_x(2); % lyticklabel_y = get(lytick_hdl(j),'ydata'); % lyticklabel_y = lyticklabel_y(2); % lyticklabel_hdl(j) = text(lyticklabel_x, lyticklabel_y, ... % num2str(y_interval2*j), ... % 'horizon', 'left', 'vertical', 'middle', ... % 'fontsize',12,'clipping','on'); end else lytick_hdl = []; end % create legend axis label % %-------------------------------- % north % if start_timepoint < 0 ltext_x = eta * ((-start_timepoint)/(timepoint-1)) + offset; else ltext_x = offset; end ltext_y = eta + offset2 + offset3; legend_text_hdl(1) = text(ltext_x,ltext_y, ... deblank(north), ... 'horizon', 'center', 'vertical', 'bottom', ... 'fontsize',font_size,'clipping','on'); % south % if start_timepoint < 0 ltext_x = eta * ((-start_timepoint)/(timepoint-1)) + offset; else ltext_x = offset; end ltext_y = offset2 - offset3; legend_text_hdl(2) = text(ltext_x,ltext_y, ... deblank(south), ... 'horizon', 'center', 'vertical', 'top', ... 'fontsize',font_size,'clipping','on'); % east % ltext_x = eta + offset + offset3; ltext_y = eta * abs(wave_bottom)/(abs(wave_bottom) ... + abs(wave_top)) + offset2; legend_text_hdl(3) = text(ltext_x,ltext_y, ... deblank(east), ... 'horizon', 'left', 'vertical', 'middle', ... 'fontsize',font_size,'clipping','on'); % west % ltext_x = offset - offset3; ltext_y = eta * abs(wave_bottom)/(abs(wave_bottom) ... + abs(wave_top)) + offset2; legend_text_hdl(4) = text(ltext_x,ltext_y, ... deblank(west), ... 'horizon', 'right', 'vertical', 'middle', ... 'fontsize',font_size,'clipping','on'); axis_hdl = []; xtick_hdl = []; ytick_hdl = []; selected_subjects = getappdata(fig,'selected_subjects'); selected_conditions = getappdata(fig,'selected_conditions'); % create new waves, and save the handles % for i=selected_chan_idx % draw sub axes % axis_hdl(:,i) = copyobj(axis_template, ha); for j=1:size(axis_hdl,1) % for all the components on axis set(axis_hdl(j,i), ... 'xdata',eta*get(axis_hdl(j,i),'xdata')+x(i), ... 'ydata',eta*get(axis_hdl(j,i),'ydata')+y(i)); end if ~isempty(xtick_template) xtick_hdl(:,i) = copyobj(xtick_template, ha); for j=1:size(xtick_hdl,1) % for all the components on xtick set(xtick_hdl(j,i), ... 'xdata',eta*get(xtick_hdl(j,i),'xdata')+x(i), ... 'ydata',eta*get(xtick_hdl(j,i),'ydata')+y(i)); if strcmp(get(xtick_hdl(j,i),'user'),'origin') set(xtick_hdl(j,i), 'visible', 'off'); end end end if ~isempty(ytick_template) ytick_hdl(:,i) = copyobj(ytick_template, ha); for j=1:size(ytick_hdl,1) % for all the components on ytick set(ytick_hdl(j,i), ... 'xdata',eta*get(ytick_hdl(j,i),'xdata')+x(i), ... 'ydata',eta*get(ytick_hdl(j,i),'ydata')+y(i)); if strcmp(get(ytick_hdl(j,i),'user'),'origin') set(ytick_hdl(j,i), 'visible', 'off'); end end end x_wave = eta * (0:timepoint-1)/(timepoint-1) + x(i); if num_avg ~= 0 seq_idx = 1; for j = selected_avg_idx y_wave = eta * ... avg_amplitude(:,i,j)/ ... (abs(wave_bottom)+abs(wave_top))+y(i); avg_wave_hdl(i,j) = plot(x_wave, y_wave, ... color_code1(seq_idx,:), ... 'linewidth', 2, ... 'userdata',[seq_idx,i]); seq_idx = seq_idx + 1; end % for selected_avg_idx else avg_wave_hdl = []; end if num_wave ~= 0 seq_idx = 1; for j = selected_wave_idx y_wave = eta * ... wave_amplitude(:,i,j)/ ... (abs(wave_bottom)+abs(wave_top))+y(i); if view_option == 1 condition_num = selected_wave_info(j,1); subject_num = selected_wave_info(j,2); else condition_num = []; subject_num = []; end % wave_hdl(channel, wave) % userdata [wave_sequence, channel, condition, subject] % wave_hdl(i,j) = plot(x_wave, y_wave, ... color_code2(seq_idx,:), ... 'buttondown','erp_plot_ui(''select_wave'');',... 'userdata', [seq_idx, i, condition_num, subject_num]); seq_idx = seq_idx + 1; end % for num_wave else wave_hdl = []; end if num_bs ~= 0 seq_idx = 1; for j = selected_bs_idx x_mask = bs_x_amplitude(:,i,j); x_mask = find(x_mask)'; if ~isempty(x_mask) bs_x_wave = x_wave(x_mask); bs_y_wave = bs_y_amplitude(:,i,j); bs_y_wave = eta * (bs_y_wave(x_mask) - abs(wave_bottom) ... / (abs(wave_bottom) + abs(wave_top))) ... + y(i); bs_wave_hdl(i,j) = plot(bs_x_wave, bs_y_wave, ... color_code3(seq_idx,:), ... 'markersize',6, ... 'userdata',[seq_idx,i]); else bs_wave_hdl(i,j) = 99999; end seq_idx = seq_idx + 1; end % for selected_bs_idx else bs_wave_hdl = []; end % create new chan_names, and save the handles % if start_timepoint < 0 text_x = x + eta*((-start_timepoint)/(timepoint-1)); else text_x = x; end text_y = y + eta*wave_top/(abs(wave_bottom)+abs(wave_top)); chan_name_hdl(i) = text(text_x(i),text_y(i), ... deblank(chan_nam(i,:)), 'userdata', [0, i], ... 'buttondown','erp_plot_ui(''select_chan_name'');', ... 'horizon', 'center', 'vertical', 'bottom', ... 'fontsize',font_size,'clipping','on'); end % for selected_chan_idx hold off; % create plot title % switch view_option case {1} h_title=text(.03,.99,'ERP Amplitude', ... 'fontsize',font_size, 'fontweight','bold'); case {2} h_title=text(.03,.99,'Grand Average Amplitude', ... 'fontsize',font_size, 'fontweight','bold'); case {3} h_title=text(.03,.99,'Electrode Salience', ... 'fontsize',font_size, 'fontweight','bold'); case {4} h_title=text(.03,.99,'Group Subject Amplitude', ... 'fontsize',font_size, 'fontweight','bold'); case {5} h_title=text(.03,.99,'Spatiotemporal Correlations', ... 'fontsize',font_size, 'fontweight','bold'); end setappdata(fig,'axis_template',axis_template); setappdata(fig,'xtick_template',xtick_template); setappdata(fig,'ytick_template',ytick_template); setappdata(fig,'laxis_hdl',laxis_hdl); setappdata(fig,'lxtick_hdl',lxtick_hdl); setappdata(fig,'lytick_hdl',lytick_hdl); setappdata(fig,'ltext_x',ltext_x); setappdata(fig,'ltext_y',ltext_y); setappdata(fig,'legend_text_hdl',legend_text_hdl); setappdata(fig,'offset',offset); setappdata(fig,'offset2',offset2); setappdata(fig,'offset3',offset3); setappdata(fig,'axis_hdl',axis_hdl); setappdata(fig,'xtick_hdl',xtick_hdl); setappdata(fig,'ytick_hdl',ytick_hdl); setappdata(fig,'wave_hdl',wave_hdl); setappdata(fig,'avg_wave_hdl',avg_wave_hdl); setappdata(fig,'bs_wave_hdl',bs_wave_hdl); setappdata(fig,'chan_name_hdl',chan_name_hdl); setappdata(fig,'h_title',h_title); setappdata(fig,'x_interval',x_interval); setappdata(fig,'y_interval',y_interval); setappdata(fig,'avg_amplitude',avg_amplitude); setappdata(fig,'wave_amplitude',wave_amplitude); setappdata(fig,'bs_x_amplitude',bs_x_amplitude); setappdata(fig,'bs_y_amplitude',bs_y_amplitude); setappdata(fig,'selected_wave_name',selected_wave_name); setappdata(fig,'selected_avg_name',selected_avg_name); setappdata(fig,'selected_bs_name',selected_bs_name); setappdata(fig,'selected_wave_idx',selected_wave_idx); setappdata(fig,'selected_avg_idx',selected_avg_idx); setappdata(fig,'selected_bs_idx',selected_bs_idx); setappdata(fig,'selected_chan_idx',selected_chan_idx); else % if there is NO wave setappdata(fig,'wave_hdl',[]); setappdata(fig,'avg_wave_hdl',[]); setappdata(fig,'bs_wave_hdl',[]); setappdata(fig,'legend_hdl',[]); setappdata(fig,'chan_name_hdl',[]); return; end % end of ~isempty(wave_disp) | ~isempty(avg_disp) init_option.color_code1 = color_code1; init_option.color_code2 = color_code2; init_option.color_code3 = color_code3; init_option.avg_max = avg_max; init_option.avg_min = avg_min; init_option.wave_max = wave_max; init_option.wave_min = wave_min; init_option.wave_top = wave_top; init_option.wave_bottom = wave_bottom; init_option.timepoint = timepoint; init_option.start_timepoint = start_timepoint; init_option.digit_interval = digit_interval; init_option.x = x; init_option.y = y; setappdata(fig,'init_option',init_option); else %%%%%%%%%%%%%%%%%%%%% if NOT first time run showplot %++++++++++++++++++++++++++++++++++++++++++++++++ init_option = getappdata(fig,'init_option'); avg_max = init_option.avg_max; avg_min = init_option.avg_min; wave_max = init_option.wave_max; wave_min = init_option.wave_min; wave_top = init_option.wave_top; wave_bottom = init_option.wave_bottom; timepoint = init_option.timepoint; start_timepoint = init_option.start_timepoint; x = init_option.x; y = init_option.y; orig_axis_hdl = getappdata(fig,'axis_template'); orig_xtick_hdl = getappdata(fig,'xtick_template'); orig_ytick_hdl = getappdata(fig,'ytick_template'); laxis_hdl = getappdata(fig,'laxis_hdl'); lxtick_hdl = getappdata(fig,'lxtick_hdl'); lytick_hdl = getappdata(fig,'lytick_hdl'); ltext_x = getappdata(fig,'ltext_x'); ltext_y = getappdata(fig,'ltext_y'); legend_text_hdl = getappdata(fig,'legend_text_hdl'); offset = getappdata(fig,'offset'); offset2 = getappdata(fig,'offset2'); offset3 = getappdata(fig,'offset3'); axis_hdl = getappdata(fig,'axis_hdl'); xtick_hdl = getappdata(fig,'xtick_hdl'); ytick_hdl = getappdata(fig,'ytick_hdl'); wave_hdl = getappdata(fig,'wave_hdl'); avg_wave_hdl = getappdata(fig,'avg_wave_hdl'); bs_wave_hdl = getappdata(fig,'bs_wave_hdl'); chan_name_hdl = getappdata(fig,'chan_name_hdl'); h_title = getappdata(fig,'h_title'); wave_amplitude = getappdata(fig,'wave_amplitude'); avg_amplitude = getappdata(fig,'avg_amplitude'); bs_x_amplitude = getappdata(fig,'bs_x_amplitude'); bs_y_amplitude = getappdata(fig,'bs_y_amplitude'); selected_wave_name = getappdata(fig,'selected_wave_name'); selected_avg_name = getappdata(fig,'selected_avg_name'); selected_bs_name = getappdata(fig,'selected_bs_name'); selected_wave_idx = getappdata(fig,'selected_wave_idx'); selected_avg_idx = getappdata(fig,'selected_avg_idx'); selected_bs_idx = getappdata(fig,'selected_bs_idx'); selected_chan_idx = getappdata(fig,'selected_chan_idx'); eta = getappdata(fig,'eta'); font_size = font_size_list(getappdata(fig,'font_size_selection')); num_wave = size(selected_wave_name,1); num_avg = size(selected_avg_name,1); num_bs = size(selected_bs_name,1); % update legend axis % for j=1:size(laxis_hdl,1) % for all the components on legend axis set(laxis_hdl(j), ... 'xdata',eta*get(orig_axis_hdl(j),'xdata') ... +offset, ... 'ydata',eta*(get(orig_axis_hdl(j),'ydata') ... +abs(wave_bottom)/(abs(wave_bottom)+abs(wave_top))) ... +offset2); end % update legend tick marks % if ~isempty(orig_xtick_hdl) for j=1:size(lxtick_hdl,1) % for all the components on legend xtick set(lxtick_hdl(j), ... 'xdata',eta*get(orig_xtick_hdl(j),'xdata') ... +offset, ... 'ydata',eta*(get(orig_xtick_hdl(j),'ydata') ... +abs(wave_bottom)/(abs(wave_bottom)+abs(wave_top))) ... +offset2); end end if ~isempty(orig_ytick_hdl) for j=1:size(lytick_hdl,1) % for all the components on legend ytick set(lytick_hdl(j), ... 'xdata',eta*get(orig_ytick_hdl(j),'xdata') ... +offset, ... 'ydata',eta*(get(orig_ytick_hdl(j),'ydata') ... +abs(wave_bottom)/(abs(wave_bottom)+abs(wave_top))) ... +offset2); end end % update legend axis label % %-------------------------------- % north % if start_timepoint < 0 ltext_x = eta * ((-start_timepoint)/(timepoint-1)) + offset; else ltext_x = offset; end ltext_y = eta + offset2 + offset3; set(legend_text_hdl(1),'position',[ltext_x,ltext_y,0], ... 'fontsize',font_size); % south % if start_timepoint < 0 ltext_x = eta * ((-start_timepoint)/(timepoint-1)) + offset; else ltext_x = offset; end ltext_y = offset2 - offset3; set(legend_text_hdl(2),'position',[ltext_x,ltext_y,0], ... 'fontsize',font_size); % east % ltext_x = eta + offset + offset3; ltext_y = eta * abs(wave_bottom)/(abs(wave_bottom) ... + abs(wave_top)) + offset2; set(legend_text_hdl(3),'position',[ltext_x,ltext_y,0], ... 'fontsize',font_size); % west % ltext_x = offset - offset3; ltext_y = eta * abs(wave_bottom)/(abs(wave_bottom) ... + abs(wave_top)) + offset2; set(legend_text_hdl(4),'position',[ltext_x,ltext_y,0], ... 'fontsize',font_size); for i=selected_chan_idx for j=1:size(axis_hdl,1) % for all the components on axis set(axis_hdl(j,i), ... 'xdata',eta*get(orig_axis_hdl(j),'xdata')+x(i), ... 'ydata',eta*get(orig_axis_hdl(j),'ydata')+y(i)); end if ~isempty(orig_xtick_hdl) for j=1:size(xtick_hdl,1) % for all the components on xtick set(xtick_hdl(j,i), ... 'xdata',eta*get(orig_xtick_hdl(j),'xdata')+x(i), ... 'ydata',eta*get(orig_xtick_hdl(j),'ydata')+y(i)); end end if ~isempty(orig_ytick_hdl) for j=1:size(ytick_hdl,1) % for all the components on ytick set(ytick_hdl(j,i), ... 'xdata',eta*get(orig_ytick_hdl(j),'xdata')+x(i), ... 'ydata',eta*get(orig_ytick_hdl(j),'ydata')+y(i)); end end x_wave = eta * (0:timepoint-1)/(timepoint-1) + x(i); if num_avg ~= 0 for j=selected_avg_idx y_wave = eta * ... avg_amplitude(:,i,j)/(abs(wave_bottom)+abs(wave_top))+y(i); set(avg_wave_hdl(i,j),'xdata',x_wave,'ydata',y_wave); end % for selected_avg_idx end if num_wave ~= 0 for j=selected_wave_idx y_wave = eta * ... wave_amplitude(:,i,j)/(abs(wave_bottom)+abs(wave_top))+y(i); set(wave_hdl(i,j),'xdata',x_wave,'ydata',y_wave); end end if num_bs ~= 0 for j = selected_bs_idx x_mask = bs_x_amplitude(:,i,j); x_mask = find(x_mask)'; bs_x_wave = x_wave(x_mask); bs_y_wave = bs_y_amplitude(:,i,j); bs_y_wave = eta * (bs_y_wave(x_mask) - abs(wave_bottom) ... / (abs(wave_bottom) + abs(wave_top))) ... + y(i); if bs_wave_hdl(i,j) ~= 99999 set(bs_wave_hdl(i,j),'xdata',bs_x_wave,'ydata',bs_y_wave); end end end if start_timepoint < 0 text_x = x + eta*((-start_timepoint)/(timepoint-1)); else text_x = x; end text_y = y + eta*wave_top/(abs(wave_bottom)+abs(wave_top)); set(chan_name_hdl(i),'position',[text_x(i),text_y(i)], ... 'fontsize',font_size); end % for selected_chan_idx set(h_title,'fontsize',font_size); % update title font size end %% if init_flag end if need_chan_name erp_plot_ui('display_chan_name','on'); else erp_plot_ui('display_chan_name','off'); end if need_chan_axes erp_plot_ui('display_chan_axes','on'); else erp_plot_ui('display_chan_axes','off'); end if need_chan_tick & need_chan_axes erp_plot_ui('display_chan_tick','on'); else erp_plot_ui('display_chan_tick','off'); end % below is the test routine, don't remove it % if(0) % set 1 for test, 0 for normal switch system.class case 1 type_str = 'BESAThetaPhi|EGI'; switch system.type case 1 load('erp_chan148'); case 2 load('erp_egi_chan129'); end case 2 type_str = 'CTF-150'; switch system.type case 1 load('meg_ctf_chan150'); end end % convert all channel location from [theta,phi] to [x,y] % for i=1:size(chan_loc,1) [x(i) y(i)] = erp_getxy(chan_loc(i,2),chan_loc(i,1)); end min_x = min(x); max_x = max(x); min_y = min(y); max_y = max(y); % normalize & shift x & y array % x = (x-min_x) / (max_x-min_x); y = (y-min_y) / (max_y-min_y); for i=1:size(chan_loc,1) text(x(i), y(i), chan_nam(i,:)); % text(x(i), y(i), num2str(i)); end end return; % erp_showplot_ui