基于matlab的区域生长算法实现

利用区域生长算法进行图像分割

算法步骤

1. 鼠标点击选取生长点
2. 3*3 8-邻域生长
3. 生长准则：图像与生长点灰度差小于阈值
4. 生长点压入栈中，直到栈中没有元素时停止

I=imread('mountain.png');

if isinteger(I)
    I=im2double(I);
end
I = rgb2gray(I);

figure 
imshow(I)
[M,N]=size(I);
[y,x]=getpts; %单击取点后，按enter结束
x1=round(x);
y1=round(y);
seed=I(x1,y1); %获取中心像素灰度值

J=zeros(M,N);
J(x1,y1)=1;

count=1; %待处理点个数
threshold=0.15;
while count>0
    count=0;
    for i=1:M %遍历整幅图像
    for j=1:N
        if J(i,j)==1 %点在“栈”内
        if (i-1)>1&(i+1)<M&(j-1)>1&(j+1)<N %3*3邻域在图像范围内
            for u=-1:1 %8-邻域生长
            for v=-1:1
                if J(i+u,j+v)==0&abs(I(i+u,j+v)-seed)<=3.75&abs(I(i+u,j+v)-seed)<=threshold
                    J(i+u,j+v)=1;
                    count=count+1;  %记录此次新生长的点个数
                end
            end
            end
        end
        end
    end
    end
end

subplot(1,2,1),imshow(I);
title("original image")
subplot(1,2,2),imshow(J);
title("segmented image")

处理结果

参考
https://blog.csdn.net/qq_35608277/article/details/78598016