preprocess.m

很全的som工具箱 四个demo可为初学者提够帮助

%%% arg's first "word" is 4 letters long and it can be:
%%%
%%% 'key '
%%% 'down'
%%% 'drag'
%%% 'up  '

if strcmp(arg,'key ') %string is 'key' + 1 space!!!
  if ~LOG
    key=get(gcf,'CurrentCharacter');
   else 
    key=arg2
  end
  if ~strcmp(key,'<') & ~strcmp(key,'>')
    return;
  end
  data=get(gcf,'UserData');
  sel=data.selected_vects;
  if length(sel) == 1
    if strcmp(key,'<') & sel ~= 1 
      data.selected_vects=sel-1;
      set(gcf,'UserData',data);
     elseif strcmp(key,'>') & sel ~= length(data.sData.data(:,1))
      data.selected_vects = sel + 1;
      set(gcf,'UserData',data);
     end
  else
    if strcmp(key,'<') & sel(1) ~= 1
      data.selected_vects=cat(2,sel(1)-1,sel);
      set(gcf,'UserData',data);
     elseif strcmp(key,'>') & sel(length(sel)) ~= length(sel)
      data.selected_vects=cat(2,sel,sel(length(sel))+1);
      set(gcf,'UserData',data);
     end
  end
  cplot_mimema;
  pro_tools('plot_hist');
  pro_tools('c_stat');
  vects=zeros(1,length(data.sData.data(:,1)));
  vects(data.selected_vects)=1;
  draw_vectors(vects,data.vector_h);
 
  if ~LOG
    data=get(gcf,'UserData');
    data.LOG{length(data.LOG)+1}=...
    sprintf('preprocess(''vector_bdf'',''key  %s'');',key);
                                                %string is 'key'+2spaces+%s
    set(gcf,'UserData',data);
  end
  return;
end

switch arg
  case 'down'
   set(gcf,'WindowButtonUpFcn','preprocess(''vector_bdf'',''up  '')');
   set(gcf,'WindowButtonMotionFcn','preprocess(''vector_bdf'',''drag'')');
   switch get(gcf,'SelectionType')
     case 'normal'
      data.lims1=round(getfield(get(gca,'CurrentPoint'),{1,1}));
      data.lims2=[];
     case 'alt'
      tmp=round(getfield(get(gca,'CurrentPoint'),{1,1}));
      if isempty(get(gca,'UserData'))
        data.lims1=tmp;
        data.lims2=[];
      else
        data.lims1=cat(2,getfield(get(gca,'UserData'),'lims1'),tmp);
        data.lims2=getfield(get(gca,'UserData'),'lims2');
      end
   end
   coords=get(gca,'CurrentPoint');
   h=line([coords(1),coords(1)],get(gca,'YLim'),'EraseMode','xor');
   set(h,'Color','red');
   h2=line([coords(1),coords(1)],get(gca,'YLim'),'EraseMode','xor');
   set(h2,'Color','red');
   data.h=h;
   data.h2=h2;
   set(gca,'UserData',data);

  case 'drag'
   coords=get(gca,'CurrentPoint');
   lim=get(gca,'XLim');
   h2=getfield(get(gca,'UserData'),'h2');
   if lim(1) >= coords(1)
     set(h2,'XData',[lim(1) lim(1)]);
   elseif lim(2) <= coords(2)
     set(h2,'XData',[lim(2) lim(2)]);
   else
     set(h2,'XData',[coords(1) coords(1)]);
   end
  case 'up  '   % string is 'up' + 2 spaces!!! 
   set(gcf,'WindowButtonUpFcn','');
   set(gcf,'WindowButtonMotionFcn','');
   if ~LOG
     data=get(gca,'UserData');
     delete(data.h);
     delete(data.h2);
     tmp=round(getfield(get(gca,'CurrentPoint'),{1,1}));
     data.lims2=cat(2,data.lims2,tmp);
     tmp_data=sort(cat(1,data.lims1,data.lims2));
     high=getfield(get(gca,'XLim'),{2})-1;
     vectors=zeros(1,high);
     tmp_data(find(tmp_data<1))=1;
     tmp_data(find(tmp_data>high))=high;

     for i=1:getfield(size(tmp_data),{2})
       vectors(tmp_data(1,i):tmp_data(2,i))=1;
     end
     selected_vects=find(vectors);
   else
     pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
     len=size(getfield(getfield(get(pre_h,'UserData'),'sData'),'data'));
     vectors=zeros(1,len(1));
     i=1;
     while i <= length(arg2) & (isspace(arg2(i)) | ~isletter(arg2(i)))
      i=i+1;
     end
     arg3=arg2(i:length(arg2));
     selected_vects=str2num(arg2(1:i-1));
     if ~isempty(arg3) & ~all(isspace(arg3))
       selected_vects=unique(cat(2,selected_vects,...
                            getfield(get(pre_h,'UserData'),'selected_vects')));
     end           
     vectors(selected_vects)=1;  
     set(pre_h,'CurrentAxes',getfield(get(pre_h,'UserData'),'vector_h'));
     set(0,'CurrentFigure',pre_h);
   end
   draw_vectors(vectors,gca);
   sData=getfield(get(gcf,'UserData'),'sData');
   h=getfield(get(gcf,'UserData'),'vect_mean_h');
   vect_means(sData,h,selected_vects);
   if ~LOG
     set(gca,'UserData',data);
   end
   data=get(gcf,'UserData');  
   data.undo.sData=data.sData;
   data.undo.selected=data.selected_vects;
   data.selected_vects=selected_vects;
   if ~LOG
     data.LOG{length(data.LOG)+1}='% Vector selection by using the mouse...';
     tmp=sprintf('preprocess(''vector_bdf'',''up   %s'');',...
                num2str(data.selected_vects));
     if length(tmp) > 500
       tmp=textwrap({tmp},500);
       data.LOG{length(data.LOG)+1}=cat(2,tmp{1},''');');
       for i=2:length(tmp)-1
         data.LOG{length(data.LOG)+1}=...
               cat(2,sprintf('preprocess(''vector_bdf'',''up   %s',...
                     tmp{i}),'add'');');
       end
       data.LOG{length(data.LOG)+1}=...
             cat(2,sprintf('preprocess(''vector_bdf'',''up   %s',...
                   tmp{length(tmp)}(1:length(tmp{length(tmp)})-3)),' add'');');
     else
       data.LOG{length(data.LOG)+1}=tmp;
     end
   end   
   set(gcf,'UserData',data);
   cplot_mimema;
   sel_comp;
end

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

function sel_button(varargin) 

%SEL_BUTTON     A Callback function. It performs the operations needed
%               when vector components are selected.
%
% See also SEL_COMP
%

if nargin == 1
  LOG=1;
  pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
  string=getfield(get(pre_h,'UserData'),'comp_names_h');
  string=getfield(get(string,'String'),{str2num(varargin{1})});
  set(0,'CurrentFigure',pre_h);
else
  LOG=0;
  val=get(getfield(get(gcf,'UserData'),'comp_names_h'),'Value');
end

sel_button_h=getfield(get(gcf,'UserData'),'sel_button_h');
sel_comps_h=getfield(get(gcf,'UserData'),'sel_comps_h');
comp_names_h=getfield(get(gcf,'UserData'),'comp_names_h');
if ~LOG
  string=getfield(get(comp_names_h,'String'),{get(comp_names_h,'Value')});
end
tmp_string=get(sel_comps_h,'String');

if iscell(tmp_string)

  for i=1:length(string)
    if ~any(strcmp(string{i},tmp_string))
      tmp_string=cat(1,tmp_string,string(i));
    end
  end
  string=tmp_string;
end

set(sel_comps_h,'String',string);
set(comp_names_h,'Value',[]);
sel_comp;
if ~LOG
  data=get(gcf,'UserData');
  data.LOG{length(data.LOG)+1}='% Select components';
  data.LOG{length(data.LOG)+1}=sprintf('preprocess(''sel_button'',''%s'');',...
                                        num2str(val));
  set(gcf,'UserData',data);
end
 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function clear_button(varargin)

%CLEAR_BUTTON  Function callback evaluated when a 'Clear'-button is
%              pressed. It removes texts from the 'selected components' 
%              -window and the 'selected component data' -window and
%              clears the 'histogram' -axis.
%
%

if nargin==1
  LOG=1;
  pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
  set(0,'CurrentFigure',pre_h);
else
  LOG=0;
end

sel_comp_h=getfield(get(gcf,'UserData'),'sel_comp_h');
sel_cdata_h=getfield(get(gcf,'UserData'),'sel_cdata_h');
sel_cplot_h=getfield(get(gcf,'UserData'),'sel_cplot_h');
sel_chist_h=getfield(get(gcf,'UserData'),'sel_chist_h');
vector_h=getfield(get(gcf,'UserData'),'vector_h');

set(sel_comp_h,'Value',1);
set(sel_cdata_h,'String',' ');
subplot(sel_chist_h);
hold off;
cla;

selected=getfield(get(gcf,'UserData'),'selected_vects');
dims=size(getfield(getfield(get(gcf,'UserData'),'sData'),'data'));
vectors=zeros(1,dims(1));
vectors(selected)=1;
subplot(vector_h);
draw_vectors(vectors,vector_h);
if ~LOG
  data=get(gcf,'UserData');
  data.LOG{length(data.LOG)+1}='% Remove components from the selected list.';
  data.LOG{length(data.LOG)+1}='preprocess(''clear_button'',''foo'');';
  set(gcf,'UserData',data);
end

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

function sel_comp(varargin)

%SEL_COMP  performs the operations needed when vector components are
%          chosen. It writes maxima, minima, mean and standard deviation
%          of the chosen component to a text box window and draws a
%          histogram of the chosen component of selected vectors' 
%
%

pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
set(0,'CurrentFigure',pre_h);
sel_comp_h=getfield(get(pre_h,'UserData'),'sel_comp_h');

if nargin == 1
  set(sel_comp_h,'Value',str2num(varargin{1}));
elseif ~isempty(gcbo)
  no=get(sel_comp_h,'Value');
  data=get(gcf,'UserData');
  data.LOG{length(data.LOG)+1}='% Select one component';
  data.LOG{length(data.LOG)+1}=cat(2,'preprocess(''sel_comp'',''',...
                                      num2str(no),''');');
  set(gcf,'UserData',data);
end

pro_tools('c_stat');
pro_tools('plot_hist');
data=get(gcf,'UserData');
sData=data.sData;  
vector_h=data.vector_h;
len=length(sData.data(:,1));
vects=zeros(1,len);
vects(data.selected_vects)=1;
draw_vectors(vects,vector_h);

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

function cplot_mimema

global no_of_sc

sData=getfield(get(gcf,'UserData'),'sData');
sel_cplot_h=getfield(get(gcf,'UserData'),'sel_cplot_h');
selected=getfield(get(gcf,'UserData'),'selected_vects');

set(sel_cplot_h,'YLim',[0 length(sData.data(1,:))+1]);

subplot(sel_cplot_h);
hold off;
cla;
hold on;

for i=1:length(sData.data(1,:))
  Max=max(sData.data(:,i));
  Min=min(sData.data(:,i));
  tmp=sData.data(selected,i);

  selMax=max(tmp);
  selMin=min(tmp);
  Mean=abs(mean(tmp(find(~isnan(tmp)))));
  Median=abs(median(tmp(find(~isnan(tmp)))));
  
  if Max ~= Min & ~all(isnan(sData.data(:,i)))

    if rem(i,no_of_sc)   % no_of_sc is defined in the beginning of this file...

      line([abs(selMin-Min)/(Max-Min) (selMax-Min)/(Max-Min)],...
           [i i],'Color','black');
      plot(abs(Mean-Min)/(Max-Min),i,'oblack');
      plot(abs(Median-Min)/(Max-Min),i,'xblack');
    else
      line([abs(selMin-Min)/(Max-Min) (selMax-Min)/(Max-Min)],...
           [i i],'Color','black','LineWidth',2);
      plot(abs(Mean-Min)/(Max-Min),i,'oblack','LineWidth',2);
      plot(abs(Median-Min)/(Max-Min),i,'xblack','LineWidth',2);
    end
  else

    if rem(i,no_of_sc)     % N is defined in the beginning of this file.

      plot(mean(get(gca,'XLim')),i,'oblack');
      plot(mean(get(gca,'XLim')),i,'xblack');
    else
      plot(mean(get(gca,'XLim')),i,'oblack','LineWidth',2);
      plot(mean(get(gca,'XLim')),i,'xblack','LineWidth',2);
    end
  end
end

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


function bool=set_sD_stats

%SET_SD_STATS Writes the data set names to popup menu.
%
%

bool=1;
data=get(gcf,'UserData');

for i=1:length(data.sD_set)
 % if ~isvalid_var_name({data.sD_set(i).name})
 %   close_preprocess;
 %   bool=0;
 %   return;
 % end
  string{i}=cat(2,sprintf('#%d: ',i),data.sD_set(i).name);
end

set(data.sD_set_h,'String',string);
data.sData=data.sD_set(get(data.sD_set_h,'Value'));
data.sData.MODIFIED=0;
data.sData.INDEX=1;
set(gcf,'UserData',data);
write_sD_stats;

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

function write_sD_stats

%WRITE_SD_STATS  writes data's name, length and dimension to text box.
%
%

pre_h=findobj(get(0,'Children'),'Tag','Preprocess');


sD_name_h=getfield(get(pre_h,'UserData'),'sD_name_h');
sData=getfield(get(pre_h,'UserData'),'sData');
dims=size(sData.data);
string{1}=cat(2,'Name:   ',sData.name);
string{2}=cat(2,'Length: ',sprintf('%d',dims(1)));
string{3}=cat(2,'Dim:     ',sprintf('%d',dims(2)));

set(sD_name_h,'String',string);
set(sD_name_h,'HorizontalAlignment','left');

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



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

function sel_sD(varargin)

%SEL_SD  sets new data to UserData's 'sData'.
%        
%

if nargin==1
  LOG=1;
  index=str2num(varargin{1});
  pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
  set(0,'CurrentFigure',pre_h);
else
  LOG=0;
end

sD_set_h=getfield(get(gcf,'UserData'),'sD_set_h');
comp_names_h=getfield(get(gcf,'UserData'),'comp_names_h');
vector_h=getfield(get(gcf,'UserData'),'vector_h');
vect_mean_h=getfield(get(gcf,'UserData'),'vect_mean_h');

if ~LOG
  index=get(sD_set_h,'Value');
end
data=get(gcf,'UserData');
data.undo = [];
INDEX=data.sData.INDEX;
data.sData=rmfield(data.sData,'MODIFIED'); 
data.sData=rmfield(data.sData,'INDEX');

tmp=data.sD_set(index);
tmp.MODIFIED=0;
tmp.INDEX=index;
data.sD_set(INDEX)=data.sData;
data.sData=tmp;

len=getfield(size(tmp.data),{1});

data.selected_vects=find(ones(1,len));
if ~LOG
  data.LOG{length(data.LOG)+1}='% Select a new data set.';
  data.LOG{length(data.LOG)+1}=sprintf('preprocess(''sel_sD'',''%s'');',...
                                        num2str(index));
end
set(gcf,'UserData',data);
write_sD_stats;
set_compnames(tmp,comp_names_h);
draw_vectors(ones(1,len),vector_h);
vect_means(tmp,vect_mean_h,data.selected_vects);
clear_button;
sel_comp;
cplot_mimema;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

function indices=get_indices

pre_h=findobj(get(0,'Children'),'Tag','Preprocess');

comp_names_h=getfield(get(pre_h,'UserData'),'comp_names_h');
indices = get(comp_names_h,'Value');

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



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

function sD_management(varargin)

if nargin ~= 1
  pre_h=findobj(get(0,'Children'),'Tag','Preprocess');
  preh_udata=get(pre_h,'UserData');
  preh_udata.LOG{length(preh_udata.LOG)+1}=...
           '% Starting the ''Data Set Management'' -window...';
  preh_udata.LOG{length(preh_udata.LOG)+1}=...
                 'preprocess(''sD_management'',''foo'');';