NDVI 数据处理,及估算植被覆盖度(一)

若是研究区在中国,则在地理空间数据云有中国区合成好的数据,直接下载用即可。数据格式为TIF。坐标为WGS84,不用再转换投影。下载网址为:http://www.gscloud.cn/sources/?cdataid=265&pdataid=10

若研究区在国外,需从nasa下载modis标准数据产品,进行投影转换等工作。后续有时间会补充如何转换投影等内容,可参见《modis AMSR 雪当量数据处理》这篇博客。


下载完成后,会得到一系列TIFF文件,由于一般研究区都是在某个流域或者某个省或行政区,需要裁剪TIF,拿到需要的范围。


以NDVI 月数据为例,本次下载数据为中国2005-2015年的逐月NDVI。需要的结果:

1:逐年逐月NDVI的时间序列变化。

2:逐年逐月NDVI的空间变化。

3:根据NDVI估算区域植被覆盖度,并给出空间图。

 

处理步骤

1:准备研究区的shp文件,在arcis裁剪出需要的范围。

在ARCGIS中选择cilp(裁剪选项,注意是Datamanagement中的裁剪,不是analysis中的)。若数据不多,则手动点击即可,但本例数据共132副图像,需要批量处理。

批量处理的脚本如下,打开arcgis的python控制器,出现如下界面:

 

代码在下方,每个都有注释,把目录改成自己需要的,粘贴到控制器,点确定,即可运行出结果。

 

import arcpy
import os
ras_path = 'F:/ZPATR-TIME-JOB/ndvi-china-/'      # TIF数据的路径
shp_path = 'F:/ZPATR-TIME-JOB/ndvi-china-/shp/xinjiang.shp'   # 用于裁剪的shp文件
outpath = 'F:/ZPATR-TIME-JOB/ndvi-china-/xinjiang/'          # 输出文件的路径

desc = arcpy.Describe(shp_path)
extent = desc.extent

files = [f for f in os.listdir(ras_path) if f.endswith('TIF')]   
for file in files:
    arcpy.Clip_management(ras_path+file,str(extent), outpath+file , shp_path, "#","ClippingGeometry", "NO_MAINTAIN_EXTENT")

 处理完成后得到裁剪后的研究区NDVI逐月的图像。

 

 

 

 

后续处理请见下篇。

展开阅读全文

没有更多推荐了,返回首页