[发明专利]一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法

权利要求书

1.一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法，其特征在于，包括如下步骤：

步骤1，基于现有叶面积指数和地表覆盖分类产品，利用聚类分析和最小差规则创建可代表全球主要土地覆盖类型的训练样本；

步骤2，通过训练长短时记忆LSTM、门控递归单元GRU、和双向LSTM深度学习模型，确定最佳叶面积指数估算模型为双向LSTM深度学习模型，即BiLSTM模型；

步骤3，通过将BiLSTM模型分别应用于MODIS 500米和250米地表反射率，生成分辨率为500米和250米的叶面积指数中间产品；

步骤4，利用时空加权平均后处理合并250米和500米分辨率LAI中间产品，从而得到全球250米分辨率时空连续的叶面积指数产品；

步骤4的具体实现方式如下；

以任意一年为例，计算该年前一年及该年的LAI序列1，以及该年及后一年的LAI序列2，对这2个LAI序列乘上时间加权函数，并在该年处相加，得到该年500米分辨率的LAI时间序列，其中时间加权函数w设计为以下形式，其中t代表时间序列的步数：

由于500米LAI的估算模型精度更高，为保持500米LAI和250米LAI的一致性，对250米LAI中间产品采用空间加权，即在每500米像素内，四个250米LAI像元被标准化以匹配500米LAI值：

最终得到空间分辨率为250米，时间分辨率8天，时间跨度从2000年至今的全球时空连续叶面积指数。

2.如权利要求1所述的一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法，其特征在于：步骤1的具体实现方式如下，

首先，通过对2000年之后的任意一年的时间序列GLASS LAI曲线进行K均值聚类分析，识别出不同类型的LAI时间曲线，基于MODIS土地覆盖产品分别对每种土地覆盖类型进行聚类；

为确保每个聚类代表真实的LAI时间序列，利用三个LAI产品进行比较，即GLASS第五版、MODIS第六版和PROBA-V第一版，以生成代表时间序列的连续LAI样本，时间采样为K1年，每8天时间间隔，时间步为：K1*365/8向上取整；对于每个小类，选择这三个LAI产品中差异最小的像元作为该小类的代表像元；由于MODIS LAI比GLASS和PROBA-V LAI波动更大，在每个时间点，如果GLASS和PROBA-V LAI的差值小于一个单位，则取其平均值作为样本值，否则取MODIS、GLASS和PROBA-V LAI的中值；在此之后，在全球每个4°×4°的窗口中排查，如果窗口中没有选中的像元，则按最小差规则添加一个代表性像元；

选取代表性像元以后，提取其对应的MODIS时间序列地表反射率数据作为控制变量，融合后的时间序列LAI作为目标变量，组成训练样本；随机将样本分为三组，分别是用于获得深度学习模型的训练数据集、用于选择最优模型的验证数据集和用于评估最终模型的测试数据集。

3.如权利要求2所述的一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法，其特征在于：利用最小差规则来选择代表性像元的具体实现方式如下，

由于MODIS和PROBA-V时间序列LAI在时间上可能是不连续的，因此这三种产品的LAI序列的共同长度对于每个像素是不同的；这里应用一个最小差规则来选择代表性像元：当纬度低于50°，共有LAI时间序列长度不小于70，当纬度高于50°且低于55°，共有LAI序列长度不小于60，其他情况下共有长度不小于45。

4.如权利要求1所述的一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法，其特征在于：步骤2中，对于三个深度学习模型：LSTM、GRU和BiLSTM，保持相同的训练算法和参数：使用Adam优化器，初始学习率为0.0001，批次大小为100，最大训练期数为200。