function wt = dwtreshape(c,s)
% function wt = dwtreshape(c,s)
% converts output of wavedec2 to "nested" decomposition image
% assumes dwtmode('per')

level=size(s,1)-2;
M1 = prod(s(1,:));
wt = reshape(c(1:M1),s(1,:));
 
for i = 2:level+1,
    d_size = s(i,:);
    M = prod(d_size);
    h_coef = c(M1+1:(M1+M));
    v_coef = c((M1+M+1):(M1+2*M));
    d_coef = c((M1+2*M+1):(M1+3*M));
    M1 = M1 + 3*M;
    wt = [wt reshape(h_coef,d_size); reshape(v_coef,d_size) reshape(d_coef,d_size)];
end