Menetapkan ID titik untuk masing-masing atribut awal dan akhir polyline?

Saya memiliki dua shapefile, Point dan Polyline. Saya mencoba untuk mengisi lapisan polyline saya dengan ID data titik (systemID) di mana garis mulai ke tempat mereka berakhir; misalnya. buat dua kolom (dari dan ke). Garis sudah memiliki arah aliran yang ditugaskan padanya. Saya menggunakan alat ArcHydro tetapi Arcpy juga baik-baik saja.

Di ArcHydro saya telah mencoba menggunakan Alat Atribut> Hasilkan Dari / Ke Node untuk baris tetapi ini menciptakan nilai untuk dan dari yang bukan dari data titik saya. Alternatif saya yang lain adalah menggunakan Feature to Vertices dan menggabungkan hasilnya ke tabel data poin saya untuk mendapatkan systemID. Tetapi kemudian saya harus memperbarui lapisan polyline juga. Saya yakin ada cara yang lebih mudah untuk melakukannya.

Catatan: tidak semua garis memiliki titik di simpul, nilai kosong baik untuk mereka.

Terkadang lebih baik tidak menggunakan solusi out of the box. Inilah sebabnya saya sarankan

1. Populasikan bidang X dan Y di lapisan simpul, konversikan ke bilangan bulat, katakan cm. Buat bidang string dan isi dengan merangkai representasi string bilangan bulat.
2. Lakukan hal serupa di tabel tautan untuk titik pertama dalam bentuk.
3. Gabung tabel node ke tautan menggunakan bidang yang dibuat terbaru dan transfer nodeID ke bidang FROM.

Goto 2, tetapi gunakan titik terakhir untuk menetapkan TO nodeID

Saya melakukan hal yang sama beberapa bulan yang lalu. Saya menggunakan arcpy tetapi skripnya terlalu panjang untuk dikirim di sini jadi saya akan memberi Anda garis besar tentang apa yang saya lakukan.

1. Saya menggunakan gabungan spasial untuk mendeteksi titik / node mana dalam fitur garis tertentu.

2. Karena sambungan spasial tidak memperhitungkan arah aliran, saya menggunakan arcpy untuk memastikan titik awal dan titik akhir mana. Saya bisa melakukan ini dengan menggunakan fungsi menjelaskan busur untuk mengekstrak informasi koordinat untuk simpul mulai / akhir dari fitur garis dan membandingkannya dengan nilai koordinat titik-titik yang terhubung.

3. Akhirnya, setelah saya menentukan titik to / from mana yang saya gunakan fungsi setValue untuk mengisi bidang ke / dari dalam dataset polyline asli.

Jelas, ada sedikit lebih dari ini tetapi saya telah menguraikan poin utama.

Saya terinspirasi oleh @FelixIP, tetapi saya ingin menulis solusi tanpa bergabung atau membuat file tambahan, karena jaringan saya cukup besar dengan 400K + pipa dan 500K + node.

Membangun jaringan geometrik memaksa X, Y dari node dan ujung pipa menjadi bertepatan. Anda dapat mengakses lokasi ini dengan token bentuk dalam kursor busur dan mencocokkannya. Token bentuk untuk garis mengembalikan array simpul sesuai urutan penggambarannya. Di jaringan saya, urutan undian pipa sangat QA'd karena kami menggunakan ini untuk mengatur arah aliran. Jadi, simpul pertama adalah awal pipa, dan simpul terakhir adalah ujung pipa.

Referensi: ASSETID = id pipa, UNITID = simpul id di awal pipa, UNITID2 = simpul id di ujung pipa.

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)