surface op• 데이터 종류: depth → image2d
• 호출: import fullseye as fs; fs.ledger.dem_slope(dem, cell_size, method='horn', units='degrees')(구현을 직접 호출하려면 import demops; demops.dem_slope(dem, cell_size, method='horn', units='degrees'), 원장에서 가져오려면 opsdem.get("dem_slope"))
경사각. units는 "degrees" / "radians" / "percent".
평면이면 닫힌 식 atan(|grad|)과 일치한다(실측 1e-13 이하).
절차: 3x3 근방에서 dz/dx(동쪽 방향)와 dz/dy(북쪽 방향 = 행이 줄어드는 방향)를취하고, g = hypot(dz/dx, dz/dy). "degrees"는 atan(g)를 도로, "radians"는 라디안, "percent"는 100 g(45도 = 100 %). 가장자리는 끝 값을 복제해 채우므로(pad mode="edge"), 바깥 둘레 1셀은 안쪽보다완만하게 나온다.
• dem: (H, W)의 표고 [m], 3x3 이상, 실수, inf 불가. 결측은 nan으로전달한다 ―― -9999 같은 센티넬 값이 그대로 들어 있으면(<= -9000이 있고 nan이 없으면) 거부한다. nan은 3x3 범위로 전파된다.
• cell_size: 셀 한 변 길이 [m], 양의 유한값. bool / 문자열은 거부. 위경도 격자는대신 dem_geodetic_slope 쪽을 쓴다.
• method: "horn"(Horn 1981의 3x3 가중 차분, 기본값) / "central" (Zevenbergen–Thorne의 중앙 차분).
• 반환값: (H, W) float64. degrees는 [0, 90).
• 실패는 모두 ValueError(형상, 센티넬 값, cell_size, 선택지).
dem_aspect와 짝으로 쓴다. dem_hillshade는 이 2개로 음영을 만든다.
• dem_terrain_analysis 패밀리 가이드
• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL URL.
• 연산자의 내력·참고문헌 —— 이 연산자 족의 바탕이 된 연구/기법의 출처.
• 알고리즘의 정전(저자·연도)과 용도는 위의 패밀리 사용 가이드에 적혀 있습니다.
• dem_geodesy_tour — py -3.11 examples/dem_geodesy_tour.py
• poc_crop_phenotyping — py -3.11 examples/poc_crop_phenotyping.py
• poc_dem_terrain — py -3.11 examples/poc_dem_terrain.py
• poc_geodetic_height_frames — py -3.11 examples/poc_geodetic_height_frames.py
• poc_lidar_terrain_change — py -3.11 examples/poc_lidar_terrain_change.py
• poc_multibeam_bathymetry — py -3.11 examples/poc_multibeam_bathymetry.py
• poc_stockpile_volume — py -3.11 examples/poc_stockpile_volume.py
image2d 를 입력으로 받는 것)—
surface)dem_aspect · dem_curvature · dem_roughness · dem_tpi
*Provenance: demops.py — DEM 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.