Python批量裁切影像

转载链接：http://blog.youkuaiyun.com/heyubingzju/article/details/5845139

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

[python]  view plain copy

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