A= imread('d:\ronghepc.tif');
[M,N]=size(A);
pa=zeros(1,256);
pa=double(pa);
pc=zeros(1,256);
A=uint8(A);
for m=1:M;
for n=1:N;
if A(m,n)==0;
k1=1;
else
k1=A(m,n);
end
pa(k1)=pa(k1)+1;
end
end
pa=pa/(M*N);
ENTROPHY1=0;
M=256;
for i=1:M;
if pa(i)==0;
ENTROPHY1=ENTROPHY1;
else
ENTROPHY1=ENTROPHY1-pa(i)*log2(pa(i));
end
end
h=ENTROPHY1;