ARCGIS中Python实现批量裁剪
合集下载
相关主题
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
程序代码：
[python]view plaincopyprint?
1.import arcpy
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:\\Vegetation Change\\Data\\Bound\\bound_Buffer_Polygon.shp"
7.for raster in rasters:
8.print(raster)
9. out= "F:\\Vegetation Change\\Data\\GIMMS Data\\new\\"+"ma_"+raster[0:8]
10. arcpy.gp.ExtractByMask_sa(raster, mask, out)
11.print("ma_"+raster[0:8]+" has done")
12.print("All done")
二、注意事项：
1.arcpy.gp.overwriteOutput=1即覆盖之前的文件；
2.输入的是.tif文件，输出的是Grid文件；
3.raster[0:8]表示从第0个开始取8个字符串；
4.bound_Buffer_Plygon.shp后缀名不可或缺。
三、运行情况：