draw1points1back16.m

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

function draw1pointsback16(n,band)
%函数作用为对k-l变换结果进行小波零树编码及非线性量化编码

band16_save=zeros(32*n,32);
v16_my1=zeros(32,32);
v16_my2=zeros(32,32);
v16_my3=zeros(32,32);
v16_my4=zeros(32,32);
%根据band(谱段）恢复(与draw1pionts16的抽取顺序相反）
i=fix((band-28)/14)+1;
j=rem((band-28),14)+1;

load v16_savetemp1;
load v16_savetemp2;
load v16_savetemp3;
load v16_savetemp4;

%对第一点进行小波变换
vv=test(v16_savetemp1);
v16_savetemp1=vv;
%恢复
v16_my1=v16_savetemp1((i-1)*32+1:i*32,(j-1)*32+1:j*32);%=v_my1;%reshape(v_save,16,16);
v16_my2=v16_savetemp2((i-1)*32+1:i*32,(j-1)*32+1:j*32);%=v_my2;%reshape(v_save,16,16);
v16_my3=v16_savetemp3((i-1)*32+1:i*32,(j-1)*32+1:j*32);%=v_my3;%reshape(v_save,16,16);
v16_my4=v16_savetemp4((i-1)*32+1:i*32,(j-1)*32+1:j*32);%=v_my4;%reshape(v_save,16,16);
%根据n进行存储
if(n==1)
for m=1:32
   band16_save((m-1)*1+1,:)=v16_my1(m,:);
   %16band_save((m-1)*4+2,:)=v_my2(m,:);
   %16band_save((m-1)*4+3,:)=v_my3(m,:);
   %16band_save((m-1)*4+4,:)=v_my4(m,:);
end
%whos
%return
end
if(n==2)
for m=1:32
   band16_save((m-1)*2+1,:)=v16_my1(m,:);
   band16_save((m-1)*2+2,:)=v16_my2(m,:);
   %16band_save((m-1)*4+3,:)=v_my3(m,:);
   %16band_save((m-1)*4+4,:)=v_my4(m,:);
end
end
if(n==3)
for m=1:32
   band16_save((m-1)*3+1,:)=v16_my1(m,:);
   band16_save((m-1)*3+2,:)=v16_my2(m,:);
   band16_save((m-1)*3+3,:)=v16_my3(m,:);
   %16band_save((m-1)*4+4,:)=v_my4(m,:);
end
end

if(n==4)
for m=1:32
   band16_save((m-1)*4+1,:)=v16_my1(m,:);
   band16_save((m-1)*4+2,:)=v16_my2(m,:);
   band16_save((m-1)*4+3,:)=v16_my3(m,:);
   band16_save((m-1)*4+4,:)=v16_my4(m,:);
end
end


%band16_save存储经小波零树及非线性量化的结果
save band16_save band16_save
%进行k-l重建
drawbacktemp4(n,band);
return
%----------------------------------