地图相关坐标转换及说明

已于 2025-02-26 09:59:39 修改
阅读量436 收藏
点赞数
文章标签： 地图坐标转换及说明 javascript
于 2023-10-11 11:59:27 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/lwx931449660/article/details/133764810
版权
本文详细介绍了地理坐标系如WGS-84、GCJ-02（中国的火星坐标系）和BD-09（百度坐标系）的概念，以及它们之间的转换方法，包括WGS-84到GCJ-02、GCJ-02到BD-09等，并提供了JavaScript代码示例。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

关于坐标系说明

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

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

/**
 * 计算两个经纬度点之间的距离
 * @param {number} lon1 第一个点的经度
 * @param {number} lat1 第一个点的纬度
 * @param {number} lon2 第二个点的经度
 * @param {number} lat2 第二个点的纬度
 * @returns {number} 两点之间的距离（单位：米）
 */
export function 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 * 1000
  } catch (ex) {
    console.log(ex)
    return 0
  }
}



/**
     * 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]
    }

验证经纬度数据是否符合格式

// 经度正则表达式，整数部分为0-180,小数部分为0到15位
const lngReg =
  /^(\-|\+)?(((\d|[1-9]\d|1[0-7]\d|0{1,3})\.\d{0,15})|(\d|[1-9]\d|1[0-7]\d|0{1,3})|180\.0{0,15}|180)$/

// 纬度正则表达式，整数部分为0-90小数部分为0到15位
const latReg = /^(\-|\+)?([0-8]?\d{1}\.\d{0,15}|90\.0{0,15}|[0-8]?\d{1}|90)$/

/**
 * 检查给定的经纬度是否是有效的坐标
 * @param {number} longitude 经度
 * @param {number} latitude 纬度
 * @returns {boolean} 是否是有效的坐标
 */
export function isValidCoordinates(longitude, latitude) {
  return lngReg.test(longitude) && latReg.test(latitude)
}

地球坐标转高德坐标 


const PI = 3.1415926535897932384626;
const a = 6378245.0;
const ee = 0.00669342162296594323;
 
function 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;
}
function 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;
}
 
function outOfChina(lng, lat) {
  return lng < 72.004 || lng > 137.8347 || lat < 0.8293 || lat > 55.8271 || false;
}
 
export default function (lng, lat) {
  if (outOfChina(lng, lat)) {
    return [lng, lat];
  } else {
    let dlat = transformLat(lng - 105.0, lat - 35.0);
    let dlng = transformLng(lng - 105.0, lat - 35.0);
    const radlat = (lat / 180.0) * PI;
    let magic = Math.sin(radlat);
    magic = 1 - ee * magic * magic;
    const sqrtmagic = Math.sqrt(magic);
    dlat = (dlat * 180.0) / (((a * (1 - ee)) / (magic * sqrtmagic)) * PI);
    dlng = (dlng * 180.0) / ((a / sqrtmagic) * Math.cos(radlat) * PI);
    const mglat = lat + dlat;
    const mglng = lng + dlng;
    return [+mglng.toFixed(6), +mglat.toFixed(6)];
  }
}
 
// 使用
import wgs84ToGCJ02 from '/@/utils/map/wgs84ToGCJ02';
const [longitude, latitude] = wgs84ToGCJ02(+item.longitude, +item.latitude);

丿刘先森
关注
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值