基于区域生长法的图像分割matlab程序

image=imread（‘mri1.jpg‘）;
I=rgb2gray（image）;
figureimshow（I）title（‘原始图像‘）
I=double（I）/255;
[yx]=getpts;                %获得区域生长起始点
y1=round（x）;                 %横坐标取整
x1=round（y）;                 %纵坐标取整

[MN]=size（I）;               %获取图像大小
stack=[y1x1];               %将生长起始点灰度值存入stack中
suit=1;                      %储存符合区域生长条件的点的个数
Y=zeros（MN）;                %作一个全零与原图像等大的图像矩阵Y，作为输出图像矩阵
Y（y1x1）=1;                  %将种子点的灰度值置1
count=1;                     %记录每次判断一点周围八点符合条件的新点的数目
threshold=0.05555;           %阈值
sum=I（y1x1）;                %存灰度值
adaptM=stack（11）;
adaptN=stack（12）;
greyvalue=I（adaptMadaptN）;
while suit>0
    adaptM=stack（11）;
    adaptN=stack（12）;
    
    for u=-1:1                                       %在围围八点找符合条件的点
          for v=-1:1
              if adaptM+u<（M+1） & ada