matlab entropy 函数,Matlab entropyfilt函数 解析

J = entropyfilt(I,NHOOD)

该函数用来计算灰度影像的局部熵。其中NHOOD为邻域，默认邻域为9*9。

邻域设置越大，计算时间越久。图像熵表示为图像灰度级集合的比特平均数，单位比特/像素，也描述了图像信源的平均信息量。

H(p)=-∑i,jp(i.j)lnp(i,j), 其中p(i,j)=x(i,j)∑i,jx(i,j),x(i,j)为图像的像元。

For pixels on the borders of I, ENTROPYFILT uses symmetric padding. 在边界区域ENTROPYFILT函数使用对称填充。下面的matlab代码没有在边界区域使用对称填充。

[m,n]=size(I);

In=zeros(size(I));%In为影像局部熵

k=25;%k为搜索半径

for i=k+1:m-k

for j=k+1:n-k

his=zeros(1,256);

for p=i-k:i+k

for q=j-k:j+k

his(I(p,q)+1)=his(I(p,q)+1)+1;%计算2k+1块的灰度直方图

end

end

his=his/sum(his);

for g=1:256

if his(g)~=0

In(i,j)=In(i,j)-his(g)*log(his(g));

end

end

end

end