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

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

波变换

开发平台：

Matlab

utthrw1d.m：源码内容
  1. function varargout = utthrw1d(option,fig,varargin)
  2. %UTTHRW1D Utilities for wavelet thresholding 1-D.
  3. %   VARARGOUT = UTTHRW1D(OPTION,FIG,VARARGIN)
  5. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 01-Jun-98.
  6. %   Last Revision: 22-May-2003.
  7. %   Copyright 1995-2004 The MathWorks, Inc.
  8. %   $Revision: 1.11.4.2 $  $Date: 2004/03/15 22:42:11 $
  10. % MemBloc of stored values.
  11. %--------------------------
  12. n_memblocTHR   = 'MB_ThrStruct';
  13. ind_thr_struct = 1;
  14. ind_int_thr    = 2;
  16. % Default values.
  17. %----------------
  18. max_lev_anal = 12;
  19. def_lev_anal = 5;
  21. % Tag property of objects.
  22. %-------------------------
  23. tag_fra_tool   = ['Fra_' mfilename];
  24. tag_lineH_up   = 'LH_u';
  25. tag_lineH_down = 'LH_d';
  26. tag_lineV      = 'LV';
  28. switch option
  29.   case {'create'}
  31.   otherwise
  32.     % ud.handlesUIC = ...
  33.     %   [fra_utl;txt_top;pop_met; ...
  34.     %    rad_sof;rad_har;txt_noi;pop_noi; ...
  35.     %    txt_BMS;sli_BMS;txt_tit(1:4),...
  36.     %    txt_nor;edi_nor;txt_npc; ...
  37.     %    txt_zer;edi_zer;txt_zpc; ...
  38.     %    tog_thr;tog_res;pus_est];
  39.     %-----------------------------------------
  40.     if ~ishandle(fig) , varargout{1} = []; return; end
  41.     uic = findobj(get(fig,'Children'),'flat','type','uicontrol');
  42.     fra = findobj(uic,'style','frame','tag',tag_fra_tool);
  43.     if isempty(fra) , return; end
  44.     calledFUN = wfigmngr('getWinPROP',fig,'calledFUN');
  45.     ud  = get(fra,'Userdata');
  46.     toolOPT = ud.toolOPT;
  47.     toolStatus = ud.status;
  48.     handlesUIC = ud.handlesUIC;
  49.     h_CMD_LVL  = ud.h_CMD_LVL;
  50.     h_GRA_LVL  = ud.h_GRA_LVL;
  51.     switch option
  52.       case 'handles'
  53.         handles = [handlesUIC(:);h_CMD_LVL(:)];
  54.         varargout{1} = handles(ishandle(handles));
  55.         return;
  57.       case {'handlesUIC','handlesuic'}
  58.         varargout{1} = handlesUIC;
  59.         return;
  61.       case {'h_CMD_LVL','h_cmd_lvl'}
  62.         varargout{1} = h_CMD_LVL;
  63.         return;
  64.     end
  65.     ind = 2;
  66.     txt_top = handlesUIC(ind); ind = ind+1;
  67.     pop_met = handlesUIC(ind); ind = ind+1;
  68.     rad_sof = handlesUIC(ind); ind = ind+1;
  69.     rad_har = handlesUIC(ind); ind = ind+1;
  70.     txt_noi = handlesUIC(ind); ind = ind+1;
  71.     pop_noi = handlesUIC(ind); ind = ind+1;
  72.     txt_BMS = handlesUIC(ind); ind = ind+1;  
  73.     sli_BMS = handlesUIC(ind); ind = ind+1;
  74.     txt_tit(1:4) = handlesUIC(ind:ind+3); ind = ind+4;
  75.     txt_nor = handlesUIC(ind); ind = ind+1;
  76.     edi_nor = handlesUIC(ind); ind = ind+1;
  77.     txt_npc = handlesUIC(ind); ind = ind+1;
  78.     txt_zer = handlesUIC(ind); ind = ind+1;
  79.     edi_zer = handlesUIC(ind); ind = ind+1;
  80.     txt_zpc = handlesUIC(ind); ind = ind+1;
  81.     tog_thr = handlesUIC(ind); ind = ind+1;
  82.     tog_res = handlesUIC(ind); ind = ind+1;
  83.     pus_est = handlesUIC(ind); ind = ind+1;
  84. end
  86. switch option
  87.     case 'create'
  88.         % Get Globals.
  89.         %--------------
  90.         [Def_Txt_Height,Def_Btn_Height,Def_Btn_Width,Pop_Min_Width, ...
  91.          X_Spacing,Y_Spacing,sliYProp,...
  92.          Def_FraBkColor,Def_EdiBkColor] = ...
  93.             mextglob('get',...
  94.                 'Def_Txt_Height','Def_Btn_Height','Def_Btn_Width',   ...
  95.                 'Pop_Min_Width','X_Spacing','Y_Spacing','Sli_YProp', ...
  96.                 'Def_FraBkColor','Def_EdiBkColor'   ...
  97.                 );
  99.         % Defaults.
  100.         %----------
  101.         xleft = Inf; xright  = Inf; xloc = Inf;
  102.         ytop  = Inf; ybottom = Inf; yloc = Inf;
  103.         bkColor = Def_FraBkColor;
  104.         enaVal  = 'Off';
  105.         %------------------------
  106.         ydir   = -1;
  107.         levmin = 1;
  108.         levmax = def_lev_anal;
  109.         levmaxMAX = max_lev_anal;
  110.         levANAL = def_lev_anal;
  111.         visVal = 'on';
  112.         isbior = 0;
  113.         statusINI = 'On';
  114.         toolOPT   = 'deno';
  115.         %------------------------
  117.         % Inputs.
  118.         %--------        
  119.         nbarg = length(varargin);
  120.         for k=1:2:nbarg
  121.             arg = lower(varargin{k});
  122.             switch arg
  123.               case 'left'     , xleft     = varargin{k+1};
  124.               case 'right'    , xright    = varargin{k+1};
  125.               case 'xloc'     , xloc      = varargin{k+1};
  126.               case 'bottom'   , ybottom   = varargin{k+1};
  127.               case 'top'      , ytop      = varargin{k+1};
  128.               case 'yloc'     , yloc      = varargin{k+1};
  129.               case 'bkcolor'  , bkColor   = varargin{k+1};
  130.               case 'visible'  , visVal    = varargin{k+1};
  131.               case 'enable'   , enaVal    = varargin{k+1};
  132.               case 'isbior'   , isbior    = varargin{k+1};
  133.               case 'ydir'     , ydir      = varargin{k+1};
  134.               case 'levmin'   , levmin    = varargin{k+1};
  135.               case 'levmax'   , levmax    = varargin{k+1};
  136.               case 'levmaxmax', levmaxMAX = varargin{k+1};
  137.               case 'levanal'  , levANAL   = varargin{k+1};
  138.               case 'status'   , statusINI = varargin{k+1};
  139.               case 'toolopt'  , toolOPT   = varargin{k+1};
  140.             end 
  141.         end
  143.         % Structure initialization.
  144.         %--------------------------
  145.         % h_CMD_LVL: [txt_lev ; pop_int; sli_lev ; edi_lev] x Level
  146.         % h_GRA_LVL: [axe_thr ; lin_min ; lin_max] x Level
  147.         h_CMD_LVL = NaN*ones(4,levmaxMAX);        
  148.         h_GRA_LVL = NaN*ones(3,levmaxMAX);
  149.         thrStruct = struct(...
  150.                         'thrParams',cell(levmaxMAX,1), ...
  151.                         'hdlLines',cell(levmaxMAX,1)   ...                        '
  152.                         );
  153.         wmemtool('wmb',fig,n_memblocTHR, ...
  154.                   ind_thr_struct,thrStruct,ind_int_thr,[]);
  155.         ud = struct(...
  156.                 'toolOPT',toolOPT, ...
  157.                 'status',statusINI,...
  158.                 'levmin',levmin, ...
  159.                 'levmax',levmax, ...
  160.                 'levmaxMAX',levmaxMAX, ...
  161.                 'levANAL',levANAL, ...
  162.                 'visible',lower(visVal),...
  163.                 'ydir', ydir,    ...
  164.                 'isbior',isbior, ...
  165.                 'handlesUIC',[], ...
  166.                 'h_CMD_LVL',h_CMD_LVL, ...
  167.                 'h_GRA_LVL',h_GRA_LVL, ...
  168.                 'handleORI' ,[], ... 
  169.                 'handleTHR',[],  ... 
  170.                 'handleRES' ,[]  ... 
  171.                 );
  173.         % Figure units.
  174.         %--------------
  175.         str_numfig = num2mstr(fig);
  176.         old_units  = get(fig,'units');
  177.         fig_units  = 'pixels';
  178.         if ~isequal(old_units,fig_units), set(fig,'units',fig_units); end       
  180.         % Positions utilities.
  181.         %---------------------
  182.         nblevs = abs(levmax-levmin)+1;        
  183.         dx = X_Spacing; bdx = 3;
  184.         dy = Y_Spacing; 
  185.         d_txt  = (Def_Btn_Height-Def_Txt_Height);
  186.         deltaY = (Def_Btn_Height+dy);
  187.         [bdy,d_lev] = get_DLEV(nblevs,toolOPT);
  188.         sli_hi = Def_Btn_Height*sliYProp;
  189.         sli_dy = 0.5*Def_Btn_Height*(1-sliYProp);
  191.         % Setting frame position.
  192.         %------------------------
  193.         if ~isequal(toolOPT,'esti_DEN') , multiY = 1; else , multiY = 1; end
  194.         switch  toolOPT
  195.           case {'deno','esti'} , NB_Height = 6;
  196.           case {'esti_REG'}    , NB_Height = 6;
  197.           case {'esti_DEN'}    , NB_Height = 4+2*multiY;
  198.           case {'comp'}        , NB_Height = 7;
  199.         end
  200.         w_fra   = wfigmngr('get',fig,'fra_width');
  201.         h_fra   = (levmaxMAX+NB_Height)*Def_Btn_Height+...
  202.                    levmaxMAX*d_lev+ ...
  203.                    depOfMachine(Def_Btn_Height)+(NB_Height-1)*bdy;
  204.         xleft   = utposfra(xleft,xright,xloc,w_fra);
  205.         ybottom = utposfra(ybottom,ytop,yloc,h_fra);
  206.         pos_fra = [xleft,ybottom,w_fra,h_fra];
  208.         % String properties.
  209.         %-------------------
  210.         str_txt_top = 'Select thresholding method';
  211.         str_txt_BMS = '-       Sparsity       +';
  212.         str_txt_tit = strvcat('Lev','Int','Select','Thresh');
  213.         str_tog_thr = ['Int. dependent threshold settings'];
  214.         str_tog_res = 'Residuals';
  215.         str_pop_met = wthrmeth(toolOPT,'names');
  216.         switch  toolOPT
  217.           case {'deno','esti','esti_REG'}
  218.             str_rad_sof = 'soft';
  219.             str_rad_har = 'hard';
  220.             str_txt_noi = 'Select noise structure';
  221.             str_pop_noi = strvcat(...
  222.                 'Unscaled white noise', ...
  223.                 'Scaled white noise',   ...
  224.                 'Non-white noise'       ...
  225.                 );
  227.           case {'esti_DEN'}
  228.             str_rad_sof = 'soft';
  229.             str_rad_har = 'hard';
  230.             str_txt_noi = 'Select noise structure';
  231.             str_pop_noi = strvcat(...
  232.                 'Unscaled white noise', ...
  233.                 'Scaled white noise',   ...
  234.                 'Non-white noise'       ...
  235.                 );
  236.           case 'comp'
  237.             if ud.isbior
  238.                 str_txt_nor = 'Norm cfs recovery';
  239.             else
  240.                 str_txt_nor = 'Retained energy';
  241.             end
  242.             str_txt_zer = 'Number of zeros';
  243.         end
  244.         switch  toolOPT
  245.           case 'deno' , str_pus_est = 'De-noise'; estOPT = 'denoise';
  246.           case 'comp' , str_pus_est = 'Compress'; estOPT = 'compress';
  247.           case {'esti','esti_REG','esti_DEN'} ,
  248.                         str_pus_est = 'Estimate'; estOPT = 'estimate';
  249.         end
  251.         % Position properties.
  252.         %---------------------
  253.         txt_width   = Def_Btn_Width;
  254.         dy_lev      = Def_Btn_Height+d_lev;
  255.         xleft       = xleft+bdx;
  256.         w_rem       = w_fra-2*bdx;
  257.         ylow        = ybottom+h_fra-Def_Btn_Height-bdy;
  259.         w_uic       = (5*txt_width)/2;
  260.         x_uic       = xleft+(w_rem-w_uic)/2;
  261.         y_uic       = ylow;
  262.         pos_txt_top = [x_uic, y_uic+d_txt/2, w_uic, Def_Txt_Height];
  263.         y_uic       = y_uic-Def_Btn_Height;
  264.         pos_pop_met = [x_uic, y_uic, w_uic, Def_Btn_Height];
  265.         y_uic       = y_uic-Def_Btn_Height;
  266.         switch toolOPT
  267.           case {'deno','esti','esti_REG','esti_DEN'}
  268.             y_uic       = y_uic-bdy;
  269.             w_rad       = Pop_Min_Width;
  270.             w_sep       = (w_uic-2*w_rad)/3;
  271.             x_rad       = x_uic+w_sep;
  272.             pos_rad_sof = [x_rad, y_uic, w_rad, Def_Btn_Height];
  273.             x_rad       = x_rad+w_rad+w_sep;
  274.             pos_rad_har = [x_rad, y_uic, w_rad, Def_Btn_Height];
  276.             y_uic       = y_uic-multiY*Def_Btn_Height;
  277.             y_BMS       = y_uic;
  278.             pos_txt_noi = [x_uic, y_uic, w_uic, Def_Txt_Height];
  279.             y_uic       = y_uic-multiY*Def_Btn_Height;
  280.             pos_pop_noi = [x_uic, y_uic, w_uic, Def_Btn_Height];
  282.           case 'comp'
  283.             y_BMS       = y_uic;
  284.             y_uic       = y_uic-Def_Btn_Height;
  285.         end
  287.         pos_txt_BMS = [x_uic, y_BMS, w_uic, Def_Txt_Height];
  288.         y_BMS       = y_BMS-multiY*Def_Btn_Height;
  289.         w_BMS       = (w_uic-bdx)/3;
  290.         pos_sli_BMS = [x_uic+w_BMS/2, y_BMS+sli_dy, 2*w_BMS, sli_hi];
  292.         wx          = 2;
  293.         wbase       = 2*(w_rem-5*wx)/5;
  294.         w_lev       = [4*wbase ; 7*wbase ; 12*wbase ; 7*wbase]/12;
  295.         x_uic       = xleft+wx;
  296.         y_uic       = y_uic-Def_Btn_Height;
  297.         pos_lev_tit = [x_uic, y_uic, w_lev(1), Def_Txt_Height];
  298.         pos_lev_tit = pos_lev_tit(ones(1,4),:);
  299.         pos_lev_tit(:,3) = w_lev; 
  300.         for k=1:3 , pos_lev_tit(k+1,1) = pos_lev_tit(k,1)+pos_lev_tit(k,3); end
  301.         y_uic = pos_lev_tit(1,2)-levmaxMAX*(Def_Btn_Height+d_lev);
  303.         switch  toolOPT
  304.           case {'deno','esti','esti_REG','esti_DEN'}
  305.           case {'comp'}
  306.             wid1 = (15*w_rem)/26;
  307.             wid2 = (8*w_rem)/26;
  308.             wid3 = (2*w_rem)/26;
  309.             wx   = (w_rem-wid1-wid2-wid3)/4;
  310.             y_uic       = y_uic-Def_Btn_Height-bdy;
  311.             pos_txt_nor = [xleft, y_uic+d_txt/2, wid1, Def_Txt_Height];
  312.             x_uic       = pos_txt_nor(1)+pos_txt_nor(3)+wx;
  313.             pos_edi_nor = [x_uic, y_uic, wid2, Def_Btn_Height];
  314.             x_uic       = pos_edi_nor(1)+pos_edi_nor(3)+wx;
  315.             pos_txt_npc = [x_uic, y_uic+d_txt/2, wid3, Def_Txt_Height];
  317.             y_uic       = y_uic-Def_Btn_Height-bdy;
  318.             pos_txt_zer = [xleft, y_uic+d_txt/2, wid1, Def_Txt_Height];
  319.             x_uic       = pos_txt_zer(1)+pos_txt_zer(3)+wx;
  320.             pos_edi_zer = [x_uic, y_uic, wid2, Def_Btn_Height];
  321.             x_uic       = pos_edi_zer(1)+pos_edi_zer(3)+wx;
  322.             pos_txt_zpc = [x_uic, y_uic+d_txt/2, wid3, Def_Txt_Height];
  323.         end
  325.         w_uic = w_rem-2*bdx;
  326.         x_uic = xleft+bdx;
  327.         h_uic = depOfMachine(Def_Btn_Height);        
  328.         y_uic = y_uic-2*bdy-h_uic;
  329.         pos_tog_thr = [x_uic, y_uic, w_uic, h_uic];
  330.         w_uic = w_fra/2-bdx;
  331.         x_uic = pos_fra(1);
  332.         h_uic = (3*Def_Btn_Height)/2;
  333.         y_uic = pos_fra(2)-h_uic-Def_Btn_Height/2;     
  334.         pos_pus_est = [x_uic, y_uic, w_uic, h_uic];
  335.         x_uic = x_uic+w_uic+2*bdx;
  336.         pos_tog_res = [x_uic, y_uic, w_uic, h_uic];
  338.         % Create UIC.
  339.         %------------
  340.         comProp = {...
  341.            'Parent',fig,    ...
  342.            'Unit',fig_units ...
  343.            'Visible','Off'  ...
  344.            };
  345.         comTxtProp = {comProp{:}, ...
  346.            'Style','Text',...
  347.            'HorizontalAlignment','center', ...
  348.            'Backgroundcolor',bkColor ...
  349.            };
  351.         fra_utl = uicontrol(comProp{:}, ...
  352.                             'Style','frame', ...
  353.                             'Position',pos_fra, ...
  354.                             'Backgroundcolor',bkColor, ...
  355.                             'Tag',tag_fra_tool ...
  356.                             );
  358.         txt_top = uicontrol(comProp{:}, ...
  359.                             'Style','Text', ...
  360.                             'Position',pos_txt_top,   ...
  361.                             'String',str_txt_top,     ...
  362.                             'Backgroundcolor',bkColor ...
  363.                             );
  365.         cba = [mfilename '(''update_methName'',' str_numfig ');'];
  366.         pop_met = uicontrol(comProp{:}, ...
  367.                             'Style','Popup',...
  368.                             'Position',pos_pop_met,...
  369.                             'Enable',statusINI, ...
  370.                             'String',str_pop_met,...
  371.                             'HorizontalAlignment','center',...
  372.                             'Userdata',1,...
  373.                             'Callback',cba ...
  374.                             );
  376.         switch toolOPT
  377.           case {'deno','esti','esti_REG','esti_DEN'}
  378.             rad_sof = uicontrol(comProp{:}, ...
  379.                                 'Style','RadioButton',...
  380.                                 'Position',pos_rad_sof,...
  381.                                 'Enable',statusINI, ...
  382.                                 'HorizontalAlignment','center',...
  383.                                 'String',str_rad_sof,...
  384.                                 'Value',1,'Userdata',1 ...
  385.                                 );
  387.             rad_har = uicontrol(comProp{:}, ...
  388.                                 'Style','RadioButton',...
  389.                                 'Position',pos_rad_har,...
  390.                                 'Enable',statusINI, ...
  391.                                 'HorizontalAlignment','center',...
  392.                                 'String',str_rad_har,...
  393.                                 'Value',0,'Userdata',0 ...
  394.                                 );
  395.             cba = [mfilename '(''update_thrType'',' str_numfig ');'];
  396.             set(rad_sof,'Callback',cba);
  397.             set(rad_har,'Callback',cba);
  399.             txt_noi = uicontrol(comProp{:}, ...
  400.                                 'Style','Text',...
  401.                                 'Position',pos_txt_noi,...
  402.                                 'Backgroundcolor',bkColor,...
  403.                                 'String',str_txt_noi...
  404.                                 );
  406.             cba = [mfilename '(''update_by_Caller'',' str_numfig ');'];
  407.             pop_noi = uicontrol(comProp{:}, ...
  408.                                 'Style','Popup',...
  409.                                 'Position',pos_pop_noi,...
  410.                                 'Enable',statusINI, ...
  411.                                 'String',str_pop_noi,...
  412.                                 'HorizontalAlignment','center',...
  413.                                 'Userdata',1,...
  414.                                 'Callback',cba ...
  415.                                 );
  416.           case {'comp'}
  417.         end
  418.         txt_BMS = uicontrol(comProp{:}, ...
  419.                             'Style','Text',...
  420.                             'Position',pos_txt_BMS,...
  421.                             'Backgroundcolor',bkColor,...
  422.                             'String',str_txt_BMS...
  423.                             );
  425.         cba = [mfilename '(''update_by_Caller'',' str_numfig ');'];
  426.         sli_BMS = uicontrol(comProp{:}, ...
  427.                             'Style','Slider',...
  428.                             'Position',pos_sli_BMS,...
  429.                             'Enable',statusINI,    ...
  430.                             'Min',1+sqrt(eps),     ...
  431.                             'Max',5-sqrt(eps),     ...
  432.                             'Value',1.5,           ...
  433.                             'Backgroundcolor',bkColor, ...
  434.                             'Callback',cba         ...
  435.                             );
  436.         txt_tit = zeros(4,1);
  437.         for k=1:4
  438.             txt_tit(k) = uicontrol(...
  439.                                    comTxtProp{:}, ...
  440.                                    'Position',pos_lev_tit(k,:), ...
  441.                                    'String',deblank(str_txt_tit(k,:))...
  442.                                    );
  443.         end
  444.         xbtn0 = xleft;
  445.         ybtn0 = pos_lev_tit(1,2)-Def_Btn_Height;
  446.         xbtn  = xbtn0;
  447.         ybtn  = ybtn0;
  448.         if ud.ydir==1
  449.             index = [1:levmaxMAX];
  450.         else
  451.             index = [levmaxMAX:-1:1];
  452.             ybtn  = ybtn0+(levmaxMAX-levmax)*dy_lev;
  453.         end
  454.         for j=1:length(index)
  455.             i = index(j);
  456.             max_lev = 1;
  457.             val_lev = 0.5;
  458.             pos_lev = [xbtn ybtn+d_txt/2 w_lev(1) Def_Txt_Height];
  459.             str_lev = sprintf('%.0f',i);
  460.             txt_lev = uicontrol(...
  461.                          comTxtProp{:},     ...
  462.                          'Position',pos_lev,...
  463.                          'String',str_lev,  ...
  464.                          'Userdata',i       ...
  465.                          );
  467.             xbtn      = xbtn+w_lev(1)+wx;
  468.             pos_lev = [xbtn ybtn w_lev(2) Def_Btn_Height];
  469.             pop_lev = uicontrol(...
  470.                          comProp{:}, ...
  471.                          'Style','popupmenu',...
  472.                          'Enable',enaVal,  ...
  473.                          'Position',pos_lev,...
  474.                          'String','1',...
  475.                          'Backgroundcolor',bkColor, ...
  476.                          'Userdata',i ...
  477.                          );
  479.             xbtn    = xbtn+w_lev(2)+wx;
  480.             pos_lev = [xbtn, ybtn+sli_dy, w_lev(3), sli_hi];
  481.             sli_lev = uicontrol(...
  482.                          comProp{:},         ...
  483.                          'Style','Slider',   ...
  484.                          'Enable',enaVal,    ...
  485.                          'Position',pos_lev, ...
  486.                          'Min',0,            ...
  487.                          'Max',max_lev,      ...
  488.                          'Value',val_lev,    ...
  489.                          'Userdata',i        ...
  490.                          );
  492.             xbtn    = xbtn+w_lev(3)+wx;
  493.             pos_lev = [xbtn ybtn w_lev(4) Def_Btn_Height];
  494.             str_val = sprintfLOC('%7.3f',val_lev);
  495.             edi_lev = uicontrol(...
  496.                          comProp{:},         ...
  497.                          'Style','Edit',     ...
  498.                          'Enable',enaVal,    ...
  499.                          'Position',pos_lev, ...
  500.                          'String',str_val,   ...
  501.                          'HorizontalAlignment','center',...
  502.                          'Backgroundcolor',Def_EdiBkColor,...
  503.                          'Userdata',i          ...
  504.                          );
  506.             beg_cba = [mfilename '(''update_by_UIC'',' str_numfig ',' str_lev];
  507.             cba_pop = [beg_cba ',''pop'');'];
  508.             cba_sli = [beg_cba ',''sli'');'];
  509.             cba_edi = [beg_cba ',''edi'');'];
  510.             set(pop_lev,'Callback',cba_pop);
  511.             set(sli_lev,'Callback',cba_sli);
  512.             set(edi_lev,'Callback',cba_edi);
  514.             h_CMD_LVL(:,i) = [txt_lev;pop_lev;sli_lev;edi_lev];
  515.             xbtn = xbtn0;
  516.             ybtn = ybtn-dy_lev;
  517.         end
  518.         switch  toolOPT
  519.           case {'deno','esti','esti_REG','esti_DEN'}
  520.           case {'comp'}
  521.             comEdiProp = {comProp{:}, ...
  522.                 'Style','Edit',...
  523.                 'String','',...
  524.                 'Enable','Inactive', ...
  525.                 'Backgroundcolor',bkColor,...
  526.                 'HorizontalAlignment','center'...
  527.                 };
  528.             txt_nor = uicontrol(comTxtProp{:}, ...
  529.                                 'Position',pos_txt_nor,...
  530.                                 'HorizontalAlignment','left',...
  531.                                 'String',str_txt_nor...
  532.                                 );
  534.             cba_nor = [mfilename '(''updateTHR'',' str_numfig ',''nor'');'];
  535.             edi_nor = uicontrol(comEdiProp{:}, ...
  536.                                 'Position',pos_edi_nor,...
  537.                                 'Callback',cba_nor ...
  538.                                 );
  540.             txt_npc = uicontrol(comTxtProp{:}, ...
  541.                                 'Position',pos_txt_npc,...
  542.                                 'String','%'...
  543.                                 );
  545.             txt_zer = uicontrol(comTxtProp{:}, ...
  546.                                 'Position',pos_txt_zer,...
  547.                                 'HorizontalAlignment','left',...
  548.                                 'String',str_txt_zer...
  549.                                 );
  551.             cba_zer = [mfilename '(''updateTHR'',' str_numfig ',''zer'');'];
  552.             edi_zer = uicontrol(comEdiProp{:}, ...
  553.                                 'Position',pos_edi_zer,...
  554.                                 'Callback',cba_zer ...
  555.                                 );
  557.             txt_zpc = uicontrol(comTxtProp{:}, ...
  558.                                 'Position',pos_txt_zpc,...
  559.                                 'String','%'...
  560.                                 );
  561.         end
  562.         cba = [mfilename '(''init_SetThr'',' str_numfig ');'];
  563.         tip = 'Interval dependent threshold Settings';
  564.         tog_thr = uicontrol(...
  565.                             comProp{:},             ...
  566.                             'Style','Togglebutton', ...
  567.                             'Position',pos_tog_thr, ...
  568.                             'String',str_tog_thr,   ...
  569.                             'Enable',enaVal,        ...
  570.                             'Callback',cba,         ...
  571.                             'tooltip',tip,          ...
  572.                             'Interruptible','Off'   ...
  573.                             );
  574.         cba = [mfilename '(''residuals'',' str_numfig ');'];
  575.         tip = 'More on Residuals';
  576.         tog_res = uicontrol(...
  577.                             comProp{:},             ...
  578.                             'Style','Togglebutton', ...
  579.                             'Position',pos_tog_res, ...
  580.                             'String',str_tog_res,   ...
  581.                             'Enable','off',         ...
  582.                             'Callback',cba,         ...
  583.                             'tooltip',tip,...
  584.                             'Interruptible','Off'   ...
  585.                             );
  587.         if isequal(toolOPT,'esti_DEN')
  588.            set(tog_res,'Visible','Off')
  589.            pos_pus_est(1) =  pos_pus_est(1)+pos_pus_est(3)/2;
  590.         end
  592.         cba_pus_est = [mfilename '(''' estOPT ''',' str_numfig ');'];
  593.         pus_est = uicontrol(comProp{:},             ...
  594.                             'Style','Pushbutton',   ...
  595.                             'Position',pos_pus_est, ...
  596.                             'String',xlate(str_pus_est),   ...
  597.                             'Enable',enaVal,        ...
  598.                             'callback',cba_pus_est  ...
  599.                             );
  601. % Add Context Sensitive Help (CSHelp).
  602. %-------------------------------------
  603.         switch  toolOPT
  604.           case {'deno','esti','esti_REG','esti_DEN'}
  605. hdl_DENO_SOFTHARD   = [rad_sof,rad_har];
  606. hdl_DENO_NOISSTRUCT = [txt_noi,pop_noi];
  607. wfighelp('add_ContextMenu',fig,hdl_DENO_SOFTHARD,'DENO_SOFTHARD');
  608. wfighelp('add_ContextMenu',fig,hdl_DENO_NOISSTRUCT,'DENO_NOISSTRUCT');
  610.           case {'comp'}
  611. hdl_COMP_SCORES = [...
  612. txt_nor,edi_nor,txt_npc,txt_zer,edi_zer,txt_zpc ...
  613. ];
  614. wfighelp('add_ContextMenu',fig,hdl_COMP_SCORES,'COMP_SCORES');
  615.         end
  616.         hdl_COMP_DENO_STRA = [...
  617. fra_utl,txt_top,pop_met,txt_BMS,sli_BMS,  ...
  618. ];
  619. hdl_IDTS_GUI = tog_thr;
  620. wfighelp('add_ContextMenu',fig,hdl_COMP_DENO_STRA,'COMP_DENO_STRA');
  621. wfighelp('add_ContextMenu',fig,hdl_IDTS_GUI,'IDTS_GUI');
  622. %-------------------------------------
  624. % Store handles.
  625. %--------------
  626. switch  toolOPT
  627.           case {'deno','esti','esti_REG'}
  628.             ud.handlesUIC = ...
  629.                 [fra_utl;txt_top;pop_met;...
  630.                 rad_sof;rad_har;txt_noi;pop_noi; ...
  631.                 txt_BMS;sli_BMS;txt_tit(1:4); ...
  632.                 NaN;NaN;NaN;NaN;NaN;NaN; ...
  633.                 tog_thr;tog_res;pus_est];
  635.           case {'esti_DEN'}
  636.             ud.handlesUIC = ...
  637.                 [fra_utl;txt_top;pop_met;...
  638.                 rad_sof;rad_har;txt_noi;pop_noi; ...
  639.                 txt_BMS;sli_BMS;txt_tit(1:4); ...
  640.                 NaN;NaN;NaN;NaN;NaN;NaN; ...
  641.                 tog_thr;tog_res;pus_est];
  643.           case {'comp'}
  644.             ud.handlesUIC = ...
  645.                 [fra_utl;txt_top;pop_met; ...
  646.                 NaN;NaN;NaN;NaN; ...
  647.                 txt_BMS;sli_BMS;txt_tit(1:4);...
  648.                 txt_nor;edi_nor;txt_npc;txt_zer;edi_zer;txt_zpc; ...
  649.                 tog_thr;tog_res;pus_est];
  650.         end
  651.         ud.h_CMD_LVL = h_CMD_LVL;
  652.         ud.h_GRA_LVL = h_GRA_LVL;
  653.         set(fra_utl,'Userdata',ud);
  654.         varargout{1} = utthrw1d('set',fig,'position',{levmin,levmax});
  656.     case {'denoise','compress','estimate'}
  657.         feval(calledFUN,option,fig);
  659.     case 'get_LVL_par'
  660.         numMeth = get(pop_met,'value');
  661.         meth    = wthrmeth(toolOPT,'shortnames',numMeth);
  662.         switch  toolOPT
  663.           case {'deno','esti','esti_REG'}
  664.              valType = get(rad_sof,'value');
  665.              if valType==1 , sorh = 's'; else , sorh = 'h'; end
  666.              switch numMeth
  667.                case {1,2,3,4}
  668.                  valNoise = get(pop_noi,'value');
  669.                  switch valNoise
  670.                    case 1 , scal = 'sln';
  671.                    case 2 , scal = 'one';
  672.                    case 3 , scal = 'mln';
  673.                  end
  674.                case {5,6,7}  , scal = get(sli_BMS,'value');
  675.              end
  676.              varargout = {numMeth,meth,scal,sorh};
  678.           case {'esti_DEN'}
  679.              valType = get(rad_sof,'value');
  680.              if valType==1 , sorh = 's'; else , sorh = 'h'; end
  681.              switch numMeth
  682.                case {1,2,3} , scal = [];
  683.                case {4}  , scal = get(sli_BMS,'value');
  684.              end
  685.              varargout = {numMeth,meth,scal,sorh};
  687.           case {'comp'}
  688.              sorh = 'h'; 
  689.              switch numMeth
  690.                case {1,2,3} , scal = get(sli_BMS,'value');
  691.                case {4,5}   , scal = NaN;
  692.              end
  693.              varargout = {numMeth,meth,scal,sorh};
  694.         end
  696.     case 'update_methName'
  697.         numMeth = get(pop_met,'value');
  698.         switch toolOPT
  699.           case {'deno','esti','esti_REG'}
  700.             HDL_1 = [txt_BMS;sli_BMS];
  701.             HDL_2 = [txt_noi;pop_noi];
  702.             switch numMeth
  703.               case {1,2,3,4}
  704.                 invis  = HDL_1;   vis = HDL_2;
  705.                 radDef = rad_sof; radNoDef = rad_har;
  706.               case {5,6,7}
  707.                 invis  = HDL_2;   vis = HDL_1;
  708.                 radDef = rad_har; radNoDef = rad_sof;
  709.             end
  710.             set(sli_BMS,'value',3)
  711.             set(invis,'Visible','off')
  712.             set(vis,'Visible','on')
  713.             set(radDef,'Value',1,'Userdata',1);
  714.             set(radNoDef,'Value',0,'Userdata',0);
  716.           case {'esti_DEN'}
  717.             switch numMeth
  718.               case {1,2,3}
  719.                 set([txt_noi;pop_noi;txt_BMS;sli_BMS],'Visible','off')
  720.                 radDef = rad_sof; radNoDef = rad_har;
  721.               case {4}
  722.                 set([txt_noi;pop_noi],'Visible','off')
  723.                 set([txt_BMS;sli_BMS],'Visible','on')
  724.                 radDef = rad_sof; radNoDef = rad_har;
  725.             end
  726.             set(sli_BMS,'value',3)
  727.             set(radDef,'Value',1,'Userdata',1);
  728.             set(radNoDef,'Value',0,'Userdata',0);
  730.           case {'comp'}
  731.             HDL_1 = [txt_BMS;sli_BMS];
  732.             HDL_2 = [];
  733.             switch numMeth
  734.               case {1,2,3} , invis = HDL_2; vis = HDL_1;                
  735.               case {4,5}   , invis = HDL_1; vis = HDL_2;
  736.             end
  737.             set(sli_BMS,'value',1.5);
  738.             set(invis,'Visible','off')
  739.             set(vis,'Visible','on')
  740.         end
  741.         utthrw1d('update_by_Caller',fig)
  743.     case 'update_thrType'
  744.         rad = gcbo;
  745.         old = get(rad,'userdata');
  746.         if isequal(old,1) , set(rad,'Value',1); return; end
  747.         if isequal(rad,rad_sof)
  748.            type = 's'; other = rad_har;
  749.         else
  750.            type = 'h'; other = rad_sof;           
  751.         end
  752.         set(other,'Value',0,'Userdata',0);
  753.         set(rad,'Value',1,'Userdata',1);
  754.         feval(calledFUN,'clear_GRAPHICS',fig);
  756.     case 'update_hdlDEN'
  757.         thrStruct = varargin{1};
  758.         NB_lev  = size(thrStruct,1);
  759.         pop_int = h_CMD_LVL(2,1:NB_lev);
  760.         sli_lev = h_CMD_LVL(3,1:NB_lev);
  761.         edi_lev = h_CMD_LVL(4,1:NB_lev);
  762.         for k = 1:NB_lev
  763.             thrParams = thrStruct(k).thrParams;
  764.             if ~isempty(thrParams)
  765.                 linvalHdl = thrStruct(k).hdlLines;
  766.                 if ishandle(linvalHdl)
  767.                    NB_int  = size(thrParams,1);
  768.                    act_int = 1;
  769.                    set(pop_int(k),'String',int2str([1:NB_int]'),'Value',act_int);
  770.                    thr     = thrParams(act_int,3);
  771.                    oldMax  = get(sli_lev(k),'Max');
  772.                    newMax  = max([thr oldMax]);
  773.                    set(sli_lev(k),'Max',newMax,'Value',thr);
  774.                    set(edi_lev(k),'string',sprintfLOC('%7.3f',thr));
  776.                    axeHdl  = get(linvalHdl,'Parent');
  777.                    lv      = findobj(axeHdl,'tag',tag_lineV);
  778.                    if ~isempty(lv) , delete(lv); end
  779.                    lu      = findobj(axeHdl,'tag',tag_lineH_up);
  780.                    ld      = findobj(axeHdl,'tag',tag_lineH_down);
  781.                    [xHOR,yHOR] = getxy(thrParams);
  782.                    yHOR    = abs(yHOR);
  783.                    set(lu,'Xdata',xHOR,'Ydata',yHOR);
  784.                    ind     = find(yHOR<sqrt(eps));
  785.                    tmp     = yHOR;
  786.                    if ~isempty(ind) , tmp(ind) = NaN; end
  787.                    set(ld,'Xdata',xHOR,'Ydata',-tmp);
  788.                    yVMin   = 2*max(abs(get(linvalHdl,'Ydata')));
  789.                    cbthrw1d('plotLV',[fig ; lu ; ld],[xHOR ; yHOR],yVMin);
  790.                 end
  791.             end
  792.         end
  793.         utthrw1d('set',fig,'thrstruct',thrStruct);
  795.     case 'update_by_UIC'
  796.         level    = varargin{1};
  797.         type_hdl = varargin{2};
  798.         pop      = h_CMD_LVL(2,level);
  799.         sli      = h_CMD_LVL(3,level);
  800.         edi      = h_CMD_LVL(4,level);
  801.         lHu      = h_GRA_LVL(2,level);
  802.         lHd      = h_GRA_LVL(3,level);
  804.         % Updating threshold.
  805.         %---------------------
  806.         ok = 0;
  807.         switch type_hdl
  808.             case 'pop'
  810.             case 'sli'
  811.               thresh = get(sli,'value');
  812.               set(edi,'string',sprintfLOC('%7.3f',thresh));
  814.             case 'edi'
  815.               valstr = get(edi,'string');
  816.               [thresh,count,err] = sscanf(valstr,'%f');
  817.               if (count~=1) || ~isempty(err)
  818.                   thresh = get(sli,'Value');
  819.                   set(edi,'string',sprintfLOC('%7.3f',thresh));
  820.                   return
  821.               else
  822.                   mi = get(sli,'Min');
  823.                   ma = get(sli,'Max');
  824.                   if     thresh<mi , thresh = mi;
  825.                   elseif thresh>ma , thresh = ma;
  826.                   end
  827.                   set(sli,'Value',thresh);
  828.                   set(edi,'string',sprintfLOC('%7.3f',thresh));
  829.               end
  830.         end
  831.         num_int = get(pop,'Value');
  832.         ydata   = get(lHu,'Ydata');
  833.         i_beg   = 3*(num_int-1)+1;
  834.         switch type_hdl
  835.             case 'pop'
  836.               thresh  = ydata(i_beg);
  837.               set(sli,'Value',thresh);
  838.               set(edi,'string',sprintfLOC('%7.3f',thresh));
  840.             case {'sli','edi'}
  841.               i_end   = i_beg+1;
  842.               ydata([i_beg i_end]) = [thresh thresh];
  843.               set(lHu,'Ydata', ydata);
  844.               if thresh<sqrt(eps)
  845.                   ydata([i_beg i_end]) = [NaN NaN];
  846.               end
  847.               set(lHd,'Ydata',-ydata);
  848.               cbthrw1d('upd_thrStruct',fig,lHu);
  849.               ok = 1;
  850.         end
  851.         feval(calledFUN,'clear_GRAPHICS',fig);
  852.         if ok , utthrw1d('show_LVL_perfos',fig); end
  854.     case 'update_by_Caller'
  855.         feval(calledFUN,'update_LVL_meth',fig);
  857.     case 'update_LVL_meth'
  858.         % called by : calledFUN('update_LVL_meth', ...)
  859.         %----------------------------------------------
  860.         valTHR = varargin{1};
  861.         NB_lev = length(valTHR);
  862.         maxTHR = utthrw1d('get',fig,'maxTHR');
  863.         maxTHR = maxTHR(ud.levmin:ud.levmin+NB_lev-1);
  864.         valTHR = min(valTHR,maxTHR);
  865.         thrStruct = wmemtool('rmb',fig,n_memblocTHR,ind_thr_struct);
  866.         xmin = thrStruct(ud.levmin).thrParams(1,1);
  867.         xmax = thrStruct(ud.levmin).thrParams(end,2);
  868.         nb_MAX = min([ud.levmax,ud.levmin+NB_lev-1]);
  869.         for k=ud.levmin:nb_MAX
  870.             thrStruct(k).thrParams = [xmin xmax valTHR(k)];
  871.         end
  872.         utthrw1d('update_hdlDEN',fig,thrStruct);
  874.     case 'show_LVL_perfos'
  875.         if isequal(toolOPT,'comp')
  876.             feval(calledFUN,'show_LVL_perfos',fig);
  877.         end
  879.     case 'init_SetThr'
  880.         utthrset('init');
  882.     case 'return_SetThr'
  883.         if isempty(varargin{1}) , return; end        
  884.         feval(calledFUN,'clear_GRAPHICS',fig);
  885.         utthrw1d('update_hdlDEN',fig,varargin{:})
  886.         utthrw1d('show_LVL_perfos',fig);
  888.     case 'status'
  889.       if length(varargin)>0
  890.          toolStatus = varargin{1};
  891.          ud.status  = toolStatus;
  892.          set(fra,'Userdata',ud);
  893.          set([pop_met;rad_sof;rad_har;pop_noi;sli_BMS],'Enable',toolStatus)
  894.          if isequal(lower(toolStatus),'off')
  895.              utthrw1d('enable',fig,'off')
  896.          end       
  897.       end
  898.       varargout{1} = toolStatus;
  899.               
  900.     case 'enable'
  901.         enaVal = varargin{1};
  902.         if length(varargin)>1
  903.             levs = varargin{2};
  904.         else
  905.             levs = [1:size(h_GRA_LVL,2)];
  906.         end
  907.         uic = h_CMD_LVL(2:4,:);
  908.         set([uic(:);tog_thr;tog_res;pus_est],'Enable','off');
  909.         if isequal(lower(enaVal),'on')
  910.             uic = h_CMD_LVL(2:4,levs);
  911.             set([uic(:);tog_thr;pus_est],'Enable',enaVal);
  912.         end
  914.     case 'enable_tog_res'
  915.         enaVal = varargin{1};
  916.         set(tog_res,'Enable',enaVal);
  918.     case 'visible'
  919.         visVal     = lower(varargin{1});
  920.         ud.visible = visVal;
  921.         handlesAXE = h_GRA_LVL(1,ud.levmin:ud.levmax);
  922.         handlesAXE = findobj(handlesAXE(ishandle(handlesAXE(:))));
  923.         if isequal(visVal,'on')
  924.             h_CMD_LVL = h_CMD_LVL(1:4,ud.levmin:ud.levmax);
  925.             numMeth = get(pop_met,'value');
  926.             switch toolOPT
  927.               case {'deno','esti','esti_REG'}
  928.                 switch numMeth
  929.                   case {1,2,3,4} , invis = [txt_BMS;sli_BMS];
  930.                   case {5,6,7}   , invis = [txt_noi;pop_noi];
  931.                 end
  933.               case {'esti_DEN'}
  934.                 switch numMeth
  935.                   case {1,2,3} , invis = [txt_noi;pop_noi;txt_BMS;sli_BMS];
  936.                   case {4}     , invis = [txt_BMS;sli_BMS];
  937.                 end
  939.               case {'comp'}
  940.                 switch numMeth
  941.                    case {1,2,3} , invis = [];
  942.                    case {4,5}   , invis = [txt_BMS;sli_BMS];
  943.                 end 
  944.             end
  945.             handlesUIC = setdiff(handlesUIC,invis);
  946.         end
  947.         handles = [h_CMD_LVL(:);handlesAXE(:);handlesUIC(:)];
  948.         set(handles(ishandle(handles)),'visible',visVal);
  949.         if isequal(toolOPT,'esti_DEN')
  950.            set(tog_res,'Visible','Off')
  951.         end
  952.         
  953.     case 'set'
  954.         nbarg = length(varargin);
  955.         if nbarg<1 , return; end
  956.         for k = 1:2:nbarg
  957.            argType = lower(varargin{k});
  958.            argVal  = varargin{k+1};
  959.            switch argType
  960.              case 'position'
  961.                [levmin,levmax] = deal(argVal{:});
  962.                nblevs = levmax-levmin+1;
  963.                if ud.ydir==1
  964.                    dnum_lev = (levmin-ud.levmin);
  965.                else
  966.                    dnum_lev = (ud.levmax-levmax);
  967.                end
  968.                ud.levmin = levmin;
  969.                ud.levmax = levmax;
  970.                set(fra,'Userdata',ud);
  971.                old_units = get(fig,'units');
  972.                tmpHandles = [h_CMD_LVL(:);handlesUIC(:)];
  973.                tmpHandles = tmpHandles(ishandle(tmpHandles));
  974.                set(tmpHandles,'visible','off');
  975.                set([fig;tmpHandles],'units','pixels');
  976.                [pos_fra,dy_lev,y_thr,y_res,y_est] = ...
  977.                        getPosFraThr(fra,nblevs,toolOPT);
  978.                set(fra,'Position',pos_fra);
  979.                ytrans = dnum_lev*dy_lev;
  980.                for j=1:size(h_CMD_LVL,2)
  981.                    for k = 1:4
  982.                        p = get(h_CMD_LVL(k,j),'Position');
  983.                        set(h_CMD_LVL(k,j),'Position',[p(1),p(2)+ytrans,p(3:4)]);
  984.                    end
  985.                end
  986.                p = get(tog_thr,'Position');
  987.                set(tog_thr,'Position',[p(1),y_thr,p(3:4)]);
  988.                ytrans = y_thr-p(2);
  989.                p = get(tog_res,'Position');
  990.                set(tog_res,'Position',[p(1),y_res,p(3:4)]);
  991.                p = get(pus_est,'Position');
  992.                set(pus_est,'Position',[p(1),y_est,p(3:4)]);
  993.                switch toolOPT
  994.                  case {'comp'}
  995.                    tmpHDL = [txt_nor;edi_nor;txt_npc;txt_zer;edi_zer;txt_zpc];
  996.                    for k = 1:length(tmpHDL)
  997.                        p = get(tmpHDL(k),'Position');
  998.                        set(tmpHDL(k),'Position',[p(1),p(2)+ytrans,p(3:4)]);
  999.                    end
  1000.                end
  1001.                set([fig;tmpHandles],'units',old_units);
  1002.                utthrw1d('visible',fig,ud.visible);
  1003.                if nargout>0
  1004.                    varargout{1} = [pos_fra(1) , y_est , pos_fra([3 4])];
  1005.                end
  1007.              case 'thrstruct'
  1008.                if iscell(argVal)
  1009.                    [xmin,xmax,thr,hdl_lines] = deal(argVal{:});
  1010.                    first = ud.levmin;
  1011.                    last  = ud.levmax;
  1012.                    argVal = struct(...
  1013.                                 'thrParams',cell(ud.levmaxMAX,1), ...
  1014.                                 'hdlLines',cell(ud.levmaxMAX,1)   ...
  1015.                                 );
  1016.                    for k=first:last
  1017.                        argVal(k).thrParams = [xmin xmax thr(k)];
  1018.                        argVal(k).hdlLines  = hdl_lines(k);
  1019.                    end
  1020.                  end
  1021.                  wmemtool('wmb',fig,n_memblocTHR,ind_thr_struct,argVal);
  1023.              case 'intdepthr'
  1024.                  wmemtool('wmb',fig,n_memblocTHR,ind_int_thr,argVal);
  1026.              case 'axes'
  1027.                h_GRA_LVL(1,1:length(argVal)) = argVal;
  1028.                ud.h_GRA_LVL = h_GRA_LVL;
  1029.                set(fra,'Userdata',ud);
  1031.              case 'handleori' , ud.handleORI = argVal; set(fra,'Userdata',ud);
  1032.              case 'handlethr' , ud.handleTHR = argVal; set(fra,'Userdata',ud);
  1033.              case 'handleres' , ud.handleRES = argVal; set(fra,'Userdata',ud);
  1034.              case 'perfos'
  1035.                switch toolOPT
  1036.                  case {'comp'}
  1037.                    set(edi_nor,'string',sprintf('%5.2f',argVal{1}));
  1038.                    set(edi_zer,'string',sprintf('%5.2f',argVal{2}));
  1039.                end
  1040.            end
  1041.         end
  1043.     case 'get'
  1044.         nbarg = length(varargin);
  1045.         if nbarg<1 , return; end
  1046.         for k = 1:nbarg
  1047.            outType = lower(varargin{k});
  1048.            switch outType
  1049.              case 'position'
  1050.                pos_fra = get(fra,'Position');
  1051.                pos_est = get(pus_est,'Position');
  1052.                varargout{k} = [pos_fra(1) , pos_est(2) , pos_fra([3 4])];
  1054.              case 'thrstruct'
  1055.                varargout{k} = wmemtool('rmb',fig,n_memblocTHR,ind_thr_struct);
  1057.              case 'intdepthr'
  1058.                varargout{k} = wmemtool('rmb',fig,n_memblocTHR,ind_int_thr);
  1060.              case 'maxthr'
  1061.                val = get(h_CMD_LVL(3,:),'Max');
  1062.                val = cat(2,val{:});
  1063.                varargout{k} = val;
  1065.              case 'valthr'
  1066.                val = get(h_CMD_LVL(3,:),'Value');
  1067.                val = cat(2,val{ud.levmin:ud.levmax});
  1068.                varargout{k} = val;
  1070.              case 'tog_thr'   , varargout{k} = tog_thr;
  1071.              case 'tog_res'   , varargout{k} = tog_res;
  1072.              case {'pus_den','pus_est','pus_com'} , varargout{k} = pus_est;
  1073.              case 'handleori' , varargout{k} = ud.handleORI;
  1074.              case 'handlethr' , varargout{k} = ud.handleTHR;
  1075.              case 'handleres' , varargout{k} = ud.handleRES;
  1076.            end
  1077.         end
  1079.     case 'clean_thr'
  1080.         set(h_CMD_LVL(2,:),'String','1','Value',1);       % pop_int;
  1081.         set(h_CMD_LVL(3,:),'Min',0,'Max',1,'Value',0.5);  % sli_lev;
  1082.         set(h_CMD_LVL(4,:),'String','');                  % edi_lev
  1084.     case 'plot_dec'
  1085.         lev = varargin{1};
  1086.         [max_det,val_lev,xmin,xmax,k] = deal(varargin{2}{:});
  1087.         set(h_CMD_LVL(3,lev),'Min',0,'Max',max_det,'Value',val_lev);
  1088.         set(h_CMD_LVL(4,lev),'String',sprintfLOC('%7.3f',val_lev));
  1089.         xHOR    = [xmin xmax];
  1090.         yHOR    = [val_lev val_lev];
  1091.         [lu,ld] = cbthrw1d('plotLH',[h_CMD_LVL(2:4,lev);h_GRA_LVL(1,lev)], ...
  1092.                            xHOR,yHOR,k,max_det);
  1093.         h_GRA_LVL(2:3,lev) = [lu;ld];
  1094.         ud.h_GRA_LVL = h_GRA_LVL;
  1095.         set(fra,'Userdata',ud);
  1096.  
  1097.     case 'residuals'
  1098.         [handleORI,handleTHR,handleRES] = ...
  1099.             utthrw1d('get',fig,'handleORI','handleTHR','handleRES');
  1100.         wmoreres('create',fig,tog_res,handleRES,handleORI,handleTHR,'blocPAR');
  1102.     case 'den_M1'
  1103.         cfs = varargin{1};
  1104.         len = varargin{2};
  1105.         thrStruct = utthrw1d('get',fig,'thrstruct');
  1106.         switch toolOPT
  1107.           case {'deno','esti','esti_REG','esti_DEN'}
  1108.             v = get(rad_sof,'Value');
  1109.             if v==1 , sorh = 's'; else , sorh = 'h'; end
  1110.           case {'comp'} , sorh = 'h';
  1111.         end
  1112.         NB_lev = size(thrStruct,1);
  1113.         for k = 1:NB_lev
  1114.             thr_par = thrStruct(k).thrParams;
  1115.             if ~isempty(thr_par)
  1116.                 NB_int = size(thr_par,1);
  1117.                 x      = [thr_par(:,1) ; thr_par(NB_int,2)];
  1118.                 x      = round(x);
  1119.                 x(x<1) = 1;
  1120.                 x(x>len) = len;
  1121.                 thr = thr_par(:,3);
  1122.                 for j = 1:NB_int
  1123.                     if j==1 , d_beg = 0; else , d_beg = 1; end
  1124.                     j_beg = x(j)+d_beg;
  1125.                     j_end = x(j+1);
  1126.                     j_ind = [j_beg:j_end];
  1127.                     cfs(k,j_ind) = wthresh(cfs(k,j_ind),sorh,thr(j));
  1128.                 end
  1129.             end
  1130.         end
  1131.         varargout{1} = cfs;
  1133.     case 'den_M2'
  1134.         coefs = varargin{1};
  1135.         longs = varargin{2};
  1136.         len   = longs(end);
  1137.         first = cumsum(longs)+1;
  1138.         first = first(end-2:-1:1);
  1139.         tmp   = longs(end-1:-1:2);
  1140.         last  = first+tmp-1;
  1141.         longs = longs(end-1:-1:2);
  1142.         thrStruct = utthrw1d('get',fig,'thrstruct');
  1143.         switch toolOPT
  1144.           case {'deno','esti','esti_REG','esti_DEN'}
  1145.             v = get(rad_sof,'Value');
  1146.             if v==1 , sorh = 's'; else , sorh = 'h'; end
  1148.           case {'comp'} , sorh = 'h';
  1149.         end
  1150.         NB_lev  = size(thrStruct,1);        
  1151.         for k = 1:NB_lev
  1152.             thr_par = thrStruct(k).thrParams;
  1153.             if ~isempty(thr_par)
  1154.                 cfs = coefs(first(k):last(k));
  1155.                 nbCFS = longs(k);            
  1156.                 NB_int = size(thr_par,1);
  1157.                 x   = [thr_par(:,1) ; thr_par(NB_int,2)];
  1158.                 alf = (nbCFS-1)/(x(end)-x(1));
  1159.                 bet = 1 - alf*x(1);
  1160.                 x   = round(alf*x+bet);
  1161.                 x(x<1) = 1;
  1162.                 x(x>nbCFS) = nbCFS;
  1163.                 thr = thr_par(:,3);
  1164.                 for j = 1:NB_int
  1165.                     if j==1 , d_beg = 0; else , d_beg = 1; end
  1166.                     j_beg = x(j)+d_beg;
  1167.                     j_end = x(j+1);
  1168.                     j_ind = [j_beg:j_end];
  1169.                     cfs(j_ind) = wthresh(cfs(j_ind),sorh,thr(j));
  1170.                 end
  1171.                 coefs(first(k):last(k)) = cfs;
  1172.             end
  1173.         end
  1174.         varargout{1} = coefs;
  1176.     case 'demo'
  1177.     % SPECIAL for DEMOS
  1178.     %------------------
  1179.     [tool,parDemo] = deal(varargin{:});
  1180.     switch tool
  1181.        case {'sw1d','wdre'}
  1182.          if ischar(parDemo{1})
  1183.              shortnames = wthrmeth(toolOPT,'shortnames');
  1184.              ind = strmatch(deblank(parDemo{1}),shortnames);
  1185.              if isempty(ind) , return; end
  1186.              set(pop_met,'Value',ind)
  1187.              utthrw1d('update_methName',fig)
  1189.          elseif isnumeric(parDemo{1})
  1190.              utthrset('demo',tog_thr,parDemo{1})
  1191.          end
  1192.     end
  1195.     otherwise
  1196.         errargt(mfilename,'Unknown Option','msg');
  1197.         error('*');
  1198. end
  1200. %=============================================================================%
  1201. % INTERNAL FUNCTIONS
  1202. %=============================================================================%
  1203. %-----------------------------------------------------------------------------%
  1204. function varargout = wthrmeth(toolOPT,varargin)
  1206. switch toolOPT
  1207.   case {'deno','esti','esti_REG'}
  1208.     thrMethods = {...
  1209.         'Fixed form threshold',  'sqtwolog',    1; ...
  1210.         'Rigorous SURE',         'rigrsure',    2; ...
  1211.         'Heuristic SURE',        'heursure',    3; ...
  1212.         'Minimax',               'minimaxi',    4; ...
  1213.         'Penalize high',         'penalhi',     5; ...
  1214.         'Penalize medium',       'penalme',     6; ...
  1215.         'Penalize low',          'penallo',     7  ...
  1216.         };
  1218.   case {'esti_DEN'}
  1219.     thrMethods = {...
  1220.         'Global threshold',      'globalth',    1; ...
  1221.         'By level threshold 1',  'bylevth1',    2; ...
  1222.         'By level threshold 2',  'bylevth2',    3; ...
  1223.         'By level threshold 3',  'bylevsth',    4  ...
  1224.         };
  1226.   case 'comp'
  1227.     thrMethods = {...
  1228.         'Scarce high',           'scarcehi',     1; ...
  1229.         'Scarce medium',         'scarceme',     2; ...
  1230.         'Scarce low',            'scarcelo',     3; ...
  1231.         'Balance sparsity-norm', 'bal_sn',       4; ...
  1232.         'Remove near 0',         'rem_n0',       5  ...
  1233.         };
  1234. end
  1235. nbin = length(varargin);
  1236. if nbin==0 , varargout{1} = thrMethods; return; end
  1238. option = varargin{1};
  1239. switch option
  1240.   case 'names'
  1241.      out = strvcat(1,thrMethods{:,1});
  1242.      varargout{1} = out(2:end,:);
  1243.      if nbin==2
  1244.          num = varargin{2};
  1245.          varargout{1} = deblank(varargout{1}(num,:));
  1246.      end
  1248.   case 'shortnames'
  1249.      out = strvcat(1,thrMethods{:,2});
  1250.      varargout{1} = out(2:end,:);
  1251.      if nbin==2
  1252.          num = varargin{2};
  1253.          varargout{1} = deblank(varargout{1}(num,:));
  1254.      end
  1256.   case 'nums'
  1257.      varargout{1} = cat(1,thrMethods{:,3});
  1258. end
  1259. %-----------------------------------------------------------------------------%
  1260. function [pos_fra,dy_lev,y_thr,y_res,y_est] = getPosFraThr(fra,nblevs,toolOPT)
  1262. Def_Btn_Height = mextglob('get','Def_Btn_Height');
  1263. [bdy,d_lev] = get_DLEV(nblevs,toolOPT);
  1264. pos_fra = get(fra,'Position');
  1265. top_fra = pos_fra(2)+pos_fra(4);
  1266. if ~isequal(toolOPT,'esti_DEN') , multiY = 1; else , multiY = 1; end
  1267. switch  toolOPT
  1268.   case {'deno','esti'} , NB_Height = 6;
  1269.   case {'esti_REG'}    , NB_Height = 6;
  1270.   case {'esti_DEN'}    , NB_Height = 4+2*multiY;  
  1271.   case {'comp'}        , NB_Height = 7;
  1272. end
  1273. h_ini   = (NB_Height-1)*bdy+NB_Height*Def_Btn_Height;
  1274. h_fra   = h_ini+ nblevs*(Def_Btn_Height+d_lev)+ ...
  1275.           depOfMachine(Def_Btn_Height);
  1276. pos_fra(2) = top_fra-h_fra;
  1277. pos_fra(4) = h_fra;
  1278. dy_lev = d_lev+Def_Btn_Height;
  1279. y_thr  = pos_fra(2)+bdy;
  1280. y_est = pos_fra(2)-(3*Def_Btn_Height)/2-Def_Btn_Height/2;
  1281. y_res = y_est;
  1282. %-----------------------------------------------------------------------------%
  1283. function [x,y] = getxy(arg)
  1284.         
  1285. NB_int  = size(arg,1);
  1286. if NB_int>1
  1287.     x = [arg(:,1:2) NaN*ones(NB_int,1)]';
  1288.     x = x(:)';
  1289.     l = 3*NB_int-1;
  1290.     x = x(1:l);
  1291.     y = [arg(:,[3 3]) NaN*ones(NB_int,1)]';
  1292.     y = y(:)';
  1293.     y = y(1:l);
  1294. else
  1295.     x = arg(1,1:2); y = arg(1,[3 3]);
  1296. end
  1297. %-----------------------------------------------------------------------------%
  1298. function varargout = depOfMachine(varargin)
  1300. btn_height = varargin{1};
  1301. scrSize = get(0,'ScreenSize');
  1302. if scrSize(4)<600
  1303.     height = btn_height;
  1304. elseif scrSize(4)<800
  1305.     height = 3*btn_height/2;
  1306. else
  1307.     height = 3*btn_height/2;
  1308. end
  1309. varargout = {height};
  1310. %-----------------------------------------------------------------------------%
  1311. function s = sprintfLOC(frm,x);
  1313. s = sprintf(frm,x);
  1314. s(s==' ') = [];
  1315. %-----------------------------------------------------------------------------%
  1316. function [bdy,d_lev] = get_DLEV(nblevs,toolOPT)
  1318. if isequal(toolOPT,'comp') && (nblevs>10)
  1319.     bdy = 2; d_lev = 0;    
  1320. else
  1321.     bdy = 4; d_lev = 2;
  1322. end
  1323. %=============================================================================%