关于arcpy中批处理栅格及矢量文件的总结

总所周知，GIS是处理栅格数据和矢量数据的极其有效的工具，但是当我们有大量数据时，用GIS自带的TOOL一个一个选择添加然后处理就会非常繁琐，这时候Arc GIS中集成的工具arcpy的优势就显现出来了，可以使用Geoprocessing工具条中的Python直接输入语句运行
当然更常用的方式是在Arctoolbox中添加Toolbox，然后导入.py脚本，直接做成类似工具条的模式以后可以自有调用，具体方式很多文章都有讲到，就不一一说明了。
另外arc gis desktop的网站上详细讲解了arcpy的相关语句和命令，可在上面查找相关命令，另外也可直接在Arctoolbox上直接右击命令条，选择item description可以看到命令的源Python代码。下面就我的工作谈谈具体批量处理数据，理论上只要是重复的工作都是可以用代码实现的。
批量处理首先涉及到读取文件列表和循环的问题，arcpy中可以直接用ListRasters或者ListFeatrues工具把文件夹中所有的矢量和栅格文件罗列出来，但是只能罗列一层文件夹中的数据，如果要罗列多层文件夹中的数据可以用语句：

for root, dirs, folders in os.walk(path):
    for folder in folders:
        if (folder.endswith('.shp')):

首先我要把二值化的栅格图像批量转换为shp数据，我将栅格图像放在同一个文件夹中，代码如下：

import arcpy
from arcpy import env

env.workspace = arcpy.GetParameterAsText(0)
field = "VALUE"
outputworkspace = arcpy.GetParameterAsText(1)
files = arcpy.ListRasters("*","tif")
for file in files:
    pos=file.index('.')
    outfile=file[0:pos]
    year=file[6:10]
    outPolygons =outputworkspace + "/" + year + "/" + outfile
    arcpy.RasterToPolygon_conversion(file, outPolygons, "SIMPLIFY", field)
print "finish"

随后我要通过读取shp文件的属性表，对shp文件的字段进行删除、添加、计算等操作，主要语句如下：

for root, dirs, folders in os.walk(path):
    for folder in folders:
        if (folder.endswith('.shp')):
            out_feature= outpath + "/" + folder #输出数据
            upcur = arcpy.da.UpdateCursor(folder, “id”) #遍历数据字段“id”的所有行
            for row in upcur:
                if (row[0] == 0 or row[0] == 32767): #找出“id”中数字为0和32767的行
                    upcur.deleteRow() #删除行，即删除要素
            arcpy.Clip_analysis(folder, outline_feature, out_feature, "") #剪切数据
            arcpy.Delete_management(outline_feature) #删除矢量数据，根据自己需要选择这个功能
        ##添加面积字段
        arcpy.AddField_management("folder",'Areakm2',"FLOAT","9","4")

随后可以检查一下数据剪切得怎样，随后可以进行计算字段和将计算的面积写入TXT文档的操作，也可以写入EXCEL，但是比较麻烦，可以按一定间隔写入TXT中，后直接在excel中打开，主要代码如下：

for root, dirs, folders in os.walk(path):
    for folder in folders:
        if (folder.endswith('.shp')):
            #计算多边形面积
            #arcpy.CalculateField_management(folder, "Areakm2", "!shape.geodesicArea@SQUAREKILOMETERS!", "PYTHON_9.3")
            #将属性表内容读取并输出至TXT文档中
            cursor2 = arcpy.da.SearchCursor(folder, ["name", "Areakm2"])
            for row in cursor2:
                f.write(str(row[0]) + " ")
                f.write(str(row[1]) + '\r\n')

以上为学习arcpy批处理数据的小小总结，欢迎提出问题及讨论啊。