ÖÂл
ͼÏñÐÔÄÜÆÀ¼ÛÖ¸±ê´úÂ룺
[X]=imread('D:\\MATLAB\\work\\image\\¾í»ýÔͼ.jpg');%×°ÔØÔͼµÄͼÏñ [Y]=imread('D:\\MATLAB\\work\\image\\¾í»ýÈ¥Ôëͼ.jpg');%×°ÔØ´¦ÀíºóµÄͼÏñ
mseij=double(zeros(240,320));%¼Ù¶¨Ëù´¦ÀíµÄͼÏñΪ8±ÈÌØµÄ»Ò¶ÈͼÏñ£¬L=255 for i=1:240 for j=1:320
mseij(i,j)=mseij(i,j)+(X(i,j)-Y(i,j))*(X(i,j)-Y(i,j)); mseij(i,j)=mseij(i,j)+(X(i,j)^2-Y(i,j)^2) end end
mse=(sum((mseij(:)))/(240*320))
psnr=10*log10(255*255/mse) %·åÖµÐÅÔë±È t(1)=0;%±£´æ¸÷¸ö»Ò¶ÈÖµ³öÏÖ¸ÅÂÊ for i=1:256 t(i)=0; end
for i=1:240 for j=1:320
t(X(i,j)+1)=t(X(i,j)+1)+1; end end
%¼ÆËã¸÷¸ö»Ò¶ÈÖµ³öÏֵĸÅÂÊ for i=1:256
t(i)=t(i)/(240*320); end
entropy=0; for i=1:256 if t(i)>0
entropy=entropy-t(i)*log2(t(i)); end
end %Ôö¼Ó entropy %
39
Ïà¹ØÍÆ¼ö£º