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

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

波变换

开发平台：

Matlab

dw1ddeno.m：源码内容
  1. function varargout = dw1ddeno(option,varargin)
  2. %DW1DDENO Wavelet 1-D de-noising.
  3. %   VARARGOUT = DW1DDENO(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.19.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. nbLOC_1_stored = 7;
  48. % MB2 (local).
  49. %-------------
  50. n_thrDATA = 'thrDATA';
  51. ind_value = 1;
  52. nbLOC_2_stored = 1;
  54. if ~isequal(option,'create') , win_denoise = varargin{1}; end
  55. switch option
  56.     case 'create'
  58.         % Get Globals.
  59.         %-------------
  60.         [Def_Btn_Height,Y_Spacing,Def_FraBkColor] = ...
  61.             mextglob('get','Def_Btn_Height','Y_Spacing','Def_FraBkColor');
  63.         % Calling figure.
  64.         %----------------
  65.         win_caller     = varargin{1};
  66.         str_win_caller = sprintf('%.0f',win_caller);
  68.         % Window initialization.
  69.         %----------------------
  70.         win_name = 'Wavelet 1-D  --  De-noising';
  71.         [win_denoise,pos_win,win_units,str_win_denoise,...
  72.                 frame0,pos_frame0,Pos_Graphic_Area,pus_close] = ...
  73.                     wfigmngr('create',win_name,'','ExtFig_CompDeno', ...
  74.                         strvcat(mfilename,'cond'),1,1,0);
  75.         set(win_denoise,'userdata',win_caller);
  76.         varargout{1} = win_denoise;
  78. % Add Help for Tool.
  79. %------------------
  80. wfighelp('addHelpTool',win_denoise,'Signal De-noising','DW1D_DENO_GUI');
  82. % Add Help Item.
  83. %----------------
  84. wfighelp('addHelpItem',win_denoise,'De-noising Procedure','DENO_PROCEDURE');
  85. wfighelp('addHelpItem',win_denoise,'Available Methods','COMP_DENO_METHODS');
  86. wfighelp('addHelpItem',win_denoise,'Variance Adaptive Thresholding','VARTHR');
  88.         % Menu construction for current figure.
  89.         %--------------------------------------
  90. m_save  = wfigmngr('getmenus',win_denoise,'save');
  91.         sav_menus(1) = uimenu(m_save,...
  92.                                 'Label','De-noised &Signal ',...
  93.                                 'Position',1,                    ...
  94.                                 'Enable','Off',                  ...
  95.                                 'Callback',                      ...
  96.                                 [mfilename '(''save_synt'','     ...
  97.                                         str_win_denoise ');']    ...
  98.                                 );
  99.         sav_menus(2) = uimenu(m_save,...
  100.                                'Label','&Coefficients ',  ...
  101.                                'Position',2,                  ...
  102.                                'Enable','Off',                ...
  103.                                'Callback',                    ...
  104.                                [mfilename '(''save_cfs'','    ...
  105.                                        str_win_denoise ');']  ...
  106.                                );
  107.         sav_menus(3) = uimenu(m_save,...
  108.                                 'Label','&Decomposition ', ...
  109.                                 'Position',3,                  ...
  110.                                 'Enable','Off',                ...
  111.                                 'Callback',                    ...
  112.                                 [mfilename '(''save_dec'','    ...
  113.                                         str_win_denoise ');']  ...
  114.                                 );
  116.         % Begin waiting.
  117.         %---------------
  118.         wwaiting('msg',win_denoise,'Wait ... initialization');
  120.         % Getting variables from dw1dtool figure memory block.
  121.         %-----------------------------------------------------
  122.         [Sig_Name,Wav_Name,Lev_Anal,Sig_Size] = ...
  123.                         wmemtool('rmb',win_caller,n_param_anal, ...
  124.                                 ind_sig_name,ind_wav_name,      ...
  125.                                 ind_lev_anal,ind_sig_size);
  126.         Wav_Fam = wavemngr('fam_num',Wav_Name);
  127.         [coefs,longs] = wmemtool('rmb',win_caller,n_coefs_longs, ...
  128.                                        ind_coefs,ind_longs);
  130.         % Parameters initialization.
  131.         %---------------------------
  132.         dy = Y_Spacing;
  134.         % Command part of the window.
  135.         %============================
  136.         
  137.         % Position property of objects.
  138.         %------------------------------
  139.         xlocINI = pos_frame0([1 3]);
  140.         ytopINI = pos_win(4)-dy;
  142.         % Data, Wavelet and Level parameters.
  143.         %------------------------------------
  144.         toolPos = utanapar('create_copy',win_denoise, ...
  145.                     {'xloc',xlocINI,'top',ytopINI},...
  146.                     {'n_s',{Sig_Name,Sig_Size},'wav',Wav_Name,'lev',Lev_Anal} ...
  147.                     );
  149.         % Threshold tool.
  150.         %----------------
  151.         ytopTHR = toolPos(2)-4*dy;
  152.         utthrw1d('create',win_denoise, ...
  153.                  'xloc',xlocINI,'top',ytopTHR,...
  154.                  'ydir',-1, ...
  155.                  'levmax',Lev_Anal, ...
  156.                  'levmaxMAX',Lev_Anal, ...
  157.                  'caller',mfilename, ...
  158.                  'toolOPT','deno' ...
  159.                  );
  161.         % Adding colormap GUI.
  162.         %---------------------
  163.         [viewType,hdl_cfs] = dw1dvdrv('get_imgcfs',win_caller);
  164.         if isequal(viewType,'image') & (Lev_Anal<8)
  165.             [pop_pal_caller,mapName,nbColors] = ...
  166.                 cbcolmap('get',win_caller,'pop_pal','mapName','nbColors');
  167.             utcolmap('create',win_denoise, ...
  168.                      'xloc',xlocINI, ...
  169.                      'bkcolor',Def_FraBkColor, ...
  170.                      'briflag',0, ...
  171.                      'enable','on');
  172.             pop_pal_loc = cbcolmap('get',win_denoise,'pop_pal');
  173.             set(pop_pal_loc,'Userdata',get(pop_pal_caller,'Userdata'));
  174.             cbcolmap('set',win_denoise,'pal',{mapName,nbColors});
  175.         end
  177.         % General graphical parameters initialization.
  178.         %--------------------------------------------
  179.         bdx      = 0.08*pos_win(3);
  180.         bdy      = 0.06*pos_win(4);
  181.         ecy      = 0.03*pos_win(4);
  182.         y_graph  = 2*Def_Btn_Height+dy;
  183.         h_graph  = pos_frame0(4)-y_graph;
  184.         w_graph  = pos_frame0(1);
  185.         fontsize = wmachdep('fontsize','normal',9,Lev_Anal);
  187.         % Axes construction parameters.
  188.         %------------------------------
  189.         w_left     = (w_graph-3*bdx)/2;
  190.         x_left     = bdx;
  191.         w_right    = w_left;
  192.         x_right    = x_left+w_left+5*bdx/4;
  193.         n_axeleft  = Lev_Anal;
  194.         n_axeright = 3;
  195.         ind_left   = n_axeleft;
  196.         ind_right  = n_axeright;
  198.         % Vertical separation.
  199.         %---------------------
  200.         w_fra = 0.01*pos_win(3);
  201.         x_fra = (w_graph-w_fra)/2;
  202.         uicontrol('Parent',win_denoise,...
  203.                   'Style','frame',...
  204.                   'Unit',win_units,...
  205.                   'Position',[x_fra,y_graph,w_fra,h_graph],...
  206.                   'Backgroundcolor',Def_FraBkColor ...
  207.                   );
  209.         % Building axes on the right part.
  210.         %---------------------------------
  211.         ecy_right = 2*ecy;
  212.         h_right =(h_graph-2*bdy-(n_axeright-1)*ecy_right)/n_axeright;
  213.         y_right = y_graph+bdy;
  214.         axe_datas = zeros(1,n_axeright);
  215.         pos_right = [x_right y_right w_right h_right];
  216.         for k = 1:n_axeright
  217.             axe_datas(k) = axes(...
  218.                                 'Parent',win_denoise, ...
  219.                                 'Units',win_units,     ...
  220.                                 'Position',pos_right,  ...
  221.                                 'Drawmode','fast',     ...
  222.                                 'Box','On'             ...
  223.                                 );
  224.             pos_right(2) = pos_right(2)+pos_right(4)+ecy_right;
  225.         end
  226.         set(axe_datas(1),'visible','off');
  228.         % Displaying the signal.
  229.         %-----------------------
  230.         sig = dw1dfile('sig',win_caller,1);
  231.         axeAct = axe_datas(3);
  232.         axes(axeAct)
  233.         curr_color = wtbutils('colors','sig');
  234.         lin_sig = line(...
  235.                        'Xdata',1:length(sig),...
  236.                        'Ydata',sig,...
  237.                        'erasemode','none',...
  238.                        'Color',curr_color,...
  239.                        'Parent',axeAct);
  240.         wtitle('Original signal','Parent',axeAct);
  241.         xlim = [1         Sig_Size];
  242.         ylim = [min(sig)  max(sig)];
  243.         if xlim(1)==xlim(2) , xlim = xlim+0.01*[-1 1]; end
  244.         if ylim(1)==ylim(2) , ylim = ylim+0.01*[-1 1]; end
  245.         set(axeAct,'Xlim',xlim,'Ylim',ylim);
  246.         utthrw1d('set',win_denoise,'handleORI',lin_sig);
  248.         % Displaying original details coefficients.
  249.         %------------------------------------------
  250.         axeAct = axe_datas(2);
  251.         [viewType,hdl_cfs] = dw1dvdrv('get_imgcfs',win_caller);
  252.         [details,set_ylim,ymin,ymax] = dw1dfile('cfs_beg',win_caller,...
  253.                             [1:Lev_Anal],1);
  254.         
  255.         if isequal(viewType,'image')
  256.             flagType = 1;
  257.             cfsMode  = [];
  258.             set(win_denoise,'Colormap',get(win_caller,'Colormap'));
  259.             [nul,i_min] = min(abs(details(:)));
  260.             if ~isempty(hdl_cfs)
  261.                 col_cfs = flipud(get(hdl_cfs,'Cdata'));
  262.             else
  263.                 col_cfs = wcodemat(details,def_nbCodeOfColors,'row',1);
  264.             end
  265.             col_min = col_cfs(i_min);
  266.             col_cfs = flipud(col_cfs);
  267.             axes(axeAct);
  268.             cfs_ori = image(col_cfs,'Parent',axeAct,'Userdata',col_min);
  269.             clear col_cfs
  270.             levlab = int2str([Lev_Anal:-1:1]');
  271.         else
  272.             flagType = -1;
  273.             axes(axeAct);
  274.             hdl_stem = copyobj(hdl_cfs,axeAct);
  275.             set(hdl_stem,'Visible','on');
  276.             cfsMode = get(hdl_stem(1),'Userdata');
  277.             levlab  = int2str([1:Lev_Anal]');
  278.             cfs_ori = '';
  279.         end
  280.         set(axeAct,...
  281.               'Userdata',cfs_ori,       ...
  282.               'Xlim',[1 Sig_Size],      ...
  283.               'YTicklabelMode','manual',...
  284.               'YTick',[1:Lev_Anal],     ...
  285.               'YTickLabel',levlab,      ...
  286.               'Ylim',[0.5 Lev_Anal+0.5] ...
  287.               );
  289.         wtitle('Original coefficients','Parent',axeAct);
  290.         wylabel('Level number','Parent',axeAct);
  292.         xylim = get(axeAct,{'Xlim','Ylim'});        
  293.         set(axe_datas(1),'Xlim',xylim{1},'Ylim',xylim{2});
  295.         % Building axes on the left part.
  296.         %--------------------------------
  297.         ecy_left = ecy/2;
  298.         h_left   = (h_graph-2*bdy-(n_axeleft-1)*ecy_left)/n_axeleft;
  299.         y_left   = y_graph+0.75*bdy;
  301.         axe_left = zeros(1,n_axeleft);
  302.         txt_left = zeros(1,n_axeleft);
  303.         pos_left = [x_left y_left w_left h_left];
  304.         commonProp = {...
  305.            'Parent',win_denoise,...
  306.            'Units',win_units,...
  307.            'Drawmode','fast',...
  308.            'Box','On'...
  309.            };
  310.         for k = 1:n_axeleft
  311.             if k~=1
  312.                 axe_left(k) = axes(commonProp{:}, ...
  313.                                   'Position',pos_left,...
  314.                                   'XTicklabelMode','manual','XTickLabel',[]);
  315.             else
  316.                 axe_left(k) = axes(commonProp{:},'Position',pos_left);
  317.             end
  318.             pos_left(2) = pos_left(2)+pos_left(4)+ecy_left;
  320.             txt_left(k) = txtinaxe('create',['d' wnsubstr(k)],...
  321.                                     axe_left(k),'left',...
  322.                                     'on','bold',fontsize);
  323.         end
  324.         utthrw1d('set',win_denoise,'axes',axe_left);
  326.         % Initializing by level threshold.
  327.         %---------------------------------
  328.         maxTHR = zeros(1,Lev_Anal);
  329.         for k = 1:Lev_Anal , maxTHR(k) = max(abs(details(k,:))); end
  330.         valTHR = dw1ddeno('compute_LVL_THR',win_denoise,win_caller);
  331.         valTHR = min(valTHR,maxTHR);
  333.         % Displaying details.
  334.         %-------------------
  335.         col_det = wtbutils('colors','det',Lev_Anal);
  336.         for k = Lev_Anal:-1:1
  337.             axeAct  = axe_left(ind_left);
  338.             axes(axeAct);
  339.             lin_cfs(k) = line(...
  340.                              'Parent',axeAct,...
  341.                              'Xdata',1:Sig_Size,...
  342.                              'Ydata',details(k,:),...
  343.                              'Color',col_det(k,:));
  344.             utthrw1d('plot_dec',win_denoise,k, ...
  345.                      {maxTHR(k),valTHR(k),1,Sig_Size,k})
  346.             maxi = max([abs(ymax(k)),abs(ymin(k))]);
  347.             if abs(maxi)<eps , maxi = maxi+0.01; end;
  348.             ylim = 1.1*[-maxi maxi];
  349.             set(axe_left(ind_left),'Xlim',xlim,'Ylim',ylim);
  350.             ind_left = ind_left-1;
  351.         end
  352.         axeAct = axe_left(Lev_Anal);
  353.         wtitle('Original details coefficients','Parent',axeAct);
  355.         % Axes attachment.
  356.         %-----------------
  357.         if flagType~=-1
  358.             axe_cmd = [axe_datas(1:3) axe_left(1:n_axeleft)];
  359.             axe_act = [];
  360.             axe_cfs = [axe_datas(1:2)];
  361.         else
  362.             axe_cmd = [axe_datas(1:3) axe_left(1:n_axeleft)];
  363.             axe_act = [];
  364.             axe_cfs = [axe_datas(1:2)];
  365.         end
  366.         dynvtool('init',win_denoise,[],axe_cmd,axe_act,[1 0], ...
  367.                         '','','dw1dcoor',[win_caller,axe_cfs,flagType*Lev_Anal]);
  369. %        dynvtool('init',win_denoise,[],axe_cmd,axe_act,[1 0], ...
  370. %                        '','','dw1dcoor',[axe_cfs flagType*Lev_Anal],...
  371. %                        'cbthrw1d',[axe_left(1:n_axeleft)]);
  373.         % Initialization of  Denoising structure.
  374.         %----------------------------------------
  375.         xmin = 1; xmax = Sig_Size;
  376.         utthrw1d('set',win_denoise,...
  377.                        'thrstruct',{xmin,xmax,valTHR,lin_cfs},...
  378.                        'intdepthr',[]);
  381.         % Memory blocks update.
  382.         %----------------------
  383.         wmemtool('ini',win_denoise,n_misc_loc,nbLOC_1_stored);
  384.         wmemtool('wmb',win_denoise,n_misc_loc,    ...
  385.                        ind_sav_menus,sav_menus,   ...
  386.                        ind_status,0,              ...
  387.                        ind_win_caller,win_caller, ...
  388.                        ind_axe_datas,axe_datas,   ...
  389.                        ind_cfsMode,cfsMode,       ...
  390.                        ind_lin_cfs,lin_cfs        ...
  391.                        );
  392.         wmemtool('ini',win_denoise,n_thrDATA,nbLOC_2_stored);
  394.         % Setting units to normalized.
  395.         %-----------------------------
  396.         wfigmngr('normalize',win_denoise);
  398.         % End waiting.
  399.         %-------------
  400.         utthrw1d('enable',win_denoise,'on',[1:Lev_Anal]);
  401.         wwaiting('off',win_denoise);
  403.     case 'denoise'
  405.         % Waiting message.
  406.         %-----------------
  407.         wwaiting('msg',win_denoise,'Wait ... computing');
  409.         % Clear & Get Handles.
  410.         %----------------------
  411.         dw1ddeno('clear_GRAPHICS',win_denoise);
  412.         win_caller = wmemtool('rmb',win_denoise,n_misc_loc,ind_win_caller);
  413.         axe_datas = wmemtool('rmb',win_denoise,n_misc_loc,ind_axe_datas);
  415.         % Getting memory blocks.
  416.         %-----------------------
  417.         [Wav_Name,Lev_Anal,Sig_Size] = ...
  418.                         wmemtool('rmb',win_caller,n_param_anal, ...
  419.                                        ind_wav_name,ind_lev_anal,ind_sig_size);
  420.         [coefs,longs] = wmemtool('rmb',win_caller,n_coefs_longs, ...
  421.                                        ind_coefs,ind_longs);
  424.         % De-noising depending on the selected thresholding mode.
  425.         %--------------------------------------------------------
  426.         cxc = utthrw1d('den_M2',win_denoise,coefs,longs);
  427.         lxc = longs;
  428.         xc  = waverec(cxc,longs,Wav_Name);
  430.         % Displaying denoised signal.
  431.         %----------------------------
  432.         lin_den = utthrw1d('get',win_denoise,'handleTHR');
  433.         if ishandle(lin_den)
  434.             set(lin_den,'Ydata',xc,'Visible','on');
  435.         else
  436.             curr_color = wtbutils('colors','ssig');
  437.             lin_den = line(...
  438.                            'Parent',axe_datas(3), ...
  439.                            'Xdata',1:length(xc),  ...
  440.                            'Ydata',xc,            ...
  441.                            'color',curr_color     ...
  442.                            );
  443.             utthrw1d('set',win_denoise,'handleTHR',lin_den);
  444.         end     
  445.         wtitle('Original and de-noised signals','Parent',axe_datas(3));
  447.         % Displaying thresholded details coefficients.
  448.         %---------------------------------------------
  449.         cfsMode = wmemtool('rmb',win_denoise,n_misc_loc,ind_cfsMode);
  450.         if isempty(cfsMode)
  451.             col_cfs   = wrepcoef(cxc,lxc);
  452.             nz_cfs    = find(col_cfs~=0);
  453.             [nbr,nbc] = size(col_cfs);
  454.             img_cfs   = get(axe_datas(2),'Userdata');
  455.             col_min   = get(img_cfs,'Userdata');
  456.             cfs_ori   = flipud(get(img_cfs,'Cdata'));
  457.             col_cfs   = col_min*ones(nbr,nbc);
  458.             col_cfs(nz_cfs) = cfs_ori(nz_cfs);
  459.             image(flipud(col_cfs),'Parent',axe_datas(1));
  460.             levlab = int2str([Lev_Anal:-1:1]');
  461.         else
  462.             dw1dstem(axe_datas(1),cxc,lxc,'mode',cfsMode,'colors','WTBX');
  463.             levlab = int2str([1:Lev_Anal]');
  464.         end
  465.         set(axe_datas(1),...
  466.                 'clipping','on',            ...
  467.                 'Xlim',get(axe_datas(2),'Xlim'),...                    
  468.                 'YTicklabelMode','manual',  ...
  469.                 'YTick',[1:Lev_Anal],       ...
  470.                 'YTickLabel',levlab,        ...
  471.                 'Ylim',[0.5 Lev_Anal+0.5]   ...
  472.                 );
  473.         wtitle('Thresholded coefficients','Parent',axe_datas(1));
  474.         wylabel('Level number','Parent',axe_datas(1));
  475.         set(findobj(axe_datas(1)),'Visible','on');
  477.         % Memory blocks update.
  478.         %----------------------
  479.         thrStruct = utthrw1d('get',win_denoise,'thrstruct');
  480.         thrParams = {thrStruct(1:Lev_Anal).thrParams};
  481.         wmemtool('wmb',win_denoise,n_thrDATA,ind_value,{xc,cxc,lxc,thrParams});
  482.         dw1ddeno('enable_menus',win_denoise,'on');
  484.         % End waiting.
  485.         %-------------
  486.         wwaiting('off',win_denoise);
  488.     case 'compute_LVL_THR'
  489.         win_caller = varargin{2};
  490.         [numMeth,meth,alfa,sorh] = utthrw1d('get_LVL_par',win_denoise);
  491.         [coefs,longs] = wmemtool('rmb',win_caller,n_coefs_longs,...
  492.                                        ind_coefs,ind_longs);
  493.         varargout{1} = wthrmngr('dw1ddenoLVL',meth,coefs,longs,alfa);
  495.     case 'update_LVL_meth'
  496.         dw1ddeno('clear_GRAPHICS',win_denoise);
  497.         win_caller = wmemtool('rmb',win_denoise,n_misc_loc,ind_win_caller);
  498.         valTHR = dw1ddeno('compute_LVL_THR',win_denoise,win_caller);
  499.         utthrw1d('update_LVL_meth',win_denoise,valTHR);
  501.     case 'clear_GRAPHICS'
  502.         status = wmemtool('rmb',win_denoise,n_misc_loc,ind_status);
  503.         if status == 0 , return; end
  505.         % Diseable Toggle and Menus.
  506.         %---------------------------
  507.         dw1ddeno('enable_menus',win_denoise,'off');
  509.         % Get Handles.
  510.         %-------------
  511.         axe_datas = wmemtool('rmb',win_denoise,n_misc_loc,ind_axe_datas);
  513.         % Setting the de-noised coefs axes invisible.
  514.         %--------------------------------------------
  515.         lin_den = utthrw1d('get',win_denoise,'handleTHR');
  516.         if ~isempty(lin_den)
  517.            vis = get(lin_den,'Visible');
  518.            if vis(1:2)=='on'
  519.                set(findobj(axe_datas(1)),'Visible','off');
  520.                axes(axe_datas(3));
  521.                wtitle('Original signal','Parent',axe_datas(3));
  522.                set(lin_den,'Visible','off');
  523.            end
  524.         end
  526.     case 'enable_menus'
  527.         enaVal = varargin{2};
  528.         sav_menus = wmemtool('rmb',win_denoise,n_misc_loc,ind_sav_menus);
  529.         set(sav_menus,'Enable',enaVal);
  530.         utthrw1d('enable_tog_res',win_denoise,enaVal);     
  531.         if strncmpi(enaVal,'on',2) , status = 1; else , status = 0; end
  532.         wmemtool('wmb',win_denoise,n_misc_loc,ind_status,status);
  534.     case 'save_synt'
  536.         % Testing file.
  537.         %--------------
  538.         [filename,pathname,ok] = utguidiv('test_save',win_denoise, ...
  539.                                      '*.mat','Save De-noised Signal');
  540.         if ~ok, return; end
  542.         % Begin waiting.
  543.         %--------------
  544.         wwaiting('msg',win_denoise,'Wait ... saving');
  546.         % Getting Analysis values.
  547.         %-------------------------
  548.         win_caller = wmemtool('rmb',win_denoise,n_misc_loc,ind_win_caller);
  549.         wname = wmemtool('rmb',win_caller,n_param_anal,ind_wav_name);
  550.         thrDATA = wmemtool('rmb',win_denoise,n_thrDATA,ind_value);
  551.         xc = thrDATA{1};
  552.         thrParams = thrDATA{4};
  554.         % Saving file.
  555.         %--------------
  556.         [name,ext] = strtok(filename,'.');
  557.         if isempty(ext) | isequal(ext,'.')
  558.             ext = '.mat'; filename = [name ext];
  559.         end
  560.         saveStr = name;
  561.         eval([saveStr '= xc ;']);
  562.         wwaiting('off',win_denoise);
  563.         try
  564.           save([pathname filename],saveStr,'thrParams','wname');
  565.         catch          
  566.           errargt(mfilename,'Save FAILED !','msg');
  567.         end
  569.     case 'save_cfs'
  571.         % Testing file.
  572.         %--------------
  573.         [filename,pathname,ok] = utguidiv('test_save',win_denoise, ...
  574.                                      '*.mat','Save Coefficients (1D)');
  575.         if ~ok, return; end
  577.         % Begin waiting.
  578.         %--------------
  579.         wwaiting('msg',win_denoise,'Wait ... saving coefficients');
  581.         % Getting Analysis values.
  582.         %-------------------------
  583.         win_caller = wmemtool('rmb',win_denoise,n_misc_loc,ind_win_caller);
  584.         wname = wmemtool('rmb',win_caller,n_param_anal,ind_wav_name);
  585.         thrDATA = wmemtool('rmb',win_denoise,n_thrDATA,ind_value);
  586.         coefs = thrDATA{2};
  587.         longs = thrDATA{3};
  588.         thrParams = thrDATA{4};
  590.         % Saving file.
  591.         %--------------
  592.         [name,ext] = strtok(filename,'.');
  593.         if isempty(ext) | isequal(ext,'.')
  594.             ext = '.mat'; filename = [name ext];
  595.         end
  596.         saveStr = {'coefs','longs','thrParams','wname'};
  597.         wwaiting('off',win_denoise);
  598.         try
  599.           save([pathname filename],saveStr{:});
  600.         catch          
  601.           errargt(mfilename,'Save FAILED !','msg');
  602.         end
  604.     case 'save_dec'
  606.         % Testing file.
  607.         %--------------
  608.         [filename,pathname,ok] = utguidiv('test_save',win_denoise, ...
  609.                                      '*.wa1','Save Wavelet Analysis (1D)');
  610.         if ~ok, return; end
  612.         % Begin waiting.
  613.         %--------------
  614.         wwaiting('msg',win_denoise,'Wait ... saving decomposition');
  616.         % Getting Analysis values.
  617.         %-------------------------
  618.         win_caller = wmemtool('rmb',win_denoise,n_misc_loc,ind_win_caller);
  619.         [wave_name,data_name] = wmemtool('rmb',win_caller,n_param_anal, ...
  620.                                           ind_wav_name,ind_sig_name);
  621.         thrDATA = wmemtool('rmb',win_denoise,n_thrDATA,ind_value);
  622.         coefs = thrDATA{2};
  623.         longs = thrDATA{3};
  624.         thrParams = thrDATA{4};
  626.         % Saving file.
  627.         %--------------
  628.         [name,ext] = strtok(filename,'.');
  629.         if isempty(ext) | isequal(ext,'.')
  630.             ext = '.wa1'; filename = [name ext];
  631.         end
  632.         saveStr = {'coefs','longs','thrParams','wave_name','data_name'};
  633.         wwaiting('off',win_denoise);
  634.         try
  635.           save([pathname filename],saveStr{:});
  636.         catch          
  637.           errargt(mfilename,'Save FAILED !','msg');
  638.         end
  640.     case 'close'
  641.         [status,win_caller] = wmemtool('rmb',win_denoise,n_misc_loc, ...
  642.                                              ind_status,ind_win_caller);
  643.         if status==1
  644.             % Test for Updating.
  645.             %--------------------
  646.             status = wwaitans(win_denoise,...
  647.                               'Update the synthesized signal ?',2,'cond');
  648.         end
  649.         switch status
  650.           case 1
  651.             wwaiting('msg',win_denoise,'Wait ... computing');
  652.             thrDATA = wmemtool('rmb',win_denoise,n_thrDATA,ind_value);
  653.             valTHR  = thrDATA{4};
  654.             lin_den = utthrw1d('get',win_denoise,'handleTHR');
  655.             wmemtool('wmb',win_caller,n_param_anal,...
  656.                      ind_ssig_type,'ds',ind_thr_val,valTHR);
  657.             dw1dmngr('return_deno',win_caller,status,lin_den);
  658.             wwaiting('off',win_denoise);
  660.           case 0 , dw1dmngr('return_deno',win_caller,status);
  661.         end
  662.         varargout{1} = status;
  664.     otherwise
  665.         errargt(mfilename,'Unknown Option','msg');
  666.         error('*');
  667. end
  670. %-------------------------------------------------
  671. function varargout = depOfMachine(varargin)
  673. btn_height = varargin{1};
  674. scrSize = get(0,'ScreenSize');
  675. if scrSize(4)<600
  676.     hvisu = btn_height;
  677. elseif scrSize(4)<700
  678.     hvisu = 3*btn_height/2;
  679. else
  680.     hvisu = 2*btn_height;
  681. end
  682. varargout{1} = hvisu;
  683. %-------------------------------------------------