开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 图形图像 > 波变换 > 查看源码
cf1dtool.m
资源名称：wavelet_toolbox.rar　[点击查看]
上传用户：haiyisale
上传日期：2013-01-09
资源大小：3246k
文件大小：35k
源码类别：

波变换

开发平台：

Matlab

cf1dtool.m：源码内容
  1. function varargout = cf1dtool(option,varargin)
  2. %CF1DTOOL Wavelet Coefficients Selection 1-D tool.
  3. %   VARARGOUT = CF1DTOOL(OPTION,VARARGIN)
  5. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
  6. %   Last Revision: 28-Jun-2001.
  7. %   Copyright 1995-2002 The MathWorks, Inc.
  8. %   $Revision: 1.12 $  $Date: 2002/06/17 12:19:28 $
  10. % Test inputs.
  11. %-------------
  12. if nargin==0 , option = 'create'; end
  13. [option,winAttrb] = utguidiv('ini',option,varargin{:});
  15. % Default values.
  16. %----------------
  17. max_lev_anal = 9;
  19. % Stem parameters.
  20. %-----------------
  21. absMode = 0;
  22. appView = 1;
  24. % Memory Blocs of stored values.
  25. %===============================
  26. % MB0.
  27. %-----
  28. n_membloc0 = 'MB0';
  29. ind_sig    = 1;
  30. ind_coefs  = 2;
  31. ind_longs  = 3;
  32. ind_first  = 4;
  33. ind_last   = 5;
  34. ind_sort   = 6;
  35. ind_By_Lev = 7;
  36. ind_sizes  = 8;  % Dummy
  37. nb0_stored = 8;
  39. % MB1.
  40. %-----
  41. n_param_anal = 'MB1';
  42. ind_sig_name =  1;
  43. ind_sig_size =  2;
  44. ind_wav_name =  3;
  45. ind_lev_anal =  4;
  46. nb1_stored   =  4;
  48. % MB2.
  49. %-----
  50. n_InfoInit   = 'MB2';
  51. ind_filename =  1;
  52. ind_pathname =  2;
  53. nb2_stored   =  2;
  55. % MB3.
  56. %-----
  57. n_synt_sig = 'MB3';
  58. ind_ssig   =  1;
  59. nb3_stored =  1;
  61. % MB4.
  62. %-----
  63. n_miscella     = 'MB4';
  64. ind_graph_area =  1;
  65. ind_axe_hdl    =  2;
  66. ind_lin_hdl    =  3;
  67. nb4_stored     =  3;
  69. if ~isequal(option,'create') , win_tool = varargin{1}; end
  70. switch option
  71.   case {'create','close'} ,
  72.   otherwise
  73.     toolATTR = wfigmngr('getValue',win_tool,'ToolATTR');
  74.     hdl_UIC  = toolATTR.hdl_UIC;
  75.     hdl_MEN  = toolATTR.hdl_MEN;
  76.     pus_ana  = hdl_UIC.pus_ana;
  77.     chk_sho  = hdl_UIC.chk_sho;
  78. end
  79. switch option
  80.   case 'create'
  81.     % Get Globals.
  82.     %--------------
  83.     [Def_Btn_Height,Def_Btn_Width,Pop_Min_Width,  ...
  84.      Def_Txt_Height,X_Spacing,Y_Spacing,Def_FraBkColor ] = ...
  85.         mextglob('get',...
  86.             'Def_Btn_Height','Def_Btn_Width','Pop_Min_Width',  ...
  87.             'Def_Txt_Height','X_Spacing','Y_Spacing','Def_FraBkColor' ...
  88.             );
  90.     % Window initialization.
  91.     %----------------------
  92.     win_title = 'Wavelet Coefficients Selection 1-D';
  93.     [win_tool,pos_win,win_units,str_numwin,...
  94.         frame0,pos_frame0,Pos_Graphic_Area,pus_close] = ...
  95.            wfigmngr('create',win_title,winAttrb,'ExtFig_Tool_3',mfilename,1,1,0);
  96.     if nargout> 0 , varargout{1} = win_tool; end
  98. % Add Help for Tool.
  99. %------------------
  100. wfighelp('addHelpTool',win_tool,'One-Dimensional &Selection','CF1D_GUI');
  102.     % Menu construction.
  103.     %-------------------
  104.     m_files = wfigmngr('getmenus',win_tool,'file');
  105.     m_load = uimenu(m_files,...
  106.                     'Label','&Load Signal', ...
  107.                     'Position',1,           ...
  108.                     'Callback',             ...
  109.                     [mfilename '(''load'',' str_numwin ');'] ...
  110.                     );
  111.     m_save = uimenu(m_files,...
  112.                     'Label','&Save Synthesized Signal', ...
  113.                     'Position',2,     ...
  114.                     'Enable','Off',   ...
  115.                     'Callback',       ...
  116.                     [mfilename '(''save'',' str_numwin ');'] ...
  117.                     );
  118.     m_demo = uimenu(m_files,...
  119.                     'Label','&Example Analysis ','Position',3);
  120.     m_demo_1 = uimenu(m_demo,'Label','Basic Signals');
  121.     m_demo_2 = uimenu(m_demo,'Label','Noisy Signals');
  122.     m_demo_3 = uimenu(m_demo,'Label','Noisy Signals - Movie');
  124.     % Submenu of test signals.
  125.     %-------------------------
  126.     names(1,:)  = 'Sum of sines               ';
  127.     names(2,:)  = 'Frequency breakdown        ';
  128.     names(3,:)  = 'Uniform white noise        ';
  129.     names(4,:)  = 'AR(3) noise                ';
  130.     names(5,:)  = 'Noisy polynomial           ';
  131.     names(6,:)  = 'Noisy polynomial           ';
  132.     names(7,:)  = 'Step signal                ';
  133.     names(8,:)  = 'Two nearby discontinuities ';
  134.     names(9,:)  = 'Two nearby discontinuities ';
  135.     names(10,:) = 'Second derivative breakdown';
  136.     names(11,:) = 'Second derivative breakdown';
  137.     names(12,:) = 'Ramp + white noise         ';
  138.     names(13,:) = 'Ramp + colored noise       ';
  139.     names(14,:) = 'Sine + white noise         ';
  140.     names(15,:) = 'Triangle + sine            ';
  141.     names(16,:) = 'Triangle + sine + noise    ';
  142.     names(17,:) = 'Electrical consumption     ';
  143.     names(18,:) = 'Cantor curve               ';
  144.     names(19,:) = 'Koch curve                 ';
  146.     files = [ 'sumsin  ' ; 'freqbrk ' ; 'whitnois' ; 'warma   ' ; ...
  147.               'noispol ' ; 'noispol ' ; 'wstep   ' ; 'nearbrk ' ; ...
  148.               'nearbrk ' ; 'scddvbrk' ; 'scddvbrk' ; 'wnoislop' ; ...
  149.               'cnoislop' ; 'noissin ' ; 'trsin   ' ; 'wntrsin ' ; ...
  150.               'leleccum' ; 'wcantor ' ; 'vonkoch '                    ];
  152.     waves = ['db3' ; 'db5' ; 'db3' ; 'db3' ; 'db2' ; 'db3' ; 'db2' ; ...
  153.              'db2' ; 'db7' ; 'db1' ; 'db4' ; 'db3' ; 'db3' ; 'db5' ; ...
  154.              'db5' ; 'db5' ; 'db3' ; 'db1' ; 'db1'                      ];
  156.     levels = ['5';'5';'5';'5';'4';'4';'5';'5';'5';'2';'2';'6';'6';'5';...
  157.                     '6';'7';'5';'5';'5'];
  159.     beg_call_str = [mfilename '(''demo'',' str_numwin ','''];
  160.     for i=1:size(files,1)
  161.             libel = ['with ' waves(i,:) ' at level ' levels(i,:) ...
  162.                             '  --->  ' names(i,:)];
  163.             action = [beg_call_str files(i,:) ''',''' ...
  164.                             waves(i,:) ''',' levels(i,:) ');'];
  165.             uimenu(m_demo_1,'Label',libel,'Callback',action);
  166.     end
  168.     names = strvcat(...
  169.               'Noisy blocks','Noisy bumps','Noisy heavysin',     ...
  170.               'Noisy Doppler','Noisy quadchirp','Noisy mishmash' ...
  171.               );
  172.     files      = [ 'noisbloc' ; 'noisbump' ; 'heavysin' ; ...
  173.                    'noisdopp' ; 'noischir' ; 'noismima'      ];
  174.     waves  = ['sym8';'sym4';'sym8';'sym4';'db1 ';'db3 '];
  175.     levels = ['5';'5';'5';'5';'5';'5'];
  176.     beg_call_str = [mfilename '(''demo'',' str_numwin ','''];
  177.     for i=1:size(files,1)
  178.         libel = ['with ' waves(i,:) ' at level ' levels(i,:) ...
  179.                         '  --->  ' names(i,:)];
  180.         action = [beg_call_str files(i,:) ''',''' ...
  181.                         waves(i,:) ''',' levels(i,:) ');'];
  182.         uimenu(m_demo_2,'Label',libel,'Callback',action);
  183.     end
  185.     for i=1:size(files,1)
  186.         libel = ['with ' waves(i,:) ' at level ' levels(i,:) ...
  187.                         '  --->  ' names(i,:)];
  188.         action = [beg_call_str files(i,:) ''',''' ...
  189.                                waves(i,:) ''',' levels(i,:) ',' ...
  190.                                '{''Stepwise''}' ');'];
  191.         uimenu(m_demo_3,'Label',libel,'Callback',action);
  192.     end
  194.     % Begin waiting.
  195.     %---------------
  196.     wwaiting('msg',win_tool,'Wait ... initialization');
  198.     % General parameters initialization.
  199.     %-----------------------------------
  200.     dx = X_Spacing;   dx2 = 2*dx;
  201.     dy = Y_Spacing;   dy2 = 2*dy;
  202.     d_txt = (Def_Btn_Height-Def_Txt_Height);
  203.  
  204.     % Command part of the window.
  205.     %============================
  206.     % Data, Wavelet and Level parameters.
  207.     %------------------------------------
  208.     xlocINI = pos_frame0([1 3]);
  209.     ytopINI = pos_win(4)-dy;
  210.     toolPos = utanapar('create',win_tool, ...
  211.                   'xloc',xlocINI,'top',ytopINI,...
  212.                   'enable','off', ...
  213.                   'wtype','dwt'   ...
  214.                   );
  215.  
  216.     w_uic   = 1.5*Def_Btn_Width;
  217.     h_uic   = 1.5*Def_Btn_Height;
  218.     bdx     = (pos_frame0(3)-w_uic)/2;
  219.     x_left  = pos_frame0(1)+bdx;
  220.     y_low   = toolPos(2)-1.5*Def_Btn_Height-2*dy;
  221.     pos_ana = [x_left, y_low, w_uic, h_uic];
  223.     commonProp = {...
  224.         'Parent',win_tool, ...
  225.         'Unit',win_units,  ...
  226.         'Enable','off'     ...
  227.         };
  229.     str_ana = xlate('Analyze');
  230.     cba_ana = [mfilename '(''anal'',' str_numwin ');'];
  231.     pus_ana = uicontrol(commonProp{:},...
  232.                          'Style','Pushbutton', ...
  233.                          'Position',pos_ana,   ...
  234.                          'String',str_ana,     ...
  235.                          'Callback',cba_ana,   ...
  236.                          'Interruptible','On'  ...
  237.                          );
  239.     % Create coefficients tool.
  240.     %--------------------------
  241.     ytopCFS = pos_ana(2)-4*dy;
  242.     toolPos = utnbcfs('create',win_tool,...
  243.                       'toolOPT','cf1d',  ...
  244.                       'xloc',xlocINI,'top',ytopCFS);
  246.     % Create show checkbox.
  247.     %----------------------
  248.     w_uic = (3*pos_frame0(3))/4;
  249.     x_uic = pos_frame0(1)+(pos_frame0(3)-w_uic)/2;
  250.     h_uic = Def_Btn_Height;
  251.     y_uic = toolPos(2)-Def_Btn_Height/2-h_uic;
  252.     pos_chk_sho = [x_uic, y_uic, w_uic, h_uic];
  253.     str_chk_sho = 'Show Original Signal';
  254.     chk_sho = uicontrol(commonProp{:},...
  255.                         'Style','checkbox',     ...
  256.                         'Visible','on',         ...
  257.                         'Position',pos_chk_sho, ...
  258.                         'String',str_chk_sho    ...
  259.                         );
  261.     %  Normalisation.
  262.     %----------------
  263.     Pos_Graphic_Area = wfigmngr('normalize',win_tool,Pos_Graphic_Area);
  264.  
  265.     % Axes contruction.
  266.     %------------------
  267.     ax     = zeros(4,1);
  268.     pos_ax = zeros(4,4);
  269.     bdx = 0.05;
  270.      ecy_top = 0.04;
  271.     ecy_bot = 0.04;
  272.     ecy_mid = 0.06;
  273.     w_ax = (Pos_Graphic_Area(3)-3*bdx)/2;
  274.     h_ax = (Pos_Graphic_Area(4)-ecy_top-ecy_mid-ecy_bot)/3;
  275.     x_ax = bdx;
  276.     y_ax = Pos_Graphic_Area(2)+Pos_Graphic_Area(4)-ecy_top-h_ax;
  277.     pos_ax(1,:) = [x_ax y_ax w_ax h_ax];
  278.     x_ax = x_ax+w_ax+bdx;
  279.     pos_ax(4,:) = [x_ax y_ax w_ax h_ax];
  280.     x_ax = bdx;
  281.     y_ax = Pos_Graphic_Area(2)+ecy_bot;
  282.     pos_ax(2,:) = [x_ax y_ax w_ax 2*h_ax];
  283.     x_ax = x_ax+w_ax+bdx;
  284.     pos_ax(3,:) = [x_ax y_ax w_ax 2*h_ax];
  285.     for k = 1:4
  286.         ax(k) = axes(...
  287.                      'Parent',win_tool,      ...
  288.                      'Unit','normalized',    ...
  289.                      'Position',pos_ax(k,:), ...
  290.                      'Xtick',[],'Ytick',[],  ...
  291.                      'Box','on',             ...
  292.                      'Visible','off'         ...
  293.                      );
  294.     end
  296.     % Callbacks update.
  297.     %------------------
  298.     hdl_den = utnbcfs('handles',win_tool);
  299.     utanapar('set_cba_num',win_tool,[m_files;hdl_den(:)]);
  300.     pop_lev = utanapar('handles',win_tool,'lev');
  301.     tmp     = num2mstr([pop_lev chk_sho]);
  302.     end_cba = [str_numwin ',' tmp ');'];
  303.     cba_pop_lev = [mfilename '(''update_level'',' end_cba];
  304.     cba_chk_sho = [mfilename '(''show_ori_sig'',' str_numwin ');'];
  305.     set(pop_lev,'Callback',cba_pop_lev);
  306.     set(chk_sho,'Callback',cba_chk_sho);
  308.     % Memory for stored values.
  309.     %--------------------------
  310.     hdl_UIC  = struct('pus_ana',pus_ana,'chk_sho',chk_sho);
  311.     hdl_MEN  = struct('m_load',m_load,'m_save',m_save,'m_demo',m_demo);
  312.     toolATTR = struct('hdl_UIC',hdl_UIC,'hdl_MEN',hdl_MEN);
  313.     wfigmngr('storeValue',win_tool,'ToolATTR',toolATTR);
  314.     hdl_STEM = struct(...
  315.                       'Hstems_O',[], ...
  316.                       'H_vert_O',[], ...
  317.                       'H_stem_O',[], ...
  318.                       'H_vert_O_Copy',[], ...
  319.                       'H_stem_O_Copy',[], ...
  320.                       'Hstems_M',[], ...
  321.                       'H_vert_M',[], ...
  322.                       'H_stem_M',[], ...
  323.                       'H_vert_M_Copy',[], ...
  324.                       'H_stem_M_Copy',[]  ...
  325.                       );
  326.     wfigmngr('storeValue',win_tool,'Stems_struct',hdl_STEM);
  327.     wmemtool('ini',win_tool,n_InfoInit,nb0_stored);
  328.     wmemtool('ini',win_tool,n_param_anal,nb1_stored);
  329.     wmemtool('ini',win_tool,n_membloc0,nb2_stored);
  330.     wmemtool('ini',win_tool,n_synt_sig,nb3_stored);
  331.     wmemtool('ini',win_tool,n_miscella,nb4_stored);
  332.     wmemtool('wmb',win_tool,n_miscella,...
  333.                    ind_graph_area,Pos_Graphic_Area,ind_axe_hdl,ax);
  335.     % End waiting.
  336.     %---------------
  337.     wwaiting('off',win_tool);
  339.   case 'load'
  340.     % Loading file.
  341.     %-------------
  342.     if length(varargin)<2
  343.        [sigInfos,sig_Anal,ok] = ...
  344.             utguidiv('load_sig',win_tool,'Signal_Mask','Load Signal');
  345.         demoFlag = 0;
  346.     else
  347.         sig_Name = deblank(varargin{2});
  348.         wav_Name = deblank(varargin{3});
  349.         lev_Anal = varargin{4};
  350.         filename = [sig_Name '.mat'];
  351.         pathname = utguidiv('WTB_DemoPath',filename);
  352.         [sigInfos,sig_Anal,ok] = ...
  353.             utguidiv('load_dem1D',win_tool,pathname,filename);
  354.         demoFlag = 1;
  355.     end
  356.     if ~ok, return; end
  358.     % Begin waiting.
  359.     %---------------
  360.     wwaiting('msg',win_tool,'Wait ... loading');
  362.     % Get Values.
  363.     %------------
  364.     axe_hdl = wmemtool('rmb',win_tool,n_miscella,ind_axe_hdl);
  366.     % Cleaning.
  367.     %----------
  368.     dynvtool('stop',win_tool);
  369.     utnbcfs('clean',win_tool)
  370.     set(hdl_MEN.m_save,'Enable','Off');
  371.     set(axe_hdl(2:end),'Visible','Off');
  372.     children = allchild(axe_hdl);
  373.     delete(children{:});
  374.     set(axe_hdl,'Xtick',[],'Ytick',[],'Box','on');
  376.     % Setting GUI values.
  377.     %--------------------
  378.     sig_Name = sigInfos.name;
  379.     sig_Size = sigInfos.size;
  380.     sig_Size = max(sig_Size);
  381.     levm     = wmaxlev(sig_Size,'haar');
  382.     levmax   = min(levm,max_lev_anal);
  383.     lev      = min(levmax,5);
  384.     str_lev_data = int2str([1:levmax]');
  385.     if ~demoFlag
  386.         cbanapar('set',win_tool, ...
  387.                  'n_s',{sig_Name,sig_Size}, ...
  388.                  'lev',{'String',str_lev_data,'Value',lev});
  389.     else
  390.         cbanapar('set',win_tool, ...
  391.                  'n_s',{sig_Name,sig_Size}, ...
  392.                  'wav',wav_Name, ...
  393.                  'lev',{'String',str_lev_data,'Value',lev_Anal});
  394.         lev = lev_Anal;
  395.     end
  396.     set(chk_sho,'Value',0)
  397.     cf1dtool('position',win_tool,lev,chk_sho);
  399.     % Drawing.
  400.     %---------
  401.     axeAct = axe_hdl(1);
  402.     lsig   = length(sig_Anal);
  403.     wtitle('Original Signal','Parent',axeAct);
  404.     col_s = wtbutils('colors','sig');
  405.     lin_hdl(1) = line(...
  406.       'Parent',axeAct,  ...
  407.       'Xdata',[1:lsig], ...
  408.       'Ydata',sig_Anal, ...
  409.       'Color',col_s,'Visible','on'...
  410.       );
  411.     ymin = min(sig_Anal);
  412.     ymax = max(sig_Anal);
  413.     dy   = (ymax-ymin)/20;
  414.     set(axeAct,...
  415.         'Xlim',[1 lsig],'Ylim',[ymin-dy ymax+dy], ...
  416.         'XtickMode','auto','YtickMode','auto','Visible','on' ...
  417.         );
  418.     axeAct = axe_hdl(4);
  419.     wtitle('Synthesized Signal','Parent',axeAct);
  420.     lin_hdl(2) = line(...
  421.       'Parent',axeAct,  ...
  422.       'Xdata',[1:lsig], ...
  423.       'Ydata',sig_Anal, ...
  424.       'Color',col_s,'Visible','off'...
  425.       );
  426.     col_ss = wtbutils('colors','ssig');
  427.     lin_hdl(3) = line(...
  428.       'Parent',axeAct,  ...
  429.       'Xdata',[1:lsig], ...
  430.       'Ydata',sig_Anal, ...
  431.       'Color',col_ss,'Visible','off'...
  432.       );
  433.     set(axeAct,...
  434.         'Xlim',[1 lsig],'Ylim',[ymin-dy ymax+dy], ...
  435.         'XtickMode','auto','YtickMode','auto'     ...
  436.         );
  438.     % Setting Analysis parameters.
  439.     %-----------------------------
  440.     wmemtool('wmb',win_tool,n_membloc0,ind_sig,sig_Anal);
  441.     wmemtool('wmb',win_tool,n_param_anal, ...
  442.                    ind_sig_name,sigInfos.name,...
  443.                    ind_sig_size,sigInfos.size ...
  444.                    );
  445.     wmemtool('wmb',win_tool,n_InfoInit, ...
  446.                    ind_filename,sigInfos.filename, ...
  447.                    ind_pathname,sigInfos.pathname  ...
  448.                    );
  450.     % Store Values.
  451.     %--------------
  452.     wmemtool('wmb',win_tool,n_miscella,ind_lin_hdl,lin_hdl);
  454.     % Setting enabled values.
  455.     %------------------------
  456.     utnbcfs('set',win_tool,'handleORI',lin_hdl(1),'handleTHR',lin_hdl(3))
  457.     cbanapar('enable',win_tool,'on');
  458.     set(pus_ana,'Enable','On' );
  459.  
  460.     % End waiting.
  461.     %-------------
  462.     wwaiting('off',win_tool);
  464.   case 'demo'
  465.     cf1dtool('load',varargin{:})
  466.     if length(varargin)>4 
  467.         parDEMO = varargin{5};
  468.     else
  469.         parDEMO = {'Global'};
  470.     end
  472.     % Begin waiting.
  473.     %---------------
  474.     wwaiting('msg',win_tool,'Wait ... computing');
  476.     % Computing.
  477.     %-----------
  478.     cf1dtool('anal',win_tool);
  479.     pause(1)
  480.     utnbcfs('demo',win_tool,parDEMO);
  482.     % End waiting.
  483.     %-------------
  484.     wwaiting('off',win_tool);
  486.   case 'save'
  487.     % Testing file.
  488.     %--------------
  489.     [filename,pathname,ok] = utguidiv('test_save',win_tool, ...
  490.                                  '*.mat','Save Synthesized Signal');
  491.     if ~ok, return; end
  493.     % Begin waiting.
  494.     %--------------
  495.     wwaiting('msg',win_tool,'Wait ... saving');
  497.     % Getting Synthesized Signal.
  498.     %---------------------------
  499.     wname   = wmemtool('rmb',win_tool,n_param_anal,ind_wav_name);
  500.     lin_hdl = wmemtool('rmb',win_tool,n_miscella,ind_lin_hdl);
  501.     lin_hdl = lin_hdl(3);
  502.     sig     = get(lin_hdl,'Ydata');
  504.     % Saving file.
  505.     %--------------
  506.     [name,ext] = strtok(filename,'.');
  507.     if isempty(ext) | isequal(ext,'.')
  508.         ext = '.mat'; filename = [name ext];
  509.     end
  510.     eval([name ' = sig ;']);
  511.     saveStr = {name,'wname'};
  512.     wwaiting('off',win_tool);
  513.     try
  514.       save([pathname filename],saveStr{:});
  515.     catch
  516.       errargt(mfilename,'Save FAILED !','msg');
  517.     end
  519.   case 'anal'
  520.     % Waiting message.
  521.     %-----------------
  522.     wwaiting('msg',win_tool,'Wait ... computing');
  523.  
  524.     % Reading Analysis Parameters.
  525.     %-----------------------------
  526.     sig_Anal = wmemtool('rmb',win_tool,n_membloc0,ind_sig);
  527.     [wav_Name,lev_Anal] = cbanapar('get',win_tool,'wav','lev');
  529.     % Setting Analysis parameters
  530.     %-----------------------------
  531.     wmemtool('wmb',win_tool,n_param_anal, ...
  532.                    ind_wav_name,wav_Name, ...
  533.                    ind_lev_anal,lev_Anal ...
  534.                    );
  535.     % Get Values.
  536.     %------------
  537.     [axe_hdl,lin_hdl] = wmemtool('rmb',win_tool,n_miscella,...
  538.                                   ind_axe_hdl,ind_lin_hdl);
  540.     % Analyzing.
  541.     %-----------
  542.     [coefs,longs] = wavedec(sig_Anal,lev_Anal,wav_Name);
  543.     [tmp,idxsort] = sort(abs(coefs));
  544.     last  = cumsum(longs(1:end-1));
  545.     first = ones(size(last));
  546.     first(2:end) = last(1:end-1)+1;
  547.     len = length(last);
  548.     idxByLev = cell(1,len);
  549.     for k=1:len
  550.         idxByLev{k} = find((first(k)<=idxsort) & (idxsort<=last(k)));
  551.     end
  553.     % Writing coefficients.
  554.     %----------------------
  555.     wmemtool('wmb',win_tool,n_membloc0,...
  556.              ind_coefs,coefs,ind_longs,longs, ...
  557.              ind_first,first,ind_last,last, ...
  558.              ind_sort,idxsort,ind_By_Lev,idxByLev,...
  559.              ind_sizes,[]);
  560.  
  561.     % Clean axes and reset dynvtool.
  562.     %-------------------------------
  563.     hdls_all = get(axe_hdl(2:3),'Children');
  564.     delete(hdls_all{:});
  565.     set(axe_hdl(2:3),'YTickLabel',[],'YTick',[]);
  566.     dynvtool('ini_his',win_tool,'reset')
  568.     % Plot original decomposition.
  569.     %-----------------------------
  570.     xlim = [1,length(sig_Anal)];
  571.     set(axe_hdl(1:4),'Xlim',xlim);
  572.     ax_prop = {'Xlim',xlim,'box','on','XtickMode','auto','Visible','On'};
  573.     axeAct = axe_hdl(2);
  574.     axes(axeAct);
  575.     Hstems_O = dw1dstem(axeAct,coefs,longs,absMode,appView,'WTBX');
  576.     set(axeAct,ax_prop{:});
  577.     wtitle('Original Coefficients','Parent',axeAct);
  579.     % Plot modified decomposition.
  580.     %-----------------------------
  581.     axeAct = axe_hdl(3);
  582.     axes(axeAct);
  583.     Hstems_M = dw1dstem(axeAct,coefs,longs,absMode,appView,'WTBX');
  584.     set(axeAct,ax_prop{:});
  585.     wtitle('Selected Coefficients','Parent',axeAct);
  587.     % Plot signal and synthezised signal.
  588.     %------------------------------------
  589.     axeAct = axe_hdl(4);
  590.     set(axeAct,'Visible','on');
  591.     set(lin_hdl(3),'Ydata',sig_Anal,'Visible','on');
  593.     % Reset tool coefficients.
  594.     %-------------------------
  595.     utnbcfs('update_NbCfs',win_tool,'anal');
  596.     utnbcfs('update_methode',win_tool,'anal');
  597.     utnbcfs('enable',win_tool,'anal');
  598.     set(hdl_MEN.m_save,'Enable','On');
  600.     % Construction of the invisible Stems.
  601.     %-------------------------------------
  602.     cf1dtool('set_Stems_HDL',win_tool,'anal',Hstems_O,Hstems_M);
  604.     % Connect dynvtool.
  605.     %------------------
  606.     params = [axe_hdl(2:3)' , -lev_Anal];
  607.     dynvtool('init',win_tool,[],axe_hdl,[],[1 0], ...
  608.             '','','cf1dcoor',params,'cf1dselc',params);
  610.     % End waiting.
  611.     %-------------
  612.     wwaiting('off',win_tool);
  613.         
  614.   case 'apply'
  615.     % Waiting message.
  616.     %-----------------
  617.     wwaiting('msg',win_tool,'Wait ... computing');
  618.  
  619.     % Analysis Parameters.
  620.     %--------------------
  621.     [first,idxsort,idxByLev] = ...
  622.         wmemtool('rmb',win_tool,n_membloc0,ind_first,ind_sort,ind_By_Lev);
  623.     [nameMeth,nbkept] = utnbcfs('get',win_tool,'nameMeth','nbkept');
  624.     len = length(idxByLev);
  626.     switch nameMeth
  627.       case {'Global','ByLevel'}
  628.         ind = [];
  629.         for k=1:len
  630.             ind = [ind , idxByLev{k}(end-nbkept(k)+1:end)];
  631.         end
  632.         idx_Cfs = idxsort(ind);
  634.         % Computing & Drawing.
  635.         %---------------------
  636.         Hstems_M = cf1dtool('plot_NewDec',win_tool,idx_Cfs,nameMeth);
  638.         % Construction of the invisible Stems.
  639.         %-------------------------------------
  640.         cf1dtool('set_Stems_HDL',win_tool,'apply',Hstems_M);
  642.       case {'Manual'}
  643.         [H_stem_O,H_stem_O_Copy] = ...
  644.             cf1dtool('get_Stems_HDL',win_tool,'Manual');
  645.         ind = [];
  646.         for k=1:len
  647.             y = len+1-k;
  648.             x_stem = get(H_stem_O(y),'Xdata');
  649.             x_stem_Copy = get(H_stem_O_Copy(y),'Xdata');
  650.             Idx = find(ismember(x_stem,x_stem_Copy));
  651.             ind = [ind , Idx+first(k)-1];
  652.         end
  653.         idx_Cfs = ind;
  655.         % Computing & Drawing.
  656.         %---------------------
  657.         cf1dtool('plot_NewDec',win_tool,idx_Cfs,nameMeth);
  658.     end 
  660.     % End waiting.
  661.     %-------------
  662.     wwaiting('off',win_tool);
  664.   case 'Apply_Movie'
  665.     movieSET = varargin{2};
  666.     if isempty(movieSET)
  667.         Hstems_M = cf1dtool('plot_NewDec',win_tool,[],'Stepwise');
  668.         return
  669.     end
  670.     nbInSet = length(movieSET);  
  671.     appFlag = varargin{3};
  672.     popStop = varargin{4};
  674.     % Waiting message.
  675.     %-----------------
  676.     if nbInSet>1
  677.         txt_msg = wwaiting('msg',win_tool,'Wait ... computing');
  678.     end
  680.     % Get Analysis Parameters.
  681.     %-------------------------
  682.     [first,last,idxsort,idxByLev] = ...
  683.         wmemtool('rmb',win_tool,n_membloc0, ...
  684.                        ind_first,ind_last,ind_sort,ind_By_Lev);
  686.     % Computing.
  687.     %-----------
  688.     len = length(last);
  689.     nbKept = zeros(1,len+1);
  690.     switch appFlag
  691.       case 1
  692.         idx_App = idxsort(idxByLev{1});
  693.         App_Len = length(idx_App);
  694.         idxsort(idxByLev{1}) = [];
  696.       case 2
  697.         idx_App = [];
  698.         App_Len = 0;
  699.        
  700.       case 3
  701.         idx_App = [];
  702.         App_Len = 0;       
  703.         idxsort(idxByLev{1}) = [];
  704.     end
  705.     for jj = 1:nbInSet
  706.         nbcfs = movieSET(jj);
  707.         nbcfs  = nbcfs-App_Len;
  708.         idx_Cfs = [idx_App , idxsort(end-nbcfs+1:end)];
  709.         if nbInSet>1 , 
  710.             for k=1:len
  711.               dummy  = find((first(k)<=idx_Cfs) & (idx_Cfs<=last(k)));
  712.               nbKept(k) = length(dummy);
  713.             end
  714.             nbKept(end) = sum(nbKept(1:end-1));
  715.             msg2 = [int2str(nbKept(end)) '  = [' int2str(nbKept(1:end-1)) ']'];
  716.             msg  = strvcat(' ', sprintf('Number of kept coefficients:  %s', msg2)); 
  717.             set(txt_msg,'String',msg);
  718.         end
  720.         % Computing & Drawing.
  721.         %---------------------
  722.         Hstems_M = cf1dtool('plot_NewDec',win_tool,idx_Cfs,'Stepwise');
  724.         if nbInSet>1 , 
  725.             % Test for stopping.
  726.             %-------------------
  727.             user = get(popStop,'Userdata');
  728.             if isequal(user,1)
  729.                set(popStop,'Userdata',[]);
  730.                break
  731.             end
  732.             pause(0.1);
  733.         end
  734.     end
  736.     % Construction of the invisible Stems.
  737.     %-------------------------------------
  738.     cf1dtool('set_Stems_HDL',win_tool,'apply',Hstems_M);
  740.     % End waiting.
  741.     %-------------
  742.     if nbInSet>1 , wwaiting('off',win_tool); end
  744.   case {'select','unselect'}
  745.      OK_Select = isequal(option,'select');
  747.      % Find Select Box.
  748.      %-----------------
  749.      [X,Y] = mngmbtn('getbox',win_tool);
  750.      xmin = ceil(min(X));
  751.      xmax = floor(max(X));
  752.      ymin = min(Y);
  753.      ymax = max(Y);
  755.      % Get stored Stems.
  756.      %------------------
  757.      [H_vert_O,H_stem_O,H_vert_O_Copy,H_stem_O_Copy,...
  758.       H_vert_M,H_stem_M,H_vert_M_Copy,H_stem_M_Copy] = ...
  759.           cf1dtool('get_Stems_HDL',win_tool,'allComponents');
  760.      nb_Stems = length(H_stem_O);
  762.      % Find points.
  763.      %-------------         
  764.      nbKept = utnbcfs('get',win_tool,'nbKept');
  765.      ylow = max(1,floor(ymin));
  766.      ytop = min(ceil(ymax),nb_Stems);
  767.      for y = ylow:ytop
  768.         xy_stem      = get(H_stem_O(y),{'Xdata','Ydata'});
  769.         xy_stem_Copy = get(H_stem_O_Copy(y),{'Xdata','Ydata'});
  770.         Idx = find(xmin<=xy_stem{1} & xy_stem{1}<=xmax & ...
  771.                    ymin<=xy_stem{2} & xy_stem{2}<=ymax);
  772.         if OK_Select
  773.             Idx = Idx(~ismember(xy_stem{1}(Idx),xy_stem_Copy{1}));
  774.         else
  775.             Idx = find(ismember(xy_stem_Copy{1},xy_stem{1}(Idx)));
  776.         end
  777.         if ~isempty(Idx)
  778.             xy_vert_Copy = get(H_vert_O_Copy(y),{'Xdata','Ydata'});
  779.             if OK_Select
  780.                 xy_stem_Copy{1} = [xy_stem_Copy{1} , xy_stem{1}(Idx)];
  781.                 xy_stem_Copy{2} = [xy_stem_Copy{2} , xy_stem{2}(Idx)];
  782.                 tmp = [xy_stem{1}(Idx); xy_stem{1}(Idx) ; xy_stem{1}(Idx)];
  783.                 xy_vert_Copy{1} = [xy_vert_Copy{1} , tmp(:)'];
  784.                 nbIdx = length(Idx);
  785.                 tmp = [y*ones(1,nbIdx); xy_stem{2}(Idx) ; NaN*ones(1,nbIdx)];
  786.                 xy_vert_Copy{2} = [xy_vert_Copy{2} , tmp(:)'];
  787.             else
  788.                 xy_stem_Copy{1}(Idx) = [];
  789.                 xy_stem_Copy{2}(Idx) = [];
  790.                 Idx = 3*Idx-4;
  791.                 Idx = [Idx,Idx+1,Idx+2];
  792.                 xy_vert_Copy{1}(Idx) = [];
  793.                 xy_vert_Copy{2}(Idx) = [];
  794.             end
  795.             set([H_stem_O_Copy(y),H_stem_M_Copy(y)],...
  796.                 'Xdata',xy_stem_Copy{1},...
  797.                 'Ydata',xy_stem_Copy{2} ...
  798.                 );
  799.             set([H_vert_O_Copy(y),H_vert_M_Copy(y)],...
  800.                 'Xdata',xy_vert_Copy{1},...
  801.                 'Ydata',xy_vert_Copy{2} ...
  802.                 );        
  803.             nbInd = length(xy_stem_Copy{1})-1;
  804.             nbKept(nb_Stems+1-y) = nbInd;
  805.         end               
  806.      end
  807.      nbKept(end) = sum(nbKept(1:end-1));
  808.      utnbcfs('set',win_tool,'nbKept',nbKept);
  810.   case 'plot_NewDec'
  811.     % Indices of preserved coefficients & Methode.
  812.     %---------------------------------------------
  813.     idx_Cfs  = varargin{2};
  814.     nameMeth = varargin{3};
  815.     
  816.     % Get Handles.
  817.     %-------------
  818.     [axe_hdl,lin_hdl] = wmemtool('rmb',win_tool,n_miscella,...
  819.                                   ind_axe_hdl,ind_lin_hdl);
  821.     % Get Analysis Parameters.
  822.     %-------------------------
  823.     wav_Name = wmemtool('rmb',win_tool,n_param_anal,ind_wav_name);
  824.     [coefs,longs] = wmemtool('rmb',win_tool,n_membloc0,ind_coefs,ind_longs);
  826.     % Compute synthezized signal.
  827.     %---------------------------
  828.     Cnew = zeros(size(coefs));
  829.     Cnew(idx_Cfs) = coefs(idx_Cfs);
  830.     SS  = waverec(Cnew,longs,wav_Name);
  831.     thr = min(abs(Cnew));
  833.     % Plot modified decomposition.
  834.     %-----------------------------
  835.     xlim    = get(axe_hdl(1),'Xlim');
  836.     ax_prop = {'Xlim',xlim,'box','on'};
  837.     axeAct = axe_hdl(3);
  838.     if ~isequal(nameMeth,'Manual')
  839.         axes(axeAct);
  840.         varargout{1} = dw1dstem(axeAct,Cnew,longs,absMode,appView,'WTBX');
  841.         set(axeAct,ax_prop{:});
  842.     else
  843.         varargout{1} = [];
  844.     end
  846.     % Plot synthezised signal.
  847.     %-------------------------
  848.     axeAct = axe_hdl(4);
  849.     axes(axeAct);
  850.     set(lin_hdl(3),'Ydata',SS);
  851.     ymin = min(SS);
  852.     ymax = max(SS);
  853.     dy   = (ymax-ymin)/20;
  854.     cf1dtool('show_ori_sig',win_tool)
  855.     % wtitle('Synthesized Signal','Parent',axeAct);
  856.     set(axeAct,...
  857.         'Xlim',xlim,  ...
  858.         'XtickMode','auto','YtickMode','auto', ...
  859.         'Box','on'  ...
  860.         );
  861.     set(hdl_MEN.m_save,'Enable','On');   
  863.   case 'get_Stems_HDL'
  864.     % Output parameter.
  865.     %------------------
  866.     mode = varargin{2};
  868.     % Get stored Stems.
  869.     %------------------
  870.     hdl_STEM  = wfigmngr('getValue',win_tool,'Stems_struct');
  871.     varargout = struct2cell(hdl_STEM);
  872.     if isequal(mode,'All') , return; end
  873.     
  874.     [...
  875.      Hstems_O,H_vert_O,H_stem_O,H_vert_O_Copy,H_stem_O_Copy, ...
  876.      Hstems_M,H_vert_M,H_stem_M,H_vert_M_Copy,H_stem_M_Copy  ...
  877.      ] = deal(varargout{:});
  878.      switch mode
  879.        case 'allComponents'
  880.          varargout = {... 
  881.            H_vert_O,H_stem_O,H_vert_O_Copy,H_stem_O_Copy, ...
  882.            H_vert_M,H_stem_M,H_vert_M_Copy,H_stem_M_Copy};
  884.        case 'Manual'
  885.          varargout = {H_stem_O,H_stem_O_Copy};
  886.      end
  888.   case 'set_Stems_HDL'
  889.     mode = varargin{2};
  890.     axe_hdl  = wmemtool('rmb',win_tool,n_miscella,ind_axe_hdl);
  891.     lev_Anal = wmemtool('rmb',win_tool,n_param_anal,ind_lev_anal);
  892.     nb_STEMS = lev_Anal+1;
  893.     hdl_STEM = wfigmngr('getValue',win_tool,'Stems_struct'); 
  894.     switch mode
  895.       case 'anal'
  896.         Hstems_O = varargin{3};
  897.         Hstems_M = varargin{4};
  899.         % Construction of the invisible duplicated coefficients axes.
  900.         %------------------------------------------------------------
  901.         [H_vert_O,H_stem_O,Hstems_O] = extractSTEMS(Hstems_O);
  902.         [H_vert_M,H_stem_M,Hstems_M] = extractSTEMS(Hstems_M);
  904.         % Store values.
  905.         %--------------
  906.         hdl_STEM.Hstems_O = Hstems_O;
  907.         hdl_STEM.H_vert_O = H_vert_O;
  908.         hdl_STEM.H_stem_O = H_stem_O;
  909.         hdl_STEM.Hstems_M = Hstems_M;
  910.         hdl_STEM.H_vert_M = H_vert_M;
  911.         hdl_STEM.H_stem_M = H_stem_M;
  913.       case 'apply'
  914.         Hstems_M = varargin{3};
  916.         % Modification of Stems.
  917.         %-----------------------
  918.         [H_vert_M,H_stem_M,Hstems_M] = extractSTEMS(Hstems_M);
  919.         hdl_STEM.Hstems_M = Hstems_M;
  920.         hdl_STEM.H_vert_M = H_vert_M;
  921.         hdl_STEM.H_stem_M = H_stem_M;
  923.       case 'reset'
  924.         nameMeth = varargin{3};
  925.         hdl_DEL = [hdl_STEM.H_vert_O_Copy,hdl_STEM.H_stem_O_Copy,...
  926.                    hdl_STEM.H_vert_M_Copy,hdl_STEM.H_stem_M_Copy];
  928.         hdl_DEL = hdl_DEL(ishandle(hdl_DEL));
  929.         delete(hdl_DEL)
  930.         if ~isempty(hdl_STEM.Hstems_M)
  931.             HDL_VIS = hdl_STEM.Hstems_M(:,[2:4]);
  932.             HDL_VIS = HDL_VIS(ishandle(HDL_VIS(:)));
  933.         else
  934.             HDL_VIS = [];
  935.         end
  936.         switch nameMeth
  937.           case {'Global','ByLevel','Stepwise'}
  938.             H_vert_O_Copy = [];  H_stem_O_Copy = [];
  939.             H_vert_M_Copy = [];  H_stem_M_Copy = [];
  940.             if ~isequal(nameMeth,'Stepwise')
  941.                 vis_VIS = 'On';
  942.             else
  943.                 vis_VIS = 'Off';
  944.             end
  946.           case {'Manual'}
  947.             [H_vert_O_Copy,H_stem_O_Copy] = initSTEMS(axe_hdl(2),nb_STEMS);
  948.             [H_vert_M_Copy,H_stem_M_Copy] = initSTEMS(axe_hdl(3),nb_STEMS);
  949.             HDL_Copy = [H_vert_O_Copy(:);H_stem_O_Copy(:);...
  950.                         H_vert_M_Copy(:);H_stem_M_Copy(:)];
  951.             vis_VIS = 'Off';
  952.             set(HDL_Copy,'Visible','On');
  953.         end
  954.         set(HDL_VIS,'Visible',vis_VIS);
  955.         hdl_STEM.H_vert_O_Copy = H_vert_O_Copy;
  956.         hdl_STEM.H_stem_O_Copy = H_stem_O_Copy;
  957.         hdl_STEM.H_vert_M_Copy = H_vert_M_Copy;
  958.         hdl_STEM.H_stem_M_Copy = H_stem_M_Copy;
  959.     end
  960.     wfigmngr('storeValue',win_tool,'Stems_struct',hdl_STEM);
  962.   case 'show_ori_sig'
  963.     lin_hdl = wmemtool('rmb',win_tool,n_miscella,ind_lin_hdl);
  964.     vis = getonoff(get(chk_sho,'Value'));
  965.     set(lin_hdl(2),'Visible',vis);
  966.     strTitle = 'Synthezised Signal';
  967.     if isequal(vis(1:2),'on')
  968.         strTitle = [strTitle ' and Original Signal'];
  969.     end
  970.     wtitle(strTitle,'Parent',get(lin_hdl(2),'Parent'))
  972.   case 'position'
  973.     lev_view = varargin{2};
  974.     chk_sho  = varargin{3};
  975.     set(chk_sho,'Visible','off');
  976.     pos_old  = utnbcfs('get',win_tool,'position');
  977.     utnbcfs('set',win_tool,'position',{1,lev_view})
  978.     pos_new  = utnbcfs('get',win_tool,'position');
  979.     ytrans   = pos_new(2)-pos_old(2);
  980.     pos_chk  = get(chk_sho,'Position');
  981.     pos_chk(2) = pos_chk(2)+ytrans;
  982.     set(chk_sho,'Position',pos_chk,'Visible','on');
  984.   case 'update_level'
  985.     pop_lev = varargin{2}(1);
  986.     chk_sho = varargin{2}(2);
  987.     levmax  = get(pop_lev,'value');
  989.     % Get Values.
  990.     %------------
  991.     [axe_hdl,lin_hdl] = wmemtool('rmb',win_tool,n_miscella,...
  992.                                   ind_axe_hdl,ind_lin_hdl);
  993.     % Clean axes.
  994.     %------------
  995.     hdls_all = get([axe_hdl(2:3)],'Children');
  996.     delete(hdls_all{:});
  997.     set(axe_hdl(2:3),'Visible','Off','YTickLabel',[],'YTick',[]);
  999.     % Hide synthezised signal.
  1000.     %-------------------------
  1001.     set(chk_sho,'Value',0);
  1002.     set(wfindobj(axe_hdl(4)),'Visible','Off');
  1003.     set(hdl_MEN.m_save,'Enable','Off');
  1005.     % Reset coefficients tool and dynvtool.
  1006.     %--------------------------------------
  1007.     utnbcfs('clean',win_tool)
  1008.     cf1dtool('position',win_tool,levmax,chk_sho);
  1009.     dynvtool('ini_his',win_tool,'reset')
  1011.   case 'handles'
  1013.   case 'close'
  1014.  
  1015.   otherwise
  1016.     errargt(mfilename,'Unknown Option','msg');
  1017.     error('*');
  1018. end
  1021. %=============================================================================%
  1022. % INTERNAL FUNCTIONS
  1023. %=============================================================================%
  1024. %-----------------------------------------------------------------------------%
  1025. function varargout = copyManyObj(axe,varargin)
  1027. nbin = length(varargin);
  1028. varargout = cell(1,nbin);
  1029. for k = 1:nbin , varargout{k} = copyobj(varargin{k},axe); end
  1030. %-----------------------------------------------------------------------------%
  1031. function varargout = extractSTEMS(HDL_Stems)
  1033. nbrow = 4;
  1034. nbcol = length(HDL_Stems)/nbrow;
  1035. HDL_Stems = reshape(HDL_Stems(:),nbrow,nbcol)';
  1036. H_vert = HDL_Stems(:,2);
  1037. H_stem = HDL_Stems(:,3);
  1038. varargout = {H_vert,H_stem,HDL_Stems};
  1039. %-----------------------------------------------------------------------------%
  1040. function varargout = initSTEMS(axe,nbSTEMS)
  1042. stemColor = 'y';
  1043. stemColor = wtbutils('colors','stem');
  1044. linProp   = {...
  1045.              'Visible','Off',             ...
  1046.              'Xdata',NaN,'Ydata',NaN,     ...
  1047.              'MarkerEdgeColor',stemColor, ...
  1048.              'MarkerFaceColor',stemColor  ...
  1049.              };
  1050. linPropVert = {linProp{:},'linestyle','-','Color',stemColor};
  1051. linPropStem = {linProp{:},'linestyle','none','Marker','o','MarkerSize',3};
  1053. linTmpVert = line(linPropVert{:},'Parent',axe);
  1054. linTmpStem = line(linPropStem{:},'Parent',axe);
  1055. dupli      = ones(nbSTEMS,1);
  1056. HDL_vert = copyobj(linTmpVert,axe(dupli));
  1057. HDL_stem = copyobj(linTmpStem,axe(dupli));
  1058. delete([linTmpVert,linTmpStem]);
  1059. varargout = {HDL_vert,HDL_stem};
  1060. %-----------------------------------------------------------------------------%
  1061. %=============================================================================%