广告

扫码关注官方微信

扫码下载APP
返回顶部
首页 > 资讯 > 后端开发 > 其他教程 >详解基于Matlab的空心散点检测
  • 298

0
分享到

详解基于Matlab的空心散点检测

2024-04-02 19:04:59 298人浏览 独家记忆
摘要

目录问题描述过程展示图像导入依据RGB值图像二值化图像腐蚀图像边缘清理联通区域查找与坐标均值计算圈查找完整代码其它形状空心散点检测问题描述 有一张这样的图片，如何提取里面的红色圈圈坐

问题描述

有一张这样的图片，如何提取里面的红色圈圈坐标，并且连接这些坐标形成两个封闭的环路？

过程展示

图像导入

oriPic=imread('test1.png');

subplot(2,2,1)
imshow(oriPic)

依据RGB值图像二值化

原理就是图中颜色种类比较少，只有红黑白，而红色和白色都是R通道数值较大，因此我们可以利用这一点进行图像分割

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

%subplot(2,2,2)
figure
imshow(grayPic)

图像腐蚀

对于白色来说是腐蚀，对于黑色来说是膨胀，这一步是为了让那些有缺口的小圆圈将缺口补起来

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

figure
imshow(bwPic)

图像边缘清理

就是把和边缘连接的不被黑色包围的区域变成黑色：

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

联通区域查找与坐标均值计算

现在每一个白点都是一个坐标区域，我们检测所有联通区域并计算各个区域的重心即可：

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

% 计算每一个联通区域 坐标均值
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)
figure
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

这效果就很美滋滋：

完整代码

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

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

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

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

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

% 计算每一个联通区域 坐标均值
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)
figure
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

其它形状空心散点检测

来波正方形试试：

可以看出效果还是很棒的，当然大家可以根据实际情况自行更改图像腐蚀模板形状，如果散点是其它颜色请自行更改第一步的图像分割条件。

后注：

若是因为点较为密集而导致圈形路径内部白色区域没被清除，可能会将内部区域也算作散点造成错误，解决方法是计算每个联通区域面积并剔除远远大于区域面积中位数的联通区域：

问题出现原因的图片描述：

如图所示种间那一大片区域也被算作散点

更改后代码如下：

function redPnt
oriPic=imread('test2.png');
figure
imshow(oriPic)

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

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

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

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

% 筛掉超大区域
pointSizeSet=zeros(1,labelNum);
for i=1:labelNum
    pointSizeSet(i)=sum(sum(LPic==i));
end
[~,ind]=find(pointSizeSet>10*median(pointSizeSet));

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


figure
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:min(5,length(ind)));
        [~,~,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

注：

2016版本及以前可能这句：

disSet=sqrt(sum((pointSet-pointSet(circleSetInd(end),:)).^2,2));

会出现数组大小不匹配问题，可以将其改为：

tempMat=repmat(pointSet(circleSetInd(end),:),[size(pointSet,1),1]);
disSet=sqrt(sum((pointSet-tempMat).^2,2));

以上就是详解基于Matlab的空心散点检测的详细内容，更多关于Matlab空心散点检测的资料请关注编程网其它相关文章！

点击免费下载>>软考高级考试备考技巧/历年真题/备考精华资料

--结束END--

本文标题: 详解基于Matlab的空心散点检测

本文链接: https://www.lsjlt.com/news/140784.html(转载时请注明来源链接)

有问题或投稿请发送至: 邮箱/279061341@qq.com    QQ/279061341

本篇文章演示代码以及资料文档资料下载

下载Word文档到电脑，方便收藏和打印～

下载Word文档
去做题
猜你喜欢
  • 详解基于Matlab的空心散点检测
    目录问题描述过程展示图像导入依据RGB值图像二值化图像腐蚀图像边缘清理联通区域查找与坐标均值计算圈查找完整代码其它形状空心散点检测问题描述 有一张这样的图片，如何提取里面的红色圈圈坐...
    99+
    2022-11-13
  • Matlab空心散点检测的示例分析
    这篇文章主要介绍Matlab空心散点检测的示例分析，文中介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们一定要看完！问题描述有一张这样的图片，如何提取里面的红色圈圈坐标，并且连接这些坐标形成两个封闭的环路？过程展示图像导入oriPic=...
    99+
    2023-06-29
  • 基于Python的人脸检测与分类过程详解
    目录人脸识别算法简介人脸检测简述数据集介绍算法介绍测试网络结果预览人脸识别 算法简介 我们的算法可以分成两个部分，识别人脸位置和确定人脸分类。这两个部分可以看成：1.检测出人脸之间相...
    99+
    2022-11-11
软考高级职称资格查询
软考职称历年真题下载
软考职称资料下载
热门wiki
近期文章
推荐阅读
热门问答
热门标签
编程网，编程工程师的家园，是目前国内优秀的开源技术社区之一，形成了由开源软件库、代码分享、资讯、协作翻译、讨论区和博客等几大频道内容，为IT开发者提供了一个发现、使用、并交流开源技术的平台。

  • 官方手机版

  • 微信公众号

  • 商务合作