CN107369161A - 一种基于改进欧式聚类的散乱工件点云分割方法 - Google Patents
一种基于改进欧式聚类的散乱工件点云分割方法 Download PDFInfo
- Publication number
- CN107369161A CN107369161A CN201710591494.9A CN201710591494A CN107369161A CN 107369161 A CN107369161 A CN 107369161A CN 201710591494 A CN201710591494 A CN 201710591494A CN 107369161 A CN107369161 A CN 107369161A
- Authority
- CN
- China
- Prior art keywords
- point
- mrow
- msub
- point cloud
- cloud
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10028—Range image; Depth image; 3D point clouds
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30108—Industrial image inspection
- G06T2207/30164—Workpiece; Machine component
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Probability & Statistics with Applications (AREA)
- Image Analysis (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
本发明提供了一种基于改进欧式聚类的散乱工件点云分割方法,涉及点云分割领域,该方法考虑到散乱工件点云固有的杂乱性和无序性,提出相应的场景分割方案,具体步骤为:对点云进行预处理,包括使用RANSAC法去除背景点、迭代半径滤波法去除离群点;采用线下模板点云的信息注册方法为线上分割提供参数选取依据,从而提高了线上分割的速度;提出先去除边缘点、再聚类分割、最后补齐边缘点的思想,避免了聚类过程中出现欠分割或过分割的现象,在聚类分割时,提出基于自适应邻域搜索半径的聚类方法,大大提高了分割速度,边缘点补齐则保留了工件的表面特征,有利于提高后续位姿定位准确性。
Description
技术领域
本发明涉及点云分割领域,尤其涉及一种基于改进欧式聚类的散乱工件点云分割方法。
背景技术
近年来随着三维扫描设备精度的提高和成本的降低,研究者能够快速精确地获取对象表面的三维点云信息。相对于二维图像,点云包含了对象的深度信息,在目标识别和定位方面具有独特的优势和潜力,因此该项技术在随机箱体抓取(Random Bin Picking,RBP)领域引起了广泛的关注。使用三维扫描设备获取箱体内散乱工件表面点云,结合点云处理算法计算单个工件位姿,引导工业机器人进行抓取,效率更高、速度更快、智能化程度更好。但是,受箱体内工件杂乱性,物体间的重叠、碰撞等因素的影响,如何将单个工件从复杂场景中准确的分离出来,仍是有效抓取面临的主要问题,分割结果的好坏,会直接影响后续位姿定位的准确性。
常见的分割算法有基于区域生长的分割、基于模型的分割、基于属性的分割。基于区域生长的点云分割因为效果好,被广泛应用在三维场景重建,逆向工程等领域,但是时间消耗大,不适合解决RBP问题;基于模型的分割方法采用几何图元进行归类,将具有相同数学表示的点归为一类,常用的方法为RANSAC(随机采样一致性)算法。该方法使用纯粹的数学方法,速度快、对于离群点的稳健性好,但是在处理不同来源点云时精度不高;基于属性的分割是将具有相同属性的点归为一类,也称聚类分割,分割结果灵活准确,但是在处理大场景多维点云时的时间消耗大。
本发明采用编码结构光3D视觉测量***获取点云,综合考虑该***获取的点云的特点、分割效果和实时性,基于属性的分割方法效果最好。聚类分割常用的属性如边缘、欧式距离、表面法线等,其中基于欧式距离的聚类分割速度最快,满足工业生产的实时性要求。欧式聚类算法根据点与点之间的欧氏距离进行归类,将距离小于阈值的点归为当前类,但该方法需要手动设置点与点之间的距离阈值,阈值过大过小会导致不同程度的欠分割和过分割,对分割的效果造成影响。
发明内容
本发明为实现散乱工件点云的分割,提供了一种基于改进欧式聚类的散乱工件点云分割方法,该方法在满足实时性要求的同时,能很好的分割出单个目标工件,且最大程度上保留了工件表面的特征。
为了达到此目的,本发明提供的技术方案实现如下:
一种基于改进欧式聚类的散乱工件点云分割方法,步骤如下:
(1)对目标点云和模板点云进行相同预处理,通过迭代半径滤波方法去除聚集分布的离群点;
(2)对步骤(1)中预处理后的模板点云进行线下信息注册,获取单个工件点云的信息,所述信息包括所述单个工件点云在不同k邻域下点与邻近点之间的平均距离和标准差以及所述单个工件点云在不同r邻域下点与邻近点之间的法线夹角均值;
(3)将步骤(2)中所述单个工件点云在不同r邻域下点与邻近点之间的法线夹角均值设置为阀值,将局部邻域法线夹角均值小于阈值的点判定为边缘点,从而分离出目标点云的边缘点和非边缘点,并分别保存成两幅点云;
(4)对步骤(3)得到的所述非边缘点点云进行分割,操作如下:将步骤(2)中所述k邻域下点与邻近点之间的平均距离和标准差设定初始邻域搜索半径阈值,将目标点云中点和其k个邻近点的平均距离和所设阈值相比较,根据两者差异确定自适应搜索半径,然后进行聚类,最终将散乱工件点云分割成多个包含单个工件的点云子集。
(5)对步骤(4)中所述多个包含单个工件的点云子集筛选聚类,滤除数目过少的子集,得到新的点云子集;
(6)对步骤(3)中得到的所述边缘点点云集中搜索属于它的边缘点并归入步骤(5)所述新的点云子集进行边缘点补齐,最终完成散乱工件点云分割。
进一步地,所述步骤(1)中通过迭代半径滤波方法去除聚集分布的离群点,其过程如下:
第一步,采用kd树建立无序点云点与点之间的拓扑关系,实现单个点邻近点的快速查找;
第二步,对点云中任意一点pi,搜索其r邻域半径内的邻近点集合将k值大小和预先设定的点云密度阈值nth比较,若k<nth,则点pi判定为离群点,并从点云中移除;对于少数聚集分布的离群点,在不改变邻域半径大小的情况下进行多次重复滤波,直到滤波前后点云的数目相等终止滤波。
进一步地,所述步骤(2)中进行模板点云进行线下信息注册,其过程如下:
第一步,设模板点云有n个数据点,构建模板点云kd树,对任意一点qi∈Q,搜索其r邻域半径内的所有邻近点,得到邻近点集合将集合内点个数记作单个点密度,遍历点云中的所有点,取平均值,得到模板点云的平均密度;
第二步,对任意一点qi∈Q,搜索距离其最近的k个邻近点,记作{qi1,qi2,…qik},计算每个邻近点距qi点的平均距离如下式所示;
其中dij表示第j个邻近点到qi的距离,遍历点云Q中所有的点,得到n个平均距离值,根据下式求得点云Q在k邻域下的平均距离均值和标准差;
第三步,采用PCA(主成分分析法)求取点云中所有点的法线,具体步骤是:对任意一点qi∈Q,构建其k邻近点{qi1,qi2,…qik}对应的协方差矩阵M:
其中表示最邻近元素的三维质心,λi和vi分别代表协方差矩阵的特征值和特征向量,λ0≤λ1≤λ2。特征向量所在的平面互相正交;具有最小特征值的特征向量即qi点的法线,记作点qi的k邻近点对应法线分别为根据如下公式:
计算得到qi点与k邻近点法线夹角的均值、点云Q的k邻域法线夹角均值。
本发明的有益效果:提供了一种基于改进欧式聚类的散乱工件点云分割方法,该方法能够快速高效地将箱体内散乱堆放工件点云分割成多个包含单个工件的点云子集,同时尽可能地保留工件表面的特征,为后续工件的定位和抓取奠定良好的基础。为了去除点云中的离群点,在基于邻域半径内点云密度的基础上,引入迭代思想,提出迭代半径滤波方法,去除少数聚集分布的离群点,提高了离群点去除的效率;线下处理阶段,对模板点云进行信息注册,为在线分割提供参数选取依据,提高了分割结果的准确性;在线点云分割阶段,通过边缘点去除和补齐避免欧式聚类容易出现欠分割或过分割的现象,保证了分割的效率,自适应搜索半径的欧式聚类方法,大大提高了分割速度。
附图说明
图1是本发明的整体流程图
图2是搜索半径自适应的聚类分割算法流程图
具体实施方式
为了更清楚的说明本发明的技术方案和优点,下面结合具体实施例,并参附图,阐述本发明的具体实施方式。
本发明的目的是将箱体内散乱堆放工件点云分割成多个包含单个工件的点云子集,主要流程分成以下五个部分:点云预处理、模板点云线下信息注册、目标点云边缘点提取、基于自适应邻域搜索半径的聚类分割、边缘点补齐,如图1所示。
具体实现步骤为:
(1)点云预处理(以目标点云P为例)
(1.1)采用随机采样一致性算法计算箱体底部平面方程,将点云中接近箱底平面的点去除。
(1.2)迭代半径滤波方法去除离群点的思想是:首先采用kd树邻近点搜索算法建立无序点云P点与点之间的拓扑关系,实现邻近点的快速查找;然后,对任意一点pi∈P,搜索其r邻域半径内的邻近点集合将k值大小和预先设定的点云密度阈值nth比较。若k<nth,则点pi判定为离群点,并从点云P中移除,在不改变邻域半径大小的情况下进行多次重复滤波,直到滤波前后点云的数目相等终止滤波。
(1.3)对去除离群点后的点云使用栅格法进行降采样,根据点云在三个维度上最小包围盒和点云分辨率定义单个体素网格大小,用距离重心最近的点代替体素网格内所有的点,达到点云精简的目的。
(2)线下模板点云Q的信息注册
(2.1)设模板点云有n个数据点,构建模板点云kd树,对任意一点qi∈Q,搜索其r邻域半径内的所有邻近点,得到邻近点集合将集合内点个数记作单个点密度,遍历点云中的所有点,取平均值,得到模板点云的平均密度。
(2.2)对任意一点qi∈Q,搜索距离其最近的k个邻近点,记作{qi1,qi2,…qik},计算每个邻近点距qi点的平均距离如下式所示。
其中dij表示第j个邻近点到qi的距离,遍历点云Q中所有的点,得到n个平均距离值,根据下式求得点云Q在k邻域下的平均距离均值和标准差。
(2.3)采用PCA(主成分分析法)求取点云中所有点的法线,具体步骤是:
对任意一点qi∈Q,构建其k邻近点{qi1,qi2,…qik}对应的协方差矩阵M:
其中表示最邻近元素的三维质心,λi和vi分别代表协方差矩阵的特征值和特征向量,λ0≤λ1≤λ2。特征向量所在的平面互相正交。具有最小特征值的特征向量即qi点的法线,记作点qi的k邻近点对应法线分别为根据如下公式计算qi点与k邻近点法线夹角的均值、点云Q的k邻域法线夹角均值。
(3)目标点云边缘点提取
根据(2.3)中得到的法线夹角均值设置阈值,对任意一点pi∈P,计算该点与k个邻近点法线夹角的均值将该值与所设阈值相比较,若大于阈值,则pi点为边缘点,反之则为非边缘点。遍历P中所有的元素,最终提取所有的边缘点记作点云P1,非边缘点记作点云P2,P1,P2∈P且P1∪P2=P。
(4)基于自适应邻域搜索半径的聚类分割
欧式聚类分割是根据点与点之间的欧式距离是否满足所设阈值进行归类的算法,其缺点在于邻域搜索半径是固定的,半径较大或较小容易出现欠分割和过分割的现象。而基于自适应半径阈值的聚类分割方法可以根据点云局部邻域密度自适应确定搜索半径,避免了欠分割或过分割现象的产生,同时提高了算法的速度,其流程如图2所示。具体实现步骤如下:
1)构建非边缘点云P2的kd树,以(2.1)节中求得的点云密度为参考初始化邻域点个数k,由(2.2)求得k邻近点距离均值mean和标准差Devia,初始化最小聚类数目Min和最大聚类数目Max。
2)新建一个空的点云索引向量cluster_indices用于存储聚类结果;新建一个整型向量queue并初始化为空,用于存储单次聚类结果;新建一个bool型向量processed,尺寸和点云P2相同,初始化向量为false,用于表示点云是否被处理过。
3)初始化i=0,对任意一点pi∈P2,执行以下步骤:
a)初始化c=0,将pi点的索引值加入当前序列queue[c],标记该点为已分类,即processed[i]=true;
b)根据序列queue中第c个元素存储的点索引找到其在点云P2中的对应点pc,在kd树中,搜索距离该点最近的(k+1)个点,记作集合求k个邻近点与pi点的距离{Dist[j],j=1,2...k}和均值AverSinglePoint,根据如下公式判断pi点是否为非边缘点,如果满足,则自适应搜索半径Rth根据如下公式求得。
AverSinglePoint≤(mean+μ*Devia)
Rth=Max{Dist[j],j=1,2...k}
反之,自适应搜索半径Rth根据如下公式求得。
c)根据b)求得的搜索半径Rth,将位于该半径范围内邻近点的索引加入序列queue,并标记邻近点为已分类,c=c+1;
d)重复步骤b)、c),直到queue中对应的所有点都被处理过,根据如下公式判断该聚类结果是否满足数目要求,若满足,则将queue中所有的点索引放入点云索引向量cluster_indices中,清空队列queue,反之,则舍弃并清空序列queue,i=i+1。
Min≤queue.size()≤Max
4)重复步骤3),直到点云P2中所有的点都已处理完毕,聚类结果存放在cluster_indices[0],cluster_indices[1]...cluster_indices[num-1]中,将每组点云索引向量中索引在P2中对应的点保存成点云文件输出。
(5)边缘点补齐
经过步骤(4),目标点云P2被分割成num个点云子集,每个子集都代表一个目标工件,表示为{cluster1,cluster2...clusterM},但此时所有子集所表示的工件缺乏边缘部分,边缘的缺失会增大位姿估计的误差,因此需要进行边缘点的修补。根据边缘点与其在点云子集中最邻近点的距离是否满足阈值进行边缘点补齐,过程如下:
1)初始化j=0,初始化阈值THk。
2)构建子集clusterj的kd树,对边缘点云P1中任意一点P1i,计算该点到子集clusterj中最邻近点的欧式距离dj,比较dj和THk大小,若dj<THk,则将P1点归入子集clusterj,并标记该点已处理,反之取下一个边缘点。
3)遍历点云P1中的数据点,检测所有属于子集clusterj的边缘点,完成边缘点补齐。
4)重复步骤2)3),直到所有子集边缘点补齐完全,输出补齐边缘点后的点云子集,完成点云分割。
以上所述仅为说明本发明的实施方式,并不用于限制本发明,对于本领域的技术人员来说,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种基于改进欧式聚类的散乱工件点云分割方法,其特征在于,步骤如下:
(1)对目标点云和模板点云进行相同预处理,通过迭代半径滤波方法去除聚集分布的离群点;
(2)对步骤(1)中预处理后的模板点云进行线下信息注册,获取单个工件点云的信息,所述信息包括所述单个工件点云在不同k邻域下点与邻近点之间的平均距离和标准差以及所述单个工件点云在不同r邻域下点与邻近点之间的法线夹角均值;
(3)将步骤(2)中所述单个工件点云在不同r邻域下点与邻近点之间的法线夹角均值设置为阀值,将局部邻域法线夹角均值小于阈值的点判定为边缘点,从而分离出目标点云的边缘点和非边缘点,并分别保存成两幅点云;
(4)对步骤(3)得到的所述非边缘点点云进行分割,操作如下:将步骤(2)中所述k邻域下点与邻近点之间的平均距离和标准差设定初始邻域搜索半径阈值,将目标点云中点和其k个邻近点的平均距离和所设阈值相比较,根据两者差异确定自适应搜索半径,然后进行聚类,最终将散乱工件点云分割成多个包含单个工件的点云子集。
(5)对步骤(4)中所述多个包含单个工件的点云子集筛选聚类,滤除数目过少的子集,得到新的点云子集;
(6)对步骤(3)中得到的所述边缘点点云集中搜索属于它的边缘点并归入步骤(5)所述新的点云子集进行边缘点补齐,最终完成散乱工件点云分割。
2.根据权利要求1所述的一种基于改进欧式聚类的散乱工件点云分割方法,其特征在于,所述步骤(1)中通过迭代半径滤波方法去除聚集分布的离群点,其过程如下:
第一步,采用kd树建立无序点云点与点之间的拓扑关系,实现单个点邻近点的快速查找;
第二步,对点云中任意一点pi,搜索其r邻域半径内的邻近点集合将k值大小和预先设定的点云密度阈值nth比较,若k<nth,则点pi判定为离群点,并从点云中移除;对于少数聚集分布的离群点,在不改变邻域半径大小的情况下进行多次重复滤波,直到滤波前后点云的数目相等终止滤波。
3.根据权利要求1所述的一种基于改进欧式聚类的散乱工件点云分割方法,其特征在于,所述步骤(2)中进行模板点云进行线下信息注册,其过程如下:
第一步,设模板点云有n个数据点,构建模板点云kd树,对任意一点qi∈Q,搜索其r邻域半径内的所有邻近点,得到邻近点集合将集合内点个数记作单个点密度,遍历点云中的所有点,取平均值,得到模板点云的平均密度;
第二步,对任意一点qi∈Q,搜索距离其最近的k个邻近点,记作{qi1,qi2,…qik},计算每个邻近点距qi点的平均距离如下式所示;
<mrow>
<msub>
<mover>
<mi>d</mi>
<mo>&OverBar;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>k</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<msub>
<mi>d</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
其中dij表示第j个邻近点到qi的距离,遍历点云Q中所有的点,得到n个平均距离值,根据下式求得点云Q在k邻域下的平均距离均值和标准差;
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>d</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&sigma;</mi>
<mo>=</mo>
<msqrt>
<mrow>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mover>
<mi>d</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
第三步,采用PCA(主成分分析法)求取点云中所有点的法线,具体步骤是:对任意一点qi∈Q,构建其k邻近点{qi1,qi2,…qik}对应的协方差矩阵M:
<mrow>
<mi>M</mi>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>k</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<mover>
<mi>q</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<mover>
<mi>q</mi>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
<mrow>
<mi>M</mi>
<mo>&CenterDot;</mo>
<mover>
<msub>
<mi>v</mi>
<mi>j</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
<mo>=</mo>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mo>&CenterDot;</mo>
<mover>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
<mo>,</mo>
<mi>j</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>}</mo>
</mrow>
其中表示最邻近元素的三维质心,λi和vi分别代表协方差矩阵的特征值和特征向量,λ0≤λ1≤λ2。特征向量所在的平面互相正交;具有最小特征值的特征向量即qi点的法线,记作点qi的k邻近点对应法线分别为根据如下公式:
<mrow>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>k</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
<mi>k</mi>
</munderover>
<msup>
<mi>cos</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<mover>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
<mo>&CenterDot;</mo>
<mover>
<msub>
<mi>v</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&RightArrow;</mo>
</mover>
</mrow>
<mrow>
<mo>|</mo>
<mover>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mo>&RightArrow;</mo>
</mover>
<mo>|</mo>
<mo>|</mo>
<mover>
<msub>
<mi>v</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&RightArrow;</mo>
</mover>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
<mrow>
<mo>&part;</mo>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>n</mi>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<msub>
<mi>&alpha;</mi>
<mi>i</mi>
</msub>
</mrow>
计算得到qi点与k邻近点法线夹角的均值、点云Q的k邻域法线夹角均值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710591494.9A CN107369161B (zh) | 2017-07-19 | 2017-07-19 | 一种基于改进欧式聚类的散乱工件点云分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710591494.9A CN107369161B (zh) | 2017-07-19 | 2017-07-19 | 一种基于改进欧式聚类的散乱工件点云分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107369161A true CN107369161A (zh) | 2017-11-21 |
CN107369161B CN107369161B (zh) | 2020-09-11 |
Family
ID=60308327
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710591494.9A Active CN107369161B (zh) | 2017-07-19 | 2017-07-19 | 一种基于改进欧式聚类的散乱工件点云分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107369161B (zh) |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108171780A (zh) * | 2017-12-28 | 2018-06-15 | 电子科技大学 | 一种基于激光雷达构建室内真实三维地图的方法 |
CN108830902A (zh) * | 2018-04-19 | 2018-11-16 | 江南大学 | 一种基于点云处理的散乱工件识别与定位方法 |
CN108872991A (zh) * | 2018-05-04 | 2018-11-23 | 上海西井信息科技有限公司 | 目标物检测与识别方法、装置、电子设备、存储介质 |
CN109559346A (zh) * | 2018-11-07 | 2019-04-02 | 西安电子科技大学 | 一种3d点云测量中的待测部位的定位及分割方法、扫描仪 |
CN110136072A (zh) * | 2018-02-08 | 2019-08-16 | 北京大学 | 点云噪声的去除方法、去噪***、计算机设备及存储介质 |
CN110188810A (zh) * | 2019-05-23 | 2019-08-30 | 西北工业大学 | 一种基于密度聚类的在线野值检测、识别与修正方法 |
CN110599506A (zh) * | 2019-10-16 | 2019-12-20 | 湖南大学 | 一种复杂异形曲面机器人三维测量的点云分割方法 |
CN110853044A (zh) * | 2019-04-25 | 2020-02-28 | 华中农业大学 | 基于条件欧式聚类的盆栽玉米点云快速分割方法 |
CN111152229A (zh) * | 2020-02-24 | 2020-05-15 | 深圳群宾精密工业有限公司 | 3d机械视觉的机械手引导方法和装置 |
CN111210500A (zh) * | 2018-11-22 | 2020-05-29 | 北京欣奕华科技有限公司 | 一种三维点云处理方法及装置 |
CN111275810A (zh) * | 2020-01-17 | 2020-06-12 | 五邑大学 | 基于图像处理的k近邻点云滤波方法、装置和存储介质 |
CN111553938A (zh) * | 2020-04-29 | 2020-08-18 | 南京航空航天大学 | 一种基于图优化的多站位扫描点云全局配准方法 |
CN111815706A (zh) * | 2020-06-23 | 2020-10-23 | 熵智科技(深圳)有限公司 | 面向单品类拆垛的视觉识别方法、装置、设备及介质 |
CN111932570A (zh) * | 2020-09-10 | 2020-11-13 | 熵智科技(深圳)有限公司 | 一种基于网格数据的边缘检测方法、装置、介质及设备 |
CN111932688A (zh) * | 2020-09-10 | 2020-11-13 | 深圳大学 | 一种基于三维点云的室内平面要素提取方法、***及设备 |
CN113012161A (zh) * | 2021-03-05 | 2021-06-22 | 华南理工大学 | 一种基于凸性区域生长的堆叠散乱目标点云分割方法 |
CN113139721A (zh) * | 2021-04-16 | 2021-07-20 | 深圳市艾赛克科技有限公司 | 骨料堆场管理***及方法 |
CN113313137A (zh) * | 2021-04-01 | 2021-08-27 | 杭州兰亭视觉科技有限责任公司 | 涂胶障碍物识别及去除方法 |
CN113344992A (zh) * | 2021-05-31 | 2021-09-03 | 山东大学 | 一种全局点云配准方法、***、存储介质及设备 |
CN114462493A (zh) * | 2021-12-29 | 2022-05-10 | 浙江大华技术股份有限公司 | 聚类方法、聚类装置及计算机可读存储介质 |
CN114492619A (zh) * | 2022-01-22 | 2022-05-13 | 电子科技大学 | 一种基于统计和凹凸性的点云数据集构建方法及装置 |
CN114612598A (zh) * | 2022-02-16 | 2022-06-10 | 苏州一径科技有限公司 | 一种点云的处理方法、装置及激光雷达 |
CN114754741A (zh) * | 2022-06-15 | 2022-07-15 | 湖南大学 | 边坡断面坡度测量方法 |
CN115760855A (zh) * | 2023-01-09 | 2023-03-07 | 中建科技集团有限公司 | 工件检查方法及相关设备 |
CN116612132A (zh) * | 2023-05-25 | 2023-08-18 | 重庆茂侨科技有限公司 | 基于集料特征的3d点云目标分割方法 |
CN117455928A (zh) * | 2023-12-25 | 2024-01-26 | 珠海市格努科技有限公司 | 无序抓取过程中抓取对象的分割方法、装置和电子设备 |
CN117808703A (zh) * | 2024-02-29 | 2024-04-02 | 南京航空航天大学 | 一种多尺度大型部件装配间隙点云滤波方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104143194A (zh) * | 2014-08-20 | 2014-11-12 | 清华大学 | 一种点云分割方法及装置 |
CN104700398A (zh) * | 2014-12-31 | 2015-06-10 | 西安理工大学 | 一种点云场景物体提取方法 |
CN105740798A (zh) * | 2016-01-27 | 2016-07-06 | 西安理工大学 | 一种基于结构分析的点云场景物体识别方法 |
CN105844629A (zh) * | 2016-03-21 | 2016-08-10 | 河南理工大学 | 一种大场景城市建筑物立面点云自动分割方法 |
CN106600622A (zh) * | 2016-12-06 | 2017-04-26 | 西安电子科技大学 | 一种基于超体素的点云数据分割方法 |
-
2017
- 2017-07-19 CN CN201710591494.9A patent/CN107369161B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104143194A (zh) * | 2014-08-20 | 2014-11-12 | 清华大学 | 一种点云分割方法及装置 |
CN104700398A (zh) * | 2014-12-31 | 2015-06-10 | 西安理工大学 | 一种点云场景物体提取方法 |
CN105740798A (zh) * | 2016-01-27 | 2016-07-06 | 西安理工大学 | 一种基于结构分析的点云场景物体识别方法 |
CN105844629A (zh) * | 2016-03-21 | 2016-08-10 | 河南理工大学 | 一种大场景城市建筑物立面点云自动分割方法 |
CN106600622A (zh) * | 2016-12-06 | 2017-04-26 | 西安电子科技大学 | 一种基于超体素的点云数据分割方法 |
Non-Patent Citations (1)
Title |
---|
LANLAN LI,ETC: "Point Cloud Simplification Based on an Affinity Propagation Clustering Algorithm", <2009 INTERNATIONAL CONFERENCE ON ARTIFICIAL INTELLIGENCE AND COMPUTATIONAL INTELLIGENCE> * |
Cited By (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108171780A (zh) * | 2017-12-28 | 2018-06-15 | 电子科技大学 | 一种基于激光雷达构建室内真实三维地图的方法 |
CN110136072B (zh) * | 2018-02-08 | 2021-02-09 | 北京大学 | 点云噪声的去除方法、去噪***、计算机设备及存储介质 |
CN110136072A (zh) * | 2018-02-08 | 2019-08-16 | 北京大学 | 点云噪声的去除方法、去噪***、计算机设备及存储介质 |
CN108830902A (zh) * | 2018-04-19 | 2018-11-16 | 江南大学 | 一种基于点云处理的散乱工件识别与定位方法 |
CN108872991A (zh) * | 2018-05-04 | 2018-11-23 | 上海西井信息科技有限公司 | 目标物检测与识别方法、装置、电子设备、存储介质 |
CN109559346A (zh) * | 2018-11-07 | 2019-04-02 | 西安电子科技大学 | 一种3d点云测量中的待测部位的定位及分割方法、扫描仪 |
CN109559346B (zh) * | 2018-11-07 | 2021-12-14 | 西安电子科技大学 | 一种3d点云测量中的待测部位的定位及分割方法、扫描仪 |
CN111210500A (zh) * | 2018-11-22 | 2020-05-29 | 北京欣奕华科技有限公司 | 一种三维点云处理方法及装置 |
CN111210500B (zh) * | 2018-11-22 | 2023-08-29 | 浙江欣奕华智能科技有限公司 | 一种三维点云处理方法及装置 |
CN110853044A (zh) * | 2019-04-25 | 2020-02-28 | 华中农业大学 | 基于条件欧式聚类的盆栽玉米点云快速分割方法 |
CN110188810A (zh) * | 2019-05-23 | 2019-08-30 | 西北工业大学 | 一种基于密度聚类的在线野值检测、识别与修正方法 |
CN110599506B (zh) * | 2019-10-16 | 2023-03-24 | 湖南大学 | 一种复杂异形曲面机器人三维测量的点云分割方法 |
CN110599506A (zh) * | 2019-10-16 | 2019-12-20 | 湖南大学 | 一种复杂异形曲面机器人三维测量的点云分割方法 |
CN111275810A (zh) * | 2020-01-17 | 2020-06-12 | 五邑大学 | 基于图像处理的k近邻点云滤波方法、装置和存储介质 |
CN111275810B (zh) * | 2020-01-17 | 2022-06-24 | 五邑大学 | 基于图像处理的k近邻点云滤波方法、装置和存储介质 |
CN111152229B (zh) * | 2020-02-24 | 2023-08-25 | 群滨智造科技(苏州)有限公司 | 3d机械视觉的机械手引导方法和装置 |
CN111152229A (zh) * | 2020-02-24 | 2020-05-15 | 深圳群宾精密工业有限公司 | 3d机械视觉的机械手引导方法和装置 |
CN111553938A (zh) * | 2020-04-29 | 2020-08-18 | 南京航空航天大学 | 一种基于图优化的多站位扫描点云全局配准方法 |
CN111815706B (zh) * | 2020-06-23 | 2023-10-27 | 熵智科技(深圳)有限公司 | 面向单品类拆垛的视觉识别方法、装置、设备及介质 |
CN111815706A (zh) * | 2020-06-23 | 2020-10-23 | 熵智科技(深圳)有限公司 | 面向单品类拆垛的视觉识别方法、装置、设备及介质 |
CN111932688A (zh) * | 2020-09-10 | 2020-11-13 | 深圳大学 | 一种基于三维点云的室内平面要素提取方法、***及设备 |
CN111932570B (zh) * | 2020-09-10 | 2021-01-19 | 熵智科技(深圳)有限公司 | 一种基于网格数据的边缘检测方法、装置、介质及设备 |
CN111932570A (zh) * | 2020-09-10 | 2020-11-13 | 熵智科技(深圳)有限公司 | 一种基于网格数据的边缘检测方法、装置、介质及设备 |
CN113012161B (zh) * | 2021-03-05 | 2023-04-28 | 华南理工大学 | 一种基于凸性区域生长的堆叠散乱目标点云分割方法 |
CN113012161A (zh) * | 2021-03-05 | 2021-06-22 | 华南理工大学 | 一种基于凸性区域生长的堆叠散乱目标点云分割方法 |
CN113313137A (zh) * | 2021-04-01 | 2021-08-27 | 杭州兰亭视觉科技有限责任公司 | 涂胶障碍物识别及去除方法 |
CN113139721B (zh) * | 2021-04-16 | 2023-12-19 | 深圳市艾赛克科技有限公司 | 骨料堆场管理***及方法 |
CN113139721A (zh) * | 2021-04-16 | 2021-07-20 | 深圳市艾赛克科技有限公司 | 骨料堆场管理***及方法 |
CN113344992B (zh) * | 2021-05-31 | 2022-06-28 | 山东大学 | 一种全局点云配准方法、***、存储介质及设备 |
CN113344992A (zh) * | 2021-05-31 | 2021-09-03 | 山东大学 | 一种全局点云配准方法、***、存储介质及设备 |
CN114462493B (zh) * | 2021-12-29 | 2024-05-10 | 浙江大华技术股份有限公司 | 聚类方法、聚类装置及计算机可读存储介质 |
CN114462493A (zh) * | 2021-12-29 | 2022-05-10 | 浙江大华技术股份有限公司 | 聚类方法、聚类装置及计算机可读存储介质 |
CN114492619B (zh) * | 2022-01-22 | 2023-08-01 | 电子科技大学 | 一种基于统计和凹凸性的点云数据集构建方法及装置 |
CN114492619A (zh) * | 2022-01-22 | 2022-05-13 | 电子科技大学 | 一种基于统计和凹凸性的点云数据集构建方法及装置 |
CN114612598A (zh) * | 2022-02-16 | 2022-06-10 | 苏州一径科技有限公司 | 一种点云的处理方法、装置及激光雷达 |
CN114754741A (zh) * | 2022-06-15 | 2022-07-15 | 湖南大学 | 边坡断面坡度测量方法 |
CN115760855B (zh) * | 2023-01-09 | 2023-05-23 | 中建科技集团有限公司 | 工件检查方法及相关设备 |
CN115760855A (zh) * | 2023-01-09 | 2023-03-07 | 中建科技集团有限公司 | 工件检查方法及相关设备 |
CN116612132A (zh) * | 2023-05-25 | 2023-08-18 | 重庆茂侨科技有限公司 | 基于集料特征的3d点云目标分割方法 |
CN117455928A (zh) * | 2023-12-25 | 2024-01-26 | 珠海市格努科技有限公司 | 无序抓取过程中抓取对象的分割方法、装置和电子设备 |
CN117455928B (zh) * | 2023-12-25 | 2024-04-02 | 珠海市格努科技有限公司 | 无序抓取过程中抓取对象的分割方法、装置和电子设备 |
CN117808703A (zh) * | 2024-02-29 | 2024-04-02 | 南京航空航天大学 | 一种多尺度大型部件装配间隙点云滤波方法 |
CN117808703B (zh) * | 2024-02-29 | 2024-05-10 | 南京航空航天大学 | 一种多尺度大型部件装配间隙点云滤波方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107369161B (zh) | 2020-09-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107369161A (zh) | 一种基于改进欧式聚类的散乱工件点云分割方法 | |
CN111299815B (zh) | 一种用于低灰度橡胶垫视觉检测与激光切割轨迹规划方法 | |
Su et al. | A deeper look at 3D shape classifiers | |
CN109887015B (zh) | 一种基于局部曲面特征直方图的点云自动配准方法 | |
Lu et al. | Fast 3D line segment detection from unorganized point cloud | |
Zhong | Intrinsic shape signatures: A shape descriptor for 3D object recognition | |
CN106447676B (zh) | 一种基于快速密度聚类算法的图像分割方法 | |
CN108830902A (zh) | 一种基于点云处理的散乱工件识别与定位方法 | |
US7995055B1 (en) | Classifying objects in a scene | |
CN114170279A (zh) | 一种基于激光扫描的点云配准方法 | |
CN106780524A (zh) | 一种三维点云道路边界自动提取方法 | |
CN110363861A (zh) | 基于激光雷达点云的田地作物三维重构方法 | |
CN108510516A (zh) | 一种散乱点云的三维线段提取方法及*** | |
CN113628263A (zh) | 一种基于局部曲率及其近邻特征的点云配准方法 | |
CN106123812A (zh) | 基于遥感影像获取起伏地表甘蔗种植面积的方法及装置 | |
CN111145129A (zh) | 基于超体素的点云去噪方法 | |
CN115187803A (zh) | 一种用于名优茶嫩芽采摘过程的定位方法 | |
CN108564116A (zh) | 一种摄像头场景图像的成分智能分析方法 | |
CN107123138A (zh) | 基于vanilla‑R点对剔除策略的点云配准算法 | |
CN109035311A (zh) | 一种弯骨骨折自动配准及内固定钢板预弯建模方法 | |
Lalonde et al. | Data structures for efficient dynamic processing in 3-d | |
Zhang et al. | Object-based classification of urban airborne LiDAR point clouds with multiple echoes using SVM | |
CN106980878A (zh) | 三维模型几何风格的确定方法及装置 | |
CN106611418A (zh) | 一种图像分割算法 | |
Sun et al. | A study on the classification of vegetation point cloud based on random forest in the straw checkerboard barriers area |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |