Assegnare ID punto ai rispettivi attributi di inizio e fine della polilinea?

Ho due shapefile, Point e Polyline. Sto cercando di popolare il mio livello di polilinea con l'ID dati punto (systemID) dove le linee iniziano a dove finiscono; per esempio. creare due colonne (da e a). Alle linee è già stata assegnata la direzione del flusso. Sto usando gli strumenti ArcHydro ma anche Arcpy va bene.

In ArcHydro ho provato a utilizzare Strumenti di attributo> Genera da / a nodi per le linee ma crea valori per A e Da che non provengono dai miei dati punto. Un'altra mia alternativa era usare Feature to Vertices e unire il risultato alla mia tabella dei dati dei punti per ottenere il systemID. Ma poi dovrei aggiornare anche il livello della polilinea. Sono sicuro che esiste un modo più semplice per farlo.

Nota: non tutte le linee hanno punti sui vertici, i valori vuoti vanno bene per loro.

A volte è meglio non usare una soluzione pronta all'uso. Questo è il motivo per cui suggerisco

1. Compilare i campi X e Y nel livello del nodo, convertirli in numeri interi, ad esempio cm. Crea un campo stringa e popolalo concatenando rappresentazioni di stringa di numeri interi.
2. Fai qualcosa di simile nella tabella dei collegamenti per il primo punto della forma.
3. Unisci la tabella dei nodi ai collegamenti utilizzando gli ultimi campi creati e trasferisci nodeID nel campo FROM.

Vai a 2, ma usa l'ultimo punto per assegnare TO nodeID

Ho fatto la stessa cosa qualche mese fa. Ho usato arcpy ma la sceneggiatura è troppo lunga per essere pubblicata qui, quindi ti darò una descrizione di ciò che ho fatto.

1. Ho usato un join spaziale per rilevare quali punti / nodi erano all'interno di una particolare funzione di linea.

2. Poiché l'unione spaziale non tiene conto della direzione del flusso, ho usato arcpy per accertare quale fosse il punto iniziale e quale il punto finale. Sono stato in grado di farlo utilizzando la funzione di descrizione arcpy per estrarre le informazioni sulle coordinate per i vertici di inizio / fine della funzione linea e confrontandole con i valori delle coordinate dei punti collegati.

3. Alla fine, una volta capito quali dei punti da / a erano quali ho usato la funzione setValue per popolare i campi da / a nel set di dati della polilinea originale.

Ovviamente, c'è un po 'più di questo, ma ho delineato i punti principali.

Sono stato ispirato da @FelixIP, ma volevo scrivere una soluzione senza join o la creazione di file extra, poiché la mia rete è abbastanza grande con 400K + pipe e 500K + nodi.

La costruzione della rete geometrica forza la X, Y dei nodi e le estremità del tubo coincidono. Puoi accedere a queste posizioni con i token forma in cursori arcpy e abbinarli. I token di forma per le linee restituiscono una matrice di vertici nell'ordine in cui sono stati disegnati. Nella mia rete, l'ordine di disegno dei tubi è fortemente QA perché lo usiamo per impostare le direzioni del flusso. Quindi, il primo vertice è l'inizio della pipe e l'ultimo vertice è la fine della pipe.

Riferimento: ASSETID = ID del tubo, UNITID = ID nodo all'inizio del tubo, UNITID2 = ID nodo alla fine del tubo.

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)