draw1points8.m

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

function draw1points8
%函数作用为产生k-l变换结果并整理，抽取各小方块的第一点，第二点。。。
%64点取8点

v_savetemp1=zeros(16*14,16*14);
v_savetemp2=zeros(16*14,16*14);
v_savetemp3=zeros(16*14,16*14);
v_savetemp4=zeros(16*14,16*14);
v_savetemp5=zeros(16*14,16*14);
v_savetemp6=zeros(16*14,16*14);
v_savetemp7=zeros(16*14,16*14);
v_savetemp8=zeros(16*14,16*14);

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);
%将从28谱段开始的14*14个谱段的k-l变换结果凑成一个大矩阵，以利于小波零树编码
for i=1:14
   for j=1:14
      %drawtemp为产生保留n点的k-l变换结果的函数
      drawtemp(8,28+(i-1)*14+j-1);
%v_save保存n点的k-l变换结果
load v_save;

for m=1:16
   v_my1(m,:)=v_save((m-1)*8+1,:);
   v_my2(m,:)=v_save((m-1)*8+2,:);
   v_my3(m,:)=v_save((m-1)*8+3,:);
   v_my4(m,:)=v_save((m-1)*8+4,:);
   v_my5(m,:)=v_save((m-1)*8+5,:);
   v_my6(m,:)=v_save((m-1)*8+6,:);
   v_my7(m,:)=v_save((m-1)*8+7,:);
   v_my8(m,:)=v_save((m-1)*8+8,:);

end   
v_savetemp1((i-1)*16+1:i*16,(j-1)*16+1:j*16)=v_my1;
v_savetemp2((i-1)*16+1:i*16,(j-1)*16+1:j*16)=v_my2;
v_savetemp3((i-1)*16+1:i*16,(j-1)*16+1:j*16)=v_my3;
v_savetemp4((i-1)*16+1:i*16,(j-1)*16+1:j*16)=v_my4;
v_savetemp5((i-1)*16+1:i*16,(j-1)*16+1:j*16)=v_my5;
v_savetemp6((i-1)*16+1:i*16,(j-1)*16+1:j*16)=v_my6;
v_savetemp7((i-1)*16+1:i*16,(j-1)*16+1:j*16)=v_my7;
v_savetemp8((i-1)*16+1:i*16,(j-1)*16+1:j*16)=v_my8;
end
end
%分别存储保留8点的第一点，第二点。。。
save v_savetemp1 v_savetemp1;
save v_savetemp2 v_savetemp2;
save v_savetemp3 v_savetemp3;
save v_savetemp4 v_savetemp4;
save v_savetemp5 v_savetemp5;
save v_savetemp6 v_savetemp6;
save v_savetemp7 v_savetemp7;
save v_savetemp8 v_savetemp8;

return