wavedec2_spiht.m

输入：为需要压缩图象的名称,该主程序仅能构处理256灰度图,读者可以自行改编为RGB处理；ratio为压缩比率；level为小波分解的级数 输出:是解压完毕的图象数据矩阵

function pic=wavdec2_SPIHT(pic,wavename,scale)
%function pic=wav_de3(img,scale)
%用法: pic=wav_de3(pic,scale)
%功能：将输入的参数pic进行scale级的小波分解并输出结果
%pic=imread(img);
%M,N为输入参数pic的尺寸（行数和列数）
pic=double(pic);
[M,N]=size(pic);
%一共作scale级小波分解，因此变量i从1一直增加到scale
for i=1:scale
   %X,Y为当前级待分解块的尺寸
   X=M/2^i;
   Y=N/2^i;
   [a,h,v,d]=dwt2(pic(1:2*X,1:2*Y),wavename);
   pic(1:X,1:Y)=a(10:X+9,10:Y+9);
   pic(1:X,Y+1:2*Y)=[h(6:X+4,10:Y+9);zeros(1,Y)];
   pic(X+1:2*X,1:Y)=[v(10:X+9,6:Y+4) zeros(X,1)];
   pic(X+1:2*X,Y+1:2*Y)=d(7:X+6,7:Y+6);
end