matlab实现最大熵法图像分割程序源代码.doc

《matlab实现最大熵法图像分割程序源代码.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、clearalla=imread('moon.tif');figure,imshow(a)count=imhist(a);[m,n]=size(a);N=m*n;L=256;count=count/N;%%每一个像素的分布概率countfori=1:Lifcount(i)~=0st=i-1;break;endendstfori=L:-1:1ifcount(i)~=0nd=i-1;break;endendndf=count(st+1:nd+1);%f是每个灰度出现的概率size(f)E=[];forTh=st:nd-1%%%设定初始分割阈值

2、为Thav1=0;av2=0;Pth=sum(count(1:Th+1));%%%第一类的平均相对熵为fori=0:Thav1=av1-count(i+1)/Pth*log(count(i+1)/Pth+0.00001);end%%%第二类的平均相对熵为fori=Th+1:L-1av2=av2-count(i+1)/(1-Pth)*log(count(i+1)/(1-Pth)+0.00001);endE(Th-st+1)=av1+av2;endposition=find(E==(max(E)));th=st+position-1fori=

3、1:mforj=1:nifa(i,j)>tha(i,j)=255;elsea(i,j)=0;endendendfigure,imshow(a);