MODIS系列之NDVI（MOD13Q1）九：单一土地类型NDVI及时序图（以耕地类型为例）

单一土地类型耕地NDVI及耕地时序图

在 MODIS系列之NDVI（MOD13Q1）五：NDVI处理流程 和 MODIS系列之NDVI（MOD13Q1）七：时间序列S-G滤波之Python  处理的都是地表各土地类型的NDVI综合。

 本次博文则以单一土地类型（以耕地类型为例）得出耕地的NDVI及时序图

一、所需数据

1.研究区域和时间（年份）的土地利用类型数据（又叫地表覆盖数据）下载。下载详见该系列博文：MODIS系列之NDVI（MOD13Q1）八：土地利用数据下载

本文以2010年河南省地表覆盖数据为例（.tif数据）

2.研究区矢量数据准备（常见矢量数据下载请见ArcMap 2 （常见矢量数据下载））

本文以河南省省界矢量数据为例（.shp数据）

二、ArcMap战场

1.将下载好的2010年河南省地表覆盖数据通过ArcMap打开。（注意：因为网站提供的地表覆盖数据受投影及条带的影响，不是按照行政区划分幅。所以下载的地表覆盖数据要包含研究区域。）

 2.打开ArcToolbox

若研究区域较大，需对几幅地表覆盖数据“合并”（非术语解释）。若研究区域较小，一副地表覆盖数据就能覆盖研究区域，请省略"合并"步骤。由于研究区为河南，故所需“合并”。

 3.找到“数据管理工具”下的“镶嵌”和“镶嵌至新栅格”

 “镶嵌”和“镶嵌至新栅格”根据需要，进行二选一。

镶嵌：将多个输入栅格镶嵌到现有栅格数据集。

 镶嵌至新栅格：将多个栅格数据集镶嵌到一个新的栅格数据集中。

 4.借助矢量数据，通过“裁剪”工具将“合并”后的地表覆盖数据进行裁剪

(注：此处用的“裁剪”工具为按掩膜提取。选择原因请参考ArcMap 3 “裁剪”)

 得到2010年河南省土地分类数据如图：

 下载的土地分类数据有空隙，始终解决不了。如果知道问题，请通过评论区告知，非常感谢。

5.先打开2010年河南省土地分类数据的属性表

 右击河南土地利用分类.tif。打开属性表。看到Value为10。

6.按属性提取。因为耕地分类为数字“10”。MODIS系列之NDVI（MOD13Q1）八：土地利用数据（又称地表覆盖数据）下载

 7.在点开的按属性提取页面，点击SQL栏

 8.进入“查询构建器”

 

 9.通过给出的按键选项，选出"Value" =10

 

 

 得到如图：

 10.将之前得到的综合NDVI.tif数据，通过按掩膜提取工具提取耕地分类数据。得到耕地NDVI.tif数据。

 由于通常所需数据量大，该步骤可以通过python批量提取。也可通过ArcToolbox中按掩膜提取的批处理（右击按掩膜提取）。

 

 

 11.得到各耕地NDVI.tif数据。后通过参考MODIS系列之NDVI（MOD13Q1）七：时间序列S-G滤波之Python 的代码。得到耕地NDVI时序图。（该NDVI时间范围内，如果种植的是某种作物。那么得到的耕地时序图，也是作物时序图。由于该数据选择的是3、4、5月MODIS数据。在该时间范围内耕地种植的是小麦。所以本耕地时序图，又叫小麦时序图。（不考虑误差））

 

 至此，NDVI基本系类处理结束。该系列后续会对之前的博文进行不断更进，以及适当补充。

对本系列（包括所有博文）有建议的，请您留言，与君共勉，非常感谢。

由于最近忙于其它，更新较慢，不足之处，还请见谅。

posted @ 2020-05-10 22:36  青灯常伴古佛  阅读(4145)  评论(0)  编辑  收藏  举报