关于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批处理数据的小小总结,欢迎提出问题及讨论啊。

  • 3
    点赞
  • 17
    收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
©️2022 CSDN 皮肤主题:大白 设计师:CSDN官方博客 返回首页
评论

打赏作者

weixin_42055454

你的鼓励将是我创作的最大动力

¥2 ¥4 ¥6 ¥10 ¥20
输入1-500的整数
余额支付 (余额:-- )
扫码支付
扫码支付:¥2
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、C币套餐、付费专栏及课程。

余额充值