python批处理栅格转点_三种利用Python批量处理地理数据的方法——以栅格数据投影转换为例...

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

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_origin',0.0],UNIT['Meter',1.0]]", "NEAREST", "9660.84811813477", "", "", "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0]