滇西北森林郁闭度估测模型——基于全球生态系统动力学调查多波束激光雷达数据  

Construction of forest canopy closure estimation model in the northwestern Yunnan based on global ecosystem dynamics investigation multi⁃beam LiDAR data

作　　者：周文武 舒清态 胥丽 杨正道 高应群 吴再昆 夏翠芬 顾纯僖 李华 ZHOU Wenwu;SHU Qingtai;XU Li;YANG Zhengdao;GAO Yingqun;WU Zaikun;XIA Cuifen;GU Chunxi;LI Hua(College of Forestry,Southwest Forestry University,Kunming 650224,China)

机构地区：[1]西南林业大学林学院,昆明650224

出　　处：《生态学报》2024年第8期3525-3539,共15页Acta Ecologica Sinica

基　　金：国家自然科学基金项目(31860205,31460194);云南省农业联合专项⁃重点项目(202301BD070001⁃002);云南省教育厅科学研究基金项目(2023Y0728)。

摘　　要：探究全球生态系统动力学调查(GEDI)多波束激光雷达数据估测区域森林郁闭度(FCC)的潜力,对于评估森林生态系统状态和林分环境具有重要作用。选取滇西北典型生态脆弱区香格里拉为研究区,以GEDI波形数据为信息源,提取46245个有林地光斑参数,使用经验贝叶斯克里金法(EBK)获取光斑参数在研究区未知空间的连续分布,结合54块实测样地数据,采用支持向量机的递归特征消除法(SVM⁃RFE)、随机森林(RF)和Pearson分析分别优选特征变量,基于贝叶斯优化(BO)随机森林回归模型(BO⁃RFR)、贝叶斯优化梯度回归模型(BO⁃GBRT)和偏最小二乘法(PLSR)研建森林郁闭度最佳估测模型。结果表明:(1)EBK法预测精度高,估测结果可靠,R2:0.20—0.92,RMSE:0.004—2812.912,MAE:0.003—1996.258,MRE:0.007—4.423;(2)基于不同特征优选方法筛选的特征变量和数量略有差异,SVM⁃RFE法优选出6个参数(cover、pai、sensitivity、rv_a1、rv_a4、rg_a4)的平均交叉验证精度达0.84,RF法以贡献度5%为阈值筛选出5个参数(cover、pai、pgap_theta_error、modis_treecover、modis_nonvegetated),Pearson法以相关性大于0.3且在0.01水平显著优选出5个参数(cover、pai、rv_a5、rg_a5、pgap_theta_error);(3)不同特征变量优选方法筛选的建模参数研建估测模型精度差异性较大,以SVM⁃RFE和RF方法优选参数构建估测模型的精度更佳,SVM⁃RFE方法优选参数研建估测模型精度变化相对稳定,以RF方法中的BO⁃GBRT模型为最佳FCC估测模型(R2=0.85、RMSE=0.069,P=86.5%);(4)采用BO⁃GBRT模型估测研究区森林郁闭度和空间制图,与GEDI pai参数预测的FCC具有较高空间相关性达0.53,FCC均值分别为0.58、0.61,主要分布在0.4—0.7,分别占比65.45%、51.79%。研究区森林郁闭度主要处于中度郁闭,北部区域主要为高度郁闭区,与研究区植被覆盖度的空间分布具有一致性,说明使用GEDI数据估测森林郁闭度的方法具有可行性Exploring the potential of global ecosystem dynamics investigation(GEDI)multi⁃beam LiDAR data to estimate regional forest canopy closure(FCC)plays an important role in assessing forest ecosystem status and stand environment.The typically ecologically fragile area,Shangri⁃la,was selected as the study area in the northwestern Yunnan.GEDI waveform data was used as the information source to extract 46245 forest footprints parameters.The empirical Bayesian kriging(EBK)method was used to obtain the continuous distribution of footprints parameters in the unknown space of the study area.Then,combined with 54 measured samples data,the recursive feature elimination method of support vector machine(SVM⁃RFE),random forest(RF)and Pearson analysis were chosen to optimize the characteristic variables,respectively.The best estimation model of forest canopy closure was studied and constructed by Bayesian optimal random forest regression model(BO⁃RFR),Bayesian optimal gradient regression model(BO⁃GBRT),and partial least square method(PLSR).The results showed that:(1)The EBK method had high prediction accuracy and reliable estimation results.Its R2 was 0.20—0.92,RMSE was 0.004—2812.912,MAE was 0.003—1996.258,and MRE was 0.007—4.423.(2)There were slight differences in the method selection of characteristic variables and number based on different characteristic optimization methods.Among the three methods,the SVM⁃RFE method selected six parameters(cover,pai,sensitivity,rv_a1,rv_a4,rg_a4)with an average cross⁃validation accuracy of 0.84.The RF method selected five parameters(cover,pai,pgap_theta_error,modis_treecover,modis_nonvegetated)with a contribution of 5%as the threshold.The Pearson method significantly selected five parameters(cover,pai,rv_a5,rg_a5,pgap_theta_error)with a correlation greater than 0.3 and at the 0.01 level.(3)The modeling parameters selected by different characteristic variable optimization methods had great differences in the prediction accuracy of the estimation model.Among them,the accur

关 键 词：全球生态系统动力学调查(GEDI) 贝叶斯优化算法 机器学习 特征变量优选 经验贝叶斯克里金 大空间尺度 

分 类 号：S718.5[农业科学—林学]