findmy16.m

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

function ratefind=findmy
load v16_savetemp2;
v_temp=v16_savetemp2;
pic=wav_de3(v_temp,3);
themax=max(max(pic));
v_t=pic./themax;
x=0;
y=0;
z=0;
for i=1:224
   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