Adding vertex to line at point location using ArcPy?

Question

I have a line layer and point layer on my ArcGIS project. The points are on the lines. I want to create vertex on lines at this point coordinates.

import arcpy

point_rows = arcpy.da.SearchCursor("D:\Test.gdb\MyPoints", ['OID@', 'SHAPE'])

rows = arcpy.UpdateCursor("D:\Test.gdb\MyLines")

for row in rows:
    geom = row.SHAPE
    for part in geom.getPart():
        for pnt in point_rows:
                row.X = pnt.X
                row.Y = pnt.Y
                rows.updateRow(row)

I could not add vertex on points, if point is on the line. How can I do this using arcpy?

Answer 1

Below is some python code that will take an un-ordered set of polylines with zero or more intersecting points and insert a vertex into the polyline at the point position.

Sample data

This code is designed to be dropped into the Python console in ArcMap, simply change lineLayer and pointLayer.

import arcpy

# Layers in MXD
lineLayer = "fcLine"
pointLayer = "fcPoint"

with arcpy.da.UpdateCursor(lineLayer,['SHAPE@']) as lineCursor:
    for row in lineCursor:
        geom = row[0]

        # Get selected points along line, this could be zero, 1 or more
        arcpy.SelectLayerByLocation_management(pointLayer,"INTERSECT",geom,"#","NEW_SELECTION")
        res = arcpy.GetCount_management(pointLayer)
        n = int(res.getOutput(0))
        if n > 0:
            # We have at least 1 point intersecting line
            # Store selected point(s) in list for further processing
            pointList = []
            with arcpy.da.SearchCursor(pointLayer,['SHAPE@']) as pointCursor:
                for pointRow in pointCursor:
                    pointList.append(pointRow[0])
            print(str(len(pointList)) + " intersect current line.")

            # Convert polyline to points
            # Assumes line is simple geometry, i,e a single part
            lineArr = geom.getPart(0) # Returns an array
            for p in lineArr:
                pointList.append(arcpy.PointGeometry(p))

            # PointList is a list of PointGeometries composed of points on the line and points from the line itself
            # Now need to work out position of points along the line
            posDict = {}
            for p in pointList:
                pos = geom.measureOnLine(p,True)
                posDict[pos] = p

            # Reconstruct polyline by using sorted keys in posDict
            newArray = arcpy.Array()
            for key in sorted(posDict.iterkeys()):
                newArray.add(posDict[key].centroid) # Array takes points not PointGeometry, hence using centroid property
            newLine = arcpy.Polyline(newArray)

            # Insert updated geometry line back into dataset
            row[0] = newLine
            lineCursor.updateRow(row)

Answer 2

I think I have a solution! I ve checked the following code and it works.

1) This example creates a Point Feature of the Vertices of your line with the help of Numpy Arrays (Numpy Needs to be installed)

2) Then it merges These Vertices as Points together with your 'points' on the line -> Theres one drawback: On the Output Pointfile of this merge you Need to create a new field and sort the Points in the right order (1-9 in this example) Otherwise the final line would look like a 'grinding by numbers' (but Ive found also examples to automate this too..

3) Finally a new line is created including the original line and the Points on the line

Feel free to ask me, in case you've got questions or some parts of the code dont work on your system

import arcpy, numpy, os

sr = arcpy.SpatialReference(25832)


### The Input: Your Points on the line (PointsLTW) // The line feature itself (InFcLine)

PointsLTW = r"C:\lineVerts\PointsLTW.shp"
InFcLine  = r"C:\\lineVerts\LineLTW.shp" # input feature class

OutFcLine = r"C:\lineVerts\LineLTW_out2.shp" # output feature class


print "Inputs accepted.."

# Creates a Point Array of your Line

arrayLine = arcpy.da.FeatureClassToNumPyArray(InFcLine,["SHAPE@XY"], "", sr, explode_to_points=True)

print "Numpy Array created.."

# ..in case no Array was created

if arrayLine.size == 0:
    arcpy.AddError(InFcLine + " contains no features!")

# From Array to Point Feature of your Input Line
else:
    arcpy.da.NumPyArrayToFeatureClass(arrayLine, OutFcLine, ['SHAPE@XY'], sr)
    print "New Multipoint Feature (LineToVertPoints) from Numpy Array.."

# Now you just merge your original PointsLTW with your original Line, now written as a Point Feature out of the lines Vertices

arcpy.Merge_management(inputs="PointsLTW;LineLTW_out2",output="C:/lineVerts/FinalPointsxx.shp",field_mappings="""Id "Id" true true false 6 Long 0 6 ,First,#,PointsLTW,Id,-1,-1,LineLTW_out2,Id,-1,-1""")

# Finally you create a line out of the merged points, the only drawback is that you need the Points to be sorted after the line Direction (From East to West or VV), I did this with a new field sorting the points by hand

arcpy.PointsToLine_management(r"C:/lineVerts/FinalPointsxx.shp" , r"D:\FinalLinexx.shp", "","sortedLTW")

Answer 3

As the easiest way (not too much coding, but questionable performance) you might try to do this (ArcGIS Desktop Advanced license required):

1. Use Split Line at Point tool to split the lines.
2. Use Unsplit Line tool to merge the lines back.

As described in the description of the tool:

If Search Radius (search_radius in Python) is unspecified, the nearest point will be used to split the line feature. If a Search Radius is specified, all points within the search radius are used to split the line.

So by default Split Line at Point only works for ONE point per line. Other points on the line are ignored. But if you set huge Search Radius all points will split the line.

________Before_____________________ Splitted_____________________ UNsplitted

_________ _________

This code should do the job (not tested):

import arcpy

points = r"D:\Test.gdb\MyPoints"
lines = r"D:\Test.gdb\MyLines"
splitted = r"D:\Test.gdb\MySplittedLines"
unsplitted = r"D:\Test.gdb\MyUnSplittedLines"
searchRadius = "9999999999999999999 Meters"

arcpy.SplitLineAtPoint_management(lines, points, splitted, searchRadius)
arcpy.UnsplitLine_management(splitted, unsplitted)

Answer 4

I think you had your geometry token for your points wrong. It should have come as "Shape@".

The code can be written as :

    Import arcpy 
    # Setting path 
    arcpy.env.workspace = "path" 
    # This line of code returns x, and y as tuple. 
    with arcpy.da.searchCursor("Feature/Table",['OID@','Shape@']) as /
    point_cursor
    for row in point_cursor: 
       pnt.X,pnt.Y = row[1]
       print('{},{}'.format(pnt.X,pnt.Y)) 
       ## The first code reads the x,y of the point layer 
       ### The next code is to write the point x,y on the line 
       with arcpy.InsertCursor('Feature/Table',"Shape@") as line_cursor 
         array = arcpy.array([pnt.X,pnt.Y]) 
         Line_ = arcpy.polyline(array) 
         line_cursor.insertRow([Line_])