三种利用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['False_Easting',0.0],PARAMETER['False_Northing',0.0],PARAMETER['central_meridian',105.0],PARAMETER['Standard_Parallel_1',25.0],PARAMETER['Standard_Parallel_2',47.0],PARAMETER['latitude_of_

你可能感兴趣的:(《Python》,ArcGIS,空间数据,批处理,Python,投影)