[发明专利]一种长时间序列的NDVI数据产品融合方法

说明书

本发明提供一种长时间序列的NDVI数据产品融合方法，能够得到有效的NDVI时间序列B的延长序列。所述方法包括：确定NDVI时间序列A的插值序列；确定第一时间段的插值序列与第一时间段的NDVI时间序列B的第一相关系数和第一混合距离，若第一相关系数没有通过显著性水平的检验或第一混合距离不小于预设的第一距离阈值，则进行基于傅里叶理论的变换与重构，重构插值序列；确定重构后的第一时间段的插值序列与第一时间段的NDVI时间序列B的第二相关系数和第二混合距离，若第二相关系数通过显著性水平的检验且第二混合距离小于预设的第二距离阈值，则将重构后的第二时间段的插值序列作为NDVI时间序列B的延长序列。本发明涉及遥感技术领域。

技术领域

本发明涉及遥感技术领域，特别是指一种长时间序列的NDVI数据产品融合方法。

背景技术

归一化植被指数(Normalized Difference Vegetation Index，NDVI)已经被广泛用于植被长势、物候和气象研究以及森林火灾等研究中。例如，就当前的NDVI数据产品来看，AVHRR GIMMS NDVI产品具有较长的时间序列，AVHRR GIMMS NDVI产品起始于1981年7月，而起始于1998年4月的SPOT VEGETATION NDVI数据的精度以及空间分辨率更高，为融合二者的优势，已有少数学者试图通过AVHRR GIMMS NDVI产品对SPOT VEGETATION NDVI数据进行延长。

现有技术中，仅线性插值法被用于延长SPOT VEGETATION NDVI数据，由于AVHRRGIMMS NDVI与SPOT VEGETATION NDVI的相关性具有显著的区域差异性，因此，单一的线性插值法限制了SPOT VEGETATION NDVI数据延长的有效性。

发明内容

本发明要解决的技术问题是提供一种长时间序列的NDVI数据产品融合方法，以解决现有技术所存在的单一的线性插值法限制了SPOT VEGETATION NDVI数据延长的有效性的问题。

为解决上述技术问题，本发明实施例提供一种长时间序列的NDVI数据产品融合方法，包括：

S101，获取NDVI时间序列A和NDVI时间序列B，其中，与NDVI时间序列A相比，NDVI时间序列B的精度以及空间分辨率更高，NDVI时间序列A与NDVI时间序列B具有的相同时间段称为第一时间段，NDVI时间序列A独有的时间段称为第二时间段；

S102，确定NDVI时间序列A的插值序列；

S103，确定第一时间段的插值序列与第一时间段的NDVI时间序列B的第一相关系数和第一混合距离，若第一相关系数没有通过显著性水平的检验或第一混合距离不小于预设的第一距离阈值，则利用插值序列与第一时间段的NDVI时间序列B进行基于傅里叶理论的变换与重构，重构插值序列；

S104，确定重构后的第一时间段的插值序列与第一时间段的NDVI时间序列B的第二相关系数和第二混合距离，若第二相关系数通过显著性水平的检验且第二混合距离小于预设的第二距离阈值，则将重构后的第二时间段的插值序列作为NDVI时间序列B的延长序列。

进一步地，S101中获取的原始NDVI时间序列A的空间分辨率的值大于原始NDVI时间序列B的空间分辨率的值；

S101中获取的原始NDVI时间序列A的时间分辨率的值大于原始NDVI时间序列B的时间分辨率的值。

进一步地，所述确定NDVI时间序列A的插值序列包括：

利用NDVI时间序列A与NDVI时间序列B的线性关系对NDVI时间序列A进行线性插值，得到NDVI时间序列A的插值序列。

进一步地，所述利用NDVI时间序列A与NDVI时间序列B的线性关系对NDVI时间序列A进行线性插值，得到NDVI时间序列A的插值序列包括：