trans=sum(A01);
a1=recon(trans,basic_info,a_fromlevel,a_tolevel,dimx,dimy);
surf(t1,t2,a1)
surf(t1,t2,A/2)
surf(t1,t2,A/2-a1)

trans=sum(A02);
b1=recon(trans,basic_info,a_fromlevel,a_tolevel,dimx,dimy);
surf(t1,t2,b1)
surf(t1,t2,X)
surf(t1,t2,X-b1)

B=(5*t1-6*t2-5);
[C2,S2] = wavedec2(B,b_tolevel,'bior1.3');
B1=waverec2(C2,S2,'bior1.3');