[Python知识库]python（arcpy）对矢量线进行批量计算长度和筛选（按属性筛选，查找最大最小值，选出最长线）

如果不想批量处理，直接把dataname替换成需要处理的数据，把遍历的循环删掉。
需要替换的参数较多，没有写成函数…初学，写的不太好，供参考。

整体的思路，就是挑选出属性符合的线，计算这些线的长度，然后查询最大最小值，当满足长度等于最大值时，被挑选出来。最终的结果呢就是：选出符合条件的线里的最长的那个

# -*- coding: utf-8 -*-
import os
import arcpy
import np
from arcpy import env
from arcpy.sa import *


#按属性选择800mm
arcpy.env.workspace =""#这里填入矢量数据所在的文件夹
path = ""#矢量数据位置
datanames = os.listdir(path)
for dataname in datanames:
    if dataname.endswith(".shp"):
        outpolyline = ""#输出位置及文件名.shp
        print(dataname+"_800mmTABLE SQL")
        arcpy.Select_analysis(dataname, outpolyline,'"CONTOUR" = 800')#这里的参数可以修改

#添加字段并计算长度
arcpy.env.workspace =""#这里填入矢量数据所在的文件夹
path =""#矢量数据位置
datanames = os.listdir(path)
for dataname in datanames:
    if dataname.endswith(".shp"):
        arcpy.AddField_management(dataname, "length", "SHORT", "", "", "#", "#", "#", "#")
        print(dataname[-8:-4] + "_ADD_FIELD")
        arcpy.CalculateField_management(dataname, "length", "!shape.geodesicLength@KILOMETERS!", "PYTHON_9.3")

#查找属性中的最大值最小值
field = "length"
max_value = []; k = 0
datanames = os.listdir(path)#注意path，如果路径有变需要重新定义
for dataname in datanames:
    if dataname.endswith(".shp"):
        datapath = path + "/" + dataname
        valuelist = []; i = 0
        for row in arcpy.da.SearchCursor(datapath, "length"):
            valuelist.append(row)
        max_value.append(max(valuelist))
        k += 1
max_value =np.array(max_value, dtype=int)
max_value = str([int(x) for item in max_value for x in item])
max_value = max_value[1:-1].split(',', k)
for i in range(20):
    max_value[i]=max_value[i].strip()
print(max_value,len(max_value))

#按属性选择最长线
arcpy.env.workspace =""#这里填入矢量数据所在的文件夹
path = ""#矢量数据位置
i=0
datanames = os.listdir(path)
for dataname in datanames:
    if dataname.endswith(".shp"):
        outpolyline = ""#输出位置及文件名.shp
        print(dataname[-8:-4]+"_MAXTABLE SQL")
        arcpy.Select_analysis(dataname, outpolyline,'"length"='+ max_value[i])
        os.remove(""+dataname)#这一步是把原数据删除，不需要该操作可以删掉
        i+=1


print ("成功")

以以上！OwO