数字图像处理实验十二迭代法阈值分割

利用迭代法对图像进行阈值分割。

clear;
clc;
%----------------------------图像的灰度化处理-------------------------------
JPG=imread('1.jpg');
JPG=rgb2gray(JPG);
figure();
imshow(JPG);title('原始图像')

%-----------------------------迭代法阈值分割--------------------------------
T=mean2(JPG);   %取均值作为初始阈值
i=0;    %记录迭代次数
% while循环进行迭代
while 1
    i=i+1;  %记录迭代次数
    R1=find(JPG<=T);  %小于阈值的部分
    R2=find(JPG>T);   %大于阈值的部分
    T1=(mean(JPG(R1))+mean(JPG(R2)))/2;  %计算分割后两部分的阈值均值的均值
    if abs(T1-T)<0.1  %两个平均灰度基本无差别时退出
        T=T1;
        break;
    end    
    T=T1;      %如果差别较大,则将分割后的均值作为新的阈值再次进行运算
end
JPG(R1)=0;   %将小于阈值的部分赋值为0
JPG(R2)=1;   %将大于阈值的部分赋值为1   这两步是将图像转换成二值图像
figure();
imshow(JPG,[]);title('迭代法阈值分割图像')

效果图展示