一种林分生长与经营交互可视化模拟方法
【专利摘要】本发明涉及一种林分生长与经营交互可视化模拟方法，包括如下步骤：步骤一建立全林分生长模型、径阶生长模型、单木生长模型和测树因子间估计模型所需的变量；步骤二构建自定义林分生长模型解析器；步骤三林分现状与生长变化后各变量的自动动态存储与提取计算；步骤四构建基于自定义林分生长模型的林木个体可视化模型参数计算法，实现自定义生长模型下的交互式林分生长可视化模拟；步骤五采用GDI与绘图技术以及MOGRE渲染引擎技术，对林分经营前后的林分结构、林分二维状态及三维场景进行可视化模拟，实现林分生长与经营交互可视化模拟。
【专利说明】一种林分生长与经营交互可视化模拟方法
【技术领域】
[0001]本发明属于林分生长可视化模拟与林业经营管理可视化模拟两个【技术领域】，具体涉及一种林分生长与经营交互可视化模拟方法。
【背景技术】
[0002]现有林分生长与林分经营模拟，均以单一的生长模型为模拟对象，不能适应不同林分生长特点的需求，同时林分生长与林分经营之间缺乏有效的交互，可视化模拟程度低等缺点，缺乏一种既能满足不同林分生长特点，又能直观反映林分生长与经营交互的可视化模拟方法，以提高林分经营管理效率和水平。
[0003]本发明综合全林分生长模型、径阶分布模型与单木生长模型所需的变量，利用编译原理中的词法、语法及语义分析等理论，构建自定义林分生长模型解析器；结合数据库技术，建立林分现状与生长变化后各变量的自动动态存储与提取算法，参考林分(林木)各调查因子间、林木个体可视化模型参数间关系模型，构建基于自定义林分生长模型的林木个体可视化模型参数估计算法，实现自定义生长模型下的交互式林分生长可视化模拟；采用GDI+绘图技术与MOGRE渲染引擎技术，对林分经营前后的林分结构、林分二维状态及三维场景进行可视化模拟，实现林分生长与经营交互可视化模拟。
【发明内容】
[0004]针对现有技术中的不足，本发明的目的在于提供一种林分生长与经营交互可视化模拟方法，实现林 分生长与经营的交互可视化模拟。
[0005]为了实现上述目的，本发明采取的技术方案是:
[0006]一种林分生长与经营交互可视化模拟方法，包括如下步骤:
[0007]步骤一建立全林分生长模型、径阶生长模型、单木生长模型和测树因子间估计模型所需的变量；
[0008]步骤二构建自定义林分生长模型解析器；
[0009]步骤三林分现状与生长变化后各变量的自动动态存储与提取计算；
[0010]步骤四构建基于自定义林分生长模型的林木个体可视化模型参数计算法，实现自定义生长模型下的交互式林分生长可视化模拟；
[0011]步骤五采用⑶I与绘图技术以及MOGRE渲染引擎技术，对林分经营前后的林分结构、林分二维状态及三维场景进行可视化模拟，实现林分生长与经营交互可视化模拟。
[0012]所述全林分生长模型建立步骤为
[0013]建立林分断面积生长模型
[0014]G = 26.9899S10-7686 {1-exp [-0.0006802 (SDI/100) 18675 X (t-2.6)]}04803(I)
[0015]其中，G为林分断面积；SI为地位指数；SDI为林分密度指数；t为预估林龄；
[0016]SDI = NX (D/D0) 0(2)
[0017]其中，N和D为现实林分每公顷株数和平均胸径，Dtl为基准平均胸径，取DpOcm，3为最大密度线，取值为1.2248 ；
[0018]建立优势木高预估模型
[0019]林分优势木高的生长与林分密度关系相对较小，主要依赖于立地质量，立地质量
越闻，优势木闻生长越快，可表达为:
[0020]
【权利要求】
1.一种林分生长与经营交互可视化模拟方法，其特征在于，包括如下步骤: 步骤一建立全林分生长模型、径阶生长模型、单木生长模型和测树因子间估计模型所需的变量； 步骤二构建自定义林分生长模型解析器； 步骤三林分现状与生长变化后各变量的自动动态存储与提取计算； 步骤四构建基于自定义林分生长模型的林木个体可视化模型参数计算法，实现自定义生长模型下的交互式林分生长可视化模拟； 步骤五采用GDI+与绘图技术以及MOGRE渲染引擎技术，对林分经营前后的林分结构、林分二维状态及三维场景进行可视化模拟，实现林分生长与经营交互可视化模拟。
2.根据权利要求1所述的方法，其特征在于:所述全林分生长模型建立步骤为 建立林分断面积生长模型G = 26.9899S10-7686 {1-exp [-0.0006802 (SDI/100) 18675 X (t_2.6)]}a 4803(I) 其中，G为林分断面积；SI为地位指数；SDI为林分密度指数；t为预估林龄； SDI = NX (D/D0)0(2) 其中，N和D为现实林分每公顷株数和平均胸径，Dtl为基准平均胸径，取^=20(^, ^为最大密度线，取值为1.2248 ； 建立优势木高预估模型 林分优势木高的生长与林分密度关系相对较小，主要依赖于立地质量，立地质量越高，优势木高生长越快，可表达为:H 二 a沿-e-(3)
其中，SI 为地位指数；A 为林分年龄；&1=47.784 ；a2=0.00102 ；a3=l.124 ；a4=l.077 ； 建立林分株数自稀疏模型 依赖于优势木高和前期密度的枯损方程，表达为: N = (Na0+H!rn-H才(4) 其中，N为期末林分密度；氏为林分期末对应的优势木高为林分初始密度即单位面积株数为林分期初对应的优势木高；a=-0.36785 ；b=-0.0000000895 ；c=3.56889 ；建立林分平均胸径预估模型 在未进行人为干预的林分中，林分平均胸径与株数、优势木高之间存在幂函数关系，所建模型表达为: Dg=axHb'N--Hb2(5) 其中，Dg为林分平均胸径；H为林分优势木高；N为林分单位面积株数；&1=1.6601 ；Id1=L 6744 ；a2=0.0862 ；b2=0.4668 ； 建立人为干预后林分平均胸径和株数估计模型 当经营措施实施后，人为地改变了林分的生长过程，需用动态模型进行估计:
D2 = D1 + 0.3451 *?0 7153 (1- C^u3346am)"0.618 *￡^.083244( 6 )
N2 = JV1 -2782.601/邊1582(1-e—u2472uWaj)211.81 *eM322妈(7 )其中，D1为期初林分平均胸径；D2为期末预估林分平均胸径A1为期初林分单位面积株数；N2为期末林分单位面积株数A1为期初林分年龄；AA为预估期长度；SI为地位指数。
3.根据权利要求1所述的方法，其特征在于:所述径阶分布模型建立步骤为 建立Weibull分布概率密度函数: bK b 1 (8) 式中，a为位置参数；b为尺度参数；c为形状参数； Weibull分布函数的一阶原点距为林分算术平均胸径D，二阶原点距为林分平均胸径的平方Dg2，即: D = a+b r (l+1/c)(10) /.),:+ 2, t ) + 2c")r(i + 1 c)-\-cr(11 ) a可定义为林分最小直径，根据同龄纯林直径结构规律，用林分平均胸径的0.45倍作为a的近似估计值，在得知D和Dg的情况下，两式反复迭代即可求得参数b和c ; 林分算术平均胸径和林分平均胸径间有着密切的线性关系，表达如下: D = -0.36756+1.01301Dg(12) 利用上式即完成对林分算术平均胸径的估计，各径阶的株数按下式求得:' V b(13) 其中，叫为第i径阶内理论株数；N为林木总株数;W为径阶距；Xi为第i径阶中值。
4.根据权利要求1所述的方法，其特征在于:所述单木生长模型的建立步骤为 建立胸径连年生长量模型 基于生长量修正法建立的杉木人工林单木胸径连年生长量模型来对中林分经营前后的林木生长状况进行预估，表达式如下:
dD/dt = 1.56364RS0-51489RD0-02685 (0.13 2 7 8S10-88556Da 2302-0.0158SI。.73315D)(14) 式中，SI为地位指数；D为胸径；RS为相对植距；RD为相对优势度； 建立胸径生长量模型 将胸径生长量作为胸径、年龄和竞争指数的函数，表达如下: LnZc = a0+a1D+a2A+a3CI+a4/A(15) a,=Jix 丄(16)
>=i xi lH 其中，CIi为简单竞争指数；X1、&分别表示竞争木i和对象木j的测树因子，Lu表示竞争木i和对象木j的的距离；ZD为胸径生长量；D为胸径;A为林龄； 建立对象木3年胸径生长量模型 将胸径3生长量作为竞争指数的函数，表达如下: Zd = 5.4137ei 3733CI(17) 其中，Zd代表林木3年胸径生长量；CI表达式如式(23)所示，其中，n=8。
5.根据权利要求1所述的方法，其特征在于:所述测树因子间估计模型的建立步骤为 对林木树高分布估计林木胸径和树闻具有一定的关系，以直角坐标系中的X轴表不胸径，y轴表不树闻，将林分中的林木绘制到坐标图上，并依据散点的分布绘制一条平滑曲线，这条曲线称为树高曲线，树高曲线能够明显的反映出树高随着胸径的变化规律， 常用的胸径-树高曲线方程有: H = ao+a*log (D)(18) H = ao+af (D) +a2*D2(19) .11 =^d"' (20) 其中，H代表树高；D代表胸径；％，&1，a2为参数。 依据林分调查资料，绘制H-D散点图，根据散点图分布选择几个树高曲线进行拟合，从中挑选拟合度最优的树高曲线方程作为研究区域的H-D树高曲线方程；采用胸径-树高方程: H = 0.4921*S°-6309^D273 (r = 0.9734)(21) 对林木冠幅分布估计 冠幅的模拟采用张铁林得出的杉木人工林胸径冠幅关系理论公式: 大径材: Cff = 0.83008D0-47306(22) 中径材: Cff = 0.75387D0-49336( 23) 小径材: Cff = 0.81540D0 46136(24) 杉木冠幅面积估计模型表达式如下: CWa = D/(2.8617-0.03498D)(25) 其中，D为胸径；CW为诞幅；CWA为诞幅面积； 建立材积估计模型 采用农林部标准LY208-77中的杉木二元材积公式: V=0.000058777042*Dl 9699831^Ha 89646157( 26 ) 式中:D为胸径；H为树高；V为材积。
6.根据权利要求1所述的方法，其特征在于:所述自定义林分生长模型解析器包括C语言编程模块、词法分析模块、语法分析模块与语义分析模块，其中执行过程如下 首先，通过定义模型因变量与自变量，对用户定义的林分生长模型进行词法分析，确保字符表达正确； 其次，执行语法分析，确保模型表达无误； 再次，执行语义分析，实现类型审查，并将模型表达为便于执行运算的后缀表达式形式，结合动态实时获取变量值的林分数据表，进行运算； 最后，用执行结果更新林分数据表，完成数据更新。
7.根据权利要求1所述的算法，其特征在于:以WF技术为基础，在构建林分生长活动、林分结构分析活动与林分经营措施活动的基础上，通过设计林分生长活动、林分结构分析活动及林分经营措施活动与林分数据表间的接口，实现各个活动与林分数据表间的数据传输，以图形化、流程化方式， 实现模拟三者间信息交互的可视化模拟，其中实现林分结构分析、林分经营措施活动、3维场景模拟需要的是林木个体信息，当林分生长模型是全林分模型或径阶模型时，在设计林分生长活动与林分数据表接口的时候，涉及到林分整体信息到林木个体信息含林木3D可视化模型构建所需参数的估计的算法，以实现任意类型生长模型下的林分数据表更新，接口设计采用数据库技术，以实现林分生长、林分结构分析与林分经营措施所需 变量的自动动态存储与提取。
【文档编号】G06F19/00GK103678870SQ201310438526
【公开日】2014年3月26日 申请日期:2013年9月24日 优先权日:2013年9月24日
【发明者】张怀清, 鞠洪波, 李永亮, 刘闽, 蒋娴 申请人:中国林业科学研究院资源信息研究所