[发明专利]一种基于重构植被指数时间序列影像的湿地植被固碳速率遥感估算方法

权利要求书

1.一种基于重构植被指数时间序列影像的湿地植被固碳速率遥感估算方法，其特征在于，包括如下步骤：

S1、重构湿地植被分布区具有高时空连续性的增强型植被指数EVI影像时间序列；

S2、基于实测数据与同期对应像元的EVI观测值建立最优湿地植被碳存量遥感反演模型；具体为：以多种统计方法建立多个实测湿地植被碳存量Bio.C与同期对应坐标位置EVI值的关系模型，以决定系数/拟合优度R²以及均方根误差/效值RMSE为选择标准对比各个模型的拟合效果，选择R²最高且RMSE最低的模型作为最优湿地植被碳存量遥感反演模型；

S3、结合最优湿地植被碳存量遥感反演模型与重构的EVI时序影像重建湿地植被碳存量年内动态变化过程；

S4、根据累积存量法逐像元计算湿地植被碳存量年累积增长量，得到湿地植被年固碳速率。

2.如权利要求1所述的基于重构植被指数时间序列影像的湿地植被固碳速率遥感估算方法，其特征在于，步骤S1中，重构湿地植被分布区具有高时空连续性的增强型植被指数EVI影像时间序列具体为：以EVI影像的概括性质量图层对EVI时序曲线中的异常值进行筛选并剔除；通过S-G滤波对剔除的异常值进行插补，得到具有高时空连续性的重构EVI时间序列影像。

3.如权利要求2所述的基于重构植被指数时间序列影像的湿地植被固碳速率遥感估算方法，其特征在于，S-G滤波方法基本运算过程为：选取EVI影像点j附近m个临近点的值拟合d阶多项式，获得各点多项式系数C_i，以预测点j处的拟合值，S-G滤波对EVI时序数据的处理过程可用如下公式描述：

式中，EVI_j^*为EVI时序曲线拟合值，EVI_j+i为原始EVI时序曲线上的观测值，C_i为S-G滤波系数，m为滑动窗口宽度，N为卷积运算数目，N＝2m+1。

4.如权利要求1所述的基于重构植被指数时间序列影像的湿地植被固碳速率遥感估算方法，其特征在于，步骤S3中，结合最优湿地植被碳存量遥感反演模型与重构的EVI时序影像重建湿地植被碳存量年内动态变化过程具体为：以步骤S2中建立的最优湿地植被碳存量遥感反演模型，对S1中重构的高质量EVI时序影像进行逐幅计算，有效重建研究时段的湿地植被碳存量年动态变化过程。

5.如权利要求1所述的基于重构植被指数时间序列影像的湿地植被固碳速率遥感估算方法，其特征在于，步骤S4中，根据累积存量法逐像元计算湿地植被碳存量年累积增长量，得到湿地植被年固碳速率具体为：在重建湿地植被地表碳存量年内动态变化过程的基础上，以累积存量法逐像元计算湿地植被年固碳速率CSQ：

其中，CSQ指湿地植被年固碳速率；k_t是t时刻与t-1时刻间的湿地植被碳存量累积增长量；Bio.C_t是t时刻的湿地植被瞬时碳存量，n为EVI时序影像每年的获取影像数量。