买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：重庆邮电大学

摘要：本发明公开了一种基于时空Kriging改进的多云地区MODISNDVI时间序列重构方法，该方法包括：结合质量信息的预插值流程、引入叠加窗口筛选机制的时空Kriging插值、和基于质量权重的S‑G迭代滤波。首先提取MOD13Q1提供的像元可信度和产品可用性数据作为质量信息标注，可以通过预插值生成参考数据。其次根据NDVI时序数据建立时空随机场，叠加窗口筛选有效样本构建时空变异函数开展时空Kriging插值补充缺失的时空信息。最后可在此基础之上开展加权S‑G迭代滤波以重构高质量的MODISNDVI时间序列数据。该方法能够有效应对多云地区NDVI时序数据连续缺失问题，能为NDVI时间序列重构提供重要参考，在一定程度上拓宽了NDVI的应用场景。

主权项：1.基于时空Kriging改进的多云地区MODISNDVI时间序列重构方法，其特征在于，包括以下步骤：1对NDVI时间序列数据进行投影转换与图像裁剪预处理；2根据质量信息标注对时间序列进行预插值；3将预插值后的NDVI时间序列数据集视为一个时空随机场，引入叠加窗口机制，借助空间探针对其附近搜索窗口内的所有像元做多时段筛选，得到满足拟合时间和空间变异函数所需数量的有效样本点；4对数据进行STL时间序列分解以满足时空区域变化量二阶平稳的要求，将数据分解为趋势项T、季节项S和随机项R三部分并剔除季节项S，在此基础上开展变异函数构建与时空Kriging插值，在插值过后加入季节项S得到最终的插值结果；所述时空Kriging插值，以邻近时空有效像元提供参考信息，插值补充缺失像元，时空Kriging的插值公式定义如下： 其中si,ti代表某样本点i在时空随机场中的坐标，Zsi,ti表示该样本点的观测值，λi代表其权重系数，且需满足条件n为样本点总数，Z*s,t代表时空点s,t处的估计值；所述时空变异函数公式为： 其中s,t代表样本点在时空随机场中的位置，h＝hs,ht代表相应的时空距离，γhs,ht为变异函数，式中E、Var表示数学期望和方差运算；5进行基于质量权重的S-G迭代滤波，平滑时间序列曲线，得到高质量的NDVI时间序列曲线；所述S-G迭代滤波的基本公式为： 其中Yj+i，表示原始时间序列值和滤波结果，wi表示相应的S-G滤波权重，N为滤波器长度，m为半窗口宽度；S-G滤波权重中对于wi＝1的点，两次迭代中值均不改变，其余值按下式更新权重： 其中wj代表更新后的权值，yj代表原始值，代表拟合过后的值。

全文数据：

权利要求：

百度查询： 重庆邮电大学 基于时空Kriging改进的多云地区MODIS NDVI时间序列重构方法