迭代法求双峰灰度阈值

function x = fenlimili(y)
I1=imread(y);
figure;
subplot(1,2,1);
imshow(I1);%显示原图
title('灰度图像')

T=0.5*(double(min(I1(:)))+double(max(I1(:))));%对阈值进行初始化
d=2;%设置一个初始标度值以便实现循环
%通过迭代求最佳阈值
while d>0.5%设置标度值为0.5
     g=I1>=T;%利用布尔数组将图像分位大于阈值的部分和小于阈值的部分
     Tn=0.5*(mean(I1(g))+mean(I1(~g)));%取大于和小于阈值的部分的平均值的平均值作为新的阈值
     d=abs(T-Tn);%与标度值进行比较，新的阈值是否已经差值很小
     T=Tn;%阈值赋值
end
% 根据最佳阈值进行图像分割
level=Tn/255;
A=im2bw(I1,level);%分割图像
subplot(1,2,2),imshow(A),title(['迭代法分割后阈值为',num2str(Tn)])%同时利用num2str在title中显示数字
end