坐标转换工具：WGS84 / GCJ02 / 百度 BD09 经纬度在线互转

• 十进制度 DD

  116.40337249, 39.91793075

  经度在前，GIS 与多数 API 默认顺序

• 纬度, 经度

  39.91793075, 116.40337249

  谷歌地图搜索框格式

• 度分秒 DMS

  116°24′12.141″E 39°55′04.551″N

  测绘成果、航海常用

• 度分 DM

  116°24.20235′E 39°55.07584′N

  NMEA / 航空常用

• WKT

  POINT(116.40337249 39.91793075)

  PostGIS / OGC 标准文本

• GeoJSON

  {"type":"Point","coordinates":[116.40337249,39.91793075]}

  RFC 7946，坐标顺序为 [经度, 纬度]

# GCJ-02 / BD-09 无官方 EPSG,以下实现与在线工具一致
import math

_A = 6378245.0
_EE = 0.00669342162296594323

def _out_of_china(lng, lat):
    return not (72.004 <= lng <= 137.8347 and 0.8293 <= lat <= 55.8271)

def _tr_lat(x, y):
    r = -100.0 + 2.0*x + 3.0*y + 0.2*y*y + 0.1*x*y + 0.2*math.sqrt(abs(x))
    r += (20.0*math.sin(6.0*x*math.pi) + 20.0*math.sin(2.0*x*math.pi)) * 2.0/3.0
    r += (20.0*math.sin(y*math.pi) + 40.0*math.sin(y/3.0*math.pi)) * 2.0/3.0
    r += (160.0*math.sin(y/12.0*math.pi) + 320*math.sin(y*math.pi/30.0)) * 2.0/3.0
    return r

def _tr_lng(x, y):
    r = 300.0 + x + 2.0*y + 0.1*x*x + 0.1*x*y + 0.1*math.sqrt(abs(x))
    r += (20.0*math.sin(6.0*x*math.pi) + 20.0*math.sin(2.0*x*math.pi)) * 2.0/3.0
    r += (20.0*math.sin(x*math.pi) + 40.0*math.sin(x/3.0*math.pi)) * 2.0/3.0
    r += (150.0*math.sin(x/12.0*math.pi) + 300.0*math.sin(x/30.0*math.pi)) * 2.0/3.0
    return r

def wgs84_to_gcj02(lng, lat):
    if _out_of_china(lng, lat):
        return lng, lat
    d_lat = _tr_lat(lng - 105.0, lat - 35.0)
    d_lng = _tr_lng(lng - 105.0, lat - 35.0)
    rad = lat / 180.0 * math.pi
    magic = 1 - _EE * math.sin(rad) ** 2
    sqrt_magic = math.sqrt(magic)
    d_lat = (d_lat * 180.0) / ((_A * (1 - _EE)) / (magic * sqrt_magic) * math.pi)
    d_lng = (d_lng * 180.0) / (_A / sqrt_magic * math.cos(rad) * math.pi)
    return lng + d_lng, lat + d_lat

def gcj02_to_wgs84(lng, lat, eps=1e-11, max_iter=30):
    """迭代精确反算,残差 < 1e-11°(约 1 μm),优于常见的 2*P - f(P) 近似"""
    if _out_of_china(lng, lat):
        return lng, lat
    w_lng, w_lat = lng, lat
    for _ in range(max_iter):
        t_lng, t_lat = wgs84_to_gcj02(w_lng, w_lat)
        d_lng, d_lat = lng - t_lng, lat - t_lat
        w_lng += d_lng
        w_lat += d_lat
        if abs(d_lng) < eps and abs(d_lat) < eps:
            break
    return w_lng, w_lat

_BD = math.pi * 3000.0 / 180.0

def gcj02_to_bd09(lng, lat):
    z = math.sqrt(lng*lng + lat*lat) + 0.00002 * math.sin(lat * _BD)
    theta = math.atan2(lat, lng) + 0.000003 * math.cos(lng * _BD)
    return z * math.cos(theta) + 0.0065, z * math.sin(theta) + 0.006

def bd09_to_gcj02(lng, lat):
    # 近似反算作初值,再用正算残差迭代,残差 < 1e-11°
    x, y = lng - 0.0065, lat - 0.006
    z = math.sqrt(x*x + y*y) - 0.00002 * math.sin(y * _BD)
    theta = math.atan2(y, x) - 0.000003 * math.cos(x * _BD)
    g_lng, g_lat = z * math.cos(theta), z * math.sin(theta)
    for _ in range(10):
        f_lng, f_lat = gcj02_to_bd09(g_lng, g_lat)
        d_lng, d_lat = lng - f_lng, lat - f_lat
        g_lng, g_lat = g_lng + d_lng, g_lat + d_lat
        if abs(d_lng) < 1e-12 and abs(d_lat) < 1e-12:
            break
    return g_lng, g_lat

# WGS84 → GCJ-02
lng, lat = 116.397128, 39.916527
lng, lat = wgs84_to_gcj02(lng, lat)
print(lng, lat)