[发明专利]一种基于立方体技术的Landsat遥感植被指数修复方法

说明书

本发明公开了一种基于立方体技术的Landsat遥感植被指数修复方法，该方法首先对每个自然月内的Landsat数据计算成NDVI产品，并以月为单位采用最大合成法，将计算得到的Landsat NDVI轨道产品在同一自然月内的全部数据合成为一期月产品，该产品分辨率为30m，简称L30产品；然后，使用MODIS 250米分辨率16天NDVI产品，以月为单位采用最大合成法，合成NDVI月产品，简称M250产品，最后使用M250数据对L30数据进行修补。本发明方法能够修复由于云造成的Landsat数据缺失问题，采用基于前后一定时间周期内的Landsat NDVI有效数据，配合MODIS NDVI产品进行无效区域的估计和修补。

技术领域

本发明属于环境如植被监测领域，具体涉及一种基于立方体技术的Landsat遥感植被指数修复方法。

背景技术

卫星遥感在进行地表植被监测上具有覆盖范围广的特点，其中NDVI产品被广泛用于评价植被健康状况与长势等信息。Landsat系列卫星遥感影像具有30m空间分辨率，具有红光与近红外波段，适合于对地表植被情况进行检测。Landsat卫星的重访周期为16天，每个月对同一地理位置一般可以观测两次。这两次观测难免会受到云的影响，而有云覆盖的区域是无法获取地表有效信息的。因此当制作逐月的Landsat NDVI序列时由于云造成的数据缺失十分严重，这个序列产品难以直接使用。

为了修复由于云造成的Landsat数据缺失问题，本发明方法采用基于前后一定时间周期内的Landsat NDVI有效数据，配合MODIS NDVI产品进行无效区域的估计和修补。

发明内容

针对上述问题，本发明提供一种基于立方体技术的Landsat遥感植被指数修复方法，该方法基于前后一定时间周期内的Landsat NDVI有效数据，配合MODIS NDVI产品进行无效区域的估计和修补。

因此，本发明提供的一种基于立方体技术的Landsat遥感植被指数修复方法，首先对每个自然月内的Landsat数据计算成NDVI产品，并以月为单位采用最大合成法，将计算得到的Landsat NDVI轨道产品在同一自然月内的全部数据合成为一期月产品，该产品分辨率为30m，简称L30产品；然后，使用MODIS 250米分辨率 16天NDVI产品，以月为单位采用最大合成法，合成NDVI月产品，简称M250产品，最后使用M250数据对L30数据进行修补。

本发明所述的方法依次包括Landsat轨道NDVI产品处理流程、Landsat NDVI原数据月合成产品处理流程、MODIS 250m NDVI 16天产品处理流程、MODIS 250m NDVI月合成产品处理流程、Landsat NDVI数据立方体修补处理流程；其中，

所述Landsat轨道NDVI产品处理流程如下：

（1）将原始Landsat T1轨道数据产品解压缩，形成geoTiff格式的波段数据和质量控制码文件；

（2）计算轨道NDVI产品；

（3）对NDVI产品进行质量控制，对于云阴影和云象元使用无效值替换；

（4）对NDVI产品进行投影变换；

（5）NDVI产品以GeoTiff格式文件保存，产品信息入库；

所述Landsat NDVI原数据月合成产品处理流程如下：

（1）查询Landsat NDVI月产品数据库是否有某一月份产品；

（2）从数据库读取该月份全部轨道NDVI产品文件信息；

（3）进行数据拼接，重叠区域使用最大值合成；

（4）月合成产品保存为GeoTiff格式，产品信息入库；