阈值分割

一、迭代法
通过迭代法来求出分割的最佳阈值，具有一定的自适应性，这种方法适用于物体的和背景有较大的区分，即其直方图有相当清楚的波谷；

具体算法步骤如下：
1.为全局阈值选择一个初始估计值T。
2.用T分割图像。产生两组像素：G1有灰度值大于T的像素组成，G2有小于等于T像素组成。
3.计算G1和G2像素的平均灰度值m1和m2；
4.计算一个新的阈值:T = (m1 + m2) / 2;
5.重复步骤2和4,直到连续迭代中的T值间的差小于一个预定义参数为止。

实现代码：

clear;
clc;
close all;
I=imread('Fig1013(a).tif');
T=0.01;
I=im2double(I);
% T1=(min(I(:))+max(I(:)))/2;
T1=mean2(I);
index=find(I<T1);
index2=find(I>&#