THESE/codes/wave/allcode/den3.m

   1
   2
   3 s1 = double(imread('st.tif'));
   4 s = s1(:,:,3);
   5 x = s + 20*randn(size(s));
   6 t = 0:5:50;         % threshold range(0~50),increment by 5
   7 e = den2(s,x,t);    % using standard method
   8 re = rden2(s,x,t);  % using reduced 2-D dual-tree method
   9 ce = cden2(s,x,t);  % using complex 2-D dual-tree method
  10
  11 figure(1)
  12 plot(t,e)
  13 hold on
  14 plot(t,re,'r')
  15 hold on
  16 plot(t,ce,'k')
  17 title('RMS error V.S. Threshold Pt.')
  18 xlabel('Threshold pt.');
  19 ylabel('RMS error');
  20 legend('Standard 2D','Reduced 2D dual','Complx 2D dual',0);
  21
  22
  23 [emin,k] = min(e);
  24 T = t(k);
  25 y = denS2D(x,T);
  26 figure(2)
  27 clf
  28 image(y/4)
  29 colormap(gray)
  30 axis image
  31 axis off
  32 print -deps den3_A
  33
  34 ax = [280*[1 1] 250*[1 1]]+[0 1 0 1]*120;
  35
  36 figure(3)
  37 clf
  38 image(y/4)
  39 colormap(gray)
  40 axis(ax)
  41 axis square
  42 axis off
  43 print -deps den3_B
  44
  45 [emin,k] = min(ce);
  46 T = t(k);
  47 y = denC2D(x,T);
  48 figure(4)
  49 clf
  50 image(y/4)
  51 colormap(gray)
  52 axis image
  53 axis off
  54 print -deps cden3_A
  55
  56 figure(5)
  57 clf
  58 image(y/4)
  59 colormap(gray)
  60 axis(ax)
  61 axis square
  62 axis off
  63 print -deps cden3_B
  64