THESE/codes/wave/allcode/denC2D.m

   1 function y = denC2D(x,T);
   2
   3 % % Example
   4 % s1 = double(imread('st.tif'));
   5 % s = s1(:,:,3);
   6 % x = s + 20*randn(size(s));
   7 % T = 40;
   8 % y = denC2D(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 = cplxdual2D(x,J,Faf,af);
  18 I = sqrt(-1);
  19 % loop thru scales:
  20 for j = 1:J
  21     % loop thru subbands
  22     for s1 = 1:2
  23         for s2 = 1:3
  24             C = w{j}{1}{s1}{s2} + I*w{j}{2}{s1}{s2};
  25             C = soft(C,T);
  26             w{j}{1}{s1}{s2} = real(C);
  27             w{j}{2}{s1}{s2} = imag(C);
  28         end
  29     end
  30 end
  31 y = icplxdual2D(w,J,Fsf,sf);
  32