I?=?medfilt2（imread（‘1.jpg‘））;
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;???????????