dem_stream_network — DEM `hydrology` op

• 資料種類:depth → image2d

• 呼叫: import fullseye as fs; fs.ledger.dem_stream_network(dem, cell_size, threshold_cells=100.0, fill=True, nodata='error')(要直接呼叫實作,import demops; demops.dem_stream_network(dem, cell_size, threshold_cells=100.0, fill=True, nodata='error');從台帳取用則 opsdem.get("dem_stream_network"))

用法

> 該運算子的說明尚無譯文,以下照原文給出。

集水量が閾値を超えたセルを河道とみなす二値マスク。

手順: `dem_flow_accumulation(dem, cell_size, fill, nodata)` で各セルの集水

セル数(自セルを 1 と数える)を取り、`acc >= threshold_cells` を 1、それ以外を 0

にする。欠測セルは `nan` のまま。

• `threshold_cells: 河道とみなす集水セル数。正の有限値(ValueError`)。

面積で決めたいなら `面積 [m^2] / cell_size^2` に換算して渡す(この op は

セル数しか受けない)。

• `fill=True(既定): 先に priority-flood(dem_fill_sinks、epsilon=1e-6`)で

窪地を埋める。埋めないと窪地で流れが止まり、下流の集水量が小さく出て河道が

途切れる。

• `nodata: "error"(既定、nan があれば拒否)/ "outlet"`(欠測 = 流出口、

水域向き)/ `"barrier"`(欠測 = 壁)。定数で穴埋めする選択肢は無い。

• `dem / cell_size の契約は dem_slope` と同じ(3x3 以上、番兵値拒否、[m])。

• 返り値: `(H, W) float64 の 0 / 1(欠測は nan`)。bool ではないので

`blob_label 等に渡すなら == 1` で二値化する。

• 計算量: 集水量のトポロジカル走査がセル数ぶんの Python ループ。

D8 流向(`dem_flow_direction`)に基づくので、平坦な埋立地では流れが一箇所に

集まりやすい ―― 閾値を下げすぎると格子状の偽河道が出る。

詳細使用指南

• dem_terrain_analysis 族使用指南

參考(範例資料・文獻)

• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。

• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。

• 演算法的正典(作者・年份)與用途見上面的族使用指南。

可執行的範例(實際呼叫該運算子並已驗證的樣例)

• dem_terrain_analysis_tour — py -3.11 examples/dem_terrain_analysis_tour.py

型別可銜接的下一個運算子(可接受 `image2d` 作為輸入)

—

同類別(`hydrology`)

dem_fill_sinks · dem_flow_direction · dem_flow_accumulation

*Provenance: demops.py — DEM 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*

dem_stream_network — DEM hydrology op

用法