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

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

波变换

开发平台：

Matlab

sw2dtool.m：源码内容
  1. function varargout = sw2dtool(option,varargin)
  2. %SW2DTOOL Stationary Wavelet Transform 2-D tool.
  3. %   VARARGOUT = SW2DTOOL(OPTION,VARARGIN)
  5. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 02-Mar-1998.
  6. %   Last Revision: 05-Feb-2004.
  7. %   Copyright 1995-2004 The MathWorks, Inc.
  8. %   $Revision: 1.9.4.2 $  $Date: 2004/03/15 22:41:46 $
  10. % Test inputs.
  11. %-------------
  12. if nargin==0 , option = 'create'; end
  13. [option,winAttrb] = utguidiv('ini',option,varargin{:});
  15. % Default values.
  16. %----------------
  17. max_lev_anal = 5;
  18. def_lev_anal = 2;
  19. def_nbcolors = 255;
  20. str_dir_det  = strvcat('Horizontal','Diagonal','Vertical');
  22. % Memory Blocks of stored values.
  23. %================================
  24. % MB1.
  25. %-----
  26. n_membloc1    = 'MB_1';
  27. ind_status    = 1;
  28. ind_sav_menu  = 2;
  29. ind_filename  = 3;
  30. ind_pathname  = 4;
  31. ind_img_name  = 5;
  32. ind_img_t_nam = 6;
  33. ind_NB_lev    = 7;
  34. ind_wave      = 8;
  35. nb1_stored    = 8;
  37. % MB2.
  38. %-----
  39. n_membloc2   = 'MB_2';
  40. ind_pus_dec  = 1;
  41. ind_axe_hdl  = 2;
  42. ind_txt_hdl  = 3;
  43. ind_gra_area = 4;
  44. nb2_stored   = 4;
  46. % MB3.
  47. %-----
  48. n_membloc3 = 'MB_3';
  49. ind_coefs  = 1;
  50. nb3_stored = 1;
  52. if ~isequal(option,'create') , win_tool = varargin{1}; end
  53. switch option
  54.     case 'create'
  55.         % Get Globals.
  56.         %-------------
  57.         [Def_Btn_Height,Y_Spacing] = ...
  58.             mextglob('get','Def_Btn_Height','Y_Spacing');
  60.         % Window initialization.
  61.         %-----------------------
  62.         win_title = 'Stationary Wavelet Transform De-noising 2-D';
  63.         [win_tool,pos_win,win_units,str_numwin,...
  64.            frame0,pos_frame0,Pos_Graphic_Area,pus_close] = ...
  65.              wfigmngr('create',win_title,winAttrb,'ExtFig_Tool_3',mfilename,1,1,0);
  66.         if nargout>0 , varargout{1} = win_tool; end
  68. % Add Help for Tool.
  69. %------------------
  70. wfighelp('addHelpTool',win_tool, ...
  71. 'Two-Dimensional Analysis for De-&noising','SW2D_GUI');
  73. % Add Help Item.
  74. %----------------
  75. wfighelp('addHelpItem',win_tool,'Stationary Wavelet Transform','SWT');
  76. wfighelp('addHelpItem',win_tool,'Available Methods','COMP_DENO_METHODS');
  77. wfighelp('addHelpItem',win_tool,'Loading and Saving','SW2D_LOADSAVE');
  79.         % Menu construction for current figure.
  80.         %--------------------------------------
  81.         m_files = wfigmngr('getmenus',win_tool,'file');
  82.         cba_menu = [mfilename '(''load'','  str_numwin ');'];
  83.         uimenu(m_files,...
  84.                       'Label','&Load Image ', ...
  85.                       'Position',1,           ...
  86.                       'Callback', cba_menu    ...
  87.                        );
  88.         cba_menu = [mfilename '(''save'','  str_numwin ');'];
  89.         m_save = uimenu(m_files,...
  90.                         'Label','&Save De-Noised Image ',...
  91.                         'Position',2,        ...
  92.                         'Enable','Off',      ...
  93.                         'Callback', cba_menu ...
  94.                         );
  95.         demoSET = {...
  96.           'Noisy Woman'   , 'noiswom'  , 'haar', 3 , '{''penallo'',46.12}'  ; ...
  97.           'Noisy Woman'   , 'noiswom'  , 'haar', 5 , '{''penallo'',48.62}'  ; ...
  98.           'Noisy Woman'   , 'noiswom'  , 'db3' , 4 , '{''penallo'',NaN}'  ; ...
  99.           'Noisy Barbara' , 'nbarb1'   , 'db1' , 4 , '{}'  ; ...
  100.           'Noisy Sinsin'  , 'noissi2d' , 'db1' , 2 , '{}'   ...
  101.           };
  103.         m_demo = uimenu(m_files,'Label','&Example Analysis','Position',3);
  104.         beg_call_str = [mfilename '(''demo'',' str_numwin ','''];
  105.         nbDEM = size(demoSET,1);
  106.         for k=1:nbDEM
  107.             nam = demoSET{k,1};
  108.             fil = demoSET{k,2};
  109.             wav = demoSET{k,3};
  110.             lev = int2str(demoSET{k,4});
  111.             par = demoSET{k,5};       
  112.             libel = ['with ' wav ' at level ' lev  '  --->  ' nam];
  113.             action = [beg_call_str fil ''',''' wav ''',' lev ',' par ');'];
  114.             uimenu(m_demo,'Label',libel,'Callback',action);
  115.         end
  117.         % Begin waiting.
  118.         %---------------
  119.         wwaiting('msg',win_tool,'Wait ... initialization');
  121.         % General parameters initialization.
  122.         %-----------------------------------
  123.         dy = Y_Spacing;
  124.         str_pus_dec = 'Decompose Image';
  126.         % Command part of the window.
  127.         %============================
  128.         comFigProp = {'Parent',win_tool,'Unit',win_units};
  130.         % Data, Wavelet and Level parameters.
  131.         %------------------------------------
  132.         xlocINI = pos_frame0([1 3]);
  133.         ytopINI = pos_win(4)-dy;
  134.         toolPos = utanapar('create',win_tool, ...
  135.                     'xloc',xlocINI,'top',ytopINI,...
  136.                     'enable','off',        ...
  137.                     'wtype','dwt',         ...
  138.                     'deflev',def_lev_anal, ...
  139.                     'maxlev',max_lev_anal  ...
  140.                     );
  142.         % Decompose pushbutton.
  143.         %----------------------
  144.         h_uic = 3*Def_Btn_Height/2;
  145.         y_uic = toolPos(2)-h_uic-2*dy;
  146.         w_uic = (3*pos_frame0(3))/4;
  147.         x_uic = pos_frame0(1)+(pos_frame0(3)-w_uic)/2;
  148.         pos_pus_dec = [x_uic, y_uic, w_uic, h_uic];
  149.         pus_dec = uicontrol(...
  150.                             'Parent',win_tool,      ...
  151.                             'Style','Pushbutton',   ...
  152.                             'Unit',win_units,       ...
  153.                             'Position',pos_pus_dec, ...
  154.                             'String',xlate(str_pus_dec),   ...
  155.                             'Enable','off',         ...
  156.                             'Interruptible','On'    ...
  157.                             );
  159.         % De-noising tool.
  160.         %-----------------
  161.         ytopTHR = pos_pus_dec(2)-4*dy;
  162.         toolPos = utthrw2d('create',win_tool, ...
  163.                     'xloc',xlocINI,'top',ytopTHR,...
  164.                     'ydir',-1, ...
  165.                     'levmax',def_lev_anal,    ...
  166.                     'levmaxMAX',max_lev_anal, ...
  167.                     'status','Off',  ...
  168.                     'toolOPT','deno' ...
  169.                     );
  171.         % Adding colormap GUI.
  172.         %---------------------
  173.         briflag = (max_lev_anal<6); 
  174.         utcolmap('create',win_tool, ...
  175.                  'xloc',xlocINI, ...
  176.                  'briflag',briflag, ...
  177.                  'enable','off');
  179.         % Callbacks update.
  180.         %------------------
  181.         hdl_den = utthrw2d('handles',win_tool);
  182.         utanapar('set_cba_num',win_tool,[m_files;hdl_den(:)]);
  183.         pop_lev = utanapar('handles',win_tool,'lev');
  184.         end_cba = [str_numwin ',' num2mstr([pop_lev]) ');'];
  185.         cba_pop_lev = [mfilename '(''update_level'',' end_cba];
  186.         cba_pus_dec = [mfilename '(''decompose'','  end_cba];
  187.         set(pop_lev,'Callback',cba_pop_lev);
  188.         set(pus_dec,'Callback',cba_pus_dec);
  190.         % General graphical parameters initialization.
  191.         %--------------------------------------------
  192.         fontsize    = wmachdep('fontsize','normal',9,max_lev_anal);
  193.         txtfontsize = 14;
  195.         % Axes construction parameters.
  196.         %------------------------------
  197.         NB_lev    = max_lev_anal;
  198.         bdx       = 0.08*pos_win(3);
  199.         ecy_up    = 0.06*pos_win(4);
  200.         ecy_mid_1 = 0.08*pos_win(4);
  201.         ecy_mid_2 = ecy_up;
  202.         ecy_down  = ecy_up;
  203.         ecy_det   = (0.04*pos_win(4))/1.4;
  204.         w_gra_rem = Pos_Graphic_Area(3);
  205.         h_gra_rem = Pos_Graphic_Area(4);
  206.         w_left    = (w_gra_rem-3*bdx)/2;
  207.         w_right   = w_left;
  208.         w_medium  = w_left;
  209.         w_detail  = (w_gra_rem-4*bdx)/3+bdx/2;
  210.         x_left    = bdx;
  211.         x_right   = x_left+w_left+5*bdx/4;
  212.         x_medium  = (w_gra_rem-w_medium)/2;
  213.         h_min     = h_gra_rem/4;
  214.         h_max     = h_gra_rem/3;
  215.         h_image   = (h_min*(NB_lev-1)+h_max*(max_lev_anal-NB_lev))/(max_lev_anal-1);
  216.         h_histo   = h_gra_rem/8;
  217.         h_gra_rem = h_gra_rem-h_histo-h_image-ecy_up-ecy_mid_1-ecy_mid_2;
  218.         h_detail  = (h_gra_rem-ecy_down-(NB_lev-1)*ecy_det)/NB_lev;
  219.         y_low_ini = pos_win(4);
  221.         % Building data axes.
  222.         %--------------------
  223.         commonProp   = {...
  224.            'Parent',win_tool,...
  225.            'Visible','off',  ...
  226.            'Units',win_units,...
  227.            'Drawmode','fast',...
  228.            'Box','On'...
  229.            };
  230.         y_low_ini   = y_low_ini-h_image-ecy_up;
  231.         pos_left    = [x_left y_low_ini w_left h_image];
  232.         axe_left_1  = axes(commonProp{:},'Position',pos_left);
  233.         pos_right   = [x_right y_low_ini w_right h_image];
  234.         axe_right_1 = axes(commonProp{:},'Position',pos_right);
  236.         y_low_ini   = y_low_ini-h_histo-ecy_mid_1;
  237.         pos_medium  = [x_medium y_low_ini w_medium h_histo];
  238.         axe_medium  = axes(commonProp{:},'Position',pos_medium);
  239.         y_low_ini   = y_low_ini-ecy_mid_2+ecy_det;
  241.         % Building histograms axes.
  242.         %==========================
  243.         axe_hist = zeros(3,NB_lev);
  244.         commonProp = {commonProp{:},'XTicklabelMode','manual','XTickLabel',[]};
  246.         % Building axes on the left part.
  247.         %-------------------------------    
  248.         txt_left = zeros(NB_lev,1);
  249.         x_left   = bdx;
  250.         y_left   = y_low_ini;
  251.         pos_left = [x_left y_left w_detail h_detail];
  252.         for j = 1:NB_lev
  253.             k = NB_lev-j+1;
  254.             pos_left(2) = pos_left(2)-pos_left(4)-ecy_det;
  255.             axe_hist(1,k) = axes(commonProp{:}, ...
  256.                                 'Position',pos_left,'NextPlot','add');
  257.             str_txt     = ['L' wnsubstr(k)];
  258.             txt_left(k) = txtinaxe('create',str_txt, ...
  259.                              axe_hist(1,k),'l','off','bold',txtfontsize);
  260.         end
  262.         % Building details axes on the middle part.
  263.         %----------------------------------------
  264.         x_mid   = x_left+w_detail+bdx/2;
  265.         y_mid   = y_low_ini;
  266.         pos_mid = [x_mid y_mid w_detail h_detail];
  267.         for k = 1:NB_lev
  268.             j = NB_lev-j+1;
  269.             pos_mid(2)    = pos_mid(2)-pos_mid(4)-ecy_det;
  270.             axe_hist(2,k) = axes(commonProp{:},'Position',pos_mid);
  271.         end
  272.         
  273.         % Building details axes on the right part.
  274.         %-----------------------------------------
  275.         x_right   = x_mid+w_detail+bdx/2;
  276.         y_right   = y_low_ini;
  277.         pos_right = [x_right y_right w_detail h_detail];
  278.         for j = 1:NB_lev
  279.             k = NB_lev-j+1;
  280.             pos_right(2)  = pos_right(2)-pos_right(4)-ecy_det;
  281.             axe_hist(3,k) = axes(commonProp{:},'Position',pos_right);
  282.         end
  284.         %  Normalization.
  285.         %----------------
  286.         Pos_Graphic_Area = wfigmngr('normalize',win_tool,Pos_Graphic_Area);
  287.         drawnow
  289.         % Memory blocks update.
  290.         %----------------------
  291.         axes_hdl = {axe_left_1,axe_right_1,axe_medium,axe_hist};
  292.         utthrw2d('set',win_tool,'axes',axe_hist);
  293.         wmemtool('ini',win_tool,n_membloc1,nb1_stored);
  294.         wmemtool('ini',win_tool,n_membloc2,nb2_stored);
  295.         wmemtool('wmb',win_tool,n_membloc1, ...
  296.                        ind_status,0,        ...
  297.                        ind_sav_menu,m_save  ...
  298.                        );
  299.         wmemtool('wmb',win_tool,n_membloc2, ...
  300.                        ind_pus_dec,pus_dec, ...
  301.                        ind_axe_hdl,axes_hdl,...
  302.                        ind_gra_area,Pos_Graphic_Area, ...
  303.                        ind_txt_hdl,txt_left ...
  304.                        );
  306.         % End waiting.
  307.         %---------------
  308.         wwaiting('off',win_tool);
  310.     case {'load','demo'}
  311.         % Loading file.
  312.         %--------------
  313.         switch option
  314.             case 'load'
  315.                 imgFileType = ['*.mat;*.bmp;*.hdf;*.jpg;' ...
  316.                         '*.jpeg;*.pcx;*.tif;*.tiff;*.gif'];
  317.                 [imgInfos,img_ori,map,ok] = utguidiv('load_img',win_tool, ...
  318.                     imgFileType,'Load Image',def_nbcolors);
  319.                 
  320.             case 'demo'
  321.                 img_Name = deblank(varargin{2});
  322.                 wav_Name = deblank(varargin{3});
  323.                 lev_Anal = varargin{4};
  324.                 if length(varargin)>4 & ~isempty(varargin{5})
  325.                     par_Demo = varargin{5};
  326.                 else
  327.                     par_Demo = '';
  328.                 end
  329.                 filename = [img_Name '.mat'];
  330.                 pathname = utguidiv('WTB_DemoPath',filename);
  331.                 [imgInfos,img_ori,map,ok] = utguidiv('load_dem2D',win_tool, ...
  332.                     pathname,filename,def_nbcolors);
  333.         end
  334.         if ~ok, return; end
  336.         % Begin waiting.
  337.         %---------------
  338.         wwaiting('msg',win_tool,'Wait ... cleaning');
  340.         % Storing values. 
  341.         %-----------------
  342.         wmemtool('wmb',win_tool,n_membloc1, ...
  343.                        ind_status,0,        ...
  344.                        ind_filename,imgInfos.filename,  ...
  345.                        ind_pathname,imgInfos.pathname,  ...
  346.                        ind_img_name,imgInfos.name,      ...
  347.                        ind_img_t_nam,imgInfos.true_name ...
  348.                        );
  350.         % Cleaning and setting GUI. 
  351.         %--------------------------
  352.         cbanapar('enable',win_tool,'Off');
  353.         dynvtool('stop',win_tool);
  354.         ax_hdl  = wmemtool('rmb',win_tool,n_membloc2,ind_axe_hdl);
  355.         ax_hdl  = [ax_hdl{1:3},ax_hdl{4}(:)'];
  356.         obj2del = [findobj(ax_hdl,'type','line');  ...
  357.                    findobj(ax_hdl,'type','patch'); ...
  358.                    findobj(ax_hdl,'type','image')];
  359.         delete(obj2del)
  360.         utthrw2d('clean_thr',win_tool);
  362.         % Setting analysis  & GUI values.
  363.         %--------------------------------
  364.         levm   = wmaxlev(imgInfos.size,'haar');
  365.         levmax = min(levm,max_lev_anal);
  366.         if isequal(imgInfos.true_name,'X')
  367.             img_Name = imgInfos.name;
  368.         else
  369.             img_Name = imgInfos.true_name;
  370.         end
  371.         if isequal(option,'demo')
  372.             anaPar = {'wav',wav_Name};
  373.         else
  374.             lev_Anal = def_lev_anal;
  375.             anaPar = {};
  376.         end
  377.         strlev = int2str([1:levmax]');
  378.         anaPar = {anaPar{:},'n_s',{img_Name,imgInfos.size}, ...
  379.                   'lev',{'String',strlev,'Value',lev_Anal}};
  380.         cbanapar('set',win_tool,anaPar{:});
  381.         NB_ColorsInPal = size(map,1);
  382.         if imgInfos.self_map , arg = map; else , arg = []; end
  383.         cbcolmap('set',win_tool,'pal',{'pink',NB_ColorsInPal,'self',arg});
  384.         utthrw2d('set',win_tool,'position',{1,lev_Anal})
  386.         % Enabling HDLG.
  387.         %---------------
  388.         sw2dtool('enable',win_tool,'ini','on');
  390.         % Setting axes. 
  391.         %--------------
  392.         sw2dtool('set_axes',win_tool);
  394.         % Initial drawing
  395.         %----------------
  396.         axe_hdl = wmemtool('rmb',win_tool,n_membloc2,ind_axe_hdl);
  397.         axe_ori = axe_hdl{1};
  398.         hdl_ori = image([1 imgInfos.size(1)],[1 imgInfos.size(2)],img_ori, ...
  399.                         'Parent',axe_ori);
  400.         wtitle('Original Image (I)','Parent',axe_ori)
  401.         utthrw2d('set',win_tool,'handleORI',hdl_ori);
  403.         % Analyze and De-noise.
  404.         %----------------------
  405.         if isequal(option,'demo')
  406.             sw2dtool('decompose',win_tool);
  407.             if ~isempty(par_Demo)
  408.                 den_Meth = par_Demo{1};
  409.                 if length(par_Demo)>1
  410.                     thr_Val  = par_Demo{2};
  411.                 else
  412.                     thr_Val  = NaN;
  413.                 end
  414.                 if isequal(den_Meth,'penallo')
  415.                     utthrw2d('demo',win_tool,'swt2',den_Meth,thr_Val);
  416.                 end
  417.             end
  418.             sw2dtool('denoise',win_tool);
  419.         end
  420.         cbanapar('enable',win_tool,'On');
  422.         % End waiting.
  423.         %-------------
  424.         wwaiting('off',win_tool);
  426.     case 'save'
  427.         % Testing file.
  428.         %--------------
  429.         [filename,pathname,ok] = utguidiv('test_save',win_tool, ...
  430.                                      '*.mat','Save De-Noised Image');
  431.         if ~ok, return; end
  433.         % Begin waiting.
  434.         %--------------
  435.         wwaiting('msg',win_tool,'Wait ... saving');
  437.         % Getting colormap.
  438.         %------------------
  439.         map = cbcolmap('get',win_tool,'self_pal');
  440.         if isempty(map) , map = pink(def_nbcolors); end
  442.         % Get de-noised signal.
  443.         %----------------------
  444.         hdl_den = utthrw2d('get',win_tool,'handleTHR');
  445.         X = get(hdl_den,'Cdata');
  447.         % Saving file.
  448.         %--------------
  449.         [wname,Lev_Anal] = wmemtool('rmb',win_tool,n_membloc1, ...
  450.                                           ind_wave,ind_NB_lev);        
  451.         valTHR = utthrw2d('get',win_tool,'allvalthr');
  452.         valTHR = valTHR(:,1:Lev_Anal);
  453.         [name,ext] = strtok(filename,'.');
  454.         if isempty(ext) | isequal(ext,'.')
  455.             ext = '.mat'; filename = [name ext];
  456.         end
  457.         saveStr = {'X','map','wname','valTHR'};
  458.         wwaiting('off',win_tool);        
  459.         try
  460.           save([pathname filename],saveStr{:});
  461.         catch
  462.           errargt(mfilename,'Save FAILED !','msg');
  463.         end
  464.    
  465.         % End waiting.
  466.         %-------------
  467.         wwaiting('off',win_tool);
  469.     case 'decompose'
  470.         % Compute decomposition and plot.
  471.         %--------------------------------
  472.         wwaiting('msg',win_tool,'Wait ... computing');
  474.         % Get handles analysis parameters.
  475.         %---------------------------------
  476.         [wname,Lev_Anal] = cbanapar('get',win_tool,'wav','lev');
  477.         hdl_ori = utthrw2d('get',win_tool,'handleORI');
  478.         img_ori = get(hdl_ori,'Cdata');
  479.         siz_ori = size(img_ori);
  480.         pow = 2^Lev_Anal;
  481.         if any(rem(siz_ori,pow))
  482.             siz_sug = ceil(siz_ori/pow)*pow;
  483.             wwaiting('off',win_tool);
  484.             oriStr = ['(' int2str(siz_ori(1)) ',' int2str(siz_ori(2)) ')'];
  485.             sugStr = ['(' int2str(siz_sug(1)) ',' int2str(siz_sug(2)) ')'];
  486.             msg = strvcat(...
  487.             ['The level of decomposition ' int2str(Lev_Anal)],...
  488.             ['and the size of the image ' oriStr],...
  489.             'are not compatible.',...
  490.             ['Suggested size: ' sugStr],...
  491.             '(see Image Extension Tool)', ...
  492.             ' ', ...
  493.             ['2^Level has to divide the size of the image.'] ...
  494.             );
  495.             errargt(mfilename,msg,'msg');
  496.             return
  497.         end        
  500.         % Clean HDLG.
  501.         %---------------
  502.         utthrw2d('clean_thr',win_tool);
  504.         % Get Handles.
  505.         %--------------------------------------------------------
  506.         % ax_hdl = {axe_left_1 axe_right_1 axe_medium  axe_hist}
  507.         %--------------------------------------------------------
  508.         ax_hdl = wmemtool('rmb',win_tool,n_membloc2,ind_axe_hdl);
  510.         % Clean axes.
  511.         %------------
  512.         axe_hist = ax_hdl{4};
  513.         axes2clean = [ax_hdl{2};ax_hdl{3};axe_hist(:)];
  514.         obj2del = [findobj(axes2clean,'type','line');   ...
  515.                    findobj(axes2clean,'type','image');  ...
  516.                    findobj(axes2clean,'type','patch')   ...
  517.                    ];
  518.         delete(obj2del)
  519.         set(axe_hist(:),'Nextplot','add')
  521.         % Setting prog status.
  522.         %----------------------
  523.         wmemtool('wmb',win_tool,n_membloc1, ...
  524.                        ind_status,1,ind_NB_lev,Lev_Anal,ind_wave,wname);
  526.         % Clean reset axes limits.
  527.         %-------------------------
  528.         set([ax_hdl{1:2}],'Xlim',[1 siz_ori(2)],'Ylim',[1 siz_ori(1)])
  530.         % Compute Decomposition.
  531.         %------------------------
  532.         wDEC = swt2(img_ori,Lev_Anal,wname);
  533.         wmemtool('wmb',win_tool,n_membloc3,ind_coefs,wDEC);
  535.         % Initializing by level threshold.
  536.         %---------------------------------
  537.         maxTHR = zeros(3,Lev_Anal);
  538.         valTHR = sw2dtool('compute_LVL_THR',win_tool);
  539.         permDir = [1 3 2];
  540.         index   = zeros(3,Lev_Anal);
  541.         for k = 1:Lev_Anal
  542.           for d=1:3 , index(d,k) = (permDir(d)-1)*Lev_Anal+k; end
  543.         end
  544.         for k=1:Lev_Anal
  545.             for d=1:3
  546.                 ind = index(d,k);
  547.                 curr_img = wDEC(:,:,ind);
  548.                 curr_img = curr_img(:);
  549.                 tmp = max(abs(curr_img));
  550.                 if tmp<eps , maxTHR(d,k) = 1; else , maxTHR(d,k) = 1.1*tmp; end
  551.             end
  552.         end
  553.         valTHR = min(maxTHR,valTHR);
  555.         % Displaying details coefficients histograms.
  556.         %--------------------------------------------
  557.         dirDef   = 1;
  558.         fontsize = wmachdep('fontsize','normal');
  559.         col_det  = wtbutils('colors','det',Lev_Anal);
  560.         nb_bins  = 50;
  561.         axeXColor = get(win_tool,'DefaultAxesXColor');        
  562.         for level = 1:Lev_Anal
  563.             for direct=1:3
  564.                 ind = index(direct,level);
  565.                 curr_img = wDEC(:,:,ind);
  566.                 curr_img = curr_img(:);
  567.                 axeAct   = axe_hist(direct,level);
  568.                 axes(axeAct);
  569.                 curr_color = col_det(level,:);
  570.                 his        = wgethist(curr_img,nb_bins);
  571.                 his(2,:)   = his(2,:)/length(curr_img);
  572.                 hdl_hist   = wplothis(axeAct,his,curr_color);
  573.                 if level>1
  574.                     wxlabel('','Parent',axeAct);
  575.                 else
  576.                     wxlabel([deblank(str_dir_det(direct,:)) ' Details'],...
  577.                             'color',axeXColor,'Parent',axeAct);
  578.                 end
  579.                 thr_val = valTHR(direct,level);
  580.                 thr_max = maxTHR(direct,level);
  581.                 ylim    = get(axeAct,'Ylim');
  582.                 utthrw2d('plot_dec',win_tool,dirDef, ...
  583.                           {thr_max,thr_val,ylim,direct,level,axeAct})
  584.                 xmax = 1.1*max([thr_max, max(abs(his(1,:)))]);
  585.                 set(axeAct,'Xlim',[-xmax xmax]);
  586.                 set(findall(axeAct),'Visible','on');
  587.             end
  588.         end
  589.         drawnow
  591.         % Initialization of denoising structure.
  592.         %----------------------------------------
  593.         utthrw2d('set',win_tool,'valthr',valTHR,'maxthr',maxTHR);
  595.         % Dynvtool Attachement.
  596.         %---------------------
  597.         axe_cmd = [ax_hdl{1:2}];
  598.         dynvtool('init',win_tool,[],axe_cmd,[],[1 1],'','','','int')
  600.         % Enabling HDLG.
  601.         %---------------
  602.         sw2dtool('enable',win_tool,'dec','on');
  604.         % End waiting.
  605.         %-------------
  606.         wwaiting('off',win_tool);
  608.     case 'denoise'
  609.         % Compute decomposition and plot.
  610.         %--------------------------------
  611.         wwaiting('msg',win_tool,'Wait ... computing');
  613.         % Diseable De-noising Tool.
  614.         %---------------------------
  615.         utthrw2d('enable',win_tool,'off');
  617.         % Get decomposition.
  618.         %-------------------
  619.         [wname,Lev_Anal] = wmemtool('rmb',win_tool,n_membloc1, ...
  620.                                           ind_wave,ind_NB_lev);        
  621.         wDEC = wmemtool('rmb',win_tool,n_membloc3,ind_coefs);
  622.         valTHR = utthrw2d('get',win_tool,'allvalthr');
  623.         valTHR = valTHR(:,1:Lev_Anal);
  624.         permDir = [1 3 2];
  625.         for level = 1:Lev_Anal
  626.             for kk=1:3
  627.                ind = (permDir(kk)-1)*Lev_Anal+level;
  628.                thr = valTHR(kk,level);
  629.                wDEC(:,:,ind) = wthresh(wDEC(:,:,ind),'s',thr);
  630.             end
  631.         end
  633.         % Plotting de-noised Image.
  634.         %---------------------------
  635.         img_den = iswt2(wDEC,wname);
  636.         axe_hdl = wmemtool('rmb',win_tool,n_membloc2,ind_axe_hdl);
  637.         axe_ori = axe_hdl{1};
  638.         axe_den = axe_hdl{2};
  639.         axe_dif = axe_hdl{3};
  640.         hdl_den = image(img_den,'Parent',axe_den);
  641.         xylim   = get(axe_ori,{'Xlim','Ylim'});
  642.         set(axe_den,'Xlim',xylim{1},'Ylim',xylim{2},'Visible','on');
  643.         wtitle('De-Noised image (DI)','Parent',axe_den);
  644.         hdl_ori  = utthrw2d('get',win_tool,'handleORI');
  645.         img_dif  = get(hdl_ori,'Cdata')-img_den;
  646.         img_dif  = img_dif(:);
  647.         nb_bins  = 80;
  648.         color    = wtbutils('colors','sw2d','histRES');
  649.         his      = wgethist(img_dif,nb_bins);
  650.         his(2,:) = his(2,:)/length(img_dif);
  651.         hdl_hist = wplothis(axe_dif,his,color);
  652.         wtitle('Histogram of residuals: (I) - (DI)','Parent',axe_dif);
  654.         % Memory blocks update.
  655.         %----------------------
  656.         utthrw2d('set',win_tool,'handleTHR',hdl_den);
  657.         wmemtool('wmb',win_tool,n_membloc1,ind_status,1);
  659.         % Dynvtool Attachement.
  660.         %---------------------
  661.         dynvtool('ini_his',win_tool,0);
  662.         dynvtool('put',win_tool)
  664.         % Enabling HDLG.
  665.         %---------------
  666.         utthrw2d('enable',win_tool,'on');
  667.         sw2dtool('enable',win_tool,'den','on');
  669.         % End waiting.
  670.         %---------------
  671.         wwaiting('off',win_tool);
  673.     case 'update_level'
  674.         pop_lev = varargin{2}(1);
  675.         lev_New  = get(pop_lev,'value');
  676.         [status,lev_Anal] = wmemtool('rmb',win_tool,n_membloc1,...
  677.             ind_status,ind_NB_lev);
  678.         utthrw2d('set',win_tool,'position',{1,lev_New});
  679.         sw2dtool('set_axes',win_tool);
  681.         vis_InAxes = 'off';
  682.         if isequal(lev_New,lev_Anal)
  683.             switch status
  684.                 case -1 , 
  685.                     sw2dtool('enable',win_tool,'ini');
  686.                 case  0 ,
  687.                     sw2dtool('enable',win_tool,'dec');
  688.                 case  1 ,
  689.                     vis_InAxes = 'on';
  690.                     sw2dtool('enable',win_tool,'dec','on');
  691.                     sw2dtool('enable',win_tool,'den','on');
  692.             end
  693.         else
  694.             sw2dtool('enable',win_tool,'ini');
  695.         end
  696.         
  697.         % Get Handles.
  698.         %--------------------------------------------------------
  699.         % ax_hdl = {axe_left_1 axe_right_1 axe_medium  axe_hist}
  700.         %--------------------------------------------------------
  701.         ax_hdl  = wmemtool('rmb',win_tool,n_membloc2,ind_axe_hdl);
  702.         axe_Off = [ax_hdl{2};ax_hdl{3}];
  703.         obj_Off = [findobj(axe_Off,'type','line');  ...
  704.                    findobj(axe_Off,'type','patch'); ...
  705.                    findobj(axe_Off,'type','image')];
  706.         set(obj_Off,'Visible',vis_InAxes);
  707.                 
  708.     case 'compute_LVL_THR'
  709.         [numMeth,meth,alfa,sorh] = utthrw2d('get_LVL_par',win_tool);
  710.         wDEC = wmemtool('rmb',win_tool,n_membloc3,ind_coefs);
  711.         varargout{1} = wthrmngr('sw2ddenoLVL',meth,wDEC,alfa);
  713.     case 'update_LVL_meth'
  714.         sw2dtool('clear_GRAPHICS',win_tool);
  715.         valTHR = sw2dtool('compute_LVL_THR',win_tool);
  716.         utthrw2d('update_LVL_meth',win_tool,valTHR);
  718.     case 'clear_GRAPHICS'
  719.         status = wmemtool('rmb',win_tool,n_membloc1,ind_status);
  720.         if status<1 , return; end
  722.         % Diseable Toggle and Menus.
  723.         %---------------------------
  724.         sw2dtool('enable',win_tool,'den','off');
  726.         % Get Handles.
  727.         %-------------
  728.         axe_hdl = wmemtool('rmb',win_tool,n_membloc2,ind_axe_hdl);
  729.         axe_Off = [axe_hdl{2};axe_hdl{3}];
  730.         obj_Off = [findobj(axe_Off,'type','line');  ...
  731.                    findobj(axe_Off,'type','patch'); ...
  732.                    findobj(axe_Off,'type','image')];
  733.         set(obj_Off,'Visible','off');
  735.     case 'enable'
  736.         type  = varargin{2};
  737.         m_sav = wmemtool('rmb',win_tool,n_membloc1,ind_sav_menu);
  738.         switch type 
  739.           case 'ini'
  740.             set(m_sav,'Enable','off');
  741.             pus_dec = wmemtool('rmb',win_tool,n_membloc2,ind_pus_dec);
  742.             col_hdl = utcolmap('handles',win_tool,'act','true');
  743.             utthrw2d('status',win_tool,'off');
  744.             set([pus_dec;col_hdl],'Enable','on');
  746.           case 'dec'
  747.             nb_lev = wmemtool('rmb',win_tool,n_membloc1,ind_NB_lev);
  748.             set(m_sav,'Enable','off');
  749.             utthrw2d('status',win_tool,'on');
  750.             utthrw2d('enable',win_tool,'on',[1:nb_lev]);
  751.             wmemtool('wmb',win_tool,n_membloc1,ind_status,0);
  753.           case 'den'
  754.             enaVal = varargin{3};
  755.             set(m_sav,'Enable',enaVal);
  756.             utthrw2d('enable_tog_res',win_tool,enaVal);
  757.             if strncmpi(enaVal,'on',2) , status = 1; else , status = 0; end
  758.             wmemtool('wmb',win_tool,n_membloc1,ind_status,status);
  759.         end
  761.     case 'set_axes'
  762.         %*************************************************************%
  763.         %** OPTION = 'set_axes' - Set axes positions and visibility **%
  764.         %*************************************************************%
  765.         pos_win = get(win_tool,'Position');
  766.         [ax_hdl,Pos_Graphic_Area] = wmemtool('rmb',win_tool,n_membloc2,...
  767.                                                    ind_axe_hdl,ind_gra_area);
  768.         [wname,NB_lev] = cbanapar('get',win_tool,'wav','lev');
  770.         % Hide axes
  771.         %-----------
  772.         ax_2_clean  = [ax_hdl{1:3},ax_hdl{4}(:)'];
  773.         obj_in_axes = findobj(ax_2_clean);
  774.         set(obj_in_axes,'Visible','off');
  776.         % Plots.
  777.         %--------------------------------------------------------
  778.         % ax_hdl = {axe_left_1 axe_right_1 axe_medium  axe_hist} 
  779.         %--------------------------------------------------------
  780.         NBaxes  = length(ax_2_clean);
  781.         ax_l_1  = ax_hdl{1};
  782.         ax_r_1  = ax_hdl{2};
  783.         ax_med  = ax_hdl{3};
  784.         ax_hist = ax_hdl{4};
  786.         % General graphical parameters initialization.
  787.         %---------------------------------------------
  788.         ecy_up    = 0.06*pos_win(4);
  789.         ecy_mid_1 = 0.08*pos_win(4);
  790.         ecy_mid_2 = ecy_up;
  791.         ecy_down  = ecy_up;
  792.         ecy_det   = (0.04*pos_win(4))/1.4;
  793.         h_gra_rem = Pos_Graphic_Area(4);
  794.         h_min     = h_gra_rem/4;
  795.         h_max     = h_gra_rem/3;
  796.         h_image   = (h_min*(NB_lev-1)+h_max*(max_lev_anal-NB_lev))/(max_lev_anal-1);
  797.         h_histo   = h_gra_rem/8;
  798.         h_gra_rem = h_gra_rem-h_histo-h_image-ecy_up-ecy_mid_1-ecy_mid_2;
  799.         h_detail  = (h_gra_rem-ecy_down-(NB_lev-1)*ecy_det)/NB_lev;
  800.         y_low_ini = 1;
  802.         % Building data axes.
  803.         %--------------------
  804.         y_low_ini = y_low_ini-h_image-ecy_up;
  805.         pos_axes  = get(ax_l_1,'Position');
  806.         pos_axes([2 4]) = [y_low_ini h_image];
  807.         set(ax_l_1,'Position',pos_axes);
  808.         axe_vis   = [ax_l_1];
  810.         pos_axes = get(ax_r_1,'Position');
  811.         pos_axes([2 4]) = [y_low_ini h_image];
  812.         set(ax_r_1,'Position',pos_axes)
  813.         axe_vis  = [axe_vis ax_r_1];
  815.         y_low_ini = y_low_ini-h_histo-ecy_mid_1;
  816.         pos_axes  = get(ax_med,'Position');
  817.         pos_axes([2 4]) = [y_low_ini h_histo];
  818.         set(ax_med,'Position',pos_axes)
  819.         axe_vis   = [axe_vis ax_med];
  820.         y_low_ini = y_low_ini-ecy_mid_2+ecy_det;
  822.         % Position for histograms axes.
  823.         %------------------------------   
  824.         axeXColor = get(win_tool,'DefaultAxesXColor');        
  825.         pos_y    = [y_low_ini , h_detail];
  826.         for level = NB_lev:-1:1
  827.             pos_y(1) = pos_y(1)-h_detail-ecy_det;
  828.             for direct=1:3
  829.                 ax_act   = ax_hist(direct,level);
  830.                 pos_axes = get(ax_act,'Position');           
  831.                 pos_axes([2 4]) = pos_y;
  832.                 set(ax_act,'Position',pos_axes);
  833.                 axe_vis = [axe_vis ax_act];
  834.                 if level>1
  835.                     wxlabel('','Parent',ax_act);
  836.                 else
  837.                     wxlabel([deblank(str_dir_det(direct,:)) ' Details'],...
  838.                             'color',axeXColor,...
  839.                             'Parent',ax_act);
  840.                 end
  841.             end           
  842.         end
  844.         % Setting axes visibility & title.
  845.         %---------------------------------
  846.         obj_in_axes_vis = findobj(axe_vis);
  847.         set(obj_in_axes_vis,'Visible','on');
  848.         wtitle('Original Image (I)','Parent',ax_l_1);
  849.         wtitle('Histogram of residuals: (I) - (DI)','Parent',ax_med);
  850.         wtitle('De-Noised Image (DI)','Parent',ax_r_1);
  851.                     
  852.     case 'close'
  854.     otherwise
  855.         errargt(mfilename,'Unknown Option','msg');
  856.         error('*');
  858. end