程序代码：
[python]view plaincopyprint?
1. importarcpy
2. arcpy.CheckOutExtension("spatial")
3. arcpy.gp.overwriteOutput=1
4. arcpy.env.workspace="F:\Modis_16\1Moasic"
5. rasters=arcpy.ListRasters("*","tif")
6. mask="F:\VegetationChange\Data\Bound\bound_Buffer_Polygon.shp"
7. forrasterinrasters:
8. print(raster)
9. out="F:\VegetationChange\Data\GIMMSData\
new\"+"ma_"+raster[0:8]
10. arcpy.gp.ExtractByMask_sa(raster,mask,out)
11. print("ma_"+raster[0:8]+"hasdone")
12. print("Alldone")
二、注意事项：
1.arcpy.gp.overwriteOutput=1 即覆盖之前的文件；
2.输入的是.tif 文件，输出的是 Grid 文件；
3.raster[0:8]表示从第 0 个开始取 8 个字符串；
4.bound_Bu&er_Plygon.shp 后缀名不可或缺。
三、运行情况：
评论0