可燃物载量间接估测模型

9.3.5　可燃物载量间接估测模型

为了寻找高精度的估测模型，本节介绍以林分平均年龄、平均胸径、郁闭度等为中间因子的云南松可燃物估测模型。

(1)以平均年龄及平均胸径为中间因子估测模型结果

云南松样地的平均胸径数据(因变量y)和TM相关因子及坡度、坡向、坡向指数、海拔等地形因子所对应灰度值(自变量)进行回归，得到云南松平均胸径的估测模型为:

y=11.933+0.300X₁-3.466X₂+0.353X₃-0.276X₄-0.065X₅+0.065X₆–0.114X₇+ 0.216X₈-0.653X₉+0.584X₁₀+0.040X₁₁-0.008X₁₂-0.021X₁₃-0.024X₁₄-0.206X₁₅-0.034X₁₆+0.462X₁₇-0.002X₁₈

式中y表示云南松平均胸径(cm)。X₁～X₁₈代表的变量见下表，F=1.227，F_0.05=1.69，F＜F_0.05，显著水平为0.05;双尾T检验显著性概率sig=2.52＞α_0.05=2.326，说明方程不显著;复相关系数为0.401，判定系数R²为0.161，调整后判定系数R²为0.030，模型效果不理想。

表9.26　模型参选自变量

同时，以平均树高为因变量和TM相关因子及地形因子所对应灰度值(自变量)进行回归，得F=1.06，F_0.01=1.69，F＜F_0.05，显著水平为0.05，双尾T检验显著性概率sig=3.57＞α_0.01=2.576，说明方程不显著;复相关系数为0.373，判定系数R²为0.139，调整后判定系数R²为0.013，模型效果也不理想。

(2)以平均郁闭度为中间因子估测模型结果

以郁闭度为因变量和TM相关因子及地形因子所对应灰度值(自变量)进行回归，回归方程如下:

y=0.569+0.727X₁-0.060X₂+0.003X₃+0.004X₄-0.002X₅-0.001X₆+0.002X₇-0.001X₈-0.001X₉-0.001X₁₀+0.001X₁₁-0.001X₁₅-0.002X₁₆-0.002X₁₇

式中y表示云南松郁闭度。X₁～X₁₈代表的变量和上表相同，但是由于坡度、坡向、坡向指数、海拔等变量系数非常小，可以忽略不计。F=52.644，F_0.01=2.08，F＞F_0.01，显著水平为0.01，双尾T检验显著性概率sig=0.00＜α_0.01=2.576，说明方程显著;复相关系数R为0.944，判定系数R²为0.892，调整后判定系数R²为0.875。

用30个样点的实测值和预测值对比图如下所示，其精度为88.3%。

图9.45　回归模型郁闭度预报值与实测值的比较

(3)建立胸径及年龄与可燃物负荷量的模型

由样地到林分平均年龄，并以林分平均年龄与林分平均胸径为自变量，拟合云南松林可燃物载量线性回归模型:

y=862.099-91.782X₁+33114.882X₂

式中y表示云南松可燃物载量(kg.hm^-2)。X₁代表年龄，X₂代表郁闭度。F=51.822，F_0.01=4.78，F＞F_0.01，显著水平为0.01，说明方程显著;复相关系数R为0.665，判定系数R²为0.442，调整后判定系数R²为0.433，双尾T检验显著性概率sig=0.00＜α_0.01=2.576，说明方程显著。

用30个样点的实测值和预测值对比图如下所示，其精度为65%左右。

图9.46　回归模型载量预报值与实测值的比较