[filename?pathname]?=?uigetfile（‘*.bmp;*.pgm;*.tif;*.jpg;*.*‘?‘Load?input?image‘）;
?filename?=?[pathname?filename];?
?im_original=imread（filename）;
?[nxnyl]=size（im_original）;?
?if?l>1
?????image_I=rgb2gray（im_original）;%將原圖象轉(zhuǎn)換為灰度圖象
?else
?????image_I=?im_original;
?end
?figure（1）;imshow（image_I）;title（‘原圖像‘）
?



?image_I=double（image_I）;
?image_I12=image_I;
??for?i=1:3
?image_I1=imfilter（image_Ifspecial（‘gaussian‘102））;
??end
?figure（3）;imshow（uint8（image_I1））;title（‘lvbo圖像‘）?
?

im_graylap=image_I1;
T1=0.5*（double（min（im_graylap（:）））+double（max（im_graylap（:））））;
done=0