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

权利要求书

1.一种基于无效值修正的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时序数据。

2.根据权利要求1所述的一种基于无效值修正的NDVI数据重建方法，其特征在于，步 骤二中所述时间序列上最近邻的不含无效数据影像为含无效数据的影像的下一个或前一个时 间的无云影像，其中，采用的无云影像为下一个时间的无云影像ImgNext或前一个时间的无 云影像ImgPre的判断方法如下：

1)获得和含无效数据的影像Img时间上最邻近的前一个时间的无云影像ImgPre和下一 个时间的无云影像ImgNext；

2)查找Img上的无效数据区A区和环状缓冲区B区，然后分别计算ImgPre和ImgNext 上对应的环状缓冲区B’和B”；

3)取出B’和B”中的像素，分别计算与B区中像素之间的相关系数R’和R”；

4)若R’>R”，则采用ImgPre对Img进行修复，否则，采用ImgNext对Img进行修复。

3.根据权利要求1所述的一种基于无效值修正的NDVI数据重建方法，其特征在于，步 骤三中所述一定长度的水平线状结构元的长度为5。

4.根据权利要求1所述的一种基于无效值修正的NDVI数据重建方法，其特征在于，步 骤一中所述无效值包括：(1)由于地表反射率数据缺失，导致计算出的光谱指数值无效，产 生的缺失无效值；(2)由于冰雪覆盖，光谱指数不能正确反应应有地物信息，产生的雪覆盖 无效值；(3)由于云遮蔽，光谱指数不能正确反应地物信息，产生的云遮蔽无效值。