clc;?clear?all;
I=dicomread（‘D:\Image\CT\Z00‘）;
%I=double（I）;
figureimshow（I[]）title（‘原始圖像‘）;
I1=imcrop（I）;?????????????%選取要分割的區(qū)域
figureimshow（I1[]）title（‘已選取‘）;
[yx]=getpts;????????????????%獲得區(qū)域生長(zhǎng)起始點(diǎn)
y1=round（x）;?????????????????%橫坐標(biāo)取整
x1=round（y）;?????????????????%縱坐標(biāo)取整

[MN]=size（I1）;???????????????%獲取圖像大小
stack=[y1x1];???????????????%將生長(zhǎng)起始點(diǎn)灰度值存入stack中
suit=1;??????????????????????%儲(chǔ)存符合區(qū)域生長(zhǎng)條件的點(diǎn)的個(gè)數(shù)
Y=zeros（MN）;????????????????%作一個(gè)全零與原圖像等大的圖像矩陣Y，作為輸出圖像矩陣
Y（y1x1）=1;??????????????????%將種子點(diǎn)的灰度值置1
count=1;?????????????????????%記錄每次判斷一點(diǎn)周圍八點(diǎn)符合條件的新點(diǎn)的數(shù)目
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