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

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

波变换

开发平台：

Matlab

dw1dcomp.m：源码内容
  1. function varargout = dw1dcomp(option,varargin)
  2. %DW1DCOMP Discrete wavelet 1-D compression.
  3. %   VARARGOUT = DW1DCOMP(OPTION,VARARGIN)
  5. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
  6. %   Last Revision: 03-Feb-2003.
  7. %   Copyright 1995-2004 The MathWorks, Inc.
  8. %   $Revision: 1.22.4.2 $
  10. % Default Value(s).
  11. %------------------
  12. def_nbCodeOfColors = 128;
  14. % Memory Blocks of stored values.
  15. %================================
  16. % MB1.
  17. %-----
  18. n_param_anal   = 'DWAn1d_Par_Anal';
  19. ind_sig_name   = 1;
  20. ind_sig_size   = 2;
  21. ind_wav_name   = 3;
  22. ind_lev_anal   = 4;
  23. ind_axe_ref    = 5;
  24. ind_act_option = 6;
  25. ind_ssig_type  = 7;
  26. ind_thr_val    = 8;
  27. nb1_stored     = 8;
  29. % MB2.
  30. %-----
  31. n_coefs_longs = 'Coefs_and_Longs';
  32. ind_coefs     = 1;
  33. ind_longs     = 2;
  34. nb2_stored    = 2;
  36. % MB1 (local).
  37. %-------------
  38. n_misc_loc = ['MB1_' mfilename];
  39. ind_sav_menus  = 1;
  40. ind_status     = 2;
  41. ind_win_caller = 3;
  42. ind_axe_datas  = 4;
  43. ind_hdl_datas  = 5;
  44. ind_cfsMode    = 6;
  45. ind_lin_cfs    = 7;
  46. ind_pop_mod    = 8;
  47. nbLOC_1_stored = 8;
  49. % MB2 (local).
  50. %-------------
  51. n_thrDATA = 'thrDATA';
  52. ind_value = 1;
  53. nbLOC_2_stored = 1;
  55. % Tag property.
  56. %--------------
  57. tag_axetxt_perf = 'Txt_Perf';
  59. if ~isequal(option,'create') , win_compress = varargin{1}; end
  60. switch option
  61.     case 'create'
  62.         % Get Globals.
  63.         %-------------
  64.         [Def_Btn_Height,Y_Spacing,Def_FraBkColor] = ...
  65.             mextglob('get','Def_Btn_Height','Y_Spacing','Def_FraBkColor' );
  67.         % Calling figure.
  68.         %----------------
  69.         win_caller = varargin{1};
  70.         str_win_caller = sprintf('%.0f',win_caller);
  72.         % Window initialization.
  73.         %----------------------
  74.         win_name = 'Wavelet 1-D  --  Compression';
  75.         [win_compress,pos_win,win_units,str_win_compress,...
  76.                 frame0,pos_frame0,Pos_Graphic_Area,pus_close] = ...
  77.                     wfigmngr('create',win_name,'','ExtFig_CompDeno', ...
  78.                                 strvcat(mfilename,'cond'),1,1,0);
  79.         set(win_compress,'userdata',win_caller);
  80.         varargout{1} = win_compress;
  82. % Add Help for Tool.
  83. %------------------
  84. wfighelp('addHelpTool',win_compress,'Signal Compression','DW1D_COMP_GUI');
  86. % Add Help Item.
  87. %----------------
  88. wfighelp('addHelpItem',win_compress,'Compression Procedure','COMP_PROCEDURE');
  89. wfighelp('addHelpItem',win_compress,'Available Methods','COMP_DENO_METHODS');
  91.         % Menu construction for current figure.
  92.         %--------------------------------------
  93. m_save  = wfigmngr('getmenus',win_compress,'save');
  94.         sav_menus(1) = uimenu(m_save,...
  95.                                 'Label','Compressed &Signal ',...
  96.                                 'Position',1,                   ...
  97.                                 'Enable','Off',                 ...
  98.                                 'Callback',                     ...
  99.                                 [mfilename '(''save_synt'','    ...
  100.                                         str_win_compress ');']  ...
  101.                                 );
  102.         sav_menus(2) = uimenu(m_save,...
  103.                                 'Label','&Coefficients ',       ...
  104.                                 'Position',2,                   ...
  105.                                 'Enable','Off',                 ...
  106.                                 'Callback',                     ...
  107.                                 [mfilename '(''save_cfs'','     ...
  108.                                         str_win_compress ');']  ...
  109.                                 );
  110.         sav_menus(3) = uimenu(m_save,...
  111.                                 'Label','&Decomposition ',      ...
  112.                                 'Position',3,                   ...
  113.                                 'Enable','Off',                 ...
  114.                                 'Callback',                     ...
  115.                                 [mfilename '(''save_dec'','     ...
  116.                                         str_win_compress ');']  ...
  117.                                 );
  119.         % Begin waiting.
  120.         %---------------
  121.         wwaiting('msg',win_compress,'Wait ... initialization');
  123.         % Getting Analysis parameters.
  124.         %-----------------------------
  125.         [Sig_Name,Wav_Name,Lev_Anal,Sig_Size] = ...
  126.                 wmemtool('rmb',win_caller,n_param_anal, ...
  127.                                ind_sig_name,ind_wav_name,...
  128.                                ind_lev_anal,ind_sig_size);
  129.         Wav_Fam  = wavemngr('fam_num',Wav_Name);
  130.         isBior   = wavemngr('isbior',Wav_Fam);
  131.         Sig_Anal = dw1dfile('sig',win_caller);
  133.         % General parameters initialization.
  134.         %-----------------------------------
  135.         dy = Y_Spacing;
  136.         if Lev_Anal>8 , ySpace = 2*dy; else , ySpace = 4*dy; end
  137.         str_pop_mod = strvcat('Global thresholding','By Level thresholding');
  139.         % Command & Graphic parts (common & global thresholding).
  140.         %========================================================
  141.         comFigProp = {'Parent',win_compress,'Unit',win_units};
  143.         % Data, Wavelet and Level parameters.
  144.         %------------------------------------
  145.         xlocINI = pos_frame0([1 3]);
  146.         ytopINI = pos_win(4)-dy;
  147.         toolPos = utanapar('create_copy',win_compress, ...
  148.                     {'xloc',xlocINI,'top',ytopINI},...
  149.                     {'n_s',{Sig_Name,Sig_Size},'wav',Wav_Name,'lev',Lev_Anal}...
  150.                     );
  152.         % Popup for mode.
  153.         %----------------
  154.         w_uic = (3*pos_frame0(3))/4;
  155.         h_uic = Def_Btn_Height;
  156.         y_uic = toolPos(2)-ySpace-h_uic;
  157.         x_uic = pos_frame0(1)+(pos_frame0(3)-w_uic)/2;
  158.         pos_pop_mod = [x_uic, y_uic, w_uic, h_uic];            
  159.         pop_mod = uicontrol(comFigProp{:},...
  160.                             'Style','Popup',...
  161.                             'Position',pos_pop_mod,...
  162.                             'Userdata',1,...
  163.                             'String',str_pop_mod...
  164.                             );
  165.         cba_pop_mod = [mfilename '(''change_mode'',' ...
  166.                           str_win_compress ',' num2mstr(pop_mod) ');'];
  167.         set(pop_mod,'Callback',cba_pop_mod);
  169.         % Global Compression tool.
  170.         %-------------------------
  171.         ytopTHR = pos_pop_mod(2)-ySpace;
  172.         utthrgbl('create',win_compress,'toolOPT','dw1dcomp', ...
  173.                  'xloc',xlocINI,'top',ytopTHR, ...
  174.                  'isbior',isBior ...
  175.                  );
  177.         % Adding colormap GUI.
  178.         %---------------------
  179.         [viewType,hdl_cfs] = dw1dvdrv('get_imgcfs',win_caller);
  180.         if isequal(viewType,'image')
  181.             [pop_pal_caller,mapName,nbColors] = ...
  182.                 cbcolmap('get',win_caller,'pop_pal','mapName','nbColors');
  183.             utcolmap('create',win_compress, ...
  184.                      'xloc',xlocINI, ...
  185.                      'bkcolor',Def_FraBkColor, ...
  186.                      'briflag',0, ...
  187.                      'enable','on');
  188.             pop_pal_loc = cbcolmap('get',win_compress,'pop_pal');
  189.             set(pop_pal_loc,'Userdata',get(pop_pal_caller,'Userdata'));
  190.             cbcolmap('set',win_compress,'pal',{mapName,nbColors});
  191.         end
  193.         % General graphical parameters initialization.
  194.         %--------------------------------------------
  195.         bdx     = 0.08*pos_win(3);
  196.         bdy     = 0.06*pos_win(4);
  197.         ecy     = 0.03*pos_win(4);
  198.         y_graph = 2*Def_Btn_Height+dy;
  199.         h_graph = pos_frame0(4)-y_graph;
  200.         w_graph = pos_frame0(1);
  201.         fontsize = wmachdep('fontsize','normal',9,Lev_Anal);
  203.         % Axes construction parameters.
  204.         %------------------------------
  205.         w_left     = (w_graph-3*bdx)/2;
  206.         x_left     = bdx;
  207.         w_right    = w_left;
  208.         x_right    = x_left+w_left+5*bdx/4;
  209.         n_axeleft  = Lev_Anal;
  210.         n_axeright = 3;
  211.         ind_left   = n_axeleft;
  212.         ind_right  = n_axeright;
  214.         % Vertical separation.
  215.         %---------------------
  216.         w_fra = 0.01*pos_win(3);
  217.         x_fra = (w_graph-w_fra)/2;
  218.         uicontrol('Parent',win_compress,...
  219.                   'Style','frame',...
  220.                   'Unit',win_units,...
  221.                   'Position',[x_fra,y_graph,w_fra,h_graph],...
  222.                   'Backgroundcolor',Def_FraBkColor ...
  223.                   );
  225.         % Building axes on the right part.
  226.         %---------------------------------
  227.         ecy_right  = 2*ecy;
  228.         h_right    =(h_graph-2*bdy-(n_axeright-1)*ecy_right)/n_axeright;
  229.         y_right    = y_graph+2*bdy/3;
  230.         axe_datas  = zeros(1,n_axeright);
  231.         pos_right  = [x_right y_right w_right h_right];
  232.         for k = 1:n_axeright
  233.             axe_datas(k) = axes(...
  234.                                 'Parent',win_compress,  ...
  235.                                 'Units',win_units,...
  236.                                 'Position',pos_right,...
  237.                                 'Drawmode','fast',...
  238.                                 'Box','On' ...
  239.                                 );
  240.             pos_right(2) = pos_right(2)+pos_right(4)+ecy_right;
  241.         end
  242.         set(axe_datas(1),'visible','off');
  244.         % Displaying the signal.
  245.         %-----------------------
  246.         hdl_datas = [NaN ; NaN];
  247.         axeAct = axe_datas(3);
  248.         axes(axeAct)
  249.         curr_color   = wtbutils('colors','sig');
  250.         hdl_datas(1) = line(...
  251.                            'Xdata',1:length(Sig_Anal),...
  252.                            'Ydata',Sig_Anal,...
  253.                            'erasemode','none',...
  254.                            'Color',curr_color,...
  255.                            'Parent',axeAct);
  256.         wtitle('Original signal','Parent',axeAct);
  257.         xlim = [1              Sig_Size];
  258.         ylim = [min(Sig_Anal)  max(Sig_Anal)];
  259.         if xlim(1)==xlim(2) , xlim = xlim+[-0.01 0.01]; end
  260.         if ylim(1)==ylim(2) , ylim = ylim+[-0.01 0.01]; end
  261.         set(axeAct,'Xlim',xlim,'Ylim',ylim);
  263.         % Displaying original details coefficients.
  264.         %------------------------------------------
  265.         axeAct = axe_datas(2);
  266.         [details,set_ylim,ymin,ymax] = dw1dfile('cfs_beg',win_caller,...
  267.                             [1:Lev_Anal],1);
  268.         
  269.         if isequal(viewType,'image')
  270.             flagType = 1;
  271.             cfsMode  = [];
  272.             [nul,i_min] = min(abs(details(:)));
  273.             if ~isempty(hdl_cfs)
  274.                 col_cfs = flipud(get(hdl_cfs,'Cdata'));
  275.             else
  276.                 col_cfs = wcodemat(details,def_nbCodeOfColors,'row',1);
  277.             end
  278.             col_min = col_cfs(i_min);
  279.             col_cfs = flipud(col_cfs);
  280.             axes(axeAct);
  281.             hdl_cfs_ori = image(col_cfs,'Parent',axeAct,'Userdata',col_min);
  282.             clear col_cfs
  283.             levlab = int2str([Lev_Anal:-1:1]');
  284.         else
  285.             flagType = -1;
  286.             axes(axeAct);
  287.             hdl_stem = copyobj(hdl_cfs,axeAct);
  288.             set(hdl_stem,'Visible','on');
  289.             cfsMode  = get(hdl_stem(1),'Userdata');
  290.             levlab = int2str([1:Lev_Anal]');
  291.             hdl_cfs_ori = [];
  292.         end
  293.         set(axeAct,...
  294.               'Userdata',hdl_cfs_ori,...
  295.               'Xlim',[1 Sig_Size],...
  296.               'YTicklabelMode','manual',...
  297.               'YTick',[1:Lev_Anal],...
  298.               'YTickLabel',levlab,...
  299.               'Ylim',[0.5 Lev_Anal+0.5]...
  300.               );
  301.         wtitle('Original coefficients','Parent',axeAct);
  302.         wylabel('Level number','Parent',axeAct);
  303.         xylim = get(axeAct,{'Xlim','Ylim'});
  304.         set(axe_datas(1),'Xlim',xylim{1},'Ylim',xylim{2});
  306.         % Initializing global threshold.
  307.         %-------------------------------
  308.         [valTHR,maxTHR,thresVALUES,rl2SCR,n0SCR] = ...
  309.             dw1dcomp('compute_GBL_THR',win_compress,win_caller);
  310.         utthrgbl('set',win_compress,'thrBOUNDS',[0,valTHR,maxTHR]);
  312.        % Displaying perfos & legend.
  313.         %--------------------------
  314.         y_axe = y_graph+2*bdy/3+h_right+ecy_right;
  315.         h_axe = 2*h_right+ecy_right;
  316.         pos_axe_perfo = [x_left y_axe w_left h_axe];
  317.         y_axe = y_graph+h_right/2;
  318.         h_axe = h_right/2;
  319.         pos_axe_legend = [x_left y_axe w_left h_axe];
  320.         utthrgbl('displayPerf',win_compress, ...
  321.                   pos_axe_perfo,pos_axe_legend,thresVALUES,n0SCR,rl2SCR,valTHR);
  322.         [perfl2,perf0] = utthrgbl('getPerfo',win_compress);
  323.         utthrgbl('set',win_compress,'perfo',[perfl2,perf0]);
  324.         drawnow
  326.         % Command & Graphic parts (by Level thresholding).
  327.         %=================================================
  328.         utthrw1d('create',win_compress, ...
  329.                  'xloc',xlocINI,'top',ytopTHR,...
  330.                  'ydir',-1, ...
  331.                  'visible','off', ...
  332.                  'enable','on', ...
  333.                  'levmax',Lev_Anal, ...
  334.                  'levmaxMAX',Lev_Anal, ...
  335.                  'isbior',isBior,  ...
  336.                  'toolOPT','comp' ...
  337.                  );
  339.         % Building axes on the left part.
  340.         %--------------------------------
  341.         comAxeProp = {comFigProp{:}, ...
  342.                       'Visible','Off',...
  343.                       'Drawmode','normal',...
  344.                       'Box','On'...
  345.                      };
  346.         ecy_left  = ecy/2;
  347.         h_left    = (h_graph-2*bdy-(n_axeleft-1)*ecy_left)/n_axeleft;
  348.         y_left    = y_graph+bdy;
  349.         axe_left  = zeros(1,n_axeleft);
  350.         pos_left  = [x_left y_left w_left h_left];
  351.         for k = 1:n_axeleft
  352.             if k~=1
  353.                 axe_left(k) = axes(comAxeProp{:},'Position',pos_left,...
  354.                                    'XTicklabelMode','manual','XTickLabel',[]);
  355.             else
  356.                 axe_left(k) = axes(comAxeProp{:},'Position',pos_left);
  357.             end
  358.             pos_left(2) = pos_left(2)+pos_left(4)+ecy_left;
  359.         end
  360.         wtitle('Original details coefficients','Parent',axe_left(Lev_Anal));
  361.         utthrw1d('set',win_compress,'axes',axe_left);
  363.         % Initializing by level threshold.
  364.         %---------------------------------
  365.         maxTHR = zeros(1,Lev_Anal);
  366.         for k = 1:Lev_Anal , maxTHR(k) = max(abs(details(k,:))); end
  367.         [valTHR,perfl2,perf0] = ...
  368.             dw1dcomp('compute_LVL_THR',win_compress,win_caller);
  369.         valTHR = min(valTHR,maxTHR);
  371.         % Displaying details.
  372.         %-------------------
  373.         col_det  = wtbutils('colors','det',Lev_Anal);
  374.         txt_left = ones(Lev_Anal,1);
  375.         for k = Lev_Anal:-1:1
  376.             axeAct  = axe_left(ind_left);
  377.             axes(axeAct);
  378.             lin_cfs(k) = line(...
  379.                            'Parent',axeAct,...
  380.                            'Visible','off', ...
  381.                            'Xdata',1:Sig_Size,...
  382.                            'Ydata',details(k,:),...
  383.                            'Color',col_det(k,:));
  384.             txt_left(k) = txtinaxe('create',['d' wnsubstr(k)],...
  385.                                     axe_left(k),'left',...
  386.                                     'off','bold',fontsize);
  387.             utthrw1d('plot_dec',win_compress,k, ...
  388.                      {maxTHR(k),valTHR(k),1,Sig_Size,k})
  390.             maxi = max([abs(ymax(k)),abs(ymin(k))]);
  391.             if abs(maxi)<eps , maxi = maxi+0.01; end;
  392.             ylim = 1.1*[-maxi maxi];
  393.             set(axe_left(ind_left),'Xlim',xlim,'Ylim',ylim);
  394.             ind_left = ind_left-1;
  395.         end
  397.         % Initialization of Compression structure.
  398.         %----------------------------------------
  399.         xmin = 1; xmax = Sig_Size;
  400.         utthrw1d('set',win_compress,...
  401.                        'thrstruct',{xmin,xmax,valTHR,lin_cfs},...
  402.                        'intdepthr',[]);
  404. % Add Context Sensitive Help (CSHelp).
  405. %-------------------------------------
  406. wfighelp('add_ContextMenu',win_compress,pop_mod,'DW1D_COMP_GUI');
  407. %-------------------------------------
  409.         % Memory blocks update.
  410.         %----------------------
  411.         utthrgbl('set',win_compress,'handleORI',hdl_datas(1));
  412.         utthrw1d('set',win_compress,'handleORI',hdl_datas(1));
  413.         wmemtool('ini',win_compress,n_misc_loc,nbLOC_1_stored);
  414.         wmemtool('wmb',win_compress,n_misc_loc,  ...
  415.                        ind_sav_menus,sav_menus,  ...
  416.                        ind_status,0,             ...
  417.                        ind_win_caller,win_caller,...
  418.                        ind_axe_datas,axe_datas,  ...
  419.                        ind_hdl_datas,hdl_datas,  ...
  420.                        ind_cfsMode,cfsMode,      ...
  421.                        ind_lin_cfs,lin_cfs,      ...
  422.                        ind_pop_mod,pop_mod       ...
  423.                        );
  424.         wmemtool('ini',win_compress,n_thrDATA,nbLOC_2_stored);
  426.         % Axes attachment.
  427.         %-----------------
  428.         if flagType~=-1
  429.             axe_cmd = [axe_datas(1:3) axe_left(1:n_axeleft)];
  430.             axe_act = [];
  431.             axe_cfs = [axe_datas(1:2)];
  432.         else
  433.             axe_cmd = [axe_datas(1:3) axe_left(1:n_axeleft)];
  434.             axe_act = [];
  435.             axe_cfs = [axe_datas(1:2)];
  436.         end
  437.         dynvtool('init',win_compress,[],axe_cmd,axe_act,[1 0], ...
  438.                         '','','dw1dcoor',[win_caller,axe_cfs,flagType*Lev_Anal]);
  440.         % Setting units to normalized.
  441.         %-----------------------------
  442.         wfigmngr('normalize',win_compress);
  444.         % End waiting.
  445.         %-------------
  446.         wwaiting('off',win_compress);
  448.     case 'compress'
  450.         % Waiting message.
  451.         %-----------------
  452.         wwaiting('msg',win_compress,'Wait ... computing');
  454.         % Clear & Get Handles.
  455.         %----------------------
  456.         dw1dcomp('clear_GRAPHICS',win_compress);
  457.         [win_caller,pop_mod]  = wmemtool('rmb',win_compress,n_misc_loc, ...
  458.                                                ind_win_caller,ind_pop_mod);
  459.         [axe_datas,hdl_datas] = wmemtool('rmb',win_compress,n_misc_loc, ...
  460.                                                ind_axe_datas,ind_hdl_datas);
  462.         % Getting memory blocks.
  463.         %-----------------------
  464.         [Wav_Name,Lev_Anal,Sig_Size] = ...
  465.                         wmemtool('rmb',win_caller,n_param_anal,       ...
  466.                                        ind_wav_name,ind_lev_anal,ind_sig_size);
  467.         [coefs,longs] = wmemtool('rmb',win_caller,n_coefs_longs,...
  468.                                        ind_coefs,ind_longs);
  469.         isBior = wavemngr('isbior',Wav_Name);
  471.         % Loading the original signal.
  472.         %-----------------------------
  473.         Sig_Anal = dw1dfile('sig',win_caller);
  475.         mode_val = get(pop_mod,'value');
  476.         if isBior
  477.             topTitle = 'Energy ratio ';
  478.         else
  479.             topTitle = 'Retained energy ';
  480.         end
  481.         switch mode_val
  482.           case 1    % Global thresholding
  483.             valTHR = utthrgbl('get',win_compress,'valthr');
  484.             thrParams = {'gbl',coefs,longs,Wav_Name,Lev_Anal,valTHR,'h',1};
  485.             if isBior
  486.                 [xc,cxc,lxc,perf0] = wdencmp(thrParams{:});
  487.                 normsig = norm(Sig_Anal);
  488.                 perfl2  = 100;
  489.                 if normsig>eps , perfl2 = perfl2*(norm(xc)/normsig)^2; end
  490.             else
  491.                 [xc,cxc,lxc]   = wdencmp(thrParams{:});
  492.                 [perfl2,perf0] = utthrgbl('getPerfo',win_compress);
  493.             end
  495.           case 2    % By level thresholding
  496.             cxc = utthrw1d('den_M2',win_compress,coefs,longs);
  497.             lxc = longs;
  498.             xc  = waverec(cxc,longs,Wav_Name);            
  499.             normsig = norm(Sig_Anal);
  500.             perfl2  = 100;
  501.             if normsig>eps , perfl2 = perfl2*(norm(xc)/normsig)^2; end
  502.             perf0 = 100*(length(find(cxc==0))/length(cxc));
  503.         end
  505.         % Displaying compressed signal.
  506.         %------------------------------
  507.         hdl_comp = hdl_datas(2);
  508.         if ishandle(hdl_comp)
  509.             set(hdl_comp,'Ydata',xc,'Visible','on');
  510.         else
  511.             curr_color = wtbutils('colors','ssig');
  512.             hdl_comp = line('Xdata',1:length(xc),...
  513.                             'Ydata',xc,...
  514.                             'color',curr_color,...
  515.                             'Parent',axe_datas(3));
  516.             hdl_datas(2) = hdl_comp;
  517.             wmemtool('wmb',win_compress,n_misc_loc,...
  518.                            ind_hdl_datas,hdl_datas);
  519.             utthrgbl('set',win_compress,'handleTHR',hdl_comp);
  520.             utthrw1d('set',win_compress,'handleTHR',hdl_comp);
  521.         end     
  523.         % Set a text as a super title.
  524.         %-----------------------------
  525.         wtitle('Original and compressed signals','Parent',axe_datas(3));
  526.         strPerfo = [topTitle num2str(perfl2,'%5.2f') ...
  527.                     ' % -- Zeros ' num2str(perf0,'%5.2f') ' %'];
  528.         wtxttitl(axe_datas(3),strPerfo,tag_axetxt_perf);
  529.         if mode_val==2
  530.             utthrw1d('set',win_compress,'perfos',{perfl2,perf0});
  531.         end 
  533.         % Displaying thresholded details coefficients.
  534.         %---------------------------------------------
  535.         cfsMode = wmemtool('rmb',win_compress,n_misc_loc,ind_cfsMode);
  536.         if isempty(cfsMode)
  537.             col_cfs   = wrepcoef(cxc,lxc);
  538.             nz_cfs    = find(col_cfs~=0);
  539.             [nbr,nbc] = size(col_cfs);
  540.             img_cfs   = get(axe_datas(2),'Userdata');
  541.             col_min   = get(img_cfs,'Userdata');
  542.             cfs_ori   = flipud(get(img_cfs,'Cdata'));
  543.             col_cfs   = col_min*ones(nbr,nbc);
  544.             col_cfs(nz_cfs) = cfs_ori(nz_cfs);
  545.             image(flipud(col_cfs),'Parent',axe_datas(1));
  546.             levlab = int2str([Lev_Anal:-1:1]');
  547.         else
  548.             dw1dstem(axe_datas(1),cxc,lxc,'mode',cfsMode,'colors','WTBX');
  549.             levlab = int2str([1:Lev_Anal]');
  550.         end
  552.         set(axe_datas(1),...
  553.             'clipping','on',            ...
  554.             'Xlim',get(axe_datas(2),'Xlim'),...
  555.             'YTicklabelMode','manual',  ...
  556.             'YTick',[1:Lev_Anal],       ...
  557.             'YTickLabel',levlab,        ...
  558.             'Ylim',[0.5 Lev_Anal+0.5]   ...
  559.             );
  561.         wtitle('Thresholded coefficients','Parent',axe_datas(1));
  562.         wylabel('Level number','Parent',axe_datas(1));
  563.         set(findobj(axe_datas(1)),'Visible','on');
  565.         % Memory blocks & HG update.
  566.         %---------------------------
  567.         switch mode_val
  568.           case 1
  569.             thrParams = utthrgbl('get',win_compress,'valthr');
  570.           case 2
  571.             thrStruct = utthrw1d('get',win_compress,'thrstruct');
  572.             thrParams = {thrStruct(1:Lev_Anal).thrParams};
  573.         end
  574.         wmemtool('wmb',win_compress,n_thrDATA,ind_value,{xc,cxc,lxc,thrParams});
  575.         dw1dcomp('enable_menus',win_compress,'on');
  577.         % End waiting.
  578.         %-------------
  579.         wwaiting('off',win_compress);
  581.     case 'change_mode'
  582.         pop_mod = varargin{2}(1);
  583.         mod_val = get(pop_mod,'value');
  584.         old_mod = get(pop_mod,'userdata');
  585.         if isequal(mod_val,old_mod) , return; end
  586.         set(pop_mod,'userdata',mod_val);
  587.         dw1dcomp('clear_GRAPHICS',win_compress);
  588.         switch mod_val
  589.           case 1 , visGBL = 'on';  visLVL = 'off';
  590.           case 2 , visGBL = 'off'; visLVL = 'on';
  591.         end
  592.         win_caller = wmemtool('rmb',win_compress,n_misc_loc,ind_win_caller);
  593.         viewType = dw1dvdrv('get_imgcfs',win_caller);
  594.         if isequal(viewType,'image')
  595.             visMAP = 'on';
  596.             if mod_val==2
  597.                Lev_Anal = wmemtool('rmb',win_caller,n_param_anal,ind_lev_anal);
  598.                if Lev_Anal>5 , visMAP = 'off'; else , visMAP = 'on'; end
  599.             end
  600.             cbcolmap('visible',win_compress,visMAP);
  601.         end            
  602.         utthrgbl('visible',win_compress,visGBL);
  603.         utthrw1d('visible',win_compress,visLVL);
  605.     case 'compute_GBL_THR'
  606.         win_caller = varargin{2};
  607.         [numMeth,meth,sliVal] = utthrgbl('get_GBL_par',win_compress);
  608.         [coefs,longs] = wmemtool('rmb',win_caller,n_coefs_longs,...
  609.                                        ind_coefs,ind_longs);
  610.         thrFLAGS = 'dw1dcompGBL';
  611.         switch numMeth
  612.           case 1
  613.             [valTHR,maxTHR,thresVALUES,rl2SCR,n0SCR] = ...
  614.                 wthrmngr(thrFLAGS,meth,coefs,longs);
  615.             if nargout==1 , varargout = {valTHR};
  616.             else          , varargout = {valTHR,maxTHR,thresVALUES,rl2SCR,n0SCR};
  617.             end
  619.           case 2
  620.             sig = dw1dfile('sig',win_caller);
  621.             valTHR = wthrmngr(thrFLAGS,meth,sig);
  622.             maxTHR = max(coefs);
  623.             valTHR = min(valTHR,maxTHR);
  624.             varargout = {valTHR};
  625.         end
  627.     case 'update_GBL_meth'
  628.         dw1dcomp('clear_GRAPHICS',win_compress);
  629.         win_caller = wmemtool('rmb',win_compress,n_misc_loc,ind_win_caller);
  630.         valTHR = dw1dcomp('compute_GBL_THR',win_compress,win_caller);
  631.         utthrgbl('update_GBL_meth',win_compress,valTHR);
  633.     case 'show_LVL_perfos'
  634.         win_caller = wmemtool('rmb',win_compress,n_misc_loc,ind_win_caller);
  635.         [coefs,longs] = wmemtool('rmb',win_caller,n_coefs_longs,...
  636.                                        ind_coefs,ind_longs);
  637.         lev_anal = wmemtool('rmb',win_caller,n_param_anal,ind_lev_anal);
  638.         [numMeth,meth,scal,sorh] = utthrw1d('get_LVL_par',win_compress);
  639.         valTHR = utthrw1d('get',win_compress,'valTHR');
  640.         [perfl2,perf0] = wscrupd(coefs,longs,lev_anal,valTHR,sorh);      
  641.         utthrw1d('set',win_compress,'perfos',{perfl2,perf0}); 
  643.     case 'compute_LVL_THR'
  644.         win_caller = varargin{2};
  645.         [numMeth,meth,alfa,sorh] = utthrw1d('get_LVL_par',win_compress);
  646.        [coefs,longs] = wmemtool('rmb',win_caller,n_coefs_longs,...
  647.                                        ind_coefs,ind_longs);
  648.         level = wmemtool('rmb',win_caller,n_param_anal,ind_lev_anal);
  650.         thrFLAGS = 'dw1dcompLVL';
  651.         switch numMeth
  652.           case {1,2,3,4} , valTHR = wthrmngr(thrFLAGS,meth,coefs,longs,alfa);
  653.           case 5          
  654.             sig = dw1dfile('sig',win_caller);
  655.             valTHR = wthrmngr(thrFLAGS,meth,sig,level);
  656.         end
  657.         [perfl2,perf0] = wscrupd(coefs,longs,level,valTHR,sorh);
  658.         varargout = {valTHR,perfl2,perf0};
  659.         utthrw1d('set',win_compress,'perfos',{perfl2,perf0}); 
  661.     case 'update_LVL_meth'
  662.         dw1dcomp('clear_GRAPHICS',win_compress);
  663.         win_caller = wmemtool('rmb',win_compress,n_misc_loc,ind_win_caller);
  664.         [valTHR,perfl2,perf0] = ...
  665.             dw1dcomp('compute_LVL_THR',win_compress,win_caller);
  666.         utthrw1d('update_LVL_meth',win_compress,valTHR);
  668.     case 'clear_GRAPHICS'
  669.         status = wmemtool('rmb',win_compress,n_misc_loc,ind_status);
  670.         if status == 0 , return; end
  672.         % Diseable Toggles and Menus.
  673.         %----------------------------
  674.         dw1dcomp('enable_menus',win_compress,'off');
  676.         % Get Handles.
  677.         %-------------
  678.         [axe_datas,hdl_datas] = wmemtool('rmb',win_compress,n_misc_loc, ...
  679.                                                ind_axe_datas,ind_hdl_datas);
  681.         % Setting the compressed coefs axes invisible.
  682.         %---------------------------------------------
  683.         hdl_comp = hdl_datas(2);
  684.         if ishandle(hdl_comp)
  685.             vis = get(hdl_comp,'Visible');
  686.             if vis(1:2)=='on'
  687.                 txt_perf = findobj(axe_datas(3),'tag',tag_axetxt_perf);
  688.                 if ishandle(txt_perf) , set(txt_perf,'Visible','off'); end
  689.                 set(findobj(axe_datas(1)),'Visible','off');
  690.                 wtitle('Original signal','Parent',axe_datas(3));
  691.                 set(hdl_comp,'Visible','off');
  692.             end
  693.         end
  694.         
  695.     case 'enable_menus'
  696.         enaVal = varargin{2};
  697.         sav_menus = wmemtool('rmb',win_compress,n_misc_loc,ind_sav_menus);
  698.         set(sav_menus,'Enable',enaVal);
  699.         utthrgbl('enable_tog_res',win_compress,enaVal);
  700.         utthrw1d('enable_tog_res',win_compress,enaVal);
  701.         if strncmpi(enaVal,'on',2) , status = 1; else , status = 0; end
  702.         wmemtool('wmb',win_compress,n_misc_loc,ind_status,status);
  704.     case 'save_synt'
  705.         win_compress = varargin{1};
  707.         % Testing file.
  708.         %--------------
  709.         [filename,pathname,ok] = utguidiv('test_save',win_compress, ...
  710.                                      '*.mat','Save Compressed Signal');
  711.         if ~ok, return; end
  713.         % Begin waiting.
  714.         %--------------
  715.         wwaiting('msg',win_compress,'Wait ... saving');
  717.         % Getting Analysis values.
  718.         %-------------------------
  719.         win_caller = wmemtool('rmb',win_compress,n_misc_loc,ind_win_caller);
  720.         wname = wmemtool('rmb',win_caller,n_param_anal,ind_wav_name);
  721.         thrDATA = wmemtool('rmb',win_compress,n_thrDATA,ind_value);
  722.         xc = thrDATA{1};
  723.         thrParams = thrDATA{4};
  724.         if length(thrParams)==1      % Global Mode
  725.             valTHR = thrParams;
  726.             thrName = 'valTHR';
  727.         else
  728.             thrName = 'thrParams';  % By Level Mode
  729.         end
  731.         % Saving file.
  732.         %--------------
  733.         [name,ext] = strtok(filename,'.');
  734.         if isempty(ext) | isequal(ext,'.')
  735.             ext = '.mat'; filename = [name ext];
  736.         end
  737.         saveStr = name;
  738.         eval([saveStr '= xc ;']);
  739.         wwaiting('off',win_compress);
  740.         try
  741.           save([pathname filename],saveStr,thrName,'wname');
  742.         catch
  743.           errargt(mfilename,'Save FAILED !','msg');
  744.         end
  746.     case 'save_cfs'
  747.         win_compress = varargin{1};
  749.         % Testing file.
  750.         %--------------
  751.         [filename,pathname,ok] = utguidiv('test_save',win_compress, ...
  752.                                      '*.mat','Save Coefficients (1D)');
  753.         if ~ok, return; end
  755.         % Begin waiting.
  756.         %--------------
  757.         wwaiting('msg',win_compress,'Wait ... saving coefficients');
  759.         % Getting Analysis values.
  760.         %-------------------------
  761.         win_caller = wmemtool('rmb',win_compress,n_misc_loc,ind_win_caller);
  762.         wname = wmemtool('rmb',win_caller,n_param_anal,ind_wav_name);
  763.         thrDATA = wmemtool('rmb',win_compress,n_thrDATA,ind_value);
  764.         coefs = thrDATA{2};
  765.         longs = thrDATA{3};
  766.         thrParams = thrDATA{4};
  767.         if length(thrParams)==1      % Global Mode
  768.             valTHR = thrParams;
  769.             thrName = 'valTHR';
  770.         else
  771.             thrName = 'thrParams';  % By Level Mode
  772.         end
  774.         % Saving file.
  775.         %--------------
  776.         [name,ext] = strtok(filename,'.');
  777.         if isempty(ext) | isequal(ext,'.')
  778.             ext = '.mat'; filename = [name ext];
  779.         end
  780.         saveStr = {'coefs','longs',thrName,'wname'};
  781.         wwaiting('off',win_compress);
  782.         try
  783.           save([pathname filename],saveStr{:});
  784.         catch
  785.           errargt(mfilename,'Save FAILED !','msg');
  786.         end
  788.     case 'save_dec'
  789.         % Testing file.
  790.         %--------------
  791.         [filename,pathname,ok] = utguidiv('test_save',win_compress, ...
  792.                                      '*.wa1','Save Wavelet Analysis (1D)');
  793.         if ~ok, return; end
  795.         % Begin waiting.
  796.         %--------------
  797.         wwaiting('msg',win_compress,'Wait ... saving decomposition');
  799.         % Getting Analysis parameters.
  800.         %-----------------------------
  801.         win_caller = wmemtool('rmb',win_compress,n_misc_loc,ind_win_caller);
  802.         [wave_name,data_name] = wmemtool('rmb',win_caller,n_param_anal,  ...
  803.                                                ind_wav_name,ind_sig_name);
  805.         % Getting Analysis values.
  806.         %-------------------------
  807.         thrDATA = wmemtool('rmb',win_compress,n_thrDATA,ind_value);
  808.         coefs = thrDATA{2};
  809.         longs = thrDATA{3};
  810.         thrParams = thrDATA{4};
  811.         if length(thrParams)==1      % Global Mode
  812.             valTHR = thrParams;
  813.             thrName = 'valTHR';
  814.         else
  815.             thrName = 'thrParams';  % By Level Mode
  816.         end
  818.         % Saving file.
  819.         %--------------
  820.         [name,ext] = strtok(filename,'.');
  821.         if isempty(ext) | isequal(ext,'.')
  822.             ext = '.wa1'; filename = [name ext];
  823.         end
  824.         saveStr = {'coefs','longs',thrName,'wave_name','data_name'};
  825.         wwaiting('off',win_compress);
  826.         try
  827.           save([pathname filename],saveStr{:});
  828.         catch
  829.           errargt(mfilename,'Save FAILED !','msg');
  830.         end
  832.     case 'close'
  833.          [status,win_caller] = wmemtool('rmb',win_compress,n_misc_loc,...
  834.                                              ind_status,ind_win_caller);
  835.         if status==1
  836.             % Test for Updating.
  837.             %--------------------
  838.             status = wwaitans(win_compress,...
  839.                               'Update the synthesized signal ?',2,'cancel');
  840.         end
  841.         switch status
  842.             case 1
  843.               wwaiting('msg',win_compress,'Wait ... computing');
  844.               thrDATA = wmemtool('rmb',win_compress,n_thrDATA,ind_value);
  845.               valTHR  = thrDATA{4};
  846.               wmemtool('wmb',win_caller,n_param_anal,...
  847.                        ind_ssig_type,'cs',ind_thr_val,valTHR);
  848.               hdl_datas = wmemtool('rmb',win_compress,n_misc_loc,ind_hdl_datas);
  849.               hdl_comp  = hdl_datas(2);
  850.               dw1dmngr('return_comp',win_caller,status,hdl_comp);
  851.               wwaiting('off',win_compress);
  853.             case 0 , dw1dmngr('return_comp',win_caller,status);
  854.         end
  855.         varargout{1} = status;
  857.     otherwise
  858.         errargt(mfilename,'Unknown Option','msg');
  859.         error('*');
  860. end