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 패밀리 가이드
• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL 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 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.