【ArcGIS遇上Python】三种利用Python批量处理地理数据的方法——以栅格数据投影转换为例

时至今日，笔者已经总结了三种用Python语言结合ArcGIS10.x提供的接口去批量处理地理空间数据的方法。即：

      1. 用IDLE、Python Tools for Visual Studio等去编写独立的py脚本文件；

      2. 用Python提供的tKinter模块去构建可视化窗体，并将其打包成exe可执行文件；

      3. 用Python脚本工具自定义工具箱。

     下面对三种方法分别做介绍。

一、用IDLE、Python Tools for Visual Studio或PyCharm等去编写py文件；

arcpy.env.workspace = "F:\\GIMMS 3g\\3Mask"
rasters = arcpy.ListRasters("*", "tif")
for raster in rasters:
    out= "F:\\GIMMS 3g\\4Project\\"+"pro_"+raster[4:12]+".tif"
    arcpy.ProjectRaster_management(raster,out, "PROJCS['MyAlbers',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Albers'],PARAMETER['Fals