東京測地系から世界測地系への変換 - eggshell blue最近新幹線に乗って地方へ出かける機会が多く、 GPS トラッキングとかして遊んでるのですが、 GPS ユニットの出力ログが日本測地系なので Google Maps 上で表示させるには測地系を変換する必要が出てきました。 なので、測地系パラメタを使う方法 を参考に、三次元直交座標を経由して変換する方法で実装してみました。ページにあったのは Perl だったので、今回は Python で書き直してみた。ていうか Python 3.0 系って 2.x 系から文法まで変わってるのな。びっくらこいた。 # coding: utf-8 import math # Radian RAD = math.pi / 180 # 日本測地系における赤道半径、扁平率、第一離心率 TKY = 6377397.155 tireT = 1/299.152813 ecT = tireT * (2 - tireT) # 世
