CN109872394B - 基于最小二乘支持向量机的狭长三角形网格优化方法 - Google Patents
基于最小二乘支持向量机的狭长三角形网格优化方法 Download PDFInfo
- Publication number
- CN109872394B CN109872394B CN201910023996.0A CN201910023996A CN109872394B CN 109872394 B CN109872394 B CN 109872394B CN 201910023996 A CN201910023996 A CN 201910023996A CN 109872394 B CN109872394 B CN 109872394B
- Authority
- CN
- China
- Prior art keywords
- triangular mesh
- vertex
- long
- vector machine
- support vector
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Generation (AREA)
Abstract
本发明公开了基于最小二乘支持向量机的狭长三角形网格优化方法,主要步骤为:1)利用工业CT技术,得到待测对象的CT切片。2)将CT切片重建为三维点云STL模型,并读取三维点云STL模型。3)优化三维点云STL模型。4)对优化后的三维点云STL模型进行狭长三角形网格识别。5)利用最小二乘支持向量机和拉普拉斯算子对狭长三角形网格的顶点坐标进行优化。本发明提出了一种基于支持向量机以最小二乘法作为向量机改进形式的拉普拉斯算子的三角形网格优化方法。本发明可以使狭长三角形网格尽可能的接近正三角形的网格,并使得优化的网格曲面更加光顺。
Description
技术领域
本发明涉及三维重建领域,具体是基于最小二乘支持向量机的狭长三角形网格优化方法。
背景技术
目前从工业CT三维重建到STL三角形网格模型在经过网格简化、优化操作之后,虽能够解决孔洞、面片重叠和多边共线等缺陷,但是网格中往往还存在大量细小或者狭长的三角形网格(如图1)。在计算机图形学中,三角形网格的三角单元越接近正三角形则质量越高,狭长三角形网格质量低而降低渲染效果,影响有限元分析等。因此需要对狭长三角形网格进行优化以提高网格的质量。
发明内容
本发明的目的是解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的,基于最小二乘支持向量机的狭长三角形网格优化方法,主要包括以下步骤:
1)利用工业CT技术,对待测对象进行处理,得到待测对象的CT切片。
2)将CT切片重建为三维点云STL模型,并读取三维点云STL模型。
所述三维点云STL模型为ASCII格式。
将CT切片重建为三维点云STL模型的方法为MC算法。
3)优化三维点云STL模型。
4)对优化后的三维点云STL模型进行狭长三角形网格识别。
对优化后的三维点云STL模型进行狭长三角形网格识别的主要步骤如下:
4.1)建立STL模型的几何元素间的拓扑关系。
4.2)根据三角形面积法识别狭长三角形网格,并找到狭长三角形网格的n阶邻域三角形。n为正整数。
狭长三角形网格面积Q满足下式:
式中,h1、h2和h3分别为三角形网格的三条边长。A为狭长三角网格面积。
4.3)在STL模型中标识出狭长三角形网格。
5)利用最小二乘支持向量机和拉普拉斯算子对狭长三角形网格的顶点坐标进行优化。
对狭长三角形顶点坐标进行优化的主要步骤如下:
5.1)利用最小二乘法拟合狭长三角形网格的一阶邻域三角形曲面,并将一阶邻域三角形曲面作为局部曲面。
最小二乘法拟合函数如下所示:
式中,(x,y)为顶点坐标。a00、a01、a02、a03、a10、a11、a12、a20、a21和a30为计算系数。f(x,y)为顶点坐标函数。
5.2)利用拉普拉斯算子δi分别计算狭长三角形网格的拉普拉斯坐标和一阶邻域三角形顶点的拉普拉斯坐标。
式中,wij为权重。vi和vj为三角形网格的顶点。N(i)为顶点集合。
5.3)利用二阶加权拉普拉斯算子对狭长三角形网格一阶邻域三角形进行优化处理,得到更新后的三角形网格顶点的拉普拉斯坐标。
对狭长三角形网格一阶邻域三角形进行优化处理的主要步骤如下:
5.3.1)计算一阶邻域三角形中每个三角形顶点沿投影梯度方向在切平面上的位移T(v)。
T(v)=-[U2(v)-(U2(v)·Nv)·Nv]。 (4)
式中,U2(v)为二阶加权拉普拉斯算子。Nv为顶点v的法向量。
二阶加权拉普拉斯算子U2(v)如下所示:
式中,wi为权重。i为顶点。U(v)为一阶加权拉普拉斯算子。U(si)为顶点v的一阶邻域顶点。
5.3.2)基于位移T(v),更新各顶点的拉普拉斯坐标。
顶点v更新后的拉普拉斯坐标v′如下所示:
v′=v+τ·T(v)。 (6)
式中,τ为顶点v的调整步长。0<τ≤1。
5.3.3)重复步骤5.3.1和步骤5.3.2,直至T(v)趋近于0。
当T(v)≈0时,所有三角形网格顶点的拉普拉斯坐标更新完成,并将更新后的三角形网格顶点的拉普拉斯坐标作为最小二乘支持向量机的学习样本。
5.4)将更新后的三角形网格顶点的拉普拉斯坐标作为最小二乘支持向量机的学习样本。将拉普拉斯坐标在同一坐标上的分量作为最小二乘支持向量机的学习训练集。
拉普拉斯坐标在同一坐标上的分量S如下所示:
式中,δis为在y坐标上的拉普拉斯坐标分量。i为任意顶点序号。l为顶点个数。Vis为在x坐标上的拉普拉斯坐标分量。R表示整个坐标维度。
5.5)将学习样本输入到最小二乘支持向量机中,并利用学习训练集对最小二乘支持向量机进行训练,得到回归函数式子d(x),即:
式中,αi为拉格朗日乘子。b为偏差。K(x,xi)为核函数。
回归函数式子分解得到三个回归函数fx(x),fy(x),fz(x),即优化的狭长三角形网格顶点坐标函数。
本发明的技术效果是毋庸置疑的。本发明提出了一种基于支持向量机以最小二乘法作为向量机改进形式的拉普拉斯算子的三角形网格优化方法。本发明可以使狭长三角形网格尽可能的接近正三角形的网格,并使得优化的网格曲面更加光顺。
本发明结合了机器学习里的支持向量机采集狭长三角形网格的三个顶点及其n阶邻域三角形网格的顶点作为学习样本训练支持向量机,使其能够表示狭长三角形网格周围局部拓扑结构的表面函数形式拟合作为局部曲面。进而得到最优的拉普拉斯算子相关参数进而优化狭长三角形网格。
附图说明
图1为茶壶模型中存在的狭长三角形;
图2为基于支持向量机的三角网格优化方法流程图;
图3为T(v)与顶点v及其一阶邻域顶点的关系;
图4为网格的局部几何关系;
图5为茶壶优化后的局部狭长三角网格;
图6为化油器优化后的局部狭长三角网格。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图2至图4,基于最小二乘支持向量机(LSSVM)的狭长三角形网格优化方法,主要包括以下步骤:
1)利用工业CT技术,对待测对象进行处理,得到待测对象的CT切片。待测对象可以为工件。
2)将CT切片重建为三维点云STL模型,并读取三维点云STL模型。
所述三维点云STL模型为ASCII格式。若三维点云STL模型不是ASCII格式,则转换为ASCII格式。
将CT切片重建为三维点云STL模型的方法为MC(Marching Cubes)算法。MC算法(移动立方体算法)是一种三维规则数据场等值面生成的经典算法。
3)优化三维点云STL模型。本实施例是针对已经简化的STL模型在一系列优化操作之后,比如孔洞修补等。在优化完成之后三维点云STL模型往往还存在一些狭长的三角面片,本实施例针对这些狭长的三角面片提出优化算法。
4)对优化后的三维点云STL模型进行狭长三角形网格识别。
对优化后的三维点云STL模型进行狭长三角形网格识别的主要步骤如下:
4.1)建立STL模型的几何元素间的拓扑关系,包括STL模型点、线、面片之间的数据结构关系。
4.2)根据三角形面积法识别狭长三角形网格,并找到狭长三角形网格的n阶邻域三角形。n为正整数。本实施例中n=1。
狭长三角形网格面积Q满足下式:
式中,h1、h2和h3分别为三角形网格的三条边长。A为狭长三角网格面积。
4.3)在STL模型中标识出狭长三角形网格,即利用OpenGL对狭长三角形网格进行可视化处理。
OpenGL全称Open Graphics Library,是一个跨平台的第三方图形绘制库。
5)利用最小二乘支持向量机和拉普拉斯算子对狭长三角形网格的顶点坐标进行优化。
对狭长三角形顶点坐标进行优化的主要步骤如下:
5.1)利用最小二乘法拟合狭长三角形网格的一阶邻域三角形曲面,并将一阶邻域三角形曲面作为局部曲面。
最小二乘法拟合函数如下所示:
式中,(x,y)为顶点坐标。a00、a01、a02、a03、a10、a11、a12、a20、a21和a30为计算系数。f(x,y)为顶点坐标函数。上标0、1、2、3表示次方。
5.2)利用拉普拉斯算子δi分别计算狭长三角形网格的拉普拉斯坐标和一阶邻域三角形顶点的拉普拉斯坐标。
式中,wij为权重。vi和vj为三角形网格的顶点。N(i)为顶点集合。
5.3)利用二阶加权拉普拉斯算子对狭长三角形网格一阶邻域三角形进行优化处理,得到更新后的三角形网格顶点的拉普拉斯坐标。
对狭长三角形网格一阶邻域三角形进行优化处理的主要步骤如下:
5.3.1)计算一阶邻域三角形中每个三角形顶点沿投影梯度方向在切平面上的位移T(v)。
T(v)=-[U2(v)-(U2(v)·Nv)·Nv]。 (4)
式中,U2(v)为二阶加权拉普拉斯算子。Nv为顶点v的法向量。
二阶加权拉普拉斯算子U2(v)如下所示:
式中,wi为权重。i为顶点。U(v)为一阶加权拉普拉斯算子。U(si)为顶点v的一阶邻域顶点。
5.3.2)基于位移T(v),更新各顶点的拉普拉斯坐标。
顶点v更新后的拉普拉斯坐标v′如下所示:
v′=v+τ·T(v)。 (6)
式中,τ为顶点v的调整步长。0<τ≤1。
5.3.3)重复步骤5.3.1和步骤5.3.2,直至T(v)趋近于0。
当T(v)≈0时,所有三角形网格顶点的拉普拉斯坐标更新完成,并将更新后的三角形网格顶点的拉普拉斯坐标作为最小二乘支持向量机的学习样本。
趋近于0是指误差小于ζ。本实施例设定ζ=0.001,也可以根据检测对象的实际情况及STL模型进行设置。
5.4)将更新后的三角形网格顶点的拉普拉斯坐标作为最小二乘支持向量机的学习样本。将拉普拉斯坐标在同一坐标上的分量作为最小二乘支持向量机的学习训练集。
拉普拉斯坐标在同一坐标上的分量S如下所示:
式中,δis为在y坐标上的拉普拉斯坐标分量。i为任意顶点序号。l为顶点个数,也即输入学习样本数。Vis为在x坐标上的拉普拉斯坐标分量。R表示整个坐标维度,也即整个坐标平面。
5.5)将学习样本输入到最小二乘支持向量机中,并利用学习训练集对最小二乘支持向量机进行训练,得到回归函数式子d(x),即:
式中,αi为拉格朗日乘子。b为偏差。K(x,xi)为核函数。
回归函数式子分解得到三个回归函数fx(x),fy(x),fz(x),即优化的狭长三角形网格顶点坐标函数。fx(x)为优化的狭长三角形网格顶点x坐标函数,fy(x)为优化的狭长三角形网格顶点y坐标函数,fz(x)为优化的狭长三角形网格顶点z坐标函数。x表示任意顶点;实施例2:
利用基于最小二乘支持向量机的狭长三角形网格优化方法对茶壶狭长三角形网格进行优化的实验,主要包括以下步骤:
1)利用工业CT技术,得到茶壶的CT切片。
2)将CT切片重建为三维点云STL模型,并读取三维点云STL模型。
3)优化三维点云STL模型。
4)对优化后的三维点云STL模型进行狭长三角形网格识别,如图1所示。
5)利用最小二乘支持向量机和拉普拉斯算子对狭长三角形网格的顶点坐标进行优化,如图5所示。
实施例3:
利用基于最小二乘支持向量机的狭长三角形网格优化方法对化油器狭长三角形网格进行优化的实验,主要包括以下步骤:
1)利用工业CT技术,得到化油器的CT切片。
2)将CT切片重建为三维点云STL模型,并读取三维点云STL模型。
3)优化三维点云STL模型。
4)对优化后的三维点云STL模型进行狭长三角形网格识别,如图1所示。
5)利用最小二乘支持向量机和拉普拉斯算子对狭长三角形网格的顶点坐标进行优化,如图6所示。
Claims (5)
1.基于最小二乘支持向量机的狭长三角形网格优化方法,其特征在于,主要包括以下步骤:
1)利用工业CT技术,对待测对象进行处理,得到待测对象的CT切片;
2)将CT切片重建为三维点云STL模型,并读取三维点云STL模型;
3)优化三维STL模型;
4)对优化后的三维点云STL模型进行狭长三角形网格识别;
5)利用最小二乘支持向量机和拉普拉斯算子对狭长三角形网格的顶点坐标进行优化;
对狭长三角形顶点坐标进行优化的主要步骤如下:
5.1)利用最小二乘法拟合狭长三角形网格的一阶邻域三角形曲面,并将一阶邻域三角形曲面作为局部曲面;
最小二乘法拟合函数如下所示:
式中,(x,y)为顶点坐标;a00、a01、a02、a03、a10、a11、a12、a20、a21和a30为计算系数;f(x,y)为顶点坐标函数;
5.2)利用拉普拉斯算子δi分别计算狭长三角形网格的拉普拉斯坐标和一阶邻域三角形顶点的拉普拉斯坐标;
式中,wij为权重;vi和vj为三角形网格的顶点;N(i)为顶点集合;
5.3)利用二阶加权拉普拉斯算子对狭长三角形网格一阶邻域三角形进行优化处理,得到更新后的三角形网格顶点的拉普拉斯坐标;
5.4)将更新后的三角形网格顶点的拉普拉斯坐标作为最小二乘支持向量机的学习样本;将拉普拉斯坐标在同一坐标上的分量作为最小二乘支持向量机的学习训练集;
拉普拉斯坐标在同一坐标上的分量S如下所示:
式中,δis为在y坐标上的拉普拉斯坐标分量;i为任意顶点序号;l为顶点个数;Vis为在x坐标上的拉普拉斯坐标分量;R表示整个坐标维度;
5.5)将学习样本输入到最小二乘支持向量机中,并利用学习训练集对最小二乘支持向量机进行训练,得到回归函数式子d(x),即:
式中,αi为拉格朗日乘子;b为偏差;K(x,xi)为核函数;
回归函数式子分解得到三个回归函数fx(x),fy(x),fz(x),即优化的狭长三角形网格顶点坐标函数。
2.根据权利要求1所述的基于最小二乘支持向量机的狭长三角形网格优化方法,其特征在于,所述三维点云STL模型为ASCII格式。
3.根据权利要求1或2所述的基于最小二乘支持向量机的狭长三角形网格优化方法,其特征在于,将CT切片重建为三维点云STL模型的方法为MC算法。
5.根据权利要求4所述的基于最小二乘支持向量机的狭长三角形网格优化方法,其特征在于,对狭长三角形网格一阶邻域三角形进行优化处理的主要步骤如下:
I)计算一阶邻域三角形中每个三角形顶点沿投影梯度方向在切平面上的位移T(v);
T(v)=-[U2(v)-(U2(v)·Nv)·Nv]; (5)
式中,U2(v)为二阶加权拉普拉斯算子;Nv为顶点v的法向量;
二阶加权拉普拉斯算子U2(v)如下所示:
式中,wi为权重;i为顶点;U(v)为一阶加权拉普拉斯算子;U(si)为顶点v的一阶邻域顶点;
II)基于位移T(v),更新各顶点的拉普拉斯坐标;
顶点v更新后的拉普拉斯坐标v′如下所示:
v′=v+τ·T(v); (7)
式中,τ为顶点v的调整步长;0<τ≤1;
III)重复步骤I)和步骤II),直至T(v)趋近于0;
当T(v)≈0时,所有三角形网格顶点的拉普拉斯坐标更新完成,并将更新后的三角形网格顶点的拉普拉斯坐标作为最小二乘支持向量机的学习样本。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910023996.0A CN109872394B (zh) | 2019-01-10 | 2019-01-10 | 基于最小二乘支持向量机的狭长三角形网格优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910023996.0A CN109872394B (zh) | 2019-01-10 | 2019-01-10 | 基于最小二乘支持向量机的狭长三角形网格优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109872394A CN109872394A (zh) | 2019-06-11 |
CN109872394B true CN109872394B (zh) | 2021-04-27 |
Family
ID=66917662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910023996.0A Expired - Fee Related CN109872394B (zh) | 2019-01-10 | 2019-01-10 | 基于最小二乘支持向量机的狭长三角形网格优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109872394B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110279476B (zh) * | 2019-07-02 | 2021-01-22 | 中国人民解放军国防科技大学 | 基于支持向量机的虚拟牙龈变形处理制作矫治器的方法 |
CN110796693B (zh) * | 2019-09-11 | 2023-03-21 | 重庆大学 | 一种工业ct切片图像直接生成二维有限元模型的方法 |
CN111291793B (zh) * | 2020-01-20 | 2023-11-14 | 北京大学口腔医学院 | 一种网格曲面的元素分类方法、装置及存储介质 |
CN113221403B (zh) * | 2021-04-28 | 2022-10-14 | 中汽数据(天津)有限公司 | 基于深度强化学习的有限元网格优化方法、设备和介质 |
CN113689566B (zh) * | 2021-07-16 | 2022-05-27 | 淮阴工学院 | 一种基于特征约束的三维模型三角网格优化方法 |
CN113792459B (zh) * | 2021-09-09 | 2023-11-14 | 中国航天科工集团第二研究院 | 一种有限元三角形网格的网格光顺优化方法及装置 |
CN117456110B (zh) * | 2023-12-25 | 2024-03-19 | 北京飞渡科技股份有限公司 | 一种三维场景变态网格数据轻量化处理方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106815481A (zh) * | 2017-01-19 | 2017-06-09 | 中国科学院深圳先进技术研究院 | 一种基于影像组学的生存期预测方法及装置 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101655992B (zh) * | 2009-09-15 | 2011-12-21 | 中国科学院软件研究所 | 一种三角网格模型的重建方法 |
EP2422741B1 (en) * | 2010-08-31 | 2013-01-23 | Straumann Holding AG | Computer-implemented method for digitally designing a dental restoration |
CN102930091B (zh) * | 2012-10-24 | 2015-05-20 | 浙江大学 | 一种面向牙齿牙龈的高效网格修补方法 |
CN105405162B (zh) * | 2015-10-16 | 2017-12-05 | 北京师范大学 | 基于局部结构和方向感知的树点云三维重建方法 |
CN105243687B (zh) * | 2015-10-23 | 2017-11-03 | 佛山市南海区广工大数控装备协同创新研究院 | 一种义齿模型三角网格优化方法 |
CN106157370B (zh) * | 2016-03-03 | 2019-04-02 | 重庆大学 | 一种基于粒子群算法的三角网格规范化方法 |
CN106683186A (zh) * | 2016-11-16 | 2017-05-17 | 浙江工业大学 | 一种曲面细节保持的三维模型修复方法 |
CN108230452B (zh) * | 2017-12-18 | 2021-06-29 | 华侨大学 | 一种基于纹理合成的模型补洞方法 |
-
2019
- 2019-01-10 CN CN201910023996.0A patent/CN109872394B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106815481A (zh) * | 2017-01-19 | 2017-06-09 | 中国科学院深圳先进技术研究院 | 一种基于影像组学的生存期预测方法及装置 |
Non-Patent Citations (1)
Title |
---|
基于动态误差控制和PSO的三角网格模型简化优化方法;王武礼 等;《计算机集成制造系统》;20180131;第24卷(第1期);第63-71页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109872394A (zh) | 2019-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109872394B (zh) | 基于最小二乘支持向量机的狭长三角形网格优化方法 | |
CN107767457B (zh) | 一种基于点云快速重建的stl数模生成方法 | |
Song et al. | A progressive point cloud simplification algorithm with preserved sharp edge data | |
US8175734B2 (en) | Methods and system for enabling printing three-dimensional object models | |
CN107610230B (zh) | 一种3d打印数字模型剖切及接口设计方法及应用 | |
CN107358636B (zh) | 一种基于纹理合成的疏松缺陷图像生成方法 | |
KR101555426B1 (ko) | 지형 렌더링 방법 및 장치 | |
CN107610061B (zh) | 一种基于二维投影的保边点云孔洞修补方法 | |
CN107622530B (zh) | 一种高效鲁棒的三角网切割方法 | |
CN102306397A (zh) | 点云数据网格化的方法 | |
Qiu et al. | Direct slicing of cloud data with guaranteed topology for rapid prototyping | |
JPWO2007083602A1 (ja) | 補間処理方法、補間処理装置、形状評価方法、および形状評価装置 | |
JP2014500484A (ja) | 連続b‐スプライン変形を利用してタイヤ表面の三次元画像を前処理する方法 | |
CN111932669A (zh) | 一种基于边坡岩体特征对象的变形监测方法 | |
CN106384386A (zh) | Lod模型生成中网格处理方法、系统以及3d重建方法和系统 | |
CN108230452B (zh) | 一种基于纹理合成的模型补洞方法 | |
CN107610228A (zh) | 基于海量点云的曲面增量拓扑重建方法 | |
CN114494641B (zh) | 一种三维模型轻量化方法及装置 | |
Wang | Computing length-preserved free boundary for quasi-developable mesh segmentation | |
CN108010114B (zh) | 基本图元点云曲面的几何形状识别方法以及特征识别方法 | |
JP2007102595A (ja) | 解析メッシュ生成装置 | |
CN112687010B (zh) | 一种端框钻模的数字化计量方法 | |
JP2007156604A (ja) | 曲面作成方法及び曲面作成プログラム並びに3次元形状処理装置 | |
JP2006277712A (ja) | 3次元モデルのフィッティング装置、方法及びプログラム | |
Azariadis et al. | Product design using point-cloud surfaces: A recursive subdivision technique for point parameterization |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210427 Termination date: 20220110 |
|
CF01 | Termination of patent right due to non-payment of annual fee |