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

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

波变换

开发平台：

Matlab

dcmdwpck.m：源码内容
  1. function varargout = dcmdwpck(varargin)
  2. %DCMDWPCK Demonstrates wavelet packet tools in the Wavelet Toolbox. 
  3. %
  4. % This is a slideshow file for use with wshowdrv.m
  5. % To see it run, type 'wshowdrv dcmdwpck', 
  6. %
  7. % See also DDENCMP, WCODEMAT, WDEN, WDENCMP, WNOISE, WPDENCMP.
  9. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
  10. %   Last Revision: 14-May-2003.
  11. %   Copyright 1995-2004 The MathWorks, Inc.
  12. % $Revision: 1.14.4.2 $
  14. % Initialization and Local functions if necessary.
  15. if nargin>0
  16. action = varargin{1};
  17. switch action
  18.       case 'addHelp'
  19. % Add Help Item.
  20. %---------------
  21.         hdlFig = varargin{2};
  22. wfighelp('addHelpItem',hdlFig,'Wavelet Packets','WP_PACKETS');
  23. wfighelp('addHelpItem',hdlFig,'Compression Procedure','COMP_PROCEDURE');
  24.         wfighelp('addHelpItem',hdlFig,'De-noising Procedure','DENO_PROCEDURE');
  25.         
  26.       case 'auto'    , wshowdrv('#autoMode',mfilename,'close');
  27.   case 'gr_auto' , wshowdrv('#gr_autoMode',mfilename,'close');
  29.   case 'getFigParam'
  30. figName  = 'Wavelet packets demo';
  31. showType = 'mix7';
  32. varargout = {figName,showType};
  33. end
  34. return
  35. end
  37. if nargout<1,
  38.   wshowdrv(mfilename)
  39. else
  40.   idx = 0;
  41.   %========== Slide 1 ==========
  42.   idx = idx+1;
  43.   slide(idx).code = {
  44.   'figHandle = gcf;';
  45.   'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow; ',   
  46.   '' };
  48.   slide(idx).text = {
  49.   '',
  50.   ' Press the "Start" button to see a demonstration of',
  51.   ' wavelet packet tools in the Wavelet Toolbox.',
  52.   '',
  53.   ' This demo uses Wavelet Toolbox functions.',
  54.   ''};
  56.   %========== Slide 2 ==========
  57.   idx = idx+1;
  58.   slide(idx).code = {
  59.   'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow; ',   
  60.   ''};
  62.   slide(idx).text = {
  63.   ' % For wavelet packets analysis demos',
  64.   ' % see the GUI mode demos.',
  65.   ' % Let us give some examples of de-noising and compression',
  66.   ' % using wavelet packets.',
  67.   ''};
  69.   %========== Slide 3 ==========
  70.   idx = idx+1;
  71.   slide(idx).code = {
  72.   ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  73.   'if (idxPREV<idxSlide) | (idxPREV>idxSlide+1)',
  74.   '   init = 2055615866; randn(''seed'',init);', 
  75.   '   [xref,x] = wnoise(5,11,7,init);',
  76.   '   wshowdrv(''#set_axes'',figHandle,[4,1]);'
  77.   '   h = wshowdrv(''#get_axes'',figHandle);',
  78.   '   set(h(1),''Visible'',''on'');',    
  79.   '   axes(h(1)); plot(xref,''r''); title(''Original signal'')',
  80.   '   axis([1 length(xref) min(xref)-eps max(xref)+eps]);',
  81.   '   set(h(2),''Visible'',''on'');',
  82.   '   axes(h(2)); plot(x,''g''), title(''Noisy Signal'')',
  83.   '   axis([1 length(x) min(x)-eps max(x)+eps]);',
  84.   '   p3 = get(h(3),''Position''); p4 = get(h(4),''Position'');',
  85.   '   p5 = [p4(1:3), p3(2)+p3(4)-p4(2)];',
  86.   '   try , delete(h(5)); end',
  87.   '   h(5) = axes(''Position'',p5,''Visible'',''Off'');',
  88.   '   wtbxappdata(''set'',figHandle,''local_AXES'',h);'
  89.   'else'
  90.   '   h = wtbxappdata(''get'',figHandle,''local_AXES'');'
  91.   '   set(findobj(h(3)),''Visible'',''off'');',
  92.   'end'
  93.   ''};
  95.   slide(idx).text = {
  96.   '% Generate heavy sine and a noisy version of it.',
  97.   ' ',
  98.   '        init = 2055615866; randn(''seed'',init);',
  99.   '        [xref,x] = wnoise(5,11,7,init);', 
  100.   ''};
  102.   slide(idx).info = 'wnoise';
  104.   %========== Slide 4 ==========
  105.   idx = idx+1;
  106.   slide(idx).code = {
  107.   ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'], 
  108.   'if idxPREV<idxSlide',
  109.   '   xwpd = wtbxappdata(''get'',h(3),''xwpd'');',
  110.   '   if isempty(xwpd)',
  111.   '      n = length(x); thr = sqrt(2*log(n*log(n)/log(2)));'
  112.   '      xwpd = wpdencmp(x,''s'',4,''sym4'',''sure'',thr,1);',
  113.   '      wtbxappdata(''set'',h(3),''xwpd'',xwpd);',
  114.   '   end',
  115.   '   set(h(3),''Visible'',''on'');',
  116.   '   axes(h(3)); plot(xwpd); title(''Denoised Signal using Wavelet Packets'')',
  117.   '   axis([1 length(xwpd) min(xwpd)-eps max(xwpd)+eps]);',
  118.   'else'
  119.   '   set(findobj(h(4)),''Visible'',''off'');',
  120.   'end'
  121.   ''};
  123.   slide(idx).text = {
  124.   ' % Use wpdencmp for signal de-noising.',
  125.   ' ',
  126.   '        n = length(x);',
  127.   '        thr = sqrt(2*log(n*log(n)/log(2)));',
  128.   '        xwpd = wpdencmp(x,''s'',4,''sym4'',''sure'',thr,1);',
  129.   ''};
  131.   slide(idx).info = 'wpdencmp';
  133.   %========== Slide 5 ==========
  134.   idx = idx+1;
  135.   slide(idx).code = {
  136.   ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  137.   'if idxPREV<idxSlide',
  138.   '   xwd = wtbxappdata(''get'',h(4),''xwd'');',
  139.   '   if isempty(xwd)',
  140.   '      xwd = wden(x,''rigrsure'',''s'',''one'',4,''sym4'');',
  141.   '      wtbxappdata(''set'',h(4),''xwd'',xwd);',
  142.   '   end',
  143.   '   set(h(4),''Visible'',''on'');',
  144.   '   axes(h(4)); plot(xwd); title(''Denoised Signal using Wavelets'')',
  145.   '   axis([1 length(xwd) min(xwd)-eps max(xwd)+eps]);',
  146.   'else',
  147.   '   set(findobj(h(5)),''Visible'',''Off'');',
  148.   '   set(findobj(h(3:4)),''Visible'',''On'');',
  149.   '   xlim = get(h(3),''Xlim'');',
  150.   '   set(h([1,2,5]),''Xlim'',xlim);',
  151.   'end',
  152.   ''};
  154.   slide(idx).text = {
  155.   ' % Compare with wavelets based denoising result.',
  156.   ' ',
  157.   '        xwd = wden(x,''rigrsure'',''s'',''one'',4,''sym4'');',
  158.   ''};
  160.   slide(idx).info = 'wden';
  162.   %========== Slide 6 ==========
  163.   idx = idx+1;
  164.   slide(idx).code = {
  165.   'xwpd = wtbxappdata(''get'',h(3),''xwpd'');',
  166.   'xwd  = wtbxappdata(''get'',h(4),''xwd'');',
  167.   'set(findobj(h(3:4)),''Visible'',''Off'')',
  168.   'set(h(1:2),''Xticklabel'',[]);',
  169.   'set(h(5),''Visible'',''on'');',
  170.   'ind = 1:800;',
  171.   'axes(h(5)); plot(ind,xwd(ind),''y'',ind,xwpd(ind),''m'')',
  172.   'title(''Denoised Signals : Wavelet Packets (m) , Wavelets (y)'')',
  173.   'set(h([1,2,5]),''Xlim'',[1 800]);',
  174.   ''};
  176.   slide(idx).text = {
  177.   ' % Zoom on the initial part of the signal,' ,
  178.   ' % the noisy signal and the two denoised signals.',
  179.   ''};
  181.   %========== Slide 7 ==========
  182.   idx = idx+1;
  183.   slide(idx).code = {
  184.   'load detfingr',
  185.   '[r,c] = size(X);  X = X(1:2:r,1:2:c);', 
  186.   'sm = size(map,1); colormap(pink(sm))',
  187.   ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  188.   'if (idxPREV<idxSlide) | (idxPREV>idxSlide+1)',
  189.   '   ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow; ',
  190.   '   tmp = subplot(1,2,1); image(wcodemat(X,sm)), title(''Original Image'')',
  191.   '   axis(''square'');',
  192.   '   set(tmp,''Xtick'',[],''Ytick'',[]);',
  193.   'else',
  194.   '   set(findobj(subplot(1,2,2)),''Visible'',''Off'')',
  195.   'end'
  196.   ''};
  198.   slide(idx).text = {
  199.   ' load detfingr ',
  200.   ' % X contains the loaded image.',
  201.   ' ',
  202.   '  [r,c] = size(X);', 
  203.   '  X = X(1:2:r,1:2:c);    % To reduce the size of X.',
  204.   ''};
  206.   slide(idx).idxPrev = 3;
  208.   %========== Slide 8 ==========
  209.   idx = idx+1;
  210.   slide(idx).code = {
  211.   'set(findobj(subplot(1,2,2)),''Visible'',''Off'')',
  212.   ''};
  214.   slide(idx).text = {
  215.   ' % Use wpdencmp for image compression.',
  216.   '        sorh = ''h''; lev = 3;',
  217.   '        crit = ''shannon'';',
  218.   '        thr = 30;',
  219.   '        keepapp = 1;',
  220.   '        [xd,t,perf0,perfl2] = ... ',
  221.   '                wpdencmp(X,sorh,lev,''bior6.8'',crit,thr,keepapp);',
  222.   ''};
  224.   slide(idx).info = 'wpdencmp';
  226.   %========== Slide 9 ==========
  227.   idx = idx+1;
  228.   slide(idx).code = {
  229.   'ax = subplot(1,2,1);',
  230.   'denFLAG_1 = wtbxappdata(''get'',ax,''denFLAG_1'');',
  231.   'if isempty(denFLAG_1)',
  232.   '   lev = 3; sorh = ''h''; crit = ''shannon''; thr = 30; keepapp = 1;',
  233.   '   [xd,t,perf0,perfl2] = wpdencmp(X,sorh,lev,''bior6.8'',crit,thr,keepapp);',
  234.   '   tmp = subplot(1,2,2); image(wcodemat(xd,sm)), title(''Compressed Image'')',
  235.   '   N2Str = [''norm cfs. rec.: '',num2str(perfl2,''%5.2f''),'' %''];',
  236.   '   N0Str = [''nul cfs: '',num2str(perf0,''%5.2f''),'' %''];',
  237.   '   strxlab = strvcat(N2Str,N0Str); xlabel(strxlab);',
  238.   '   axis(''square'');',
  239.   '   set(tmp,''Xtick'',[],''Ytick'',[]);',
  240.   '   wtbxappdata(''set'',ax,''denFLAG_1'',1);',
  241.   'else',
  242.   '  set(findobj(subplot(1,2,2)),''Visible'',''On'')',
  243.   'end'
  244.   ''};
  246.   slide(idx).text = slide(idx-1).text;
  248.   slide(idx).info = 'wpdencmp';
  250.   %========== Slide 10 ==========
  251.   idx = idx+1;
  252.   slide(idx).code = {
  253.   'load  sinsin;  X = X(1:64,1:64);',
  254.   'sm = size(map,1); colormap(pink(sm))',
  255.   ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  256.   'if idxPREV<idxSlide',
  257.   '   ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow; ',   
  258.   '   tmp = subplot(2,2,1); image(wcodemat(X,sm)), title(''Original Image'')',
  259.   '   set(tmp,''Xtick'',[],''Ytick'',[]);',
  260.   'else'
  261.   '   set(findobj(subplot(2,2,2)),''Visible'',''Off'')',
  262.   'end'
  263.   ''};
  265.   slide(idx).text = {
  266.   ' % Load original image.',
  267.   ' ',
  268.   '        load  sinsin',
  269.   '        % X contains the loaded image.',
  270.   ' ',
  271.   '        X = X(1:64,1:64); ',
  272.   ''};
  274.   slide(idx).idxPrev = 7;
  276.   %========== Slide 11 ==========
  277.   idx = idx+1;
  278.   slide(idx).code = {
  279.   'init = 2055615866; randn(''seed'',init);',
  280.   'x = X+15*randn(size(X));',
  281.   ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  282.   'if idxPREV<idxSlide',
  283.   '   tmp = subplot(2,2,2); image(wcodemat(x,sm)), title(''Noisy Image'')',
  284.   '   set(tmp,''Xtick'',[],''Ytick'',[]);',
  285.   'else'
  286.   '   set(findobj(subplot(2,2,3)),''Visible'',''Off'')',
  287.   'end'
  288.   ''};
  289.   
  290.   slide(idx).text = {
  291.   ' % Compute and plot the Noisy Image.',
  292.   ' ',
  293.   '        init = 2055615866; randn(''seed'',init);',
  294.   '        x = X+15*randn(size(X));',
  295.   ''};
  296.   
  297.   %========== Slide 12 ==========
  298.   idx = idx+1;
  299.   slide(idx).code = {
  300.   'ax = subplot(2,2,1);',
  301.   'denFLAG_2 = wtbxappdata(''get'',ax,''denFLAG_2'');',
  302.   'if isempty(denFLAG_2)',
  303.   '   [thr,sorh,keepapp,crit] = ddencmp(''den'',''wp'',x);',
  304.   '   n  = length(x(:)); thr = sqrt(2*log(n*log(n)/log(2)))*15;',
  305.   '   xd = wpdencmp(x,sorh,3,''sym4'',crit,thr,keepapp);',
  306.   '   tmp = subplot(2,2,3); image(wcodemat(xd,sm)), title(''Denoised Image'')',
  307.   '   set(tmp,''Xtick'',[],''Ytick'',[]);',
  308.   '   wtbxappdata(''set'',ax,''denFLAG_2'',1);',
  309.   'else',
  310.   '   set(findobj(subplot(2,2,3)),''Visible'',''On'')',
  311.   'end'
  312.   ''};
  313.   
  314.   slide(idx).text = {
  315.   ' % Compute and plot the Denoised Image.',
  316.   ' ',
  317.   ''};
  319.   slide(idx).info = 'wpdencmp';
  321.   varargout{1} = slide;
  322. end