Logo AND Algorithmique Numérique Distribuée

Private GIT Repository
final avant rapport
[these_gilles.git] / 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