qiushang.m

我自己编的图像求熵的程序。运行经过了优化。

A=imread('LENA_G.bmp'); 

[M,N]=size(A); 
temp=zeros(1,256); 

%对图像的灰度值在[0,255]上做统计 
for m=1:M; 
for n=1:N; 


i=A(m,n)+1; 

temp(i)=temp(i)+1; 
end 
end 
temp=temp./(M*N); 

%由熵的定义做计算 
result=0; 


for i=1:256 
if temp(i)==0; 
result=result; 
else 
result=result-temp(i).*log2(temp(i)); 
end 
end 
result