trans.m

另ㄧ筆SPIHT的寫法~~有一點複雜就是了

function ll4=trans(x,wname)

x=x-100;
dwtmode('per','nodisp');
[a b c d]=dwt2(x,wname);
[a1 b1 c1 d1]=dwt2(a,wname);
[a2 b2 c2 d2]=dwt2(a1,wname);
[a3 b3 c3 d3]=dwt2(a2,wname);
ll=[a3 b3;c3 d3];
ll2=[ll b2; c2 d2];
ll3=[ll2 b1;c1 d1];
ll4=[ll3 b;c d];
imshow(mat2gray(ll4))