som_show_add.m

此源码提供了ssvm的完整的matlab源代码

	       mode_ '''.']);
  else 
    Property=setfield(Property,identifier,value);
  end
end

% set default subplot
if isempty(Property.subplot)
  % search the subplot number for current axis
  value=find(gca==handles);    
  if isempty(value) | value>length(handles) 
    error('SubPlot default value setting: current axis is not in the figure!');
  else
    Property.subplot=value;
  end
end
  
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%% Main switch: select the right subfunction %%%%%%%%%%%%%%%%%%%

switch mode
 case 'hit'
  h_=hit(D, lattice, msize, Property);  
 case 'mhit'
  h_=mhit(D, lattice, msize, Property);  
 case 'label'
  h_=label(D, lattice, msize, Property);
 case 'traj'
  h_=traj(D, lattice, msize, Property);
 case 'comet'
  %error('Not yet implemented.'); 
  h_=comet(D, lattice, msize, Property);
 otherwise
  error('Whoops! Internal error: unknown mode!');
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Build output if necessary %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

if nargout>0
  h=h_;
end

%%%% SUBFUNCTIONS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function h_=hit(Hits, lattice, msize, Property);  

% number of map units
munits=prod(msize);

% subplots
p=Property.subplot;
handles=Property.handles;

% Set default marker
if isempty(Property.marker),
  if strcmp(Property.text,'on')
    Property.marker='none';
  else
    Property.marker='lattice';
  end
end

% Set default markersize
if isempty(Property.markersize)
  if strcmp(Property.marker,'none'),
    warning('MarkerSize is not meaningful since Marker is set to ''none''.');
  elseif strcmp(Property.marker,'lattice'),
    Property.markersize=1; % normalized size
  else
    Property.markersize=12; % points
  end
end

% Set default colors
if ~isempty(Property.markercolor),
  if strcmp(Property.marker,'none')
    warning('MarkerColor is not used since Marker is set to ''none''.');
    Property.markercolor=[]; % not used
  else
    ; % ok
  end
elseif any(strcmp(Property.marker,{'+','*','.','x'})),
  % these don't use fill color: 'none' will cause default
  % edgecolor to be 'k'.
  Property.markercolor='none'; 
else
  Property.markercolor='k';
end

if ~isempty(Property.edgecolor),
  if strcmp(Property.marker,'none')
    warning(['EdgeColor is not used since Marker is set to' ...
	     ' ''none''.']);
  else
    ; %ok 
  end
elseif ~strcmp(Property.markercolor,'none'),
  Property.edgecolor='none';
else
  Property.edgecolor='k';
end

% Set default text
if isempty(Property.text),
  Property.text='off';
end

% Set default textsize
if isempty(Property.textsize)
  Property.textsize=10;
elseif strcmp(Property.text,'off')  
  warning('TextSize not used as hits are not set to be shown as numbers.');
end

% Set default textcolor
if isempty(Property.textcolor)
  Property.textcolor='w';
elseif strcmp(Property.text,'off')  
  warning('TextColor not used as hits are not set to be shown as numbers.');
end

%% Action %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

h_=[];          % this variable is for collecting the object handles

% Select the drawing mode 

if ~strcmp(Property.marker,'none') 
  
  %%%%% Draw spots %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  
  % unit coordinates 
  coord=som_vis_coords(lattice,msize);
  
  % Calculate the size of the spots
  mx=max(Hits);
  
  if mx==0,
    % nothing to draw!
    h_=[]; 
    return
  else
    Size=sqrt(Hits./mx);      
  end
  % coordinates for non-zero hits (only those are drawn)
  coord=coord(Size~=0,:);
  Size=Size(Size~=0);
  N=size(Size,1);
  
  % som_cplane can't draw one unit with arbitrary
  % coordinates as it its mixed with msize:
  if size(coord,1)==1 & strcmp(Property.marker,'lattice'),
    Size=[Size;Size];
    coord=[coord;coord];
  end
  
  for i=1:length(p),
    % Set axes
    axes(handles(p(i)));
    % Get hold state and caxis
    memhold=ishold; cax=caxis;     
    hold on;
    
    switch Property.marker
     case 'lattice'
      h_(i,1)=som_cplane(lattice, coord, Property.markercolor, ...
			 Property.markersize*Size);
     otherwise  
      [S,m]=som_grid(lattice, [N 1],...
		     'Coord',coord, ...
		     'Line','none',...
		     'Marker',Property.marker,...
		     'MarkerColor',Property.markercolor,...
		     'MarkerSize', Size*Property.markersize);
      h_=[h_;m(:)];
    end
    
    % Restore hold state
    if ~memhold         
      hold off;
    end
  end
  
  % Set edgecolor
  if strcmp(Property.marker,'lattice')
    set(h_,'edgecolor',Property.edgecolor);
  else
    set(h_,'markeredgecolor',Property.edgecolor);
  end
end
  
if strcmp(Property.text,'on'),
  %%%%% Draw numbers %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%    

  % Do numbers
  Hits=reshape(Hits,[munits 1]);
  labels=cell([munits 1]);
  for i=1:length(Hits)
    if Hits(i)              % zero hit won't be shown
      labels(i)={num2str(Hits(i))};
    end
  end

  for i=1:length(p),
    axes(handles(p(i)));    % Set axes
    memhold=ishold;         % Get hold state
    hold on;
    [S,m,l,t]=som_grid(lattice, msize, ...
		       'Line','none',...
		       'Marker','none', ...
		       'Label',labels, ...
		       'LabelColor', Property.textcolor, ...
		       'LabelSize', Property.textsize);
    % Get handles
    h_=[h_;t(:)];
    
    % Restore hold state and caxis
    if ~memhold     
      hold off;
    end
    caxis(cax);
  end

  % Remove zero object handles (missing objects)
  h_=setdiff(h_,0);
end

%% Set object tags (for som_show_clear) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%

set(h_,'Tag','Hit')

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

function h_=mhit(Hits, lattice, msize, Property);  

% number of map units
munits=prod(msize);

% subplots
p=Property.subplot;
handles=Property.handles;


% Set default marker
if isempty(Property.marker),
  Property.marker=lattice;
end

% variable 'mode' indicates which kind of markers are used:

if iscell(Property.marker),
  mode='marker';
elseif vis_valuetype(Property.marker,{'markerstyle'}),
  mode='marker';
elseif strcmp(Property.marker,'pie'),
  mode='pie';
else
  mode='lattice';
end

% Set default size scaling
if isempty(Property.sizefactor)
  Property.sizefactor='separate';
end

% Set default markersize 
if isempty(Property.markersize)
  if any(strcmp(mode,{'lattice','pie'})),
    Property.markersize=1; % normalized
  else
    Property.markersize=12;  % points
  end
end

% Set default colors

if isempty(Property.markercolor),
  Property.markercolor=hsv(size(Hits,2));
end

if isempty(Property.edgecolor),
  if vis_valuetype(Property.markercolor,{'none'}),
    Property.edgecolor='k';
  else
    Property.edgecolor='none';
  end
end

% Set default shift
if isempty(Property.shift)
  Property.shift=0;
end

%% Action %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

h_=[];          % this variable is for collecting the object handles

switch mode
  case {'marker','lattice'}
   % Number of hits histograms
   n_Hits=size(Hits,2);
   % Calculate the size of the spots
   
   if strcmp(Property.sizefactor,'common')
     mx=max(max(Hits));
     if mx==0 % nothing to draw!
       h_=[]; return
     end
     spotSize=sqrt(Hits./mx);
   else
     mx=repmat(max(Hits),munits,1);
     mx(mx==0)=1; % Prevent division by zero
     spotSize=sqrt(Hits./mx);
   end
   
   %%% Make spotSize
   
   %reshape Size to a vector [spotSizeforHist(:,1); spotSizeforHist(:,2);...]
   spotSize=spotSize(:);
   
   % indices for non-zero hits (only those are drawn)
   notZero=find(spotSize ~= 0);
   
   % Drop zeros away from spotSize
   spotSize=spotSize(notZero);
   
   % Order spots so that bigger will be drawn first, so that they 
   % won't hide smaller ones
   [dummy, sizeOrder]=sort(spotSize); sizeOrder=sizeOrder(end:-1:1);
   spotSize=spotSize(sizeOrder);
   
   %%% Make unit coordinates 
   coord=som_vis_coords(lattice,msize);
   
   move=repmat(linspace(-.1,.1,n_Hits),size(coord,1),1)*Property.shift;
   move=repmat(move(:),1,2);
   
   % do n_Hits copies of unit coordinates so that they match spotSize
   coord=repmat(coord,n_Hits,1)+move;
   
   % Drop zeros away from coords and order
   coord=coord(notZero,:);
   coord=coord(sizeOrder,:);
   
   %%% Make unit colors
   
   if vis_valuetype(Property.markercolor,{'nx3'}),
     % If multiple colors Copy unit colors so that they match spotSize
     color=Property.markercolor(reshape(repmat([1:n_Hits]',1,munits)',...
					munits*n_Hits,1),:);
     % drop zeros away & order
     color=color(notZero,:);
     color=color(sizeOrder,:);
   else
     % only on color
     color=Property.markercolor;
   end
   
   %%% Make unit markers
   
   if iscell(Property.marker),
     %marker shows class: 
     marker=char(Property.marker);
     marker=marker(reshape(repmat([1:n_Hits]',1,munits)',...
			   munits*n_Hits,1),:);
     % Drop zeros, order & make to cell array (for som_grid)
     marker=marker(notZero,:);
     marker=cellstr(marker(sizeOrder,:));
   else
     marker=Property.marker;
   end

   % som_cplane can't draw one unit with arbitrary
   % coordinates as it its mixed with msize:
   if size(coord,1)==1 & strcmp(mode,'lattice'),
     spotSize = [spotSize; spotSize];
     coord = [coord; coord];
   end

   N=length(notZero); % for som_grid visuzalization routine
 case 'pie'
  % marker 'pie' requires size parameter totHits
  if strcmp(mode,'pie')
    coord=som_vis_coords(lattice, msize);
    notZero=sum(Hits,2)>0;
    Hits=Hits(notZero,:);
    coord=coord(notZero,:);
    N=size(notZero,1);
    totHits=sqrt(sum(Hits,2)./max(sum(Hits,2)));
  end
  
  % som_pieplane can't draw one unit with arbitrary
  % coordinates as it its mixed with msize:
  if size(coord,1)==1,
    Hits= [Hits; Hits];
    coord = [coord; coord];
  end
 otherwise
  error('Whoops: internal error. Bad mode in subfunction mhit');
end

for i=1:length(p),    %%% Main loop begins