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