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

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

波变换

开发平台：

Matlab

dcmddw2d.m：源码内容
  1. function varargout = dcmddw2d(varargin)
  2. %DCMDDW2D Demonstrates discrete 2-D wavelet tools in the Wavelet Toolbox. 
  3. %
  4. % This is a slideshow file for use with wshowdrv.m
  5. % To see it run, type 'wshowdrv dcmddw2d', 
  6. %
  7. % See also APPCOEF2, DETCOEF2, DWT2, IDWT2, UPCOEF2, 
  8. %          WAVEDEC2, WAVEREC2, WRCOEF2, WCODEMAT.
  10. %   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
  11. %   Last Revision: 18-Apr-2001.
  12. %   Copyright 1995-2002 The MathWorks, Inc.
  13. % $Revision: 1.15 $
  15. % Initialization and Local functions if necessary.
  16. if nargin>0
  17. action = varargin{1};
  18. switch action
  19.   case 'auto'    , wshowdrv('#autoMode',mfilename,'close');
  20.   case 'gr_auto' , wshowdrv('#gr_autoMode',mfilename,'close');
  22.   case 'getFigParam'
  23. figName  = 'Orthogonal Wavelets 2-D';
  24. showType = 'mix9';
  25. varargout = {figName,showType};
  27.   case 'initShowViewer'
  28. figHandle = varargin{2};
  29. slideData = get(figHandle,'Userdata');
  30. autoHndl = slideData.autoHndl;
  31. propAuto = get(autoHndl,{'Units','Position','BackGroundColor'});
  32. txtBkColor = get(slideData.slitxtHndl,'BackGroundColor');
  33. pos  = propAuto{2};
  34. hBtn = pos(4);
  35. pos(2) = pos(2)-2*hBtn;
  36.         popstr = strvcat(...
  37.                         'woman2','sinsin','detfingr', ...
  38.                         'wbarb','detail','geometry', ...
  39.                         'tire','mandrill','woman',...
  40.                         'wifs','facets','tartan');
  41.         txtstr = 'Image';
  42. propUIC = {'Units',propAuto{1},'Position',pos,'BackGroundColor',txtBkColor};
  43. slideData.txtLocHandle = uicontrol('Parent',figHandle,...
  44. 'style','text','String',txtstr,propUIC{:});
  45. pos(2) = pos(2)-hBtn+hBtn/3;
  46. propUIC = {'Units',propAuto{1},'Position',pos,'BackGroundColor',propAuto{3}};
  47. propAuto{2} = pos;
  48. propUIC{4} = propAuto{2};
  49. slideData.popLocHandle = uicontrol('Parent',figHandle,...
  50. 'style','popupmenu','String',popstr,propUIC{:});
  51. set(figHandle,'Userdata',slideData);
  52. set(figHandle,'DefaultAxesXtick',[],'DefaultAxesYtick',[]);
  53. end
  54. return
  55. end
  57. if nargout<1,
  58.   wshowdrv(mfilename)
  59. else
  60.   idx = 0;
  61.   %========== Slide 1 ==========
  62.   idx = idx+1;
  63.   slide(idx).code = {
  64. 'figHandle = gcf;',
  65. 'ax  = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',
  66. 'popTMP = findobj(figHandle,''style'',''popupmenu'');',
  67. 'set(popTMP,''enable'',''on'');',
  68. '' };
  70.   slide(idx).text = {
  71. '',
  72. ' Press the "Start" button to see a demonstration of',
  73. ' 2-D wavelet tools in the Wavelet Toolbox.',
  74. '',
  75. ' This demo uses Wavelet Toolbox functions.',
  76. ''};
  78.   %========== Slide 2 ==========
  79.   idx = idx+1;
  80.   slide(idx).code = {
  81. 'pop = findobj(figHandle,''style'',''popupmenu'');',
  82. 'lstIMG = get(pop,''String''); idxIMG = get(pop,''Value'');',
  83. 'set(pop,''enable'',''off'');',
  84. 'imgName = deblankl(lstIMG(idxIMG,:));',
  85. 'feval(''load'',imgName)',
  86. 'newTXT = [''         load '',imgName];';
  87. 'wshowdrv(''#modify_Comment'',figHandle,2,newTXT);',
  88. '[r,c] = size(X);'
  89. 'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',
  91. 'if strcmp(imgName,''detail'')',
  92. ' c = r; X = X(1:2:r,1:2:c);',
  93. 'elseif strcmp(imgName,''mandrill'')',
  94. ' X = X(1:4:r,1:5:c);',
  95. 'elseif strcmp(imgName,''tire'')',
  96. ' X = X(1:204,:);',
  97. 'elseif strcmp(imgName,''detfingr'')',
  98. ' r = 296; c = 296; X = X(1:2:r,1:2:c);',
  99. 'elseif strcmp(imgName,''durer'')',
  100. ' r = 640; X = X(1:4:r,1:4:c);',
  101. 'elseif strcmp(imgName,''wifs'')',
  102. ' ',
  103. 'end',
  104. '[r,c] = size(X);',
  105. 'if rem(r,4)>0 | rem(c,4)>0',
  106. ' X = X(1:r-rem(r,4),1:c-rem(c,4));',
  107. 'end',
  108. 'sX = size(X);',
  109. 'row = sX(1); col = sX(2);',
  110. 'nbcol = size(map,1);',
  111. 'cod_X = wcodemat(X,nbcol);',     
  112. 'row_1 = 1; col_1 = 2;',
  113. 'colormap(pink(nbcol));',
  114. 'h(1) = subplot(row_1,col_1,1);',
  115. 'image(cod_X); axis(''square'');',
  116. 'title(''Original image X.'');',
  117. '' };
  119.   slide(idx).text = {
  120. ' % Load original image.',
  121. '         load  imgName',
  122. '         % X is the loaded image.',
  123. '         % and map the loaded colormap.',
  124. ' ',
  125. ' % Image coding.',
  126. '         nbcol = size(map,1);',
  127. '         cod_X = wcodemat(X,nbcol);',
  128. ''};
  130.   slide(idx).info = 'wcodemat';
  132.   %========== Slide 3 ==========
  133.   idx = idx+1;
  134.   slide(idx).code = {
  135.   '[ca1,ch1,cv1,cd1] = dwt2(X,''db1'');',
  136.   'cod_ca1 = wcodemat(ca1,nbcol);',
  137.   'cod_ch1 = wcodemat(ch1,nbcol);',
  138.   'cod_cv1 = wcodemat(cv1,nbcol);',
  139.   'cod_cd1 = wcodemat(cd1,nbcol);',
  140.   'ax = subplot(row_1,col_1,2);',
  141.   'lineCOL = ''r'';',
  142.   'image([cod_ca1,cod_ch1;cod_cv1,cod_cd1]); axis(''square'');',
  143.   'line(''Xdata'',[1 col],''Ydata'',[row/2 row/2],''Color'',lineCOL);',
  144.   'line(''Xdata'',[col/2 col/2],''Ydata'',[1 row],''Color'',lineCOL);',
  145.   'title(''Decomposition at level 1'');',
  146.   '' };
  148.   slide(idx).text = {
  149.   ' % Perform one step decomposition.',
  150.   '        [ca1,ch1,cv1,cd1] = dwt2(X,''db1'');',
  151.   ' '
  152.   ' % Image coding.'
  153.   '        cod_ca1 = wcodemat(ca1,nbcol);',
  154.   '        cod_ch1 = wcodemat(ch1,nbcol);',
  155.   '        cod_cv1 = wcodemat(cv1,nbcol);',
  156.   '        cod_cd1 = wcodemat(cd1,nbcol);',
  157.   '        image([cod_ca1,cod_ch1;cod_cv1,cod_cd1]);',
  158.   ''};
  160.   slide(idx).info = 'dwt2';
  162.   %========== Slide 4 ==========
  163.   idx = idx+1;
  164.   slide(idx).code = {
  165.   'nul = zeros(size(cod_ca1));',
  166.   'subplot(row_1,col_1,2);',
  167.   'image([cod_ca1,nul;nul,nul]); axis(''square'');',
  168.   'line(''Xdata'',[1 col],''Ydata'',[row/2 row/2],''Color'',lineCOL);',
  169.   'line(''Xdata'',[col/2 col/2],''Ydata'',[1 row],''Color'',lineCOL);',
  170.   'title(''Coef. of Approximation at level 1'');',
  171.   ''};
  173.   slide(idx).text = slide(idx-1).text;
  175.   slide(idx).info = 'dwt2';
  177.   %========== Slide 5 ==========
  178.   idx = idx+1;
  179.   slide(idx).code = {
  180.   'subplot(row_1,col_1,2);',
  181.   'image([nul,cod_ch1;nul,nul]); axis(''square'');',
  182.   'line(''Xdata'',[1 col],''Ydata'',[row/2 row/2],''Color'',lineCOL);',
  183.   'line(''Xdata'',[col/2 col/2],''Ydata'',[1 row],''Color'',lineCOL);',
  184.   'title(''Coef. of Horizontal Detail at level 1'');',
  185.   '' };
  187.   slide(idx).text = slide(idx-1).text;
  189.   slide(idx).info = 'dwt2';
  191.   %========== Slide 6 ==========
  192.   idx = idx+1;
  193.   slide(idx).code = {
  194.   'subplot(row_1,col_1,2);',
  195.   'image([nul,nul;cod_cv1,nul]); axis(''square'');',
  196.   'line(''Xdata'',[1 col],''Ydata'',[row/2 row/2],''Color'',lineCOL);',
  197.   'line(''Xdata'',[col/2 col/2],''Ydata'',[1 row],''Color'',lineCOL);',
  198.   'title(''Coef. of Vertical Detail at level 1'');',
  199.   '' };
  201.   slide(idx).text = slide(idx-1).text;
  203.   %========== Slide 7 ==========
  204.   idx = idx+1;
  205.   slide(idx).code = {
  206.   'subplot(row_1,col_1,2);',
  207.   'image([nul,nul;nul,cod_cd1]); axis(''square'');',
  208.   'line(''Xdata'',[1 col],''Ydata'',[row/2 row/2],''Color'',lineCOL);',
  209.   'line(''Xdata'',[col/2 col/2],''Ydata'',[1 row],''Color'',lineCOL);',
  210.   'title(''Coef. of Diagonal Detail at level 1'');',
  211.   '' };
  213.   slide(idx).text = slide(idx-1).text;
  215.   slide(idx).info = 'dwt2';
  217.   %========== Slide 8 ==========
  218.   idx = idx+1;
  219.   slide(idx).code = {
  220.   'subplot(row_1,col_1,2);',
  221.   'image([cod_ca1,cod_ch1;cod_cv1,cod_cd1]);axis(''square'');',
  222.   'line(''Xdata'',[1 col],''Ydata'',[row/2 row/2],''Color'',lineCOL);',
  223.   'line(''Xdata'',[col/2 col/2],''Ydata'',[1 row],''Color'',lineCOL);',
  224.   'title(''Decomposition at level 1'');',
  225.   'clear nul',
  226.   '' };
  228.   slide(idx).text = slide(idx-1).text;
  230.   slide(idx).info = 'dwt2';
  232.   %========== Slide 9 ==========
  233.   idx = idx+1;
  234.   slide(idx).code = {
  235.  ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  236.   'if idxPREV>idxSlide',
  237.   '   subplot(row_1,col_1,1);',
  238.   '   image(cod_X); axis(''square'');',
  239.   '   title(''Original image X.'');',
  240.   'end',
  241.   '[ca2,ch2,cv2,cd2] = dwt2(ca1,''db1'');',
  242.   'cod_ca2 = wcodemat(ca2,nbcol);        clear ca2',
  243.   'cod_ch2 = wcodemat(ch2,nbcol);        clear ch2',
  244.   'cod_cv2 = wcodemat(cv2,nbcol);        clear cv2',
  245.   'cod_cd2 = wcodemat(cd2,nbcol);        clear cd2',
  246.   'cod_ca1 = [cod_ca2,cod_ch2;cod_cv2,cod_cd2];',
  247.   'clear cod_ca2 cod_ch2 cod_cv2 cod_cd2',
  248.   'subplot(row_1,col_1,2)',
  249.   'image([cod_ca1,cod_ch1;cod_cv1,cod_cd1]); axis(''square'');',
  250.   'line(''Xdata'',[1 col/2],''Ydata'',[row/4 row/4],''Color'',lineCOL);',
  251.   'line(''Xdata'',[col/4 col/4],''Ydata'',[1 row/2],''Color'',lineCOL);',
  252.   'line(''Xdata'',[1 col],''Ydata'',[row/2 row/2],''Color'',lineCOL);',
  253.   'line(''Xdata'',[col/2 col/2],''Ydata'',[1 row],''Color'',lineCOL);',
  254.   'title(''Decomposition at level 2'');',
  255.   '' };
  257.   slide(idx).text = {
  258.   ' % Perform second step decomposition :'
  259.   ' % decompose approx. cfs of level 1.'
  260.   '        [ca2,ch2,cv2,cd2] = dwt2(ca1,''db1'');'
  261.   ' ' 
  262.   ' % Code ca2, ch2,cv2 and cd2.'
  263.   '        cod_ca2 = wcodemat(ca2,nbcol);'
  264.   '        cod_ch2 = wcodemat(ch2,nbcol);'
  265.   '        cod_cv2 = wcodemat(cv2,nbcol);'
  266.   '        cod_cd2 = wcodemat(cd2,nbcol);',
  267.   '' };
  269.   slide(idx).info = 'dwt2';
  271.   %========== Slide 10 ==========
  272.   idx = idx+1;
  273.   slide(idx).code = {
  274.   'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  275.   'image([cod_ca1,cod_ch1;cod_cv1,cod_cd1]); axis(''square'');',
  276.   'line(''Xdata'',[1 col/2],''Ydata'',[row/4 row/4],''Color'',lineCOL);',
  277.   'line(''Xdata'',[col/4 col/4],''Ydata'',[1 row/2],''Color'',lineCOL);',
  278.   'line(''Xdata'',[1 col],''Ydata'',[row/2 row/2],''Color'',lineCOL);',
  279.   'line(''Xdata'',[col/2 col/2],''Ydata'',[1 row],''Color'',lineCOL);',
  280.   'title(''Decomposition at level 2'');',
  281.   '' };
  283.   slide(idx).text = {
  284.   ' % Zoom the decomposition at level 2.', 
  285.   '' };
  287.   slide(idx).info = 'dwt2';
  289.   %========== Slide 11 ==========
  290.   idx = idx+1;
  291.   slide(idx).code = {
  292.   'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',   
  293.   'row_2 = 1; col_2 = 3;',
  294.   'h(1) = subplot(row_2,col_2,1); image(cod_X);',
  295.   'axis(''square'');',
  296.   'title(''Original image X.'');',
  297.   'drawnow;',
  298.   'cod_ca1  = wcodemat(ca1,nbcol);',
  299.   'h(3) = subplot(row_2,col_2,3);',
  300.   'image([cod_ca1,cod_ch1;cod_cv1,cod_cd1]);',
  301.   'axis(''square'');',
  302.   'clear cod_ca1 cod_ch1 cod_cv1 cod_cd1',
  303.   'line(''Xdata'',[1 col],''Ydata'',[row/2 row/2],''Color'',lineCOL);',
  304.   'line(''Xdata'',[col/2 col/2],''Ydata'',[1 row],''Color'',lineCOL);',
  305.   'title(''Decomp. at level 1'');',
  306.   'a0 = idwt2(ca1,ch1,cv1,cd1,''db1'',sX);',
  307.   'h(2) = subplot(row_2,col_2,2); image(wcodemat(a0,nbcol));',
  308.   'axis(''square'');',
  309.   'title(''Recons. image : a0.'');',
  310.   '' };
  312.   slide(idx).text = {
  313.   ' % Invert directly decomposition of X' ,
  314.   ' % using coefficients at level 1.', 
  315.   ' ',
  316.   '         a0 = idwt2(ca1,ch1,cv1,cd1,''db1'',size(X));', 
  317.   ' ',
  318.   '' };
  320.   slide(idx).info = 'idwt2';
  322.   %========== Slide 12 ==========
  323.   idx = idx+1;
  324.   slide(idx).code = {
  325.   'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = []; drawnow;',
  326.   'h(1) = subplot(row_2,col_2,1); image(cod_X);',
  327.   'axis(''square'');',
  328.   'title(''Original image X.'');',
  329.   '[c,s] = wavedec2(X,2,''db1'');',   
  330.   '' };
  332.   slide(idx).text = {
  333.   ' % Perform decomposition at level 2',
  334.   ' % of X using db1.',
  335.   ' ',
  336.   '         [c,s] = wavedec2(X,2,''db1'');',
  337.   '' };
  339.   slide(idx).info = 'wavedec2';
  341.   %========== Slide 13 ==========
  342.   idx = idx+1;
  343.   slide(idx).code = {
  344.  ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  345.   'if idxPREV>idxSlide',
  346.   '   delete(subplot(row_2,col_2,3));',
  347.   'end',
  348.   'ca2 = appcoef2(c,s,''db1'',2);',
  349.   'ch2 = detcoef2(''h'',c,s,2);',
  350.   'cv2 = detcoef2(''v'',c,s,2);',
  351.   'cd2 = detcoef2(''d'',c,s,2);',
  352.   'cod_ca2 = wcodemat(ca2,nbcol); clear ca2',
  353.   'cod_ch2 = wcodemat(ch2,nbcol); clear ch2',
  354.   'cod_cv2 = wcodemat(cv2,nbcol); clear cv2',
  355.   'cod_cd2 = wcodemat(cd2,nbcol); clear cd2',
  356.   'cod_ca1 = [cod_ca2,cod_ch2;cod_cv2,cod_cd2];',
  357.   'clear cod_ca2 cod_ch2 cod_cv2 cod_cd2',
  358.   'nul = zeros(size(cod_ca1));',
  359.   'h(2) = subplot(row_2,col_2,2); image([cod_ca1,nul;nul,nul]);',
  360.   'axis(''square'');',
  361.   'clear nul',
  362.   'line(''Xdata'',[1 col/2],''Ydata'',[row/4 row/4],''Color'',lineCOL);',
  363.   'line(''Xdata'',[col/4 col/4],''Ydata'',[1 row/2],''Color'',lineCOL);',
  364.   'line(''Xdata'',[1 col],''Ydata'',[row/2 row/2],''Color'',lineCOL);',
  365.   'line(''Xdata'',[col/2 col/2],''Ydata'',[1 row],''Color'',lineCOL);',
  366.   'title(''Coef. at level 2.'');',
  367.   '' };
  369.   slide(idx).text = {
  370.   ' % Extract coefficients at level 2 from',
  371.   ' % the wavelet decomposition structure [c,s].',
  372.   ' ',
  373.   '         ca2 = appcoef2(c,s,''db1'',2);',
  374.   '         ch2 = detcoef2(''h'',c,s,2);',
  375.   '         cv2 = detcoef2(''v'',c,s,2);',
  376.   '         cd2 = detcoef2(''d'',c,s,2);',
  377.   '' };
  379.   slide(idx).info = 'appcoef2';
  381.   %========== Slide 14 ==========
  382.   idx = idx+1;
  383.   slide(idx).code = {
  384.  ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  385.   'if idxPREV>idxSlide',
  386.   '   subplot(row_2,col_2,1);',
  387.   '   image(cod_X); axis(''square'');',
  388.   '   title(''Original image X.'');',
  389.   'end',
  390.   'ca1 = appcoef2(c,s,''db1'',1);',
  391.   'ch1 = detcoef2(''h'',c,s,1);',
  392.   'cv1 = detcoef2(''v'',c,s,1);',
  393.   'cd1 = detcoef2(''d'',c,s,1);',
  394.   'cod_ca1 = wcodemat(ca1,nbcol);       clear ca1',
  395.   'cod_ch1 = wcodemat(ch1,nbcol);       clear ch1',
  396.   'cod_cv1 = wcodemat(cv1,nbcol);       clear cv1',
  397.   'cod_cd1 = wcodemat(cd1,nbcol);       clear cd1',
  398.   'h(3) = subplot(row_2,col_2,3); image([cod_ca1,cod_ch1;cod_cv1,cod_cd1]);',
  399.   'axis(''square'');',
  400.   'clear cod_ca1 cod_ch1 cod_cv1 cod_cd1',
  401.   'line(''Xdata'',[1 col],''Ydata'',[row/2 row/2],''Color'',lineCOL);',
  402.   'line(''Xdata'',[col/2 col/2],''Ydata'',[1 row],''Color'',lineCOL);',
  403.   'title(''Coef. at level 1.'');',
  404.   '' };
  406.   slide(idx).text = {
  407.   ' % Extract coefficients at level 1,',
  408.   ' % from wavelet decomposition structure [c,s].',
  409.   ' '
  410.   '        ca1 = appcoef2(c,s,''db1'',1);',
  411.   '        ch1 = detcoef2(''h'',c,s,1);',
  412.   '        cv1 = detcoef2(''v'',c,s,1);',
  413.   '        cd1 = detcoef2(''d'',c,s,1);',
  414.   '' };
  416.   slide(idx).info = 'detcoef2';
  418.   %========== Slide 15 ==========
  419.   idx = idx+1;
  420.   slide(idx).code = {
  421.   'ax = findall(figHandle,''type'',''axes''); delete(ax); drawnow;  h = [];',
  422.   'row_3 = 3; col_3 = 4;';
  423.   'ax = subplot(row_3,col_3,1); image(cod_X); axis(''square'');',
  424.   'title(''X'');',
  425.   '' };
  427.   slide(idx).text = {
  428.   ' % Plot the original image.',
  429.   '' };
  431.   slide(idx).info = 'detcoef2';
  433.   slide(idx).idxPrev = 12;
  435.   %========== Slide 16 ==========
  436.   idx = idx+1;
  437.   slide(idx).code = {
  438.  ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  439.   'if idxPREV>idxSlide',
  440.   '   for j = 1:3 , tmp(j) = subplot(row_3,col_3,j+9); end;',
  441.   '   delete(tmp); clear tmp',
  442.   'end',
  443.   'a2 = wrcoef2(''a'',c,s,''db1'',2);',
  444.   'ax = subplot(row_3,col_3,9); image(wcodemat(a2,nbcol)); axis(''square'');',
  445.   'title(''a2'')',
  446.   'clear a2',
  447.   '' };
  449.   slide(idx).text = {
  450.   ' % Reconstruct approximation at level 2 from',
  451.   ' % the wavelet decomposition structure [c,s].',
  452.   ' ',
  453.   '        a2 = wrcoef2(''a'',c,s,''db1'',2);',
  454.   '' };
  456.   slide(idx).info = 'wrcoef2';
  458.   %========== Slide 17 ==========
  459.   idx = idx+1;
  460.   slide(idx).code = {
  461.  ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  462.   'if idxPREV<idxSlide',
  463.   '    h2 = wrcoef2(''h'',c,s,''db1'',2);',
  464.   '    v2 = wrcoef2(''v'',c,s,''db1'',2);',
  465.   '    d2 = wrcoef2(''d'',c,s,''db1'',2);',
  466.   '    ax = subplot(row_3,col_3,10); image(wcodemat(h2,nbcol)); axis(''square'');',
  467.   '    title(''h2'');',
  468.   '    clear h2',
  469.   '    ax = subplot(row_3,col_3,11); image(wcodemat(v2,nbcol)); axis(''square'');',
  470.   '    title(''v2'');',
  471.   '    clear v2',
  472.   '    ax = subplot(row_3,col_3,12); image(wcodemat(d2,nbcol)); axis(''square'');',
  473.   '    title(''d2'');',
  474.   '    clear d2',
  475.   'end',
  476.   '' };
  478.   slide(idx).text = {
  479.   ' % Reconstruct detail at level 2 from',
  480.   ' % the wavelet decomposition structure [c,s].',
  481.   '         h2 = wrcoef2(''h'',c,s,''db1'',2);',
  482.   '         v2 = wrcoef2(''v'',c,s,''db1'',2);',
  483.   '         d2 = wrcoef2(''d'',c,s,''db1'',2);',
  484.   '' };
  486.   slide(idx).info = 'wrcoef2';
  488.   %========== Slide 18 ==========
  489.   idx = idx+1;
  490.   slide(idx).code = {
  491.   'sc = size(c);',
  492.   'val_s = s;',
  493.   '[c,s] = upwlev2(c,s,''db1'');',
  494.   'sc = size(c);',
  495.   'val_s = s;',
  496.   '' };
  498.   slide(idx).text = {
  499.   ' % One step reconstruction of wavelet',
  500.   ' % decomposition structure [c,s].',
  501.   '        [c,s] = upwlev2(c,s,''db1'');',
  502.   '' };
  504.   slide(idx).info = 'upwlev2';
  506.   %========== Slide 19 ==========
  507.   idx = idx+1;
  508.   slide(idx).code = {
  509.  ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  510.   'if idxPREV>idxSlide',
  511.   '   for j = 1:4 , tmp(j) = subplot(row_3,col_3,j+4); end; ',
  512.   '   delete(tmp); clear tmp',
  513.   'end',   
  514.   '' };
  516.   slide(idx).text = {
  517.   ' % Reconstruct approximation and details',
  518.   ' % at level 1, from coefficients.',
  519.   ' %'
  520.   ' % step 1 : extract coefficients',
  521.   ' % decomposition structure [c,s].',
  522.   ' %',
  523.   ' % step 2 : reconstruct.',
  524.   '' };
  526.   slide(idx).info = '';
  528.   %========== Slide 20 ==========
  529.   idx = idx+1;
  530.   slide(idx).code = {
  531.  ['idxPREV = wshowdrv(''#get_idxSlide'',figHandle); idxSlide = ' int2str(idx) ';'],
  532.   'if idxPREV>idxSlide',
  533.   '   delete(subplot(row_3,col_3,3));',
  534.   'else',
  535.   '   siz = s(size(s,1),:);',
  536.   '   ca1 = appcoef2(c,s,''db1'',1);',
  537.   '   a1  = upcoef2(''a'',ca1,''db1'',1,siz);',
  538.   '   clear ca1',
  539.   '   ch1 = detcoef2(''h'',c,s,1);',
  540.   '   hd1 = upcoef2(''h'',ch1,''db1'',1,siz);',
  541.   '   clear ch1',
  542.   '   cv1 = detcoef2(''v'',c,s,1);',
  543.   '   vd1 = upcoef2(''v'',cv1,''db1'',1,siz);',
  544.   '   clear cv1',
  545.   '   cd1 = detcoef2(''d'',c,s,1);',
  546.   '   dd1 = upcoef2(''d'',cd1,''db1'',1,siz);',
  547.   '   clear cd1',
  548.   '   ax = subplot(row_3,col_3,5); image(wcodemat(a1,nbcol)); axis(''square'');',
  549.   '   title(''a1'');',
  550.   '   clear a1',
  551.   '   ax = subplot(row_3,col_3,6); image(wcodemat(hd1,nbcol)); axis(''square'');',
  552.   '   title(''h1'');',
  553.   '   clear hd1',
  554.   '   ax = subplot(row_3,col_3,7); image(wcodemat(vd1,nbcol)); axis(''square'');',
  555.   '   title(''v1'');',
  556.   '   clear vd1',
  557.   '   ax = subplot(row_3,col_3,8); image(wcodemat(dd1,nbcol)); axis(''square'');',
  558.   '   title(''d1'');',
  559.   '   clear dd1',
  560.   'end'
  561.   '' };
  563.   slide(idx).text = {
  564.   '    siz = s(size(s,1),:);',
  565.   '    ca1 = appcoef2(c,s,''db1'',1);',
  566.   '    ch1 = detcoef2(''h'',c,s,1);',
  567.   '    cv1 = detcoef2(''v'',c,s,1);',
  568.   '    cd1 = detcoef2(''d'',c,s,1);',
  569.   '    a1  = upcoef2(''a'',ca1,''db1'',1,siz);',
  570.   '    hd1 = upcoef2(''h'',ch1,''db1'',1,siz);',
  571.   '    vd1 = upcoef2(''v'',cv1,''db1'',1,siz);',
  572.   '    dd1 = upcoef2(''d'',cd1,''db1'',1,siz);',
  573.   '' };
  575.   slide(idx).info = 'upcoef2';
  577.   %========== Slide 21 ==========
  578.   idx = idx+1;
  579.   slide(idx).code = {
  580.   'a0 = waverec2(c,s,''db1'');',
  581.   'subplot(row_3,col_3,3); image(wcodemat(a0,nbcol)); axis(''square'');',
  582.   'title(''Reconstructed image a0.'');',
  583.   '' };
  585.   slide(idx).text = {
  586.   ' % Reconstruct X from the wavelet',
  587.   ' % decomposition structure [c,s].',
  588.   ' ',
  589.   '         a0 = waverec2(c,s,''db1'');',
  590.   '' };
  592.   slide(idx).info = 'waverec2';
  594.   varargout{1} = slide;
  595. end