How to delete vertices from a line when inside a selecting polygon

LindsayRaabe_FPCWA

Hi Brains Trust. I'm doing some cleanup work on a road network to make it suitable for Network Analysis, and have come across the following issue. One of the steps I have is usin the Integrate tool to align spatially different datasets. Unfortunately, a side effect of this is that lines passing across each other which normally don't intersect (as in no co-located vertices which Network Analyst uses to permit turns) now do. I've identified a lot of overpass/underpass areas which need to have these no-turn intersections cleaned, but can't find a way to bulk delete vertices when identified by an intersecting polygon. Any ideas out there on how to achieve this?

 

Below screenshot shows where a route turns off a bridge instead of following along the off ramp. 

 

LindsayRaabe_FPCWA

Hi. Do you have an advanced license? In a nut shell, I use one of these to:

1. Convert line vertices to points (Feature Vertices To Points (Data Management)—ArcGIS Pro | Documentation - ADVANCED license)
2. Erase the points that fall within the polygon (Erase (Analysis)—ArcGIS Pro | Documentation - ADVANCED license)
3. Reconstruct the lines from the remaining points (Points To Line (Data Management)—ArcGIS Pro | Documentation - BASIC license)

At the time I didn't figure out how to do it without an advanced license, but with some playing around with python, now I have! 

Starting data:

1. Convert lines to points using below python code:

import arcpy

# Set paths
input_fc = r"C:\YourProject\YourGDB.gdb\YourPolylineLayer"
output_fc = r"C:\YourProject\YourGDB.gdb\PolylineVertices"

# Use the same spatial reference as the input
spatial_ref = arcpy.Describe(input_fc).spatialReference

# Create output point feature class
arcpy.CreateFeatureclass_management(
    out_path=r"C:\YourProject\YourGDB.gdb",
    out_name="PolylineVertices",
    geometry_type="POINT",
    spatial_reference=spatial_ref
)

# Add a field to track original feature ID (optional)
arcpy.AddField_management(output_fc, "SourceID", "LONG")

# Insert points
with arcpy.da.SearchCursor(input_fc, ["OID@", "SHAPE@"]) as search_cursor, \
     arcpy.da.InsertCursor(output_fc, ["SHAPE@", "SourceID"]) as insert_cursor:

    for oid, polyline in search_cursor:
        for part in polyline:
            for point in part:
                if point:
                    insert_cursor.insertRow((point, oid))

 

2. Use Select by Location to select undesirable points

3. Use Delete Rows to delete selected points

4. Use Points to Line to reconstruct lines from remaining points

Turn that python code into a python tool and you can string that into a model with the rest of the tools to repeat as needed.