基于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://www.geek-share.com/detail/2722275890.html