THESE/codes/wave/allcode/afb2D.m

   1 function [lo, hi] = afb2D(x, af1, af2)
   2
   3 % 2D Analysis Filter Bank
   4 %
   5 % USAGE:
   6 %   [lo, hi] = afb2D(x, af1, af2);
   7 % INPUT:
   8 %   x - N by M matrix
   9 %       1) M, N are both even
  10 %       2) M >= 2*length(af1)
  11 %       3) N >= 2*length(af2)
  12 %   af1 - analysis filters for columns
  13 %   af2 - analysis filters for rows
  14 % OUTPUT:
  15 %    lo - lowpass subband
  16 %    hi{1} - 'lohi' subband
  17 %    hi{2} - 'hilo' subband
  18 %    hi{3} - 'hihi' subband
  19 % EXAMPLE:
  20 %   x = rand(32,64);
  21 %   [af, sf] = farras;
  22 %   [lo, hi] = afb2D(x, af, af);
  23 %   y = sfb2D(lo, hi, sf, sf);
  24 %   err = x - y;
  25 %   max(max(abs(err)))
  26 %
  27 % WAVELET SOFTWARE AT POLYTECHNIC UNIVERSITY, BROOKLYN, NY
  28 % http://taco.poly.edu/WaveletSoftware/
  29
  30 if nargin < 3
  31    af2 = af1;
  32 end
  33
  34 % filter along columns
  35 [L, H] = afb2D_A(x, af1, 1);
  36
  37 % filter along rows
  38 [lo,    hi{1}] = afb2D_A(L, af2, 2);
  39 [hi{2}, hi{3}] = afb2D_A(H, af2, 2);
  40