reconstructdwttile.m

SPIHT coding implementation. using wavelet transforms as a key tool. compression and decoding is don

function [rec] = reconstructDwtTile(tile, Lo_R, Hi_R, level)

% reconstructed image assembly
%
% parameters:   tile - decomposition image
%               level - transform depth
%               Lo_R, Hi_R - reconstruction filter coeficients
% output:       rec - reconstructed image



im_size = size(tile, 1);
chng = 2^(log2(im_size)-level);

cA = tile(1:chng,1:chng);

for i=1:level
    low = 0;
    med = chng;
    high = 2*chng;
    
    cH = tile(low+1:med, med+1:high);
    cV = tile(med+1:high, low+1:med );
    cD = tile(med+1:high, med+1:high);
    
    cA = idwt2(cA, cH, cV, cD, Lo_R, Hi_R);
    
    chng = chng*2;
end

rec = cA;

rec = min(255,rec);
rec = max(0,rec);