Question

怎么用matlab计算图像熵

Answer 1

%---------------------------------
%求一幅数字图像的熵值
%---------------------------------
I=imread('lena.bmp');
%I=double(I);
[C,R]=size(I); %求图像的规格
Img_size=C*R; %图像像素点的总个数
L=256; %图像的灰度级
H_img=0;
nk=zeros(L,1);
for i=1:C
for j=1:R
Img_level=I(i,j)+1; %获取图像的灰度级
nk(Img_level)=nk(Img_level)+1; %统计每个灰度级像素的点数
end
end
for k=1:L
Ps(k)=nk(k)/Img_size; %计算每一个灰度级像素点所占的概率
if Ps(k)~=0; %去掉概率为0的像素点
H_img=-Ps(k)*log2(Ps(k))+H_img; %求熵值的公式
end
end
H_img
%entropy(I)

Answer 2

　　function entr=yentropy(a)
　　a=uint8(a); %这里a为8位的单色图像或24为的RGB彩色图像
　　[m n l]=size(a);
　　entr=0;
　　for k=1:l
　　hi=zeros(1,256);
　　for i=1:m
　　for j=1:n
　　hi(a(i,j,k)+1)=hi(a(i,j,k)+1)+1; %求每种值的在图像中出现的次数
　　end
　　end
　　hi=sort(hi,'descend');
　　hi=hi./m./n; %求概率
　　en=0.0;
　　for i=1:256
　　if hi(i)>0
　　en=en-hi(i).*log2(hi(i)); %概率不为0 累加求熵
　　else
　　i=257; %否则停止
　　end
　　end
　　entr=entr+en;
　　end
　　entr=entr/l; %当l=1时a为单色图像；
　　当l=3时a为彩色图像，三个页面的熵平均.