entropy_calc.m

matlab 关于二值图像熵和交叉熵的计算程序,这一程序与结合其它方法对图像进行分割,作为一个判断标准.

%计算一副图像的熵 
A=imread('lena1.bmp');
[M,N]=size(A); 
temp=zeros(1,256); 

%对图像的灰度值在[0,255]上做统计 
for m=1:M; 
    for n=1:N; 
        if A(m,n)==0; 
            i=1; 
        else 
            i=A(m,n); 
        end 
            temp(i)=temp(i)+1; 
    end 
end 
temp=temp./(M*N); 

%由熵的定义做计算 
result=0; 
for i=1:length(temp) 
    if temp(i)==0; 
        result=result; 
    else 
        result=result-temp(i)*log2(temp(i)); 
    end 
end 
result