学术 | 基于无人机航测数据的森林郁闭度和蓄积量估测

 提示：点击上方"达北云"↑免费订阅本刊

森林

森林是地球生态系统中重要的组成部分，其碳汇功能在维持全球碳平衡方面发挥着重要作用。森林蓄积量是生物量和碳储量研究的重要参考因素，是评价森林固碳能力的重要指标，也是森林资源调查的重要因素。
传统的一、二类森林资源调查是测算森林蓄积量的重要方法，这两种方法虽然准确度较高，但在时间与空间上存在较大的局限性，调查过程会耗费大量的时间和人力。

如何在准确获取森林蓄积量信息的同时，快速、高效地为区域生态状况评估提供参考已成为研究热点。

近年来，随着遥感技术的发展，基于遥感因子的森林蓄积量估测方法因其可以快速、高效获得蓄积量模型而引起了国内外学者的广泛关注。回归建模估测法因其工作量相对较少、估测频率高、估测覆盖面积大等优点逐步成为基于遥感因子蓄积量估测研究的主要方向。
按照自变量数量可以将蓄积量回归模型分为一元和多元回归，按照方程表现形式可以分为线性和非线性回归。研究森林蓄积量的线性回归方法主要有最小二乘法、逐步回归法、主成分回归法、偏最小二乘法。其中，偏最小二乘法结合了多种回归方法的功能和优点，有效地改善了模型自变量间的多重共线性问题，目前已广泛地应用于各领域的估测研究中。
带有点云数据的航测数据与其他遥感数据相比，具有分辨率高、细节更完善、时效性强、成本低等优点，可以更准确地提取蓄积量的相关特征因子，从而提高蓄积量估测精度。

东北林业大学苏迪，高心丹利用无人机航摄影像点云数据估测树高和胸径因子，使用正射影像通过分水岭算法提取树冠个数，结合提取的坡向等因子，用主成分回归方法估测郁闭度，结合全部因子通过偏最小二乘回归方法估测森林蓄积量。

1 研究区域与数据

1.1　研究区域

1.2　研究数据

1.3　样地筛选

2 特征参数提取与结果

利用GIS技术提取与森林蓄积量相关的树冠个数、平均树高、平均胸径、坡度、坡向、海拔、小班面积、郁闭度这8个GIS因子建立蓄积量模型，通过该方法提高蓄积量的估测精度。本研究先使用ArvGIS10.2软件将DOM图像和点云数据图像按照小班边界进行裁剪，为进行特征因子提取做准备。本研究结合点云数据、DOM图像和研究区小班矢量数据，以二类调查小班为单位对树冠个数、平均树高等因子进行估测，对海拔、坡向等因子进行提取，获得建立蓄积量模型所需的特征因子，特征因子提取流程图见图3。

2.1　特征参数的提取

2.1.1　分水岭树冠个数提取

分水岭算法适合提取边缘微弱较敏感、分辨率高、细节多的影像轮廓，但传统分水岭算法存在过分割问题，为了解决该问题，本研究将120个小班正射影像作为分水岭分割的原始图像，首先使用支持向量机的分类方法将其分为针叶林和阔叶林；其次采用最大熵法分别获得针叶林和阔叶林的最优阈值，比较针、阔叶林的最大熵最优阈值，发现针叶林最优阈值大于阔叶林最优阈值；最后对形态学开闭重建处理后的针阔混交林图像进行基于阈值标记的分水岭分割，第1次将小于阔叶林阈值区域进行标记和分水岭分割，第2次将大于阔叶林最优阈值，小于针叶林最优阈值进行分割得到树冠轮廓图像。

2.1.2　树高和胸径的估测

使用LIDAR360工具对点云数据进行去噪处理，生成0.5m分辨率的DEM(digital elevation model， 数字高程模型） 和DSM(digital surface model，数字地表模型）。由于DEM模型中只有地形的高程信息，不包括DSM模型中的森林树木高度等地表信息的高程，因此从DSM中减去DEM即可获得CHM(canopy height model， 冠层高度模型）。CHM中每个点的像素值即为该点的高程值，通过MATLAB 2013b软件编程，读取各个小班样本的CHM模型平均高程，只保留高程值不为0的点进行平均高程的计算。由于CHM为冠层模型不能等同于树高，因此需要结合读取的CHM平均高程利用一元线性回归方法建立平均树高模型和平均胸径模型。

2.1.3　海拔、坡向、坡度、小班面积因子的提取

使用LIDAR360软件在DEM基础上生成等高线图（图4a）。其中比例尺为1∶1000、三角形最大边为30m、间曲线等高距为2.5m、首曲线等高距为5m、计曲线等高线为25m。利用DEM数据生成坡向图像和以10°为间隔进行处理生成坡度图像，如图4b、c所示。使用Python 2.0编程，用双线性插值法提取等高线、坡向和坡度平均值。使用ArcGIS10.2软件结合小班矢量图计算小班面积，提取各项因子为接下来建立郁闭度和蓄积量模型提供特征因子。

2.1.4　郁闭度估测

主成分回归是主成分分析和回归估计结合的方法，主成分分析是一种变量降维技术，是将多个相关变量重新组合成几个主成分变量，减少变量个数的变量分析方法。该方法能够简化分析过程，提高分析效率。同理主成分回归是用合成后的主成分变量进行回归分析，不仅保留了原自变量对因变量的影响，还消除了各个变量相互间的共线性问题。本研究通过主成分回归将树冠个数、平均树高、平均胸径、坡度、坡向、海拔、小班面积这7个特征因子组成若干相互独立的主成分，根据各个主成分的贡献率大小选取模型最终的主成分变量，建立郁闭度回归模型。使用MATLAB 2013b编程，分析1～70号小班数据，得出了各主成分的特征根、主成分得分、主成分累积贡献率。按照贡献率累积80％以上的选取标准，根据贡献率的大小顺序选取主成分，主成分与特征因子间的关系方程，再根据主成分的构成，还原郁闭度方程，得到郁闭度与７个特征因子的关系模型。

2.2　特征参数的提取结果

2.2.1　树冠个数提取结果实验使用

2.1.1部分改进后分水岭分割算法提取树冠轮廓信息，将树冠分割结果图像进行二值化处理计算连通区域得到树冠个数，通过公式（2）计算分割准确率和平均分割准确率来显示分割效果。

式中： A_d为准确率；A_d为平均准确率；N_c为分割的树冠个数；N_d为小班树木株数总数（实测小班株数）；n为图像个数。通过计算得到120个样地DOM图像的传统分水岭平均分割精度为52.89％，改进后分水岭平均分割精度为80.03％，能够达到基本分割要求，改善了传统分水岭存在的过分割现象，分割后得到的树冠个数可以用于建立蓄积量模型。经计算，23号样地树冠分割精度与平均分割精度最相近，所以选择该样地图像进行显示（图5）。

2.2.2　平均树高、平均胸径估测结果

根据2.1.2部分得到平均树高和平均胸径模型结果如下：

式中:R_DBH为估算的平均胸径；H_CHM为CHM提取小班的平均高程。平均树高模型（公式3）精度达到97.34％，平均胸径模型（公式4）精度达到91.27％。用71～120号这50个检验样本对平均树高和平均胸径的估测值与实测值建立一元线性回归方程，50个检验样本的平均树高和平均胸径估测值与实测值的拟合曲线见图6。由图6和平均树高、平均胸径模型的平均精度可知，估测值与实测值的拟合程度较好，且都是由CHM直接进行估测而成，充分利用点云数据所生成的CHM数据。

2.2.3　郁闭度估测结果

由2.1.4部分将7个变量进行主成分分析，其中第1主成分贡献率32.087％，第2主成分贡献率24.375％，第3主成分贡献率18.005％，第4主成分贡献率13.769％。前4个主成分贡献率已达到88.235％，符合主成分选取要求，选取前4个主成分因子为新变量。则公式（5）～（8）为4个主成分与特征因子间的关系方程。

方程（9）和（10）中X₁为郁闭度，通过71～120号50个检验数据进行精度检验，郁闭度模型（公式10）精度为83.18％，可以达到精度要求，再结合前文所有提取和估测的特征因子建立蓄积量模型，进一步研究各个特征因子与蓄积量之间的关系。

3 蓄积量估测模型

3.1　偏最小二乘回归方法

偏最小二乘回归方法是一种集典型相关分析和主成分分析优点以及功能于一身的多元线性分析方法。该方法由于精度高，稳健性实用性好而被广泛地应用于各种估测研究中。设已知单个因变量Y和自变量X＝［x₁，x₂，…，x_p］，样本个数n，偏最小二乘法从X和Y的相关矩阵中提取主成分t，使用Y和X对t进行回归，对于回归分析的需要，偏最小二乘回归在提取成分时添加以下两个目标：１）尽可能多地在t中携带X矩阵数据表中的变异信息；２）t和Y的相关程度能够达到最大。在对X和Y进行一次成分提取后，分别进行X对t的回归和Y对t的回归。如果回归方程达到了满意的精度，就停止计算，否则将X关于t回归后的残差矩阵与Y关于t回归后的残差矩阵进行新一轮的成分提取。反复重复这个过程，直至满足交叉有效性原则规定的条件，最终确定提取主变量成分的个数，建立偏最小二乘回归方程。

3.2　蓄积量模型的建立

借助SPSS软件，通过相关性分析得到与蓄积量存在显著关系的8个自变量如表1所示，表中X1—X8同前文。表1表明各因子间存在多重相关性，一般回归模型已不适用，偏最小二乘法对自变量的选择门槛不高，不需要选择最优因子，且较多变量有益于对提取的主成分进行累计解释分析，因此本研究采用偏最小二乘法建立森林蓄积量估测模型。在MATLAB 2013b软件中按照偏最小二乘原理编写PLS回归代码，对1～120号样本用偏最小二乘法建立蓄积量模型，对模型进行交叉有效检验，得到原始变量的回归方程系数，据此可以拟合出基于偏最小二乘回归方法蓄积量方程。

式中：V为蓄积量，方程包含的1—X8与上文一致，将30个固定样地检验样本数据代入方程（公式11），得到蓄积量估测值，为模型评价和精度检验做准备。

3.3　模型评价和精度检验

结合预留的121～150号30个验证样本实测数据和蓄积量方程（公式11）估算的森林蓄积量得到图7实测值与估测值折线图。

从图7中可以看出估测的蓄积量值整体偏大，通过对蓄积量实测值和估测值进行配对T检验，得到T＝-0.374，P＝0.657＞0.05，说明实测值和估测值相比虽然整体偏小，但不存在显著性差异。使用实测值与估测值对模型进行评价，根据模型的决定系数R²，均方根误差R_MSE，总体相对误差RS，平均相对误差M_RE，预估精度P_r为评价指标。根据公式（12）～（16）计算：

式中：y_i为蓄积量实测值；＾y_i为蓄积量模型预测值；i表示第i个样本。本研究蓄积量模型（公式11）的R²＝0.7386，R_MSE＝5.1353m³/hm²，说明模型拟合效果较好。利用模型偏差统计量进行比较并评价模型的预测能力，R_S＝-1/2858、M_RE＝ -0.2630，说明估测蓄积量与实测蓄积量偏离不大，可以很好地进行蓄积量估测；精度达到88.43％，估测能力较好，并且有一定的可推广性。以蓄积量实测值为横坐标，估测值为纵坐标，建立实测值和估测值散点图，建立的散点图趋势线拟合效果较好，实测值和估测值比较吻合（图8）。

结论

利用无人机影像生成的点云数据和正射影像估测蓄积量的方法，能够在充分利用无人机影像的同时最大限度地保留无人机影像的细节和各种因子的特征，提高特征参数的提取精度。该方法不仅提高了实施效率、降低了估测成本，还为代替人工野外实测森林调查信息提供了可能。无人机数据生成的正射影像与高分辨率卫星影像相比具有重叠度大、分辨率高等优点；其生成的点云数据不仅可以代替机载雷达获得高程信息，还具有航线灵活、操作方便、时效性好等优点。无人机航测数据弥补了高分辨率卫星影像无法同时获得影像和高程信息的缺点，在估测森林生物量和蓄积量方面应用前景广阔。本研究利用无人机航摄方法采集了黑龙江帽儿山林场老山实验区航摄影像，利用三维点云数据提取CHM模型高程，估测了平均树高和平均胸径；用改进的分水岭算法获得了树冠个数；用DEM数据获得坡度等特征信息，估测了森林蓄积量。使用帽儿山林场固定样地数据作为检验样本，结果表明，平均树高提取精度达到97.34％，平均胸径提取精度达到91.27％，树冠个数提取精度达到80.03％，郁闭度模型精度为83.18％，森林蓄积量估测精度达到88.43％，满足森林调查精度要求。研究得到如下结论：１）在提取CHM平均高程值时不计算高程值为0的点，减少了点云数据空白点对平均树高、胸径模型的影响；２）通过遥感技术提取特征因子，用主成分回归的方法可有效地估测郁闭度；３）模型中所有因子都是通过遥感与回归方法提取或估测生成，使用偏最小二乘回归方法建立蓄积量模型，为信息化、自动化提取特征因子估测蓄积量提出可能。由于时间和试验条件限制，仍存在很多不足和改进的地方，如波段信息对蓄积量估测的影响等方面有待进一步研究。

该文发表于《林业工程学报》2020年第1期。

往期推荐