vis_trajgui.m

it is matlab code , som(slef organizing map) tool for matlab

function vis_trajgui(trajStruct,arg)

% VIS_TRAJGUI subfuntion for SOM_TRAJECTORY
%
% This function is the actual GUI called by SOM_TRAJECTORY
% function. 
%
% See also SOM_TRAJECTORY.

% Contributed code to SOM Toolbox 2.0, February 11th, 2000 by Juha Parhankangas
% Copyright (c) by Juha Parhankangas.
% http://www.cis.hut.fi/projects/somtoolbox/        

% Version 2.0beta juha 180699

if nargin == 1

  sM_h=trajStruct.figure;

  if size(trajStruct.bmus,1) ~= 1 & size(trajStruct.bmus,2) ~= 1
    fuzzy_traj(trajStruct,[]);
    return;
  end
  
  
  if size(trajStruct.bmus,1) == 1 | size(trajStruct.bmus,2) == 1

    udata.bmus = trajStruct.bmus;
    udata.a_h=[findobj(get(sM_h,'Children'),'Tag','Uplane');...
	       findobj(get(sM_h,'Children'),'Tag','Cplane')];
    udata.sM_h=trajStruct.figure;
    udata.traj=[];
    data1 = trajStruct.primary_data;
    if ~isempty(trajStruct.primary_names)
      names=trajStruct.primary_names;
    else
      for i=1:size(data1,2)
	names{i,1}=sprintf('Var%d',i);
      end
    end

    udata.lattice=trajStruct.lattice;
    form = 0.7*vis_patch(udata.lattice);
    udata.msize = trajStruct.msize;


    %%%%%%%%%%%%%%%%%%%%%%%%
    %
    % forming a patch object, which is placed above every component plane
    %
    
    
    l = size(form,1);
    
    nx = repmat(form(:,1),1,prod(udata.msize));
    ny = repmat(form(:,2),1,prod(udata.msize));
    
    x=reshape(repmat(1:udata.msize(2),l*udata.msize(1),1),l,prod(udata.msize));
    y=repmat(repmat(1:udata.msize(1),l,1),1,udata.msize(2));
    
    if strcmp(udata.lattice,'hexa')
      t = find(~rem(y(1,:),2));
      x(:,t)=x(:,t)+.5;
    end
    x=x+nx;
    y=y+ny;
    
    colors=reshape(ones(prod(udata.msize),1)*[NaN NaN NaN],...
		   [1 prod(udata.msize) 3]);
    
    set(0,'CurrentFigure',udata.sM_h);
    
    %%%%%%%%%%%%%%%%%%%%%%
    %
    % drawing patch
    %
    % caxis -commands keep the colormap of the original patch unchanged.
    %
    
    for i=1:length(udata.a_h)
      udata.real_patch(i)=get(udata.a_h(i),'Children');
      set(udata.real_patch(i),'ButtonDownFcn',...
			'vis_trajgui([],''click'')');
      subplot(udata.a_h(i));
      v=caxis;
      udata.tmp_patch(i)=patch(x,y,colors,'EdgeColor','none',...
			       'ButtonDownFcn',...
			       'vis_trajgui([],''click'')',...
			       'Tag','TmpPatch');
      caxis(v);
    end
    
    %%%%%%%%%%%%%%%%%%%%
    
    
    
    udata.length_of_traj=length(trajStruct.size);
    udata.size=trajStruct.size;
    udata.color=trajStruct.color;
    udata.poly.x=[];
    udata.poly.y=[];
    udata.poly.h=[];
    udata.new_marks=[];
    udata.all_marks=[];
    udata.d_mark2=[];
    udata.fig1 = figure;
    set(udata.fig1,'KeyPressFcn','vis_trajgui([],''key'')',...
		   'Name','Primary Data');
    
    
    %%%%%%%%%%%%%%%%%%%%
    %
    % making the 'Tools' -menu
    %
    
    udata.m_i=uimenu(udata.fig1,'Label','Trajectoy T&ools');
    udata.m_i(2)=uimenu(udata.m_i,'Label','&Remove Trajectory',...
			'Callback',...
			'vis_trajgui([],''remove_traj'')');
    udata.m_i(3)=uimenu(udata.m_i(1),'Label','&Dye Nodes',...
			'Callback',...
			'vis_trajgui([],''dye_gui'')');
    udata.m_i(4)=uimenu(udata.m_i(1),'Label','&Clear Markers',...
			'Callback',...
			'vis_trajgui([],''clear'')');
    udata.m_i(5)=uimenu(udata.m_i(1),'Label','&Save',...
			'Callback',...
			'vis_trajgui([],''save'')');
    udata.m_i(6)=uimenu(udata.m_i(1),'Label','&Load',...
			'Callback',...
			'vis_trajgui([],''load'')');

    %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
    %
    % drawing data components to the figure ....
    %
    % 
    
    if nargin < 5 | isempty(comps) | (isstr(comps) & strcmp(comps,'all'))
      comps = 1:size(data1,2);
    end
    
    x=1:size(data1,1);
    
    for i=1:length(comps)
      subplot(length(comps),1,i);
      udata.h(i)=gca;
      udata.d_mark(i).h=[];
      
      udata.d(i)=plot(x,data1(:,comps(i)),...
		      'ButtonDownFcn',...
		      'vis_trajgui([],''line_down'')');        
      set(gca,'XLim',[1 size(data1,1)],...
	      'XTick',[],...	      
	      'ButtonDownFcn','vis_trajgui([],''line_down'')'); %,...
	      %'YLim',[min(data1(:,comps(i))) max(data1(:,comps(i)))]);
	      
      ylabel(names{comps(i)});
      hold on;
      ymin=get(udata.h(i),'YLim');
      pos=mean(get(udata.h(i),'XLim'));
      udata.l(i) = line([pos pos],[ymin(1) ymin(2)],...
			'Color','red',...
			'ButtonDownFcn',...
			'vis_trajgui([],''down'')');
    end
    udata.text1=[];
    
    udata.fig2=[];
    udata.h2=[];
    udata.l2=[];
    udata.text2=[];
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%%
    %
    % ... and to the figure 2.
    %
    
    
    if ~isempty(trajStruct.secondary_data)
      data2=trajStruct.secondary_data;
      if isempty(trajStruct.secondary_names)
	for i=1:size(data1,2)
	  names2{i,1}=sprintf('Var%d',i);
	end
      else
	names2=trajStruct.secondary_names;
      end
      
      udata.fig2 = figure;
      set(udata.fig2,'Name','Secondary Data');
      set(udata.fig2,'KeyPressFcn',...
		     'vis_trajgui([],''key'')');
      for i=1:size(data2,2)
	subplot(size(data2,2),1,i);
	udata.h2(i) = gca;
	udata.d_mark2(i).h=[];
	udata.d2(i) = plot(x,data2(:,i),...
			   'ButtonDownFcn',...
			   'vis_trajgui([],''line_down'')');
	set(gca,'XLim',[1 size(data1,1)],'XTick',[],'ButtonDownFcn',...
		'vis_trajgui([],[],[],[],[],[],''line_down'')');
	ylabel(names2{i});
	hold on;
	ymin = get(udata.h2(i),'YLim');
	pos = mean(get(udata.h2(i),'XLim'));
	udata.l2(i) = line([pos pos],ymin,'Color','red',...
			   'ButtonDownFcn','vis_trajgui([],''down'')');
      end
    end  
    
    %%%%%%%%%%%%%%%%%%%%%%%%%%
    
    set(udata.fig1,'UserData',udata);
    if ~isempty(udata.fig2);
      tmp.fig1=udata.fig1;
      set(udata.fig2,'UserData',tmp);
    end
    tmp=get(udata.sM_h,'UserData');
    tmp.fig1=udata.fig1;
    set(udata.sM_h,'UserData',tmp);  
    set_numbers(round(pos)); 
    return;
  end
  
end

%%%%%%%%%%%%%%%%%
%
% if figures have been drawn, the only function calls that may exist
% are the ones that change the state of the application.
%
udata=get(gcf,'UserData');
udata=get(udata.fig1,'UserData');

switch arg
 case 'fuzzy'
  fuzzy_traj(tS,[]);
  return;
 case 'move_fuzzy'
  fuzzy_traj([],'move');
  return;
 case 'remove_traj'
  remove_traj;
  return;
 case 'line_down'
  line_bdf('down');
  return;
 case 'line_drag'
  line_bdf('drag');
  return;
 case 'line_up'
  line_bdf('up');
  return;
 case 'dye_gui';
  color_gui(udata.fig1);
  return;
 case {'dye','cyan','magenta','yellow','red','green','blue','white','grey'}
  dye_nodes(arg);
  return;
 case 'clear'
  clear_markers;
  return;
 case 'key'
  key_bdf;
  return;   
 case 'click'
  click;
  return;
 case 'save'
  save_data;
  return;
 case 'load'
  load_data;
  return;
end


%%%%%%%%%%
%
% lines in the data figure(s) are dragged ...
%

udata=get(gcf,'UserData');
udata=get(udata.fig1,'UserData');


lims=get(gca,'XLim');
x = getfield(get(gca,'CurrentPoint'),{1});  % the location of the line

if x < lims(1)
  x=lims(1);
elseif x > lims(2)
  x=lims(2);
end

old = gcf;

switch arg
 case 'down',...
      
  % mouse button is pressed down above the line
  
  set(gcf,'WindowButtonMotionFcn','vis_trajgui([],''drag'')');
  set(gcf,'WindowButtonUpFcn','vis_trajgui([],''up'')');
  set(udata.l,'EraseMode','xor');
  delete(udata.text1);
  if ~isempty(udata.l2);
    set(udata.l2,'EraseMode','xor');
    delete(udata.text2);
  end
  set(gcf,'Pointer','crosshair');
  
 case 'drag'
  % change the location of the lines
  
  set(0,'CurrentFigure',udata.fig1);
  set(udata.l,'XData',[x x]);
  if ~isempty(udata.fig2)
    set(0,'CurrentFigure',udata.fig2);
    set(udata.l2,'XData',[x x]);
  end
  draw_traj(round(x));
  set(0,'CurrentFigure',old);
 case 'up'
  
  % draw trajectory and set figure to the normal state.
  
  set(udata.l,'EraseMode','normal');
  set(gcf,'Pointer','arrow','WindowButtonMotionFcn','',...
	  'WindowButtonUpFcn','');
  draw_traj(round(x));
  set_numbers(round(x));
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function draw_traj(point)

udata=get(gcf,'UserData');
udata=get(udata.fig1,'UserData');
color=udata.color;
eMode='normal';
if isstr(udata.color) & strcmp(udata.color,'xor')
  eMode='xor';
  color='black';
end
ind=udata.bmus(point);
[i j] = ind2sub(udata.msize,ind);
if ~mod(i,2)
  j=j+0.5;
end
old = gcf;
set(0,'CurrentFigure',udata.sM_h);
hold on;
if isempty(udata.traj) | length(udata.traj.h) ~= length(udata.a_h)
  
  % trajectory does not exist
  
  for i=1:length(udata.a_h)
    subplot(udata.a_h(i));
    hold on;
    new.h = plot(j,i,'Color',color,'EraseMode',eMode,'LineStyle','none');
    udata.traj.h(i)=new;
    udata.traj.j=j;
    udata.traj.i=i;
  end
else
  if length(udata.traj.j) == udata.length_of_traj
    % if the length of trajectory == ..., 
    udata.traj.j(1) = [];         % the first (the oldest) coordinate pair
    udata.traj.i(1) = [];         % is removed.
  end
  udata.traj.j=[udata.traj.j;j]; % the new point is added to the
  udata.traj.i=[udata.traj.i;i]; % end of coordinate vectors (i and j)
  for i=1:length(udata.a_h)
    subplot(udata.a_h(i));            % remove the existing trajectory
    delete(udata.traj.h(i).h);          % and plot the new one.
    for j=1:length(udata.traj.j)
      udata.traj.h(i).h(j)=plot(udata.traj.j(j),udata.traj.i(j),...
				'Color',color,...
				'EraseMode',eMode,'Marker','o','LineWidth',2,...
				'MarkerSize',udata.size(udata.length_of_traj-j+1),...
				'LineStyle','none');
    end
  end
end 
set(0,'CurrentFigure',udata.fig1);
set(udata.fig1,'UserData',udata);


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function set_numbers(x);

% This function writes the numbers beside of the pointer lines


udata=get(gcf,'UserData');
udata=get(udata.fig1,'UserData');

xlim = get(gca,'XLim');
ylim = get(gca,'YLim');
p = ylim(1) + 0.9*(ylim(2)-ylim(1));

old = gcf;
set(0,'CurrentFigure',udata.fig1);


for i=1:length(udata.h)
  subplot(udata.h(i));
  
  % check if the text is placed to the left side of the line...
  
  if abs(x-xlim(1)) > (abs(x-xlim(2)))
    udata.text1(i)=text(x-1,p,sprintf('%d ->',x),...
                        'VerticalAlignment','top',...
                        'HorizontalAlignment','right',...
                        'FontWeight','demi');
  else
    
    %  or to the right side.
    
    udata.text1(i)=text(x+1,p,sprintf('<- %d',x),...
			'VerticalAlignment','top',...
			'FontWeight','demi');
  end
end

if ~isempty(udata.fig2)
  set(0,'CurrentFigure',udata.fig2);
  
  for i=1:length(udata.h2)
    subplot(udata.h2(i));
    
    if abs(x-xlim(1)) > (abs(x-xlim(2)))
      udata.text2(i)=text(x-1,p,sprintf('%d ->',x),...
			  'VerticalAlignment','top',...
			  'HorizontalAlignment','right',...
			  'FontWeight','demi');
    else
      udata.text2(i)=text(x+1,p,sprintf('<- %d',x),...
			  'VerticalAlignment','top',...
			  'FontWeight','demi');
    end
  end
end

set(0,'CurrentFigure',old);
set(udata.fig1,'UserData',udata);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function remove_traj()

% delete trajectory -object from every component plane.

udata=get(gcf,'UserData');
udata=get(udata.fig1,'UserData');

if isempty(udata.traj)
  return;
end


for i=1:length(udata.traj.h)
  delete(udata.traj.h(i).h);
end

udata.traj=[];
set(udata.fig1,'UserData',udata);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function line_bdf(arg)

% this function takes care of action when region is selected in the
% data figure.


udata=get(gcf,'UserData');
udata=get(udata.fig1,'UserData');
xlim=get(gca,'XLim');

if ~(any(strcmp('THIS',fieldnames(udata))))
  p = getfield(get(gca,'CurrentPoint'),{1});
else
  p = getfield(get(udata.THIS,'CurrentPoint'),{1});
end

if p < xlim(1)
  p = xlim(1);
elseif p > xlim(2)
  p = xlim(2);
end



switch arg
 case 'down'
  
  % the mouse button is pressed down, the pointer lines are drawn.
  % and the state of the figure is changed.
  
  udata.THIS=gca;
  set(gcf,'WindowButtonMotionFcn',...
	  'vis_trajgui([],''line_drag'')',...
	  'WindowButtonUpFcn','vis_trajgui([],''line_up'')');
  udata.start_p=p;
  
  old = gcf;
  set(0,'CurrentFigure',udata.fig1);
  
  for i=1:length(udata.h)
    subplot(udata.h(i));
    udata.t_line.h(i)=line([p p],get(gca,'YLim'),'Color','red');
    udata.t_line.h2(i)=line([p p],get(gca,'YLim'),'Color','red',...
			    'EraseMode','xor');
  end
  if ~isempty(udata.h2)
    set(0,'CurrentFigure',udata.fig2);
    for i=1:length(udata.h2)
      subplot(udata.h2(i));
      udata.t_line2.h(i)=line([p p],get(gca,'YLim'),'Color','red');
      udata.t_line2.h2(i)=line([p p],get(gca,'YLim'),'Color','red',...
			       'EraseMode','xor');
    end
  end
  
 case 'drag'
  
  % change the position of the pointer lines
  
  old = gcf;
  set(0,'CurrentFigure',udata.fig1);
  set(udata.t_line.h2,'XData',[p p]);
  if ~isempty(udata.fig2)
    set(0,'CurrentFigure',udata.fig2);
    set(udata.t_line2.h2,'XData',[p p]);
  end
  set(0,'CurrentFigure',old);
 case 'up'
  
  
  % sort the 'points' -vector and draw the markers to the data and nodes
  
  points=sort([round(udata.start_p) round(p)]);
  draw_markers(points(1):points(2));
  udata=get(udata.fig1,'UserData');
  udata.new_marks=unique([udata.new_marks ;(points(1):points(2))']);