arcpy_intersectmerge

主要想实现两条直线（其中一个线shp是固定的，另一个线shp是多个）之间交点的批量处理，并对结果进行整合，方便后续加载数据和增加经纬度坐标信息。

1 . 利用arcpy 中Intersect_analysis函数实现两直线相交功能，具体代码：
输入数据：
smallshape 为需要进行相交处理的多个线shp；
region 为某个固定的线shp.
目的：将smallshape 中每个线shp与region 进行intersect处理，得到交点shp.

import arcpy,os,datetimefrom arcpy.sa import *import pandas as pdimport refrom arcpy import envpath = r"存放需要处理的shp文件路径"arcpy.env.workspace = patharcpy.env.overwriteOutput = Trueprint("strart")smallshapes  = [f for f in arcpy.ListFeatureClasses() if f.split('.')[0].isdigit()]#List a:print(smallshapes)#infeature = ["big.shp", "small1.shp"]join_attributes = "ALL"#cluster_tolerance = "-1 Unknown"output_type = "POINT"region = r"固定的线shp存放路径" #一个固定的shpi = 0for smallshape in smallshapes:#outfeature = r"输出文件路径" + "{0}.shp".format(i)outfeature = r"输出文件路径名称\intersect_" + str(i) #该程序输出路径命名为intersect_1,intersect_2，如果有其它要求，可以修改代码。print(smallshape)arcpy.Intersect_analysis([smallshape, region], outfeature, join_attributes, "", output_type)i+=1print("end")

参考资料：
Using ArcPy with Intersect and loop：https://gis.stackexchange.com/questions/358430/using-arcpy-with-intersect-and-loop

ListFeatureClasses：https://desktop.arcgis.com/zh-cn/arcmap/latest/analyze/arcpy-functions/listfeatureclasses.htm

利用Intersect_analysis进行分区统计：https://zhuanlan.zhihu.com/p/61445743

相交 (分析)：https://pro.arcgis.com/zh-cn/pro-app/latest/tool-reference/analysis/intersect.htm

2 . 将上面的结果整合merge成一个shp
注意这个步骤，需要控制输入点shp的个数，如果数量超出设定的范围，就会报错：

但是没有找到更好的方法，只能暂时分批次整合（每次小于1000进行操作）。

具体代码：

import arcpy,osout_path =r"整合后文件路径" #整合后的文件路径arcpy.env.workspace = r"需要整合的文件路径" #需要整合的shp路径shplist = arcpy.ListFeatureClasses('*.shp')print(shplist)arcpy.Merge_management(shplist,os.path.join(out_path,'Merged_intersection.shp'))print("done")

参考资料：
Simple Merge of Shapefiles in a single folder：https://community.esri.com/t5/python-questions/simple-merge-of-shapefiles-in-a-single-folder/td-p/337113

arcpy批量合并shp:https://blog.csdn.net/imxiezy/article/details/88605946

arcgis merge：https://desktop.arcgis.com/zh-cn/arcmap/10.3/tools/data-management-toolbox/merge.htm

arcpy_intersectmerge

最新文章

热门文章