[发明专利]一种点云主曲率计算方法

说明书

本发明提供一种点云主曲率计算方法，包括以下步骤：构造自由曲面，在曲面上采集点云得到Q1并计算Q1中每个点的法矢；根据得到的法矢，在Q1中选定一个点p，在其k邻域内计算质心，并由质心计算协方差矩阵；根据得到的协方差矩阵，计算其特征值和特征向量；比较得到最大和次大特征值及其对应的特征向量，通过计算得到主曲率；根据得到的主曲率构造各向异性高斯核函数，得到权函数；根据权函数构造局部加权法矢，由加权法矢得到局部各向异性最小二乘曲面；循环操作至最后一个点的计算。本发明提供的一种点云主曲率计算方法，通过使用两个主曲率作为高斯核函数的参数，增加利用自适应最小二乘曲面拟合点云在平坦区域的光滑度，有效突出了细节部分。

技术领域

本发明涉及工程零件表面点云测量领域，更具体的，涉及一种点云主曲率计算方法。

背景技术

三维数据测量现已广泛应用于复杂型面零件的加工质量检测、运行状态监测及快速检修过程中，以及时进行外形评价并做出决策，避免产生大量加工废品或运行安全事故。为此，工程中常通过反求工程方法，利用激光扫描仪或三坐标测量机等对工件外形进行测量，并对测量数据进行对齐与曲面重建，通过与工件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为半径创建球域，设球域坐标系为(β,α,γ)，球域内的点p_j可表示为：

其中，α∈[0,2π],u∈[-1,1],r∈[0,R]，在整个球域内生成n_b个均匀分布点；

S22：将得到的n_b个点p_j投影到由p点及其邻域构建的移动最小二乘曲面上，得到投影点p′_j，若则p′_j在区域A内部，即球域与移动最小二乘曲面的相交区域的内部，否则在外部；