WGS-84 转高德、百度、腾讯地图坐标系

正文与代码来自原稿。GPS 设备多为 WGS-84，国内在线底图需加密坐标系。

1、坐标系简介

坐标系	说明	常见用途
WGS-84	GPS 原始	设备上报、国际地图
GCJ-02	国测局加密	高德、腾讯
BD-09	百度再加密	百度地图

示例：WGS-84 39.908823, 116.397470，转到百度后偏移通常更大（在 GCJ-02 上再偏一次）。

关系：

• WGS-84 → GCJ-02：公开加密算法
• GCJ-02 → BD-09：百度算法
• 反向也有公开实现，但厂商不鼓励随意逆变换

大陆展示地图一般不允许直接用未加密 WGS-84 叠官方瓦片。

2、转换代码（JavaScript）

WGS-84 → GCJ-02（高德 / 腾讯）

function wgs84ToGcj02(lng, lat) {
  if (outOfChina(lng, lat)) return [lng, lat];
  let a = 6378245.0;
  let ee = 0.00669342162296594323;
  let dLat = transformLat(lng - 105.0, lat - 35.0);
  let dLng = transformLng(lng - 105.0, lat - 35.0);
  let radLat = lat / 180.0 * Math.PI;
  let magic = Math.sin(radLat);
  magic = 1 - ee * magic * magic;
  let sqrtMagic = Math.sqrt(magic);
  dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * Math.PI);
  dLng = (dLng * 180.0) / (a / sqrtMagic * Math.cos(radLat) * Math.PI);
  return [lng + dLng, lat + dLat];
}

function transformLat(x, y) {
  return -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y +
    0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x)) +
    (20.0 * Math.sin(6.0 * x * Math.PI) + 20.0 * Math.sin(2.0 * x * Math.PI)) * 2.0 / 3.0 +
    (20.0 * Math.sin(y * Math.PI) + 40.0 * Math.sin(y / 3.0 * Math.PI)) * 2.0 / 3.0 +
    (160.0 * Math.sin(y / 12.0 * Math.PI) + 320 * Math.sin(y * Math.PI / 30.0)) * 2.0 / 3.0;
}

function transformLng(x, y) {
  return 300.0 + x + 2.0 * y + 0.1 * x * x +
    0.1 * x * y + 0.1 * Math.sqrt(Math.abs(x)) +
    (20.0 * Math.sin(6.0 * x * Math.PI) + 20.0 * Math.sin(2.0 * x * Math.PI)) * 2.0 / 3.0 +
    (20.0 * Math.sin(x * Math.PI) + 40.0 * Math.sin(x / 3.0 * Math.PI)) * 2.0 / 3.0 +
    (150.0 * Math.sin(x / 12.0 * Math.PI) + 300.0 * Math.sin(x / 30.0 * Math.PI)) * 2.0 / 3.0;
}

function outOfChina(lng, lat) {
  return (lng < 72.004 || lng > 137.8347 || lat < 0.8293 || lat > 55.8271);
}

GCJ-02 → BD-09

function gcj02ToBd09(lng, lat) {
  let z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * Math.PI * 3000.0 / 180.0);
  let theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * Math.PI * 3000.0 / 180.0);
  return [z * Math.cos(theta) + 0.0065, z * Math.sin(theta) + 0.006];
}

function wgs84ToBd09(lng, lat) {
  const [gcjLng, gcjLat] = wgs84ToGcj02(lng, lat);
  return gcj02ToBd09(gcjLng, gcjLat);
}

示例

const wgs84 = [116.397128, 39.916527];
console.log("GCJ-02:", wgs84ToGcj02(...wgs84));
console.log("BD-09:", wgs84ToBd09(...wgs84));
// 约: GCJ-02 [116.40337, 39.91793]  BD-09 [116.40974, 39.92427]

3、小结

规则链或自定义部件里，上报后先转换再存/再画图；与 地图绘制路线 配合时，底图厂商与转换目标保持一致。