function xdaub2 = dwtdaub2(x2, J, Hord)
% Daubechies wavelet transform i två dimensioner
%       xdaub2 = dwtdaub2(x2, J, Hord)
%
% x2    En matris.
% J     Antalet uppdelningssteg.
% Hord  Filterordning av Daubechies filter.
% xdaub2 Daubechies wavelet transform av x2 med filter av ordningen 
%       Hord och J uppdelningssteg.
%

sz = size(x2);
nrows = sz(1);
ncols = sz(2);

% kolla att nrows och ncols är dividerbara med 2^J
if ( rem(nrows, 2^J) ~= 0 ) || ( rem(ncols, 2^J) ~= 0 )
    error('The dimensions of x2 must be multiples of 2^J');
end

% beräkna Daubechies filterkoefficienter
H = hdaub(Hord);

xdaub1 = zeros(nrows,ncols);
xdaub2 = zeros(nrows,ncols);

% Transform first rowswise, then columnwise
for i=1:nrows
    xdaub1(i,:) = dwtdaub( x2(i,:), J, H );
end

for j=1:ncols
    xdaub2(:,j) = dwtdaub( xdaub1(:,j), J, H );
end