THESE/codes/wave/allcode/denS2D.m

   1 function y = denS2D(x,T)
   2
   3 % % Example
   4 % s1 = double(imread('st.tif'));
   5 % s = s1(:,:,3);
   6 % x = s + 20*randn(size(s));
   7 % T = 35;
   8 % y = denS2D(x,T);
   9 % imagesc(y)
  10 % colormap(gray)
  11 % axis image
  12
  13 [af, sf] = farras;
  14 J = 4;
  15 w = dwt2D(x,J,af);
  16 % loop thru scales:
  17 for j = 1:J
  18     % loop thru subbands
  19     for s = 1:3
  20         w{j}{s} = soft(w{j}{s},T);
  21     end
  22 end
  23 y = idwt2D(w,J,sf);
  24