hydrology opdepth → image2dimport fullseye as fs; fs.ledger.dem_flow_accumulation(dem, cell_size, fill=True, epsilon=1e-06, nodata='error') (実装を直接呼ぶなら import demops; demops.dem_flow_accumulation(dem, cell_size, fill=True, epsilon=1e-06, nodata='error')、台帳から引くなら opsdem.get("dem_flow_accumulation"))集水セル数。各セルへ流れ込む上流セルの個数(自セルを 1 と数える)。
fill=True なら先に窪地を埋める(埋めないと窪地で流れが止まり、
下流の集水量が静かに小さく出る)。欠測セルの集水量は nan。
グラフの構築はベクトル化してある。ただし効果は小さい: 1024x1024 で 2.71 秒 -> 2.39 秒。支配しているのは構築ではなく、そのあとのトポロジカル 走査(セル数ぶんの Python ループ)のほうだった。速くしたいならそこを 書き換える必要がある —— 直したつもりで直っていない、を残さないため明記する。
py -3.11 examples/dem_terrain_analysis_tour.pypy -3.11 examples/poc_dem_terrain.pyimage2d を入力に取れる)—
hydrology)dem_fill_sinks · dem_flow_direction · dem_stream_network
Provenance: demops.py — DEM operator registry. この per-op ノートは tools/opdocs.py md が自動生成(手編集しない)。
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.