測地系の変換 - Geodetic Datum Transformation - 麦酒堂

Willcom の端末が取得した位置情報を Google Map に投げてみたところ、測地系が異なることが判明。 Willcom は日本測地系（Tokyo）で、Google Map は世界測地系（WGS 1984）のようです。 こいつを変換してやりたいのだけれど、実はこれは非常に厄介な問題。 本気で取り組もうとするとこんなことになります。 Web版TKY2JGD 世界測地系移行のための座標変換ソフトウェア "TKY2JGD" (PDF) 多少大雑把にやろうとしてもこんな具合。 座標系変換パラメータ この値を足して補正してやる必要があります。 手順自体は簡単だけれど、どの値を使うかを条件式で判断するのが非常に萎える。 これらは、三角点測量の誤差の補正も行おうとしているためにこんなに複雑になっています。 次善の策として、準拠楕円体、測地座標系を数学的に変換するだけでお茶を濁すことにしました。 測定誤差は無視。 この変換にはいくつか方法があるようですが、今回は Molodensky変換というアルゴリズムを用います。 02 DATUM Molodensky.pl 計算に必要な準拠楕円体と測地座標系の情報はここから入手。 Geodetic Datum Overview Reference Ellipsoids and Geodetic Datum Transformation Parameters (Local to WGS-84) From NIMA 8350.2 4 July 1977 and MADTRAN 1 October 1996 この表を awk で変換して PHP の参照配列に突っ込んで利用しました。 で完成したライブラリがこちら。 GeodeticDatumTransformer.php 使い方はこんな感じです。