[发明专利]一种基于GPU加速的多分辨率大规模森林演替过程仿真方法

说明书

技术领域

本发明涉及森林动态模拟仿真领域，尤其是一种大规模森林演替过程仿真方法。 

背景技术

对于森林动态模拟仿真，国内外已经提出了很多生长模型，这些模型中，森林的动态仿真大部分都基于以下几个环节：种子分布阶段，萌芽阶段，生长阶段，死亡阶段。在这些阶段，其中包括光照模型的计算，环境资源计算，成年树生长计算，幼年树生长计算等。传统森林仿真不具有科学真实性，森林场景只由主观想象绘制的树木图像简单叠加而成，未充分体现出森林场景的复杂性与客观真实性。另一方面，虽然有很多基于林学数学模型的计算机模拟森林动态仿真，但很多在森林场景的规模上都未达到很大规模，在达到大规模场森林场景的动态仿真时，其方法运行效率不高，通常需要耗费高代价的时间来获得模拟结果，其未能充分发挥现在PC计算机硬件的潜力，尤其是图形硬件GPU处理器上越来越通用的并行计算能力。 

现有方法对上述这些森林动态生长阶段的加速计算主要基于GPU基础上，对于种子分布阶段，本系列方法发明所借助的模型是2000年Philip T.LePage等人提出的WeiBull Function模型，由于对于某块样地的种子数量涉及整片森林场景的树木，因此其计算量是非常庞大的，现在已有的加速方法是基于在CPU上建立四叉树数据结构实现数据优化，进行计算的方法。 

对于成年树阶段，借助的生长模型是2009年K.David Coates等人提出的温带植物资源竞争模型。其中在该成年树生长模型中，对于计算机的模拟主要耗时的是成年树领域内树木竞争的计算与光照资源竞争的计算，在该模型中，其中成年 树领域树的竞争以NCI指数所表示，光照资源以光照指数(GLI)所表示。对于成年树的光照指数(GLI)计算，现有使用基于光线与树冠求交的方法进行实现，即使在GPU上实现求交，由于场景规模的巨大，其计算量与所需线程量也相当庞大。 

现有多分辨率技术主要用于图形图像中的网格简化，通过构造原始网格模型的多个逼近表示，结合硬件资源的绘制能力和绘制误差选择最优的细节层次(level-of-detail，LOD)进行绘制，并在保证绘制速度的前提下尽可能提高绘制质量。 

随着计算机硬件的不断发展，GPU并行计算能力与通用性的提高，另外由于森林场景中树木个体之间具有相同的生长阶段，生长模型公式具有相同性，设计适用于GPU硬件架构的加速方法对于高效大规模的森林动态仿真具有很高的实用价值，对于高效的管理大片森林资源，预测未来发展趋势具有良好的应用价值。 

发明内容

为了克服已有森林动态模拟仿真方法的正确性较差、运行效率较低的不足，本发明提供一种正确性良好、提升运行效率的基于GPU加速的多分辨率大规模森林演替过程仿真方法。 

本发明解决其技术问题所采用的技术方案是： 

一种基于GPU加速的多分辨率大规模森林演替过程仿真方法，所述仿真方法包括以下步骤： 

1)，建立树木数据组织与样地数据组织； 

根据GPU架构对树木数据进行组织，形成利于GPU计算的数数据组织结构：所有计算的树木数据结构将被设计成数组形式；在多分辨率数据组成中聚类结点第k层与第k+1层的数据表示，对于第k层单位结点数据由m个k-1层结点数据累加而成，对于每个单位结点数据，由s种树种对应单位成年树结点数据构成； 

每棵成年树或者幼年树x的树冠被建模成圆柱体，它们在地表的投影都是一个圆C(x)，设样地单元边长为r，则样地单元网格的面积即r²，对于投影完全包含k个样地单元网格的树木，则其暴露面积即为kr²，为正确计算不同树木发生遮挡时暴露面积，设树木树冠上顶面离地面距离为H(x)，所有样地网格单元标记为M_ij，在理想不想交投影区域，则如果某样地网格M_ij完全包含在树木投影C(x)中，则将M_ij分配给树木x，否者不对其分配；在投影相交区域，对于样地网格M_ij被多个投影所覆盖，该样地网格将会被分配给能够完全包含它的H(x)最大的那棵树木x；