plot.m
上传用户:haiyisale
上传日期:2013-01-09
资源大小:3246k
文件大小:4k
源码类别:

波变换

开发平台:

Matlab

  1. function varargout = plot(t,varargin);
  2. %PLOT Plot WVTREE object.
  3. %   PLOT(T) plots the WVTREE object T.
  4. %   FIG = PLOT(T) returns the handle of the figure, which
  5. %   contains the tree T.
  6. %   PLOT(T,FIG) plots the tree T in the figure FIG, which
  7. %   already contains a tree.
  8. %
  9. %   PLOT is a graphical tree-management utility. The figure
  10. %   that contains the tree is a GUI tool. It lets you change
  11. %   the Node Label to Depth_Position or Index, and Node Action
  12. %   to Split-Merge or Visualize.
  13. %   The default values are Depth_Position and Visualize.
  14. %
  15. %   You can click the nodes to execute the current Node Action.
  16. %
  17. %   After some split or merge actions you can get the new tree
  18. %   using the handle of the figure, which contains it.
  19. %   You must use the following special syntax:
  20. %       NEWT = PLOT(T,'read',FIG).
  21. %   In fact, the first argument is dummy. Then the most general
  22. %   syntax for this purpose is:
  23. %       NEWT = PLOT(DUMMY,'READ',FIG);
  24. %   where DUMMY is any object parented by an NTREE object.
  25. %
  26. %   DUMMY can be any object constructor name, which returns
  27. %   an object parented by an NTREE object. For example:
  28. %      NEWT = PLOT(ntree,'read',FIG);
  29. %      NEWT = PLOT(dtree,'read',FIG);
  30. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Feb-2003.
  31. %   Last Revision: 21-May-2003.
  32. %   Copyright 1995-2004 The MathWorks, Inc.
  33. %   $Revision: 1.1.6.2 $  $Date: 2004/03/15 22:38:53 $ 
  34. nbin = length(varargin);
  35. fig_tree = NaN;
  36. switch nbin
  37.   case 0    , option = 'create';
  38.   otherwise
  39.     option = varargin{1};
  40.     if isnumeric(option)
  41.         fig_tree = option;
  42.         option = 'create';
  43.     end
  44. end
  45. switch option
  46.   case 'create'
  47.   case {'length','type'} , fig_tree = varargin{2};  % Added Node Labels
  48.   case 'Reconstruct'     , fig_tree = varargin{2};  % Added Node Action
  49.   otherwise
  50.       if ischar(option) , % Added Node Labels
  51.           fig_tree = varargin{2}; 
  52.       end
  53. end
  54. switch option
  55.   case 'create'
  56.     fig_tree = plot(t.dtree,fig_tree);
  57.     
  58.     if nargout>0 , varargout{1} = fig_tree; end
  59. set(fig_tree,'NumberTitle','Off','Name',...
  60.         ['Fig ' int2str(fig_tree) ...
  61.  ' - Wavelet Tree Object - DWT extension mode: ',get(t,'extMode')]);
  62.     % Store the WVTREE.
  63.     %------------------
  64.     plot(dtree,'write',fig_tree,t);
  65.     % Add Node Label menus.
  66.     %----------------------
  67.     plot(ntree,'addNodeLabel',fig_tree,'length');
  68.     plot(ntree,'addNodeLabel',fig_tree,'type');
  69.     % Add Node Action menu.
  70.     %----------------------
  71.     plot(ntree,'addNodeAction',fig_tree,'Reconstruct');
  72.     % Set default Node Label to 'Index'.
  73.     %-----------------------------------
  74.     plot(ntree,'setNodeLabel',fig_tree,'Index');
  75.   case {'length','type'}
  76.     t = plot(ntree,'read',fig_tree);
  77.     if nbin<3 , nodes = allnodes(t); else , nodes = varargin{3}; end
  78.     n = length(nodes);
  79.     switch option
  80.       case 'length' , labtype = 's';
  81.       case 'type'   , labtype = 't';
  82.     end
  83.     labels = tlabels(t,labtype,nodes);
  84.     err    = ~isequal(n,size(labels,1));
  85.     varargout = {labels,err};
  86.   case 'Reconstruct'
  87.     node = plot(ntree,'getNode',fig_tree);
  88.     if isempty(node) , return; end
  89.     t = plot(ntree,'read',fig_tree);
  90.     axe_vis = plot(ntree,'getValue',fig_tree,'axe_vis');
  91.     mousefrm(fig_tree,'watch')
  92.     x = rnodcoef(t,node);
  93.     if ~isempty(x)
  94.         if min(size(x))<2
  95.             plot(x,'Color','r','parent',axe_vis);
  96.             lx = length(x);
  97.             if lx> 1 , set(axe_vis,'Xlim',[1,lx]); end
  98.         else
  99.            NBC = 128;
  100.            colormap(pink(NBC))
  101.            image(wcodemat(x,NBC,'mat',0),'parent',axe_vis);
  102.         end
  103.         endTitle = '.';
  104.     else
  105.         delete(get(axe_vis,'Children'))
  106.         endTitle = ' ==> NONE.';
  107.     end
  108.     order = treeord(t);
  109.     [d,p] = ind2depo(order,node);
  110.     ldep = sprintf('(%0.f,%0.f)',d,p);
  111.     lind = sprintf('(%0.f)',node);
  112.     axeTitle = ['data for node: ' lind ' or ' ldep endTitle];
  113.     wtitle(axeTitle,'parent',axe_vis);
  114.     mousefrm(fig_tree,'arrow')
  115.   otherwise
  116.     try
  117.       nbout = nargout;
  118.       varargout{1:nbout} = plot(dtree,varargin{:});
  119.     catch
  120.       labels = tlabels(t,varargin{1});
  121.       varargout = {labels,0};
  122.     end
  123.         
  124. end