落花成冢
matlab练习程序（图像局部熵）
⼜被国内论⽂坑了，话说写论⽂要有良⼼。
今天在看局部熵⽅⾯的内容，看论⽂中介绍的内容感觉局部熵挺容易了，于是就有了实现的想法，结果效果⾮常糟糕。你若离去后会无期
得到的⼏乎是⼀张空⽩的图⽚，就像下⾯⼀样：
究其原因是各种论⽂上都写了这样⼀个公式：
这⾥f(i,j)就是在m*n这个局部的像素，这个没问题，不过这⾥的p是什么东西，按这⾥的定义p是当前像素灰度占局部总灰度的概率，⽽p原本应该是局部直⽅图，也就是当前像素灰度个数占局部像素总个数的概率，所以这⾥的意义我也不明⽩了，结果按公式计算就得到了上图，和各种论⽂中的图都不⼀样。
⽽真正的公式应该是这样：
仓鼠的习性熵的公式本来不就是这样么，这⾥的p是归⼀化后的局部直⽅图。
这样运⾏的结果就正确了，和matlab系统提供的entropyfilt函数结果是⼀样的。
⾃⼰函数处理的
系统函数处理的
ppic9k
matlab代码如下：
clear all;
clo all;
clc;
img=imread('lena.jpg');
[m n]=size(img);
w=3; %模板半径
imgn=zeros(m,n);
for i=1+w:m-w
for j=1+w:n-w
Hist=zeros(1,256);
for p=i-w:i+w
for q=j-w:j+w
Hist(img(p,q)+1)=Hist(img(p,q)+1)+1; %统计局部直⽅图end
end
Hist=Hist/sum(Hist);
报春儿歌如何成为党员for k=1:256
if Hist(k)~=0
imgn(i,j)=imgn(i,j)+Hist(k)*log(1/Hist(k)); %局部熵end
仙女湖景区
end
%{
p=sum(sum(img(i-w:i+2,j-w:j+w))); %这⾥是按第⼀个公式写的 s=img(i-w:i+w,j-w:j+w)/p;
imgn(i,j)=-sum(sum(s.*log(s)));
%}
end
end
文明礼仪班会教案imshow(imgn,[])
imgn=entropyfilt(img); %系统的局部熵函数
figure;
imshow(imgn,[])
