话不多说,直接开撸
#!/usr/bin/env python3
# -*- coding: utf-8 -*-
#转换经纬度
#高德地图,腾讯地图用的GCJ02坐标
#百度地图用的BD09坐标
import math
#----------------GCJ02坐标 --->BD09坐标
def gcj2_to_bd09(tx_jd, tx_wd):
x_pi = math.pi*3000.0/180.0
x = tx_jd
y = tx_wd
z = math.sqrt(x*x+y*y)+0.00002*math.sin(y*x_pi)
theta = math.atan2(y, x)+0.000003*math.cos(x*x_pi)
jd_new = z*math.cos(theta)+0.0065
wd_new = z*math.sin(theta)+0.006
print("百度地图纬度:"+'%.6f' % wd_new)
print("百度地图经度:"+'%.6f' % jd_new)
#----------------BD09坐标--->GCJ02坐标
def bd09_to_gcj2(bd_jd, bd_wd):
x_pi = math.pi*3000.0/180.0
x = bd_jd-0.0065
y = bd_wd-0.006
z = math.sqrt(x*x+y*y)-0.00002*math.sin(y*x_pi)
theta = math.atan2(y, x)-0.000003*math.cos(x*x_pi)
jd_new = z*math.cos(theta)
wd_new = z*math.sin(theta)
print("高德/腾讯地图纬度:"+'%.6f' % wd_new)
print("高德/腾讯地图经度:"+'%.6f' % jd_new)
wd = 22.569111
jd = 114.117024
bd09_to_gcj2(jd, wd)
#gcj2_to_bd09(jd, wd)