image=imread（‘mri1.bmp‘）;
I=rgb2gray（image）;
figureimshow（I）title（‘原始圖像‘）
I=double（I）;
[MN]=size（I）;
[yx]=getpts;?????????????%獲得區域生長起始點
x1=round（x）;????????????%橫坐標取整
y1=round（y）;????????????%縱坐標取整
seed=I（x1y1）;???????????%將生長起始點灰度值存入seed中
Y=zeros（MN）;??????????%作一個全零與原圖像等大的圖像矩陣Y，作為輸出圖像矩陣
Y（x1y1）=1;?????????????%將Y中與所取點相對應位置的點設置為白場
sum=seed;??????????????%儲存符合區域生長條件的點的灰度值的和
suit=1;?????????????????%儲存符合區域生長條件的點的個數
count=1;???????????????%記錄每次判斷一點周圍八點符合條件的新點的數目
threshold=15;???????%域值
while?count>0
?s=0;??????????????????%記錄判斷一點周圍八點時，符合條件的新點的灰度值之和
?count=0;
?for?i=1:M
???for?j=1:N
?????if?Y（ij）==1
??????if?（i-1）>0?&?（i+1）<（M+1）?&?（j-1）>0?&?（j+1）<（N+1）?%判斷此點是否為圖像邊界上的點
???????for?u=?-1:1??????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1758??2008-06-11?18:49??quyushengzhangfa.m

-----------?---------??----------?-----??----

?????????????????1758????????????????????1