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 {Number} lat

# */

def gcj02tobd09(lng, lat):

    z = math.sqrt(lng * lng + lat * lat) + 0.00002 * math.sin(lat * x_PI)

    theta = math.atan2(lat, lng) + 0.000003 * math.cos(lng * x_PI)

    bd_lng = z * math.cos(theta) + 0.0065

    bd_lat = z * math.sin(theta) + 0.006

    return [bd_lng, bd_lat]

 

# /**

# * WGS84坐标系转火星坐标系GCj02 / 即WGS84 转谷歌、高德

# * @ param {Number} lng

# * @ param {Number} lat

# */

def wgs84togcj02(lng, lat):

    if (outOfChina(lng, lat)):

        return [lng, lat]

    else:

        dlat = transformlat(lng - 105.0, lat - 35.0)

        dlng = transformlng(lng - 105.0, lat - 35.0)

        radlat = lat / 180.0 * PI

        magic = math.sin(radlat)

        magic = 1 - ee * magic * magic

        sqrtmagic = math.sqrt(magic)

        dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI)

        dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * PI)

        mglat = lat + dlat

        mglng = lng + dlng

        return [mglng, mglat]

 

# /**

# * GCJ02（火星坐标系） 转换为 WGS84 / 即谷歌高德转WGS84

# * @ param {Number} lng

# * @ param {Number} lat

# */

def gcj02towgs84(lng, lat):

    if (outOfChina(lng, lat)):

        return [lng, lat]

    else:

        dlat = transformlat(lng - 105.0, lat - 35.0)

        dlng = transformlng(lng - 105.0, lat - 35.0)

        radlat = lat / 180.0 * PI

        magic = math.sin(radlat)

        magic = 1 - ee * magic * magic

        sqrtmagic = math.sqrt(magic)

        dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI)

        dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * PI)

        mglat = lat + dlat

        mglng = lng + dlng

        return [lng * 2 - mglng, lat * 2 - mglat]

 

# /**

# * 百度坐标系转wgs84坐标系

# * @ param {*} lng

# * @ param {*} lat

# */

def bd09towgs84(lng, lat):

    # // 百度坐标系先转为火星坐标系

    gcj02 = bd09togcj02(lng, lat)

    # // 火星坐标系转wgs84坐标系

    result = gcj02towgs84(gcj02[0], gcj02[1])

    return result

 

# /**

# * wgs84坐标系转百度坐标系

# * @ param {*} lng

# * @ param {*} lat

# */

def wgs84tobd09(lng, lat):

    # // wgs84先转为火星坐标系

    gcj02 = wgs84togcj02(lng, lat)

    # // 火星坐标系转百度坐标系

    result = gcj02tobd09(gcj02[0], gcj02[1])

    return result

 

# /**

# * 经度转换

# * @ param {Number} lng

# * @ param {Number} lat

# */

def transformlat(lng, lat):

    ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * \

        lat + 0.1 * lng * lat + 0.2 * math.sqrt(abs(lng))

    ret += (20.0 * math.sin(6.0 * lng * PI) + 20.0 *

            math.sin(2.0 * lng * PI)) * 2.0 / 3.0

    ret += (20.0 * math.sin(lat * PI) + 40.0 *

            math.sin(lat / 3.0 * PI)) * 2.0 / 3.0

    ret += (160.0 * math.sin(lat / 12.0 * PI) + 320 *

            math.sin(lat * PI / 30.0)) * 2.0 / 3.0

    return ret

 

# /**

# * 纬度转换

# * @ param {Number} lng

# * @ param {Number} lat

# */

def transformlng(lng, lat):

    ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + \

        0.1 * lng * lat + 0.1 * math.sqrt(abs(lng))

    ret += (20.0 * math.sin(6.0 * lng * PI) + 20.0 *

            math.sin(2.0 * lng * PI)) * 2.0 / 3.0

    ret += (20.0 * math.sin(lng * PI) + 40.0 *

            math.sin(lng / 3.0 * PI)) * 2.0 / 3.0

    ret += (150.0 * math.sin(lng / 12.0 * PI) + 300.0 *

            math.sin(lng / 30.0 * PI)) * 2.0 / 3.0

    return ret

 

# /**

# * 判断是否在国内，不在国内则不做偏移

# * @ param {*} lng

# * @ param {*} lat

# */

def outOfChina(lng, lat):

    return (lng < 72.004 or lng > 137.8347) or ((lat < 0.8293 or lat > 55.8271) or False)