角规与样圆测树内业计算

 

本软件是我在1996年森林资源二类调查时编写的，原是用ＤＯＳ版本的Foxpro2.6编写的，现用VFP6.0,在原有的计算功能基础上进行了升级，改造成For Windows版本，并对输入界面和计算功能进行了改进和扩充，加进了样圆测树和全林分每木检尺计算功能。原有版本除了输入外业调查测树因子外，还要输入各树种不同径阶的出材率，并且只能得到单位面积（亩）的蓄积量和出材量，要再进行手工伐区汇总。升级后的版本对查表计算过程进行了重大改进，现只要输入外业调查测树因子，如树种、胸径、树高、角规断面积（样圆为每亩株数）、坡度、比例成数、伐区（小班）面积等，其余过程全部由程序自动完成，原始数据输入完毕后按“计算”，只需瞬间即可得到伐区（小班）测树因子综合值和伐区汇总表。    

一、运行界面及输入方法

    本软件将角规测树与样圆测树集成在同一个界面，全林分单独一个界面。这里给出角规与样圆测树运行界面和样例：

1、 根据外业调查种类，选择角规或样圆测树 。

　

图1：角规测树、样圆测树运行界面

　

2、选择从表输入或逐点输入 。

如果选“从表输入”，出现以下输入框：

图2 从表输入并计算得到综合结果的表样式                           

   

如果选“逐点输入”，则逐点输入框有效：

图3：逐点输入框

    逐点输入外业调查数据，每输入一点后按“添加”，把输入的数据添加到表。此时，表输入框底色变灰色，但仍然处在可编辑状态，发现上一点数据输错，可在表输入框中直接修改。
4、输入完毕并确认无误后，按“计算”，然后按“打印结果”
5、如果出材率值不适用于当地，请按“修改数据源表”，对出材率进行修改。
6、可一次性输入外业调查表数据，不论是单一树种的小班还是多个树种的小班。
7、表同时也用于存放计算结果，是共享的，计算前，如果发现表已有数据，请按“清除”以清除上次计算结果。
8、 输入框“把计算结果输出到：”，指把计算结果表另存为XXX.dbf文件，输入一个文件名，不要扩展名，然后按确定。注意：文件名不能与系统中的表文件名相同，注意文件名已存在提示信息。系统表文件名如下：

          dmjdata.dbf      用于存放角规断面积蓄积量速见表，表格式如图5：
          sjbdata.dbf       用于直接从表输入，同时用于存放计算结果的表，格式见表输入框。
     smcclb.dbf        杉木出材率表。
          mwscclb.dbf      马尾松出材率表。
          kyscclb.dbf           阔叶树出材率表

二、计算方法

      1、查表
    采用内差法查福建省角规测树断面积蓄积量速见表法，速见表格式见图5。字段：
Trk代表树种，H3代表树高为3米，H3x代表修正系数，其余类推。限于篇幅，只列出到H9。


        图4：角规断面积蓄积量速见表

算法例：
   求算杉木角规断面积7.8平方米(G/ha)，树高8.5米（H）的角规点蓄积量，先求7.8平方米（G/ha）－8米（H）、7.8平方米（G/ha）－9米（H）角规点蓄积量分别为7.8X0.318-0.05=2.4立方米／亩和7.8X0.345-0.05=2.6立方米／亩；然后，内差树高得角规点每亩蓄积量为2.5立方米／亩。

2、样圆测树立木单材积公式：

    样圆法求算蓄积量采用的单材积公式。

    主要树种（组）立木单材积公式：
    杉木立木二元材积公式
    v= 0.00005806*d^ 1.9553351*h^ 0.89403304
    马尾松立木二元材积公式
    v=0.000062341803*d^1.8551497*h^0.95682492
    阔叶树立木二元材积公式
    v=0.00005276429*d^1.8821611*h^1.0093166

    单材积乘于每亩株数即得每亩蓄积量。

3、各树种出材率表

      图5：杉木出材率表

    树高级1：指树高大于“树高级1”的值时查表范围 H＞h1。
    树高级2：指树高小于或等于“树高级1”和大于或等于“树高级3”时查表范围。h3<=H<=h1。
    树高级3：指树高小于“树高级3”的值时查表范围。H<h3。

3、综合
根据各调查点权重（比例成数），对各调查点进行加权平均。

    树高综合：
    Hbar=(h1*r1+h2*r2+... +hn*rn)/(r1+r2+... +rn)
    胸径综合：
    Dbar=(d1*r1+d2*r2+... +dn*rn)/(r1+r2+... +rn)
    胸径和树高加权值要除与比例成数之和，因为，具有某一树种的调查点数量可能比调查总点数少。
    r1+r2+...+rn<=1 rn为某点比例成数
    亩蓄积综合：
    Vpm=vpm1*r1+vpm2*r2+...+vpmn*rn

4、打印

按“打印结果”，出现以下对话框：

图6 打印任务选项

 

打印格式如下：（用空白纸张）

图7 角规测树计算表打印预览

图8 伐区汇总表打印预览

合计栏的树种、胸径、树高为优势树种的值，其余为各树种的综合值累计。表大小为16K纸横向。

如果选“打印在外业纪录表”，外业纪录表１６Ｋ竖放（见角规测树调查纪录表）