tic
I=imread（‘C:\Users\zdh\Desktop\bishe\image\9.jpg‘）;
%?figureimshow（I）title（‘原圖像‘）;
I=rgb2gray（I）;
I1=double（I）;
T=（min（I1（:））+max（I1（:）））/2;
d=false;
i=0;
while~d
????u1=find（I1<=T）;
????u2=find（I1>T）;
????Tn=（mean（I1（u1））+mean（I1（u2）））/2;
????d=abs（Tn-T）<1;
????T=Tn;
????i=i+1;
end
I1（u1）=0;
I1（u2）=1;
%?figureimshow（I1）title（‘二值圖像‘）;
se=strel（‘disk‘18）;
A=imerode（I1se）