常用的几种在线地图（天地图、百度地图、高德地图）坐标系之间的转换算法

最新推荐文章于 2024-03-26 10:40:34 发布
最新推荐文章于 2024-03-26 10:40:34 发布
阅读量8.5w 收藏 342
点赞数 72
分类专栏： 坐标系学习 文章标签： 常用在线地图之间坐标转换 天地图坐标转百度坐标 墨卡托坐标转经纬度坐标 坐标转换算法js版
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_36377037/article/details/86479796
版权

1、首先弄明白几种在线地图的坐标系；
（1）天地图：CGCS2000，2000国家大地坐标系；我们其实很多时候直接用WGS84的坐标来代替CGCS2000坐标。因为CGCS2000的定义与WGS84实质一样。采用的参考椭球非常接近。扁率差异引起椭球面上的纬度和高度变化最大达0.1mm。当前测量精度范围内，可以忽略这点差异。可以说两者相容至cm级水平，但若一点的坐标精度达不到cm水平，则不认为CGCS2000和WGS84的坐标是相容的。

（2）百度地图：bd09II坐标。首先了解一下火星坐标，它是在国际标准坐标WGS-84上进行的一次加密，由于国内的电子地图都要至少使用火星坐标进行一次加密，百度直接就任性一些，直接自己又研究了一套加密算法，来了个二次加密，这就是我们所熟知的百度坐标(BD-09)。

（3）高德地图：gcj02坐标，也称为火星坐标。火星坐标是国家测绘局为了国家安全在原始坐标的基础上进行偏移得到的坐标，基本国内的电子地图、导航设备都是采用的这一坐标系或在这一坐标的基础上进行二次加密得到的。

（4）一些注意点：非专业人士，对于墨卡托投影坐标和经纬度坐标可以理解为一个是平面坐标一个是球面坐标。还有用wgs84坐标系的一些地图，比如： osm、谷歌、arcgisonline等。

2、几种坐标系之间的转换算法；
（1）百度坐标转火星坐标：

var x_pi=3.14159265358979324 * 3000.0 / 180.0;
function baiduTomars(baidu_point){
    var mars_point={lon:0,lat:0};
    var x=baidu_point.lon-0.0065;
    var y=baidu_point.lat-0.006;
    var z=Math.sqrt(x*x+y*y)- 0.00002 * Math.sin(y * x_pi);
    var theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_pi);
    mars_point.lon=z * Math.cos(theta);
    mars_point.lat=z * Math.sin(theta);
    return mars_point;
}

（2）火星坐标转百度坐标：

x_pi=3.14159265358979324 * 3000.0 / 180.0;
function marsTobaidu(mars_point){
    var baidu_point={lon:0,lat:0};
    var x=mars_point.lon;
    var y=mars_point.lat;
    var z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * x_pi);
    var theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * x_pi);
    baidu_point.lon = z * Math.cos(theta) + 0.0065;
    baidu_point.lat = z * Math.sin(theta) + 0.006;
    return baidu_point;
}

（3）地球坐标系(WGS-84)转火星坐标系(GCJ)：

var pi = 3.14159265358979324;
var a = 6378245.0;
var ee = 0.00669342162296594323;
/*判断是否在国内，不在国内则不做偏移*/
function outOfChina(lon, lat)
{
    if ((lon < 72.004 || lon > 137.8347)&&(lat < 0.8293 || lat > 55.8271)){
        return true;
    }else {
        return false;
    }
}
function transformLat(x,y)
{
    var 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)
{
    var 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;
}

function transform(wgLat,wgLon)
{
    var mars_point={lon:0,lat:0};
    if (outOfChina(wgLat, wgLon))
    {
        mars_point.lat = wgLat;
        mars_point.lon = wgLon;
        return;
    }
    var dLat = transformLat(wgLon - 105.0, wgLat - 35.0);
    var dLon = transformLon(wgLon - 105.0, wgLat - 35.0);
    var radLat = wgLat / 180.0 * pi;
    var magic = Math.sin(radLat);
    magic = 1 - ee * magic * magic;
    var sqrtMagic = Math.sqrt(magic);
    dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);
    dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);
    mars_point.lat = wgLat + dLat;
    mars_point.lon = wgLon + dLon;
    return mars_point
}

（4）火星坐标系GCJ02转地球坐标系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
  }

（5）百度经纬度坐标转百度墨卡托坐标：

    var LLBAND =[75, 60, 45, 30, 15, 0]
    var LL2MC = [
        [-0.0015702102444, 111320.7020616939, 1704480524535203, -10338987376042340, 26112667856603880, -35149669176653700, 26595700718403920, -10725012454188240, 1800819912950474, 82.5],
        [0.0008277824516172526, 111320.7020463578, 647795574.6671607, -4082003173.641316, 10774905663.51142, -15171875531.51559, 12053065338.62167, -5124939663.577472, 913311935.9512032, 67.5],
        [0.00337398766765, 111320.7020202162, 4481351.045890365, -23393751.19931662, 79682215.47186455, -115964993.2797253, 97236711.15602145, -43661946.33752821, 8477230.501135234, 52.5],
        [0.00220636496208, 111320.7020209128, 51751.86112841131, 3796837.749470245, 992013.7397791013, -1221952.21711287, 1340652.697009075, -620943.6990984312, 144416.9293806241, 37.5],
        [-0.0003441963504368392, 111320.7020576856, 278.2353980772752, 2485758.690035394, 6070.750963243378, 54821.18345352118, 9540.606633304236, -2710.55326746645, 1405.483844121726, 22.5],
        [-0.0003218135878613132, 111320.7020701615, 0.00369383431289, 823725.6402795718, 0.46104986909093, 2351.343141331292, 1.58060784298199, 8.77738589078284, 0.37238884252424, 7.45]
    ]
    function getRange(cC, cB, T) {
        if (cB != null) {
            cC = Math.max(cC, cB);
        }
        if (T != null) {
            cC = Math.min(cC, T);
        }
        return cC;
    }
    function getLoop(cC, cB, T) {
        while (cC > T) {
            cC -= T - cB;
        }
        while (cC < cB) {
            cC += T - cB;
        }
        return cC;
    }
    function convertor(cC, cD) {
      if (!cC || !cD) {
        return null;
      }
      let T = cD[0] + cD[1] * Math.abs(cC.x);
      const cB = Math.abs(cC.y) / cD[9];
      let cE = cD[2] + cD[3] * cB + cD[4] * cB * cB +
        cD[5] * cB * cB * cB + cD[6] * cB * cB * cB * cB +
        cD[7] * cB * cB * cB * cB * cB +
        cD[8] * cB * cB * cB * cB * cB * cB;
      T *= (cC.x < 0 ? -1 : 1);
      cE *= (cC.y < 0 ? -1 : 1);
      return [T, cE];
    }
    function convertLL2MC(T) {
        let cD, cC, len;
        T.x = getLoop(T.x, -180, 180);
        T.y = getRange(T.y, -74, 74);
        const cB = T;
        for (cC = 0, len = LLBAND.length; cC < len; cC++) {
            if (cB.y >= LLBAND[cC]) {
                cD = LL2MC[cC];
                break;
            }
        }
        if (!cD) {
            for (cC = LLBAND.length - 1; cC >= 0; cC--) {
                if (cB.y <= -LLBAND[cC]) {
                    cD = LL2MC[cC];
                    break;
                }
            }
        }
        const cE = convertor(T, cD);
        return cE;
    }

（6）百度墨卡托坐标转百度经纬度坐标：

var MCBAND = [12890594.86, 8362377.87, 5591021, 3481989.83, 1678043.12, 0]
    var MC2LL = [
      [1.410526172116255e-8, 0.00000898305509648872, -1.9939833816331, 200.9824383106796, -187.2403703815547, 91.6087516669843, -23.38765649603339, 2.57121317296198, -0.03801003308653, 17337981.2],
      [-7.435856389565537e-9, 0.000008983055097726239, -0.78625201886289, 96.32687599759846, -1.85204757529826, -59.36935905485877, 47.40033549296737, -16.50741931063887, 2.28786674699375, 10260144.86],
      [-3.030883460898826e-8, 0.00000898305509983578, 0.30071316287616, 59.74293618442277, 7.357984074871, -25.38371002664745, 13.45380521110908, -3.29883767235584, 0.32710905363475, 6856817.37],
      [-1.981981304930552e-8, 0.000008983055099779535, 0.03278182852591, 40.31678527705744, 0.65659298677277, -4.44255534477492, 0.85341911805263, 0.12923347998204, -0.04625736007561, 4482777.06],
      [3.09191371068437e-9, 0.000008983055096812155, 0.00006995724062, 23.10934304144901, -0.00023663490511, -0.6321817810242, -0.00663494467273, 0.03430082397953, -0.00466043876332, 2555164.4],
      [2.890871144776878e-9, 0.000008983055095805407, -3.068298e-8, 7.47137025468032, -0.00000353937994, -0.02145144861037, -0.00001234426596, 0.00010322952773, -0.00000323890364, 826088.5]
    ]
    function convertor(cC, cD) {
      if (!cC || !cD) {
        return null;
      }
      let T = cD[0] + cD[1] * Math.abs(cC.x);
      const cB = Math.abs(cC.y) / cD[9];
      let cE = cD[2] + cD[3] * cB + cD[4] * cB * cB +
        cD[5] * cB * cB * cB + cD[6] * cB * cB * cB * cB +
        cD[7] * cB * cB * cB * cB * cB +
        cD[8] * cB * cB * cB * cB * cB * cB;
      T *= (cC.x < 0 ? -1 : 1);
      cE *= (cC.y < 0 ? -1 : 1);
      return [T, cE];
    }
    
    function convertMC2LL(cB) {
      const cC = {
        x: Math.abs(cB.x),
        y: Math.abs(cB.y)
      };
      let cE;
      for (let cD = 0, len = MCBAND.length; cD < len; cD++) {
        if (cC.y >= MCBAND[cD]) {
          cE = MC2LL[cD];
          break;
        }
      }
      const T = convertor(cB, cE);
      return T;
    }

（7）经纬度坐标转墨卡托坐标（不太精确）：

function lonlatTomercator(lonlat) {
    var mercator={x:0,y:0};
    var x = lonlat.x *20037508.34/180;
    var y = Math.log(Math.tan((90+lonlat.y)*Math.PI/360))/(Math.PI/180);
    y = y *20037508.34/180;
    mercator.x = x;
    mercator.y = y;
    return mercator ;
}

（8）墨卡托坐标转经纬度坐标（不太精确）：

function mercatorTolonlat(mercator){
    var lonlat={x:0,y:0};
    var x = mercator.x/20037508.34*180;
    var y = mercator.y/20037508.34*180;
    y= 180/Math.PI*(2*Math.atan(Math.exp(y*Math.PI/180))-Math.PI/2);
    lonlat.x = x;
    lonlat.y = y;
    return lonlat;
}

3、感谢下面两位博主的分享：
博客地址：
（1）https://blog.csdn.net/qq_35771567/article/details/84766271
（2）https://www.cnblogs.com/fwc1994/p/5884115.html

四域公子
关注
  • 72
    点赞
  • 踩
  • 342
    收藏
    觉得还不错? 一键收藏
  • 12
    评论
专栏目录
Android中GPS坐标转换为高德地图坐标详解
01-05
地图坐标大致分为几种：  1、GPS、WGS84，也就是原始坐标体系，这是国际公认的世界标准坐标体系；  2、GCJ-02，又称为“火星坐标”，国家测绘局在02年发布的坐标体系，在国内，至少得使用此坐标体系，比如：google...
高德、百度坐标系转天地图坐标系
06-06
高德、百度坐标系转天地图坐标系，亲测实际有效，误差很小。
Cesium 常用地图纠偏 百度、高德、Arcgis、腾讯地图纠偏
04-10
Cesium 常用地图纠偏 百度、高德、Arcgis、腾讯地图纠偏 资源预览:http://dongnan185.com:8083/videos/jiupian.mp4 启动:npm i npm run serve
天地图-坐标拾取
最新发布
qq_44681558的博客
03-26 971
在网页使用高德，百度地图这些时，会因为没有购买他们的使用版权，从而出现以下水印，因此需要换一个地图-天地图 天地图官方文档：http://lbs.tianditu.gov.cn/ 然而天地图并不像其他地图一样这么方便，比如没有可以通过地址搜索，直接可以获取坐标的，因此可以小改动一下，方便自己使用。
前端地图坐标系详解
Aurora_1970s的博客
08-10 2171
1.地图类型 国外 国外一般就是GPS，谷歌地图 国内 百度地图，高德地图，腾讯地图，搜狗地图，天地图… 2.坐标系类型 WGS84坐标系 最常见的GPS，用的是WGS-84坐标系，也是原始坐标体系。一般用国际标准的GPS记录仪记录下来的坐标，都是GPS的坐标。 GCJ-02坐标系 国测局02年发布的坐标体系。又称火星坐标系。火星坐标系就是在84坐标系的基础上进行了一次加密偏移所得到的坐标系。在中国，必须至少使用GCJ-02的坐标体系。比如腾讯，高德都在用这个坐标体系。GCJ-02也是国内最广泛使用的坐标体
【GIS坐标系介绍】WGS84坐标系、Web墨卡托、GCJ02坐标系、BD09坐标系
qq_16599281的博客
10-28 4834
【GIS坐标系介绍】WGS84坐标系、Web墨卡托、GCJ02坐标系、BD09坐标系
GIS-国内常见的地理坐标系
sunny05296的博客
06-14 5447
GIS-国内常见的地理坐标系WGS84(EPSG:4326) : 目前最流行的地理坐标系统，美国GPS就是使用的这个。Pseudo-Mercator(EPSG:3857) : 投影坐标系。伪墨卡托投影，也称为球体墨卡托，Web Mercator。它是基于墨卡托投影的，把 WGS84 投影到正方形。各大互联网地图公司以它为准。伪墨卡托非常适合显示数据、但不适合存储数据。一般用 WGS84 存储数据、用伪墨卡托显示数据。CGCS2000(4490) : 国家大地坐标系。天地图采用的就是这个。可以提供高精度三维坐
天地图坐标系
weixin_45188172的博客
09-13 3353
WGS84: WGS84坐标系，全称为世界大地测量系统（英语：World Geodetic System）是一种用于地图学、大地测量学和导航（包括全球定位系统）的大地测量系统标准。WGS84坐标是观测历元的动态坐标，用于导航，强调实时性、动态性。两者用途不同，特点不同，但都统一于ITRS坐标系，都对准ITRF框架。可通过历元归算、框架转换互相转换。
国内常用地图坐标系总结
热门推荐
prcyang的专栏
05-12 1万+
国内地图坐标系
《学习总结》天地图
Leporc的博客
08-05 3914
国内比较常用的地图api开发平台大概就百度、高德、谷歌、天地图了其中百度、高德、谷歌都要收费，天地图完全免费天地图用的是CGCS-2000坐标系，与WGS84坐标系都是地心坐标系，相差不大。高德采用的是GCJ-02坐标系，也就是再WGS84坐标系上进行了一次加密得到的坐标系百度采用的是BD-09坐标系，是在GCJ-02上再进行了一次加密得到的。.........
天地图、百度地图、高德地图经纬度＜=＞墨卡托坐标转换
weixin_41234005的博客
11-22 2645
首先了解一下火星坐标，它是在国际标准坐标WGS-84上进行的一次加密，由于国内的电子地图都要至少使用火星坐标进行一次加密，百度直接就任性一些，直接自己又研究了一套加密算法，来了个二次加密，这就是我们所熟知的百度坐标(BD-09)。火星坐标是国家测绘局为了国家安全在原始坐标的基础上进行偏移得到的坐标，基本国内的电子地图、导航设备都是采用的这一坐标系或在这一坐标的基础上进行二次加密得到的。（4）一些注意点：非专业人士，对于墨卡托投影坐标和经纬度坐标可以理解为一个是平面坐标一个是球面坐标。
cj_wg_地图_高德_谷歌_坐标系_百度_
10-03
谷歌、高德、百度地图坐标系互转JS源码，对电子地图开发的很有帮助。
高德天地图坐标转换.xlsm
07-03
excel中高德与天地图坐标的转换，可使用内部定义的函数实现批量转换
百度服务获取坐标放置在天地图上实例
rrrrroy_Ha的博客
04-18 5340
开始之前还是要提前写一下涉及到的几个坐标系 WGS-84：是国际标准，GPS坐标（GoogleEarth使用、或者GPS模块） GCJ-02：中国坐标偏移标准，GoogleMap、高德、腾讯使用 BD-09：百度坐标偏移标准，BaiduMap使用 CGCS2000：天地图使用 先以北京华宾公司为例，地址是“北京市西城区西安门大街107号”，百度地图上是如下位置 访问百度WEB服务...
GIS数据漫谈（五）— 地理坐标系统
ThingJS的博客
07-27 2418
2000国家大地坐标系。由于天地图的坐标系统是CGCS2000的，在天地图上拾取一个点的经纬度来观察。比如由于区域的特殊性，公开的矢量数据本身就在原始测绘数据的基础上进行了偏移，降低精度。这里用的是QGIS，通过XYZTiles图层加载的天地图卫星底图，再把shp文件叠加加载。（1）一般情况下是由于经纬度的坐标系和地图底图的坐标系不匹配导致。（腾讯地图同理），又称火星坐标系，它是一种对经纬度数据的加密算法，即加入随机的偏差。的标准三维场景的构建服务，使城市三维场景的构建速度提升至分钟级。...
openlayer 6 + 天地图之坐标系
josiecici的博客
09-06 1388
2. 什么是EPSG:3857坐标系（投影坐标） 在openlayers 中默认的坐标就是google的摩卡托坐标，也就是我们经常看到的 EPSG:3857 坐标系。 EPSG:3857 的数据一般是这种的。[12914838.35,4814529.9]，看上去相对数值较大。不利于存储，比较占内存。 3. 什么是EPSG:4326 坐标系（地理坐标） 4326 WGS-84：是国际标准，GPS坐标（Google Earth使用、或者GPS模块） EPSG:4326 的数据一般是这种的。[22.37，114.
GIS开发：天地图开放地图投影
地理信息技术杂谈 博客，欢迎技术交流
05-20 811
引言 天地图（国家地理信息公共服务平台）开放地图服务地址： https://service.tianditu.gov.cn/#/ 有全球、全国、各省的开放地图服务数据，一般提供了wmts的切片服务，可以在客户端的地图上进行调用。 地图的投影有两种：一种是经纬度，另一种是球面墨卡托。 经纬度 天地图发布的经纬度是cgcs2000坐标系的投影，wmts的切片规则可参考下列网址进行查看： http://t0.tianditu.gov.cn/vec_c/wmts?service=wmts&request=
地图坐标系以及常见的坐标转换
博客简介
05-29 1627
前端在对接硬件厂商发来的点位时需注意的事项：选择地图以及坐标系等
js百度地图坐标转高德地图坐标
11-30
要将JS百度地图坐标转换为高德地图坐标，可以使用第三方的坐标转换服务来实现。首先，需要获取到百度地图的经纬度坐标，然后通过调用第三方的坐标转换接口，将百度地图坐标转换为高德地图坐标。接着，将转换后的高德地图坐标应用到相应的地图服务上，以显示在高德地图上。 在JS中，可以使用ajax或fetch等方式来调用第三方的坐标转换接口，将百度地图坐标作为参数传入，然后获取到转换后的高德地图坐标。然后，将获得的高德地图坐标应用到相应的地图服务上，通过相应的方法将标记点或者路径在高德地图上显示出来。 需要注意的是，使用第三方的坐标转换服务可能需要申请相应的API密钥或者进行注册，并且可能需要支付一定的费用。另外，转换坐标时也需要考虑坐标系的不同以及可能存在的误差，尽量选择精准的坐标转换服务来进行转换。 总之，通过调用第三方的坐标转换服务，以及在JS中处理转换后的坐标数据，即可实现将百度地图坐标转换为高德地图坐标，并在高德地图上进行显示。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 12
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值