CN109033028A - 一种点云主曲率计算方法 - Google Patents
一种点云主曲率计算方法 Download PDFInfo
- Publication number
- CN109033028A CN109033028A CN201810682212.0A CN201810682212A CN109033028A CN 109033028 A CN109033028 A CN 109033028A CN 201810682212 A CN201810682212 A CN 201810682212A CN 109033028 A CN109033028 A CN 109033028A
- Authority
- CN
- China
- Prior art keywords
- point
- principal curvatures
- cloud
- calculation
- curved surface
- 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.)
- Pending
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims abstract description 14
- 238000010276 construction Methods 0.000 claims abstract description 8
- 230000003044 adaptive effect Effects 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012407 engineering method Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/20—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring contours or curvatures, e.g. determining profile
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明提供一种点云主曲率计算方法,包括以下步骤:构造自由曲面,在曲面上采集点云得到Q1并计算Q1中每个点的法矢;根据得到的法矢,在Q1中选定一个点p,在其k邻域内计算质心,并由质心计算协方差矩阵;根据得到的协方差矩阵,计算其特征值和特征向量;比较得到最大和次大特征值及其对应的特征向量,通过计算得到主曲率;根据得到的主曲率构造各向异性高斯核函数,得到权函数;根据权函数构造局部加权法矢,由加权法矢得到局部各向异性最小二乘曲面;循环操作至最后一个点的计算。本发明提供的一种点云主曲率计算方法,通过使用两个主曲率作为高斯核函数的参数,增加利用自适应最小二乘曲面拟合点云在平坦区域的光滑度,有效突出了细节部分。
Description
技术领域
本发明涉及工程零件表面点云测量领域,更具体的,涉及一种点云主曲率计算方法。
背景技术
三维数据测量现已广泛应用于复杂型面零件的加工质量检测、运行状态监测及快速检修过程中,以及时进行外形评价并做出决策,避免产生大量加工废品或运行安全事故。为此,工程中常通过反求工程方法,利用激光扫描仪或三坐标测量机等对工件外形进行测量,并对测量数据进行对齐与曲面重建,通过与工件CAD模型的比较分析以评价其加工质量或规划修复区域加工工艺。由于受测量工件外形轮廓尺寸及传感器测量范围限制,单一视角很难获取工件外形完整测量数据,往往需要多个视角甚至多个传感器测量才可获取完整测量数据。
传统的MLS曲面采用单一固定半径高斯核函数定义方法,在获得的拟合曲面在平坦区域振荡,在高曲率区域欠逼近,无法突出细节,且无法准确反映测量工件上任意点处沿不同方向形状变化较大的特征,此将使逼近的MLS曲面存在一定的固有逼近误差,并最终影响测量数据到MLS曲面的对齐精度。
发明内容
本发明为克服现有的MLS曲面采用单一固定半径高斯核函数定义方法,存在无法突出细节,无法准确反映测量工件上任意点处沿不同方向形状变化较大的特征的技术缺陷,提供一种点云主曲率计算方法。
为解决上述技术问题,本发明的技术方案如下:
一种点云主曲率计算方法,其特征在于,包括以下步骤:
S1:构造自由曲面,在曲面上采集点云得到Q1并计算Q1中每个点的法矢;
S2:根据得到的法矢,在Q1中选定一个点p,在其k邻域内计算质心,并由质心计算协方差矩阵;
S3:根据得到的协方差矩阵,计算其特征值和特征向量;
S4:比较得到最大和次大特征值及其对应的特征向量,通过计算得到主曲率;
S5:根据得到的主曲率构造各向异性高斯核函数,得到权函数;
S6:根据权函数构造局部加权法矢,由加权法矢得到局部各向异性最小二乘曲面;
S7:判断是否为最后一个点,若是,则结束;若否,则执行步骤S2。
其中,所述步骤S2包括以下步骤:
S21:对于点云Q1中的点p,利用KD-Tree算法获取到点p的k邻域,以点p为中心,以R为半径创建球域,设球域坐标系为(β,α,γ),球域内的点pj可表示为:
其中,α∈[0,2π],u∈[-1,1],r∈[0,R],在整个球域内生成nb个均匀分布点;
S22:将得到的nb个点pj投影到由p点及其邻域构建的移动最小二乘曲面上,得到投影点p′j,若则p′j在区域A内部,即球域与移动最小二乘曲面的相交区域的内部,否则在外部;
S23:收集区域A内部的点构成点集用nBin表示B内点的数目,则区域A的体积Vojb可表示为:其中,Vb是以点p为中心,半径为R的球的体积,大小为4πR3/3;
S24:根据得到的区域A的体积,计算质心b和协方差矩阵J(A),其计算公式为:
其中,在所述步骤S4中,根据协方差矩阵J(A)得到最大及次大的特征值,记为Mb1,Mb2以及其所对应的方向向量v1、v2,p点处的主曲率可表示为 其计算公式为:
其中,所述步骤S5构造各向异性高斯核函数的公式具体为:
其中,θ(||x-pi||)为权函数。
其中,构造局部加权法矢,其计算公式为:
其中,n(x)表示每一个点的法矢,根据加权法矢构造局部各向异性最小二乘曲面,其计算公式为:
与现有技术相比,本发明技术方案的有益效果是:
本发明提供的一种点云主曲率计算方法,通过使用两个主曲率作为高斯核函数的参数,增加利用自适应最小二乘曲面拟合点云在平坦区域的光滑度,有效突出了细节部分。
附图说明
图1为一种点云主曲率计算方法方法流程图。
图2为创建的球域示意图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
如图1、图2所示,一种点云主曲率计算方法,包括以下步骤:
S1:构造自由曲面,在曲面上采集点云得到Q1并计算Q1中每个点的法矢;
S2:根据得到的法矢,在Q1中选定一个点p,在其k邻域内计算质心,并由质心计算协方差矩阵;
S3:根据得到的协方差矩阵,计算其特征值和特征向量;
S4:比较得到最大和次大特征值及其对应的特征向量,通过计算得到主曲率;
S5:根据得到的主曲率构造各向异性高斯核函数,得到权函数;
S6:根据权函数构造局部加权法矢,由加权法矢得到局部各向异性最小二乘曲面;
S7:判断是否为最后一个点,若是,则结束;若否,则执行步骤S2。
更具体的,所述步骤S2包括以下步骤:
S21:对于点云Q1中的点p,利用KD-Tree算法获取到点p的k邻域,以点p为中心,以R为半径创建球域,设球域坐标系为(β,α,γ),球域内的点pj可表示为:
其中,α∈[0,2π],u∈[-1,1],r∈[0,R],在整个球域内生成nb个均匀分布点;
S22:将得到的nb个点pj投影到由p点及其邻域构建的移动最小二乘曲面上,得到投影点p′j,若则p′j在区域A内部,即球域与移动最小二乘曲面的相交区域的内部,否则在外部;
S23:收集区域A内部的点构成点集用nBin表示B内点的数目,则区域A的体积Vojb可表示为:其中,Vb是以点p为中心,半径为R的球的体积,大小为4πR3/3;
S24:根据得到的区域A的体积,计算质心b和协方差矩阵J(A),其计算公式为:
更具体的,在所述步骤S4中,根据协方差矩阵J(A)得到最大及次大的特征值,记为Mb1,Mb2以及其所对应的方向向量v1、v2,p点处的主曲率可表示为其计算公式为:
更具体的,所述步骤S5构造各向异性高斯核函数的公式具体为:
其中,θ(||x-pi||)为权函数。
更具体的,构造局部加权法矢,其计算公式为:
其中,n(x)表示每一个点的法矢,根据加权法矢构造局部各向异性最小二乘曲面,其计算公式为:
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (5)
1.一种点云主曲率计算方法,其特征在于,包括以下步骤:
S1:构造自由曲面,在曲面上采集点云得到Q1并计算Q1中每个点的法矢;
S2:根据得到的法矢,在Q1中选定一个点p,在其k邻域内计算质心,并由质心计算协方差矩阵;
S3:根据得到的协方差矩阵,计算其特征值和特征向量;
S4:比较得到最大和次大特征值及其对应的特征向量,通过计算得到主曲率;
S5:根据得到的主曲率构造各向异性高斯核函数,得到权函数;
S6:根据权函数构造局部加权法矢,由加权法矢得到局部各向异性最小二乘曲面;
S7:判断是否为最后一个点,若是,则结束;若否,则执行步骤S2。
2.根据权利要求1所述的一种点云主曲率计算方法,其特征在于,步骤S2包括以下步骤:
S21:对于点云Q1中的点p,利用KD-Tree算法获取到点p的k邻域,以点p为中心,以R为半径创建球域,设球域坐标系为(β,α,γ),球域内的点pj可表示为:
其中,α∈[0,2π],u∈[-1,1],r∈[0,R],在整个球域内生成nb个均匀分布点;
S22:将得到的nb个点pj投影到由p点及其邻域构建的移动最小二乘曲面上,得到投影点p′j,若则p′j在区域A内部,即球域与移动最小二乘曲面的相交区域的内部,否则在外部;
S23:收集区域A内部的点构成点集用nBin表示B内点的数目,则区域A的体积Vojb可表示为:其中,Vb是以点p为中心,半径为R的球的体积,大小为4πR3/3;
S24:根据得到的区域A的体积,计算质心b和协方差矩阵J(A),其计算公式为:
3.根据权利要求2所述的一种点云主曲率计算方法,其特征在于,在所述步骤S4中,根据协方差矩阵J(A)得到最大及次大的特征值,记为Mb1,Mb2以及其所对应的方向向量v1、v2,p点处的主曲率可表示为其计算公式为:
4.根据权利要求3所述的一种点云主曲率计算方法,其特征在于,所述步骤S5构造各向异性高斯核函数的公式具体为:
其中,θ(||x-pi||)为权函数。
5.根据权利要求4所述的一种点云主曲率计算方法,其特征在于,构造局部加权法矢,其计算公式为:
其中,n(x)表示每一个点的法矢,根据加权法矢构造局部各向异性最小二乘曲面,其计算公式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810682212.0A CN109033028A (zh) | 2018-06-27 | 2018-06-27 | 一种点云主曲率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810682212.0A CN109033028A (zh) | 2018-06-27 | 2018-06-27 | 一种点云主曲率计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109033028A true CN109033028A (zh) | 2018-12-18 |
Family
ID=65520585
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810682212.0A Pending CN109033028A (zh) | 2018-06-27 | 2018-06-27 | 一种点云主曲率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109033028A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110009029A (zh) * | 2019-03-28 | 2019-07-12 | 北京智行者科技有限公司 | 基于点云分割的特征匹配方法 |
CN111105381A (zh) * | 2020-02-07 | 2020-05-05 | 武汉玄景科技有限公司 | 一种基于球面模型的稠密点云平滑方法 |
CN113327275A (zh) * | 2021-06-18 | 2021-08-31 | 哈尔滨工业大学 | 一种基于多约束点到局部曲面投影的点云双视角精配准方法 |
-
2018
- 2018-06-27 CN CN201810682212.0A patent/CN109033028A/zh active Pending
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110009029A (zh) * | 2019-03-28 | 2019-07-12 | 北京智行者科技有限公司 | 基于点云分割的特征匹配方法 |
CN111105381A (zh) * | 2020-02-07 | 2020-05-05 | 武汉玄景科技有限公司 | 一种基于球面模型的稠密点云平滑方法 |
CN111105381B (zh) * | 2020-02-07 | 2023-03-31 | 武汉玄景科技有限公司 | 一种基于球面模型的稠密点云平滑方法 |
CN113327275A (zh) * | 2021-06-18 | 2021-08-31 | 哈尔滨工业大学 | 一种基于多约束点到局部曲面投影的点云双视角精配准方法 |
CN113327275B (zh) * | 2021-06-18 | 2022-11-08 | 哈尔滨工业大学 | 一种基于多约束点到局部曲面投影的点云双视角精配准方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107301648B (zh) | 基于重叠区域边界角度的冗余点云去除方法 | |
CN101887586B (zh) | 基于图像轮廓尖锐度的自适应角点检测方法 | |
CN109033028A (zh) | 一种点云主曲率计算方法 | |
CN102207441B (zh) | 一种静态接触角和动态接触角的自动检测方法 | |
CN104200212B (zh) | 一种基于机载LiDAR数据的建筑物外边界线提取方法 | |
CN108549087A (zh) | 一种基于激光雷达的在线检测方法 | |
CN103886593A (zh) | 一种基于三维点云曲面圆孔检测方法 | |
CN106874580A (zh) | 一种基于点云数据的弯管模型重建方法 | |
CN101645170A (zh) | 多视点云精确配准方法 | |
CN104123746A (zh) | 一种三维扫描点云中实时法向量的计算方法 | |
CN103940356A (zh) | 一种基于三维激光扫描技术的建筑物整体变形监测方法 | |
CN103617328A (zh) | 一种飞机三维姿态解算方法 | |
CN104880160B (zh) | 一种工件表面轮廓的二维激光实时探测方法 | |
CN108805143B (zh) | 一种基于切片点云的数据特征提取方法 | |
CN103033253A (zh) | 一种塔式起重机结构的非接触振动检测方法 | |
CN108416787A (zh) | 应用于机器视觉检测的工件直线边缘定位方法 | |
CN105469404A (zh) | 一种基于三维点云数据的旋转体拟合方法及装置 | |
CN113031005A (zh) | 基于激光雷达的吊车动态障碍物识别方法 | |
CN107633523A (zh) | 基于点云的提取建筑特征线方法和*** | |
CN103544728A (zh) | 一种在实景地图中基于法向量的被探面过滤方法 | |
CN109191583B (zh) | 一种基于各向异性mls曲面精确对齐方法 | |
CN104517318A (zh) | 三维量测模拟取点***及方法 | |
CN109712190A (zh) | 三维激光扫描仪与三维激光雷达的外参标定方法 | |
CN109141266A (zh) | 一种钢结构测量方法及*** | |
CN107180409A (zh) | 一种针对弯曲骨架型对象三维点云的广义圆柱投影方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181218 |
|
RJ01 | Rejection of invention patent application after publication |