¿Asignar ID de punto a los respectivos atributos de inicio y fin de la polilínea?

Tengo dos shapefiles, Point y Polyline. Estoy tratando de llenar mi capa de polilínea con el ID de datos de punto (systemID) donde las líneas comienzan donde terminan; p.ej. crear dos columnas (desde y hasta). Las líneas ya tienen asignada la dirección del flujo. Estoy usando las herramientas de ArcHydro pero Arcpy también está bien.

En ArcHydro, he intentado usar Herramientas de atributo> Generar nodos de / a para líneas, pero crea valores para To y From que no son de mis datos de puntos. Mi otra alternativa era usar Feature to Vertices y unir el resultado a mi tabla de datos de puntos para obtener el ID del sistema. Pero entonces también tendría que actualizar la capa de polilínea. Estoy seguro de que hay una manera más fácil de hacerlo.

Nota: no todas las líneas tienen puntos en los vértices, los valores en blanco están bien para ellos.

A veces es mejor no usar una solución lista para usar. Por eso sugiero

1. Rellene los campos X e Y en la capa de nodo, conviértalos a enteros, digamos cm. Cree un campo de cadena y complételo concatenando representaciones de cadenas de enteros.
2. Haga algo similar en la tabla de enlaces para el primer punto de la forma.
3. Unir la tabla de nodos a enlaces utilizando los campos creados más recientemente y transferir el ID de nodo al campo DESDE.

Ir a 2, pero usar el último punto para asignar TO nodeID

Hice lo mismo hace unos meses. Utilicé arcpy pero el guión es demasiado largo para publicarlo aquí, así que te daré un resumen de lo que hice.

1. Utilicé una unión espacial para detectar qué puntos / nodos estaban dentro de una entidad de línea particular.

2. Debido a que la unión espacial no tiene en cuenta la dirección del flujo, utilicé arcpy para determinar cuál era el punto inicial y cuál era el punto final. Pude hacer esto usando la función de descripción arcpy para extraer información de coordenadas para los vértices de inicio / fin de la entidad de línea y los comparé con los valores de coordenadas de los puntos conectados.

3. Finalmente, una vez que resolví cuáles de los puntos para / desde los cuales usé la función setValue para llenar los campos para / desde en el conjunto de datos de polilínea original.

Obviamente, hay un poco más que esto, pero he esbozado los puntos principales.

Me inspiré en @FelixIP, pero quería escribir una solución sin combinaciones o la creación de archivos adicionales, ya que mi red es bastante grande con 400K + tuberías y 500K + nodos.

La construcción de la red geométrica obliga a que X, Y de los nodos y los extremos de la tubería coincidan. Puede acceder a estas ubicaciones con los tokens de forma en cursores de arco y combinarlos. Los tokens de forma para líneas devuelven una matriz de vértices en el orden en que fueron dibujados. En mi red, el orden de extracción de las tuberías está fuertemente controlado por calidad porque lo usamos para establecer las direcciones de flujo. Entonces, el primer vértice es el comienzo de la tubería, y el último vértice es el final de la tubería.

Referencia: ASSETID = identificación de la tubería, UNITID = identificación del nodo al inicio de la tubería, UNITID2 = identificación del nodo al final de la tubería.

nodes = "mergeNodes"
pipes = "SEWER_1"

nodeDict = {}
pipeDict = {}

#populate node dictionary with X,Y as the key and node ID as the value
for node in arcpy.da.SearchCursor(nodes, ["UNITID", "SHAPE@XY"]):
    nodeDict[(node[1][0], node[1][1])] = node[0]

#populate pipe dictionary with pipe ID as the key and list of X,Y as values 
#vertices populated in the order that the line was draw
#so that [0] is the first vertex and [-1] is the final vertex
for pipe in arcpy.da.SearchCursor(pipes, ["ASSETID", "SHAPE@"]):
    for arrayOb in pipe[1]:
        for point in arrayOb:
            if pipe[0] in pipeDict:
                pipeDict[pipe[0]].append((point.X, point.Y))
            else: 
                pipeDict[pipe[0]] = [(point.X, point.Y)]

#populate UNITID with the first vertex of the line
#populate UNITID2 with the final vertex of the line
with arcpy.da.UpdateCursor(pipes, ["ASSETID", "UNITID", "UNITID2"]) as cur:
    for pipe in cur:
        if pipeDict[pipe[0]][0] in nodeDict:
            pipe[1] = nodeDict[pipeDict[pipe[0]][0]]
        if pipeDict[pipe[0]][-1] in nodeDict:
            pipe[2] = nodeDict[pipeDict[pipe[0]][-1]]
        cur.updateRow(pipe)