-
树冠外轮廓模型是指以树木树冠任意位置处半径为因变量,以冠幅、冠长、树高或胸径等因子为自变量的函数,用来模拟树冠形状变化规律的数学表达式[1]。树冠外轮廓模型不仅能够估计树冠任意位置处的树冠半径和推导计算整个树冠体积,而且可以驱动林分树木树冠三维可视化模拟。此外,树冠形态与树冠体积对于间接估算树木树冠部分生物量具有重要意义[2]。为了描述树冠外轮廓形态变化规律和推导计算树冠体积,国内外学者常利用以下方法预测树冠相关变量,包括简单几何形状模拟、经验模型方法、树冠轮廓模型积分法、激光扫描方法等。早期学者们利用规则几何体描述树冠外轮廓二维形状,常采用如下公式定义树冠形状:RCR=(1-RDH)k,其中:RCR表示相对树冠半径,RCH为从树冠基部到任意位置处树冠长度与最大树冠长度的比,k的取值决定了其几何形状,当k为0、0.5、1.0、1.5时分别表示圆柱形、抛物线形、圆锥形或凹面形等,该模型被广泛应用于道格拉斯杉Pseudotsuga douglasii[3]、西部铁杉Tsuga heterophylla[4]、火炬松Pinus taeda[5]、华北落叶松Larix principis-rupprechtii[6]等树冠形状模拟和体积计算研究中。PRETZSCH等[7]、SADONO[8]将树冠分为“阳冠”“阴冠”两部分,分别采用抛物线体和圆锥体表示。由于树冠形状并不是简单规则形式,不同树种树冠形状并不相同,树冠体积大小与冠幅、冠长、胸径、树高等林木因子之间存在相关关系,研究者针对不同树种利用经验模型方法构建幂函数、指数函数等多种形式树冠体积模型,如楠木Phoebe zhennan[9]、樟子松Pinus sylvestris var. mongolica[10]、杉木Cunninghamia lanceolata[11]、长白落叶松Larix olgensis[12]、欧洲赤松Pinus sylvestris和云杉Picea asperata[13]、黑荆树Acacia mearnsii[14],也有研究者通过构建枝条长度、角度、半径等枝条属性经验方程间接预测树冠半径和模拟树冠形态变化规律,如李凤日[15]、LI等[16]、GILMORE等[17]、姜立春等[18],但是存在树冠枝条测量困难且误差较大问题。三维激光扫描等技术的发展,不仅提供了快速数据获取方式,而且能够更加准确获取树冠数据[19-21],但由于数据获取成本高、处理困难、数据量大等限制了它的使用。除了以上方法,为了更加准确预测树冠体积,一些学者利用多项式[22]、幂函数[23]、修正Beta曲线[24]、连续分段函数[25]等曲线形式描述树冠外轮廓形态,并将树冠曲线围绕树干轴旋转得到的旋转体,通过对旋转体积分推导计算树冠体积[26-28]。该方法只需要简单林木因子就能够准确模拟树冠形态与树冠体积。目前,树冠外轮廓模型积分法计算体积研究,一般将树冠外轮廓与树冠体积模型单独进行研究,分别拟合不同模型方程得到2套模型参数;由于外轮廓与体积方程之间存在内在相关性,采用传统最小二乘方法分别求解模型参数,不能保证2个模型误差同时最小,无法满足模型参数估计的渐进无偏性、有效性和一致性。通过对树冠轮廓模型积分推导得到树冠体积方程,将2个方程联立构建一致性方程组,能够有效解决以上模型参数估计存在的问题。目前一致性方程组的研究主要集中在树干削度与材积一致性方程组研究方面[29-30],关于利用树冠外轮廓与树冠体积一致性模型方面研究还未见报道。本研究以福建地区杉木人工林为研究对象,采用似乎不相关回归方法,构建以最大树冠半径与相对冠长为自变量的树冠外轮廓和树冠体积一致性方程组模型,研究构建的模型能够用于驱动树木树冠三维可视化模拟,直观地反映不同生长条件下树冠生长活力及林木个体间树冠重叠程度,指导林分抚育间伐活动,同时为通过树冠体积间接推算树冠生物量提供新思路。
HTML
-
福建省位于中国东南沿海地区(23°33′~28°20′N,115°50′~120°40′E),地势东南低西北高,地形多丘陵少平原,不同区域气候相差较大,其中东南沿海属南亚热带气候,东北与西北区域属中亚热带气候,年平均气温为15.0~22.0 ℃,年平均降水量为1 400~2 000 mm,土壤类型主要包括红土壤、黄土壤、山地草甸土壤,主要用材树种有杉木、马尾松Pinus massoniana、巨尾桉Eucalyptus grandis × E. urophylla等。杉木为杉科Taxodiaceae乔木,幼树树冠尖塔形,大树树冠圆锥形,为中国长江流域、秦岭以南地区栽培最广、生长快、经济价值高的用材树种。
研究数据来自于福建省顺昌县大历和岚下林场布设的杉木人工纯林临时样地,选择不同龄组、林分密度和立地条件类型设置标准样地,共设置98块30 m × 20 m的样地,每个标准地内选择3~5株树木,共413株杉木。测量每株树木胸径(D)、树高(H)、冠幅(CW)、最大冠长(LCL)、枝下高(HCB),树冠长度(CHi)以及相应树冠半径(CRi),其中:i(i=0.10、0.25、0.50、0.75、0.90)表示从树冠基部到树梢顶端的相对位置。树冠因子测量如图 1所示,该装置由透明绘图板、三脚架以及照准装置组成,具体使用方法如下:①根据相似三角形原理,在距离所测树冠一定距离位置处(通常距离1倍树高),将该绘图板固定在三脚架之上保持板面垂直于地面,然后通过透明板观察树冠,往后移动三脚架,保持透明板与地面垂直,直到通过透明板能够观察到整个树冠。②用笔将树冠外轮廓绘制在透明板之上,首先绘制树冠枝下高位置处的点A、树冠顶端的点B以及最大树冠位置处的点C与点D,然后从A点开始到B、C点,将树冠轮廓绘制出来。并将透明板上的树冠轮廓草图复制到有计算网格方格的透明硫酸纸上。③利用测高仪器和皮尺分别测量该树冠的冠长L(LCL=A′B′)和冠幅(CW=E′F′),计算树冠冠长测量值A′B′与绘图纸上AB之间比值以及树冠冠幅测量值E′F′与CD之间比值,根据这2个比值在绘图纸上分别计算树冠冠长0.10、0.25、0.50、0.75、0.90位置处对应的树冠长度和半径值。为了尽可能减少树冠测量误差,应从多个方向观察树冠,取测量均值。通过以上调查方法收集到杉木树冠调查数据情况如表 1所示。
项目 t/a N/(株·hm-2) D/cm H/cm CW/m LCR/m LCL/m HCB/m CR/m 平均值 16.00 1 850.00 16.36 12.60 3.50 1.75 6.00 6.40 1.10 最小值 5.00 900.00 5.90 3.00 1.60 0.80 1.10 0.30 0.10 最大值 29.00 4 000.00 33.20 25.50 7.60 3.80 10.30 16.30 3.80 标准差 7.25 60.39 5.52 3.91 1.05 0.53 2.65 3.17 0.56 说明:t为年龄;N为株数密度 Table 1. Summary statistics of measurements of tree variables
-
树冠体积数据可采用分层切割法处理得到。本研究在对杉木树冠进行调查时,采用平均断面积求积法,将树冠从树梢到树冠基部按照相对冠长分6部分(0~0.10,0.10~0.25,0.25~0.50,0.50~0.75,0.75~0.90,0.90~1.00)。将树冠最上部分看作圆锥体近似计算体积,其他各部分采用平均断面积法。计算公式如下:
其中:Vc表示整个树冠体积(m3);Vi为第i部分树冠体积(m3);gi为第i区分段中央断面积(m2);li为第i分段长度(m);g′为梢头底端断面积(m2);l′为梢头树冠长度(m);n为分段个数。
-
利用SAS软件的proc model模块中似乎不相关回归过程SUR方法对树冠外轮廓模型和树冠体积预测模型的一致性方程组同时进行拟合,表 2给出了不同模型拟合的统计量,即决定系数(R2)、均方根误差(RMSE)。
模型 参数 R2 RMSE a0 a1 a2 CR VC CR VC 模型1 0.661 4(0.006 8)* 0.802 5 0.858 1 0.260 3 8.605 4 模型2 1.975 3(0.054 7)* -2.703 0(0.044 4)* 0.808 2 0.860 9 0.256 6 8.518 9 模型3 0.827 9(0.009 1)* 0.138 1(0.138 1)* -0.870 5(0.043 5)* 0.829 2 0.861 0 0.242 1 8.517 6 模型4 0.835 5(0.005 5)* -0.748 7(0.010 0)* 2.277 1(0.078 1)* 0.829 5 0.861 0 0.241 9 8.517 4 说明:a0、a1、a2为模型系数,参数拟合结果数值表示估计值(标准差);*表示P<0.000 1 Table 2. Fitting results of crown profile and crown volume models
由表 2拟合结果可知:树冠外轮廓模型拟合指标R2的对比结果为模型4>模型3>模型2>模型1,轮廓模型拟合指标RMSE的对比结果为模型4<模型3<模型2<模型1,树冠体积模型拟合指标R2的对比结果为模型4=模型3>模型2>模型1,体积模型拟合指标RMSE结果为模型4<模型3<模型2<模型1。根据R2最大且RMSE相对较小的最优模型选择标准,不管从树冠外轮廓还是体积模型结果来看,模型4显示了较好的拟合结果,最终利用模型4来描述福建地区杉木树冠外轮廓和树冠体积。
-
建模数据只能反映模型拟合的好坏,不能反映模型的预测性能。模型系统的独立性检验是采用建模时未使用的独立样本数据,对各模型系统的预测性能进行综合评价。基于表 3的参数估计值和检验数据,利用SAS软件计算各模型系统树冠外轮廓和体积的绝对误差、均方根误差。从表 3可以看出:各树冠外轮廓模型的检验指标R2对比结果为模型4>模型3>模型2>模型1,外轮廓模型RMSE对比结果为模型4<模型3<模型2<模型1,各树冠体积模型检验指标R2对比结果为模型3=模型4>模型2>模型1,体积模型RMSE对比结果为模型4<模型3<模型2<模型1,不管是从R2还是RMSE来看,模型4均优于其他模型。此外,模型4的检验指标MAE和MD绝对值最小,真实值与预测值之间误差最小,进一步验证了模型4作为杉木树冠外轮廓-体积模型的合理性。
模型 树冠外轮廓模型 体积模型 MAE MD R2 RMSE MAE MD R2 RMSE 模型1 0.170 9 -0.027 3 0.760 1 0.277 6 5.258 3 -0.842 7 0.853 1 8.306 9 模型2 0.197 0 0.026 3 0.781 5 0.264 9 5.156 0 0.337 4 0.855 8 8.232 2 模型3 0.177 8 -0.011 3 0.803 3 0.251 3 5.161 4 0.007 4 0.856 0 8.226 2 模型4 0.177 4 -0.010 2 0.803 9 0.250 9 5.161 2 -0.000 4 0.856 0 8.226 1 Table 3. Validation results of different crown models
-
模型拟合的总体评价反映了总体树冠外轮廓和体积的变化,不能反映各模型是否存在异方差性和无偏性,评价这2个指标最直观的方法就是利用残差分布图。为全面评价模型4效果,分别绘制模型4加权前后的残差分布图。图 2为未增加权函数时相容模型的残差图,左侧为树冠外轮廓模型残差分布,右侧为树冠体积模型残差分布,从图 2可以看出:残差均存在明显喇叭口形状,说明异方差问题显著。图 3为增加权函数后相容模型残差分布图,从图 3可以看出:加权后模型的残差散点图分布变得均匀,说明权函数明显消除了异方差。进一步说明加权后模型4显示了较高的等方差性和无偏性,效果较好。