THESE/codes/wave/allcode/sfb2D_A.m

   1 function y = sfb2D_A(lo, hi, sf, d)
   2
   3 % 2D Synthesis Filter Bank
   4 % (along single dimension only)
   5 %
   6 % y = sfb2D_A(lo, hi, sf, d);
   7 % sf - synthesis filters
   8 % d  - dimension of filtering
   9 % see afb2D_A
  10
  11
  12 lpf = sf(:, 1);     % lowpass filter
  13 hpf = sf(:, 2);     % highpass filter
  14
  15 if d == 2
  16    lo = lo';
  17    hi = hi';
  18 end
  19
  20 N = 2*size(lo,1);
  21 L = length(sf);
  22 y = upfirdn(lo, lpf, 2, 1) + upfirdn(hi, hpf, 2, 1);
  23 y(1:L-2, :) = y(1:L-2, :) + y(N+[1:L-2], :);
  24 y = y(1:N, :);
  25 y = cshift2D(y, 1-L/2);
  26
  27 if d == 2
  28    y = y';
  29 end
  30