一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法技术

本发明专利技术涉及一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法，包括如下步骤：步骤1，基于现有叶面积指数和地表覆盖分类产品，利用聚类分析和最小差规则创建可代表全球主要土地覆盖类型的训练样本；步骤2，通过训练长短时记忆、门控递归单元、和双向LSTM深度学习模型，确定最佳叶面积指数估算模型；步骤3，通过将BiLSTM模型分别应用于MODIS 500米和250米地表反射率，生成分辨率为500米和250米的叶面积指数中间产品；步骤4，利用时空加权平均后处理合并两个250米和500米分辨率LAI中间产品，从而得到全球250米分辨率时空连续的叶面积指数产品。本发明专利技术可以填补目前产品在高纬地区的空白，是目前唯一可满足全球气候观测系统模拟碳循环的叶面积指数产品。统模拟碳循环的叶面积指数产品。统模拟碳循环的叶面积指数产品。

【技术实现步骤摘要】
一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法


[0001]本专利技术属于定量遥感卫星产品生成领域，特别是涉及一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法。

技术介绍

[0002]叶面积指数(Leaf Area Index,LAI)是全球气候观测系统(GCOS)指定的陆地基本气候变量之一，广泛应用于陆地生态系统模型模拟、作物产量估算和植被变化监测等多种科学应用。卫星观测数据为LAI全球时间序列测绘提供了唯一可靠的手段。目前全球LAI产品存在着一些局限性。最突出的问题是，输入地表反射率往往受到云雾或高浓度气溶胶的污染，导致时间序列产品的波动或数据缺失。例如在多云有雪的地区，MODIS LAI数据缺失率可达40％，尽管GLASS
‑
LAI算法已经使用基于植被指数的平滑方法重建了地表反射率，但现有的V5产品在中高纬地区冬季仍存在异常值和假生长季的现象。此外，由于GLASS需要对地表反射率进行优化重建，导致其生产过程非常耗时。
[0003]其次，由于输入观测数据与反演算法的差异，现有的LAI产品存在显著差异。例如，在热带林区，不同产品之间的平均赖氨酸差异可达1个LAI标准单位。这会在植被变化分析和地表模型模拟中造成很大的不确定性。如何通过整合产品的差异，优势互补，仍是一个挑战。第三，全球气候观测系统(GCOS)需要250米分辨率的LAI产品用于碳建模(GCOS 2016)，然而，目前能够满足这种要求的LAI产品只是区域性的或短于一年的，如多伦多大学(UFFT)LAI产品。因此，目前还没有满足该要求的叶面积指数产品。
[0004]基于以后研究背景，目前迫切需要一种具有高计算效率、可生产高精度、高分辨率的长时间系列全球叶面积指数的新方法。

技术实现思路

[0005]本专利技术针对现有技术的不足，基于MODIS(Moderate Resolution Imaging Spectroradiometer，中等分辨率成像光谱仪)卫星21年长的地表观测数据，提供一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法。
[0006]为了达到上述目的，本专利技术提供的技术方案是：
[0007]步骤1，基于现有叶面积指数和地表覆盖分类产品，利用聚类分析和最小差规则创建可代表全球主要土地覆盖类型的训练样本；
[0008]步骤2，通过训练长短时记忆、门控递归单元、和双向LSTM深度学习模型，确定最佳叶面积指数估算模型；
[0009]步骤3，通过将BiLSTM模型分别应用于MODIS 500米和250米地表反射率，生成分辨率为500米和250米的叶面积指数中间产品；
[0010]步骤4，利用时空加权平均后处理合并两个250米和500米分辨率LAI中间产品，从而得到全球250米分辨率时空连续的叶面积指数产品。
[0011]进一步的，步骤1的具体实现方式如下，
[0012]首先，通过对2000年之后的任意一年的时间序列GLASS LAI曲线进行K均值聚类分析，识别出不同类型的LAI时间曲线，基于MODIS土地覆盖产品分别对每种土地覆盖类型进行聚类；
[0013]为确保每个聚类代表真实的LAI时间序列，利用三个LAI产品进行比较，即GLASS第五版、MODIS第六版和PROBA
‑
V第一版，以生成代表时间序列的连续LAI样本，时间采样为K1年，每8天时间间隔，时间步为：K1*365/8向上取整；对于每个小类，选择这三个LAI产品中差异最小的像元作为该小类的代表像元；由于MODIS LAI比GLASS和PROBA
‑
V LAI波动更大，在每个时间点，如果GLASE和PROBA
‑
V LAI的差值小于一个单位，则取其平均值作为样本值，否则取MODIS、GLASE和CGLS LAI的中值；在此之后，在全球每个4
°×4°
的窗口中排查，如果窗口中没有选中的像元，则按最小差规则添加一个代表性像元；
[0014]选取代表性像元以后，提取其对应的MODIS时间序列地表反射率数据作为控制变量，融合后的时间序列LAI作为目标变量，组成训练样本；随机将样本分为三组，分别是用于获得深度学习模型的训练数据集、用于选择最优模型的验证数据集和用于评估最终模型的测试数据集。
[0015]进一步的，利用最小差规则来选择代表性像元的具体实现方式如下，
[0016]由于MODIS和PROBA
‑
V时间序列LAI在时间上可能是不连续的，因此这三种产品的LAI序列的共同长度对于每个像素是不同的；这里应用一个最小差规则来选择代表性像元：当纬度低于50
°
，共有LAI时间序列长度不小于70，当纬度高于50
°
且低于55
°
，共有LAI序列长度不小于60，其他情况下共有长度不小于45。
[0017]进一步的，步骤2中，对于三个深度学习模型：LSTM、GRU和BiLSTM，保持相同的训练算法和参数：使用Adam优化器，初始学习率为0.0001，批次大小为100，最大训练期数为200。
[0018]进一步的，步骤2中500米LAI模型采用BiLSTM模型，利用MODIS 6个反射率波段以及3个太阳和卫星观测角作为特征序列；250m LAI模型除采用MODIS红和近红外反射率波段组合作为特征输入外，其余设置与上述500m LAI模型相同；其中MODIS 6个反射率波段中不含波段5，太阳和卫星观测角的时间长度为K1，间隔为8天。
[0019]进一步的，步骤3的具体实现方式如下；
[0020]利用步骤2中确定的BiLSTM模型和时间序列的MODIS地表反射率数据，估算出空间分辨率为500m、时间分辨率为8天的全球时间序列LAI，然后利用250m地表反射率，采用红光和近红外两波段BiLSTM模型，估算全球250m 8天LAI。
[0021]进一步的，步骤4的具体实现方式如下；
[0022]以任意一年为例，计算该年前一年及该年的LAI序列1，以及该年及后一年的LAI序列2，对这2个LAI序列乘上时间加权函数，并在该年处相加，得到该年500米分辨率的LAI时间序列，其中时间加权函数w设计为以下形式，其中t代表时间序列的步数：
[0023][0024]由于500米LAI的估算模型精度更高，为保持500米LAI和250米LAI的一致性，对250
米LAI中间产品采用空间加权，即在每500米像素内，四个250米LAI像元被标准化以匹配500米LAI值：
[0025][0026]最终得到空间分辨率为250米，时间分辨率8天，时间跨度从2000年至今的全球时空连续叶面积指数。
[0027]与现有技术相比，本专利技术的优点和有益效果如下：1.该方法可生产自2000年21年长时间序列的250米分辨率的时空连续叶面积指数，填补目前产品在高纬地区的空白；是目前唯一可满足全球气候观测系统模拟碳循环本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法，其特征在于，包括如下步骤：步骤1，基于现有叶面积指数和地表覆盖分类产品，利用聚类分析和最小差规则创建可代表全球主要土地覆盖类型的训练样本；步骤2，通过训练长短时记忆、门控递归单元、和双向LSTM深度学习模型，确定最佳叶面积指数估算模型；步骤3，通过将BiLSTM模型分别应用于MODIS 500米和250米地表反射率，生成分辨率为500米和250米的叶面积指数中间产品；步骤4，利用时空加权平均后处理合并两个250米和500米分辨率LAI中间产品，从而得到全球250米分辨率时空连续的叶面积指数产品。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波动更大，在每个时间点，如果GLASE和PROBA
‑
V LAI的差值小于一个单位，则取其平均值作为样本值，否则取MODIS、GLASE和CGLS LAI的中值；在此之后，在全球每个4
°×4°
的窗口中排查，如果窗口中没有选中的像元，则按最小差规则添加一个代表性像元；选取代表性像元以后，提取其对应的MODIS时间序列地表反射率数据作为控制变量，融合后的时间序列LAI作为目标变量，组成训练样本；随机将样本分为三组，分别是用于获得深度学习模型的训练数据集、用于选择最优模型的验证数据集和用于评估最终模型的测试数据集。3.如权利要求2所述的一种全球250米分辨率时空连续的叶面积指数卫星产品生成方法，其特征在于：利用最小差规则来选择代表性像元的具体实现方式如下，由于MODIS和PROBA
‑
V时间序列LAI...

【专利技术属性】
技术研发人员：马晗，
申请(专利权)人：武汉大学，
类型：发明
国别省市：