function x2 = idwtdaub2(xdaub2, J, Hord)
% Inversa Daubechies wavelet transform i två dimensioner
%       x2 = idwtdaub2(xdaub2, J, Hord)
%
% xdaub2 Daubechies wavelet transform av x2 med filter av ordningen 
%       Hord och J uppdelningssteg.
% J     Antalet uppdelningssteg.
% Hord  Filterordning av Daubechies filter.
% x2    Den utgående från xdaub2 rekonstruerad matris
%

sz = size(xdaub2);
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);

x1 = zeros(nrows,ncols);
x2 = zeros(nrows,ncols);

% Inverse transform first columnwise, then rowwise
for j=1:ncols
    x1(:,j) = idwtdaub( xdaub2(:,j), J, H );
end

for i=1:nrows
    x2(i,:) = idwtdaub( x1(i,:), J, H );
end