基于种子点8-邻域的区间生长算法的Matlab实现

function J = regionGrow（I）
% 区域生长，需要以交互方式设定初始种子点，具体方法为鼠标单击图像中一点后，按下回车键
%
% 输入：I - 原图像
% 输出：J - 输出图像

if isinteger（I）
    I=im2double（I）;
end
figureimshow（I）title（‘原始图像‘）
[MN]=size（I）;
[yx]=getpts;             %获得区域生长起始点
x1=round（x）;            %横坐标取整
y1=round（y）;            %纵坐标取整
seed=I（x1y1）;           %将生长起始点灰度值存入seed中
J=zeros（MN）;          %作一个全零与原图像等大的图像矩阵J，作为输出图像矩阵
J（x1y1）=1;             %将J中与所取点相对应位置的点设置为白
sum=seed;              %储存符合区域生长条件的点的灰度值的和
suit=1;                 %储存符合区域生长条件的点的个数
count=1;               %记录每次判断一点周围八点符合条件的新点的数目
threshold=0.15;         %阈值，注意需要和double类型存储的图像相符合
while count>0
    s=0;                   %记录判断一点周围八点时，符合条件的新点的灰度值之和
     count=0;
     for i=1:M
       for j=1:N
         if J（ij）==1
          if （i-1）>0 & （i+1