地图相关坐标转换及说明

关于坐标系说明

简单介绍一下常见的地理经纬度坐标系：

WGS-84 

WGS-84坐标系 是国际标准坐标系，GPS坐标（Google Earth使用、或者GPS模块）

GCJ-02  

GCJ-02坐标系，即火星坐标系，是由中国国家测绘局制订的地理信息系统的坐标系统。

它是一种对经纬度数据的加密算法，即加入随机的偏差。

国内出版的地图系统包括腾讯地图和高德地图都采用了GCJ-02火星坐标系。

出于国家安全考虑，必须至少采用GCJ-02对地理位置进行首次加密。

BD-09

BD-09，即百度坐标系，也就是在GCJ-02的基础上，进行了BD-09二次加密措施生成的一种坐标系。

关于坐标转换及理解推荐一篇文章：

国内地图坐标转换 - 掘金

地图坐标转换API

 /*
   *   计算经纬度之间的距离(公里)
   */
  this.GetDistance = (lng1, lat1, lng2, lat2) => {
    try {
      var radLat1 = lat1 * Math.PI / 180.0;
      var radLat2 = lat2 * Math.PI / 180.0;
      var a = radLat1 - radLat2;
      var b = (lng1 * Math.PI / 180.0) - (lng2 * Math.PI / 180.0);
      var s = 2 * Math.asin(Math.sqrt(Math.pow(Math.sin(a / 2), 2) + Math.cos(radLat1) * Math.cos(radLat2) * Math.pow(
        Math.sin(b / 2), 2)));
      s = s * 6378.137;
      s = Math.round(s * 10000) / 10000;
      radLat1 = null;
      radLat2 = null;
      a = null;
      b = null;
      return s;
    } catch (ex) {
      console.log(ex);
    }
  };
 
/**
     * WGS84转GCj02
     * @param lng
     * @param lat
     * @returns {*[]}
     */
    const x_PI = 3.14159265358979324 * 3000.0 / 180.0;
    const PI = 3.1415926535897932384626;
    const a = 6378245.0;
    const ee = 0.00669342162296594323;
    this.wgs84togcj02 = (lng, lat) => {
      if (this.out_of_china(lng, lat)) {
        return [lng, lat]
      } else {
        let dlat = this.transformlat(lng - 105.0, lat - 35.0);
        let dlng = this.transformlng(lng - 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);
        dlng = (dlng * 180.0) / (a / sqrtmagic * Math.cos(radlat) * PI);
        let mglat = lat + dlat;
        let mglng = lng + dlng;
        return [mglng, mglat]
      }
    }
    
    /**
     * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换
     * 即谷歌、高德 转 百度
     * @param lng
     * @param lat
     * @returns {*[]}
     */
    this.gcj02tobd09 = (lng, lat) => {
      let z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * x_PI);
      let theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * x_PI);
      let bd_lng = z * Math.cos(theta) + 0.0065;
      let bd_lat = z * Math.sin(theta) + 0.006;
      return [bd_lng, bd_lat]
    }    
 
    /**
     * 百度坐标系gcj02的转换
     */   
    this.bd09togcj02 = (bd_lon, bd_lat) => {
      bd_lon = +bd_lon;
      bd_lat = +bd_lat;
      let x = bd_lon - 0.0065;
      let 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);
      let gg_lng = z * Math.cos(theta);
      let gg_lat = z * Math.sin(theta);
      return [gg_lng, gg_lat]
    }
 
     /**
     * 将gcj02转换为wgs84
     * @param lng
     * @param lat
     * @returns []
     */
    this.gcj02towgs84 = (lng, lat) => {
      var lat = +lat;
      var lng = +lng;
        if (this.out_of_china(lng, lat)) {
            return [lng, lat]
        } else {
            var dlat = this.transformlat(lng - 105.0, lat - 35.0);
            var dlng = this.transformlng(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 [lng * 2 - mglng, lat * 2 - mglat]
      }
    }
 
     /**
     * 将百度坐标转换为wgs84
     */
    this.bd09towgs84 = (lng, lat) =>{
      let gcj02Points = this.bd09togcj02(lng,lat);
      let wpg84Points = this.gcj02towgs84(gcj02Points[0],gcj02Points[1]);
      return wpg84Points;
    }
 
     /**
     * 将wgs84转换为百度坐标 
     */
    this.wgs84tobd09 = (lng, lat) =>{
      let gcj02PointT = this.wgs84togcj02(lng,lat);
      let bdPoints = this.gcj02tobd09(gcj02PointT[0],gcj02PointT[1]);
      return bdPoints;
    }
 
    /**
     * 判断是否在国内，不在国内则不做偏移
     * @param lng
     * @param lat
     * @returns {boolean}
     */
    this.out_of_china = (lng, lat) => {
      return (lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false);
    }
    this.transformlat = (lng, lat) => {
      let 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
    }
    this.transformlng = (lng, lat) => {
      let 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
    }
 
    //百度坐标转高德（传入经度、纬度）
    this.bd_decrypt = function (bd_lng, bd_lat) {
      var X_PI = Math.PI * 3000.0 / 180.0;
      var x = bd_lng - 0.0065;
      var y = bd_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);
      var gg_lng = z * Math.cos(theta);
      var gg_lat = z * Math.sin(theta);
      return [gg_lng,gg_lat]
    }
 
    //google坐标转高德（传入经度、纬度）
    this.bd_encrypt = function (gg_lng, gg_lat) {
      var X_PI = Math.PI * 3000.0 / 180.0;
      var x = gg_lng,
        y = gg_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);
      var bd_lng = z * Math.cos(theta) + 0.0065;
      var bd_lat = z * Math.sin(theta) + 0.006;
      return [ bd_lat,bd_lng]
    }