EXECL表格软件在样地调查计算中的应用.docxVIP

? ? EXECL表格软件在样地调查计算中的应用 ? ? 苏振海 摘 要：本文采用FORECAST、LOOKUP和LN等函数编写了样地调查计算模板，并通过设置数据有效性和保护工作表等功能，设置了数据输入的有效性检验，将模板需要输入的内容降到最低，不再需要复制粘贴过程，有效提高了工作效率，降低了误输入造成的错误，通过实验取得了良好的效果。 Key：样地调查;EXECL;计算 ：S159 ? ? ? ：A DOI：10.19754/j.nyyjs.20200515023 1 研究目的 样地调查在森林规划调查、森林采伐设计、森林资产评估、征占用林地设计中是最常用的方法[1]，特别是对人工林调查设计中，比角规测树方法更为准确。在样地调查中，最重要的是调查样地的活立木蓄积，这就需要在外业计算平均胸径（平方平均数），通过平均胸径确定测高树分布，根据这2个因子利用二元材积表或材积模型进行样地蓄积计算。 2 模板基本功能 记录胸径及树高实测记录;自动计算胸径平方平均值;自动计算各径阶树高及样地平均树高;根据树种自动计算单株蓄积、样地蓄积并推算小班蓄积、小班株数等;对可能产生的错误输入或误修改进行保护、检测并提示。 3 功能实现 3.1 树种选择录入 采用模型计算蓄积时，树种选择决定了蓄积模型的参数，为方便实现自动计算功能和输入，采取下拉菜单形式设计，通过点选就可输入树种。具体方法如下。 以笔者设计的表格样式为例（下同），在A1单元格输入“树种选择”，选择B1作为显示树种名称单元格，选中B1依次选择工具栏中的“数据”-“数据有效性”，在弹出的界面中选择“设置”，在“有效性条件”中选择“允许”-“序列”在“来源”下选择录入表单的位置，如=$N$3∶$N$8（笔者只实验了6个树种）。为方便下面的编制，需在树种名称前增加序列值如“1辽东日落”或“a辽东日落”均可，并且表单按升序排列。后面会解释原因。 3.2 胸径及测高记录 在单元格A3～f3中分别输入“径阶、株数、加权平方和、径阶高、径阶单株蓄积、径阶蓄积合计”;在单元格I3～K3中分别输入“树种、胸径、树高”。A4～A36分别输入6～70的整化径阶值（2cm），如胸径径阶大于70，已经超出蓄积模型适用范围，可采取其它方法近似计算。 3.3 计算胸径平方平均值 在C4单元格输入公式：=A5*A5*B5或=A5^2*B5或POWER（A5，2）*B5均可，以下至C36复制公式即可。 C39中输入“=SUM（C4∶C36）” （合计加权平方），B39输入“=SUM （B4∶B36）” （合计株数），D39输入“= （C39/B39）^0.5” （胸径平方平均值）。 3.4 计算各径阶平均树高 各径阶的平均树高一般采用数式法，拟合树高曲线再通过曲线求算各径阶平均树高[2]。笔者没用采用散点图添加趋势线拟合的方式进行求算，因为拟合后的曲线参数每次都不一样，需要在公式中手动修改参数值才能计算径阶平均高，增加了产生错误录入的机会，也不方便操作。根据多年的制表经验，辽宁省的树高曲线一般都满足对数模型，因此直接将胸径进行对数变换，采用FORECAST函数直接求算回归值。FORECAST函数的用途是根據一条线性回归拟合线返回一个预测值。表达式为FORECAST （x，known_ys，known_xs），其中X为需要进行预测的数据点，Known_ys为因变量数组或数据区域，Known_xs为自变量数组或数据区域。 在L4单元格中输入公式“=LN（J4）”L5～L36复制公式。进行实测胸径变换。 在H4单元格中输入公式“=LN（A4）”H5～H36复制公式。进行整化径阶变换。 在D4单元格中输入公式“=FORECAST（H4，$K$4∶$K$36，$L$4∶$L$36）”D5～D36复制公式。进行径阶平均高求算。 求算结果与拟合曲线后带入整化径阶求算径阶平均高完全一致。如果后期需要打印输出树高曲线及模型，可采用散点图加趋势线方式，本文不再赘述。在E39单元格中输入公式“=FORECAST（LN（D39），K4∶K36，L4∶L36）”，计算样地平均树高。 3.5 计算蓄积 3.5.1 单株蓄积计算 辽宁省二元立木蓄积模型采用的都是山本式V=c0Dc1Hc2参数较少，使用方便。根据3.1中树种选择录入结果自动查找3个参数。方法如下： 在C40单元格中输入公式“=LOOKUP（B1，N3∶N8，O3∶O8）”查取参数c0 在C41单元格中输入公式“=LOOKUP（B1，N3∶N8，P3∶P8）”查取参数c1 在C42单元格中输入公式“=LOOKUP（B1，N3∶N8，Q3∶Q8）”查取参数c2 O3～Q8单元格为录入好的相应树种的二元材积表参数值。 这里LOOKUP函数，要查询一个明确的值的时候，查询列必须按升序排