dem_geocentric_grid — DEM `geodesy` op

• データ種: depth → pointmap

• 呼び出し: import demops; demops.dem_geocentric_grid(dem, lat0_deg, lon0_deg, cell_size, spherical=False) (または opsdem.get("dem_geocentric_grid"))

使い方

DEM の各セルを地球中心から見た座標にする。

`(H, W) の標高格子を、北西角が (lat0_deg, lon0_deg)` にある局所平面と

みなし、各セルを ECEF(既定)または地心球座標へ写す。

Args:

dem: `(H, W)` の標高 [m]。

lat0_deg / lon0_deg: 格子の北西角の緯度経度 [度]。

cell_size: セル寸法 [m]。緯度方向・経度方向とも同じとみなす

(Web メルカトルのタイルはそうなっている。

:func:dem_cell_size_webmercator を参照)。

spherical: `True なら (半径 r[m], 地心緯度[度], 経度[度])` を返す。

`False(既定)なら ECEF (x, y, z)` [m]。

Returns:

`(H, W, 3)`。

★ 地心緯度は測地緯度ではありません。`spherical=True` の返りの緯度は

地球中心から見た角度で、地図の緯度(楕円体の法線が赤道面となす角)とは

最大 0.19 度違います。地図に戻すときは ECEF 側を

:func:dem_ecef_to_geodetic に渡してください。

詳しい使い方ガイド

• dem_terrain_analysis ファミリガイド

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

• サンプルデータカタログ(DL URL / ライセンス) — 2-D は skimage.data(BSD/public)+ 合成、3-D は実データ源(Stanford/PDS 等)の DL URL。

• 演算子の来歴・参考文献 — この op 族の元になった研究/手法の出典。

• アルゴリズムの正典(著者・年)と用途は上記ファミリ使い方ガイドに記載。

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

• (まだありません)

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

—

同カテゴリ(`geodesy`)

dem_geodetic_to_ecef · dem_ecef_to_geodetic · dem_earth_curvature_drop · dem_cell_size_webmercator · dem_geodetic_slope

*Provenance: demops.py — DEM operator registry. この per-op ノートは tools/opdocs.py md が自動生成(手編集しない)。*

dem_geocentric_grid — DEM geodesy op

使い方