hydrology op• Data kinds: depth → image2d
• Call: import fullseye as fs; fs.ledger.dem_stream_network(dem, cell_size, threshold_cells=100.0, fill=True, nodata='error') (to call the implementation directly, import demops; demops.dem_stream_network(dem, cell_size, threshold_cells=100.0, fill=True, nodata='error'); from the registry, opsdem.get("dem_stream_network"))
> This operator's description has not been translated yet. The original text follows as it is.
集水量が閾値を超えたセルを河道とみなす二値マスク。
手順: `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 family guide
• Sample-data catalog (download URLs / licences) — 2-D uses skimage.data (BSD/public domain) plus synthetic images; 3-D lists download URLs for real data sources (Stanford, PDS, …).
• Operator provenance and references — the sources of the research/methods this op family came from.
• The canonical algorithm (author, year) and its uses are named in the family usage guide above.
• dem_terrain_analysis_tour — py -3.11 examples/dem_terrain_analysis_tour.py
image2d as input)—
hydrology)dem_fill_sinks · dem_flow_direction · dem_flow_accumulation
*Provenance: demops.py — DEM operator registry. This per-op note is generated by tools/opdocs.py md (do not hand-edit).*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.