draw1points16.m

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

function draw1points16
%函数作用为产生k-l变换结果并整理，抽取各小方块的第一点，第二点。。。
%16点取1,2,3点

v16_savetemp1=zeros(32*14,32*14);
v16_savetemp2=zeros(32*14,32*14);
v16_savetemp3=zeros(32*14,32*14);
v16_savetemp4=zeros(32*14,32*14);

v16_my1=zeros(32,32);
v16_my2=zeros(32,32);
v16_my3=zeros(32,32);
v16_my4=zeros(32,32);
%将从28谱段开始的14*14个谱段的k-l变换结果凑成一个大矩阵，以利于小波零树编码
for i=1:14
   for j=1:14
      %drawtemp4为产生保留n点的k-l变换结果的函数
drawtemp4(4,28+(i-1)*14+j-1);
%v_save保存n点的k-l变换结果
load v_save4;
for m=1:32
   v16_my1(m,:)=v_save4((m-1)*4+1,:);
   v16_my2(m,:)=v_save4((m-1)*4+2,:);
   v16_my3(m,:)=v_save4((m-1)*4+3,:);
   v16_my4(m,:)=v_save4((m-1)*4+4,:);
   
end   
v16_savetemp1((i-1)*32+1:i*32,(j-1)*32+1:j*32)=v16_my1;
v16_savetemp2((i-1)*32+1:i*32,(j-1)*32+1:j*32)=v16_my2;
v16_savetemp3((i-1)*32+1:i*32,(j-1)*32+1:j*32)=v16_my3;
v16_savetemp4((i-1)*32+1:i*32,(j-1)*32+1:j*32)=v16_my4;

end
end
%分别存储保留4点的第一点，第二点。。。
save v16_savetemp1 v16_savetemp1;
save v16_savetemp2 v16_savetemp2;
save v16_savetemp3 v16_savetemp3;
save v16_savetemp4 v16_savetemp4;
return