%一維最大熵法閾值分割?

I?=?imread（‘block.jpg‘）;??
E?=?rgb2gray（I）;??
vHist=imhist（E）;??????????????????????????????????????%得到灰度直方圖
[mn]=size（E）;??
p=vHist（find（vHist>0））/（m*n）;????????????????%求每一不為零的灰度值的概率??

Pt=cumsum（p）;??????????????????????????????????????%計算出選擇不同t值時，A區(qū)域的概率??

Ht=-cumsum（p.*log（p））;????????????????????????%計算出選擇不同t值時，A區(qū)域的熵??

HL=-sum（p.*log（p））;??????????????????????????????%計算出全圖的熵??

Yt=log（Pt.*（1-Pt）+eps）+Ht./（Pt+eps）+（HL-Ht）./（1-Pt+eps）;????????????????????????????%計算出選擇不同t值時，判別函數的值??

[ath]=max（Yt）;???????????????????????????????????????%?th即為最佳閾值??
segImg=（E>th）;??

subplot（121）;
imshow（I）?;
xlabel（‘原圖‘）;
subplot（122）;
imshow（segImg）?;
xlabel（‘最大熵分割法‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????807??2016-12-19?21:23??maximum_entropy.m