routine.m

等高线拟合源代码（matlab源码）

function routine(action)

% ROUTINE Handles different subfunctions of GUI.

% All the software included in this package is presented as is.
% It may be distributed freely. The author can, however, not be
% held responsible for any problems whatever.
% 
% Designed by Johan Baeten.
% Last updated: 22-03-2000
% [e-mail]

% Put stringvariables like action(1:9) at the end 
% -> ohterwise an error will occur

global MENUBARHANDLES 
global SUBMENUHANDLES
global FITHANDLES
global UICTRHANDLES
global XSIZE
global YSIZE
global IMAGENAME
global NEWIMAGENAME
global MAINFIG
global VIEWSIZE
global NUMBER   % only in this m-file
global FILENAME
global FILECHANGED
global NR_OF_FITS;


maxnrofpoints = 12;



setting('showtitle');

if strcmp(action,'selectpoints'),
   if strcmp(IMAGENAME,'noname.pgm'),
       seterror(' There is no image presently loaded. First open an image, then try again.');
   else
       set(FITHANDLES,'visible','on');
       set(SUBMENUHANDLES(6),'enable','off');
       seterror([' Make a list of points using ''Clear List'', ''Add'' and ''Delete''.',...
                 ' Then press ''Fit'' to see a fit trough these points.']);
   end


elseif strcmp(action,'done'),
   seterror(' ... Closing interactive fit controls.');
   set(FITHANDLES(1:4,:),'visible','off');
   set(SUBMENUHANDLES(6),'enable','on');


elseif strcmp(action,'new'),
   if get(FITHANDLES(1,12),'value')==0,
     seterror(' ... List of points was already empty.');     
   else 
     seterror(' ... Clearing list of points.');
     set(FITHANDLES(3:4,1:maxnrofpoints),'string','','value',-1);
     set(FITHANDLES(1,12),'value',0);
     FILENAME = '*.mat';
     FILECHANGED = 0;
     setting('showtitle');
     hpunten = findobj(get(gca,'children'),'type','line');
     delete(hpunten); 
   end


elseif strcmp(action,'add'),
   nrofpoints = get(FITHANDLES(1,12),'value');
   if (nrofpoints < maxnrofpoints), 
      setting('menubarenableoff');
      seterror(' Select a point in the image or press the number in the list');
      set(FITHANDLES(1,1:7),'enable','off');
      set(FITHANDLES(1,17:18),'enable','off');
%      set(UICTRHANDLES(2:4),'enable','off');
      set(FITHANDLES(1,8), 'string','Stop Adding', ...
                'callback','routine(''stopadding'')');
      set(FITHANDLES(2,nrofpoints+1),'enable','on');
      set(MAINFIG,'pointer','crosshair');
      set(MAINFIG,'Windowbuttondownfcn','routine(''voegpunttoe'')');
   else
     routine('stopadding');
     seterror(' Maximum number of points reached! Use ''New'' or ''Delete''.');

   end
  


elseif strcmp(action,'voegpunttoe'),
   nrofpoints = get(FITHANDLES(1,12),'value');
   [x_size,y_size]=size(get(MAINFIG,'userdata')); 
   p=get(gca,'currentpoint');
   x=p(1,1); y = p(1,2);
   xl = get(gca,'xlim'); yl = get(gca,'ylim');
   % range for pixels from 1 to 128 e.g pixels 64 lies just befor the index 64  
   if ((x >= xl(1)) & (x<= xl(2)) & (y >= yl(1)) & (y <= yl(2))),
      nrofpoints = nrofpoints + 1;
      set(FITHANDLES(2,nrofpoints),'enable','off');
      set(FITHANDLES(1,12),'value',nrofpoints);
      x=x-0.5; y=y-0.5;
      set(FITHANDLES(2,nrofpoints),'enable','off');
      set(FITHANDLES(3,nrofpoints),'string',num2str(x),'value',x);
      set(FITHANDLES(4,nrofpoints),'string',num2str(y),'value',y);
      routine('plotpunten');
      if (nrofpoints < maxnrofpoints), 
        set(FITHANDLES(2,nrofpoints+1),'enable','on');
      else
        routine('stopadding');
        seterror(' Maximum number of points reached.');
      end
      FILECHANGED = 1;
   end
   

elseif strcmp(action,'plotpunten'),
   nrofpoints = get(FITHANDLES(1,12),'value'); 
   if nrofpoints > 0,
     for i = 1:nrofpoints,
        xc(i) = get(FITHANDLES(3,i),'value')+0.5;
        yc(i) = get(FITHANDLES(4,i),'value')+0.5;
     end
     hpunten = findobj(get(gca,'children'),'marker','o');
     l= length(hpunten);
     if l == 0,
       % no plotted points found
       hold on, plot(xc,yc,'ro'), hold off;
     else
       if l > 1, 
         seterror(' Message: More handles found then needed, taking first one ...');
       end
       set(hpunten(1),'xdata',xc,'ydata',yc); 
     end
   end


elseif strcmp(action,'stopadding'),
   set(MAINFIG,'Windowbuttondownfcn','','pointer','arrow');
   setting('menubarenableon'); 
   set(FITHANDLES(2,1:maxnrofpoints),'enable','off');
   set(FITHANDLES(1,1:7),'enable','on');
   set(FITHANDLES(1,17:18),'enable','on');
%  set(UICTRHANDLES(2:4),'enable','on');
   set(FITHANDLES(1,8), 'string','Close', ...
       'callback','routine(''done'')');
%   if strcmp(get(SUBMENUHANDLES(2),'label'),'Zoom Off'),
%      setting('togglezoomfunction');
%   end


elseif strcmp(action,'delete'),
   seterror(' Point to a previous selected point in the image or press a number in the list to be deleted');
   setting('menubarenableoff');
   set(FITHANDLES(1,1:7),'enable','off');
   set(FITHANDLES(1,17:18),'enable','off');
%   set(UICTRHANDLES(2:4),'enable','off');
   set(FITHANDLES(1,8), 'string','Cancel Delete','callback','routine(''canceldelete'')');
   nrofpoints = get(FITHANDLES(1,12),'value');
   if (nrofpoints == 0),
     seterror(' There is nothing to delete!  Press ''Cancel Delete''');
   else
     set(FITHANDLES(2,1:nrofpoints),'enable','on');
     set(MAINFIG,'WindowButtonDownFcn','routine(''deletebypointing'')');
   end

elseif strcmp(action,'canceldelete'),
   routine('stopadding');


elseif strcmp(action,'deletebypointing'),
   nrofpoints = get(FITHANDLES(1,12),'value');
   p=get(gca,'currentpoint');
   x=p(1,1); y = p(1,2);
   xl = get(gca,'xlim'); yl = get(gca,'ylim');
   if ((x >= xl(1)) & (x<= xl(2)) & (y >= yl(1)) & (y <= yl(2))),
     x=x-0.5; y=y-0.5;
     for i = 1:nrofpoints,
       xp(i) = get(FITHANDLES(3,i),'value');
       yp(i) = get(FITHANDLES(4,i),'value');
     end
     marge = (xl(2)-xl(1))/50;
     nr = 0;
     for i = 1 : nrofpoints,
       if (abs(x-xp(i)) < marge) & (abs(y-yp(i)) < marge),
       nr = i;
       end
     end
     if nr == 0,
       seterror(' No point found in this area.');
     else
       hulpstr = ['point',num2str(nr)];
       routine(hulpstr);
     end
     FILECHANGED = 1;
   end
   
   
elseif strcmp(action,'move'),
   nrofpoints = get(FITHANDLES(1,12),'value');
   if (nrofpoints == 0),
     seterror(' There is nothing to edit or move!');
   else
     setting('menubarenableoff');
     set(FITHANDLES(1,1:7),'enable','off');
     set(FITHANDLES(1,17:18),'enable','off');
%     set(UICTRHANDLES(2:4),'enable','off');
     set(FITHANDLES(1,8), 'string','Stop Editing','callback','routine(''stopadding'')');
     seterror(' Edit a point in the list or move a point in the image by ''dragging''.');
     set(FITHANDLES(2,1:nrofpoints),'enable','on');
     set(MAINFIG,'pointer','arrow',...
                 'WindowbuttonUpfcn','',...
                 'Windowbuttondownfcn','routine(''beweegpunt'')');
   end


elseif strcmp(action,'beweegpunt'),
   nrofpoints = get(FITHANDLES(1,12),'value');
   p=get(gca,'currentpoint');
   x=p(1,1); y = p(1,2);
   xl = get(gca,'xlim'); yl = get(gca,'ylim');
   if ((x >= xl(1)) & (x<= xl(2)) & (y >= yl(1)) & (y <= yl(2))),
     x=x-0.5; y=y-0.5;
     for i = 1:nrofpoints,
       xp(i) = get(FITHANDLES(3,i),'value');
       yp(i) = get(FITHANDLES(4,i),'value');
     end
     marge = (xl(2)-xl(1))/50;
     nr = 0;
     for i = 1 : nrofpoints,
       if (abs(x-xp(i)) < marge) & (abs(y-yp(i)) < marge),
       nr = i;
       end
     end
     if nr == 0,
       seterror(' No point found in this area.');
     else
       NUMBER = nr;
       set(MAINFIG,'pointer','cross');
       hold on, plot(xp(nr)+0.5,yp(nr)+0.5,'rx'),hold off;
       cbstr =['routine(''wijzigpunt'');routine(''plotpunten'');routine(''move'')'];
       set(MAINFIG,'WindowButtonUpFcn',cbstr);