draw1points1back8.m

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

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

band_save=zeros(16*8,16);
v_my1=zeros(16,16);
v_my2=zeros(16,16);
v_my3=zeros(16,16);
v_my4=zeros(16,16);
v_my5=zeros(16,16);
v_my6=zeros(16,16);
v_my7=zeros(16,16);
v_my8=zeros(16,16);

%根据band(谱段）恢复(与draw1pionts8的抽取顺序相反）
i=fix((band-28)/14)+1;
j=rem((band-28),14)+1;
load v_savetemp1;
load v_savetemp2;
load v_savetemp3;
load v_savetemp4;
load v_savetemp5;
load v_savetemp6;
load v_savetemp7;
load v_savetemp8;

%对第一点进行小波变换
vv=test(v_savetemp1);
v_savetemp1=vv;
%恢复
v_my1=v_savetemp1((i-1)*16+1:i*16,(j-1)*16+1:j*16);
v_my2=v_savetemp2((i-1)*16+1:i*16,(j-1)*16+1:j*16);
v_my3=v_savetemp3((i-1)*16+1:i*16,(j-1)*16+1:j*16);
v_my4=v_savetemp4((i-1)*16+1:i*16,(j-1)*16+1:j*16);
v_my5=v_savetemp5((i-1)*16+1:i*16,(j-1)*16+1:j*16);
v_my6=v_savetemp6((i-1)*16+1:i*16,(j-1)*16+1:j*16);
v_my7=v_savetemp7((i-1)*16+1:i*16,(j-1)*16+1:j*16);
v_my8=v_savetemp8((i-1)*16+1:i*16,(j-1)*16+1:j*16);
for m=1:16
   band_save((m-1)*8+1,:)=v_my1(m,:);
   band_save((m-1)*8+2,:)=v_my2(m,:);
   band_save((m-1)*8+3,:)=v_my3(m,:);
   band_save((m-1)*8+4,:)=v_my4(m,:);
   band_save((m-1)*8+5,:)=v_my5(m,:);
   band_save((m-1)*8+6,:)=v_my6(m,:);
   band_save((m-1)*8+7,:)=v_my7(m,:);
   band_save((m-1)*8+8,:)=v_my8(m,:);

end
%band_save存储经小波零树及非线性量化的结果
save band_save band_save

%进行k-l重建
draw1backtemp(8,band);
return