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

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

波变换

开发平台：

Matlab

wdretool.m：源码内容
  1. function varargout = wdretool(option,varargin)
  2. %WDRETOOL Wavelet Density and Regression tool.
  3. %   VARARGOUT = WDRETOOL(OPTION,VARARGIN)
  5. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 05-Dec-96.
  6. %   Last Revision: 05-Feb-2004.
  7. %   Copyright 1995-2004 The MathWorks, Inc.
  8. %   $Revision: 1.12.4.2 $  $Date: 2004/03/15 22:42:38 $
  10. % Test inputs.
  11. %-------------
  12. if nargin==0 , option = 'createREG'; end
  13. [option,winAttrb] = utguidiv('ini',option,varargin{:});
  15. % Memory Blocks of stored values.
  16. %================================
  17. % MB1.
  18. %-----
  19. n_membloc1   = 'MB_1';
  20. ind_xdata    = 1;
  21. ind_ydata    = 2;
  22. ind_xbounds  = 3;
  23. ind_filename = 4;
  24. ind_pathname = 5;
  25. ind_sig_name = 6;
  26. nb2_stored   = 6;
  28. % MB2.
  29. %-----
  30. n_membloc2   = 'MB_2';
  31. ind_status   = 1;
  32. ind_lin_den  = 2;
  33. ind_gra_area = 3;
  34. nb1_stored   = 3;
  36. % MB3.
  37. %-----
  38. n_membloc3 = 'MB_3';
  39. ind_coefs  = 1;
  40. ind_longs  = 2;
  41. ind_sig    = 3;
  42. nb4_stored = 3;
  44. % Default values.
  45. %----------------
  46. def_MinBIN = 64;
  47. def_DefBIN = 256;
  48. default_wave = 'sym4';
  49. NB_max_lev = 8;
  50. NB_def_lev = 5;
  51. yLevelDir  = -1;
  53. % Tag property.
  54. %--------------
  55. tag_ori = 'Sig';
  56. tag_dat = 'Proc_Data';
  57. tag_app = 'App';
  58. tag_est = 'Est';
  60. switch option
  61.   case {'createDEN','createREG'}
  62.   case 'close' ,
  63.   otherwise
  64.     win_tool = varargin{1};
  65.     handles  = wfigmngr('getValue',win_tool,['WDRE_handles']);
  66.     hdl_UIC  = handles.hdl_UIC;
  67.     hdl_AXE  = handles.hdl_AXE;
  68.     txt_hdl  = handles.hdl_TXT;
  69.     hdl_MEN  = handles.hdl_MEN;
  70.     men_sav  = hdl_MEN(end);
  71.     dummy    = struct2cell(hdl_UIC);
  72.     [txt_bin,sli_bin,edi_bin,pus_dec,chk_den] = deal(dummy{:});
  73.     axe_hdl  = struct2cell(hdl_AXE);
  74.     [axe_L_1,axe_R_1,axe_L_2,axe_R_2,axe_cfs,axe_det,axe_app] =  ...
  75.            deal(axe_hdl{:});
  76.     axe_hdl  = cat(2,axe_hdl{:});
  77.     colors   = wfigmngr('getValue',win_tool,['WDRE_colors']);
  78.     toolATTR = wfigmngr('getValue',win_tool,['WDRE_toolATTR']);
  79.     switch option
  80.       case {'load','demo'}
  81.       otherwise , toolMode = toolATTR.toolMode;
  82.     end
  83. end
  85. switch option
  86.     case {'createDEN','createREG'}
  88.         % Parameters initialization.
  89.         %---------------------------
  90.         indic_vis_lev = getLevels(NB_max_lev,yLevelDir);
  92.         % Get Globals.
  93.         %-------------
  94.         [Def_Txt_Height,Def_Btn_Height,X_Spacing,Y_Spacing, ...
  95.          sliYProp,Def_EdiBkColor, Def_FraBkColor] = ...
  96.             mextglob('get',...
  97.                 'Def_Txt_Height','Def_Btn_Height','X_Spacing','Y_Spacing', ...
  98.                 'Sli_YProp','Def_EdiBkColor','Def_FraBkColor' ...
  99.                 );
  101.         % Window initialization.
  102.         %-----------------------
  103.         switch option
  104.           case 'createREG'
  105.             win_title = 'Regression Estimation 1-D';
  106.             estiNAME  = 'esti_REG';
  107. figExtMode = 'ExtFig_Tool_1';
  109.           case 'createDEN'
  110.             win_title = 'Density Estimation 1-D';
  111.             estiNAME  = 'esti_DEN';
  112. figExtMode = 'ExtFig_Tool_3';
  113.         end
  114.         [win_tool,pos_win,win_units,str_numwin,...
  115.             frame0,pos_frame0,Pos_Graphic_Area,pus_close] = ...
  116.                wfigmngr('create',win_title,winAttrb,figExtMode,mfilename,1,1,0);
  117.         if nargout>0 , varargout{1} = win_tool; end
  119.         % Menu construction for current figure.
  120.         %--------------------------------------
  121.         m_files = wfigmngr('getmenus',win_tool);
  122.         switch option
  123.           case 'createREG'
  125. % Add Help for Tool.
  126. %------------------
  127. wfighelp('addHelpTool',win_tool, ...
  128. 'One-Dimensional Estimation - &Random Design','REGR_GUI');
  129. wfighelp('addHelpTool',win_tool, ...
  130. 'One-Dimensional Estimation - &Fixed Design','REGF_GUI');
  132. % Add Help Item.
  133. %----------------
  134. wfighelp('addHelpItem',win_tool,'Regression Estimation','REG_EST');
  135. wfighelp('addHelpItem',win_tool,'Available Methods','COMP_DENO_METHODS');
  136. wfighelp('addHelpItem',win_tool,'Variance Adaptive Thresholding','VARTHR');
  137. wfighelp('addHelpItem',win_tool,'Loading and Saving','REG_LOADSAVE');
  139. m_load = wfigmngr('getmenus',win_tool,'load');  
  140.             lab = 'Data for &Fixed Design Regression';
  141.             men_fix = uimenu(m_load,...
  142.                          'Label',lab,'Position',1, ...
  143.                          'Callback',               ...
  144.                          [mfilename '(''load'','   ...
  145.                              str_numwin ',''fixreg'');'] ...
  146.                          );
  147.             lab = 'Data for &Stochastic Design Regression';
  148.             men_sto = uimenu(m_load,...
  149.                         'Label',lab,'Position',2, ...
  150.                         'Callback',               ...
  151.                         [mfilename '(''load'','   ...
  152.                             str_numwin ',''storeg'');'] ...
  153.                         );
  154.             men_den = [];
  155.             men_sav = uimenu(m_files,...
  156.                          'Label','&Save Estimated Function',...
  157.                          'Position',2,            ...
  158.                          'Enable','Off',          ...
  159.                          'Callback',              ...
  160.                          [mfilename '(''save'','  ...
  161.                              str_numwin ',''fun'');']...
  162.                          );
  164.             m_demo = uimenu(m_files,...
  165.                         'Label','&Example Analysis','Position',3,'Separator','Off');
  166.  
  167.             numDEM = 1;
  168.             labDEM = 'Fixed Design';
  169.             m_demo_1 = uimenu(m_demo, ...
  170.                         'Label',labDEM,'Position',numDEM,'Separator','Off');
  172.             demoSET = {...
  173.               'fixreg' , 'Example I'     , 'ex1nfix'  , 'db2'  , 5 ; ...
  174.               'fixreg' , 'Example II'    , 'ex2nfix'  , 'sym4' , 5 ; ...
  175.               'fixreg' , 'Example II'    , 'ex3nfix'  , 'db3'  , 5 ; ...
  176.               'fixreg' , 'Noisy blocks'  , 'noisbloc' , 'haar' , 5 ; ...
  177.               'fixreg' , 'Noisy Doppler' , 'noisdopp' , 'db5'  , 5 ; ...
  178.               'fixreg' , 'Noisy bumps'   , 'noisbump' , 'db5'  , 5   ...
  179.               };
  180.             setDEMOS(m_demo_1,mfilename,str_numwin,demoSET,0)
  182.             numDEM = 2;
  183.             labDEM = ['Fixed Design - Interval Dependent Noise Variance'];
  184.             m_demo_2 = uimenu(m_demo, ...
  185.                         'Label',labDEM,'Position',numDEM,'Separator','Off');
  187.             demoSET = {...
  188.               'fixreg' , 'Noisy blocks (I)  - 3 intervals' , 'nblocr1' , 'sym4', 5 , '{3}' ; ...
  189.               'fixreg' , 'Noisy blocks (II) - 3 intervals' , 'nblocr2' , 'sym4', 5 , '{3}' ; ...
  190.               'fixreg' , 'Noisy Doppler (I) - 3 intervals' , 'ndoppr1' , 'sym4', 5 , '{3}' ; ...
  191.               'fixreg' , 'Noisy bumps (I)   - 3 intervals' , 'nbumpr1' , 'sym4', 5 , '{3}' ; ...
  192.               'fixreg' , 'Noisy bumps (II)  - 2 intervals' , 'nbumpr2' , 'sym4', 5 , '{2}' ; ...
  193.               'fixreg' , 'Noisy bumps (III) - 4 intervals' , 'nbumpr3' , 'sym4', 5 , '{4}' ; ...
  194.               'fixreg' , 'Elec. consumption - 3 intervals' , 'nelec'   , 'sym4', 5 , '{3}'   ...
  195.               };
  196.             setDEMOS(m_demo_2,mfilename,str_numwin,demoSET,0)
  198.             numDEM = 3;
  199.             labDEM = 'Stochastic Design';
  200.             m_demo_3 = uimenu(m_demo,...
  201.                         'Label',labDEM,'Position',numDEM,'Separator','On');
  202.             demoSET = {...
  203.               'storeg' , 'Example I'     , 'ex1nsto'  , 'sym4' ,5 ; ...
  204.               'storeg' , 'Example II'    , 'ex2nsto'  , 'haar' ,5 ; ...
  205.               'storeg' , 'Example III'   , 'ex3nsto'  , 'db6'  ,5 ; ...
  206.               'storeg' , 'Noisy Doppler' , 'noisdopp' , 'db5'  ,5 ; ...
  207.               'storeg' , 'Noisy bumps'   , 'noisbump' , 'db5'  ,5   ...
  208.               };
  209.             setDEMOS(m_demo_3,mfilename,str_numwin,demoSET,0)
  210.  
  211.             numDEM = 4;
  212.             labDEM = ['Stochastic Design '...
  213.                 ' - Interval Dependent Noise Variance'];
  214.             m_demo_4 = uimenu(m_demo,...
  215.                         'Label',labDEM,'Position',numDEM,'Separator','Off');
  216.             demoSET = {...
  217.               'storeg' , 'Noisy blocks (I)  - 3 intervals' , 'snblocr1' , 'sym4', 5 , '{3}' ; ...
  218.               'storeg' , 'Noisy blocks (II) - 3 intervals' , 'snblocr2' , 'sym4', 5 , '{3}' ; ...
  219.               'storeg' , 'Noisy Doppler (I) - 3 intervals' , 'sndoppr1' , 'sym4', 5 , '{3}' ; ...
  220.               'storeg' , 'Noisy bumps (I)   - 3 intervals' , 'snbumpr1' , 'sym4', 5 , '{3}' ; ...
  221.               'storeg' , 'Noisy bumps (II)  - 2 intervals' , 'snbumpr2' , 'sym4', 5 , '{2}' ; ...
  222.               'storeg' , 'Noisy bumps (III) - 4 intervals' , 'snbumpr3' , 'sym4', 5 , '{4}' ; ...
  223.               'storeg' , 'Elec. consumption - 3 intervals' , 'snelec'   , 'sym4', 5 , '{3}'   ...
  224.               };
  225.             setDEMOS(m_demo_4,mfilename,str_numwin,demoSET,0)
  228.           case 'createDEN'
  230. % Add Help for Tool.
  231. %------------------
  232. wfighelp('addHelpTool',win_tool, ...
  233. '&One-Dimensional Estimation','EDEN_GUI');
  235. % Add Help Item.
  236. %----------------
  237. wfighelp('addHelpItem',win_tool,'Available Methods','COMP_DENO_METHODS');
  238. wfighelp('addHelpItem',win_tool,'Variance Adaptive Thresholding','VARTHR');
  239. wfighelp('addHelpItem',win_tool,'Loading and Saving','EDEN_LOADSAVE');
  240. wfighelp('addHelpItem',win_tool,'Density Estimation','DENS_EST');
  241.   
  242.             lab = '&Load data for Density Estimate';
  243.             men_fix = [];
  244.             men_sto = [];
  245.             men_den = uimenu(m_files,...
  246.                          'Label',lab,'Position',1, ...
  247.                          'Callback',               ...
  248.                          [mfilename '(''load'','   ...
  249.                              str_numwin ',''denest'');'] ...
  250.                          );
  252.             men_sav = uimenu(m_files,...
  253.                          'Label','&Save Density', ...
  254.                          'Position',2,           ...
  255.                          'Enable','Off',         ...
  256.                          'Callback',             ...
  257.                          [mfilename '(''save'',' ...
  258.                              str_numwin ',''den'');'] ...
  259.                          );
  260.             numDEM = 3;
  261.             demoSET = {...
  262.               'denest' , 'ex1cusp1' , 'ex1cusp1' , 'sym4'  , 5 ; ...
  263.               'denest' , 'ex2cusp1' , 'ex2cusp1' , 'sym6'  , 5 ; ...
  264.               'denest' , 'ex1cusp2' , 'ex1cusp2' , 'sym4'  , 5 ; ...
  265.               'denest' , 'ex2cusp2' , 'ex2cusp2' , 'coif1' , 5 ; ...
  266.               'denest' , 'ex1gauss' , 'ex1gauss' , 'sym3'  , 5 ; ...
  267.               'denest' , 'ex2gauss' , 'ex2gauss' , 'sym4'  , 5   ...
  268.               };
  269.             m_demo = uimenu(m_files,'Label','&Example Analysis ','Position',3);
  270.             setDEMOS(m_demo,mfilename,str_numwin,demoSET,0);
  271.         end
  273.         % Begin waiting.
  274.         %---------------
  275.         wwaiting('msg',win_tool,'Wait ... initialization');
  277.         % General parameters initialization.
  278.         %-----------------------------------
  279.         dx = X_Spacing;
  280.         dy = Y_Spacing;
  281.         d_txt  = (Def_Btn_Height-Def_Txt_Height);
  282.         deltaY = (Def_Btn_Height+dy);
  283.         sli_hi = Def_Btn_Height*sliYProp;
  284.         sli_dy = 0.5*Def_Btn_Height*(1-sliYProp);
  285.         minimum_bin = 16;
  286.         maximun_bin = 1024;
  287.         default_bin = 256;
  289.         % String property of objects.
  290.         %----------------------------
  291.         str_pus_dec = 'Decompose';
  292.         str_chk_den = 'Overlay Estimated Function';
  293.         str_txt_bin = 'Nb bins';
  294.         str_edi_bin = sprintf('%.0f',default_bin);
  296.         % Command part of the window.
  297.         %============================
  298.         comFigProp = {'Parent',win_tool,'Unit',win_units};
  300.         % Data, Wavelet and Level parameters.
  301.         %------------------------------------
  302.         xlocINI = pos_frame0([1 3]);
  303.         ytopINI = pos_win(4)-dy;
  304.         toolPos = utanapar('create',win_tool, ...
  305.                     'xloc',xlocINI,'top',ytopINI,...
  306.                     'enable','off',      ...
  307.                     'wtype','dwt',       ...
  308.                     'deflev',NB_def_lev, ...
  309.                     'maxlev',NB_max_lev  ...
  310.                     );
  311.         pop_lev = utanapar('handles',win_tool,'lev');
  313.         % Callbacks.
  314.         %-----------
  315.         cba_sli_bin = [mfilename '(''upd_bin'',' str_numwin ',''sli'');'];
  316.         cba_edi_bin = [mfilename '(''upd_bin'',' str_numwin ',''edi'');'];
  317.         cba_pus_dec = [mfilename '(''decompose'',' str_numwin ');'];
  318.         cba_chk_den = [mfilename '(''show_lin_den'',' str_numwin ');'];
  320.         % Bins settings.
  321.         %---------------
  322.         w_bas = toolPos(3)/48;        
  323.         h_uic = Def_Btn_Height;
  325.         w_uic = 11*w_bas;
  326.         x_uic = toolPos(1);
  327.         y_uic = toolPos(2)-h_uic-2*dy;
  328.         pos_txt_bin = [x_uic, y_uic+d_txt/2, w_uic, Def_Txt_Height];
  330.         x_uic = x_uic+w_uic;
  331.         w_uic = 27*w_bas;
  332.         pos_sli_bin = [x_uic, y_uic+sli_dy, w_uic, sli_hi];
  334.         x_uic = x_uic+w_uic+4;
  335.         w_uic = 10*w_bas-4;
  336.         pos_edi_bin = [x_uic, y_uic, w_uic, h_uic];
  337.        
  338.         txt_bin = uicontrol(comFigProp{:},...
  339.                             'Style','Text',...
  340.                             'Position',pos_txt_bin,...
  341.                             'HorizontalAlignment','left',...
  342.                             'Backgroundcolor',Def_FraBkColor,...
  343.                             'String',str_txt_bin...                            
  344.                             );
  345.  
  346.         sli_bin = uicontrol(comFigProp{:},...
  347.                             'Style','Slider',...
  348.                             'Position',pos_sli_bin,...
  349.                             'Min',minimum_bin,...
  350.                             'Max',maximun_bin,...
  351.                             'Value',default_bin, ...
  352.                             'Enable','off', ...
  353.                             'Callback',cba_sli_bin ...
  354.                             );
  356.         edi_bin = uicontrol(comFigProp{:},...
  357.                             'Style','Edit',...
  358.                             'Backgroundcolor',Def_EdiBkColor,...
  359.                             'Position',pos_edi_bin,...
  360.                             'String',str_edi_bin,...
  361.                             'Enable','off',...
  362.                             'Callback',cba_edi_bin ...
  363.                             );
  365.         % Decompose pushbutton.
  366.         %----------------------
  367.         h_uic = 3*Def_Btn_Height/2;
  368.         y_uic = y_uic-h_uic-2*dy;
  369.         w_uic = pos_frame0(3)/2;
  370.         x_uic = pos_frame0(1)+(pos_frame0(3)-w_uic)/2;
  371.         pos_pus_dec = [x_uic, y_uic, w_uic, h_uic];
  372.         pus_dec = uicontrol(comFigProp{:},          ...
  373.                             'Style','Pushbutton',   ...
  374.                             'Position',pos_pus_dec, ...
  375.                             'String',xlate(str_pus_dec),   ...
  376.                             'Enable','off',         ...
  377.                             'Interruptible','On',   ...
  378.                             'Callback',cba_pus_dec  ...
  379.                             );
  381.         % Thresholding tool.
  382.         %-------------------
  383.         ytopTHR = pos_pus_dec(2)-2*dy;
  384.         toolPos = utthrw1d('create',win_tool, ...
  385.                  'xloc',xlocINI,'top',ytopTHR,...
  386.                  'ydir',yLevelDir,       ...
  387.                  'levmax',NB_def_lev,    ...
  388.                  'levmaxMAX',NB_max_lev, ...
  389.                  'status','Off',         ...
  390.                  'toolOPT',estiNAME      ...
  391.                  );
  394.         % Estimated Line(s) Check.
  395.         %-------------------------
  396.         w_uic = (3*pos_frame0(3))/4;
  397.         x_uic = pos_frame0(1)+(pos_frame0(3)-w_uic)/2;
  398.         h_uic = Def_Btn_Height;
  399.         y_uic = toolPos(2)-Def_Btn_Height/2-h_uic;
  400.         pos_chk_den = [x_uic, y_uic, w_uic, h_uic];
  401.         if isequal(option,'createDEN') , vis = 'Off'; else , vis = 'On'; end  
  402.         chk_den = uicontrol(comFigProp{:},          ...
  403.                             'Style','checkbox',     ...
  404.                             'Visible',vis,          ...
  405.                             'Position',pos_chk_den, ...
  406.                             'String',str_chk_den,   ...
  407.                             'Enable','off',         ...
  408.                             'Callback',cba_chk_den  ...
  409.                             );
  411.         % Callbacks update.
  412.         %------------------
  413.         hdl_den = utthrw1d('handles',win_tool);
  414.         utanapar('set_cba_num',win_tool,[m_files;hdl_den(:)]);
  415.         pop_lev = utanapar('handles',win_tool,'lev');
  416.         tmp     = num2mstr(pop_lev );
  417.         cba_pop_lev = [mfilename '(''upd_lev'',' str_numwin ',' tmp ');'];
  418.         set(pop_lev,'Callback',cba_pop_lev);
  420.         % General graphical parameters initialization.
  421.         %--------------------------------------------
  422.         txtLRProp = {'off','bold',14};
  424.         % Axes construction parameters.
  425.         %------------------------------
  426.         NB_lev    = NB_max_lev;    % dummy
  427.         w_gra_rem = Pos_Graphic_Area(3);
  428.         h_gra_rem = Pos_Graphic_Area(4);
  429.         ecx_left  = 0.08*pos_win(3);
  430.         ecx_med   = 0.07*pos_win(3);
  431.         ecx_right = 0.06*pos_win(3);
  432.         w_axe     = (w_gra_rem-ecx_left-ecx_med-ecx_right)/2;
  433.         x_cfs     = ecx_left;
  434.         x_det     = x_cfs+w_axe+ecx_med;
  435.         ecy_up    = 0.06*pos_win(4);
  436.         ecy_mid_1 = 0.07*pos_win(4);
  437.         ecy_mid_2 = ecy_mid_1;
  438.         ecy_mid_3 = ecy_mid_1;
  439.         ecy_det   = (0.04*pos_win(4))/1.4;
  440.         ecy_down  = ecy_up;
  441.         h_min     = h_gra_rem/12;
  442.         h_max     = h_gra_rem/5;
  443.         h_axe_std = (h_min*NB_lev+h_max*(NB_max_lev-NB_lev))/NB_max_lev;
  444.         h_space   = ecy_up+ecy_mid_1+ecy_mid_2+ecy_mid_3+...
  445.                     (NB_lev-1)*ecy_det+ecy_down;
  446.         h_detail  = (h_gra_rem-2*h_axe_std-h_space)/(NB_lev+1);
  447.         y_low_ini = pos_win(4);
  449.         % Building data axes.
  450.         %--------------------
  451.         comAxeProp = {...
  452.           comFigProp{:},'Visible','off','Drawmode','fast','Box','On'};
  453.         y_low_ini = y_low_ini-h_axe_std-ecy_up;
  454.         pos_L_1   = [x_cfs y_low_ini w_axe h_axe_std];
  455.         axe_L_1   = axes(comAxeProp{:},'Position',pos_L_1);
  456.         pos_R_1   = [x_det y_low_ini w_axe h_axe_std];
  457.         axe_R_1   = axes(comAxeProp{:},'Position',pos_R_1);
  458.         y_low_ini = y_low_ini-h_axe_std-ecy_mid_1;
  459.         pos_L_2   = [x_cfs y_low_ini w_axe h_axe_std];
  460.         axe_L_2   = axes(comAxeProp{:},'Position',pos_L_2);
  461.         pos_R_2   = [x_det y_low_ini w_axe h_axe_std];
  462.         axe_R_2   = axes(comAxeProp{:},'Position',pos_R_2);
  463.         y_low_ini = y_low_ini-h_axe_std-ecy_mid_2;
  465.         % Building approximation axes on the right part.
  466.         %-----------------------------------------------
  467.         pos_app = [x_det y_low_ini w_axe h_detail];
  468.         axe_app = axes(comAxeProp{:},'Position',pos_app);
  469.         str_txt = ['a' wnsubstr(abs(NB_max_lev))];
  470.         txt_app = txtinaxe('create',str_txt,axe_app,'r',txtLRProp{:});
  471.         y_low_ini = y_low_ini-h_detail-ecy_mid_3+ecy_det;
  473.         % Building details axes on the left part.
  474.         %----------------------------------------
  475.         comAxePropMore = {...
  476.           comAxeProp{:},'XTicklabelMode','manual','XTickLabel',[]};
  477.         axe_cfs = zeros(1,NB_max_lev);
  478.         txt_cfs = zeros(1,NB_max_lev);
  479.         y_cfs   = y_low_ini;
  480.         pos_cfs = [x_cfs y_cfs w_axe h_detail];
  481.         for j = 1:NB_max_lev
  482.             k = indic_vis_lev(j);
  483.             pos_cfs(2) = pos_cfs(2)-pos_cfs(4)-ecy_det;
  484.             axe_cfs(k) = axes(comAxePropMore{:},'Position',pos_cfs);
  485.             str_txt    = ['d' wnsubstr(k)];
  486.             txt_cfs(k) = txtinaxe('create',str_txt,axe_cfs(k),'l',txtLRProp{:});
  487.             set(txt_cfs(k),'Userdata',k,'tag','');
  488.         end
  489.         utthrw1d('set',win_tool,'axes',axe_cfs);
  491.         % Building details axes on the right part.
  492.         %-----------------------------------------
  493.         axe_det = zeros(1,NB_max_lev);
  494.         txt_det = zeros(1,NB_max_lev);
  495.         y_det   = y_low_ini;
  496.         pos_det = [x_det y_det w_axe h_detail];
  497.         for j = 1:NB_max_lev
  498.             k = indic_vis_lev(j);
  499.             pos_det(2) = pos_det(2)-pos_det(4)-ecy_det;
  500.             axe_det(k) = axes(comAxePropMore{:},'Position',pos_det);
  501.             str_txt    = ['d' wnsubstr(k)];
  502.             txt_det(k) = txtinaxe('create',str_txt,axe_det(k),'r',txtLRProp{:});
  503.         end
  505.         %  Normalization.
  506.         %----------------
  507.         Pos_Graphic_Area = wfigmngr('normalize',win_tool,Pos_Graphic_Area);
  508.         drawnow
  510.         % Set default wavelet.
  511.         %---------------------
  512.         cbanapar('set',win_tool,'wav',default_wave);
  514. % Add Context Sensitive Help (CSHelp).
  515. %-------------------------------------
  516. hdl_BINS  = [txt_bin,sli_bin,edi_bin];
  517. switch option
  518. case 'createDEN' , helpName = 'EDEN_BINS';
  519. case 'createREG' , helpName = 'REG_BINS';
  520. end
  521. wfighelp('add_ContextMenu',win_tool,hdl_BINS,helpName);
  522. %-------------------------------------
  524.         % Memory blocks update.
  525.         %----------------------
  526.         wmemtool('ini',win_tool,n_membloc2,nb1_stored);
  527.         wmemtool('ini',win_tool,n_membloc1,nb2_stored);
  528.         wmemtool('ini',win_tool,n_membloc3,nb4_stored);
  529.         wmemtool('wmb',win_tool,n_membloc2,...
  530.                                 ind_status,0,        ...
  531.                                 ind_lin_den,[NaN,NaN], ...
  532.                                 ind_gra_area,Pos_Graphic_Area ...
  533.                                 );
  534.         fields = {...
  535.           'txt_bin','sli_bin','edi_bin', ...
  536.           'pus_dec','chk_den'            ...
  537.           };
  538.         values = {...
  539.            txt_bin,sli_bin,edi_bin, ...
  540.            pus_dec,chk_den          ...
  541.            };
  542.         hdl_UIC = cell2struct(values,fields,2);
  543.         hdl_MEN = [men_fix ; men_sto ; men_den ; men_sav];
  544.         fields = {...
  545.           'axe_L_1' , 'axe_R_1' , 'axe_L_2' , 'axe_R_2' ,  ...
  546.           'axe_cfs' , 'axe_det' , 'axe_app' ...
  547.           };
  548.         values = {...
  549.           axe_L_1 , axe_R_1 , axe_L_2 , axe_R_2 ,  ...
  550.           axe_cfs , axe_det , axe_app ...
  551.            };
  552.         hdl_AXE = cell2struct(values,fields,2);
  553.         hdl_TXT = [NaN  NaN  NaN  NaN  txt_cfs txt_det txt_app];
  554.         handles = struct(...
  555.             'hdl_MEN',hdl_MEN, ...
  556.             'hdl_UIC',hdl_UIC, ...
  557.             'hdl_AXE',hdl_AXE, ...
  558.             'hdl_TXT',hdl_TXT  ...
  559.             );
  560.         wfigmngr('storeValue',win_tool,['WDRE_handles'],handles);
  562.         %-----------------------------------------------------------------
  563.         % colors = struct( ...
  564.         %   'sigColor',[1 0 0],'denColor',[1 0 0],'appColor',[0 1 1], ...
  565.         %   'detColor',[0 1 0],'cfsColor',[0 1 0],'estColor',[1 1 0]  ...
  566.         %   );
  567.         %-----------------------------------------------------------------
  569.         colors = wtbutils('colors','wdre');
  570.         wfigmngr('storeValue',win_tool,['WDRE_colors'],colors);
  572.         toolATTR = struct('toolMode','','toolState','','level',[],'wname','','NBClasses',[]);
  573.         wfigmngr('storeValue',win_tool,['WDRE_toolATTR'],toolATTR);
  575.         % End waiting.
  576.         %---------------
  577.         wwaiting('off',win_tool);
  579.     case {'load','demo'}
  580.         tool_OPT = varargin{2};
  581.         switch option
  582.           case 'load'
  583.             dialName = ['Load data - '];
  584.             switch tool_OPT
  585.               case 'denest' , dialName = [dialName 'density estimation' ];
  586.               case 'fixreg' , dialName = [dialName 'fixed design'];
  587.               case 'storeg' , dialName = [dialName 'stochastic design'];
  588.             end
  589.             [filename,pathname,ok] = ...
  590.                 utguidiv('test_load',win_tool,'*.mat',dialName);
  591.             if ~ok, return; end
  593.           case 'demo'
  594.             sig_Name = deblank(varargin{3});
  595.             wav_Name = deblank(varargin{4});
  596.             lev_Anal = varargin{5};
  597.             if length(varargin)>5  & ~isempty(varargin{6})
  598.                 parDemo = varargin{6};
  599.             else
  600.                 parDemo = '';
  601.             end
  602.             filename = [sig_Name '.mat'];
  603.             pathname = utguidiv('WTB_DemoPath',filename);
  604.         end
  606.         % Loading file.
  607.         %--------------
  608.         wwaiting('msg',win_tool,'Wait ... loading');
  609.         fullName = fullfile(pathname,filename);
  610.         [fileStruct,err] = wfileinf(fullName);
  611.         if ~err
  612.             try
  613.               load(fullName,'-mat');
  614.             catch
  615.               err = 1;
  616.               msg = sprintf('File %s is not a valid file.', filename);
  617.             end
  618.         end
  619.         if ~err
  620.             % Keep only numeric values.
  621.             %--------------------------
  622.             dum = struct2cell(fileStruct);
  623.             dumClass = dum(4,:);
  624.             idxClass = find(~strcmp(dumClass,'double') & ...
  625.                        ~strcmp(dumClass,'uint8') & ...
  626.                        ~strcmp(dumClass,'sparse'));
  627.             fileStruct(idxClass) = [];
  628.             err = isempty(fileStruct);
  629.         end
  631.         if ~err
  632.             % Keep only one dim values.
  633.             %--------------------------
  634.             dum = struct2cell(fileStruct);
  635.             dumSize = dum(2,:);
  636.             dumSize = cat(1,dumSize{:});
  637.             mini    = min(dumSize,[],2);
  638.             maxi    = max(dumSize,[],2);
  639.             idxSize = find(mini~=1 | maxi<2);
  640.             fileStruct(idxSize) = [];
  641.             err = isempty(fileStruct);
  642.         end
  644.         if ~err
  645.             err = 1;
  646.             switch tool_OPT
  647.               case 'denest'
  648.                 xdata = eval(fileStruct(1).name);
  649.                 ydata = ones(size(xdata));
  650.                 err = 0;
  652.               case {'fixreg','storeg'}
  653.                 % Seek X and Y values.
  654.                 %---------------------
  655.                 dum = struct2cell(fileStruct);
  656.                 dum = lower(dum(1,:));
  657.                 idxVarSET = [1:length(dum)];
  658.                 flagX = 0;
  659.                 flagY = 0;
  660.                 idx_Xdata = find(strcmp(dum,'xdata'));
  661.                 if ~isempty(idx_Xdata)
  662.                    flagX = 1;
  663.                    idxVarSET = setdiff(idxVarSET,idx_Xdata);
  664.                 else
  665.                    idxSET = strmatch('x',dum);
  666.                    if ~isempty(idxSET)
  667.                        flagX = 1;
  668.                        idx_Xdata = idxSET(1);
  669.                        idxVarSET = setdiff(idxVarSET,idx_Xdata);
  670.                    end
  671.                 end
  672.                 idx_Ydata = find(strcmp(dum,'ydata'));
  673.                 if ~isempty(idx_Ydata)
  674.                    flagY = 1;
  675.                    idxVarSET = setdiff(idxVarSET,idx_Ydata);
  676.                 else
  677.                    idxSET = strmatch('y',dum);
  678.                    if ~isempty(idxSET)
  679.                        flagY = 1;
  680.                        idx_Ydata = idxSET(1);
  681.                        idxVarSET = setdiff(idxVarSET,idx_Ydata);
  682.                    end
  683.                 end
  684.                 if ~flagX & ~isempty(idxVarSET)
  685.                     flagX = 1;
  686.                     idx_Xdata = idxVarSET(1);
  687.                     idxVarSET(1) = [];
  688.                 end
  689.                 if ~flagY & ~isempty(idxVarSET)
  690.                     flagY = 1;
  691.                     idx_Ydata = idxVarSET(1);
  692.                     idxVarSET(1) = [];
  693.                 end
  694.                 if ~isempty(idx_Xdata)
  695.                      xdata = eval(fileStruct(idx_Xdata).name);
  696.                 end
  697.                 if ~isempty(idx_Ydata)
  698.                      ydata = eval(fileStruct(idx_Ydata).name);
  699.                 end
  700.                 fileStruct([idx_Xdata,idx_Ydata]) = [];
  701.                 if flagX & ~flagY
  702.                    flagY = 1;
  703.                    ydata = xdata;
  704.                    xdata = [1:length(ydata)];
  706.                 elseif flagY & ~flagX
  707.                    flagX = 1;
  708.                    xdata = [1:length(ydata)];
  709.                 end
  710.                 if flagX & flagY & (length(xdata)==length(ydata))
  711.                     err = 0;
  712.                 end
  713.             end
  714.         end
  715.         if err
  716.             msg = ['File not found or invalid values for variables!'];
  717.         elseif ~isreal(xdata) | ~isreal(ydata)
  718.             msg = strvcat(sprintf('File %s doesn''t contain real data.', filename),' ');
  719.             err = 1;                     
  720.         end
  721.         if err
  722.             wwaiting('off',win_tool);
  723.             errordlg(msg,'Load ERROR','modal');
  724.             return            
  725.         end
  726.         
  727.         [sig_Name,ext] = strtok(filename,'.');
  728.         sig_Size = length(xdata);
  729.         xbounds  = [min(xdata),max(xdata)];
  731.         % Begin waiting.
  732.         %---------------
  733.         wwaiting('msg',win_tool,'Wait ... cleaning');
  735.         % Tool Settings.
  736.         %---------------
  737.         levm   = wmaxlev(sig_Size,'haar');
  738.         levmax = min(levm,NB_max_lev);
  739.         if isequal(option,'demo')
  740.             anaPar = {'wav',wav_Name};
  741.         else
  742.             wav_Name = cbanapar('get',win_tool,'wav');
  743.             lev_Anal = NB_def_lev;
  744.             anaPar = {};
  745.         end
  746.         strlev = int2str([1:levmax]');
  747.         anaPar = {anaPar{:},'n_s',{sig_Name,sig_Size}, ...
  748.                   'lev',{'String',strlev,'Value',lev_Anal}};
  749.         defaultBIN = min([def_DefBIN,fix(sig_Size/4)]);
  750.         minBIN     = min([def_MinBIN,defaultBIN]);            
  751.         maxBIN     = sig_Size;     
  753.         % Store tool settings & analysis parameters.
  754.         %------------------------------------------- 
  755.         toolATTR.toolMode  = tool_OPT;
  756.         toolATTR.toolState = 'ini';
  757.         toolATTR.level     = lev_Anal;
  758.         toolATTR.wname     = wav_Name;
  759.         toolATTR.NBClasses = defaultBIN;
  760.         wfigmngr('storeValue',win_tool,['WDRE_toolATTR'],toolATTR);
  762.         % Store analysis parameters.
  763.         %--------------------------- 
  764.         wmemtool('wmb',win_tool,n_membloc2,ind_status,0);
  765.  
  766.         wmemtool('wmb',win_tool,n_membloc1,   ...
  767.                        ind_filename,filename, ...
  768.                        ind_pathname,pathname, ...
  769.                        ind_sig_name,sig_Name, ...
  770.                        ind_xdata,xdata, ...
  771.                        ind_ydata,ydata, ...
  772.                        ind_xbounds,xbounds ...
  773.                        );
  775.         % Clean , Set analysis & GUI values.
  776.         %-----------------------------------
  777.         dynvtool('stop',win_tool)
  778.         utthrset('stop',win_tool);
  779.         wdretool('clean',win_tool,option);
  780.         cbanapar('set',win_tool,anaPar{:});
  782.         % Set bins.
  783.         %-----------
  784.         set(sli_bin,'Min',minBIN,'Value',defaultBIN,'Max',maxBIN)
  785.         set(edi_bin,'String',int2str(defaultBIN));
  787.         % Setting axes and UIC. 
  788.         %----------------------
  789.         wdretool('position',win_tool,lev_Anal);
  790.         wdretool('enable',win_tool,'ini','on');
  791.         wdretool('set_axes',win_tool);
  793.         % Plot.
  794.         %------
  795.         color   = colors.sigColor;
  796.         linProp = {...
  797.             'Parent',axe_L_1, ...
  798.             'Color',color,    ...
  799.             'LineStyle','none',...
  800.             'Marker','o',     ...
  801.             'Markersize',2,   ...
  802.             'MarkerEdgeColor',color, ...
  803.             'MarkerFaceColor',color, ...
  804.             'tag',tag_ori     ...
  805.             };
  806.         
  807.         switch tool_OPT
  808.           case 'denest'
  809.             xval = linspace(xbounds(1),xbounds(2),sig_Size);
  810.             ll = line(linProp{:},'Xdata',xval,'Ydata',xdata);
  811.             ylim = getylim(xdata);
  812.             set(axe_L_1,'Xlim',xbounds,'Ylim',ylim,'Xtick',[],'XtickLabel',[])
  813.             wtitle('Data','Parent',axe_L_1)
  814.             color  = colors.denColor;
  815.             his    = wgethist(xdata,minBIN);
  816.             wplothis(axe_L_2,his,color);
  817.             strTitle = sprintf('%d bins - Histogram',minBIN);
  818.             wtitle(strTitle,'Parent',axe_L_2)
  819.             wdretool('set_Bins',win_tool);
  820.             
  821.           case {'fixreg','storeg'}
  822.             lin_ori = line(linProp{:},'Xdata',xdata,'Ydata',ydata);
  823.             utthrw1d('set',win_tool,'handleORI',lin_ori);
  824.             ylim = getylim(ydata);
  825.             set(axe_L_1,'Ylim',ylim)          
  826.             wtitle('Data (X,Y)','Parent',axe_L_1)
  827.             if tool_OPT(1)=='s' , wdretool('set_Bins',win_tool); end
  828.         end
  830.         % Plotting Processed data.
  831.         %------------------------
  832.         wdretool('plot_Processed_Data',win_tool);
  834.         % if DEMO, analyze and estimate.
  835.         %-------------------------------
  836.         if isequal(option,'demo')
  837.             wdretool('decompose',win_tool);
  838.             if ~isempty(parDemo)
  839.                  utthrw1d('demo',win_tool,'wdre',parDemo);
  840.             end
  841.             wdretool('estimate',win_tool);
  842.             wdretool('show_lin_den',win_tool,'On')
  843.         end
  844.         cbanapar('enable',win_tool,'On');
  846.         % End waiting.
  847.         %-------------
  848.         wwaiting('off',win_tool);
  850.     case 'save'
  851.         % Testing file.
  852.         %--------------
  853.         dialName = ['Save estimated '];
  854.         switch toolMode
  855.           case {'fixreg','storeg'} , dialName = [dialName 'function'];
  856.           case 'denest' , dialName = [dialName 'denest' ];
  857.         end
  858.         [filename,pathname,ok] = utguidiv('test_save',win_tool, ...
  859.                                      '*.mat',dialName);
  860.         if ~ok, return; end
  862.         % Begin waiting.
  863.         %--------------
  864.         wwaiting('msg',win_tool,'Wait ... saving');
  866.         % Saving file.
  867.         %--------------
  868.         [name,ext] = strtok(filename,'.');
  869.         if isempty(ext) | isequal(ext,'.')
  870.             ext = '.mat'; filename = [name ext];
  871.         end
  873.         % Get de-noising parameters & de-noised signal.
  874.         %-----------------------------------------------
  875.         toolATTR = wfigmngr('getValue',win_tool,'WDRE_toolATTR');
  876.         NB_lev = toolATTR.level;
  877.         wname  = toolATTR.wname;
  878.         [thrStruct,hdl_den] = utthrw1d('get',win_tool,...
  879.                                  'thrstruct','handleTHR');
  880.         thrParams = {thrStruct(1:NB_lev).thrParams};
  881.         xdata = get(hdl_den,'Xdata');
  882.         ydata = get(hdl_den,'Ydata');
  883.         try
  884.           save([pathname filename],'xdata','ydata','thrParams','wname');
  885.         catch
  886.           errargt(mfilename,'Save FAILED !','msg');
  887.         end
  888.    
  889.         % End waiting.
  890.         %-------------
  891.         wwaiting('off',win_tool);
  893.     case 'decompose'
  894.     
  895.         % Compute decomposition and plot.
  896.         %--------------------------------
  897.         wwaiting('msg',win_tool,'Wait ... computing');
  899.         % Clean HDLG.
  900.         %---------------
  901.         utthrw1d('clean_thr',win_tool);
  902.         wdretool('show_lin_den',win_tool,'off')
  904.         % Get analyzis parameters.
  905.         %-------------------------
  906.         [wname,NB_lev] = cbanapar('get',win_tool,'wav','lev');
  907.         NBClasses = round(get(sli_bin,'Value'));
  909.         % Tool Settings.
  910.         %---------------
  911.         toolATTR.toolState = 'dec';
  912.         toolATTR.level     = NB_lev;
  913.         toolATTR.wname     = wname;
  914.         toolATTR.NBClasses = NBClasses;
  915.         wfigmngr('storeValue',win_tool,['WDRE_toolATTR'],toolATTR);
  917.         % Get Handles.
  918.         %-------------
  919.         indic_vis_lev = getLevels(NB_lev,yLevelDir);
  921.         % Clean axes
  922.         %------------
  923.         axes2clean = [axe_R_2,axe_cfs,axe_det,axe_app];
  924.         obj2del = findobj(axes2clean,'type','line');
  925.         delete(obj2del)
  927.         % Get X bounds values.
  928.         %---------------------
  929.         xbounds = wmemtool('rmb',win_tool,n_membloc1,ind_xbounds);        
  930.         xmin = xbounds(1);  xmax = xbounds(2);
  932.         % Compute.
  933.         %---------
  934.         sig_pro = wmemtool('rmb',win_tool,n_membloc3,ind_sig);
  935.         [coefs,longs] = wavedec(sig_pro,NB_lev,wname);
  936.         [det,app] = getFullDecTec(coefs,longs,wname,NB_lev);
  937.         wmemtool('wmb',win_tool,n_membloc3, ...
  938.                        ind_coefs,coefs, ...
  939.                        ind_longs,longs ...
  940.                        );
  942.         % Initializing by level threshold.
  943.         %---------------------------------
  944.         maxTHR = zeros(1,NB_lev);
  945.         for k = 1:NB_lev
  946.             maxTHR(k) = max(abs(detcoef(coefs,longs,k)));
  947.         end
  948.         valTHR = wdretool('compute_LVL_THR',win_tool);
  949.         valTHR = min(valTHR,maxTHR);
  951.         % Plotting Coefficients.
  952.         %-----------------------
  953.         hdl_lines = NaN*ones(NB_max_lev,1);
  954.         color = colors.cfsColor;
  955.         len   = longs(end);
  956.         viewSTEMS = 1;
  957.         if viewSTEMS 
  958.            dummy =  wfilters(wname);
  959.            lf = length(dummy);
  960.         end
  961.         x_cfs = linspace(xbounds(1),xbounds(2),len);
  962.         for j = 1:NB_lev
  963.             k       = indic_vis_lev(j);
  964.             axe_act = axe_cfs(k);
  965.             axes(axe_act);            
  966.             cfs = detcoef(coefs,longs,k);
  967.             tag = ['cfs_' int2str(k)];
  968.             ybounds = [-valTHR(k) , valTHR(k) , -maxTHR(k) , maxTHR(k)];            
  969.             if ~viewSTEMS
  970.                 cfs = cfs(ones(1,2^k),:);
  971.                 cfs = wkeep1(cfs(:)',len);                
  972.                 hdl_lines(k) = plotline(axe_act,x_cfs,cfs,color,tag,ybounds);
  973.             else
  974.                 ld = length(cfs);
  975.                 xd = coefsLOC([1:ld],k,lf,len);                
  976.                 x_tmp = x_cfs(xd);               
  977.                 hh = plotstem(axe_act,x_tmp,cfs,color,1,tag,ybounds);
  978.                 hdl_lines(k) = hh(3);
  979.             end
  980.             utthrw1d('plot_dec',win_tool,k,{maxTHR(k),valTHR(k),xmin,xmax,k});
  981.         end
  983.         i_axe_cfs_max = indic_vis_lev(NB_lev);
  984.         xt = get(axe_app,{'Xtick','XtickLabel'});
  985.         set([axe_det,axe_cfs],'Xtick',[],'XtickLabel',[]);
  986.         set([axe_cfs(i_axe_cfs_max),axe_det(i_axe_cfs_max)], ...
  987.            'Xtick',xt{1},'XtickLabel',xt{2}, ...
  988.            'XtickMode','auto','XTickLabelMode','auto' ...
  989.            );
  990.         set(axe_hdl,'Xlim',[xmin xmax])
  991.         set(txt_hdl(end),'String',['a' wnsubstr(abs(NB_lev))]);
  992.   
  993.         % Dynvtool Attachement.
  994.         %---------------------
  995.         if ~isequal(toolMode,'denest')
  996.            axe_IND = []; axe_CMD = axe_hdl;       
  997.         else
  998.            axe_IND = axe_hdl(1); axe_CMD = axe_hdl(2:end);
  999.         end
  1000.         dynvtool('init',win_tool,axe_IND,axe_CMD,[],[1 0],'','','')
  1002.         % Initialization of Denoising structure.
  1003.         %---------------------------------------
  1004.         utthrw1d('set',win_tool,...
  1005.                        'thrstruct',{xmin,xmax,valTHR,hdl_lines},...
  1006.                        'intdepthr',[]);
  1008.         % Enabling HDLG.
  1009.         %---------------
  1010.         wdretool('enable',win_tool,'dec','on');
  1012.         % Setting prog status.
  1013.         %----------------------
  1014.         wmemtool('wmb',win_tool,n_membloc2,ind_status,1);
  1016.         % End waiting.
  1017.         %-------------
  1018.         wwaiting('off',win_tool);
  1020.     case 'estimate'
  1021.         % Compute decomposition and plot.
  1022.         %--------------------------------
  1023.         wwaiting('msg',win_tool,'Wait ... computing');
  1025.         % Diseable De-noising Tool.
  1026.         %---------------------------
  1027.         utthrw1d('enable',win_tool,'off');
  1029.         % Read Tool Memory Block.
  1030.         %-----------------------
  1031.         [coefs,longs] = wmemtool('rmb',win_tool,n_membloc3,ind_coefs,ind_longs);
  1032.         NB_lev    = toolATTR.level;
  1033.         wname     = toolATTR.wname;
  1034.         NBClasses = toolATTR.NBClasses;
  1035.         indic_vis_lev = getLevels(NB_lev,yLevelDir);
  1037.         % De-noising & Plot de-noised signal.
  1038.         %------------------------------------
  1039.         cden = utthrw1d('den_M2',win_tool,coefs,longs);
  1040.         xbounds = wmemtool('rmb',win_tool,n_membloc1,ind_xbounds);
  1041.         xmin = xbounds(1); xmax = xbounds(2);
  1042.         xval = linspace(xmin,xmax,NBClasses);
  1043.         [det,app,sig_est] = getFullDecTec(cden,longs,wname,NB_lev);
  1045.         % Reset dynvtool.
  1046.         %----------------
  1047.         dynvtool('get',win_tool,0,'force');
  1049.         % Plotting details.
  1050.         %------------------
  1051.         color = colors.detColor;
  1052.         for j = 1:NB_lev
  1053.             k   = indic_vis_lev(j);
  1054.             tag = ['det_' int2str(k)];
  1055.             ll  = plotline(axe_det(k),xval,det(k,:),color,tag);
  1056.         end
  1058.         % Plotting approximation.
  1059.         %------------------------
  1060.         color = colors.appColor;
  1061.         plotline(axe_app,xval,app,color,tag_app);
  1062.         set(txt_hdl(end),'String',['a' wnsubstr(abs(NB_lev))]);
  1064.         % Plotting Estimation.
  1065.         %---------------------
  1066.         color = colors.estColor;
  1067.         if ~isequal(toolMode,'denest')
  1068.             lin_den = wmemtool('rmb',win_tool,n_membloc2,ind_lin_den);
  1069.             set(lin_den,'Ydata',sig_est);
  1070.         else
  1071.             xdata = wmemtool('rmb',win_tool,n_membloc1,ind_xdata);
  1072.         end
  1073.         lin_den = plotline(axe_R_2,xval,sig_est,color,tag_est);
  1074.         i_axe_cfs_max = indic_vis_lev(NB_lev);
  1075.         xt = get(axe_app,{'Xtick','XtickLabel'});
  1076.         set([axe_det,axe_cfs],'Xtick',[],'XtickLabel',[]);
  1077.         set([axe_cfs(i_axe_cfs_max),axe_det(i_axe_cfs_max)], ...
  1078.            'Xtick',xt{1},'XtickLabel',xt{2}, ...
  1079.            'XtickMode','auto','XTickLabelMode','auto' ...
  1080.            );
  1082.         % Dynvtool Attachement.
  1083.         %---------------------
  1084.         if ~isequal(toolMode,'denest')
  1085.            axe_IND = []; axe_CMD = axe_hdl;       
  1086.         else
  1087.            axe_IND = axe_hdl(1); axe_CMD = axe_hdl(2:end);
  1088.         end
  1089.         dynvtool('init',win_tool,axe_IND,axe_CMD,[],[1 0],'','','')
  1091.         % Enabling HDLG.
  1092.         %---------------
  1093.         utthrw1d('set',win_tool,'handleTHR',lin_den);
  1094.         utthrw1d('enable',win_tool,'on');
  1095.         wdretool('enable',win_tool,'den','on');
  1097.         % Storing tool State.
  1098.         %--------------------
  1099.         toolATTR.toolState = 'den';
  1100.         wfigmngr('storeValue',win_tool,['WDRE_toolATTR'],toolATTR);
  1101.         
  1102.         % End waiting.
  1103.         %-------------
  1104.         wwaiting('off',win_tool);
  1106.     case 'show_lin_den'
  1107.         if isequal(toolMode,'denest') , return; end
  1108.         lin_den = wmemtool('rmb',win_tool,n_membloc2,ind_lin_den);
  1109.         ok = ishandle(lin_den(1));
  1110.         if length(varargin)>1
  1111.             vis = lower(varargin{2});
  1112.             if isequal(vis,'on') & ok , val = 1; else , val = 0; end 
  1113.             set(chk_den,'value',val);
  1114.         else
  1115.             vis = deblank(getonoff(get(chk_den,'value')));
  1116.             if ~ok , set(chk_den,'value',0); end
  1117.         end
  1118.         
  1119.         if ok
  1120.             if isequal(vis,'on')
  1121.                 strTitle1 = ['Data (X,Y) and Estimate Y = f(X)'];
  1122.                 strTitle2 = ['Processed Data (X,Y) and Estimate Y = f(X)'];
  1123.             end
  1124.             set(lin_den,'Visible',vis);
  1125.         end
  1126.         wtitle('Data (X,Y)','Parent',axe_L_1)
  1127.         wtitle('Processed Data (X,Y)','Parent',axe_R_1)
  1129.     case 'position'
  1130.         NB_lev  = varargin{2};
  1131.         set(chk_den,'Visible','off');
  1132.         pos_old = utthrw1d('get',win_tool,'position');
  1133.         utthrw1d('set',win_tool,'position',{1,NB_lev})
  1134.         pos_new = utthrw1d('get',win_tool,'position');
  1135.         ytrans  = pos_new(2)-pos_old(2);
  1136.         pos_chk = get(chk_den,'Position');
  1137.         pos_chk(2) = pos_chk(2)+ytrans;
  1138.         switch toolMode
  1139.           case 'denest' , vis_chk = 'off';
  1140.           otherwise ,     vis_chk = 'on';
  1141.         end
  1142.         set(chk_den,'Position',pos_chk,'Visible',vis_chk);
  1144.     case 'upd_lev'
  1145.         pop_lev = varargin{2};
  1146.         lev_New = get(pop_lev,'value');
  1147.         wdretool('position',win_tool,lev_New);
  1148.         wdretool('set_axes',win_tool);
  1149.         lev_Anal = toolATTR.level;
  1150.         flagView = 2;
  1151.         if isequal(lev_New,lev_Anal)
  1152.             state =  toolATTR.toolState;
  1153.             switch state
  1154.                 case 'ini' ,
  1155.                     flagView = 1;
  1156.                     wdretool('enable',win_tool,'ini');
  1157.                 case 'dec' , 
  1158.                     wdretool('enable',win_tool,'dec');
  1159.                 case 'den' ,
  1160.                     val = get(chk_den,'Value');
  1161.                     wdretool('enable',win_tool,'dec','on');
  1162.                     wdretool('enable',win_tool,'den','on');
  1163.                     set(chk_den,'Value',val);
  1164.             end
  1165.         else
  1166.             flagView = 0;
  1167.             wdretool('enable',win_tool,'ini');
  1168.         end
  1169.         if flagView<2
  1170.             wdretool('show_lin_den',win_tool,'off')
  1171.             axe_off = [axe_R_2,axe_app,axe_det];
  1172.             axe_off = axe_off(ishandle(axe_off));
  1173.             lin_Off = findobj(axe_off,'type','line');
  1174.             set(lin_Off,'Visible','off');
  1175.         end
  1177.     case 'upd_bin'
  1178.         typ_upd = varargin{2};
  1179.         switch typ_upd
  1180.           case 'sli'
  1181.             nbOld = round(wstr2num(get(edi_bin,'String')));
  1182.             nbNew = round(get(sli_bin,'Value'));
  1184.           case 'edi'
  1185.             sliVal = get(sli_bin,{'Min','Value','Max'});
  1186.             minNb = sliVal{1};
  1187.             nbOld = sliVal{2};
  1188.             maxNb = sliVal{3};
  1189.             nbNew = round(wstr2num(get(edi_bin,'String')));
  1190.             if isempty(nbNew) | (nbNew<minNb) | (nbNew>maxNb)
  1191.                nbNew = nbOld;
  1192.             end
  1193.         end
  1194.         set(edi_bin,'String',int2str(nbNew));
  1195.         set(sli_bin,'Value',nbNew)
  1196.         if isequal(nbNew,nbOld) , return; end
  1197.         wdretool('clean',win_tool,'bin');
  1199.     case 'compute_LVL_THR'
  1200.         [numMeth,meth,alfa,sorh] = utthrw1d('get_LVL_par',win_tool);
  1201.         [coefs,longs] = wmemtool('rmb',win_tool,n_membloc3,ind_coefs,ind_longs);
  1202.         switch toolMode
  1203.           case 'denest'
  1204.             varargout{1}  = wthrmngr('dw1ddenoDEN',meth,coefs,longs,alfa); 
  1205.           case {'fixreg','storeg'}
  1206.             varargout{1}  = wthrmngr('dw1ddenoLVL',meth,coefs,longs,alfa);
  1207.         end
  1209.     case 'update_LVL_meth'
  1210.         wdretool('clear_GRAPHICS',win_tool);
  1211.         valTHR = wdretool('compute_LVL_THR',win_tool);
  1212.         utthrw1d('update_LVL_meth',win_tool,valTHR);
  1213.  
  1214.     case 'clear_GRAPHICS'
  1215.         status = wmemtool('rmb',win_tool,n_membloc2,ind_status);
  1216.         if status<1 , return; end
  1218.         wdretool('enable',win_tool,'den','off');
  1219.         wdretool('show_lin_den',win_tool,'off')
  1220.         axe_off = [axe_R_2,axe_app,axe_det];
  1221.         axe_off = axe_off(ishandle(axe_off));
  1222.         lin_Off = findobj(axe_off,'type','line');
  1223.         set(lin_Off,'Visible','off');
  1225.     case 'enable'
  1226.         type = varargin{2};
  1227.         if length(varargin)>2 , enaVal = varargin{3}; else , enaVal = 'on'; end
  1228.         switch type
  1229.           case 'ini'
  1230.             set([men_sav;chk_den],'Enable','off');
  1231.             utthrw1d('status',win_tool,'off');
  1232.             set([pus_dec;sli_bin;edi_bin],'Enable','on');
  1234.           case 'dec'
  1235.             NB_lev = toolATTR.level;
  1236.             set(chk_den,'Value',0);
  1237.             set([men_sav;chk_den],'Enable','off');
  1238.             utthrw1d('status',win_tool,'on');
  1239.             utthrw1d('enable',win_tool,'on');
  1240.             utthrw1d('enable',win_tool,enaVal,[1:NB_lev]);
  1242.           case 'den'
  1243.             set([men_sav;chk_den],'Enable',enaVal);
  1244.             utthrw1d('enable_tog_res',win_tool,enaVal);
  1245.             if strncmpi(enaVal,'on',2) , status = 1; else , status = 0; end
  1246.             wmemtool('wmb',win_tool,n_membloc2,ind_status,status);
  1247.         end
  1249.     case 'clean'
  1250.         calling_opt = varargin{2};        
  1251.         wdretool('show_lin_den',win_tool,'off')
  1252.         lin_den = wmemtool('rmb',win_tool,n_membloc2,ind_lin_den);
  1253.         lin_den = lin_den(ishandle(lin_den));
  1254.         delete(lin_den);
  1255.         wmemtool('wmb',win_tool,n_membloc2,ind_lin_den,[NaN,NaN]); 
  1256.         switch calling_opt
  1257.           case {'load','demo'}
  1258.             obj2del = [findobj(axe_hdl,'type','line'); ...
  1259.                        findobj(axe_hdl,'type','patch')];
  1260.             delete(obj2del)
  1261.             switch toolMode
  1262.               case 'fixreg'
  1263.                 win_title = 'Regression Estimation 1-D - Fixed Design';
  1264.                 set(win_tool,'Name',win_title);
  1265.               case 'storeg'
  1266.                 win_title = 'Regression Estimation 1-D - Stochastic Design';
  1267.                 set(win_tool,'Name',win_title);
  1268.             end
  1269.             sig_name = wmemtool('rmb',win_tool,n_membloc1,ind_sig_name);
  1270.             cbanapar('set',win_tool,'nam',sig_name);
  1271.             utthrw1d('clean_thr',win_tool);
  1273.           case 'bin'
  1274.             switch toolMode
  1275.               case 'denest'  , axetoClean = axe_hdl([2,4:end]);
  1276.               case {'fixreg','storeg'} , axetoClean = axe_hdl(2:end);
  1277.             end
  1278.             obj2del = [findobj(axetoClean,'type','line'); ...
  1279.                        findobj(axetoClean,'type','patch')];
  1280.             delete(obj2del)
  1281.             wdretool('set_Bins',win_tool);
  1282.             wdretool('plot_Processed_Data',win_tool);
  1283.             wdretool('enable',win_tool,'ini');
  1284.         end
  1286.     case 'set_Bins'
  1287.         xdata  = wmemtool('rmb',win_tool,n_membloc1,ind_xdata);
  1288.         nbBINS = get(sli_bin,'Value');
  1289.         toolATTR.NBClasses = nbBINS;
  1290.         wfigmngr('storeValue',win_tool,['WDRE_toolATTR'],toolATTR);
  1291.         switch toolMode
  1292.           case 'denest'
  1293.             color  = colors.denColor;
  1294.             his    = wgethist(xdata,nbBINS);
  1295.             wplothis(axe_R_1,his,color);
  1296.             wtitle('Binned Data','Parent',axe_R_1)
  1298.           case 'fixreg'
  1300.           case 'storeg'
  1301.             color  = colors.denColor;
  1302.             his    = wgethist(xdata,nbBINS);
  1303.             wplothis(axe_L_2,his,color);
  1304.             wtitle('Histogram of X','Parent',axe_L_2)
  1305.         end
  1307.     case 'plot_Processed_Data'
  1308.         [xdata,ydata] = wmemtool('rmb',win_tool,n_membloc1, ...
  1309.                                        ind_xdata,ind_ydata);
  1310.         [sig_pro,xval] = wedenreg(toolATTR,xdata,ydata);
  1311.         wmemtool('wmb',win_tool,n_membloc3,ind_sig,sig_pro);
  1312.    
  1313.         % Plotting Processed data & initial Estimation(s).
  1314.         %-------------------------------------------------        
  1315.         switch toolMode
  1316.           case {'fixreg','storeg'}
  1317.             color = colors.sigColor;
  1318.             lin_processed = plotline(axe_R_1,xval,sig_pro,color,tag_dat);            
  1319.             color = colors.estColor;
  1320.             lin_den = wmemtool('rmb',win_tool,n_membloc2,ind_lin_den);
  1321.             lin_den(1) = plotline(axe_L_1,xval,sig_pro,color,'');
  1322.             lin_den(2) = plotline(axe_R_1,xval,sig_pro,color,'');
  1323.             set(lin_den(1:2),'Visible','Off');
  1324.             ylim = getylim(ydata);
  1325.             set(axe_L_1,'Ylim',ylim);
  1326.             wmemtool('wmb',win_tool,n_membloc2,ind_lin_den,lin_den);
  1327.             hdl_est = plotline(axe_R_2,xval,sig_pro,color,tag_est);
  1328.             set(hdl_est,'Visible','Off');
  1329.             set(axe_hdl,'Xlim',[xval(1) xval(end)]);
  1331.           case {'denest'}
  1332.  
  1333.         end
  1335.     case 'set_axes'
  1336.         %*************************************************************%
  1337.         %** OPTION = 'set_axes' - Set axes positions and visibility **%
  1338.         %*************************************************************%
  1339.         if strcmp(toolMode,'nul') , return; end
  1340.         Pos_Graphic_Area = wmemtool('rmb',win_tool,n_membloc2,ind_gra_area);
  1341.         
  1342.         % Hide axes
  1343.         %-----------
  1344.         if ~isequal(toolMode,'denest')
  1345.             lin_den = wmemtool('rmb',win_tool,n_membloc2,ind_lin_den);
  1346.             if ishandle(lin_den(1)) , vis_den = get(lin_den(1),'Visible'); end
  1347.         end
  1348.         obj_in_axes = findobj(axe_hdl);
  1349.         set(obj_in_axes,'Visible','off');
  1351.         % Parameters initialization.
  1352.         %---------------------------
  1353.         NB_lev = cbanapar('get',win_tool,'lev');
  1354.         indic_vis_lev = getLevels(NB_lev,yLevelDir);
  1355.  
  1356.         % General graphical parameters initialization.
  1357.         %---------------------------------------------
  1358.         pos_win   = get(win_tool,'Position');
  1359.         ecy_up    = 0.06*pos_win(4);
  1360.         ecy_mid_1 = 0.07*pos_win(4);
  1361.         ecy_mid_2 = ecy_mid_1;
  1362.         ecy_mid_3 = ecy_mid_1;
  1363.         ecy_det   = (0.04*pos_win(4))/1.4;
  1364.         ecy_down  = ecy_up;
  1365.         h_gra_rem = Pos_Graphic_Area(4);
  1366.         h_min     = h_gra_rem/12;
  1367.         h_max     = h_gra_rem/5;
  1368.         h_axe_std = (h_min*NB_lev+h_max*(NB_max_lev-NB_lev))/NB_max_lev;
  1369.         h_space   = ecy_up+ecy_mid_1+ecy_mid_2+ecy_mid_3+...
  1370.                     (NB_lev-1)*ecy_det+ecy_down;
  1371.         h_detail  = (h_gra_rem-2*h_axe_std-h_space)/(NB_lev+1);
  1372.         y_low_ini = 1;
  1374.         % Building data axes.
  1375.         %--------------------
  1376.         y_low_ini = y_low_ini-h_axe_std-ecy_up;
  1377.         pos_L_1 = get(axe_L_1,'Position');
  1378.         pos_L_1([2 4]) = [y_low_ini h_axe_std];
  1379.         set(axe_L_1,'Position',pos_L_1);
  1380.         pos_R_1 = get(axe_R_1,'Position');
  1381.         pos_R_1([2 4]) = [y_low_ini h_axe_std];
  1382.         set(axe_R_1,'Position',pos_R_1);
  1383.         axe_vis = [axe_L_1,axe_R_1];
  1384.         
  1385.         y_low_ini = y_low_ini-h_axe_std-ecy_mid_1;
  1386.         pos_L_2  = get(axe_L_2,'Position');
  1387.         pos_L_2([2 4]) = [y_low_ini h_axe_std];
  1388.         set(axe_L_2,'Position',pos_L_2)
  1389.         pos_R_2  = get(axe_R_2,'Position');
  1390.         pos_R_2([2 4]) = [y_low_ini h_axe_std];
  1391.         set(axe_R_2,'Position',pos_R_2)
  1392.         switch toolMode
  1393.           case 'denest' , axe_vis = [axe_vis,axe_L_2,axe_R_2];
  1394.           case 'fixreg' , axe_vis = [axe_vis,axe_R_2];
  1395.           case 'storeg' , axe_vis = [axe_vis,axe_L_2,axe_R_2];
  1396.         end
  1398.         % Position for approximation axes on the right part.
  1399.         %---------------------------------------------------
  1400.         y_low_ini = y_low_ini-h_detail-ecy_mid_2;
  1401.         pos_axes = pos_R_2;
  1402.         pos_y   = [y_low_ini , h_detail];
  1403.         pos_axes([2 4]) = pos_y;
  1404.         set(axe_app,'Position',pos_axes);
  1405.         axe_vis = [axe_vis,axe_app];
  1406.         y_low_ini = y_low_ini-ecy_mid_3+ecy_det;
  1408.         % Position for details axes on the left part.
  1409.         %--------------------------------------------          
  1410.         pos_y  = [y_low_ini , h_detail];
  1411.         for j = 1:NB_lev
  1412.             i_axe    = indic_vis_lev(j);
  1413.             axe_act  = axe_cfs(i_axe);
  1414.             pos_axes = get(axe_act,'Position');
  1415.             pos_y(1) = pos_y(1)-h_detail-ecy_det;
  1416.             pos_axes([2 4]) = pos_y;
  1417.             set(axe_act,'Position',pos_axes);
  1418.             axe_vis = [axe_vis axe_act];
  1419.         end
  1420.         i_axe_cfs_min = indic_vis_lev(1);
  1421.         i_axe_cfs_max = indic_vis_lev(NB_lev);
  1423.         % Position for details axes on the right part.
  1424.         %---------------------------------------------
  1425.         pos_y   = [y_low_ini , h_detail];
  1426.         for j = 1:NB_lev
  1427.             i_axe    = indic_vis_lev(j);
  1428.             axe_act  = axe_det(i_axe);
  1429.             pos_axes = get(axe_act,'Position');
  1430.             pos_y(1) = pos_y(1)-h_detail-ecy_det;
  1431.             pos_axes([2 4]) = pos_y;
  1432.             set(axe_act,'Position',pos_axes);
  1433.             axe_vis  = [axe_vis axe_act];
  1434.         end
  1435.         i_axe_det_min = indic_vis_lev(1);
  1437.         % Modification of app_text.
  1438.         %--------------------------
  1439.         status = wmemtool('rmb',win_tool,n_membloc2,ind_status);
  1440.         if status==0
  1441.             txt_app = txt_hdl(end);
  1442.             num_app = NB_lev;
  1443.             set(txt_app,'string',['a' wnsubstr(abs(num_app))]);
  1444.         end
  1446.         % Set axes.
  1447.         %-----------
  1448.         xt = get(axe_L_1,{'Xtick','XtickLabel'});
  1449.         ind_axe_cfs = [i_axe_cfs_min:i_axe_cfs_max-1];
  1450.         set(axe_cfs(ind_axe_cfs),'Xtick',[],'XtickLabel',[]);
  1451.         set([axe_cfs(i_axe_cfs_max),axe_det(i_axe_cfs_max)], ...
  1452.             'Xtick',xt{1},'XtickLabel',xt{2} , ...
  1453.             'XtickMode','auto','XTickLabelMode','auto' ...
  1454.             )
  1455.         titles = get([axe_cfs;axe_det],'title');
  1456.         titles = cat(1,titles{:});
  1457.         set(titles,'String','');
  1458.         obj_in_axes_vis = findobj(axe_vis);
  1459.         set(obj_in_axes_vis,'Visible','on');
  1460.         if ~isequal(toolMode,'denest')
  1461.             if ishandle(lin_den(1)) , set(lin_den,'Visible',vis_den); end
  1462.         end
  1463.        %  hdl_den = utthrw1d('get',win_tool,'handleTHR')
  1464.        %  set(hdl_den,'color','g')
  1465.         
  1466.         % Setting axes title
  1467.         %--------------------
  1468.         switch toolMode
  1469.           case 'fixreg'
  1470.             wtitle('Data (X,Y)','Parent',axe_L_1);
  1471.             wtitle('Regression Estimate Y = f(X)','Parent',axe_R_2);
  1473.           case 'storeg'
  1474.             wtitle('Data (X,Y)','Parent',axe_L_1)
  1475.             wtitle('Histogram of X','Parent',axe_L_2);
  1476.             wtitle('Regression Estimate Y = f(X)','Parent',axe_R_2);
  1477.             wtitle('Processed Data (X,Y)','Parent',axe_R_1);
  1479.           case 'denest'
  1480.             wtitle('Data','Parent',axe_L_1);
  1481.             wtitle(sprintf('%d bins - Histogram',def_MinBIN),'Parent',axe_L_2);
  1482.             wtitle('Binned Data','Parent',axe_R_1);
  1483.             wtitle('Density Estimate','Parent',axe_R_2);
  1484.         end
  1485.         wtitle('Details coefficients','Parent',axe_cfs(i_axe_cfs_min));
  1486.         wtitle('Details','Parent',axe_det(i_axe_det_min));
  1487.         wtitle('Approximation','Parent',axe_app);
  1489.     case 'close'
  1491.     otherwise
  1492.         errargt(mfilename,'Unknown Option','msg');
  1493.         error('*');
  1494. end
  1496. %=============================================================================%
  1497. % INTERNAL FUNCTIONS
  1498. %=============================================================================%
  1499. function setDEMOS(m_demo,funcName,str_numwin,demoSET,sepFlag)
  1501. beg_call_str = [funcName '(''demo'',' str_numwin ','''];
  1502. [nbDEM,nbVAL] = size(demoSET);
  1503. for k=1:nbDEM
  1504.     typ = demoSET{k,1};
  1505.     nam = demoSET{k,2};
  1506.     fil = demoSET{k,3};
  1507.     wav = demoSET{k,4};
  1508.     len = length(wav);
  1509.     wavDum = [wav , blanks(4-len)];
  1510.     lev = int2str(demoSET{k,5});
  1511.     if nbVAL>5 , par = demoSET{k,6}; else , par = '[]'; end
  1512.     libel = ['with  ' wavDum '   at level ' lev '  --->  ' nam];
  1513.     action = [beg_call_str ...
  1514.               typ ''',''' fil ''',''' wav  ''',' lev ',' par ');'];
  1515.     if sepFlag & (k==1) , sep = 'on'; else , sep = 'off'; end
  1516.     uimenu(m_demo,'Label',libel,'Separator',sep,'Callback',action);
  1517. end
  1518. %-----------------------------------------------------------------------------%
  1519. function indic_vis_lev = getLevels(level,yDir)
  1521. indic_vis_lev = [1:level]';
  1522. if yDir==-1 , indic_vis_lev = flipud(indic_vis_lev); end
  1523. %-----------------------------------------------------------------------------%
  1524. function [det,app,sig] = getFullDecTec(coefs,longs,wname,level)
  1526. det  = wrmcoef('d',coefs,longs,wname,[1:level]);
  1527. app  = wrcoef('a',coefs,longs,wname,level);
  1528. if nargout<3 , return; end
  1529. sig  = waverec(coefs,longs,wname);
  1530. %-----------------------------------------------------------------------------%
  1531. function ylim = getylim(sig)
  1533. mini = min(sig);
  1534. maxi = max(sig);
  1535. dy   = maxi-mini;
  1536. tol  = sqrt(eps);
  1537. if abs(dy)<tol
  1538.     maxi = maxi+tol;
  1539.     mini = mini-tol;
  1540. end;
  1541. ylim = [mini maxi]+0.05*dy*[-1 1];
  1542. %-----------------------------------------------------------------------------%
  1543. function ll = plotline(axe,x,y,color,tag,ylimplus)
  1545. vis = get(axe,'visible');
  1546. ll = findobj(axe,'type','line','tag',tag);
  1547. if isempty(ll)
  1548.      ll = line(...
  1549.             'Parent',axe,'Xdata',x,'Ydata',y, ...
  1550.             'Color',color,'Visible',vis,'tag',tag);
  1551. else
  1552.      set(ll,'Xdata',x,'Ydata',y,'Color',color,'Visible',vis,'tag',tag);
  1553. end
  1554. if nargin<6
  1555.     ylim = getylim(y);
  1556. else
  1557.     ylim = getylim([y(:) ; ylimplus(:)]);
  1558. end
  1559. set(axe,'Ylim',ylim);
  1560. %-----------------------------------------------------------------------------%
  1561. function loc = coefsLOC(idx,lev,lf,lx)
  1562. %COEFSLOC coefficient location
  1564. up  = idx;
  1565. low = idx;
  1566. for j=1:lev
  1567.     low = 2*low+1-lf;
  1568.     up  = 2*up;
  1569. end
  1570. loc = max(1,min(lx,round((low+up)/2)));
  1571. %-----------------------------------------------------------------------------%
  1572. function h = plotstem(axe,x,y,color,flgzero,tag,ylimplus)
  1573. %PLOTSTEM Plot discrete sequence data.
  1575. vis = get(axe,'visible');
  1576. xAxeColor = get(axe,'xcolor');
  1577. q = [min(x) max(x)];
  1578. h = NaN*ones(1,4);
  1579. h(1) = line(...
  1580.               'Xdata',[q(1) q(2)],'Ydata',[0 0],...
  1581.               'Parent',axe,'color',xAxeColor ...
  1582.               );
  1584. indZ = find(abs(y)<eps);
  1585. xZ   = x(indZ);
  1586. yZ   = y(indZ);
  1587. x(indZ) = [];
  1588. y(indZ) = [];
  1590. n = length(x);
  1591. if n>0
  1592.     MSize = 3; Mtype = 'o';
  1593.     MarkerEdgeColor = color;
  1594.     MarkerFaceColor = color;
  1595.     xx = [x;x;NaN*ones(size(x))];
  1596.     yy = [zeros(1,n);y;NaN*ones(size(y))];
  1597.     h(2) = line(...
  1598.                   'Xdata',xx(:),'Ydata',yy(:),...
  1599.                   'Parent',axe,'LineStyle',...
  1600.                   '-','Color',color...
  1601.                   );
  1602.     h(3) = line(...
  1603.                   'Xdata',x,'Ydata',y,...
  1604.                   'Parent',axe,...
  1605.                   'Marker',Mtype, ...
  1606.                   'MarkerEdgeColor',MarkerEdgeColor, ...
  1607.                   'MarkerFaceColor',MarkerFaceColor, ...
  1608.                   'MarkerSize',MSize, ...
  1609.                   'LineStyle','none',...
  1610.                   'Color',color ...
  1611.                   );
  1612. end
  1614. nZ = length(xZ);
  1615. if flgzero & (nZ>0)
  1616.     MSize = 3; Mtype = 'o';
  1617.     h(4) = line(...
  1618.                   'Xdata',xZ,'Ydata',yZ,...
  1619.                   'Parent',axe,...
  1620.                   'Marker',Mtype, ...
  1621.                   'MarkerEdgeColor',xAxeColor, ...
  1622.                   'MarkerFaceColor',xAxeColor, ...
  1623.                   'MarkerSize',MSize, ...
  1624.                   'LineStyle','none',...
  1625.                   'Color',xAxeColor...
  1626.                   );
  1627. end
  1629. set(h(ishandle(h)),'Visible',vis,'tag',tag);
  1630. if nargin<7
  1631.     ylim = getylim(y);
  1632. else
  1633.     ylim = getylim([y(:) ; ylimplus(:)]);
  1634. end
  1635. set(axe,'Ylim',ylim);
  1636. %-----------------------------------------------------------------------------%
  1637. %=============================================================================%
  1640. %=============================================================================%
  1641. function [ysig,xsig,coefs,longs] = wedenreg(toolATTR,x,y)
  1642. %WEDENREG Density and Regression Estimation.
  1644. %== Initialisation ==========================================================%
  1645. toolMode  = toolATTR.toolMode;
  1646. level     = toolATTR.level;
  1647. wname     = toolATTR.wname;
  1648. NBClasses = toolATTR.NBClasses;
  1649. if nargin<3 , y = ones(size(x)); end
  1650. %============================================================================%
  1652. %== Traitement sur les couples (Xm,Ym) ======================================%
  1653. interpol = 0; % pas d'interpolation.
  1654. [sx,sy] = prepxy(toolMode,x,y,NBClasses,interpol);
  1655. ind_sx  = find(sx>0);
  1656. %============================================================================%
  1658. %============================================================================%
  1659. % Calcul de la grille en x.
  1660. %-------------------------
  1661. xmin = min(x);
  1662. xmax = max(x);
  1663. xsig = [0:NBClasses-1]/(NBClasses-1);
  1664. xsig = (xmax-xmin)*xsig+xmin;
  1665. %============================================================================%
  1667. %============================================================================%
  1668. % Normalization & Processed Data.
  1669. %-------------------------------
  1670. switch toolMode
  1671.     case {'denest'}
  1672.         ysig = sy;
  1673.         delta = xsig(2)-xsig(1);
  1674.         integ = sum(ysig)*delta;
  1675.         ysig = ysig/integ;
  1676.         
  1677.     case {'fixreg','storeg'}
  1678.         ysig = zeros(size(sx));
  1679.         ysig(ind_sx) = sy(ind_sx)./sx(ind_sx);
  1680. end
  1682. % Decomposition of processed signal.
  1683. %-----------------------------------
  1684. if nargout>2
  1685.     [coefs,longs] = wavedec(ysig,level,wname);
  1686.     varargout = {ysig,xsig,coefs,longs};
  1687. else
  1688.     varargout = {ysig,xsig};
  1689. end
  1690. %============================================================================%
  1691. % INTERNAL FUNCTIONS for WEDENREG
  1692. %============================================================================%
  1693. %----------------------------------------------------------------------------%
  1694. function [sx,sy] = prepxy(option,x,y,NBClasses,interpol)
  1695. %PREPXY Traitement des couples (X,Y)
  1696. %   
  1697. %   [sx,sy] = prepxy(option,x,y,NBClasses,interpol)
  1698. %
  1699. %   Entrees:
  1700. %       option    = 'denest','fixreg','storeg'
  1701. %       (x,y)     = couples des donnees (y==1 pour 'estidens') 
  1702. %       NBClasses = Nombre de classes de la grille sur [0,1].
  1703. %                   La grille sur l'intervalle [0,1] est:
  1704. %                   xgrid = [0.5:NBClasses-0.5]/NBClasses
  1705. %
  1706. %   Sorties:
  1707. %       sx = nombre d'element dans chaque classe.
  1708. %       sy = somme des elements de chaque classe.
  1709. %               (sy == sx pour 'denest' )
  1711. %----------------------------------------%
  1712. % On change l'intervalle.                %
  1713. % On travaille sur l'intervalle [0,1]    %
  1714. %----------------------------------------%
  1715. xmin  = min(x);
  1716. xmax  = max(x);
  1717. x     = (x-xmin)/(xmax-xmin);
  1718. I1    = find(x==1);
  1719. x(I1) = x(I1)-eps;
  1721. %----------------------------------------%
  1722. % La grille sur l'intervalle [0,1] est : %
  1723. % xf = [0.5:NBClasses-0.5]/NBClasses     %
  1724. %----------------------------------------%
  1725. ex = round(NBClasses*x+0.5); 
  1727. %----------------------------------------%
  1728. % Recherche des valeurs repetees et      %
  1729. % manquantes de ex.                      %
  1730. % Calcul des sommes des valeurs de y     %
  1731. % pour chaque ex distinct.               %
  1732. %----------------------------------------%
  1733. lx = length(ex);
  1734. sx = full(sum(sparse(1:lx,ex,1,lx,NBClasses)));
  1736. switch option
  1737.    case 'denest' , sy = sx;
  1738.    case 'fixreg' , sy = full(sum(sparse(1:lx,ex,y,lx,NBClasses)));
  1739.    case 'storeg'
  1740.      sy = full(sum(sparse(1:lx,ex,y,lx,NBClasses)));
  1741.      if nargin<5 , return; end
  1743.      % PROBLEME DES VALEURS MANQUANTES: interpol = 1;
  1744.      %----------------------------------------------
  1745.      if interpol
  1746.          % Interpolation lineaire pour les trous.
  1747.          Ind_0  = find(sx==0);
  1748.          Ind_sx = (find(sx>0))';
  1749.          for k=1:length(Ind_0)
  1750.              Ik = Ind_0(k);
  1751.              av = find(Ind_sx<Ik);
  1752.              ap = find(Ind_sx>Ik);
  1753.              Iav = Ind_sx(av(length(av)));
  1754.              Iap = Ind_sx(ap(1));
  1755.              sy(Ik) = ( (sy(Iap)*(Ik-Iav))/sx(Iap)+ ...
  1756.                         (sy(Iav)*(Iap-Ik))/sx(Iav) )/(Iap-Iav);
  1757.              sx(Ik) = 1;
  1758.          end
  1759.      end
  1760. end
  1761. %----------------------------------------------------------------------------%
  1762. %============================================================================%