THESE/codes/wave/allcode/dwt2D.m

   1 function w = dwt2D(x, J, af)
   2
   3 % discrete 2-D wavelet transform
   4 %
   5 % USAGE:
   6 %   w = dwt2D(x, stages, af)
   7 % INPUT:
   8 %   x - N by M matrix
   9 %       1) M, N are both even
  10 %       2) min(M,N) >= 2^(J-1)*length(af)
  11 %   J - number of stages
  12 %   af - analysis filters
  13 % OUPUT:
  14 %   w - cell array of wavelet coefficients
  15 % EXAMPLE:
  16 %   [af, sf] = farras;
  17 %   x = rand(128,64);
  18 %   J = 3;
  19 %   w = dwt2D(x,J,af);
  20 %   y = idwt2D(w,J,sf);
  21 %   err = x - y;
  22 %   max(max(abs(err)))
  23 %
  24 % WAVELET SOFTWARE AT POLYTECHNIC UNIVERSITY, BROOKLYN, NY
  25 % http://taco.poly.edu/WaveletSoftware/
  26
  27 for k = 1:J
  28     [x w{k}] = afb2D(x, af, af);
  29 end
  30 w{J+1} = x;
  31