ArcGIS Pro脚本工具（3）——直线转弧线

OD成本矩阵求解后的生成线为直线，如下左图，略丑。一通捣鼓终于实现了将直线转为弧线，效果见右图。

实现的方法主要是参考了这片文章，在此先感谢前辈大佬。

生成流向弧线工具 - GIS知乎-新一代GIS问答社区http://zhihu.geoscene.cn/article/2811在前文的基础上，修改了添加折点的方法（为此还解了一道几何题，几乎抓破脑壳），增加了对弧线方向和弧度的控制。值得一提的是，将原直线的属性信息也一并增加到了新生成的弧线中，方便后续的符号化。

工具脚本

    # -*- coding: utf-8 -*-
     
    #引入模块
    import os, arcpy
    import arcpy.cartography as CA
     
     
    #输入原始的line
    line = arcpy.GetParameterAsText(0)
     
     
    #计算弧线的顶点坐标的函数
    k=arcpy.GetParameter(2)
    d=arcpy.GetParameterAsText(3)
     
    def addPoint(x0, y0, x2, y2):
        if d =="顺向":
            x1,y1 = (x2+x0)/2-(y2-y0)/k,(y2+y0)/2+(x2-x0)/k
        elif d=="逆向":
            x1,y1 = (x2+x0)/2+(y2-y0)/k,(y2+y0)/2-(x2-x0)/k
        return (x1, y1)  
     
      
    # 创建一个临时线要素类addVertex
    addVertex = arcpy.CreateFeatureclass_management("in_memory", "addVertex", "POLYLINE", spatial_reference = arcpy.Describe(line).spatialReference)
    arcpy.AddField_management(addVertex, "ORIGID", "LONG")  
     
     
    #创建与line对应的折线并写入addVertex
    with arcpy.da.SearchCursor(line, ["SHAPE@", "OID@"]) as sc:  
        with arcpy.da.InsertCursor(addVertex, ["SHAPE@", "ORIGID"]) as ic:  
            for row in sc:  
                shp = row[0]
                p0 = shp.firstPoint  
                p2 = shp.lastPoint
                x1,y1 = addPoint(p0.X, p0.Y, p2.X, p2.Y)  
                p1 = arcpy.Point(x1,y1)  
                new_geometry = arcpy.Polyline(arcpy.Array([p0, p1, p2]))  
                new_row = [new_geometry, row[1]]  
                ic.insertRow(new_row)
     
     
    #调用SmoothLine工具，输入折线addVertex生成弧线curveline
    curveline = arcpy.GetParameterAsText(1)
    CA.SmoothLine(addVertex, curveline, "BEZIER_INTERPOLATION", "", "", "NO_CHECK")
     
     
    #将原始line的属性信息连接到curveline中
    try:
        arcpy.management.JoinField(curveline, "ORIGID", line, "OBJECTID", )
    except:
        arcpy.management.JoinField(curveline, "ORIGID", line, "FID", )
     
     
    # Delete inmemo fc
    arcpy.Delete_management("in_memory")

 工具参数

 工具使用

 工具下载

代码请自取，完整工具请私信联系
工具说明

有网友获取了代码制作了Pro和ArcMap工具箱，这里需要说明，ArcMap10.7中使用该代码制作的工具会出现错误，弧线曲度较大的情况下生成不了预想中的弧线，原因是ArcMap中的平滑线工具不能正常工作。10.8正常，其他版本没有试过。
————————————————
版权声明：本文为CSDN博主「学学GIS」的原创文章，遵循CC 4.0 BY-SA版权协议，转载请附上原文出处链接及本声明。
原文链接：https://blog.csdn.net/baidu_28157641/article/details/124437500