wavedec2_speck.m

这是一个不错的图像压缩的Matlab程序,希望能对你有所帮助.

function pic=wavedec2_SPECK(pic0,wavename,scale)
%function pic=wavedec2_SPECK(pic0,wavename,scale)
%用法: pic=wavedec2_SPECK(pic0,wavename,scale)
%功能：将输入的参数pic0进行scale级的小波分解并输出结果
%pic=imread(img);
%M,N为输入参数pic的尺寸（行数和列数）
pic=double(pic0);
[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