load charlie2.mat
J=50;K=10
A=X(J:J+512,K:K+512);
load db1
h=db1;
HH=kron(h',h);
g=wrev(h);
g(2:2:2)=-g(2:2:2)
GH=kron(g',h);

N=512

AA1=conv2(HH,A);
A1=zeros(2,2);
A1(1:1:N/2,1:1:N/2)=AA1(2:2:N,2:2:N);
imagesc(A1)

HG=kron(h',g);
DD1=conv2(HG,A);
D1=zeros(2,2);
D1(1:1:N/2,1:1:N/2)=DD1(2:2:N,2:2:N);

GH=kron(g',h);
DD2=conv2(GH,A);
D2=zeros(2,2);
D2(1:1:N/2,1:1:N/2)=DD2(2:2:N,2:2:N);

GG=kron(g',g);
DD3=conv2(GG,A);
D3=zeros(2,2);
D3(1:1:N/2,1:1:N/2)=DD3(2:2:N,2:2:N);

%%%%%%%%%%Level 2
AA2=conv2(HH,A1);
A2(1:1:N/4,1:1:N/4)=AA2(2:2:N/2,2:2:N/2);
imagesc(A2)

DD21=conv2(HG,A1);
D21(1:1:N/4,1:1:N/4)=DD21(2:2:N/2,2:2:N/2);

DD22=conv2(GH,A1);
D22(1:1:N/4,1:1:N/4)=DD22(2:2:N/2,2:2:N/2);

DD23=conv2(GG,A1);
D23(1:1:N/4,1:1:N/4)=DD23(2:2:N/2,2:2:N/2);

A3=([A2 D22;D21 D23]);
% 1 level C=([A1 D2 ;D1 D3])
%C=([A3 D2 ;D1 D3]) % 2 level
colormap(copper)
%imagesc(C)
subplot(2,2,1);imagesc(A3)
subplot(2,2,2);imagesc(D2)
subplot(2,2,3);imagesc(D1)
subplot(2,2,4);imagesc(D3)
