wav_de3.m

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

function pic=wav_de3(pic,scale)
%用法: pic=wav_de3(pic,scale)
%功能：将输入的参数pic进行scale级的小波分解并输出结果

%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),'coif5');
   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