tb_dem_ecef_to_geodetic — 2D typed op

• データ種: points → points

• 呼び出し: fullseye.apply(img, "tb_dem_ecef_to_geodetic", a=0.5, b=0.5) (2-D は 1 画像 + 2 スカラつまみ a,b∈[0,1] のモデル)

*図なし: この op は points を入力に取る。画像から始まる Studio のプログラムでは型が届かないので、下の「実行できる例」で使い方を見ること。*

使い方

ECEF → 測地座標。返りは `(..., 3) の (緯度[度], 経度[度], 高さ[m])`。

Bowring (1976) の閉形式に近い解法。往復(測地→ECEF→測地)の誤差は

実測で緯度・経度が 1e-12 度未満、高さが 1e-7 m 未満。

地心球座標が欲しいだけなら、`r = |xyz|` と

`geocentric_lat = asin(z/r)` で足りる —— ただしそれは測地緯度ではない

(両者は最大 0.19 度、距離にして約 21 km ずれる)。この op が返すのは

地図や GPS と同じ測地緯度のほう。

2-D 進化レジストリへ橋渡しした dem の op `dem_ecef_to_geodetic。実装は同じで、呼び出し規約だけ op(v, a, b) に合わせてある。この op に調整点は無く、a も b` も使われない。

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

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

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

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

• (まだありません)

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

identity · tb_points_to_voxel · tb_estimate_point_normals · tb_iss_keypoints · tb_angle_3points · tb_project_points · tb_render_point_depth · tb_statistical_outlier_removal

同カテゴリ(typed)

tb_points_to_voxel · tb_estimate_point_normals · tb_iss_keypoints · tb_angle_3points · tb_project_points · tb_render_point_depth · tb_statistical_outlier_removal · tb_radius_outlier_removal

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

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.