第一推是谁做的网站,松江网站开发培训学校,竞价网官网,地铁网站建设特点任务分析#xff1a;在采用遥感方法研究植被物候变化时#xff0c;需要下载长时间序列#xff08;本文用到的是30年#xff09;GIMMS 3G NDVI数据。原始数据是全球的#xff0c;预处理过程包括几何校正、裁剪、投影变换、Hants时间序列谐波分析等等。在这过程中#xff0… 任务分析在采用遥感方法研究植被物候变化时需要下载长时间序列本文用到的是30年GIMMS 3G NDVI数据。原始数据是全球的预处理过程包括几何校正、裁剪、投影变换、Hants时间序列谐波分析等等。在这过程中由于是长时间序列数据需要重复操作当然ArcGIS软件提供了批量Batch操作功能但是需要我们手动添加比较费时费事此时我们首先想到的是用Python进行。 Extract by Mask Batch 既然Python很了不起那么我们就得利用它为我们服务。
一、完整程序代码
import arcpy
arcpy.CheckOutExtension(spatial)
arcpy.gp.overwriteOutput1
arcpy.env.workspace F:\\Modis_16\\1Moasic
rasters arcpy.ListRasters(*, tif)
mask F:\\Vegetation Change\\Data\\Bound\\bound_Buffer_Polygon.shp
for raster in rasters:print(raster)out F:\\Vegetation Change\\Data\\GIMMS Data\\new\\ma_raster[0:8]arcpy.gp.ExtractByMask_sa(raster, mask, out)print(ma_raster[0:8] has done)
print(All done)
二、注意事项 1.arcpy.gp.overwriteOutput1即覆盖之前的文件 2.输入的是.tif文件输出的是Grid文件 3.raster[0:8]表示从第0个开始取8个字符串 4.bound_Buffer_Plygon.shp后缀名不可或缺。
三、运行情况 刘一哥GIS专注测绘地理信息教育、探索地理奥秘、分享GIS价值
