geodesy opなし → points(引数だけで決まる op —— 画像やデータの入力を取らない)import fullseye as fs; fs.ledger.dem_geodetic_to_ecef(lat_deg, lon_deg, height_m=0.0) (実装を直接呼ぶなら import demops; demops.dem_geodetic_to_ecef(lat_deg, lon_deg, height_m=0.0)、台帳から引くなら opsdem.get("dem_geodetic_to_ecef"))測地座標(緯度・経度・楕円体高)→ 地心直交座標 ECEF [m]。
地球の中心を原点、赤道面の本初子午線方向を x、東経 90 度を y、北極を z と する右手系。地形を「地球中心から見た座標」で扱いたいときの入口。
球ではなく WGS84 楕円体で計算する —— 球近似は緯度 45 度あたりで 最大 21 km ずれる(極半径が赤道半径より短いぶん)。
Args:
lat_deg / lon_deg: 緯度・経度 [度]。配列可(同じ形)。
height_m: 楕円体高 [m]。スカラでも配列でも可。
Returns:
(..., 3) の ECEF 座標 [m]。
py -3.11 examples/dem_geodesy_tour.pypy -3.11 examples/poc_geodetic_height_frames.pypoints を入力に取れる)geodesy)dem_ecef_to_geodetic · dem_geocentric_grid · dem_earth_curvature_drop · dem_cell_size_webmercator · dem_geodetic_slope
Provenance: demops.py — DEM operator registry. この per-op ノートは tools/opdocs.py md が自動生成(手編集しない)。
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.