minimum-watershed 实现粘连分割

实现粘连分割，对有粘连的物体分割效果很好，使用了局部极值结合分水岭算法 )

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