ポリラインの開始と終了の属性にそれぞれポイントIDを割り当てますか？

ポイントとポリラインの2つのシェープファイルがあります。ポリラインレイヤーに、ラインの始点から終点までのポイントデータID（systemID）を入力しようとしています。例えば。2つの列（fromおよびto）を作成します。ラインにはすでにフローの方向が割り当てられています。私はArcHydroツールを使用していますが、Arcpyも問題ありません。

ArcHydroで、属性ツール>ラインのFrom / Toノードの生成を使用してみましたが、ポイントデータからではないToおよびFromの値が作成されます。もう1つの方法は、Feature to Verticesを使用して、その結果をポイントデータテーブルに結合し、systemIDを取得することでした。ただし、ポリラインレイヤーも更新する必要があります。もっと簡単な方法があるはずです。

注：すべての線の頂点にポイントがあるわけではありません。空白の値で十分です。

時には、すぐに使えるソリューションを使用しない方がよい場合があります。これが私が提案する理由です

1. ノードレイヤーにXフィールドとYフィールドを入力し、それらを整数（cmなど）に変換します。文字列フィールドを作成し、整数の文字列表現を連結して入力します。
2. シェイプの最初のポイントについて、リンクテーブルで同様のことを行います。
3. 最後に作成されたフィールドを使用してノードテーブルをリンクに結合し、nodeIDをFROMフィールドに転送します。

移動2、ただし最後のポイントを使用してTOノードIDを割り当てる

私は数ヶ月前に同じことをしました。私はarcpyを使用しましたが、スクリプトが長すぎてここに投稿できないので、私が行ったことの概要を説明します。

1. 空間結合を使用して、特定のラインフィーチャ内にあるポイント/ノードを検出しました。

2. 空間結合では流れの方向が考慮されていないため、開始点と終了点を確認するためにarcpyを使用しました。これを行うには、arcpy describe関数を使用してラインフィーチャの開始/終了頂点の座標情報を抽出し、これらを接続されたポイントの座標値と比較しました。

3. 最後に、どちらのto / fromポイントがどれであるかを調べたら、setValue関数を使用して、元のポリラインデータセットのto / fromフィールドに値を入力しました。

もちろん、これ以外にも少しありますが、要点をまとめました。

@FelixIPに触発されましたが、ネットワークが400K以上のパイプと500K以上のノードで非常に大きいため、結合や追加ファイルの作成なしでソリューションを書きたかったのです。

ジオメトリックネットワークの構築により、節点とパイプの両端のX、Yが一致するように強制されます。arcpyカーソルの形状トークンを使用してこれらの場所にアクセスし、それらを一致させることができます。線の形状トークンは、頂点の配列を、描画された順序で返します。私のネットワークでは、フローの方向を設定するためにこれを使用するため、パイプの描画順序は大幅にQAされます。つまり、最初の頂点はパイプの始まりであり、最後の頂点はパイプの終わりです。

参照：ASSETID =パイプのID、UNITID =パイプの開始時のノードID、UNITID2 =パイプの終了時のノードID

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)