百度服务获取坐标放置在天地图上实例

开始之前还是要提前写一下涉及到的几个坐标系

WGS-84：是国际标准，GPS坐标（Google Earth使用、或者GPS模块）
GCJ-02：中国坐标偏移标准，Google Map、高德、腾讯使用
BD-09：百度坐标偏移标准，Baidu Map使用

CGCS2000：天地图使用

---

先以北京华宾公司为例，地址是“北京市西城区西安门大街107号”，百度地图上是如下位置

访问百度WEB服务API，找到“正/逆地理编码”-“地理编码”，可以用这个输入上面公司的地址，返回坐标。

调用接口如下，得到坐标

116.38426468294915,39.928903206383797

访问百度API示例，查看定位点

可以看到是有一点偏差的。

---

下面将百度的BD-09坐标转成WGS84（CGCS2000可以看作和WGS84“毫无差别”），放到天地图上。

百度地图没有提供直转的方法，因此需要将BD-09先转成GCJ-02，再转WGS84。

回到百度WEB服务API，找到坐标转换，根据参数说明，from为5，to为3。

结果为116.37782967212564,39.92277118837159。调用结果如下：

GCJ-02坐标转WGS84网上很容易查到，也放在这里吧

var PI = 3.14159265358979324;
function transformGCJ2WGS(gcjLat, gcjLon) {
    let d = delta(gcjLat, gcjLon)
    return {
        'lat': gcjLat - d.lat,
        'lon': gcjLon - d.lon
    }
}

function delta(lat, lon) {
    let a = 6378245.0 //  a: 卫星椭球坐标投影到平面地图坐标系的投影因子。
    let ee = 0.00669342162296594323 //  ee: 椭球的偏心率。
    let dLat = transformLat(lon - 105.0, lat - 35.0)
    let dLon = transformLon(lon - 105.0, lat - 35.0)
    let radLat = lat / 180.0 * 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) * PI)
    dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * PI)
    return {
        'lat': dLat,
        'lon': dLon
    }
}
function transformLat(x, y) {
    let ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x))
    ret += (20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0 / 3.0
    ret += (20.0 * Math.sin(y * PI) + 40.0 * Math.sin(y / 3.0 * PI)) * 2.0 / 3.0
    ret += (160.0 * Math.sin(y / 12.0 * PI) + 320 * Math.sin(y * PI / 30.0)) * 2.0 / 3.0
    return ret
}
function transformLon(x, y) {
    let ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.sqrt(Math.abs(x))
    ret += (20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0 / 3.0
    ret += (20.0 * Math.sin(x * PI) + 40.0 * Math.sin(x / 3.0 * PI)) * 2.0 / 3.0
    ret += (150.0 * Math.sin(x / 12.0 * PI) + 300.0 * Math.sin(x / 30.0 * PI)) * 2.0 / 3.0
    return ret
}

调用transformGCJ2WGS（注意参数是纬度在前）方法，返回116.37159805544552,39.921379265627216。分别将上述坐标放到天地图demo中，查看如下：

对比图（天地图-百度地图）

由此可见，偏差不大。鉴于时间不够，只能选取一个点来测试，有时间再选其他点尝试下。