CN102930594A - 基于保持边界的计算机三维地形网格的绘制方法 - Google Patents
基于保持边界的计算机三维地形网格的绘制方法 Download PDFInfo
- Publication number
- CN102930594A CN102930594A CN2012103531944A CN201210353194A CN102930594A CN 102930594 A CN102930594 A CN 102930594A CN 2012103531944 A CN2012103531944 A CN 2012103531944A CN 201210353194 A CN201210353194 A CN 201210353194A CN 102930594 A CN102930594 A CN 102930594A
- Authority
- CN
- China
- Prior art keywords
- summit
- landform
- drafting
- computerized
- regular grid
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 81
- 238000012423 maintenance Methods 0.000 title claims abstract description 11
- 238000009877 rendering Methods 0.000 title abstract description 9
- 230000001788 irregular Effects 0.000 claims abstract description 30
- 230000000007 visual effect Effects 0.000 claims abstract description 10
- 230000008602 contraction Effects 0.000 claims description 44
- 230000008569 process Effects 0.000 claims description 11
- 230000001143 conditioned effect Effects 0.000 claims description 3
- 238000012952 Resampling Methods 0.000 claims description 2
- 239000011800 void material Substances 0.000 claims 2
- 238000005070 sampling Methods 0.000 abstract description 4
- 238000005192 partition Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 6
- 239000000203 mixture Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 3
- 230000000875 corresponding effect Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000013507 mapping Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000008439 repair process Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000011960 computer-aided design Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000007500 overflow downdraw method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Landscapes
- Processing Or Creating Images (AREA)
- Image Generation (AREA)
Abstract
本发明公开了一种计算机三维地形网格的绘制方法。该方法包括:对数字高程模型进行重采样生成地形规则网格数据;将地形规则网格数据转换为网格简化辅助信息;利用网格简化辅助信息对地形规则网格数据进行保持边界的不规则边折叠简化;将不规则边折叠简化结果与地形规则网格数据按层次结构表示方法组织成具有保持边界特性的地形格式。利用本发明,可以完全避免三维虚拟场景显示中地形分块间的裂缝问题,在运行时无需占用额外计算资源即可获得无缝拼接平滑过渡的视觉效果,同时它采用的不规则边折叠简化方式能够很好地保留地形原始特征,提高远景时的显示质量。
Description
技术领域
本发明涉及计算机三维图形绘制辅助设计领域,尤其涉及一种适合于大范围地形场景实时绘制的基于保持边界的计算机三维地形网格的绘制方法。
背景技术
实时可交互地渲染带有大范围地形的户外场景一直是计算机图形学中的一项重要挑战,广泛应用在模拟飞行、虚拟战场、三维游戏等虚拟现实系统中。但实时地形绘制的工作负荷很大,尤其是涉及范围较广、逼真度较高的数据。为减轻工作负荷,目前主要有两种方法:一种是层次细节(level of details,LoD)技术,另一种是地形数据分块技术。
LoD的基本思想是在不影响视觉效果的前提下,通过简化模型、降低需要绘制的三角面片数量,达到减轻硬件负担、加速绘制的目的。LoD技术主要分面向CPU设计与面向GPU设计这两种实现方式。
面向CPU设计的LoD方法通常是在显示运行时,依靠CPU计算能力对模型的三角面片等底层基本图元进行细化或者简化操作,从而获得在当前视觉条件下最合适的网格模型进行绘制。面向CPU设计的LoD有两种典型的实现算法:一种是Lindstrom在1996年《Real-time,continuous levelof detail rendering of height fields》一文中提出的基于四叉树的地形渲染算法,他将所有顶点数据自底向上地组织成四叉树三角网,根据与视点的远近关系逐层简化网格;另一种是Duchaineau在1997年《ROAMing terrain:real-time optimally adapting meshes》一文中提出的基于三角形二叉树的ROAM算法,并引入三角形分裂和合并两个优先队列来控制算法的执行效率,是应用较广泛的一种算法。但面向CPU的LoD方法由于操作过于底层和精细,需要耗费大量CPU计算资源,而且这些算法更着重对一整块地形进行视相关多分辨率优化,对地形块内的裂缝采用顶点依赖关系进行修补,但没有对多块地形之间的边界裂缝避免进行阐述,这是面向CPU的LoD方法两大缺陷。
面向GPU的LoD算法将简化和细化的操作转移到GPU上,同时操作粒度由基本图元上升到批量图元集合,在运行时对模型进行粗粒度优化,充分发挥了GPU的并行处理和高效计算能力。Boer在《Fast terrainrendering using geometrical mipmapping》一文中提出著名的GeometricalMipmapping方法,将地形块如同图片一样进行分级采样,形成金字塔状LoD结构用作切换,并通过改变边界点的连接方式来填补不同LoD级别地形块间的裂缝。Boer的方法思想简单易行,但Mipmapping是均匀采样,各处简化程度一致,导致粗略模型误差较大,而且实际运行时需要监视所有在显示地形块级别并不断调整边界连接方式,占用CPU资源。Ulrich在2002年《Rendering massive terrains using chunked level of detail control》一文中描述了ChunkLoD的方法,基本思想是对每个地形块进行预处理,根据视相关条件和误差阈值在地形块内预生成不同精度等级的瓦片,显示时根据视觉信息在这些瓦片间切换,另外给这些瓦片加上垂直裙边,这样无论是地形块内还是块间裂缝都可以被遮掩。ChunkLoD是GPU友好的批处理LoD算法,但最大的缺点是遮掩裂缝用的垂直裙边会造成几何以及纹理空间的不连贯,影响视觉效果。Hao在2007年《A method for terrainrendering real-time based on two-level model》一文中采用了高精度规则格网和简化不规则三角网相结合的地形显示方式,而用图像空间融合方法掩饰裂缝,实际上并没有真正填补裂缝问题。从上面看出,面向GPU的LoD方法平衡了CPU与GPU的负载,但无法或不便在图元层面上精确修补裂缝,更多是简单地弥补或掩饰,影响视觉效果。
地形数据分块是应对数据量规模过大,不适合一次性完全加载到内存而采取的策略,基本思路是将大地形数据分成大小适中的块,运行时只加载视点邻近相关范围的地形块用于显示。分块思想不仅减轻了计算机硬件的压力,在实际运行时也因为可以有次序地分批处理显得灵活方便。在上一段分析面向GPU的LoD方法时,地形分块思想已经体现在了诸如Mipmapping、ChunkLoD等方法中。分块的最大问题就是如何处理不同地形块间因为精细程度不同而产生的裂缝,上文中也提及了相应措施,总的来说,都因为无法精确修补而影响视觉效果。
为使LoD技术与地形分块思想更好地结合,克服占用计算资源和精确修补裂缝两者之间的矛盾,本发明提供了一种计算机三维地形网格的绘制方法,使大范围地形场景中的地形显示在保证高效流畅的同时也具有很好的视觉效果。
发明内容
(一)要解决的技术问题
有鉴于此,本发明的主要目的在于提供一种保证大范围地形场景平滑流畅显示,基于保持边界思想的地形网格简化与组织的方法。
(二)技术方案
为达到上述目的,本发明提供了一种计算机三维地形网格的绘制方法,包括:对数字高程模型进行重采样生成地形规则网格数据;将地形规则网格数据转换为网格简化辅助信息;利用网格简化辅助信息对地形规则网格数据进行保持边界的不规则边折叠简化;将不规则边折叠简化结果与地形规则网格数据按层次结构表示方法组织成具有保持边界特性的地形格式。
(三)有益效果
从上述技术方案可以看出,本发明具有以下有益效果:
1、本发明提供的这种计算机三维地形网格的绘制方法,是基于保持边界的思想来避免块间裂缝问题的,以此方法形成的地形无论处于什么精度等级,都能保持边界在几何与纹理空间的一致性,从而在显示时不用耗费额外计算资源去做裂缝修补工作。
2、本发明提供的这种计算机三维地形网格的绘制方法,是基于不规则边折叠简化的思想进行网格简化的,一方面不规则的简化方式可以用更少的顶点构成地形,减少数据量,另一方面选择好的边折叠简化规则可以更大限度地保留地形特征,降低简化误差。
3、本发明提供的这种计算机三维地形网格的绘制方法,是基于层次结构表示方法来组织地形格式的,层次结构表示方法能够支持面向GPU的LoD算法,即方便对图元集合进行批量操作,有利于快速绘制大范围地形场景。
附图说明
图1为依照本发明实施例的计算机三维地形网格的绘制方法流程图;
图2为依照本发明实施例的地形规则网格的行列坐标与顶点经纬度的示意图;
图3为依照本发明实施例的由地形规则网格中一点映射到数字高程模型中对应位置的示意图;
图4为依照本发明实施例的地形规则网格的三维坐标系统的示意图;
图5为依照本发明实施例的地形规则网格的三角面片顶点索引数组的存储方式示意图;
图6为依照本发明实施例的计算边折叠代价时两点间的曲率量几何意义示意图;
图7为依照本发明实施例的按照层次结构表示方法组织的地形格式的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
如图1所示,图1为依照本发明实施例的计算机三维地形网格的绘制方法流程图,主要步骤包括:
步骤101、对数字高程模型进行重采样生成地形规则网格数据;
步骤102、将地形规则网格数据转换为网格简化辅助信息;
步骤103、利用网格简化辅助信息对地形规则网格数据进行保持边界的不规则边折叠简化;
步骤104、将不规则边折叠简化结果与地形规则网格数据按层次结构表示方法组织成具有保持边界特性的地形格式。
步骤101所述的重采样过程具体如下:
第一,确定规则网格顶点的数目。通常规则网格一边的定点数取为2n+1(n=6或7)最为合适,因此顶点总数目VertNum=(2n+1)2,分配顶点三维坐标数组的存储空间为VertPos[VertNum],VertPos[]每一项为Vector3类型,Vector3类型的结构如下
struct Vector3{float x,y,z}.
第二,确定每个顶点的经纬度。如图2所示,假设生成规则网格的左下角经纬度为llLong、llLat,右上角经纬度为ruLong、ruLat,则可以认为网格内任一点v(c,r)(0≤c,r≤2n+1)的经纬度vLong、vLat分别是
stepLong=(ruLong-llLong)/(VertNum-1);
vLong=llLong+c*stepLong;
stepLat=(ruLat-llLat)/(VertNum-1);
vLat=llLat+r*stepLat.
第三,根据每个顶点的经纬度,将顶点映射到数字高程模型中的对应位置。如图3所示,假设数字高程模型有Nr行Nc列数据,左下角经纬度分别为demllLong、demllLat,采样间隔弧度为cellsize。则从规则网格的点v映射到数字高程模型中的点v’,其行列坐标计算方法为
cd=(vLong-demllLong)/cellsize;
rd=(vLat-demllLat)/cellsize.
此时,cd、rd不一定是整数。
第四,找出数字高程模型中与上述位置最近的四个高程数据,使用双线性插值计算得出结果,作为规则网格中对应顶点的高度值。仍如图3所示,与v’最接近的四个高程数据以及它们的坐标分别为
hll(floor(cd),floor(rd));
hlu(floor(cd),floor(rd)+1);
hru(floor(cd)+1,floor(rd)+1);
hrl(floor(cd)+1,floor(rd)).
因此使用双线性插值计算点v的高度y为
hl=hll*(floor(rd)+1-rd)+hlu*(rd-floor(rd));
hr=hrl*(floor(rd)+1-rd)+hru*(rd-floor(rd));
y=hl*(floor(cd)+1-cd)+hr*(cd-floor(cd)).
第五,对地形规则网格的所有顶点进行上述计算,生成地形规则网格数据。地形规则网格数据包括顶点三维坐标数组和顶点索引数组。顶点三维坐标数组VertPos[VertNum]中,每一项的y值已经由双线性插值重采样获得,计算x、z的方法如图4所示,X轴指向地形的正东方向,Z轴指向地形的正南方向,原点O(Oc,Or)位于地形中心点,假设点距为dist,则任一点v(c,r)的三维坐标x,z分量计算如下
Oc=2n-1;
Or=2n-1;
x=(c-Oc)*dist;
z=(Or-r)*dist.
通过上述过程,按照顶点经度由西向东、纬度由南向北的顺序计算并得到顶点三维坐标数组VertPos[VertNum]。
对于三角面片顶点索引数组,如图5所示,图5右侧是地形规则网格三角化之后的状态,三角面片数目TriNum=22n+1,因此需要的索引数目为IndexNum=3*TriNum,设三角面片顶点索引数组的存储空间为TriIndex[IndexNum],TriIndex[]每一项均为int型,记录某一顶点在三维坐标数组VertPos[VertNum]中的序号。
TriIndex[IndexNum]存储顶点索引的顺序,从三角面片层次上看,如图5右侧箭头线所示,是按z字型方式存储的,实线箭头代表了同一行中三角面片的存储顺序,虚线箭头代表一行存储结束后折返到上一行的起始处。对具体三角面片顶点进行存储的顺序,如图5左侧的局部放大图所示,首先存储三角面片1的顶点v(0,0)、v(1,0)、v(0,1),按箭头指示顺序,然后存储三角面片2的顶点v(0,1)、v(1,0)、v(1,1),按箭头指示顺序,其他以此类推。任一顶点v(c,r)在三维坐标数组VertPos[VertNum]中的序号id的推算公式为
id=r*(2n+1)+c.
通过上述计算可完全得到三角面片顶点索引数组TriIndex[IndexNum]。从而最终完成地形规则网格数据的生成。
步骤102,将地形规则网格数据转换为网格简化辅助信息。
网格简化辅助信息包括顶点信息和三角面片信息。顶点信息的具体实现为
三角面片信息具体实现为
分别设两者存储空间为Vertex[VertNum]和Triangle[TriNum]。
步骤102中,由地形规则网格数据转换到顶点信息和三角面片信息的过程如下:
第一,遍历顶点三维坐标数组VertPos[VertNum],填充顶点信息Vertex[VertNum]中的id、position、boundary等内容;
第二,遍历三角面片顶点索引数组TriIndex[IndexNum],填充三角面片信息Triangle[TriNum]中的vertex[3]、valid等内容,并对vertex[3]中的每一项,更新其adjVertices与adjFaces的内容;
第三,计算Triangle[TriNum]中的normal的值,公式如下
normal=(vertex[1]-vertex[0])×(vertex[2]-vertex[1]).
计算出的normal须做归一化处理。
步骤103,利用网格简化辅助信息对地形规则网格数据进行保持边界的不规则边折叠简化。
根据保持边界的原则,对顶点信息Vertex[VertNum]中的cost、collapse等内容进行初始化。对于边界点,cost直接置为无穷大,collapse置为NULL;对非边界点v,cost计算方法如下:
第一,计算点v折叠到与其连接的某一点u的代价cost(v,u)
其中,C为v、u两点之间曲率,Tv为包含点v的三角面片集合,Tv,u为包含点v、u的三角面片集合,nf代表三角面片f的法向,ng同理。C的几何意义如图6所示,三角面片1、2、3、4均包含点v,法向为n1、n2、n3、n4,三角面片1、2均包含点v、u,对每个包含点v的三角面片,计算其法向与包含点v、u的三角面片法向偏离值之和,再从中选择最大值作为曲率量C。
对顶点v中的相连接的顶点列表adjVertices的每一项,都按上述方法计算。
第二,从上面计算结果中选择最小值,赋给顶点v的cost变量,取最小值时所折叠到的目标顶点,赋给collapse。
初始化完成后,则可以开始进行不规则边折叠简化,过程为:
第一,从顶点信息Vertex[VertNum]中选择cost最小的那一项,设为v;
第二,设v.collapse所指向的顶点为u,则将v.adjVertices和v.adjFaces更新到u.adjVertices和u.adjFaces,对u.adjFaces列表中同时含有顶点v、u的三角面片将其valid值置为无效(false),将v.cost更新为无穷大,重新计算顶点u的cost和collapse。
第三,重复上述两步直到满足一定条件,如简化顶点数目达到设定阈值时,简化过程结束,此时简化结果由三角面片信息Triangle[TriNum]中valid变量值为有效(true)的所有三角面片构成。
步骤104,将不规则边折叠简化结果与地形规则网格数据按层次结构表示方法组织成具有保持边界特性的地形格式。
根据简化预设条件不同,不规则边折叠简化结果可以有多个,这需要相应数量的Triangle[TriNum]拷贝,每一份经过步骤103处理过的拷贝都是一个简化结果。
将不规则边折叠简化结果与地形规则网格数据按层次结构表示方法组织,如图7所示,分为两大层次,上层精细层和下层粗略层。精细层由地形规则网格数据构造。粗略层则由不规则边折叠简化结果来构造,并根据简化结果的数量细分成相应数量的若干子层,每种简化结果是一个子层,且简化结果保留的顶点数目越少,子层的级别越低。在使用地形过程中,选择哪一层次的数据作为输出显示,是根据视觉信息来决定的,通常使用视点距离信息,距离越近,使用层次级别越高的数据。由于不规则边折叠简化结果是在保持边界的前提下完成的,如图7所示,因此任何层次的地形数据在边界处都具有相同的顶点,保持很好的一致性。
按层次结构表示方法组织的地形格式,其文件记录存储的具体方式为:
第一,存储所有顶点的三维坐标数据,使用单精度浮点数存储,数据来自地形规则网格数据中的顶点三维坐标数组VertPos[VertNum];
第二,存储地形规则网格的所有三角面片,每个三角面片用三个顶点索引表示,索引以整型存储,数据来自地形规则网格数据中的三角面片顶点索引数组TriIndex[IndexNum]。
第三,存储不规则边折叠简化结果,按照保留顶点数目由多到少的顺序依次存储,对每个简化结果,存储其保留的三角面片序列,每个三角面片同样用三个顶点索引表示并以整型存储,数据来自相应的Triangle[TriNum]拷贝中有效的三角面片信息。
综上所述,本发明通过对数字高程模型进行重采样,生成地形规则网格数据;然后将地形规则网格数据转换为网格简化辅助信息;接着利用网格简化辅助信息,对地形规则网格数据进行保持边界的不规则边折叠简化;最后将不规则边折叠简化结果与地形规则网格数据按层次结构表示方法,便能组织成具有保持边界特性的地形格式,整个过程完成。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (26)
1.一种计算机三维地形网格的绘制方法,其特征在于,包括:
对数字高程模型进行重采样生成地形规则网格数据;
将地形规则网格数据转换为网格简化辅助信息;
利用网格简化辅助信息对地形规则网格数据进行保持边界的不规则边折叠简化;
将不规则边折叠简化结果与地形规则网格数据按层次结构表示方法组织成具有保持边界特性的地形格式。
2.根据权利要求1所述的计算机三维地形网格的绘制方法,其特征在于,所述对数字高程模型进行重采样生成地形规则网格数据,是采用双线性插值法进行重采样计算的。
3.根据权利要求2所述的计算机三维地形网格的绘制方法,其特征在于,其中双线性插值法重采样计算的具体步骤包括:
确定规则网格顶点的数目;
确定规则网格每个顶点的经纬度;
根据每个顶点的经纬度,将顶点映射到数字高程模型中的对应位置;
找出数字高程模型中与上述位置最近的四个高程数据,使用双线性插值计算得出结果,作为规则网格中对应顶点的高度值;
对地形规则网格的所有顶点进行上述计算,生成地形规则网格数据。
4.根据权利要求1所述的计算机三维地形网格的绘制方法,其特征在于,所述对数字高程模型进行重采样生成地形规则网格数据,其中地形规则网格数据是采用顶点三维坐标数组和三角面片顶点索引数组来表达的。
5.根据权利要求4所述的计算机三维地形网格的绘制方法,其特征在于,其中顶点三维坐标数组的每一项记录了某个顶点在三维空间中的坐标信息,三维空间采用直角坐标系,原点水平位置在地形网格中央,高度为0,X轴指向正东,Z轴指向正南,Y轴垂直向上。
6.根据权利要求4所述的计算机三维地形网格的绘制方法,其特征在于,其中顶点三维坐标数组的存储顺序,是按照顶点的经度由西向东、纬度由南向北的顺序依次记录坐标信息的。
7.根据权利要求4所述的计算机三维地形网格的绘制方法,其特征在于,其中三角面片顶点索引数组的每一项,记录三角面片某个顶点在顶点三维坐标数组中的序号。
8.根据权利要求4所述的计算机三维地形网格的绘制方法,其特征在于,其中的三角面片顶点索引数组,从起始处开始每三项隐式地表示一个三角面片,整个三角面片顶点索引数组隐式地表达了一串三角面片序列,这串三角面片序列是按照经度由西向东、纬度由南向北的顺序排列的。
9.根据权利要求1所述的计算机三维地形网格的绘制方法,其特征在于,所述将地形规则网格数据转换为网格简化辅助信息,其中网格简化辅助信息包括顶点信息和三角面片信息。
10.根据权利要求9所述的计算机三维地形网格的绘制方法,其特征在于,其中顶点信息具体内容包括:
该顶点在三维坐标数组中的序号;
该顶点的三维坐标;
该顶点是否是边界点;
与该顶点相连接的顶点列表;
与该顶点相连接的三角面片列表;
该顶点的边折叠代价;
该顶点的边折叠目标顶点。
11.根据权利要求9所述的计算机三维地形网格的绘制方法,其特征在于,其中三角面片信息具体内容包括:
该三角面片的顶点;
该三角面片的法向;
该三角面片是否有效。
12.根据权利要求9所述的计算机三维地形网格的绘制方法,其特征在于,其中网格简化辅助信息的具体转换步骤包括:
遍历地形规则网格数据中的顶点三维坐标数组,据此更新顶点信息中的序号和三维坐标内容;
遍历地形规则网格数据中的三角面片顶点索引数组,据此更新三角面片信息中的顶点内容,并对各顶点更新其包含的相连接顶点列表和相连接三角面片列表内容,最后计算三角面片法向。
13.根据权利要求1所述的计算机三维地形网格的绘制方法,其特征在于,所述利用网格简化辅助信息对地形规则网格数据进行保持边界的不规则边折叠简化,包括以下步骤:
确定保持边界的原则;
初始化所有顶点信息中的边折叠代价和边折叠目标顶点;
获取当前具有最小边折叠代价的顶点,并将其折叠到目标顶点,同时更新目标顶点的相连接顶点列表、相连接三角面片列表、边折叠代价、边折叠目标顶点;
重复上一步直到需简化的顶点数目满足预设条件。
14.根据权利要求13所述的计算机三维地形网格的绘制方法,其特征在于,其中确定保持边界的原则,是指位于网格边界上的顶点其边折叠代价始终为无穷大,边折叠目标顶点始终为空,且当计算其他顶点到边界顶点的边折叠代价时,代价值直接置为无穷大。
15.根据权利要求13所述的计算机三维地形网格的绘制方法,其特征在于,其中顶点信息中的边折叠代价,是从与该顶点相连接的顶点列表中,选取计算出的多个边折叠代价里最小的值。
16.根据权利要求14所述的计算机三维地形网格的绘制方法,其特征在于,其中的边折叠代价,是用待折叠顶点与预选目标顶点之间的距离,乘以它们之间的曲率,所得结果再进行平方来计算的,但若两个顶点中有一个是边界点,边折叠代价直接置为无穷大。
17.根据权利要求13所述的计算机三维地形网格的绘制方法,其特征在于,其中顶点信息中的边折叠目标顶点,是从与该顶点相连接的顶点列表中,选取具有最小边折叠代价值的那个顶点。
18.根据权利要求13所述的计算机三维地形网格的绘制方法,其特征在于,其中更新目标顶点的相连接顶点列表、相连接三角面片列表,是将已折叠顶点的相连接顶点列表、三角面片列表添加到目标顶点的相连接顶点列表和三角面片列表中,并去除无效项,无效项包括重复的顶点或三角面片、不再构成三角形状的面片。
19.根据权利要求13所述的计算机三维地形网格的绘制方法,其特征在于,其中更新目标顶点的边折叠代价、边折叠目标顶点,是在已更新目标顶点的相连接顶点列表后,重新计算选取。
20.根据权利要求1所述的计算机三维地形网格的绘制方法,其特征在于,所述将不规则边折叠简化结果与地形规则网格数据按层次结构表示方法组织成具有保持边界特性的地形格式,其中不规则边折叠简化结果可以是多个满足不同简化预设条件的结果。
21.根据权利要求1所述的计算机三维地形网格的绘制方法,其特征在于,所述将不规则边折叠简化结果与地形规则网格数据按层次结构表示方法组织成具有保持边界特性的地形格式,其中层次结构表示方法是指将地形分为上层精细层和下层粗略层这两个主要层次来表示。
22.根据权利要求21所述的计算机三维地形网格的绘制方法,其特征在于,其中的上层精细层,是由地形规则网格数据构造。
23.根据权利要求21所述的计算机三维地形网格的绘制方法,其特征在于,其中的下层粗略层,是由不规则边折叠简化结果来构造,粗略层根据简化结果的数量细分成相应数量的若干子层,每种简化结果是一个子层,并且简化结果保留的顶点数目越少,子层的级别越低。
24.根据权利要求21所述的计算机三维地形网格的绘制方法,其特征在于,其中的层次结构表示方法,在选择使用哪一层次的数据来表示地形时,是根据视觉信息来决定的,通常使用视点距离信息,距离越近,使用层次级别越高的数据。
25.根据权利要求1所述的计算机三维地形网格的绘制方法,其特征在于,所述将不规则边折叠简化结果与地形规则网格数据按层次结构表示方法组织成具有保持边界特性的地形格式,其中的保持边界特性是指,由于不规则边折叠简化结果是在保持边界的前提下完成的,因此在边界处与地形规则网格具有相同的顶点。
26.根据权利要求1所述的计算机三维地形网格的绘制方法,其特征在于,所述将不规则边折叠简化结果与地形规则网格数据按层次结构表示方法组织成具有保持边界特性的地形格式,其中的地形格式具体记录存储方式为:首先使用单精度浮点数存储所有顶点的三维坐标数据;其次存储地形规则网格的所有三角面片,三角面片用整型存储的顶点索引值表示;最后按照保留顶点数目由多到少的顺序存储不规则边折叠简化结果,对每个简化结果同样存储其保留的三角面片序列的顶点索引值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210353194.4A CN102930594B (zh) | 2012-09-20 | 2012-09-20 | 基于保持边界的计算机三维地形网格的绘制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210353194.4A CN102930594B (zh) | 2012-09-20 | 2012-09-20 | 基于保持边界的计算机三维地形网格的绘制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102930594A true CN102930594A (zh) | 2013-02-13 |
CN102930594B CN102930594B (zh) | 2015-06-17 |
Family
ID=47645384
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210353194.4A Active CN102930594B (zh) | 2012-09-20 | 2012-09-20 | 基于保持边界的计算机三维地形网格的绘制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102930594B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104680573A (zh) * | 2015-01-20 | 2015-06-03 | 西安科技大学 | 一种基于三角网格简化的纹理映射方法 |
CN105303610A (zh) * | 2015-12-05 | 2016-02-03 | 中国航空工业集团公司洛阳电光设备研究所 | 一种用于机载地形三维图像绘制的地形数据处理方法 |
CN106153030A (zh) * | 2015-03-25 | 2016-11-23 | 高德软件有限公司 | 一种导航电子地图数据存储方法及装置 |
CN106777816A (zh) * | 2017-01-18 | 2017-05-31 | 上海核工程研究设计院 | 一种场景动态模型管理系统 |
CN107689189A (zh) * | 2017-08-31 | 2018-02-13 | 城市生活(北京)资讯有限公司 | 一种三维电子地图的绘制方法及装置 |
CN108053475A (zh) * | 2018-02-09 | 2018-05-18 | 城市生活(北京)资讯有限公司 | 一种三维地形显示方法及装置 |
CN108986212A (zh) * | 2018-06-21 | 2018-12-11 | 东南大学 | 一种基于裂缝消除的三维虚拟地形lod模型的生成方法 |
CN109389678A (zh) * | 2017-08-08 | 2019-02-26 | 龙芯中科技术有限公司 | 确定操作焦点坐标的方法、装置及设备 |
CN109992640A (zh) * | 2019-04-11 | 2019-07-09 | 北京百度网讯科技有限公司 | 位置网格的确定方法及装置、设备及存储介质 |
CN110555085A (zh) * | 2018-03-29 | 2019-12-10 | 中国石油化工股份有限公司 | 一种三维模型加载方法与装置 |
CN110633282A (zh) * | 2019-09-18 | 2019-12-31 | 四川九洲空管科技有限责任公司 | 一种空域资源多级三维网格化方法及工具 |
CN111105492A (zh) * | 2019-12-18 | 2020-05-05 | 四川大学 | 一种地形平滑过程中瓦片高程数据快速调度方法 |
WO2022257594A1 (zh) * | 2021-06-10 | 2022-12-15 | 广东奥普特科技股份有限公司 | 一种高精度三维数据实时渐进式渲染方法和系统 |
-
2012
- 2012-09-20 CN CN201210353194.4A patent/CN102930594B/zh active Active
Non-Patent Citations (1)
Title |
---|
任远红: "基于LOD的地形算法的研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104680573A (zh) * | 2015-01-20 | 2015-06-03 | 西安科技大学 | 一种基于三角网格简化的纹理映射方法 |
CN106153030A (zh) * | 2015-03-25 | 2016-11-23 | 高德软件有限公司 | 一种导航电子地图数据存储方法及装置 |
CN105303610A (zh) * | 2015-12-05 | 2016-02-03 | 中国航空工业集团公司洛阳电光设备研究所 | 一种用于机载地形三维图像绘制的地形数据处理方法 |
CN105303610B (zh) * | 2015-12-05 | 2018-08-17 | 中国航空工业集团公司洛阳电光设备研究所 | 一种用于机载地形三维图像绘制的地形数据处理方法 |
CN106777816A (zh) * | 2017-01-18 | 2017-05-31 | 上海核工程研究设计院 | 一种场景动态模型管理系统 |
CN106777816B (zh) * | 2017-01-18 | 2023-06-02 | 上海核工程研究设计院股份有限公司 | 一种场景动态模型管理系统 |
CN109389678A (zh) * | 2017-08-08 | 2019-02-26 | 龙芯中科技术有限公司 | 确定操作焦点坐标的方法、装置及设备 |
CN107689189A (zh) * | 2017-08-31 | 2018-02-13 | 城市生活(北京)资讯有限公司 | 一种三维电子地图的绘制方法及装置 |
CN107689189B (zh) * | 2017-08-31 | 2019-09-27 | 城市生活(北京)资讯有限公司 | 一种三维电子地图的绘制方法及装置 |
CN108053475A (zh) * | 2018-02-09 | 2018-05-18 | 城市生活(北京)资讯有限公司 | 一种三维地形显示方法及装置 |
CN110555085A (zh) * | 2018-03-29 | 2019-12-10 | 中国石油化工股份有限公司 | 一种三维模型加载方法与装置 |
CN110555085B (zh) * | 2018-03-29 | 2022-01-14 | 中国石油化工股份有限公司 | 一种三维模型加载方法与装置 |
CN108986212A (zh) * | 2018-06-21 | 2018-12-11 | 东南大学 | 一种基于裂缝消除的三维虚拟地形lod模型的生成方法 |
CN108986212B (zh) * | 2018-06-21 | 2022-05-13 | 东南大学 | 一种基于裂缝消除的三维虚拟地形lod模型的生成方法 |
CN109992640A (zh) * | 2019-04-11 | 2019-07-09 | 北京百度网讯科技有限公司 | 位置网格的确定方法及装置、设备及存储介质 |
CN110633282A (zh) * | 2019-09-18 | 2019-12-31 | 四川九洲空管科技有限责任公司 | 一种空域资源多级三维网格化方法及工具 |
CN111105492A (zh) * | 2019-12-18 | 2020-05-05 | 四川大学 | 一种地形平滑过程中瓦片高程数据快速调度方法 |
CN111105492B (zh) * | 2019-12-18 | 2020-11-17 | 四川大学 | 一种地形平滑过程中瓦片高程数据快速调度方法 |
WO2022257594A1 (zh) * | 2021-06-10 | 2022-12-15 | 广东奥普特科技股份有限公司 | 一种高精度三维数据实时渐进式渲染方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102930594B (zh) | 2015-06-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102930594B (zh) | 基于保持边界的计算机三维地形网格的绘制方法 | |
CN110827402B (zh) | 一种基于栅格化技术的类建筑物三维模型简化方法及系统 | |
CN105336003B (zh) | 结合gpu技术实时流畅绘制出三维地形模型的方法 | |
CN103077549B (zh) | 一种基于kd树的实时大规模地形可视化实现方法 | |
US20140152664A1 (en) | Method of rendering a terrain stored in a massive database | |
CN103559374B (zh) | 一种多子网格模型上进行面分裂型曲面细分的方法 | |
CN102360515B (zh) | 一种三维模型渐进网格数据组织方法 | |
CN102867331A (zh) | 面向gpu的大规模地形快速绘制方法 | |
CN101976468B (zh) | 一种多分辨率动态地形可视化方法及系统 | |
CN106251400B (zh) | 一种基于多四边形网格构建地形图的方法及装置 | |
CN103268634B (zh) | 一种外存模型基于顶点聚类的快速并行自适应简化方法 | |
CN103679807A (zh) | 一种带边界约束的散乱点云重构方法 | |
CN104658027A (zh) | 一种面向不规则海洋流场数据的三维流线动态可视化算法 | |
CN104966317A (zh) | 一种基于矿体轮廓线的三维自动建模方法 | |
CN102117494A (zh) | 一种三维建模中的数据分层分块方法及装置 | |
CN107918957A (zh) | 一种保持结构和纹理特征的三维建筑模型化简方法 | |
CN112328722B (zh) | 电力变电站三维gis和模型数据快速加载方法和系统 | |
Zheng et al. | A morphologically preserved multi-resolution TIN surface modeling and visualization method for virtual globes | |
CN101169871A (zh) | 面向大规模静态场景实时渲染的hlod预处理方法 | |
CN112102486A (zh) | 一种基于合并根节点的倾斜摄影数据lod重建方法 | |
CN102096946A (zh) | 全球六边形离散格网上的矢量数据绘制方法 | |
KR102454180B1 (ko) | 3차원 지형 데이터를 생성 장치 및 방법 | |
Wang et al. | An improving algorithm for generating real sense terrain and parameter analysis based on fractal | |
CN105957142A (zh) | 一种面向TIN构建的Hilbert排序方法 | |
CN108986212A (zh) | 一种基于裂缝消除的三维虚拟地形lod模型的生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |