dem_geocentric_grid — DEM `geodesy` op

• データ種: depth → pointmap

• 呼び出し: import fullseye as fs; fs.ledger.dem_geocentric_grid(dem, lat0_deg, lon0_deg, cell_size, spherical=False) (実装を直接呼ぶなら 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 メルカトルのタイルはそうなっている。

dem_cell_size_webmercator を参照)。

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

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

Returns:

(H, W, 3)。

★ 地心緯度は測地緯度ではありません。spherical=True の返りの緯度は地球中心から見た角度で、地図の緯度(楕円体の法線が赤道面となす角)とは最大 0.19 度違います。地図に戻すときは ECEF 側を dem_ecef_to_geodetic に渡してください。

詳しい使い方ガイド

• 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(`pointmap` を入力に取れる)

—

同カテゴリ(`geodesy`)

dem_geodetic_to_ecef · dem_ecef_to_geodetic · dem_earth_curvature_drop · dem_cell_size_webmercator · dem_geodetic_slope · dem_geoid_height · dem_height_frame_convert · dem_height_frame_residual

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

dem_geocentric_grid — DEM geodesy op

使い方