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

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

波变换

开发平台：

Matlab

dcmdextm.m：源码内容
  1. function varargout = dcmdextm(varargin)
  2. %DCMDEXTM Demonstrates border distortions tools in the Wavelet Toolbox. 
  3. %
  4. % This is a slideshow file for use with wshowdrv.m
  5. % To see it run, type 'wshowdrv dcmdextm', 
  6. %
  7. % See also DWTMODE, WAVEDEC, WAVEDEC2, WRCOEF2. 
  9. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
  10. %   Last Revision: 07-May-2003.
  11. %   Copyright 1995-2004 The MathWorks, Inc.
  12. %   $Revision: 1.14.4.2 $  $Date: 2004/03/15 22:36:59 $
  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,'Border distortion','BORDER_DIST');
  23.         
  24.   case 'auto'    , wshowdrv('#autoMode',mfilename,'close');
  25.   case 'gr_auto' , wshowdrv('#gr_autoMode',mfilename,'close');
  27.   case 'getFigParam'
  28. figName  = 'Border distortion';
  29. showType = 'mix6';
  30. varargout = {figName,showType};
  32.   case 'localPROC_1'
  33. [figHandle,modeDWT,axeHDL,strTIT] = deal(varargin{2:end});
  34. dataSTRUCT = wtbxappdata('get',figHandle,'dataSTRUCT');
  35. [x,lev,wname,nbcol] = deal(dataSTRUCT{:});
  36. dwtmode(modeDWT,'silent');
  37. [c,l] = wavedec(x,lev,wname);
  38. len = length(x);
  39. cfd = zeros(lev,len);
  40. for k = 1:lev
  41. d = detcoef(c,l,k);
  42. d = d(:)';
  43. d = d(ones(1,2^k),:);
  44. cfd(k,:) = wkeep1(d(:)',len);
  45. end
  46. cfd =  cfd(:);
  47. I = find(abs(cfd)<sqrt(eps));
  48. cfd(I) = zeros(size(I));
  49. cfd    = reshape(cfd,lev,len);
  50. imgCFS = wcodemat(cfd,nbcol);
  51. tics = 1:lev; labs = int2str([1:lev]');
  52. axes(axeHDL); image(imgCFS);
  53. set(axeHDL,'YTicklabelMode','manual','YTick',tics, ...
  54. 'YTicklabel',labs,'YDir','normal','Box','On');
  55. ylabel('level');
  56. title(strTIT);
  57. colormap(pink(nbcol));
  59.   case 'localPROC_2'
  60. [figHandle,modeDWT,axeHDL,strTIT] = deal(varargin{2:end});
  61. dataSTRUCT = wtbxappdata('get',figHandle,'dataSTRUCT');
  62. [X,lev,wname,nbcol] = deal(dataSTRUCT{:});
  63. dwtmode(modeDWT,'silent');
  64. [c,s] = wavedec2(X,lev,wname);
  65.     a = wrcoef2('a',c,s,wname,lev);
  66. imgCFS = wcodemat(a,nbcol);
  67. axes(axeHDL); image(imgCFS);
  68. set(axeHDL,'XTick',[],'YTick',[]);
  69. title(strTIT);
  71. end
  72. return
  73. end
  75. if nargout<1,
  76.   wshowdrv(mfilename)
  77. else
  78.   idx = 0;
  79.   %========== Slide 1 ==========
  80.   idx = idx+1;
  81.   slide(idx).code = {
  82.   'figHandle = gcf;',
  83.   'try',
  84.   '   wtbxappdata(''del'',figHandle,''local_AXES'');',
  85.   '   wtbxappdata(''del'',figHandle,''dataSTRUCT'');',
  86.   'end',
  87.   'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow; ',   
  88.   '' };
  90.   slide(idx).text = {
  91.   '',
  92.   ' Press the "Start" button to see a demonstration on',
  93.   ' border distortions for Discrete Wavelet Transform (DWT).',
  94.   '',
  95.   ' This demo uses Wavelet Toolbox functions.',
  96.   ''};
  97.   
  98.   %========== Slide 2 ==========
  99.   idx = idx+1;
  100.   slide(idx).code = {
  101.   ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  102.   'if idxPREV<idxSlide',
  103.   '   wshowdrv(''#set_axes'',figHandle,[4,1]);'
  104.   '   h = wshowdrv(''#get_axes'',figHandle);',
  105.   '   load nearbrk; signal = nearbrk; clear nearbrk',
  106.   '   lev = 5; wname = ''db2''; nbcol = 128;',
  107.   '   dataSTRUCT = {signal,lev,wname,nbcol};',
  108.   '   colormap(pink(nbcol));',
  109.   '   axes(h(1)); plot(signal,''r'');',
  110.   '   title(''Original Signal : two near first derivative discontinuities'');',
  111.   '   wtbxappdata(''set'',figHandle,''local_AXES'',h);'
  112.   '   wtbxappdata(''set'',figHandle,''dataSTRUCT'',dataSTRUCT);'
  113.   'else ',
  114.   '   h = wtbxappdata(''get'',figHandle,''local_AXES'');'
  115.   'end',
  116.   'set(findobj(h(2:4)),''Visible'',''Off'');',
  117.   'set(findobj(h(1)),''Visible'',''On'');',
  118.   '' };
  120.   slide(idx).text = {
  121. ' % Load signal.',
  122. '        load nearbrk;',
  123. '        x = nearbrk;',
  124. ''};
  126.   %========== Slide 3 ==========
  127.   idx = idx+1;
  128.   slide(idx).code = {
  129. ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  130. 'if idxPREV>idxSlide',
  131. '   set(findobj(h(2:4)),''Visible'',''Off'');',
  132. '   set(findobj(h(1)),''Visible'',''On'');',
  133. 'end',
  134. '' };
  136.   slide(idx).text = {
  137. ' % Set DWT signal extension mode to zero padding.',
  138. ' ',
  139. '        dwtmode(''zpd'')',
  140. ''};
  142.   slide(idx).info = 'dwtmode';
  144.   %========== Slide 4 ==========
  145.   idx = idx+1;
  146.   slide(idx).code = {
  147. ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  148. 'if idxPREV<idxSlide',
  149. '   modeDWT = ''zpd'';',
  150. '   strTIT = ''Coefficients using zero padding'';',
  151. [   mfilename ,'(''localPROC_1'',figHandle,modeDWT,h(2),strTIT);'],
  152. 'end',
  153. 'set(findobj(h(3:4)),''Visible'',''Off'');',
  154. 'set(findobj(h(1:2)),''Visible'',''On'');',
  155. '' };
  157.   slide(idx).text = {
  158. ' % Perform decomposition at level 5',
  159. ' % for the signal using db2 and',
  160. ' % expand discrete wavelet coefficients for plot.',
  161. '        lev = 5;',
  162. '        [c,l] = wavedec(x,lev,''db2'');',
  163. ''};
  164.  
  165.   slide(idx).info = 'wavedec';
  167.   %========== Slide 5 ==========
  168.   idx = idx+1;
  169.   slide(idx).code = {
  170. ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  171. 'if idxPREV<idxSlide',
  172. '   modeDWT = ''sym'';',
  173. '   strTIT = ''Coefficients using symmetrization padding'';',
  174. [   mfilename ,'(''localPROC_1'',figHandle,modeDWT,h(3),strTIT);'],
  175. 'end',
  176. 'set(findobj(h(4)),''Visible'',''Off'');',
  177. 'set(findobj(h(1:3)),''Visible'',''On'');',
  178. '' };
  180.   slide(idx).text = {
  181. ' % Set DWT signal extension mode to symmetrization.',
  182. '        dwtmode(''sym'')',
  183. ' % Perform the same two previous steps.',
  184. '        [c,l] = wavedec(x,lev,''db2'');',
  185. ''};
  186.  
  187.   slide(idx).info = 'wavedec';
  189.   %========== Slide 6 ==========
  190.   idx = idx+1;
  191.   slide(idx).code = {
  192. ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  193. 'if idxPREV<idxSlide',
  194. '   modeDWT = ''spd'';',
  195. '   strTIT = ''Coefficients using smooth padding'';',
  196. [   mfilename ,'(''localPROC_1'',figHandle,modeDWT,h(4),strTIT);'],
  197. 'else'
  198. ' dataSTRUCT = wtbxappdata(''get'',figHandle,''dataSTRUCT'');',
  199. '   nbcol = dataSTRUCT{4};',
  200. '   colormap(pink(nbcol));'
  201. 'end',
  202. 'set(findobj(h(1:4)),''Visible'',''On'');',
  203. '' };
  205.   slide(idx).text = {
  206. ' % Set DWT signal extension mode to smooth padding.',
  207. '        dwtmode(''spd'')',
  208. ' % Perform the same two previous steps.',
  209. '        [c,l] = wavedec(x,lev,''db2'');',
  210. ''};
  211.  
  212.   slide(idx).info = 'wavedec';
  214.   %========== Slide 7 ==========
  215.   idx = idx+1;
  216.   slide(idx).code = {
  217. ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  218. 'colormap(cool(2))',   
  219. '' };
  221.   slide(idx).text = {
  222. ' % Change colormap.',
  223. ' ',
  224. '        [colormap(cool(2))',
  225. ''};
  226.  
  227.   slide(idx).info = 'dwtmode';
  229.   %========== Slide 8 ==========
  230.   idx = idx+1;
  231.   slide(idx).code = {
  232.   ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  233.   'if idxPREV<idxSlide',
  234.   '   wshowdrv(''#set_axes'',figHandle,[2,2]);',
  235.   '   h = wshowdrv(''#get_axes'',figHandle);',
  236.   '   load geometry; [row,col] = size(X);',
  237.   '   lev = 3; wname = ''sym4''; nbcol = size(map,1);',
  238.   '   dataSTRUCT = {X,lev,wname,nbcol};',
  239.   '   colormap(pink(nbcol));',
  240.   '   axes(h(1)); image(wcodemat(X,nbcol));',
  241.   '   title(''Original image'');',
  242.   '   wtbxappdata(''set'',figHandle,''local_AXES'',h);'
  243.   '   wtbxappdata(''set'',figHandle,''dataSTRUCT'',dataSTRUCT);'
  244.   'else',
  245.   '   h = wtbxappdata(''get'',figHandle,''local_AXES'');'
  246.   'end',
  247.   'set(findobj(h(2:4)),''Visible'',''Off'');',
  248.   'set(findobj(h(1)),''Visible'',''On'');',
  249. '' };
  251.   slide(idx).text = {
  252. ' % Load artificial image.',
  253. '        load geometry;' ,
  254. ' % Set DWT signal extension mode to zero padding.',
  255. '        dwtmode(''zpd'')', 
  256. ''};
  257.  
  258.   slide(idx).info = 'dwtmode';
  260.   slide(idx).idxPrev = 1; 
  262.   %========== Slide 9 ==========
  263.   idx = idx+1;
  264.   slide(idx).code = {
  265. ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  266. 'if idxPREV<idxSlide',
  267. '   modeDWT = ''zpd'';',
  268. '   strTIT = ''app. 3 using zero padding'';',
  269. [   mfilename ,'(''localPROC_2'',figHandle,modeDWT,h(2),strTIT);'],
  270. 'end',
  271. 'set(findobj(h(3:4)),''Visible'',''Off'');',
  272. 'set(findobj(h(1:2)),''Visible'',''On'');',
  273. '' };
  275.   slide(idx).text = {
  276. ' % Perform decomposition at level 3',
  277. ' % of X using sym4, and reconstruct',
  278. ' % approximation of level 3.',
  279. '        lev = 3; [c,s] = wavedec2(X,lev,''sym4'');',
  280. '        a = wrcoef2(''a'',c,s,''sym4'',lev);',
  281. ''};
  282.  
  283.   slide(idx).info = 'wavedec2';
  285.   %========== Slide 10 ==========
  286.   idx = idx+1;
  287.   slide(idx).code = {
  288. ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  289. 'if idxPREV<idxSlide',
  290. '   modeDWT = ''sym'';',
  291. '   strTIT = ''app. 3 using symmetrization padding'';',
  292. [   mfilename ,'(''localPROC_2'',figHandle,modeDWT,h(3),strTIT);'],
  293. 'end',
  294. 'set(findobj(h(4)),''Visible'',''Off'');',
  295. 'set(findobj(h(1:3)),''Visible'',''On'');',
  296. '' };
  298.   slide(idx).text = {
  299. ' % Set DWT signal extension mode to symmetrization.',
  300. '        dwtmode(''sym'')',
  301. ' % Perform the same two previous steps.',
  302. '        [c,s] = wavedec2(X,lev,''sym4'');',
  303. '        a = wrcoef2(''a'',c,s,''sym4'',lev);',
  304. ''};
  305.  
  306.   slide(idx).info = 'wrcoef2';
  308.   %========== Slide 11 ==========
  309.   idx = idx+1;
  310.   slide(idx).code = {
  311. ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  312. 'if idxPREV<idxSlide',
  313. '   modeDWT = ''spd'';',
  314. '   strTIT = ''app. 3 using smooth padding'';',
  315. [   mfilename ,'(''localPROC_2'',figHandle,modeDWT,h(4),strTIT);'],
  316. 'end',
  317. 'set(findobj(h(1:4)),''Visible'',''On'');',
  318. '' };
  320.   slide(idx).text = {
  321. ' % Set DWT signal extension mode to smooth.',
  322. '        dwtmode(''spd'')',
  323. ' % Perform the same two previous steps.',
  324. '        [c,s] = wavedec2(X,lev,''sym4'');',
  325. '        a = wrcoef2(''a'',c,s,''sym4'',lev);',
  326. ''};
  327.  
  328.   slide(idx).info = 'wrcoef2';
  330.   varargout{1} = slide;
  331. end