%2-d 最大熵法(递推方法)
clc;
clear all;
tic
%a=imread('d:\image\x.jpg','jpg');
a=imread('d:\image\alumgrns.tif');
%a=imread('d:\image\airplane2.tiff');
%a=imnoise(a,'gaussian',0,0.005);
%a=imread('d:\image\eight.tif');
%a=imnoise(a,'salt & pepper',0.008);figure,imshow(a);
%a=imnoise(a,'gaussian',0,0.01);
%a = imnoise(a,'speckle',0.09);figure,imshow(a);
%a=imnoise(a,'gaussian',0,0.0018);
%figure,imshow(a);
% a=imread('shot1.tif');
%a=imread('cameraman.tif');
%a=Image;
%clear Image;
subplot(1,2,1);
imshow(a);
a0=double(a);
[m,n]=size(a);
h=1;
a1=zeros(m,n);
% 计算平均领域灰度的一维灰度直方图
for i=1:m
for j=1:n
for k=-h:h
for w=-h:h;
p=i+k;
q=j+w;
if (p<=0)|( p>m)
p=i;
end
if (q<=0)|(q>n)
q=j;
end
a1(i,j)=a0(p,q)+a1(i,j);
end
end
a2(i,j)=uint8(1/9*a1(i,j));
end
end
fxy=zeros(256,256);
% 计算二维直方图
for i=1:m
for j=1:n
c1=a0(i,j);
d=double(a2(i,j));
fxy(c1+1,d+1)=fxy(c1+1,d+1)+1;
end
end
Pxy=fxy/m/n;
% figure,
% mesh(Pxy);
% title('二维灰度直方图');
%计算Hl
Hl=0;
for i=1:256
for j=1:256
if Pxy(i,j)>0.00001
Hl=Hl-Pxy(i,j)*log(Pxy(i,j));
else
Hl=Hl;
end
end
end
%计算PA,HA
PA=zeros(256,256);
HA=zeros(256,256);
PB=zeros(256,256);
PA(1,1)=Pxy(1,1);
if PA(1,1)<1e-4
HA(1,1)=0;
else
HA(1,1)=-PA(1,1)*log(PA(1,1));
end
for i=2:256
PA(i,1)=PA(i-1,1)+Pxy(i,1);
if Pxy(i,1)>0.00001
HA(i,1)=HA(i-1,1)-Pxy(i,1)*log(Pxy(i,1));
else
HA(i,1)=HA(i-1,1);
end
end
for j=2:256
PA(1,j)=PA(1,j-1)+Pxy(1,j);
if Pxy(1,j)>0.00001
HA(1,j)=HA(1,j-1)-Pxy(1,j)*log(Pxy(1,j));
else
HA(1,j)=HA(1,j-1);
end
end
for i=2:256
for j=2:256
PA(i,j)=PA(i-1,j)+PA(i,j-1)-PA(i-1,j-1)+Pxy(i,j);
if Pxy(i,j)>0.00001
HA(i,j)=HA(i-1,j)+HA(i,j-1)-HA(i-1,j-1)-Pxy(i,j)*log(Pxy(i,j));
else
HA(i,j)=HA(i-1,j)+HA(i,j-1)-HA(i-1,j-1);
end
end
end
%计算最大熵
PB=1-PA;
h=zeros(256,256);
hmax=0;
for i=1:256
for j=1:256
if abs(PA(i,j))>0.00001&abs(PB(i,j))>0.00001
h(i,j)=log(PA(i,j)*PB(i,j))+HA(i,j)/PA(i,j)+(Hl-HA(i,j))/PB(i,j);
else
h(i,j)=0;
end
if h(i,j)>hmax
hmax=h(i,j);
s=i-1;
t=j-1;
end
end
end
z=ones(m,n);
for i=1:m
for j=1:n
if a0(i,j)<=s&a2(i,j)<=t
%if double(a(i,j))+double(a2(i,j))+a3(i,j)<=s+t+q
z(i,j)=0;
else
z(i,j)=255;
end
end
end
hmax
s
t
subplot(1,2,2);
imshow(z);
toc
没有合适的资源?快使用搜索试试~ 我知道了~
entropy_max_2d.rar_2d entropy_entropy c++ builder_threshold_二维 阈
共2个文件
txt:1个
m:1个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 125 浏览量
2022-07-14
20:37:47
上传
评论
收藏 1KB RAR 举报
温馨提示
用于实现灰度图像阈值分割的二维最大熵方法,matlab源程序
资源推荐
资源详情
资源评论
收起资源包目录
entropy_max_2d.rar (2个子文件)
entropy_max_2d.m 3KB
www.pudn.com.txt 218B
共 2 条
- 1
资源评论
JonSco
- 粉丝: 66
- 资源: 1万+
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功