一种基于立方体技术的Landsat遥感植被指数修复方法技术

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

A restoration method of Landsat remote sensing vegetation index based on cube technology

【技术实现步骤摘要】
一种基于立方体技术的Landsat遥感植被指数修复方法
本专利技术属于环境如植被监测领域，具体涉及一种基于立方体技术的Landsat遥感植被指数修复方法。
技术介绍
卫星遥感在进行地表植被监测上具有覆盖范围广的特点，其中NDVI产品被广泛用于评价植被健康状况与长势等信息。Landsat系列卫星遥感影像具有30m空间分辨率，具有红光与近红外波段，适合于对地表植被情况进行检测。Landsat卫星的重访周期为16天，每个月对同一地理位置一般可以观测两次。这两次观测难免会受到云的影响，而有云覆盖的区域是无法获取地表有效信息的。因此当制作逐月的LandsatNDVI序列时由于云造成的数据缺失十分严重，这个序列产品难以直接使用。为了修复由于云造成的Landsat数据缺失问题，本专利技术方法采用基于前后一定时间周期内的LandsatNDVI有效数据，配合MODISNDVI产品进行无效区域的估计和修补。
技术实现思路
针对上述问题，本专利技术提供一种基于立方体技术的Landsat遥感植被指数修复方法，该方法基于前后一定时间周期内的L本文档来自技高网...

【技术保护点】
1.一种基于立方体技术的Landsat遥感植被指数修复方法，包括下述步骤：/n该方法首先对每个自然月内的Landsat数据计算成NDVI产品，并以月为单位采用最大合成法，将计算得到的Landsat NDVI轨道产品在同一自然月内的全部数据合成为一期月产品，该产品分辨率为30m，简称L30产品；然后，使用MODIS 250米分辨率 16天NDVI产品，以月为单位采用最大合成法，合成NDVI月产品，简称M250产品，最后使用M250数据对L30数据进行修补。/n

【技术特征摘要】
1.一种基于立方体技术的Landsat遥感植被指数修复方法，包括下述步骤：
该方法首先对每个自然月内的Landsat数据计算成NDVI产品，并以月为单位采用最大合成法，将计算得到的LandsatNDVI轨道产品在同一自然月内的全部数据合成为一期月产品，该产品分辨率为30m，简称L30产品；然后，使用MODIS250米分辨率16天NDVI产品，以月为单位采用最大合成法，合成NDVI月产品，简称M250产品，最后使用M250数据对L30数据进行修补。


2.根据权利要求1所述的方法，其特征在于：
该方法依次包括Landsat轨道NDVI产品处理流程、LandsatNDVI原数据月合成产品处理流程、MODIS250mNDVI16天产品处理流程、MODIS250mNDVI月合成产品处理流程、LandsatNDVI数据立方体修补处理流程。


3.根据权利要求2所述的方法，其特征在于：
所述Landsat轨道NDVI产品处理流程如下：
（1）将原始LandsatT1轨道数据产品解压缩，形成geoTiff格式的波段数据和质量控制码文件；
（2）计算轨道NDVI产品；
（3）对NDVI产品进行质量控制，对于云阴影和云象元使用无效值替换；
（4）对NDVI产品进行投影变换；
（5）NDVI产品以GeoTiff格式文件保存，产品信息入库；
所述LandsatNDVI原数据月合成产品处理流程如下：
（1）查询LandsatNDVI月产品数据库是否有某一月份产品；
（2）从数据库读取该月份全部轨道NDVI产品文件信息；
（3）进行数据拼接，重叠区域使用最大值合成；
（4）月合成产品保存为GeoTiff格式，产品信息入库。


4.根据权利要求2所述的方法，其特征在于：
所述MODIS250mNDVI16天产品处理流程如下：
（1）获取MOD13Q1NDVI250m16天产品原始HDF分片数据文件；
（2）检查覆盖中国区的分片数据是否完整；
（3）对分片数据进行质量控制；
（4）将分片数据拼接到一个整个影像上面；
（5）投影变换为Albers投影坐标系；
（6）形成MODIS16天NDVI产品，产品信息入库。


5.根据权利要求2所述的方法，其特征在于：
所述MODIS250mNDVI月合成产品处理流程如下：
（1）以自然月为时间区间，获取某个自然月全部MODIS16天产品；
（2）使用最大值合成将16天产品合成为月产品；
（3）月产品保存为GeoTiff格式，产品信息入库。


6.根据权利要求2所述的方法，其特征在于：
所述LandsatNDVI数据立方体修补处理流程如下：
（1）查看数据获取没有进行过数据修补的Landsat月产品；
（2）针对该Landsat月产品进行逐像素判断是否无效值；
（3）对于无效值像素，查看同期MODIS月产品是否有效，如果MODIS月产品是有效值则进入数据立方体填充步骤，否则将新生成的结果影像该像素位置填充为无效值；
（4）获取前后各12期的Landsat和MODIS在该像素位置的NDVI值；
（5）找到距离当前月份最近的两组有效值值；
（6）通过MOD...

【专利技术属性】
技术研发人员：肖风劲，刘秋锋，王杰，王峰，黄大鹏，
申请(专利权)人：国家气候中心，
类型：发明
国别省市：北京;11