该博客介绍了如何利用Python的arcpy模块和Spatial Analyst扩展功能，批量将MOD13Q1系列的NDVI、RVI、MSAVI、Red和Nir等不同波段的栅格值提取到对应的点特征上。代码遍历2011年至2020年的数据，将每个波段的栅格值作为新字段添加到点数据集，实现了自动化处理。
import arcpy
import numpy as np
from arcpy.sa import *
arcpy.CheckOutExtension("spatial")
arcpy.env.workspace = "I:\\RSinstitue\\MOD13Q1\\AGB_VI" #点图层
for y in range(2011,2020):
inPointFeatures = "I:\\RSinstitue\\MOD13Q1\\AGB_VI\\AGB_VI_shp\\mod13qAGB"+str(y)+".shp"
inRasterList = [["K:\\Wu_Data\\MOD13Q_Band\\MOD13Q1_NDVI\\MOD13Q"+str(y)+"177_NDVI.tif", "ND626_711"], #输入栅格，字段名称不超10位
["K:\\Wu_Data\\MOD13Q_Band\\MOD13Q1_NDVI\\MOD13Q"+str(y)+"193_NDVI.tif", "ND712_727"],
["K:\\Wu_Data\\MOD13Q_Band\\MOD13Q1_NDVI\\MOD13Q"+str(y)+"209_NDVI.tif", "ND728_812"],
["K:\\Wu_Data\\MOD13Q_Band\\MOD13Q1_NDVI\\MOD13Q"+str(y)+"225_NDVI.tif", "ND813_828"],
["K:\\Wu_Data\\MOD13Q_Band\\
最低0.47元/天 解锁文章
扫一扫
1155
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
