基于 UAV 和机器学习方法的甘南地区高寒草 地地上生物量遥感估测研究

杨利霞

基于 UAV 和机器学习方法的甘南地区高寒草

地地上生物量遥感估测研究

& H6 `4 b0 p; P% ^; B  F

1 q. O- t5 |  @# o6 e( o' ]) W草地地上生物量（Above Ground Biomass, AGB）是区域碳循环分析和草地
. c- @: M& _4 D2 O9 o5 j2 N6 ?生态系统健康评价、草地资源可持续开发利用的关键指标。在青藏高原地区，受
特殊气候和复杂地形条件的影响，生态环境极其脆弱，草地生物量能够反映高原
地区草地生态系统和植被物种的动态变化情况，是全球气候变化研究的重要指标
之一。因此，构建高精度的草地生物量估测模型对草地资源管理、草畜平衡分析、
草地生长状况评估和生态环境保护均具有重要的意义。
许多学者针对青藏高原地区的草地生物量估测已开展了大量的研究。然而，
现有草地生物量估测模型差异较大，适合高寒草地生物量估测的反演方法尚缺乏
# D5 d2 C) O' ~" o) ~系统全面的研究。因此，本研究以青藏高原东缘的甘南地区作为典型研究区，利
* X6 o" ?7 z9 @  [/ L' _用与草地生物量有显著相关性的地理位置、地形、植被生物物理指标、气象、土
壤和遥感植被指数等因子，分别构建了草地生物量单因素参数模型、多因素参数
模型和非参数模型，利用 10-fold 交叉验证的方法系统比较了各种模型的精度和
稳定性。在此基础上，分析了 2000-2016 年研究区草地生物量的时空动态变化特
征。研究结果表明：
1）对 MOD13Q1 植被指数数据的滤波去噪处理是降低基于 MODIS 数据高
寒草甸草地 AGB 反演误差的关键之一。基于 Gaussian 函数、Double-Logistic 曲
6 A. O0 v3 l6 Q5 D% o3 f! U8 M线和Savitzky-Golay三种方法进行滤波处理后的EVI比MOD13Q1 EVI最优AGB
0 Z; W; n0 j' a9 b$ l4 d7 g0 x估测模型的精度高，其测试集的 R2 增加 0.001~0.009，RMSE 减少 0.72 kg/ha~1.33
# _9 @2 u* O. ?' Y. j$ kkg/ha；基于以上三种方法进行滤波处理后的 NDVI 比 MOD13Q1 NDVI 最优 AGB
估测模型的精度高，其测试集的 R2 增加 0.057~0.068，RMSE 减小 19 kg/ha~34.21
& r! u; I( [" Okg/ha。
2）在地理位置、地形、植被生物物理指标、土壤和气象等 14 种因子中，经
3 v5 |. n0 W4 M4 h度（X）、纬度（Y）、海拔（DEM）、坡位（TPI）、草层高度（H）、草地盖度（C）、
8 e9 D" q7 a- `/ V& i1 Y气温（T）、降水（P）和底层土壤粘粒含量（Clay2）等 9 种因子和草地生物量之
- S3 D8 X5 C/ I7 m' c6 M  Q! O间具有显著的相关性。在所有单因子最优估测模型中，基于草地盖度的指数模型
: C2 ~3 D/ j; u8 Q# Y精度最高，其测试集的 R2 和 RMSE 分别为 0.39 和 859.368 kg/ha；其它因子最优
' a/ `' B6 {) P; u1 H5 p* v9 `+ uI II
* {1 p9 ?7 }2 ~! g; }- e估测模型的精度由高到低依次为 H、P、T、Y、DEM、X、TPI 和 Clay2。总体
来看，9 个单因子 AGB 参数模型的估测精度较低，存在较大的不确定性，仅能
8 J& Z- c! W. t  T: u解释生长季草地生物量变化的 1.5%~36.3%。
3）与单因素 AGB 模型相比，构建草地生物量多因素参数模型和非参数模
5 Q% l1 X; V3 Z4 O2 l0 q型是提高草地生物量监测精度的重要途径。在 6 大类因素中，基于植被生物物理
指标 H 和 C 构建的最优估测模型在参数和非参数两类模型中的精度均最高，其
次依次为气象（T、P）、地理位置（X、Y）和地形（DEM、TPI）等因素构建的
  K6 a  e0 W7 f" y  ]( X9 V最优草地生物量估测模型；并且以上 4 类因素所构建的非参数模型精度较参数模
" r; D* v% ?: C4 h9 L2 ]" T型高，其测试集的 R2增加 0.07~0.1，RMSE 减小 1.41kg/ha~80.65kg/ha。在基于
% R% Z* l8 i; i. M各因素组合的参数和非参数模型中，不同组合模型的精度随着构建模型时自变量
  M9 W0 v8 @2 H5 ]/ B& E的减少，两类模型的精度均呈现逐渐递减的趋势。在相同组合中，非参数模型精
: g% _6 M( r3 A0 d% ^: g% [度均高于参数模型，其测试集的 R2增加 0.04~0.12，RMSE 减小 74.91 kg/ha~138.34
kg/ha。
# C% _/ e. B' y$ c9 N/ B# s5 e4）构建多因素参数和非参数模型是改进草地生物量估测模型稳定性的关键。
在参数和非参数 2 类模型中，随着参与构建模型因素数量的增加，模型的稳定性
也呈现增加趋势。在所有模型中，基于 X、Y、DEM、TPI、H、C、T、P、Clay2、
/ z8 [1 X; p4 L2 C- i  M% b- yEVISG等 10 种因子的最优非参数模型稳定性最高，其 RMSE 的标准偏差较最优
4 Q& U9 l1 w  R- n! \9 O单因素模型减少 51.643 kg/ha ~61.286 kg/ha，较最优多因素参数模型减少 37.508
; D3 Q* C' B: J$ L# [) Wkg/ha~61.286 kg/ha。
7 M+ ~4 u6 |$ e! V0 P, R5）2000-2016 年甘南地区草地生物量总体呈现增加趋势，年均增加量约为
12.463kg/ha。从草地生物量变化趋势的空间分布特征来看，研究区大多数区域的
9 R, M  ?/ {2 x' g1 h1 s草地生物量变化趋势处于稳定和增加的状态，其中保持稳定区域占整个研究区草
地面积的 78.53%，增加趋势的区域占整个研究区草地面积的 20.50%，而生物量
呈减少趋势的区域仅占 0.98%。从草地类型来看，除低平地草甸类草地生物量有
轻微减小趋势外，其它各类型的草地生物量变化主要处于稳定和增加趋势的状态，
" H4 a1 h6 m! L/ ]/ j- y8 z其中保持稳定的区域面积占各种草地类型的面积比例介于 76.81%~97.64%，呈现
3 E$ m9 e2 X3 Y6 _0 D: {. L# F! u增加趋势的比例介于 12.55%~23.02%，而具有减小趋势的比例仅为 0.17%~6.43%。
; o- H+ d$ V1 n2 v$ h+ K& e5 @$ {