欢迎访问《应用生态学报》官方网站,今天是 分享到:

应用生态学报 ›› 2023, Vol. 34 ›› Issue (7): 1806-1816.doi: 10.13287/j.1001-9332.202307.021

• 研究论文 • 上一篇    下一篇

基于多源遥感数据协同的滇西北森林郁闭度估测

周文武, 舒清态*, 王书伟, 杨正道, 罗绍龙, 胥丽, 肖劲楠   

  1. 西南林业大学林学院, 昆明 650224
  • 收稿日期:2022-12-07 接受日期:2023-05-09 出版日期:2023-07-15 发布日期:2024-01-15
  • 通讯作者: *E-mail: shuqt@163.com
  • 作者简介:周文武, 男, 1996年生, 硕士研究生。主要从事3S 技术在林业中的应用研究。E-mail: 1723074419@qq.com
  • 基金资助:
    云南省农业联合专项-重点项目(202301BD070001-002)、国家自然科学基金项目(31860205,31460194)和云南省教育厅科学研究基金项目(2023Y0728)

Estimation of forest canopy closure in northwest Yunnan based on multi-source remote sensing data colla-boration

ZHOU Wenwu, SHU Qingtai*, WANG Shuwei, YANG Zhengdao, LUO Shaolong, XU Li, XIAO Jinnan   

  1. College of Forestry, Southwest Forestry University, Kunming 650224, China
  • Received:2022-12-07 Accepted:2023-05-09 Online:2023-07-15 Published:2024-01-15

摘要: 森林郁闭度(FCC)是评价森林资源和生物多样性的重要参数,利用多源遥感协同手段以较小成本高精度实现区域FCC反演是当前研究热点。本研究以星载激光雷达ICESat-2/ATLAS为主要信息源,结合54块实测样地数据,采用贝叶斯优化(BO)算法改进后的随机森林(RF)、K-最近邻值法(KNN)、梯度回归(GBRT)模型获取光斑尺度ATLAS光斑内FCC,协同多源遥感影像Sentinel-1/2及地形因子基于BO算法优化后的全连接深度神经网络模型(DNN)进行区域尺度的滇西北香格里拉市FCC遥感估测。结果表明: 在提取的50个ATLAS激光雷达光斑参数指标中,经RF特征变量优选后, 6个特征参数(乔木冠层百分比、冠层顶光子相对高度的标准差、冠层高度最小值、区段内98%冠层高度值与冠层高度中位数的差值、冠层顶部光子数、表观反射率)贡献率较大,可作为光斑尺度遥感估测模型变量。在BO-RF、BO-KNN、BO-GBRT模型中,以BO-GBRT模型估测的FCC结果最优,留一交叉验证的决定系数(R2)为0.65、均方根误差(RMSE)为0.10、绝对残差均值(RS)为0.079,预测精度(P)为0.792,可作为研究区有林地74808个ATLAS光斑的FCC估测模型。以有林地ATLAS光斑尺度FCC值作为区域尺度BO-DNN模型的大样本数据,联合多源遥感因子进行研究区FCC估测,十折交叉验证的BO-DNN模型验证精度为R2=0.47、RMSE=0.22、P=0.558。使用BO-DNN模型估测及普通克里格(OK)插值的研究区FCC均值分别为0.46、0.52,主要分布在0.3~0.6,分别占比77.8%、81.4%。直接通过OK插值方法获取FCC效率较高(R2=0.26),但预测精度明显低于BO-DNN模型(R2=0.49)。FCC高值区域在研究区由西北向东南贯穿分布,北部地区和东南部分别为FCC高值、低值主要分布区。基于ICESat-2/ATLAS高密集光斑进行山地FCC估测具有一定优势,以光斑尺度的小样本数据估测结果可作为区域尺度深度学习模型的大样本数据,能为光斑尺度上推至区域尺度低成本、高精度估测FCC提供一种参考。

关键词: 深度学习, ICESat-2/ATLAS, 贝叶斯优化算法, 多源遥感数据, Sentinel数据, 普通克里格插值

Abstract: Forest canopy closure (FCC) is an important parameter to evaluate forest resources and biodiversity. Using multi-source remote sensing collaborative means to achieve regional forest canopy closure inversion with low cost and high-precision is a research hotspot. Taking ICESat-2/ATLAS data as the main information source and combined with data of 54 measured plots, we estimated FCC value by the Bayesian optimization (BO) algorithm improved random forest (RF), K-nearest neighbor (KNN), and gradient boosting regression tree (GBRT) model at footprint-scale. Combined with multi-source remote sensing image Sentinel-1/2 and terrain factors, we estimated the regional-scale FCC value of Shangri-La in the northwest Yunnan based on deep neural network (DNN) optimized by BO algorithm. The results showed that six characteristic parameters (percentage of tree canopy, standard deviation of relative height of photons at the top of the canopy, minimum canopy height, difference between 98% canopy height and median canopy height in the segment, number of top canopy photons, apparent surface reflectance) out of the 50 parameters that were extracted from ATLAS lidar footprint had higher contribution rate after RF characteristic variable optimization, which could be used as model variable for footprint-scale remote sensing estimation. Among BO-RF, BO-KNN, and BO-GBRT models, the FCC results estimated by the BO-GBRT model were the best at footprint-scale. The coefficient of determination (R2) was 0.65, the root mean square error (RMSE) was 0.10, the mean absolute residual (RS) was 0.079, and the prediction accuracy (P) was 0.792 for leave-one-out cross validation. It could be used as the FCC estimation model of 74808 ATLAS footprints for forest in the study area. We used the ATLAS footprint-scale FCC value of forest as the large sample data of the regional-scale BO-DNN model and combined with multi-source remote sensing factors to estimate FCC in the study area, the accuracy of the 10-fold cross-validation BO-DNN model was R2=0.47, RMSE=0.22, P=0.558. The mean values of FCC in the study area estimated by BO-DNN model and ordinary Kriging (OK) interpolation were 0.46 and 0.52, respectively, and the values mainly distributed in 0.3-0.6, accounting for 77.8% and 81.4%, respectively. The FCC efficiency obtained directly by the OK interpolation method was higher (R2=0.26), but the prediction accuracy was significantly lower than the BO-DNN model (R2=0.49). The FCC high value was distributed from northwest to southeast in the study area, and the northern and southeastern regions were the main distribution areas of high and low FCC values, respectively. It had certain advantages to estimate mountain area FCC based on ICESat-2/ATLAS high-density footprint, and the estimation results of small sample data at footprint-scale could be used as large sample data of deep learning model at region-scale, which would provide a reference for the low-cost and high-precision to FCC estimation on the footprint-scale up to the extrapolated regional-scale.

Key words: deep learning, ICESat-2/ATLAS, Bayesian optimization algorithm, multi-source remote sensing data, Sentinel data, ordinary Kriging interpolation