数字图像处理——图像分割

数字图像处理——图像分割

迭代阈值分割

clear,clc;
Image = double(imread('迭代阈值分割.jpg'));
Image = Image(:,:,1);
[lines,samples] = size(Image);
%迭代阈值影像分割算法
T = median( reshape(Image,[1,lines*samples]));
T0 = inf;
while( abs(T0 - T) > 1e-2)
    T0 = T;
    G1 = Image>=T;
    K1 = G1.*Image;
    G2 = Image