专业坐标系转换工具

坐标系选择

坐标输入支持多种格式：十进制度、度分秒、度分

跨系统转换一律以 WGS84 为中枢。逐步核对中间值，可以直接定位是哪一段引入了偏移。

1. 1
   WGS84 → GCJ-02国测局加偏算法
   116.41024450, 39.91640428

   仅在中国境内生效,境外原值返回

在你的代码里复现这次转换

片段使用与本页完全相同的算法与参数，可直接粘贴运行。

# 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):
    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)
    return z * math.cos(theta), z * math.sin(theta)

# WGS84 → GCJ-02
lng, lat = 116.404, 39.915
lng, lat = wgs84_to_gcj02(lng, lat)
print(lng, lat)

使用您本地测绘成果的 Bursa-Wolf 七参数进行 WGS84 与地方坐标系间的高精度转换，参数可保存为预设。