THESE/codes/wave/allcode/denR2D.m

   1 function y = denR2D(x,T);
   2
   3 % % Example
   4 % s1 = double(imread('st.tif'));
   5 % s = s1(:,:,3);
   6 % x = s + 20*randn(size(s));
   7 % T = 10;
   8 % y = denR2D(x,T);
   9 % imagesc(y)
  10 % colormap(gray)
  11 % axis image
  12 % sqrt(mean(mean((y-s).^2)))
  13
  14 [Faf, Fsf] = FSfarras;
  15 [af, sf] = dualfilt1;
  16 J = 4;
  17 w = dualtree2D(x,J,Faf,af);
  18 % loop thru scales:
  19 for j = 1:J
  20     % loop thru subbands
  21     for s1 = 1:2
  22         for s2 = 1:3
  23             w{j}{s1}{s2} = soft(w{j}{s1}{s2},T);
  24         end
  25     end
  26 end
  27 y = idualtree2D(w,J,Fsf,sf);
  28