CN103268631B - 点云骨架提取方法及装置 - Google Patents
点云骨架提取方法及装置 Download PDFInfo
- Publication number
- CN103268631B CN103268631B CN201310196243.2A CN201310196243A CN103268631B CN 103268631 B CN103268631 B CN 103268631B CN 201310196243 A CN201310196243 A CN 201310196243A CN 103268631 B CN103268631 B CN 103268631B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- point cloud
- munder
- point
- 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.)
- Active
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 22
- 238000009826 distribution Methods 0.000 claims abstract description 24
- 238000005070 sampling Methods 0.000 claims description 67
- 238000000034 method Methods 0.000 claims description 26
- 239000011159 matrix material Substances 0.000 claims description 19
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 238000009499 grossing Methods 0.000 claims description 8
- 230000008602 contraction Effects 0.000 abstract description 10
- 230000008569 process Effects 0.000 description 11
- 238000005516 engineering process Methods 0.000 description 9
- 235000011449 Rosa Nutrition 0.000 description 5
- 230000001965 increasing effect Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000004590 computer program Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011960 computer-aided design Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/56—Particle system, point based geometry or rendering
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
- Image Generation (AREA)
- Image Processing (AREA)
Abstract
一种点云骨架提取方法,包括:获取输入的点云采样数据;根据迭代公式进行点云收缩得到骨架分支,所述迭代公式为:其中,其中,i,j为正整数,J为点云采样数据的点集合,qj为点集合J中的采样点,I为采样点qj的邻域点集合,xi为邻域点集合I中的邻域点,为邻域点集合I中不属于xi的邻域点,R为规整项,γi为权重系数,h为邻域点集合I的邻域半径,σi为分布系数;连接所述骨架分支得到点云骨架。此外,还提供了一种点云骨架提取装置。上述点云骨架提取方法及装置能够提高点云骨架的准确性。
Description
技术领域
本发明涉及计算机图形处理技术领域,特别是涉及一种点云骨架提取方法及装置。
背景技术
当今的生产应用中,计算机辅助设计,逆向工程,虚拟现实,三维动画和游戏等三维技术发展的一个主要瓶颈是:至今仍没有一个便利的方法快速获得用计算机存储的三维模型。如何直接从扫描点云数据出发,直接快速得到实用的点云模型,仍是一个未完近年来,三维激光扫描仪凭借其能够方便灵活地获取现实物体的三维表面数据的优势,得到了广泛应用。然而,目前在世界范围全解决的难题。点云数据处理的难点在于,点云数据一般是散乱的,有大量缺失,噪声和外点的。经过二十多年的研究,从较为完整的点云数据中重构三维模型的技术已经较为成熟。然而,制约其成为一种通用的逆向重构技术的关键在于,点云数据往往有大面积的缺失,这是现有硬件扫描设备无法解决的问题。而现有的基于拟合的重构手段都难以快速从缺失的点云数据直接重构出令人满意的三维模型。
所谓点云模型,一般指由三维扫描设备发射扫描光到被测物体表面上,接收反射光,并以此推算出物体表面的三维坐标点的集合。所谓点云三维重建,是指在一定的点云模型数据基础上,重建出能表示源模型的网格数据,便于计算机渲染和用户交互。
为了更好地处理点云数据,从点云中提取骨架是很重要的一个预处理步骤。因为骨架中包含了物体的拓扑结构信息。所以点云骨架的提取其实相当于对原物体的形状的理解过程。一旦得到了点云的骨架后,我们可以利用骨架来反补缺失和增强点云数据,为下一步的点云重构工作奠定基础。
从另外一个角度说,骨架信息的提取一直图形学领域的一个重要的研究课题。因为,不管是二维还是三维的形状,骨架都是对源数据的一个重要的描述特征。有了骨架信息后,可以大大帮助我们分析并操纵各种图形数据,比如基于骨架的形状匹配,骨架动画等。
从封闭的三维网格中提取骨架的技术已经较为成熟,但是从散乱点云中提取骨架却是个新兴的问题。从网格中提取骨架的技术并不能直接应用到点云骨架提取,因为点云中缺少网格的连接信息,并且现有的网格骨架提取技术并不能应对点云中的噪声,外点和缺失等问题。而目前能够从点云中提取曲线骨架的方法是2009年提出的方案,简称ROSA(rotational symmetry axis)。该方案提出了一个旋转对称轴的概念,并且假设输入模型的基本形状主要由圆柱体构成。假设已知一个二维的点集并且知道每个点的法向量,可以算出一个旋转对称中心点,最能表达该点集的中心位置。基于这个二维概念,ROSA方案首先为三维点云的每一个点找到一个最佳的切平面。然后,把切平面附近的输入点投影到该二维平面上,并算出该平面点集的旋转对称中心。把这些旋转对称中心点进行进一步收缩,并且对接合点区域作特殊处理后,得到最后的旋转对称中轴。另外,值得一提的是文献提出的基于拉普拉斯收缩的点云骨架提取技术。该技术假设输入的点云数据足够稠密和规整,从而可以建立基于点集的拉普拉斯连接关系,相当于构造了一个隐性网格。然后利用较为成熟的拉普拉斯网格骨架提取技术提取骨架。
然而,传统技术中的ROSA方案严重依赖点云法向量的信息,以及对物体基本形状的假设。在ROSA方案中,法向量信息是计算局部旋转对称中心的必要条件。而对基本形状是圆柱的假设,是它对缺失数据仍能较好提取中轴的关键。然而,法向量信息并不是直接从硬件设备中得到,而是由局部邻点的信息估算出来的值。对于有大量外点和噪声的情况,很可能估算出错误的法向量方向和朝向,因此,传统技术总的ROSA方案提取点云骨架的准确度不高。
发明内容
基于此,有必要提供一种能提高准确性的点云骨架提取方法。
一种点云骨架提取方法,包括:
获取输入的点云采样数据;
根据迭代公式进行点云收缩得到骨架分支,所述迭代公式为:
其中,
其中,i,j为正整数,J为点云采样数据的点集合,qj为点集合J中的采样点,I为采样点qj的邻域点集合,xi为邻域点集合I中的邻域点,xi'为邻域点集合I中不属于xi的邻域点,R为规整项,γi为权重系数,h为邻域点集合I的邻域半径,σi为分布系数;
所述根据迭代公式进行点云收缩的步骤为:
设置初始邻域半径 以hi=hi-1+h0/2的方式逐步扩大邻域半径,其中,dbb为点云包围盒对角线长度,|J|为点云采样数据的点集合中采样点的个数;
根据扩大的邻域半径以及所述迭代公式进行点云收缩;
所述分布系数σi的获取步骤为:根据PCA算法为采样点建立协方差矩阵;
根据下述公式得到分布系数:
其中,λ0i、λ1i和λ2i为所述协方差矩阵的三个特征值,且λ0i≤λ1i≤λ2i;
连接所述骨架分支得到点云骨架。
在其中一个实施例中,所述连接所述骨架分支得到点云骨架的步骤之后还包括:
对所述点云骨架进行平滑化和中心化处理。
此外,还有必要提供一种能提高准确性的点云骨架提取装置。
一种点云骨架提取装置,包括:
采样数据获取模块,用于获取输入的点云采样数据;
骨架分支生成模块,用于根据迭代公式进行点云收缩得到骨架分支,所述迭代公式为:
其中,
其中,i,j为正整数,J为点云采样数据的点集合,qj为点集合J中的采样点,I为采样点qj的邻域点集合,xi为邻域点集合I中的邻域点,xi'为邻域点集合I中不属于xi的邻域点,R为规整项,γi为权重系数,h为邻域点集合I的邻域半径,σi为分布系数;所述骨架分支生成模块还用于根据PCA算法为采样点建立协方差矩阵;
根据下述公式得到分布系数:
其中,λ0i、λ1i和λ2i为所述协方差矩阵的三个特征值,且λ0i≤λ1i≤λ2i;
所述骨架分支生成模块还用于设置初始邻域半径以hi=hi-1+h0/2的方式逐步扩大邻域半径,其中,dbb为点云包围盒对角线长度,|J|为点云采样数据的点集合中采样点的个数,根据扩大的邻域半径以及所述迭代公式进行点云收缩;
点云骨架生成模块,用于连接所述骨架分支得到点云骨架。
在其中一个实施例中,所述装置还包括点云骨架调整模块,用于对所述点云骨架进行平滑化和中心化处理。
上述点云骨架提取提取方法及装置,在传统的基于L1中值定理收缩点云骨架的基础上,增加了规整项,使得在采样点分布不均匀的情况下,能够依靠规整项的调节作用使得计算得到的中值点为接近真实环境的中值点,从而提高了准确性。
附图说明
图1为一个实施例中点云骨架提取方法的流程图;
图2为一个实施例中点云骨架提取装置的结构示意图;
图3为另一个实施例中点云骨架提取装置的结构示意图。
具体实施方式
在一个实施例中,如图1所示,一种点云骨架提取方法,该方法完全依赖于计算机程序,该计算机程序可运行于基于冯洛伊曼体系的计算机系统之上,该方法包括以下步骤:
步骤S102,获取输入的点云采样数据。
对点云数据采集设备采集到的散乱点集进行采样得到点云采样数据(散乱点集的采样点),通过对采样点的收缩从而生成点云骨架,可大大减少计算量。
在本实施例中,采样点数一般取输入的散乱点集中数据点数目的5%,采样位置随机,即随机抽取散乱点集中数据点的5%作为采样点。为了方便描述,下文统一称输入的散乱点集为Q,而下采样的点集为X。所谓点集的邻域有两个方面。一方面是指点集X自身的邻域I,即对于点集X中每个点 xi的邻域I为距离xi的长度小于h值(邻域半径)的邻点的集合。另一方面是指点集X在点集Q中的邻域J,即对于集合中每个点xi,距离xi的长度小于h值的邻点的集合。点集X的邻域的大小,由h值控制。
步骤S104,根据迭代公式进行点云收缩得到骨架分支,迭代公式为:
其中,
其中,i,j为正整数,J为点云采样数据的点集合,qj为点集合J中的采样点,I为采样点qj的邻域点集合,xi为邻域点集合I中的邻域点,xi'为邻域点集合I中不属于xi的邻域点,R为规整项,γi为权重系数,h为邻域点集合I的邻域半径,σi为分布系数。
传统技术中的L1中值定理为:
其计算得到的中值点x是到该点集中所有点的距离之和最小的点。而传统技术中的基于L1中值定理计算收缩点云的方法中,采样点可能分布不均匀,使得很多采样点都重叠汇集到附近的中值点,形成稀疏的分布,因此使得收缩得到的点云骨架不准确。例如,若采样点的分布接近一条直线,则收缩得到的骨架的则位置偏离较大
而在本实施例中,加入了规整项R,规整项R中包含有分布系数σi,使得对于分布不均匀的采样点可在规整项的调解下仍可得到准确的骨架位置。
在本实施例中,可根据PCA算法计算分布系数σi。
根据PCA算法为采样点建立协方差矩阵;
根据下述公式得到分布系数:
其中,λ0i、λ1i和λ2i为协方差矩阵的三个特征值,且λ0i≤λ1i≤λ2i。
PCA算法的目的是“降噪”和“去冗余”,协方差矩阵的主对角线上的元素是各个维度上的方差(即能量),其他元素是两两维度间的协方差(即相关性)。为了达到“降噪”目的,希望不同维度间的相关性尽可能小,即让协方差矩阵中非对角线元素都基本为零。这可利用线性代数中的矩阵对角化来实现。对角化后得到的矩阵,其对角线上是协方差矩阵的特征值。这特征值是各维度上的新方差,同时也代表各个维度本身应该拥有的能量。
可根据PCA算法对建立采样点的协方差矩阵:
其中,xi`为xi的均值。
由上述公式可看出,协方差矩阵Ci必然具有三个特征值,则可先将该三个特征值按照从小到大的顺序排列,依次为λ0、λ1和λ2。
在本实施例中,根据迭代公式进行点云收缩的步骤还可包括逐步扩大邻域半径,根据扩大的邻域半径以及迭代公式进行点云收缩。
进一步的,逐步扩大邻域半径的步骤可具体为:
设置初始邻域半径以hi=hi-1+h0/2的方式逐步扩大邻域半径,其中,dbb为点云包围盒对角线长度,|J|为点云采样数据的点集合中采样点的个数。
由于点云本身没有点与点之间的连接关系,所以邻域范围是绝大部分点云处理算法的一个重要的全局的固定的参数。但是,在本问题中,用一个固定的邻域范围是不可行的。对于稍微复杂的模型,有些地方的采样点需要考虑较小的邻域范围,否则会受不同形状部位的点的影响;有些则要考虑较大的邻域范围,这样才能把采样点收缩到正确的局部形状内部。因此,在本实施例中,可先用较小的邻域大小,即作为初始邻域大小进行收缩。收缩到一定程度后,让明显可以构成骨架分支的局部采样点集连成骨架,并且把这些采样点固定下来不再参与往后的收缩。接着,按一定的增长速率逐步增加邻域大小(默认为,其中),用新的邻域大小进行进一步的收缩和骨架提取。本实施例中,就是通过不断增加邻域大小,直到所有采样点都固定,形成最终骨架为止。
步骤S106,连接骨架分支得到点云骨架。
在本实施例中,连接骨架分支得到点云骨架的步骤可包括构建,延展和连接三个子过程。
在一定的邻域范围下,采样点集X经过多次迭代收缩后,有些局部点集逐渐形成明显的曲线骨架分布,把这些点集简单地逐个连起来即可形成一条新的骨架分支。在前述的迭代收缩的过程中,现有的骨架分支可以根据头尾的采样点分布进行延展或和其他骨架分支合并。
新骨架分支的构建。经过一定的收缩后,可先通过粗筛选,选出未被固定的采样点中有可能形成新骨架分支的候选点集,然后用约束规则确定是否能够形成新分支。粗筛选的方法是再次利用PCA算法求出每个采样点的σ值,然后取xi点的KNN(k-Nearest Neighbor algorithm,最邻近结点算法)点进行平滑,即,K=5。最后把平滑后的σ值大于0.9的,未被固定的采样点作为骨架分支构建候选点。于是,从σ值最大的候选点开始搜索,搜索方向是沿着PCA主方向(特征值最大的特征向量)的前后两端进行搜索,搜索范围是当前的邻域大小,如果发现在搜索范围内有候选点满足夹角规则(即两个相邻分段的夹角小于25°)那么把该候选点纳入当前候选分支。当候选分支的头尾两端都没有合适的候选点,搜索停止。如果搜索结果的候选分支包含的采样点数大于5,我们认为该分支可信,并把对应的采样点标记为固定点。否则把对应的采样点从候选点中排除。经过一轮搜索后,如果没有候选点剩余,则终止新骨架分支的构建,否则,继续从σ值最大的候选点进行搜索。
由于采用逐步递增的邻域大小,在不同邻域形成的骨架分支,很可能会丢失连接关系,因此,可引入辅助的桥接点。每当一个新的骨架分支形成,分支的头尾点就会连接一个桥接点。该点作为骨架分支的临时部分,但实际指向最近的未固定的采样点并参与进一步的收缩。桥接点一直跟踪最近的采样点直到出现以下情况会失效:1)跟踪的距离超出给定的阈值;2)跟踪的夹角(桥接点和对应分支的实际头尾点的连线,与该分支主方向形成的夹角)大于90°;3)局部遇到其他桥接点,并满足条件合并或连接分支。
基于桥接点的骨架分支延展过程。在更大的邻域大小的收缩下,已有骨架的头尾的采样点可能会重新满足夹角规则的约束。因此,在桥接点的帮助下,对应的采样点回纳入现有的骨架分支,并继续往前搜索直到最近的采样点不满足夹角规则为止,最后更新桥接点的位置。
基于桥接点的骨架分支连接过程。在一定的局部区域内,可能两个或多个桥接点相遇。只在增加邻域大小前,对这种情况进行处理。如果区域内有且只有两个桥接点,并且对应的两条分支的方向的夹角大于145°,则把这两条分支合并成一条,桥接点失效;否则不做任何处理。如果区域内有大于两个桥接点,求出这些桥接点的平均点并把对应分支的头尾移动到该平均点,桥接点失效。
进一步的,还可去除干扰采样点。在逐步扩大邻域半径并进行不同阶段的收缩的同时,可能会有少数采样点没有形成骨架分支却游离在已有骨架周围或头尾。因此,为提高算法鲁棒性和效率,对于这些靠近已有分支的采样点,做删除处理。
在本实施例中,连接骨架分支得到点云骨架的步骤之后还可对点云骨架进行平滑化和中心化处理。
平滑化的过程如下:
第一步,对于每条骨架分支的除头尾外的每个结点v1,可以算出它和前后节点v0与v2之间对应的夹角,如果该夹角大于50°,对该节点进行一维拉普拉斯平滑,即v1=v0/4+v1/2+v2/4。
第二步,对每条分支进行四点插值细分,直到该分支最长的分段小于一个阈值。细分的方法如下,假设v0、v1、v2、v3是分支上连续四个点,插入新点的规则如下。
1.v1,v2的新插入点v=(-v0+9v1+9v2-v3)/16。
2.如果v1是端点,即v1之前没有点:v=(3v1+6v2-v3)/8。
3.如果v2是端点,即v2之后没有点:v=(-v0+6v1+3v2)/8。
第三步,对细分后的分支进行下采样,即从取出第一个点开始,逐步取出下一个点作为采样结点,该结点刚好满足和前一个取出点的距离大于一个阈值。最终得到一个平滑的,并且分段长度基本相同的一维曲线骨架。
中心化的过程如下:
对于骨架的每一条分支的除头尾外的结点,都可以简单得到一个切平面。该平面经过结点并且垂直于该结点和下一个结点连线方向的向量。把结点对应的邻近源点集Q上的点投影到该切平面上,然后用几何椭圆拟合的方法去用一个椭圆逼近该二维点集。如果得到的椭圆的拟合误差小于默认阈值,即把该结点移动到对应椭圆的中心。如果一条分支上有多个结点被重新中心化,对该分支再次做平滑化。需要说明的是,这里的切平面邻域的半径,是由该结点的在收缩过程中第一次固定下来时候对应的邻域半径所决定。
在一个实施例中,如图2所示,一种危险信息通知装置,包括采样数据获取模块102、骨架分支生成模块104以及点云骨架生成模块106,其中:
采样数据获取模块102,用于获取输入的点云采样数据;
骨架分支生成模块104,用于根据迭代公式进行点云收缩得到骨架分支,所述迭代公式为:
其中,
其中,i,j为正整数,J为点云采样数据的点集合,qj为点集合J中的采样点,I为采样点qj的邻域点集合,xi为邻域点集合I中的邻域点,xi'为邻域点集合I中不属于xi的邻域点,R为规整项,γi为权重系数,h为邻域点集合I的邻域半径,σi为分布系数;
点云骨架生成模块106,用于连接所述骨架分支得到点云骨架。
在一个实施例中,骨架分支生成模块104还用于根据PCA算法为采样点建立协方差矩阵;
根据下述公式得到分布系数:
其中,λ0i、λ1i和λ2i为所述协方差矩阵的三个特征值,且λ0i≤λ1i≤λ2i。
在一个实施例中,骨架分支生成模块104还用于逐步扩大邻域半径,根据扩大的邻域半径以及所述迭代公式进行点云收缩。
在一个实施例中,骨架分支生成模块104还用于设置初始邻域半径以hi=hi-1+h0/2的方式逐步扩大邻域半径,其中,dbb为点云包围盒对角线长度,|J|为点云采样数据的点集合中采样点的个数。
在一个实施例中,如图3所示,点云骨架提取装置还包括点云骨架调整模块108,用于对点云骨架进行平滑化和中心化处理。
上述点云骨架提取提取方法及装置,在传统的基于L1中值定理收缩点云骨架的基础上,增加了规整项,使得在采样点分布不均匀的情况下,能够依靠规整项的调节作用使得计算得到的中值点为接近真实环境的中值点,从而提高了准确性。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random Access Memory,RAM)等。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。
Claims (4)
1.一种点云骨架提取方法,包括:
获取输入的点云采样数据;
根据迭代公式进行点云收缩得到骨架分支,所述迭代公式为:
其中,
其中,i,j为正整数,J为点云采样数据的点集合,qj为点集合J中的采样点,I为采样点qj的邻域点集合,xi为邻域点集合I中的邻域点,xi'为邻域点集合I中不属于xi的邻域点,R为规整项,γi为权重系数,h为邻域点集合I的邻域半径,σi为分布系数;
所述根据迭代公式进行点云收缩的步骤为:
设置初始邻域半径以hi=hi-1+h0/2的方式逐步扩大邻域半径,其中,dbb为点云包围盒对角线长度,|J|为点云采样数据的点集合中采样点的个数;
根据扩大的邻域半径以及所述迭代公式进行点云收缩;
所述分布系数σi的获取步骤为:根据PCA算法为采样点建立协方差矩阵;
根据下述公式得到分布系数:
其中,λ0i、λ1i和λ2i为所述协方差矩阵的三个特征值,且λ0i≤λ1i≤λ2i;
连接所述骨架分支得到点云骨架。
2.根据权利要求1所述的点云骨架提取方法,其特征在于,所述连接所述骨架分支得到点云骨架的步骤之后还包括:
对所述点云骨架进行平滑化和中心化处理。
3.一种点云骨架提取装置,其特征在于,包括:
采样数据获取模块,用于获取输入的点云采样数据;
骨架分支生成模块,用于根据迭代公式进行点云收缩得到骨架分支,所述迭代公式为:
其中,
其中,i,j为正整数,J为点云采样数据的点集合,qj为点集合J中的采样点,I为采样点qj的邻域点集合,xi为邻域点集合I中的邻域点,xi'为邻域点集合I中不属于xi的邻域点,R为规整项,γi为权重系数,h为邻域点集合I的邻域半径,σi为分布系数;
所述骨架分支生成模块还用于根据PCA算法为采样点建立协方差矩阵;
根据下述公式得到分布系数:
其中,λ0i、λ1i和λ2i为所述协方差矩阵的三个特征值,且λ0i≤λ1i≤λ2i;
所述骨架分支生成模块还用于设置初始邻域半径以hi=hi-1+h0/2的方式逐步扩大邻域半径,其中,dbb为点云包围盒对角线长度,|J|为点云采样数据的点集合中采样点的个数,根据扩大的邻域半径以及所述迭代公式进行点云收缩;点云骨架生成模块,用于连接所述骨架分支得到点云骨架。
4.根据权利要求3所述的点云骨架提取装置,其特征在于,所述装置还包括点云骨架调整模块,用于对所述点云骨架进行平滑化和中心化处理。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310196243.2A CN103268631B (zh) | 2013-05-23 | 2013-05-23 | 点云骨架提取方法及装置 |
PCT/CN2013/083441 WO2014187046A1 (zh) | 2013-05-23 | 2013-09-13 | 点云骨架提取方法及装置 |
US14/378,976 US9390552B1 (en) | 2013-05-23 | 2013-09-13 | Method and device for extracting skeleton from point cloud |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310196243.2A CN103268631B (zh) | 2013-05-23 | 2013-05-23 | 点云骨架提取方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103268631A CN103268631A (zh) | 2013-08-28 |
CN103268631B true CN103268631B (zh) | 2015-09-30 |
Family
ID=49012258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310196243.2A Active CN103268631B (zh) | 2013-05-23 | 2013-05-23 | 点云骨架提取方法及装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9390552B1 (zh) |
CN (1) | CN103268631B (zh) |
WO (1) | WO2014187046A1 (zh) |
Families Citing this family (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2816966B1 (en) | 2012-02-22 | 2023-10-25 | Veran Medical Technologies, Inc. | Steerable surgical catheter comprising a biopsy device at the distal end portion thereof |
CN103268631B (zh) | 2013-05-23 | 2015-09-30 | 中国科学院深圳先进技术研究院 | 点云骨架提取方法及装置 |
CN103778654A (zh) * | 2013-12-10 | 2014-05-07 | 深圳先进技术研究院 | 一种三维几何体特征线引导的表面光滑向量场计算方法 |
CN103745497B (zh) * | 2013-12-11 | 2017-05-10 | 中国科学院深圳先进技术研究院 | 植物生长建模方法和系统 |
CN103714555B (zh) * | 2013-12-13 | 2017-06-13 | 中国科学院深圳先进技术研究院 | 一种基于运动轨迹的四维运动点云分割与重建方法 |
CN104156730B (zh) * | 2014-07-25 | 2017-12-01 | 山东大学 | 一种基于骨架的抗噪声汉字特征提取方法 |
CN105373814B (zh) * | 2014-08-26 | 2019-04-02 | 南京林业大学 | 一种基于激光点云数据的真实阔叶树器官分类识别方法 |
CN104408765B (zh) * | 2014-11-11 | 2017-01-11 | 中国科学院深圳先进技术研究院 | 植物扫描与重建方法 |
CN104966287B (zh) * | 2015-06-08 | 2017-08-08 | 浙江大学 | 层次化的多片点云刚性配准方法 |
CN105005995B (zh) * | 2015-07-29 | 2017-07-25 | 武汉大学 | 一种计算三维点云模型骨骼的方法 |
CN105405162B (zh) * | 2015-10-16 | 2017-12-05 | 北京师范大学 | 基于局部结构和方向感知的树点云三维重建方法 |
CN106504331B (zh) * | 2016-09-27 | 2021-05-18 | 西安科技大学 | 基于三维模型检索的牙齿建模方法 |
CN106780458B (zh) * | 2016-12-09 | 2020-04-28 | 重庆邮电大学 | 一种点云骨架提取方法及装置 |
TWI657691B (zh) * | 2017-05-24 | 2019-04-21 | 鈺立微電子股份有限公司 | 矯正原始立體掃描結果的錯誤法向量的裝置及其相關方法 |
CN108280833B (zh) * | 2018-01-18 | 2021-09-24 | 华南农业大学 | 一种植物根系分叉特征的骨架提取方法 |
CN109064471B (zh) * | 2018-07-18 | 2021-09-03 | 中北大学 | 一种基于骨架的三维点云模型分割方法 |
CN109325993B (zh) * | 2018-08-10 | 2023-01-06 | 华北电力大学(保定) | 一种基于类八叉树索引的显著性特征强化采样方法 |
CN109887009B (zh) * | 2019-01-24 | 2022-12-09 | 西北大学 | 一种点云局部匹配方法 |
CN110827233B (zh) * | 2019-08-29 | 2022-06-14 | 杭州电子科技大学 | 一种牙齿三维点云数据表面窝沟区域的提取方法 |
WO2021114026A1 (zh) * | 2019-12-09 | 2021-06-17 | 深圳大学 | 一种基于局部参考坐标系的3d形状匹配方法及装置 |
CN111462202B (zh) * | 2020-04-08 | 2022-09-02 | 中国科学技术大学 | 非刚性注册方法及系统 |
CN111666946A (zh) * | 2020-05-26 | 2020-09-15 | 东华大学 | 基于点云过分割与面片生长的植株点云单叶片分割方法 |
CN112102178A (zh) * | 2020-07-29 | 2020-12-18 | 深圳市菲森科技有限公司 | 一种点云保特征去噪方法、装置、电子设备及存储介质 |
CN111968089B (zh) * | 2020-08-15 | 2024-06-21 | 晋江市博感电子科技有限公司 | 一种基于最大内切球机制的l1中值骨架提取方法 |
CN112465832B (zh) * | 2020-11-25 | 2024-04-16 | 重庆大学 | 一种基于双目视觉的单面树木点云骨架线提取方法及系统 |
CN112686799B (zh) * | 2020-12-25 | 2022-05-10 | 燕山大学 | 一种基于法向量和l1中值的环形锻件截面形线提取方法 |
CN112712509B (zh) * | 2020-12-31 | 2023-09-01 | 重庆大学 | 基于点云的树木参数获取方法、生长评估方法、装置及系统 |
CN112967333B (zh) * | 2021-02-04 | 2024-02-09 | 重庆大学 | 基于等级划分的复杂点云骨架提取方法及系统 |
CN112802089B (zh) * | 2021-02-05 | 2023-09-01 | 重庆大学 | 一种基于分叉数量自动估计的点云骨架线提取方法及系统 |
CN113298919A (zh) * | 2021-04-14 | 2021-08-24 | 江苏大学 | 一种三维植物点云模型的骨架提取方法 |
CN114419276B (zh) * | 2022-01-19 | 2024-05-03 | 厦门大学 | 基于最优传输和聚类的曲线骨架提取方法 |
CN115841581B (zh) * | 2023-02-20 | 2023-05-16 | 宜科(天津)电子有限公司 | 一种门板骨架的特征点提取方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8004517B1 (en) * | 2005-06-24 | 2011-08-23 | Geomagic, Inc. | Methods, apparatus and computer program products that model three-dimensional surface structures |
CN102467753A (zh) * | 2010-11-04 | 2012-05-23 | 中国科学院深圳先进技术研究院 | 基于骨架配准的时变点云重建方法及系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8064731B2 (en) * | 2006-11-13 | 2011-11-22 | Siemens Audiologische Technik Gmbh | Generalized rigid alignment of 3D ear impression models |
CN101866495B (zh) * | 2010-06-01 | 2012-03-28 | 中国科学院自动化研究所 | 基于骨架点云的树木建模方法 |
CN101887596B (zh) | 2010-06-01 | 2013-02-13 | 中国科学院自动化研究所 | 树木点云数据基于分割和自动生长的三维模型重建方法 |
CN103098100B (zh) * | 2010-12-03 | 2016-01-20 | 中国科学院自动化研究所 | 基于感知信息的三维模型形状分析方法 |
EP2674913B1 (en) * | 2012-06-14 | 2014-07-23 | Softkinetic Software | Three-dimensional object modelling fitting & tracking. |
CN103268631B (zh) | 2013-05-23 | 2015-09-30 | 中国科学院深圳先进技术研究院 | 点云骨架提取方法及装置 |
-
2013
- 2013-05-23 CN CN201310196243.2A patent/CN103268631B/zh active Active
- 2013-09-13 US US14/378,976 patent/US9390552B1/en active Active
- 2013-09-13 WO PCT/CN2013/083441 patent/WO2014187046A1/zh active Application Filing
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8004517B1 (en) * | 2005-06-24 | 2011-08-23 | Geomagic, Inc. | Methods, apparatus and computer program products that model three-dimensional surface structures |
CN102467753A (zh) * | 2010-11-04 | 2012-05-23 | 中国科学院深圳先进技术研究院 | 基于骨架配准的时变点云重建方法及系统 |
Non-Patent Citations (4)
Title |
---|
Consolidation of Unorganized Point Clouds for Surface Reconstruction;Hui Huang et al.;《Proceedings of ACM SIGGRAPH Asia 2009》;20091231;第1-7页 * |
Curve Skeleton Extraction from Incomplete Point Cloud;Andrea Tagliasacchi et al.;《Proceedings of ACM SIGGRAPH 2009》;20090803;第28卷(第3期);第1-9页 * |
Robust statistics in data analysis-A review Basic concepts;M. Daszykowski et al.;《Chemometrics and Intelligent Laboratory Systems》;20070215;第85卷(第2期);第203-119页 * |
基于Lobe表示的树木快速3D建模;程章林 等;《计算机辅助设计与图形学学报》;20120131;第24卷(第1期);第2-4页 * |
Also Published As
Publication number | Publication date |
---|---|
WO2014187046A1 (zh) | 2014-11-27 |
US20160203636A1 (en) | 2016-07-14 |
CN103268631A (zh) | 2013-08-28 |
US9390552B1 (en) | 2016-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103268631B (zh) | 点云骨架提取方法及装置 | |
Patil et al. | An adaptive approach for the reconstruction and modeling of as-built 3D pipelines from point clouds | |
Saha et al. | A survey on skeletonization algorithms and their applications | |
US8041141B2 (en) | Method and software for shape representation with curve skeletons | |
KR100810294B1 (ko) | 3차원 메쉬 데이터의 특징-유지 간략화 방법 | |
WO2012071688A1 (zh) | 基干感知信息的三维模型形狀分析方法 | |
CN112233249B (zh) | 基于密集点云的b样条曲面拟合方法及装置 | |
CN106780458B (zh) | 一种点云骨架提取方法及装置 | |
CN104851126B (zh) | 基于广义圆柱体的三维模型分割方法及装置 | |
Nana et al. | Automatic reconstruction of beam structures from 3D topology optimization results | |
Governi et al. | 3D geometry reconstruction from orthographic views: A method based on 3D image processing and data fitting | |
Zhu et al. | Computing a compact spline representation of the medial axis transform of a 2D shape | |
Cohen et al. | Virtual reconstruction of archeological vessels using expert priors and intrinsic differential geometry information | |
Chang et al. | Surface reconstruction from point clouds by transforming the medial scaffold | |
Hu et al. | Surface segmentation for polycube construction based on generalized centroidal Voronoi tessellation | |
Tabib et al. | Learning-based hole detection in 3D point cloud towards hole filling | |
Monga et al. | Representing geometric structures in 3D tomography soil images: Application to pore-space modeling | |
Hu et al. | Curve skeleton extraction from 3D point clouds through hybrid feature point shifting and clustering | |
Ngom et al. | 3D shape extraction segmentation and representation of soil microstructures using generalized cylinders | |
Zhu et al. | Propagation strategies for stereo image matching based on the dynamic triangle constraint | |
Verykokou et al. | A Comparative analysis of different software packages for 3D Modelling of complex geometries | |
Fu et al. | Terrestrial laser scanning assisted dimensional quality assessment for space frame components | |
Clémot et al. | Neural skeleton: Implicit neural representation away from the surface | |
Sahillioğlu et al. | Coarse-to-fine surface reconstruction from silhouettes and range data using mesh deformation | |
Cao et al. | Patch layout generation by detecting feature networks |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |