[发明专利]一种基于无效值修正的NDVI数据重建方法

摘要

本发明公开了一种基于无效值修正的NDVI数据重建方法，由于该技术方案是在对NDVI产品不同“无效值”标记基础上，采用了基于时空最近邻线性回归方法进行“无效值”的修正，并引入了数学形态学来消除噪声，充分考虑了时间序列数据的时间相关性和空间相关性，克服了传统的单纯采用S-G滤波方法进行数据重建的不足。能够使得原始NDVI时间序列曲线的大量无效点被剔除，植被分布区误差值主要分布在0-0.2之间，大于0.2的数据很少，与传统S-G滤波方法相比，具有明显优势。

主权项

一种基于无效值修正的NDVI数据重建方法，其特征在于，包括如下步骤：步骤一，NDVI产品无效值标记；根据影像的红光和近红外波段的地表反射率数据，获得NDVI时间序列曲线，将NDVI时间序列曲线上的突变噪声点标记为无效点，其对应的NDVI值称为无效值；步骤二，基于时空最近邻线性回归方法修正无效值；通过建立含无效数据影像与其在时间序列上最近邻的不含无效数据影像之间同名区域灰度值之间关系,修复无效值数据；所述时间序列上最近邻的不含无效数据影像为含无效数据的影像的下一个或前一个时间的无云影像，所述无云影像对应区域全部是有效数据；其中：A区为含无效数据的影像的无效数据区，B区为A区外面的环状缓冲区，A’区为不含无效数据影像中与A区相对应的同名区域，B’区为不含无效数据影像中与环状缓冲区B区相对应的同名区域，B区和B’区中的像素都为有效像素；根据回归方程：g_B＝αg_B'+β首先，取出B区和B’区对应的像素进行回归计算，求得系数α和β，然后将A’区的像素逐点带入回归方程，计算出A区的像素值，即完成无效数据的修复；步骤三，采用结合数学形态学的S‑G滤波方法进行数据重建；利用数学形态学中的闭运算，采用一定长度的水平线状结构元对上述修复后的NDVI时间序列曲线进行闭运算处理，对处理后的NDVI时间序列曲线进行S‑G滤波处理，得到重建NDVI时序数据。