AI 助理
备案控制台
开发者社区 开发与运维 文章 正文

84坐标系、02坐标系、百度坐标之间相互转换算法

2023-06-11 1066
版权
版权声明：
本文内容由阿里云实名注册用户自发贡献，版权归原作者所有，阿里云开发者社区不拥有其著作权，亦不承担相应法律责任。具体规则请查看《 阿里云开发者社区用户服务协议》和 《阿里云开发者社区知识产权保护指引》。如果您发现本社区中有涉嫌抄袭的内容，填写 侵权投诉表单进行举报，一经查实，本社区将立刻删除涉嫌侵权内容。
简介： 最近有同学反馈之前的坐标系转换有问题，对之前的工具类进行了修正。 一、地图坐标转换java工具类 包含84坐标系、02坐标系、百度地图、高德地图、腾讯地图坐标之间相互转换的算法wgs84ToGcj02：将 WGS84 坐标系下的经纬度转换为 GCJ02 坐标系下的经纬度。gcj02ToWgs84：将 GCJ02 坐标系下的经纬度转换为 WGS84 坐标系下的经纬度。gcj02ToBd09：将 GCJ02 坐标系下的经纬度转换为 BD09 坐标系下的经纬度。bd09ToGcj02：将 BD09 坐标系下的经纬度转换为 GCJ02 坐标系下的经纬度。

一、地图坐标转换java工具类        

包含84坐标系、02坐标系、百度地图、高德地图、腾讯地图坐标之间相互转换的算法

wgs84ToGcj02：将 WGS84 坐标系下的经纬度转换为 GCJ02 坐标系下的经纬度。

gcj02ToWgs84：将 GCJ02 坐标系下的经纬度转换为 WGS84 坐标系下的经纬度。

gcj02ToBd09：将 GCJ02 坐标系下的经纬度转换为 BD09 坐标系下的经纬度。

bd09ToGcj02：将 BD09 坐标系下的经纬度转换为 GCJ02 坐标系下的经纬度。

gcj02ToQQ：将 GCJ02 坐标系下的经纬度转换为腾讯坐标系（QQ）下的经纬度。

qqToGcj02：将腾讯坐标系（QQ）下的经纬度转换为 GCJ02 坐标系下的经纬度。

import com.alibaba.fastjson.JSONObject;
public class CoordinateConversionUtils {
    private static final double EARTH_RADIUS = 6378137; // 地球半径，单位：米
    private static final double X_PI = Math.PI * 3000.0 / 180.0;
    private static final double A = 6378245.0;
    private static final double EE = 0.00669342162296594323;
    /**
     * 将84坐标系（WGS84）转换为02坐标系（GCJ02）
     *
     * @param lng 经度
     * @param lat 纬度
     * @return 转换后的经纬度
     */
    public static JSONObject wgs84ToGcj02(double lng, double lat) {
        JSONObject result = new JSONObject();
        if (outOfChina(lng, lat)) {
            result.put("lng", lng);
            result.put("lat", lat);
            return result;
        }
        double dLat = transformLat(lng - 105.0, lat - 35.0);
        double dLng = transformLng(lng - 105.0, lat - 35.0);
        double radLat = lat / 180.0 * Math.PI;
        double magic = Math.sin(radLat);
        magic = 1 - EE * magic * magic;
        double sqrtMagic = Math.sqrt(magic);
        dLat = (dLat * 180.0) / ((EARTH_RADIUS * (1 - EE)) / (magic * sqrtMagic) * Math.PI);
        dLng = (dLng * 180.0) / (EARTH_RADIUS / sqrtMagic * Math.cos(radLat) * Math.PI);
        double mgLat = lat + dLat;
        double mgLng = lng + dLng;
        result.put("lng", mgLng);
        result.put("lat", mgLat);
        return result;
    }
    /**
     * 将02坐标系（GCJ02）转换为84坐标系（WGS84）
     *
     * @param lng 经度
     * @param lat 纬度
     * @return 转换后的经纬度
     */
    public static JSONObject gcj02ToWgs84(double lng, double lat) {
        JSONObject result = new JSONObject();
        if (outOfChina(lng, lat)) {
            result.put("lng", lng);
            result.put("lat", lat);
            return result;
        }
        double dLat = transformLat(lng - 105.0, lat - 35.0);
        double dLng = transformLng(lng - 105.0, lat - 35.0);
        double radLat = lat / 180.0 * Math.PI;
        double magic = Math.sin(radLat);
        magic = 1 - EE * magic * magic;
        double sqrtMagic = Math.sqrt(magic);
        dLat = (dLat * 180.0) / ((EARTH_RADIUS * (1 - EE)) / (magic * sqrtMagic) * Math.PI);
        dLng = (dLng * 180.0) / (EARTH_RADIUS / sqrtMagic * Math.cos(radLat) * Math.PI);
        double mgLat = lat + dLat;
        double mgLng = lng + dLng;
        result.put("lng", lng * 2 - mgLng);
        result.put("lat", lat * 2 - mgLat);
        return result;
    }
    /**
     * 将02坐标系（GCJ02）转换为百度坐标系（BD09）
     *
     * @param lng 经度
     * @param lat 纬度
     * @return 转换后的经纬度
     */
    public static JSONObject gcj02ToBd09(double lng, double lat) {
        JSONObject result = new JSONObject();
        double z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * X_PI);
        double theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * X_PI);
        double bdLng = z * Math.cos(theta) + 0.0065;
        double bdLat = z * Math.sin(theta) + 0.006;
        result.put("lng", bdLng);
        result.put("lat", bdLat);
        return result;
    }
    /**
     * 将百度坐标系（BD09）转换为02坐标系（GCJ02）
     *
     * @param lng 经度
     * @param lat 纬度
     * @return 转换后的经纬度
     */
    public static JSONObject bd09ToGcj02(double lng, double lat) {
        JSONObject result = new JSONObject();
        double x = lng - 0.0065, y = lat - 0.006;
        double z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * X_PI);
        double theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * X_PI);
        double gcjLng = z * Math.cos(theta);
        double gcjLat = z * Math.sin(theta);
        result.put("lng", gcjLng);
        result.put("lat", gcjLat);
        return result;
    }
    /**
     * 将02坐标系（GCJ02）转换为腾讯坐标系（QQ）
     *
     * @param lng 经度
     * @param lat 纬度
     * @return 转换后的经纬度
     */
    public static JSONObject gcj02ToQQ(double lng, double lat) {
        JSONObject result = new JSONObject();
        double x = lng;
        double y = lat;
        double z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * X_PI);
        double theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * X_PI);
        double qqLng = z * Math.cos(theta) + 0.00330;
        double qqLat = z * Math.sin(theta) + 0.00150;
        result.put("lng", qqLng);
        result.put("lat", qqLat);
        return result;
    }
    /**
     * 将腾讯坐标系（QQ）转换为02坐标系（GCJ02）
     *
     * @param lng 经度
     * @param lat 纬度
     * @return 转换后的经纬度
     */
    public static JSONObject qqToGcj02(double lng, double lat) {
        JSONObject result = new JSONObject();
        double x = lng - 0.00330;
        double y = lat - 0.00150;
        double z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * X_PI);
        double theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * X_PI);
        double gcjLng = z * Math.cos(theta);
        double gcjLat = z * Math.sin(theta);
        result.put("lng", gcjLng);
        result.put("lat", gcjLat);
        return result;
    }
    /**
     * 判断是否超出中国范围
     *
     * @param lng 经度
     * @param lat 纬度
     * @return 是否超出中国范围
     */
    private static boolean outOfChina(double lng, double lat) {
        return lng < 72.004 || lng > 137.8347 || lat < 0.8293 || lat > 55.8271;
    }
    /**
     * 转换经度
     *
     * @param lng 经度
     * @param lat 纬度
     * @return 转换后的经度
     */
    private static double transformLng(double lng, double lat) {
        double 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 * Math.PI) + 20.0 * Math.sin(2.0 * lng * Math.PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(lng * Math.PI) + 40.0 * Math.sin(lng / 3.0 * Math.PI)) * 2.0 / 3.0;
        ret += (150.0 * Math.sin(lng / 12.0 * Math.PI) + 300.0 * Math.sin(lng / 30.0 * Math.PI)) * 2.0 / 3.0;
        return ret;
    }
    /**
     * 转换纬度
     *
     * @param
* lat 纬度
     * @return 转换后的纬度
     */
    private static double transformLat(double lng, double lat) {
            double 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 * Math.PI) + 20.0 * Math.sin(2.0 * lng * Math.PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(lat * Math.PI) + 40.0 * Math.sin(lat / 3.0 * Math.PI)) * 2.0 / 3.0;
        ret += (160.0 * Math.sin(lat / 12.0 * Math.PI) + 320 * Math.sin(lat / 30.0 * Math.PI)) * 2.0 / 3.0;
        return ret;
    }
}

wgs84ToGcj02` 和 `gcj02ToWgs84` 方法采用了经典的坐标转换算法，即将 WGS84 坐标系下的经纬度转换为火星坐标系（GCJ02 坐标系）下的经纬度，或将火星坐标系下的经纬度转换为 WGS84 坐标系下的经纬度。这两个方法的转换精度较高，适用于绝大部分场景。

`gcj02ToBd09` 和 `bd09ToGcj02` 方法分别实现了 GCJ02 坐标系和 BD09 坐标系之间的转换，这两种坐标系都是由中国的地图服务提供商所定义的，因此在国内的地图服务中比较常见。

`gcj02ToQQ` 和 `qqToGcj02` 方法分别实现了 GCJ02 坐标系和腾讯坐标系（QQ）之间的转换，QQ 坐标系是腾讯地图特有的一种坐标系，因此在使用腾讯地图时比较有用。

以上方法可以满足大部分常见的坐标转换需求，但在一些特殊情况下可能会存在一定的误差。如需更高精度的坐标转换算法，可以参考一些开源的第三方库，如 Proj4J、Geotools 等。

二、JavaScript 实现

       84坐标系、02坐标系、百度坐标之间相互转换算法

/**
 * 提供了百度坐标（BD09）、国测局坐标（火星坐标，GCJ02）、和WGS84坐标系之间的转换
 */
//UMD魔法代码
// if the module has no dependencies, the above pattern can be simplified to
(function (root, factory) {
  if (typeof define === 'function' && define.amd) {
    // AMD. Register as an anonymous module.
    define([], factory);
  } else if (typeof module === 'object' && module.exports) {
    // Node. Does not work with strict CommonJS, but
    // only CommonJS-like environments that support module.exports,
    // like Node.
    module.exports = factory();
  } else {
    // Browser globals (root is window)
    root.coordtransform = factory();
  }
}(this, function () {
  //定义一些常量
  var x_PI = 3.14159265358979324 * 3000.0 / 180.0;
  var PI = 3.1415926535897932384626;
  var a = 6378245.0;
  var ee = 0.00669342162296594323;
  /**
   * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
   * 即 百度 转 谷歌、高德
   * @param bd_lon
   * @param bd_lat
   * @returns {*[]}
   */
  var bd09togcj02 = function bd09togcj02(bd_lon, bd_lat) {
    var bd_lon = +bd_lon;
    var bd_lat = +bd_lat;
    var x = bd_lon - 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]
  };
  /**
   * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换
   * 即谷歌、高德 转 百度
   * @param lng
   * @param lat
   * @returns {*[]}
   */
  var gcj02tobd09 = function gcj02tobd09(lng, lat) {
    var lat = +lat;
    var lng = +lng;
    var z = Math.sqrt(lng * lng + lat * lat) + 0.00002 * Math.sin(lat * x_PI);
    var theta = Math.atan2(lat, lng) + 0.000003 * Math.cos(lng * x_PI);
    var bd_lng = z * Math.cos(theta) + 0.0065;
    var bd_lat = z * Math.sin(theta) + 0.006;
    return [bd_lng, bd_lat]
  };
  /**
   * WGS84转GCj02
   * @param lng
   * @param lat
   * @returns {*[]}
   */
  var wgs84togcj02 = function wgs84togcj02(lng, lat) {
    var lat = +lat;
    var lng = +lng;
    if (out_of_china(lng, lat)) {
      return [lng, lat]
    } else {
      var dlat = transformlat(lng - 105.0, lat - 35.0);
      var dlng = 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 [mglng, mglat]
    }
  };
  /**
   * GCJ02 转换为 WGS84
   * @param lng
   * @param lat
   * @returns {*[]}
   */
  var gcj02towgs84 = function gcj02towgs84(lng, lat) {
    var lat = +lat;
    var lng = +lng;
    if (out_of_china(lng, lat)) {
      return [lng, lat]
    } else {
      var dlat = transformlat(lng - 105.0, lat - 35.0);
      var dlng = 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]
    }
  };
  var transformlat = function transformlat(lng, lat) {
    var lat = +lat;
    var lng = +lng;
    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
  };
  var transformlng = function transformlng(lng, lat) {
    var lat = +lat;
    var lng = +lng;
    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
  };
  /**
   * 判断是否在国内，不在国内则不做偏移
   * @param lng
   * @param lat
   * @returns {boolean}
   */
  var out_of_china = function out_of_china(lng, lat) {
    var lat = +lat;
    var lng = +lng;
    // 纬度3.86~53.55,经度73.66~135.05 
    return !(lng > 73.66 && lng < 135.05 && lat > 3.86 && lat < 53.55);
  };
  return {
    bd09togcj02: bd09togcj02,
    gcj02tobd09: gcj02tobd09,
    wgs84togcj02: wgs84togcj02,
    gcj02towgs84: gcj02towgs84
  }
}));

8a321a19e9aeb7463baf924e4e2bf97e.jpg


文章标签：
Java
JavaScript
前端开发
定位技术
算法
关键词：
坐标算法
百度坐标
坐标系百度
算法坐标
冰点.
目录
相关文章
秋夜雨巷
|
3月前
|
前端开发 JavaScript 定位技术
GPS坐标转百度坐标
GPS坐标转百度坐标
秋夜雨巷
34 1
zcq5jk6y5y5e4
|
3月前
|
算法 前端开发 计算机视觉
基于均值坐标(Mean-Value Coordinates)的图像融合算法的优化实现
基于均值坐标(Mean-Value Coordinates)的图像融合算法的优化实现
zcq5jk6y5y5e4
39 0
zcq5jk6y5y5e4
|
3月前
|
自然语言处理 并行计算 算法
基于均值坐标(Mean-Value Coordinates)的图像融合算法的具体实现
基于均值坐标(Mean-Value Coordinates)的图像融合算法的具体实现
zcq5jk6y5y5e4
43 0
朝阳39
|
4月前
|
JavaScript 定位技术
vue-baidu-map 百度地图检索、获取坐标
vue-baidu-map 百度地图检索、获取坐标
朝阳39
46 1
jjczqyjxgdegi
|
5月前
|
算法 Java 定位技术
心得经验总结：标准坐标系与火星坐标系（高德）百度坐标系之间互转
心得经验总结：标准坐标系与火星坐标系（高德）百度坐标系之间互转
jjczqyjxgdegi
65 1
漏刻有时
|
定位技术 开发工具 数据安全/隐私保护
GIS数据格式坐标转换（地球坐标WGS84、GCJ-02、火星坐标、百度坐标BD-09、国家大地坐标系CGCS2000）
GIS数据格式坐标转换（地球坐标WGS84、GCJ-02、火星坐标、百度坐标BD-09、国家大地坐标系CGCS2000）
漏刻有时
2103 1
fensnote
|
6月前
|
JavaScript 前端开发 定位技术
GPS坐标显示在百度地图上（Qt+百度地图）
GPS坐标显示在百度地图上（Qt+百度地图）
fensnote
79 0
猿长大人
|
11月前
|
算法 数据可视化 C#
C# | Chaikin算法 —— 计算折线对应的平滑曲线坐标点
本文将介绍一种计算折线对应的平滑曲线坐标点的算法。该算法使用Chaikin曲线平滑处理的方法，通过控制张力因子和迭代次数来调整曲线的平滑程度和精度。通过对原始点集合进行切割和插值操作，得到平滑的曲线坐标点集合。实验结果表明，该算法能够有效地平滑折线，并且具有较高的精度和可控性。
猿长大人
225 0
C# | Chaikin算法 —— 计算折线对应的平滑曲线坐标点
热烈的马
|
6月前
|
机器学习/深度学习 算法 数据挖掘
【python机器学习】K-Means算法详解及给坐标点聚类实战（附源码和数据集 超详细）
【python机器学习】K-Means算法详解及给坐标点聚类实战（附源码和数据集 超详细）
热烈的马
287 0
漏刻有时
|
定位技术 Python
百度地图webAPI：python坐标批量转换
百度地图webAPI：python坐标批量转换
漏刻有时
110 0

热门文章

最新文章

  • 1
    百度余凯：数据是极好的竞争壁垒
  • 2
    我自己开发的工具，打印出百度贴吧某用户发表过的所有帖子
  • 3
    关于百度地图api测距显示NaN的解决方案
  • 4
    [超详实]SEO百度优化实战长篇纪实，第一篇：缘起
  • 5
    百度知道发布iPhone应用
  • 6
    百度兴趣点下载工具设计和实现
  • 7
    【android-tips】如何在android应用中插入百度广告（附源码）
  • 8
    百度熊掌号SEO指南1.0，seo重点摘录
  • 9
    王海峰被任命为百度 CTO，李一男曾担任该职位
  • 10
    推送sdk （类似百度 推送 ） android mina (待优化... )
  • 1
    Machine Learning机器学习之决策树算法 Decision Tree（附Python代码）
    301
  • 2
    【数据结构】八大排序之计数排序算法
    39
  • 3
    【数据结构】八大排序之归并排序算法
    44
  • 4
    【数据结构】八大排序之快速排序算法
    91
  • 5
    深入理解InnoDB索引数据结构和算法
    312
  • 6
    Java入门高频考查算法逻辑基础知识3-编程篇（超详细18题1.8万字参考编程实现）
    84
  • 7
    sklearn应用线性回归算法
    77
  • 8
    深度思考：雪花算法snowflake分布式id生成原理详解
    1533
  • 9
    金石推荐 | 【算法数据结构专题】「延时队列算法」史上手把手教你针对层级时间轮（TimingWheel）实现延时队列的开发实战落地（下）（一）
    125
  • 10
    大厂算法指南：优选算法 ——双指针篇（下）
    59

    • 相关课程

    更多
  • 智能运维赛（复赛）：利用数据和算法，快速定位系统异常并进行根因分析
  • 智能创作赛（复赛）：相册应用中的视频故事生成算法介绍
  • 智能创作赛（初赛）：相册应用中的故事生成算法介绍
  • 相册服务中的故事生成算法介绍
  • Go语言核心编程 - 数据结构和算法
  • 神经网络概览及算法详解

    • 相关电子书

    更多
  • AIOps-百度的思考与实践
  • 百度外卖从IDC到云端服务迁移历程
  • 百度研发工具集的应用实践——打造持续交付路上的“三堵墙”

    • 相关实验场景

    更多
  • 使用Swing算法实现商品推荐
  • RSA密码算法设计与实现
  • 欧拉图的构造性证明与算法实现
  • 推荐系统入门之使用ALS算法实现打分预测
    • 下一篇
    阿里云无影云电脑免费试用，最长可试用3个月