fullseye

dem_geodetic_slope — DEM `geodesy` op

データ種: depth → image2d
呼び出し: import fullseye as fs; fs.ledger.dem_geodetic_slope(dem, lat0_deg, d_lat_deg, d_lon_deg, method='horn', units='degrees') (実装を直接呼ぶなら import demops; demops.dem_geodetic_slope(dem, lat0_deg, d_lat_deg, d_lon_deg, method='horn', units='degrees')、台帳から引くなら opsdem.get("dem_geodetic_slope"))

使い方

緯度経度の格子(等角度間隔)の 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 と同じ。

詳しい使い方ガイド

dem_terrain_analysis ファミリガイド

参考(サンプルデータ・文献)

サンプルデータカタログ(DL URL / ライセンス) — 2-D は skimage.data(BSD/public)+ 合成、3-D は実データ源(Stanford/PDS 等)の DL URL。
演算子の来歴・参考文献 — この op 族の元になった研究/手法の出典。
アルゴリズムの正典(著者・年)と用途は上記ファミリ使い方ガイドに記載。

実行できる例(この op を実際に呼ぶ検証済みサンプル)

dem_geodesy_tour — py -3.11 examples/dem_geodesy_tour.py

型が繋がる次の op(`image2d` を入力に取れる)

—

同カテゴリ(`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. この per-op ノートは tools/opdocs.py md が自動生成(手編集しない)。

This site is open source. Improve this page.