我有一张这样的图片，如何提取里面的红色圈圈坐标，并且连接这些坐标形成两个封闭的环路？我有很多类似的图片，需要批处理

Question

下面是给出的图片，我还有很多类似的图片，处理起来需要返回圆圈的像素坐标，最后连接封闭图形的顺序坐标

Answer 1

大概写了一下不知道是否有帮助：

function redPnt
oriPic=imread('test1.png');
subplot(2,2,1)
imshow(oriPic)

% 删除红色外的部分并构造二值图
grayPic=rgb2gray(oriPic);
grayPic(oriPic(:,:,1)<255)=255;
grayPic(grayPic<250)=0;
subplot(2,2,2)
imshow(grayPic)

% 图像膨胀，使未连接边缘连接
SE=[0 1 0;1 1 1;0 1 0];
bwPic=imerode(grayPic,SE);

% 边缘清理：保留圆圈联通区域
bwPic=imclearborder(bwPic);
subplot(2,2,3)
imshow(bwPic)

% 获取每一个联通区域
LPic=bwlabel(bwPic);
labelNum=max(max(LPic));

% 计算每一个联通区域 坐标均值
pointSet=zeros(labelNum,2);
for i=1:labelNum
    [X,Y]=find(LPic==i);
    Xmean=mean(X);
    Ymean=mean(Y);
    pointSet(i,:)=[Xmean,Ymean];
end

subplot(2,2,4)
imshow(bwPic)
hold on
scatter(pointSet(:,2),pointSet(:,1),'r','LineWidth',1)

n=1;
while ~isempty(pointSet)
    circleSetInd=1;
    for j=1:length(pointSet)
        disSet=sqrt(sum((pointSet-pointSet(circleSetInd(end),:)).^2,2));
        [~,ind]=sort(disSet);
        ind=ind(1:5);
        [~,~,t_ind]=intersect(circleSetInd,ind);
        ind(t_ind)=[];
        if ~isempty(ind)
            circleSetInd=[circleSetInd;ind(1)];
        else
            circleSet{n}=pointSet(circleSetInd,:);
            pointSet(circleSetInd,:)=[];
            n=n+1;
            break
        end
    end
end

figure
imshow(oriPic)
hold on
for i=1:n-1
plot(circleSet{i}(:,2),circleSet{i}(:,1),'LineWidth',2)
end

end

没有用霍夫圆形检测之类算法，只是单纯的图像膨胀，边缘清理，寻找联通区域，计算联通区域坐标均值
大概过程如下图所示：