wave_com3.m

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

function pic=wave_com3(pic,scale)
%用法: pic=wave_com(pic,scale)
%功能：将输入的参数pic进行scale级的小波合成并输出结果
%与wave_decom功能正好相反

%M,N为输入参数pic的尺寸（行数和列数）
[M,N]=size(pic);
%一共作scale级小波分解，因此变量i从1一直增加到scale
for i=scale:-1:1
   %X,Y为当前级待合成块的四个子块尺寸
   X=M/2^i;
   Y=N/2^i;
   a=zeros(X+14,Y+14);
   h=zeros(X+14,Y+14);
   v=zeros(X+14,Y+14);
   d=zeros(X+14,Y+14);
   a(10:X+9,10:Y+9)=pic(1:X,1:Y);
   h(6:X+5,10:Y+9)=pic(1:X,Y+1:2*Y);
   v(10:X+9,6:Y+5)=pic(X+1:2*X,1:Y);
   d(7:X+6,7:Y+6)=pic(X+1:2*X,Y+1:2*Y);
        %作第i级小波合成
   pic(1:2*X,1:2*Y)=idwt2(a,h,v,d,'coif5',[2*X,2*Y]);
end