python 地图坐标系转换

import math

#  定义一些常量

x_PI = 3.14159265358979324 * 3000.0 / 180.0

PI = 3.1415926535897932384626

a = 6378245.0

ee = 0.00669342162296594323

 

# /**

#  * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换 / 即百度转谷歌、高德

#  * @param { Number } bd_lon

#  * @param { Number } bd_lat

#  */

def bd09togcj02(bd_lon, bd_lat):

    x_pi = 3.14159265358979324 * 3000.0 / 180.0

    x = bd_lon - 0.0065

    y = bd_lat - 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)

    gg_lng = z * math.cos(theta)

    gg_lat = z * math.sin(theta)

    return [gg_lng, gg_lat]

 

# /**

# * 火星坐标系(GCJ-02) 与百度坐标系(BD-09) 的转换 / 即谷歌、高德 转 百度

# * @ param {Number} lng

# * @ param {Numb