J=?imread（‘RoadG1.jpg‘）;
[xy]=size（J）;
I=double（J）;
z0=max（max（I））;???????????????????%?求出圖象中最大的灰度
z1=min（min（I））;???????????????????%?最小的灰度?

T=（z0+z1）/2;??????????????????????
TT=0;
S0=0;?n0=0;
S1=0;?n1=0;
allow=0.5;???????????????????????%?新舊閾值的允許接近程度
d=abs（T-TT）;
count=0;?????????????????????????%?記錄幾次循環

while（d>=allow）?????????????????%?迭代最佳閾值分割算法
????count=count+1;
????for?i=1:x
????????for?j=1:y
????????????if?（I（ij）>=T）
????????????????S0=S0+I（ij）;
????????????????n0=n0+1;
????????????end
????????????if?（I（ij）????????????????S1=S1+I（ij）;
????????????????n1=n1+1;
????????????end
????????end
????end?
????T0=S0/n0;
????T1=S1/n1;
????TT=（T0+T1）/2;
????d=abs（T-TT）;
????T=TT;
end

Seg=zeros（xy）;
for?i=1:x
????for?j=1:y
????????if（I（ij）>=T）
????????????Seg（ij）=1;???????????????%?閾值分割的圖象
????????end
????end
end

J0=Seg;
SE=strel（‘rectangle‘?[6?4]）;?????????%?結構定義
IM=imopen（J0SE）;??????????????%??“開”運算

figureimshow（J）;????%?圖9-7（a）
figureimhist（J）;?????%?圖9-7（b）
figureimshow（Seg）;??%?圖9-7（c）
figureimshow（IM）;??%?圖9-7（d）

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

?????文件???????1169??2004-04-04?18:50??chap9\P0901.m

?????文件????????754??2004-04-04?18:52??chap9\P0902.m

?????文件???????1497??2004-04-04?18:53??chap9\P0903.m

?????文件????????719??2004-04-04?18:58??chap9\P0904.m

?????文件????????185??2004-04-04?19:04??chap9\P0905.m

?????文件????????165??2004-07-21?13:46??chap9\readme.txt

?????文件??????55560??2004-01-12?19:27??chap9\RoadG1.jpg

?????文件??????36241??2003-12-15?14:18??chap9\RoadG2.jpg

?????文件??????36241??2003-12-15?14:18??chap9\z12.jpg

?????文件???????7303??2004-01-15?14:48??chap9\z12R.jpg

????...D..R?????????0??2009-07-03?11:50??chap9

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

???????????????139834????????????????????11