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

说明书

技术领域

本发明涉及一种NDVI数据重建方法，具体涉及一种基于无效值修正的NDVI数据重建 方法。

背景技术

MODIS(中分辨率成像光谱仪)是搭载在Terra和Aqua卫星上的一个重要的传感器，是 卫星上唯一将实时观测数据通过x波段向全世界直接广播，并可以免费接收数据并无偿使用 的星载仪器，全球许多国家和地区都在接收和使用MODIS数据。MODIS遥感资料在获取数 据过程中，由于传感器的因素，会导致某些像素的地表反射率数据缺失，同时，受云雾、冰 雪的影响，也会导致该区域无法获得真实的地表反射率，这些严重偏离真实数据的值称其为 “无效值”。“无效值”使得NDVI(植被归一化指数)指数的计算结果不能反映地表的真实 信息，从而导致指数时间序列失真，将严重影响森林扰动监测的结果。从NDVI数据集本身 出发设计降噪处理方法，主要有阈值法、各种数学滤波、拟合函数拟合重建法以及利用辅助 数据的数据拟合法。目前使用较多的S-G滤波(Savitzky-Golay)方法是一种无先验信息的滤 波，即S-G滤波时对时间序列上的噪声不考虑其产生原因，“一视同仁”进行处理，另外，用 S-G滤波方法在一定程度上可以对长时间序列NDVI数据时域噪声进行去除，但没有考虑地 物的空间相关关系，因此，缺乏一定的科学性和合理性，从而影响数据重建效果。

发明内容

针对目前常用的S-G滤波方法存在未针对噪声类型和产生的原因进行处理，以及滤波过 程中没有考虑地物空间相关关系的不足，提出依据“无效值”标记，识别不同类型噪声，采 用合适的修正方法，再利用数学形态学与S-G滤波相结合滤波，对MODIS产品数据进行重 建，在高保真度下，使重建后的数据能有效地消除时间序列曲线和图像上的异常噪声点，恢 复原本的数据值，获得植被生长的趋势表征，达到提高NDVI指数时间序列可靠性目的。

为实现上述目的，本发明采用的技术方案为：一种基于无效值修正的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时序数据。

进一步的，步骤二中所述时间序列上最近邻的不含无效数据影像为含缺失数据的影像的 下一个或前一个时间的无云影像，其中，采用的无云影像为下一个时间的无云影像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进行修复。

进一步的，步骤三中所述一定长度的水平线状结构元的长度为5。