dem_aspect — DEM `surface` op

• 데이터 종류: depth → image2d

• 호출: import fullseye as fs; fs.ledger.dem_aspect(dem, cell_size, method='horn', flat_tol=1e-12)(구현을 직접 호출하려면 import demops; demops.dem_aspect(dem, cell_size, method='horn', flat_tol=1e-12), 원장에서 가져오려면 opsdem.get("dem_aspect"))

사용법

> 이 연산자의 설명은 아직 번역이 없습니다. 원문을 그대로 싣습니다.

斜面方位 [度]。北 0 度・東回り。平坦なセルは :data:ASPECT_FLAT。

平坦を 0 で返さないのは、0 が「北向き斜面」を意味してしまうため

—— 区別できない値を返すのは、例外を出さずに間違える典型。

手順: `dem_slope と同じ 3x3 勾配 (dz/dx, dz/dy)` を取り、下り勾配の

向き `atan2(-dz/dx, -dz/dy) を度にして [0, 360)` に折り返す。行 0 が北・

列が増える向きが東という前提で、北 0 / 東 90 / 南 180 / 西 270。

行 0 が南の格子(上下反転した配列)を渡すと南北が入れ替わり、例外は出ない。

• `dem, cell_size, method: dem_slope と同じ契約(3x3 以上、nan`

で欠測、番兵値は拒否、`cell_size` は正の [m])。

• `flat_tol: 勾配の大きさ hypot(dz/dx, dz/dy)` [m/m] がこれ以下なら平坦と

みなし `ASPECT_FLAT`(= -1.0)を返す。既定 1e-12 は「数値的に厳密に平坦」だけ。

整数標高の平地でも微小勾配で方位が付くことがあるので、意味のある閾

(例 1e-3)を明示するとよい。

• 返り値: `(H, W) float64。[0, 360) または -1.0。nan セルの周りは nan`。

• 方位のヒストグラムや平均を取るときは `-1 を先に除く(out >= 0`)。

자세한 사용 가이드

• dem_terrain_analysis 패밀리 가이드

참고(샘플 데이터·문헌)

• 샘플 데이터 카탈로그(DL URL / 라이선스) —— 2-D 는 skimage.data(BSD/public)+ 합성, 3-D 는 실데이터 소스(Stanford/PDS 등)의 DL URL.

• 연산자의 내력·참고문헌 —— 이 연산자 족의 바탕이 된 연구/기법의 출처.

• 알고리즘의 정전(저자·연도)과 용도는 위의 패밀리 사용 가이드에 적혀 있습니다.

실행 가능한 예제(이 연산자를 실제로 호출하는 검증된 샘플)

• poc_dem_terrain — py -3.11 examples/poc_dem_terrain.py

• poc_geodetic_height_frames — py -3.11 examples/poc_geodetic_height_frames.py

타입이 이어지는 다음 연산자(`image2d` 를 입력으로 받는 것)

—

같은 카테고리(`surface`)

dem_slope · dem_curvature · dem_roughness · dem_tpi

*Provenance: demops.py — DEM 연산자 레지스트리. 이 op 노트는 tools/opdocs.py md 가 자동 생성합니다(직접 편집하지 마세요).*

dem_aspect — DEM surface op

사용법