CN1877640A - 一种基于三角形插值曲面细分的几何数据细分方法 - Google Patents

一种基于三角形插值曲面细分的几何数据细分方法 Download PDF

Info

Publication number
CN1877640A
CN1877640A CN 200610036349 CN200610036349A CN1877640A CN 1877640 A CN1877640 A CN 1877640A CN 200610036349 CN200610036349 CN 200610036349 CN 200610036349 A CN200610036349 A CN 200610036349A CN 1877640 A CN1877640 A CN 1877640A
Authority
CN
China
Prior art keywords
triangle
summits
summit
limit
subdivision
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
Application number
CN 200610036349
Other languages
English (en)
Inventor
罗笑南
凌若天
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sun Yat Sen University
National Sun Yat Sen University
Original Assignee
National Sun Yat Sen University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National Sun Yat Sen University filed Critical National Sun Yat Sen University
Priority to CN 200610036349 priority Critical patent/CN1877640A/zh
Publication of CN1877640A publication Critical patent/CN1877640A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Image Generation (AREA)

Abstract

本发明提供一种基于三角形插值曲面细分的几何数据细分方法,其特征在于先利用拓扑规则将三角形的拓扑网格进行分裂,以增加新顶点并形成新的网格拓扑,再利用几何规则计算所有顶点的儿何位置;循环迭代上述步骤至满足循环停止条件。本发明通过在三角形拓扑结构实体上的细分规则,网格中的点将被插值,并进一步细分,逐渐形成稠密的体数据点集。

Description

一种基于三角形插值曲面细分的几何数据细分方法
                          技术领域
本发明属于几何体造型技术领域,特别是涉及一种基于三角形插值曲面细分的几何数据细分方法。
                          技术背景
细分方法是一个从离散到离散的过程,从初始多面体或多边形网格出发,递归地调用细分规则加密控制网格,最终在极限意义下,网格序列收敛到连续甚至光滑曲线、曲面,具有许多参数方法和隐式方法表示的优点,适合计算机的离散表示。2002年,Bajaj等人在论文“A smooth subdivision scheme forhexahedral meshes.The Journal of Visual Computer,2002,18:343-356”中利用三线性插值的技术和基于六面体网格的平均算法,提出了类似的体细分方法,在基于数值实验分析的基础上,分析了其极限性质,这种方法可以用于产生非流形的网格上的几何形体的生成,本质上讲是一种逼近型的体细分方法。
近年来使用细分方法迭代构造三维几何实体作为一个新的研究方向,在理论和应用中逐渐得到了重视和发展,成为一种新的体造型的手段。曲面细分方法目的在于建立了一种基于离散数据迭代细化的方法快速建立三维数据体的技术,表现出了很多新的特性和优点。其中,三角形拓扑网格是形成三维几何形体和进行有限元等计算中常用的三维形体表示的拓扑结构。
目前与基于曲面细分方法的几何造型方法相关的技术方案有:1996年,D.Zorin,P.Schrder,W.Sweldens在论文“Interpolating Subdivision for Meshes withArbitrary Topology.SIGGRAPH 96 Proceedings(1996):189-192”中提出在三角形拓扑网格上的插值型曲面细分方法。这种方法生成的曲面往往绕动过多,曲面光滑性较差,在模型的尖点处的效果往往达不到要求。针对上述方法生成曲面绕动过多、光滑行不足的这些缺陷,本发明提出了一种基于三角形网格拓扑结构上的插值型几何数据细分方法,可以运用于计算机动画、计算机可视化及三维体造型方面。
                          发明内容
本发明的目的在于克服现有技术的不足,提供一种基于三角形插值曲面细分的几何数据细分方法。
为了实现发明目的,采用的技术方案为:
一种基于三角形插值曲面细分的几何数据细分方法,先利用拓扑规则将三角形的拓扑网格进行分裂,以增加新顶点并形成新的网格拓扑,再利用几何规则计算所有顶点的几何位置;循环迭代上述步骤至满足循环停止条件。
上述技术方案中,所述拓扑规则为:
(1)在拓扑网格的每一个三角形中,对于三角形的三条边,各产生一个新的边点。原三角形每条边应包含两个顶点,这里称为A、B,现于每条边都生成了一个新的边点C,分别连接边点C与A,C与B,新生成两条边AC、BC,用新生成的两条边AC、BC代替原有边AB。在一个三角形中,生成了3个新边点C1、C2、C3,三个边点分别两两连接,得到C1C2,C1C3、C2C3三条新的边,把三条新生成的边加入网格拓扑中。
(2)原有的拓扑网格顶点保持不变,所述(1)所生成的每一个新的边点被连接到原有边的两个顶点上。
所述几何规则包括不规则三角形几何规则和规则三角形几何规则,所述不规则三角形几何规则具体如下:
(I)计算边点,以不规则顶点V0为中心,得到所有与不规则顶点相连接的顶点构成一个顶点环V1、V2、V3、...、Vn,其中W1、W2、W3、...、Wn为与顶点对应的权重,那么新顶点的计算方法为:
q = w 0 p 0 k + Σ i = 1 n w i p i k
权重的计算法方法:
w 0 = 3 4 w 1 = 1 2 n + 1 4 w k = 1 2 n cos 2 πk n , ( k = 2,3 , . . . , n )
(II)原顶点位置保持不变。
所述规则三角形几何规则,具体如下:
(I)计算边点:每条边都包含两个顶点,取出所有与该两个顶点相邻接的所有顶点,在规则情况下,这些邻接的顶点为8个,加上边的原有的2个顶点,共有10个顶点,把它们标记为P1、P2、P3、...、P10。其中P1、P2为边的两个顶点,P3、P4为与P1、P2都相邻的两个顶点。那么新顶点的计算方法为:
q = 17 32 ( P 1 + P 2 ) + 1 16 ( P 3 + P 4 ) - 1 32 ( P 5 + P 6 + P 7 + P 8 + P 9 + P 10 )
(II)原顶点位置保持不变。
上述规则三角形几何规则所涉及的拓扑意义上的规则三角形的几何约束为:
(1)所有的拓扑网格都很好连接在一起,也就是说没有孤立点或者孤立边,每一个点都在边上,每一条边都在面上;
(2)每一个面元素都是封闭体,由三角形包围而成,也就是说三角形包围的面不会形成漏缝或者洞,面不会相交和自相交;
(3)每条边都为2个三角形所共享。
(4)对于规则三角形拓扑网格,一个点仅与6条边想连。
所述循环停止条件为三角形已经达到设定的足够的计算精度或显示精度、或三角形的最长的一条边的长度小于设定的值、或三角形最大的元素小于设定的参数值、或三角形的数量大于设定的参数值。
本发明方案定义的是在三角形拓扑结构实体上的细分规则,在面上要求有分布较均匀的拓扑网格点,这些网格点将被插值,并进一步细分,逐渐形成稠密的体数据点集。
其优越性如下:
(1)对于三角形的三维网格数据,可以通过本发明,经过迭代计算生成稠密的连续面数据;
(2)插值型方法所生成的图形经过其初始的控制顶点,当需要进行几何形体变形时,可以直接操纵控制顶点,进行交互式的操作,易于构造任意形状的三维几何形体。
                          附图说明
图1为三角形面片分裂示意图
图2为规则三角形边点细分规则示意图
图3为不规则三角形边点细分规则示意图
                        具体实施方式
下面结合附图对本发明做进一步的说明。
本发明对规则三角形的细分过程如附图1所示,先对该三角形的拓扑网格进行分裂,附图1左图为原三角形,右图为分裂后三角形,具体如下:
(1)在三角形的3条边上,各生成一个边点。
(2)分别连接新生成的3个边点。
(3)按照附图1右图的规则连接新生成的边点与三角形原有的顶点。
(4)从几何数据中删除三角形原有的3条边,用附图1右图中的9条边取而代之。
(5)从几何数据中删除原有三角形,用附图1右图中的4个三角形取而代之。
对网格中的所有三角形进行分裂后,通过几何规则计算点的位置,具体如下:
(1)若三角形为规则三角形,即三角形3个顶点都为规则顶点,与三角形3个顶点相连接的边的数目均为6,那么按照图2给出的掩膜,从三角形网格数据中取出相对应的点,按照以下方法计算新插入点的位置:
q = 17 32 ( P 1 + P 2 ) + 1 16 ( P 3 + P 4 ) - 1 32 ( P 5 + P 6 + P 7 + P 8 + P 9 + P 10 )
(2)若三角形为非规则三角形,即三角形3个顶点中,至少有一个顶点满足此条件:与该顶点相邻接的边的数目不为6。对于边点的生成,取出连接的边数目不为6的顶点P0,然后再取出所有与该顶点相连接的顶点,把相连接的顶点按照逆时针或者顺时针顺序排放好,以边点所在的顶点为P1,其余点按排放顺序编号为P2、P3、P4、....,假设共有n个顶点,那么新边点的位置应该用如下方法确定:
q = w 0 p 0 k + Σ i = 1 n w i p i k
其中Wi为与Pi对应的权重,计算法方法:
w 0 = 3 4 w 1 = 1 2 n + 1 4 w k = 1 2 n cos 2 πk n , ( k = 2,3 , . . . , n )
(3)计算顶点,顶点位置在细分过程中保持不变。
通过以上的三角形插值曲面细分算法,一个三角形将被分成9个小的三角形,在曲面细分循环迭代计算中,可能的循环停止条件如下:已经达到足够的计算精度或显示精度;最长的一条边的长度小于某个预先给定的值;或者最大的面元素都小于某个参数值;面元素的数量大于某个参数值。
实验表明,一般进行3-5次曲面细分就能够得到较稠密的体数据,满足应用要求。
本发明在三角形拓扑结构实体上的细分规则,,网格中的点将被插值,并进一步细分,逐渐形成稠密的体数据点集。

Claims (5)

1、一种基于三角形插值曲面细分的几何数据细分方法,其特征在于先利用拓扑规则将三角形的拓扑网格进行分裂,以增加新顶点并形成新的网格拓扑,再利用几何规则计算所有顶点的几何位置;循环迭代上述步骤至满足循环停止条件。
2、根据权利要求1所述的基于三角形插值曲面细分的几何数据细分方法,其特征在于所述拓扑规则为:
(1)在拓扑网格的每一个三角形中,对于三角形的三条边,各产生一个新的边点。原三角形每条边应包含两个顶点,这里称为A、B,现于每条边都生成了一个新的边点C,分别连接边点C与A,C与B,新生成两条边AC、BC,用新生成的两条边AC、BC代替原有边AB。在一个三角形中,生成了3个新边点C1、C2、C3,三个边点分别两两连接,得到C1C2,C1C3、C2C3三条新的边,把三条新生成的边加入网格拓扑中;
(2)原有的拓扑网格顶点保持不变。
3、根据权利要求1或2所述的基于三角形曲面插值细分的几何数据细分方法,其特征在于所述几何规则为不规则的三角形几何规则,具体如下:
(I)计算边点,以不规则顶点V0为中心,得到所有与不规则顶点相连接的顶点构成一个顶点环V1、V2、V3、...、Vn,其中W1、W2、W3、...、Wn为与顶点对应的权重,那么新顶点的计算方法为:
q = w 0 p 0 k + Σ i = 1 n w i p i k
权重的计算法方法:
w 0 = 3 4 w 1 = 1 2 n + 1 4 w k = 1 2 n cos 2 πk n . ( k = 2,3 , . . . , n )
(II)原顶点位置保持不变。
4、根据权利要求1或2所述的基于三角形插值体细分的几何数据细分方法,其特征在于所述几何规则为规则三角形几何规则,具体如下:
(I)计算边点:每条边都包含两个顶点,取出所有与该两个顶点相邻接的所有顶点,在规则情况下,这些邻接的顶点为8个,加上边的原有的2个顶点,共有10个顶点,把它们标记为P1、P2、P3、...、P10。其中P1、P2为边的两个顶点,P3、P4为与P1、P2都相邻的两个顶点。那么新顶点的计算方法为:
q = 17 32 ( P 1 + P 2 ) + 1 16 ( P 3 + P 4 ) - 1 32 ( P 5 + P 6 + P 7 + P 8 + P 9 + P 10 )
(II)顶点保持不变。
5、根据权利要求1所述的基于三角形插值曲面细分的几何数据细分方法,其特征在于所述循环停止条件为六面体已经达到设定的足够的计算精度或显示精度、或三角形的最长的一条边的长度小于设定的值、或三角形最大的面元素小于设定的参数值、或三角形的面元素的数量大于设定的参数值。
CN 200610036349 2006-07-05 2006-07-05 一种基于三角形插值曲面细分的几何数据细分方法 Pending CN1877640A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200610036349 CN1877640A (zh) 2006-07-05 2006-07-05 一种基于三角形插值曲面细分的几何数据细分方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200610036349 CN1877640A (zh) 2006-07-05 2006-07-05 一种基于三角形插值曲面细分的几何数据细分方法

Publications (1)

Publication Number Publication Date
CN1877640A true CN1877640A (zh) 2006-12-13

Family

ID=37510059

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200610036349 Pending CN1877640A (zh) 2006-07-05 2006-07-05 一种基于三角形插值曲面细分的几何数据细分方法

Country Status (1)

Country Link
CN (1) CN1877640A (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101833778A (zh) * 2010-03-29 2010-09-15 中山大学 一种满足几何约束的细分方法
CN101833748A (zh) * 2010-03-29 2010-09-15 中山大学 一种基于细分方法的图像放大方法
CN101841641A (zh) * 2010-03-29 2010-09-22 中山大学 一种基于细分方法的视频放大方法及系统
CN102598063A (zh) * 2009-09-09 2012-07-18 先进微装置公司 曲面细分引擎及其应用
CN105389851A (zh) * 2015-10-22 2016-03-09 北京航空航天大学 一种基于顶点和新边点统一调整的Loop细分实现方法
CN106504249A (zh) * 2016-10-19 2017-03-15 青岛兰信医学科技有限公司 部分曲面细分算法
CN107194079A (zh) * 2017-05-25 2017-09-22 中国商用飞机有限责任公司 像素细分载荷传递方法和系统
CN107274488A (zh) * 2017-05-26 2017-10-20 天津工业大学 一种三维模型的生成方法
CN103824323B (zh) * 2014-02-26 2018-01-30 中国农业大学 一种基于单幅二维图像的三维重建方法
CN110322547A (zh) * 2019-05-27 2019-10-11 中国石油大学(华东) 一种储层自适应四面体剖分方法

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102598063A (zh) * 2009-09-09 2012-07-18 先进微装置公司 曲面细分引擎及其应用
CN102598063B (zh) * 2009-09-09 2016-06-08 先进微装置公司 曲面细分引擎及其应用
CN101833748A (zh) * 2010-03-29 2010-09-15 中山大学 一种基于细分方法的图像放大方法
CN101841641A (zh) * 2010-03-29 2010-09-22 中山大学 一种基于细分方法的视频放大方法及系统
CN101833748B (zh) * 2010-03-29 2013-01-02 中山大学 一种基于细分方法的图像放大方法
CN101841641B (zh) * 2010-03-29 2013-01-02 中山大学 一种基于细分方法的视频放大方法及系统
CN101833778A (zh) * 2010-03-29 2010-09-15 中山大学 一种满足几何约束的细分方法
CN103824323B (zh) * 2014-02-26 2018-01-30 中国农业大学 一种基于单幅二维图像的三维重建方法
CN105389851B (zh) * 2015-10-22 2018-07-27 北京航空航天大学 一种基于顶点和新边点统一调整的Loop细分实现方法
CN105389851A (zh) * 2015-10-22 2016-03-09 北京航空航天大学 一种基于顶点和新边点统一调整的Loop细分实现方法
CN106504249A (zh) * 2016-10-19 2017-03-15 青岛兰信医学科技有限公司 部分曲面细分算法
CN106504249B (zh) * 2016-10-19 2019-03-19 青岛达芬奇科技有限公司 部分曲面细分方法
CN107194079A (zh) * 2017-05-25 2017-09-22 中国商用飞机有限责任公司 像素细分载荷传递方法和系统
CN107274488A (zh) * 2017-05-26 2017-10-20 天津工业大学 一种三维模型的生成方法
CN107274488B (zh) * 2017-05-26 2021-03-23 天津工业大学 一种三维片上网络的三维模型的生成方法
CN110322547A (zh) * 2019-05-27 2019-10-11 中国石油大学(华东) 一种储层自适应四面体剖分方法
CN110322547B (zh) * 2019-05-27 2023-09-19 中国石油大学(华东) 一种储层自适应四面体剖分方法

Similar Documents

Publication Publication Date Title
CN1877640A (zh) 一种基于三角形插值曲面细分的几何数据细分方法
CN110110413B (zh) 一种基于材料场缩减级数展开的结构拓扑优化方法
CN102306396B (zh) 一种三维实体模型表面有限元网格自动生成方法
CN111709085B (zh) 一种约束阻尼薄板结构拓扑优化设计方法
CN109377561A (zh) 一种基于共形几何的数模表面网格生成方法
CN101110126A (zh) 一种三维模型网格重建方法
CN109159425A (zh) 三维模型的切片方法及三维打印装置
CN107038308B (zh) 一种基于线性内插的规则格网地形建模方法
CN102509339A (zh) 一种带纹理约束的三维模型顶点聚类简化方法
CN103353916B (zh) 基于工程的复合材料层合板铺层优化后处理方法
CN108427861A (zh) 一种基于mpt工具包构建材料周期性多晶结构的方法
Wang et al. Fast mesh simplification method for three-dimensional geometric models with feature-preserving efficiency
CN102831648A (zh) 一种定节点量曲线逼近和曲面网格生成优化技术
CN110188423A (zh) 一种基于有限元网格划分的线性工程结构快速bim建模方法
CN106981097A (zh) 一种基于分区局部光顺权因子的t样条曲面拟合方法
CN1828672A (zh) 一种基于六面体插值体细分的几何数据细分方法
CN101354788A (zh) 三维网格模型的一维化无损几何压缩方法
CN1975787A (zh) 一种基于高阶三角形插值曲面细分的几何数据细分方法
CN114139471A (zh) 适用于cfd的城市街区三维模型构建和网格划分方法及系统
CN104504758A (zh) 义齿冠表面曲面生成方法
CN111047684A (zh) 一种基于三维模型特征的模型简化方法
CN107644139A (zh) 一种从cad模型到cae模型的属性映射方法
CN108062433B (zh) 基于增材再制造点云模型的梯度曲面分层方法
CN113111612A (zh) 一种基于自适应空间剖分的离散点云重复点快速查找方法
CN106067192A (zh) 基于逆插值Loop的网格简化方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication