arcpy_intersect&merge

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

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

import arcpy,os,datetime

from arcpy.sa import *

import pandas as pd

import re

from arcpy import env

path = r"存放需要处理的shp文件路径"

arcpy.env.workspace = path

arcpy.env.overwriteOutput = True

print("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存放路径" #一个固定的shp

i = 0

for 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+=1

print("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,os

out_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