matlab区域生长代码修订版

clc; clear all;
I=dicomread（‘D:\Image\CT\Z00‘）;
%I=double（I）;
figureimshow（I[]）title（‘原始图像‘）;
I1=imcrop（I）;             %选取要分割的区域
figureimshow（I1[]）title（‘已选取‘）;
[yx]=getpts;                %获得区域生长起始点
y1=round（x）;                 %横坐标取整
x1=round（y）;                 %纵坐标取整

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