Like Silk: Script: cnvgps

python関連のメモ

GPSの測地系変換(日本測地系:Tokyo Datum ⇔ 世界測地系:WGS84)

日本測地系と世界測地系でGPSの緯度経度が若干異なる. 変換用のスクリプトです.変換公式は,座標(緯度経度)の表示方法と変換方法 - end0tknrのkipple - web写経開発のコピペです.理由はわかりません.

code

cnvgps.py←右クリックで保存できます.
#!/usr/bin/env python
#coding: utf-8

def cnvTokyoToWGS84(lati,long):
	lati2 = lati -lati * 0.00010695  + long * 0.000017464 + 0.0046017
	long2 = long -lati * 0.000046038 - long * 0.000083043 + 0.010040
	return lati2, long2

def cnvWGS84ToTokyo(lati,long):
	lati2 = lati +lati * 0.00010696  - long * 0.000017467 - 0.0046020
	long2 = long +latg * 0.000046047 + long * 0.000083049 - 0.010041
	return lati2, long2


########################################################################
if __name__ == '__main__':
	# Tokyo station in Tokyo Datum
	lati = 35.67808267
	long = 139.7704438
	print cnvTokyoToWGS84(lati,long)

参考URL


Like Silk