% Forward dual-tree DWT
% Either FSfarras or AntonB function can be used to compute the stage 1 filters
-%[Faf, Fsf] = FSfarras;
-[Faf, Fsf] = AntonB;
+[Faf, Fsf] = FSfarras;
+%[Faf, Fsf] = AntonB;
[af, sf] = dualfilt1;
W = cplxdual2D(x, J, Faf, af);
W = normcoef(W,J,nor);
% Extract the image
ind = 2^(J-1)+1:2^(J-1)+L;
-y = y(ind,ind);
\ No newline at end of file
+y = y(ind,ind);