ARCGIS arcpy批量裁剪命令

最新推荐文章于 2024-09-14 12:48:47 发布

hydromonkey

最新推荐文章于 2024-09-14 12:48:47 发布

阅读量7.5k

点赞数 4

本文链接：https://blog.csdn.net/hengcall/article/details/82560994

版权

全部自动化批量处理，当然右键batch也可以处理，但是一旦有上百个文件，那就运动手指了。

代码如下，注释都有，根据需要改动参数即可。

import arcpy

# Set the current workspace
arcpy.env.workspace = ("F:/ZPATR-TIME-JOB/ndvi-china-/")  # "D:\\Test2\\Grids.gdb"

# Polygon Shapefiles Workspace
clipFeatures = ('F:/ZPATR-TIME-JOB/ndvi-china-/shp/xinjiang.shp/')  # "D:\\Test2\\Shps\\"

# Output Directory
outDirectory = ('F:/ZPATR-TIME-JOB/ndvi-china-/xinjiang/')  # "C:\\Users\\sgraf\\Documents\\ArcGISData\\Test2\\Outs.gdb\\"

# Get and print a list of GRIDs from the workspace
rasters = arcpy.ListRasters("*", "TIF")

for raster in rasters:
    arcpy.AddMessage("Clipping " + raster + " with " + raster + ".shp")

    desc = arcpy.Describe(clipFeatures + "\\" + raster + ".shp")
    extent = desc.extent

    arcpy.Clip_management(raster, str(extent), outDirectory + "\\" + raster, clipFeatures + "\\" + raster + ".shp", "#",
        "ClippingGeometry", "NO_MAINTAIN_EXTENT")