Python批量赋Spatial Reference

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

作用：给指定目录下的栅格数据批量赋Spatial Reference

 

[python]  view plain copy

1. # ---------------------------------------------------------------------------  
2. # define sr.py  
3. # Created on: 周六 三月 06 2010 11:18:20 上午  
4. #   (generated by ArcGIS/ModelBuilder)  
5. # Usage: define sr <Input_Dataset_or_Feature_Class> <Coordinate_System>   
6. # ---------------------------------------------------------------------------  
7. # Import system modules  
8. import sys, string, os, arcgisscripting  
9. # Create the Geoprocessor object  
10. gp = arcgisscripting.create(9.3)  
11. # Load required toolboxes...  
12. gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")  
13. # Script arguments...  
14. Input_Raster_Folder = gp.GetParameterAsText(0)  
15. Coordinate_System = gp.GetParameterAsText(1)  
16. # Process: Define Projection...  
17. gp.Workspace = Input_Raster_Folder  
18. tiffs =gp.ListRasters("*","TIF")  
19. for tiff in tiffs:  
20.         try:  
21.                 #gp.defineprojection(tiff, Coordinate_System)  
22.                 gp.DefineProjection_management(tiff, Coordinate_System)  
23.                 gp.addwarning("Define raster dataset " + tiff + " spatial reference succeed!")  
24.         except:  
25.                 gp.addwarning(gp.getmessage(2))