im=imread('im2.jpg');
figure;imshow(im,[]);title('原图');
I = double(im);
Ir = I(:,:,1)-I(:,:,2); %红色减去绿色
figure;imshow(Ir,[]);title('红色减去绿色');
Irmask = Ir > 20; %提取二值化的红色小棍信息
Igmask = Ir < -50;
figure;imshow(Irmask,[]);title('红色小棍二值图');
figure;imshow(Igmask,[]);title('绿色小棍二值图');
bwr=bwareaopen(Irmask,500);
figure;imshow(bwr);
bwg=bwareaopen(Igmask,500,4);
[L,num]=bwlabel(bwr,4);
[L1,num1]=bwlabel(bwg,4);
aa=(L1==3);
figure;imshow(aa);
m1=num;
m=zeros(2,num);
for i=1:m1
% 计算目标区域中心，用于显示编号的位置
[p,q]=find(L==i);
temp=[p,q];
[x,y]=size(temp);
m(1,i)=sum(p)/x;
m(2,i)=sum(q)/x;
end
for i=1:m1
figure(1);
text(m(2,i),m(1,i),int2str(i),'color','green')
end
m3=num1;%ahe1(1,1);
m2=zeros(2,num1);
for i=1:m3
% 计算目标区域中心，用于显示编号的位置
[p,q]=find(L1==i);
temp=[p,q];
[x,y]=size(temp);
m2(1,i)=sum(p)/x;
m2(2,i)=sum(q)/x;
end
for i=1:m3
figure(1);
text(m2(2,i),m2(1,i),int2str(i),'color','red');
end
