本发明提供一种实物表面样点棱边特征曲面拓扑逼近重建方法，属于产品逆向工程技术领域。
背景技术：
工业产品逆向工程中，实物表面棱边特征重建结果的形位精度对后续的设计分析和再制造有重要的影响。由于曲面采样数据通常为缺乏任何特征信息的散乱点集，在基于此类数据的曲面重建过程中，目前曲面拓扑重建技术未能实现棱边特征的精确重建，因此棱边特征区域的正确重建问题一直深受关注。
曲面重建经典算法通常假设采样数据来自光滑曲面，因而难以用于重建含棱边特征的曲面。t.k.dey等基于高斯加权拉普拉斯图和reeb图相组合来识别和重建棱边特征曲线，在delaunay剖分过程中，将特征曲线上采样点和非尖锐特征点加权保护，以提高棱边特征重建形位精度。weber等通过识别特征点局部逼近特征曲线，将特征点用作修改mls投影点的局部邻域指标，并对局部点邻域进行分段和增强，使分段光滑的表面保持所有的尖锐特征。lipmany等同样对移动最小二乘法(mls)投影方法上改进，算法定义检测奇异性指标域(sif)，使用sif来指导mls框架中的局部多项式拟合来重建曲面，通过使用连续的奇异性指示域，使得算法能够对奇点精细重构。
综上所述，现有的棱边特征区域重建方法均以优化棱边特征为目标，在重建过程中对棱边特征点加以保护来提高重建形位精度。但是当曲面采样数据含有噪声且为非均匀采样数据时，易出现凹痕等错误，导致棱边形位精度难以保证。因此，提供一种能够重建棱边特征且对非均匀采样数据有良好的适应性的重建方法已成为本领域技术人员亟待解决的技术问题。
技术实现要素：
本发明要解决的技术问题是：克服现有技术不足，提供一种实物表面样点棱边特征曲面重建方法，快速、准确重建棱边曲面特征。
为解决上述技术问题，本发明所采用的技术方案是一种实物表面样点棱边特征曲面拓扑逼近重建方法，其特征在于步骤依次为：(1)读取数据点集g，以尖锐特征区域样点为界限，基于平坦区域连通性，对点集进行分割处理，获取分割后点集集合gi,i＝1,2,3，…；(2)获取gi的边界点邻域点集并进行增益优化，以边界点法向约束优化后邻域点集的均值点对点集集合gi进行扩展；(3)基于零等值面重建算法对扩展后的点集进行曲面重建并细分，得到曲面集合h(gi)；(4)提取h(gi)中各曲面相交线，将交线段端点坐标和相交三角面片内部端点作为细分点进行剖分重建棱边特征曲面；(5)输出棱边特征曲面和平坦区域曲面，曲面重建完成。
为实现发明目的，所述实物表面样点棱边特征曲面拓扑逼近重建方法，其特征在于步骤(1)中所述以尖锐特征区域样点为界限，基于平坦区域连通性，对点集进行分割处理方法，具体步骤为：(1)读取采样数据点集g，采用局部邻域拟合法计算点集法向并进行法向统一；(2)估算采样数据点云曲率；(3)设定采样点邻域高斯映射聚类算法阈值剔除采样数据尖锐特征点；(4)对剩余采样点按照曲率值大小进行排序；(5)将曲率值最小的采样点作为种子点增长算法的起始种子点；(6)查询种子点近邻点，计算近邻点的法向与种子点法向之间的角度。若角度小于算法设定的阈值，把该种子点从排序列表中删除；否则执行步骤(8)；(7)若种子点邻近点的曲率值小于设定阈值，则将其近邻点作为种子增长算法新的种子点；(8)重复执行步骤(6)和(7)直至不再有新的种子点,则完成一块联通区域的分割；(9)将识别所得平坦联通区域点集并入集合u中，对g'＝g-u中样点重复执行步骤(4)-(8)，直至中样点的最小曲率值大于设定的曲率阈值为止，则完成所有联通区域分割。
为实现发明目的，所述实物表面样点棱边特征曲面拓扑逼近重建方法，其特征在于步骤(2)中所述非特征区域点集扩展方法，具体为：设x为分块点集gi边界点，xi是x的近邻点集；非特征区域点集扩展步骤为：(1)获取边界点x邻域点集λ(x)；(2)求解邻域参考均值点t；(3)对邻域点集λ(x)进行增益优化；优化过程具体为：沿公式
反向偏移λ(x)的均值点t(x)和t的欧氏距离至tx点，重新搜索tx近邻点加入λ(x)并合并为λ'(x)，若d(t,t(λ(x)))大于d(t,t(λ'(x)))，说明λ'(x)更能反映原始曲面形状；将将偏移、搜索和合并的过程迭代计算,直至d(t,t(λ'(x)))不小于d(t,t(λ(x)))，则邻域点集增益优化结束；(4)将λ'(x)均值点向x所在平面投影得到点r(x)；(5)以x与r(x)之差为方向向量对投影点偏移得到扩展点；(6)对gi点集中所有边界点执行上述步骤完成所有非特征区域点集扩展；上述过程步骤(2)中均值参考点的求解步骤为：①查询边界点集中目标边界点近邻点集，求解其三个坐标值方差，以三个坐标轴中方差最大的轴坐标大小将目标点和近邻点集进行排序；②用三次bezier曲线插值法求解排序点集构成的曲线；③求解目标样点、实际均值点和邻域样本中距目标样点最远点构成平面的法线s；④计算过目标样点且同时垂直于其在边界曲线处的切线和s的直线方程f；⑤计算直线f上距目标点r/2的点，以距离t(x)较近的点为均值参考点。
本发明与现有技术相比，具有以下优点：
(1)对边界点邻域进行增益优化，可一定程度上弥补采样数据不均与导致的拓扑邻域信息缺失，对非均匀数据有较好的适应性，有效避免非均匀采样数据中因边界邻域缺失而导致的曲面重建错误。
(2)可有效避免直接对曲面棱边特征进行重建出现孔洞等错误，提高曲面棱边重建的形位精度。
附图说明
图1是本发明棱边特征曲面拓扑逼近重建方法的程序实现流程图；
图2是边界样点邻域缺失和增益优化后偏移结果示意图；
图3是实施例一中fandisk采样数据在点云分割、扩展和重建步骤效果图；
图4是实施例一中fandisk采样数据最终重建结果的整体和局部棱边特征效果图。
具体实施方式
下面结合附图及实施例对本发明作进一步说明。
图1是本发明棱边特征曲面拓扑逼近重建方法程序实现流程图，通过避免直接对棱边特征重建，基于点云平坦区域联通性对点云进行分割，以边界点邻域点集作为曲面局部样本反映曲面局部特征对点集向外扩展，通过重建扩展后的点云数据，提取各曲面彼此相交的特征线，将交线段端点坐标和三角形内部端点作为细分点进行剖分重建棱边特征。
图2是边界样点邻域缺失和增益优化后偏移结果，将边界点的k-近邻点集作为曲面局部样本能够反映曲面的局部特征，但是k邻域点集的获取主要基于欧氏距离，通常会偏向局部样点分布密集区域，若某一区域同时包含样点分布稀疏区和样点分布密集区，则会导致点云偏移歪斜，最终导致该处重建出现孔洞等错误；对样点邻域进行增益优化，使邻域点集向稀疏区域适度扩展，能够有效避免偏移点歪斜错误，其优化过程为：沿公式
反向偏移λ(x)的均值点t(x)和t的欧氏距离至tx点，重新搜索tx近邻点加入λ(x)并合并为λ'(x)，若d(t,t(λ(x)))大于d(t,t(λ'(x)))，说明λ'(x)更能反映原始曲面形状；将将偏移、搜索和合并的过程迭代计算,直至d(t,t(λ'(x)))不小于d(t,t(λ(x)))，则邻域点集增益优化结束。
图3是fandisk采样数据在本发明中点云分割、扩展和非特征区域重建及细分的结果。如图3所示，通过设定尖锐特征识别阈值，将更多尖锐特征附近点剔除，能够为平坦联通区域分割提供隐形的分割线，从而保证了点云分割的正确性；如[0013]所述，对边界点邻域进行增益优化后，能够适应非均匀采样数据情况，从而提高点云扩展的正确性；为保证重建的曲面足够光滑，在非特征区域曲面重建后对网格模型进行自适应细分，其细分步骤具体为：设三角面片t的法矢为
计算当前三角面片t与其邻域三角面片集合x中各面片法矢点积的方差，将其近似表示当前三角面片的曲率，f(t)的取值范围区间为[0,1]，越接近0则越平坦，否则越陡峭。三角网格细分过程如下：
(1)读取三角网格数据，采用k-d树建立三角网格曲面空间索引结构；(2)计算三角面片法矢点积方差；(3)判断三角形三条边与其公共边所在三角面片曲率值是否大于预设阈值，若至少有一条边大于预设阈值，执行步骤(4)，否则执行步骤(5)；(4)设共边两三角面片为(v1,v2,v3)、(v1,v2,v4)，其公共边(v1,v2)大于设定阈值，则在这条边插入新点ve，计算公式为：
实施例一：对fandisk采样数据应用本文所述方法进行曲面重建。fandisk点云含有棱边特征和自由曲面，并且总体分布非常不均匀。通过观察图4可知，本文方法能够较好重建棱边尖锐特征，对非均匀采样数据有较好的适应性。
以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例。但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。