Python批量裁切影像

作用：批量裁切指定目录下的影像数据

# --------------------------------------------------------------------------- # Created on: 周六 三月 06 2010 11:18:20 上午 # (generated by ArcGIS/ModelBuilder) # Usage: Batch extract rasters # --------------------------------------------------------------------------- # Import system modules import sys, string, os, arcgisscripting # Create the Geoprocessor object gp = arcgisscripting.create(9.3) # Check out any necessary licenses gp.CheckOutExtension("spatial") # Load required toolboxes... gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx") # Script arguments... Input_Raster_Folder = gp.GetParameterAsText(0) Boundary_Mask_Feature_Class = gp.GetParameterAsText(1) Output_Raster_Folder = gp.GetParameterAsText(2) # Process: extract rasters... gp.Workspace = Input_Raster_Folder tiffs =gp.ListRasters("*","TIF") for tiff in tiffs: try: Output_raster = Output_Raster_Folder + "/" + tiff gp.ExtractByMask_sa(tiff, Boundary_Mask_Feature_Class, Output_raster) gp.addwarning("Extraction of raster dataset " + tiff + " succeeded!") except: gp.addwarning(gp.getmessage(2))