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