function?J?=?regionGrow（I）
%?區(qū)域生長(zhǎng)，需要以交互方式設(shè)定初始種子點(diǎn)，具體方法為鼠標(biāo)單擊圖像中一點(diǎn)后，按下回車(chē)鍵
%
%?輸入：I?-?原圖像
%?輸出：J?-?輸出圖像
clc;
clear;
?I=imread（‘colorpic398.jpg‘）;

?I1=rgb2gray（I）;
?thresh=graythresh（I1）;
?
%???thresh?=?0.5?*?（double（min（I（:）））?+?double（max（I（:））））;?%初始閾值
%???????done?=?false;?%結(jié)束標(biāo)志
%???????while?~done
%?	??????g?=?I?>=?thresh;
%?	??????Tnext?=?0.5?*?（mean（I（g））?+?mean（I（~g）））;
%?	??????done?=?abs（thresh?-?Tnext）?%?	??????thresh?=?Tnext;
%???????end;
??????
if??isinteger（I）;???%判斷給定數(shù)組是否是整型數(shù)組
????I=im2double（I）;???%把圖像數(shù)據(jù)類(lèi)型轉(zhuǎn)換為雙精度浮點(diǎn)類(lèi)型
end
figureimshow（I）title（‘原始圖像‘）
[MN]=size（I）;
[yx]=getpts;?????????????%獲得區(qū)域生長(zhǎng)起始點(diǎn)
x1=round（x）;????????????%橫坐標(biāo)取整
y1=round（y）;????????????%縱坐標(biāo)取整
seed=I（230270）;???????????%將生長(zhǎng)起始點(diǎn)灰度值存入seed中
J=zeros（MN）;?????%作一個(gè)全零與原圖像等大的圖像矩陣