I?=?imread（‘C:\Users\Admin\Desktop\chepail.bmp‘）;
figureimshow（I）
[xy]=size（I）;????????????????%?求出圖象大小
b=double（I）;??????????????????
????????????????????????
z0=max（max（max（I）））;???????????????????%?求出圖象中最大的灰度
z1=min（min（min（I）））;???????????????????%?最小的灰度?

T=（z0+z1）/2;??????????????????????%?T賦初值，為最大值和最小值的平均值
TT=0;?????????????????????????????%TT賦初值
S0=0.0;?n0=0.0;???????????????????????%為計算灰度大于閾值的元素的灰度總值、個數賦值
S1=0.0;?n1=0.0;???????????????????????%為計算灰度小于閾值的元素的灰度總值、個數賦值
allow=50;???????????????????????%?新舊閾值的允許接近程度
d=abs（T-TT）;
count=0;?????????????????????????%?記錄幾次循環
while（d>=allow）?????????????????%?迭代最佳閾值分割算法
????count=count+1;
????for?i=1:x
????????for?j=1:y
????????????if?I（ij）>=T
????????????????S0=S0+b（ij）;?????
????????????????n0=n0+1;
????????????end
????????????if?I（ij）????????????????S1=S1+b（ij）;
????????????????n1=n1+1;
????????????end
????????end
????end?
????T0=S0/n0;
????T1=S1/n1;
????TT=（T0+T1）/2;????????????%在閾值T下，迭代閾值的計算過程
????d=abs（T-TT）;
????T=TT;
end
（T1/T0）

i1=im2bw（IT/255）;??????????%?圖像在最佳閾值下二值化
figureimshow（i1）title（‘最佳閾值‘）;
i2=bwperim（i1）;?????????%?圖像細化過程
figureimshow（i2）?????????????????????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????8210??2019-03-11?19:29??01.bmp
?????文件???????11192??2017-07-06?09:21??02.bmp
?????文件????????1352??2017-07-14?09:54??erzhihuaguocheng.m
?????文件???????20999??2017-07-14?12:52??GUI.fig
?????文件????????8548??2017-07-14?15:53??GUI.m
?????文件??????????42??2019-03-11?21:25??說明.txt