常见地图坐标系以及转换方法、转换工具

l 为什么要进行坐标系转换？

因为国内对地理坐标有特殊的政策，GPS获取的原始坐标不能在互联网上对外提供，所有电子地图必须对位置做偏移处理，这就导致从底层接口得到的坐标展示在地图上会有偏移。这种偏移不是线性的，不能简单地通过加减来校正，并且不同的地图提供商采用的算法也不一样。不同的坐标系之间可能有几十到几百米的偏移，所以在地图上绘制矢量图形、做地理数据可视化时，需要修正不同坐标系之间的偏差。

l 关于坐标系

通常GPS获取的是WGS84坐标体系，包括很多国外地图厂商提供的地图数据也是WGS84坐标体系，但是在我国，由于政策原因，在地图发布和出版的时候，坐标至少需要经过国家测绘局加密，也就是对WGS84坐标进行一次非线性加偏，得到GCJ02坐标系，俗称火星坐标。

另外，国内许多地图公司，由于各种原因，会在火星坐标的基础上进行二次非线性加偏，得到自己的坐标系统，各地图厂商使用坐标系情况如下。

坐标系	描述	地图厂商
WGS-84	地心坐标系，GPS原始坐标体系 是世界通用坐标系，也是目前广泛使用的GPS全球卫星定位系统使用的坐标系	MapBox、Cerium、OpenStreetMap
CGCS2000	我国当前最新的国家大地坐标系 原点为包括海洋和大气的整个地球的质量中心，和WGS84椭球极为相似，偏差仅有0.11mm，完全可以兼容使用	天地图
GCJ-02	又称火星坐标系 由WGS84坐标系经加密后形成，是由中国国家测绘局制订的地理信息系统的坐标系统	高德地图、腾讯地图、谷歌中国
BD-09	百度坐标系 在GCJ02坐标系基础上再次加密，从百度产品中得到的坐标都是BD-09坐标系	百度地图

l 坐标转换方法

整理了一下坐标系转换的方法，仅限单个数组转换

火星坐标转wgs84

// 火星坐标转wgs84
function gcj02towgs84(lng, lat) {
var x_PI = 3.14159265358979324 * 3000.0 / 180.0;
var PI = 3.1415926535897932384626;
var a = 6378245.0;
var ee = 0.00669342162296594323;
if ((lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false)) {
return [lng, lat]
}
else {
var dlat = ((lng, lat) => {
var ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0.1 * lng * lat + 0.2 * Math.sqrt(Math.abs(lng));
ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;
ret += (20.0 * Math.sin(lat * PI) + 40.0 * Math.sin(lat / 3.0 * PI)) * 2.0 / 3.0;
ret += (160.0 * Math.sin(lat / 12.0 * PI) + 320 * Math.sin(lat * PI / 30.0)) * 2.0 / 3.0;
return ret
})(lng - 105.0, lat - 35.0);
var dlng = ((lng, lat) => {
var ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math.sqrt(Math.abs(lng));
ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;
ret += (20.0 * Math.sin(lng * PI) + 40.0 * Math.sin(lng / 3.0 * PI)) * 2.0 / 3.0;
ret += (150.0 * Math.sin(lng / 12.0 * PI) + 300.0 * Math.sin(lng / 30.0 * PI)) * 2.0 / 3.0;
return ret
})(lng - 105.0, lat - 35.0);
var radlat = lat / 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);
dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);
mglat = lat + dlat;
mglng = lng + dlng;
return [lng * 2 - mglng, lat * 2 - mglat]
}
}

wgs84转火星坐标

// wgs84转火星坐标
function wgs84togcj02(lng, lat) {
var x_PI = 3.14159265358979324 * 3000.0 / 180.0;
var PI = 3.1415926535897932384626;
var a = 6378245.0;
var ee = 0.00669342162296594323;
if ((lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false)) {
return [lng, lat]
} else {
var dlat = ((lng, lat) => {
var ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * lat + 0.1 * lng * lat + 0.2 * Math.sqrt(Math.abs(lng));
ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;
ret += (20.0 * Math.sin(lat * PI) + 40.0 * Math.sin(lat / 3.0 * PI)) * 2.0 / 3.0;
ret += (160.0 * Math.sin(lat / 12.0 * PI) + 320 * Math.sin(lat * PI / 30.0)) * 2.0 / 3.0;
return ret
})(lng - 105.0, lat - 35.0);
var dlng = ((lng, lat) => {
var ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + 0.1 * lng * lat + 0.1 * Math.sqrt(Math.abs(lng));
ret += (20.0 * Math.sin(6.0 * lng * PI) + 20.0 * Math.sin(2.0 * lng * PI)) * 2.0 / 3.0;
ret += (20.0 * Math.sin(lng * PI) + 40.0 * Math.sin(lng / 3.0 * PI)) * 2.0 / 3.0;
ret += (150.0 * Math.sin(lng / 12.0 * PI) + 300.0 * Math.sin(lng / 30.0 * PI)) * 2.0 / 3.0;
return ret
}
)(lng - 105.0, lat - 35.0);
var radlat = lat / 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);
dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);
var mglat = lat + dlat;
var mglng = lng + dlng;
return [mglng, mglat]
}
}

百度转火星坐标

// 百度转火星坐标
const bd2xh = (bd_lng, bd_lat) => {
const x_pi = Math.PI * 3000.0 / 180.0;
const x = bd_lng - 0.0065;
const y = bd_lat - 0.006;
let z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * x_pi);
let theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_pi);
return [z * Math.cos(theta), z * Math.sin(theta)];
}

l 坐标转换工具

有的静态数据不需要动态转换也有现成的在线坐标系转换工具可以直接拿来用

EasyMaphttps://map.easyv.cloud/transform

支持常见百度地图（BD09）、火星坐标系（GCJ02）、WGS84坐标系之间的相互转换

具体支持的坐标系转换如下

1）WGS-84 → GCJ-02

2）WGS-84 → BD-09

2）GCJ-02 → WGS-84

3）GCJ-02 → BD-09

4）BD-09 → GCJ-02

5）BD-09 → WGS-84

支持单个数组、不包含嵌套关系的对象、标准GeoJSON

这个工具也还可以获取行政边界数据、地图数据模拟