基于区域生长的MATLAB实现

image=imread（‘C:\Users\yin\Desktop\毕业设计图像\1.jpg‘）; %读取图1-1（c）
I=rgb2gray（image）; %将图像转为灰度图像
figureimshow（I）; %显示灰度图像
I=double（I）/255; %将数据类型转为双精度类型
[MN]=size（I）; %计算数据大小w
[yx]=getpts; %获得区域生长起始点
x1=round（x）; %横坐标取整
y1=round（y）; %纵坐标取整
seed=I（x1y1）; %将生长起始点灰度值存入seed中
Y=zeros（MN）; %作一个全零与原图像等大的图像矩阵Y，作为输出图像矩阵
Y（x1y1）=1; %将Y中与所取点相对应位置的点设置为白场
sum=seed; %储存符合区域生长条件的点的灰度值的和
suit=1; %储存符合区域生长条件的点的个数
count=1; %记录每次判断一点周围八点符合条件的新点的数目
threshold=0.33; %全局阈值
while count>0
 s=0; %记录判断一点周围八点时，符合条件的新点的灰度值之和
 count=0;
 for i=1:M
  for j=1:N
   if Y（ij）==1
    if （i-1）>0 && （i+1）<（M+1） && （j-1）>0 && （j+1）<（N+1）; %判断此点是否为图像边界	上的点
     for u= -1:1  %判断点周围八点是否符合阈值条件
      for v= -1:1  %uv为偏移量
       if Y（i+uj+v）==0 && abs（I（i+uj+v）-