dem_geodetic_slope — DEM `geodesy` op

• Data kinds: depth → image2d

• Call: import demops; demops.dem_geodetic_slope(dem, lat0_deg, d_lat_deg, d_lon_deg, method='horn', units='degrees') (or opsdem.get("dem_geodetic_slope"))

Usage

> This operator's description has not been translated yet. The original text follows as it is.

緯度経度の格子(等角度間隔)の DEM の傾斜。セル寸法が緯度で変わる。

公開 DEM の多くは「1 秒メッシュ」のように角度で等間隔です。この格子を

そのまま :func:dem_slope に一定のセル寸法で渡すと、経度方向の実距離が

`cos(緯度)` 倍だけ短いことが無視されます。緯度方向の寸法を両軸に使うと

東西の傾斜が過小になり(北緯 60 度で半分)、経度方向の寸法を両軸に使えば

今度は南北が過大になります。どちらに転ぶかは何を定数にしたかで決まるので、

「過大/過小」を覚えるのではなく緯度ごとに寸法を計算するのが正解です。

★ 最初この docstring に「東西が過大になる」と書いたが、テストで測ったら

半分になった(過小)。符号や向きの主張は測ってから書くこと。

ここでは緯度ごとに東西のセル寸法を計算し、格子を等距離へ直してから測ります。

Args:

dem: `(H, W)`。行 0 が北。

lat0_deg: 北西角の緯度 [度]。

d_lat_deg / d_lon_deg: 1 セルあたりの緯度・経度の刻み [度]。正の値。

method / units: :func:dem_slope と同じ。

Detailed usage guide

• dem_terrain_analysis family guide

References (sample data, literature)

• 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.

Runnable examples (verified samples that actually call this op)

• (none yet)

Ops the type connects to (they accept `image2d` as input)

—

Same category (`geodesy`)

dem_geodetic_to_ecef · dem_ecef_to_geodetic · dem_geocentric_grid · dem_earth_curvature_drop · dem_cell_size_webmercator

*Provenance: demops.py — DEM operator registry. This per-op note is generated by tools/opdocs.py md (do not hand-edit).*

dem_geodetic_slope — DEM geodesy op