[发明专利]基于带惩罚项的局部二次误差测量的地物网格化简方法

说明书

技术领域

本发明涉及渐进网格化简模型和三维格网化简技术在地理空间物理三维表达领域。尤其涉及一种基于带惩罚项的局部二次误差测量的地物网格化简方法。

背景技术

随着科学技术的发展，硬件计算能力越来越高，同时人对图像拟真度的要求越来越严格。随着云计算的发展，利用大规模集群对三维模型进行离线处理成为游戏、视频甚至是电影工业的重要的场景渲染方式。移动设备随着计算机工业的发展也得到长足的进步，但是由于电池等原因，移动设备上的大场景渲染依旧困难，业界主流的方式都是通过过度化简以减少移动设备三角形计算消耗或者将渲染完全放在服务端运行，以视频的形式提供给终端用户。前者虽然很大程度减少了移动设备的计算负担，但过度失真导致所渲染场景与真实场景相差过大，后者能够提供与桌面渲染媲美的场景渲染、光影特效等，但是受限于移动网络的传输带宽，在网络状况欠佳的情况下会造成跳帧、卡顿等状况，同时也不适合交互性较高的移动三维地图。

目前对三维场景显示采用细节层次、遮挡裁剪的方式减少额外的计算消耗。细节层次的主要技术通过网络化简对三维模型进行化简分层，常见的网格简化算法一般分顶点操作和边操作。目前对于网格简化算法有：顶点删除算法、网格重新划分算法、区域合并算法、小波分解算法、顶点聚类算法、包络网格算法、边折叠算法等。顶点删除操作涉及的面较多，一般对网格形状变化比较大。网格重新划分、区域合并、小波分解、顶点聚类、包络网格等算法都比边折叠需要的计算量更大，因此在近几年的网格简化算法研究主要针对边折叠算法展开。

发明内容

本发明的目的是为克服现有技术存在的问题，提供一种基于带惩罚项的局部二次误差测量的地物网格化简方法。

基于带惩罚项的局部二次误差测量的地物网格化简方法包括如下步骤：

1)利用渐进网格算法对精细网格进行化简，确定被折叠边选择的折叠代价，根据顶点的边与边相关面的夹角以及边长度获取折叠代价；

2)对折叠边按照折叠代价从低到高进行初步排序；

3)从已排序的队列中，依次获取折叠边，得到边的局部二次误差测量度；

4)获取顶点曲率值及相关面最大曲率值，选取最大曲率值与调整系数，获得QEM的惩罚项，根据该惩罚项，确定折叠边的最终折叠代价；

5)从队列中取出未确定最终折叠代价的边，依次确定直至所有边的折叠代价确定完毕；

6)调整惩罚项的系数，利用正交投影判断不同系数下化简模型与原模型的相似率，评价系数可用性；

7)利用地物模型对系数进行训练，获取最优系数，并应用于最终地物化简中。

所述的步骤1)包括：根据渐进网格算法对地物网格进行化简，优化全局能量函数：

E(M)＝E_dist(M)+E_spring(M)+E_scalar(M)+E_disc(M)(1)

网格化简中判断一条边是否可以折叠，取决于它的边长与曲率值的乘积，从模型中选取一条边(u，v)，代价公式：

其中，T_u为顶点包含u的三角形的集合，T_uv为边包含(u，v)的三角形的集合，即和normal_s为平面f和s的法向量，计算得到边折叠代价，写入边(u，v)缓存中，循环步骤1)直到所有边代价计算完毕。

所述的步骤2)包括：在确定折叠边折叠代价之后，利用优先队列，以折叠代价作为条件因子进行升序排序，利用R*树索引对空间进行索引，分割空间，形成多优先队列。

所述的步骤3)包括：从步骤2)得到的多优先队列中获取最小折叠代价边，取其端点，计算端点到边相关面的距离平方和作为该顶点的二次误差测量度，公式：

其中为顶点u到平面f的距离平方，u＝[u_xu_yu_z 1]^T，f表示的空间中的平面:

ax+by+cz+d＝０，(a²+b²+c²＝1) (4)

因此基于局部二次误差测量的折叠边代价定义为：

其中，为顶点u的二次误差测量的4×4对称矩阵。

所述的步骤4)包括：

5.1获取顶点曲率值及相关面最大曲率值