T=0.5*(double(min(d(:)))+double(max(d(:)))); % 计算初始阈值done=false;while~done g=d>=T %阈值比较 Tnext=0.5*(mean(d(g))+mean(d(~g)));%计算当前阈值下,分割后图像均值,并作为下一次的阈值 done=abs(T-Tnext)<0.5;%判断是否可以终止分割 T=Tnext; end