findmy.m

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

function ratefind=findmy
load v_savetemp3;
load v_savetemp2;
load v_savetemp4;
load v_savetemp5;
load v_savetemp6;
load v_savetemp7;
load v_savetemp8;
v_temp=zeros(224*7,224);
v_temp(1:224,:)=v_savetemp2;
v_temp(225:224*2,:)=v_savetemp3;
v_temp(224*2+1:224*3,:)=v_savetemp4;
v_temp(224*3+1:224*4,:)=v_savetemp5;
v_temp(224*4+1:224*5,:)=v_savetemp6;
v_temp(224*5+1:224*6,:)=v_savetemp7;
v_temp(224*6+1:224*7,:)=v_savetemp8;

pic=wav_de3(v_temp,3);
themax=max(max(pic));
v_t=pic./themax;
x=0;
y=0;
z=0;
for i=1:224*3
   for j=1:224
      if(v_t(i,j)>0.5)
         x=x+1;
      end
      if(v_t(i,j)>0.25)&(v_t(i,j)<0.5)
         y=y+1;
      end
      if(v_t(i,j)<=0.25)
         z=z+1;
      end
   end
end
disp(['>0.5        ',int2str(x)]);
disp(['<0.5>0.25   ',int2str(y)]);
disp(['<0.25       ',int2str(z)]);
%x+y+z-224*224*3
format short
ratefind=(224*224*7*32)/(x*6+y*4+z*2);
return