I=imread（‘1021.jpg‘）;
imshow（I）
r=I（::1）;
g=I（::2）;
b=I（::3）;
a0=0.27*r+0.57*g+0.06*b;
amax=max（a0）;
amin=min（a0）;
%figureimhist（a0）
a1=im2bw（a069/255）;
figureimshow（a1）
a2=imfill（a1‘holes‘）;
%figureimshow（a2）

se?=?strel（‘disk‘3）;
Image?=?imopen（a2se）;
%figureimshow（Image）title（‘預處理后圖‘）

I2=medfilt2（Image[3?3]）;
%figureimshow（I2）title（‘濾波后圖‘）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
I3=imerode（I2se）;
I3=imerode（I3se）;
%I3=imerode（I3se）;

%figureimshow（I3）title（‘腐蝕后圖‘）
I4=imdilate（I3se）;
I4=imdilate（I4se）;
%I4=imdilate（I4se）;

figureimshow（I4）title（‘膨脹后圖‘）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
L0=?bwlabel（I4）;%%%
Ar=zeros（1max（L0（:）））;
for?i=1:max（L0（:））
????Ar（i）=sum（I4（L0==i））;
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-08-26?16:07??極值結合分水嶺-大螟\
?????文件??????147998??2011-09-02?13:22??極值結合分水嶺-大螟\1021.jpg
?????文件????????1773??2011-09-03?15:50??極值結合分水嶺-大螟\jubujizhi.m
?????文件???????42988??2011-09-02?13:38??極值結合分水嶺-大螟\二值化圖.jpg
?????文件???????25045??2011-09-03?15:50??極值結合分水嶺-大螟\原圖分割后結果.jpg
?????文件???????22641??2011-09-03?15:50??極值結合分水嶺-大螟\最終分割結果.jpg
?????文件???????13608??2011-09-03?15:50??極值結合分水嶺-大螟\粘連分割后結果.jpg