背景:从GPS和北斗卫星定位得到的定位数据采用的都是WGS84坐标系,即地球坐标系,但是国内不管是高德地图、百度地图采用的并不是WGS84坐标系,所以需要经过转换后才能使用,前端用百度API提供的方法转换速度较慢。通过网上搜索收集资料,然后自己编写代码测试,总结了一下有关坐标系简介与转换的知识,分享给大家!

文章目录

  • 1.各坐标系简介
  • 2.各坐标系转换
    • 2.1坐标点实体类
    • 2.2各坐标系转换工具类
  • 3.测试

1.各坐标系简介

  • WGS84坐标系
    即地球坐标系,国际上通用的坐标系。
    设备一般包含GPS芯片或者北斗芯片获取的经纬度为WGS84地理坐标系。谷歌地图采用的是WGS84地理坐标系(中国范围除外,谷歌中国地图采用的是GCJ02地理坐标系。)

  • GCJ02坐标系
    即火星坐标系,WGS84坐标系经加密后的坐标系。
    出于国家安全考虑,国内所有导航电子地图必须使用国家测绘局制定的加密坐标系统,即将一个真实的经纬度坐标加密成一个不正确的经纬度坐标。

  • BD09坐标系
    即百度坐标系,GCJ02坐标系经加密后的坐标系。搜狗坐标系、图吧坐标系等,估计也是在GCJ02基础上加密而成的。

  • 各主流地图API采用的坐标系
    高德MapABC地图API 火星坐标
    腾讯搜搜地图API 火星坐标
    阿里云地图API 火星坐标
    灵图51ditu地图API 火星坐标

    百度地图API 百度坐标
    搜狐搜狗地图API 搜狗坐标
    图吧MapBar地图API 图吧坐标

  • 这里介绍几个常用的在线坐标工具
    高德开放平台
    百度拾取坐标系统
    在线坐标转换工具

2.各坐标系转换

2.1坐标点实体类

/*** Description: 坐标点* * @author JourWon* @date Created on 2018年6月19日*/
public class Point implements Serializable {private static final long serialVersionUID = 3584864663880053897L;/*** 经度*/private double lng;/*** 纬度*/private double lat;

2.2各坐标系转换工具类

package com.jourwon.util;import com.jourwon.pojo.Point;/*** Description: 各坐标系之间的转换工具类* * @author JourWon* @date Created on 2018年6月19日*/
public class CoordinateTransformUtils {// 圆周率πprivate static final double PI = 3.1415926535897932384626D;// 火星坐标系与百度坐标系转换的中间量private static final double X_PI = 3.14159265358979324 * 3000.0 / 180.0D;// Krasovsky 1940// 长半轴a = 6378245.0, 1/f = 298.3// b = a * (1 - f)// 扁率ee = (a^2 - b^2) / a^2;// 长半轴private static final double SEMI_MAJOR = 6378245.0D;// 扁率private static final double FLATTENING = 0.00669342162296594323D;// WGS84=>GCJ02 地球坐标系=>火星坐标系public static Point wgs84ToGcj02(double lng, double lat) {if (outOfChina(lng, lat)) {return new Point(lng, lat);}double[] offset = offset(lng, lat);double mglng = lng + offset[0];double mglat = lat + offset[1];return new Point(mglng, mglat);}// GCJ02=>WGS84 火星坐标系=>地球坐标系(粗略)public static Point gcj02ToWgs84(double lng, double lat) {if (outOfChina(lng, lat)) {return new Point(lng, lat);}double[] offset = offset(lng, lat);double mglng = lng - offset[0];double mglat = lat - offset[1];return new Point(mglng, mglat);}// GCJ02=>WGS84 火星坐标系=>地球坐标系(精确)public static Point gcj02ToWgs84Exactly(double lng, double lat) {if (outOfChina(lng, lat)) {return new Point(lng, lat);}double initDelta = 0.01;double threshold = 0.000000001;double dLat = initDelta, dLon = initDelta;double mLat = lat - dLat, mLon = lng - dLon;double pLat = lat + dLat, pLon = lng + dLon;double wgsLat, wgsLng, i = 0;while (true) {wgsLat = (mLat + pLat) / 2;wgsLng = (mLon + pLon) / 2;Point point = wgs84ToGcj02(wgsLng, wgsLat);dLon = point.getLng() - lng;dLat = point.getLat() - lat;if ((Math.abs(dLat) < threshold) && (Math.abs(dLon) < threshold))break;if (dLat > 0)pLat = wgsLat;elsemLat = wgsLat;if (dLon > 0)pLon = wgsLng;elsemLon = wgsLng;if (++i > 10000)break;}return new Point(wgsLng, wgsLat);}// GCJ-02=>BD09 火星坐标系=>百度坐标系public static Point gcj02ToBd09(double lng, double lat) {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 bd_lng = z * Math.cos(theta) + 0.0065;double bd_lat = z * Math.sin(theta) + 0.006;return new Point(bd_lng, bd_lat);}// BD09=>GCJ-02 百度坐标系=>火星坐标系public static Point bd09ToGcj02(double lng, double lat) {double x = lng - 0.0065;double 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 gcj_lng = z * Math.cos(theta);double gcj_lat = z * Math.sin(theta);return new Point(gcj_lng, gcj_lat);}// WGS84=>BD09 地球坐标系=>百度坐标系public static Point wgs84ToBd09(double lng, double lat) {Point point = wgs84ToGcj02(lng, lat);return gcj02ToBd09(point.getLng(), point.getLat());}// BD09=>WGS84 百度坐标系=>地球坐标系public static Point bd09ToWgs84(double lng, double lat) {Point point = bd09ToGcj02(lng, lat);return gcj02ToWgs84(point.getLng(), point.getLat());}/*** Description: 中国境外返回true,境内返回false* @param lng    经度* @param lat 纬度* @return*/public static boolean outOfChina(double lng, double lat) {if (lng < 72.004 || lng > 137.8347)return true;if (lat < 0.8293 || lat > 55.8271)return true;return false;}// 经度偏移量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 * 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;}// 纬度偏移量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 * 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;}// 偏移量public static double[] offset(double lng, double lat) {double[] lngLat = new double[2];double dlng = transformLng(lng - 105.0, lat - 35.0);double dlat = transformLat(lng - 105.0, lat - 35.0);double radlat = lat / 180.0 * PI;double magic = Math.sin(radlat);magic = 1 - FLATTENING * magic * magic;double sqrtmagic = Math.sqrt(magic);dlng = (dlng * 180.0) / (SEMI_MAJOR / sqrtmagic * Math.cos(radlat) * PI);dlat = (dlat * 180.0) / ((SEMI_MAJOR * (1 - FLATTENING)) / (magic * sqrtmagic) * PI);lngLat[0] = dlng;lngLat[1] = dlat;return lngLat;}}

3.测试

/*** Description: 坐标系转换工具测试类* * @author JourWon* @date Created on 2018年6月19日*/
public class CoordinateTransformUtilsTest {/*** Description: 地球坐标系 =>火星坐标系、百度坐标系*/@Testpublic void test01() {// 广州市中大地铁站Point point = new Point(113.28749670783887D, 23.094783676708065D);System.out.println("地球坐标系 : " + point);Point wgs84ToGcj02 = CoordinateTransformUtils.wgs84ToGcj02(point.getLng(), point.getLat());System.out.println("火星坐标系 : " + wgs84ToGcj02);Point wgs84ToBd09 = CoordinateTransformUtils.wgs84ToBd09(point.getLng(), point.getLat());System.out.println("百度坐标系 : " + wgs84ToBd09);}/*** Description: 火星坐标系=>地球坐标系、百度坐标系*/@Testpublic void test02() {// 广州市珠江新城地铁站Point point = new Point(113.321171D, 23.119285D);System.out.println("火星坐标系 : " + point);Point wgs84ToGcj02 = CoordinateTransformUtils.gcj02ToWgs84(point.getLng(), point.getLat());System.out.println("地球坐标系 : " + wgs84ToGcj02);Point wgs84ToBd09 = CoordinateTransformUtils.gcj02ToBd09(point.getLng(), point.getLat());System.out.println("百度坐标系 : " + wgs84ToBd09);}/*** Description: 百度坐标系=>地球坐标系、火星坐标系*/@Testpublic void test03() {// 广州市体育西路地铁站Point point = new Point(113.328035D, 23.136929D);System.out.println("百度坐标系 : " + point);Point wgs84ToGcj02 = CoordinateTransformUtils.bd09ToWgs84(point.getLng(), point.getLat());System.out.println("地球坐标系 : " + wgs84ToGcj02);Point wgs84ToBd09 = CoordinateTransformUtils.bd09ToGcj02(point.getLng(), point.getLat());System.out.println("火星坐标系 : " + wgs84ToBd09);}}

完整代码下载链接

WGS84地球坐标系,GCJ02火星坐标系,BD09百度坐标系简介与转换相关推荐

  1. WGS84地球坐标系,GCJ02火星坐标系,BD09百度坐标系简介与转换,mybatis字段映射原理

    1.各坐标系简介 2.各坐标系转换 2.1坐标点实体类 2.2各坐标系转换工具类 3.测试 1.各坐标系简介 WGS84坐标系 即地球坐标系,国际上通用的坐标系. 设备一般包含GPS芯片或者北斗芯片获 ...

  2. GPS坐标系转换(标准坐标系WGS84转GCJ-02火星坐标系)

    GPS坐标系转换(标准坐标系WGS84转GCJ-02火星坐标系) 坐标系简介 WGS-84 标准坐标系 GCJ-02 - 国测局坐标 BD-09 - 百度坐标系 干货前的说明 WGS84转GCJ-02 ...

  3. 基于Postgresql和PostGIS实现火星坐标系、百度坐标系、WGS84坐标系、CGCS2000坐标系互转

    背景 最近有一个需求,需要将WGS84转成火星坐标系.个人觉得在代码中逐个点坐标进行转换,太麻烦,而且效率低.PostGIS的st_transform虽然可以进行坐标转换,但是不支持国内这些坐标系.最 ...

  4. GCJ-02火星坐标系、BD-09百度坐标系和WGS-84坐标系转换

    WGS-84:GPS坐标系 GCJ-02:火星坐标系,国测局02年发布的坐标体系,高德,腾讯等使用. BD-09:百度坐标系 GCJ-02火星坐标系和WGS-84坐标系转换关系 var pi = 3. ...

  5. 地球坐标系,火星坐标系,百度坐标系

    地球坐标系--WGS84:常见于 GPS 设备,Google 地图等国际标准的坐标体系. 火星坐标系--GCJ-02:中国国内使用的被强制加密后的坐标体系,高德坐标就属于该种坐标体系. 百度坐标系-- ...

  6. 地图坐标转换(火星坐标系、百度坐标系、谷歌、腾讯、高德等坐标)

    好文章记录 概念解释 首先搞清楚什么是火星坐标?到底有哪些坐标体系?由来是怎样的? 1.1 名词解释 坐标系统:用于定位的系统,就跟二维笛卡尔坐标系统一样,一个点使用(x,y),就能确定该点在笛卡尔坐 ...

  7. php百度坐标转腾讯坐标,火星坐标系、百度坐标系、腾讯坐标系的相互转换

    /** * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 将 GCJ-02 坐标转换成 BD-09 坐标 * * @param gg_lat * @param gg_lon ...

  8. JS火星坐标系转百度坐标系(公式法)

    使用该公式即可将获取到的火星坐标系坐标转成百度坐标系,而无需在调用百度API进行转换 但该公式会有一定的误差,本人测试过后发现误差大概在300~400m. function ggToBaidu(lng ...

  9. 火星坐标和百度坐标之间的转换算法

    在开始这个题目之前,先给大家再次扫扫盲,扫的不是坐标系统的盲,而是我们国家所使用的坐标系统.大家都知道,美国GPS使用的是WGS84的坐标系统,以经纬度的形式来表示地球平面上的某一个位置,这应该是国际 ...

  10. android gps 火星坐标,GPS真实坐标与火星地图坐标/百度地图坐标的转换

    #include #include #include static const uint32_t GPSBaud = 9600; TinyGPSPlus gps; HardwareSerial ss( ...

最新文章

  1. python入门经典例题-Python入门经典练习题
  2. EmEditor编辑器正则表达式的优点
  3. 【Matlab 图像】边缘检测算法及效果演示
  4. Java黑皮书课后题第10章:*10.5(显示素数因子)编写一个程序,提示用户输入一个正整数,然后以降序显示它的所有最小因子
  5. 分享:一款快速传输大文件、支持英文版的工具
  6. mybatis 使用merge into
  7. linux显示系统信息软件下载,linux查看系统信息软件安装信息命令学习笔记
  8. 阿里云Kubernetes服务上使用Tekton完成应用发布初体验
  9. 基于阿里云Serverless架构下函数计算的最新应用场景详解(二)
  10. C#LeetCode刷题之#59-螺旋矩阵 II(Spiral Matrix II)
  11. EL表达式JSON应用
  12. keras中文文档_Keras分词器Tokenizer
  13. Java解决找不到主类或者无法加载主类
  14. Kafka follower副本为什么不对外提供服务
  15. “属你最牛” | 树莓派百变打印机
  16. 智力题----囚犯猜帽子问题
  17. 宁芝普拉姆键盘说明书以及键盘校验软件
  18. opencv(python)------二值化阙值(threshold)、自适应阙值(adaptiveThreshold)
  19. 如何使用智能默认值来减少认知负担
  20. 哈工大计算机学院答辩,翟天临们为啥让人愤怒?看看哈工大这场博士答辩就知道了……...

热门文章

  1. TCON控制字及TMOD寄存器
  2. 华为安装gsm框架_华为Mate30Pro怎么安装谷歌服务框架?谷歌服务GMS框架安装方法...
  3. 智慧城管三维可视化决策系统平台(数字孪生)-解决方案开发案例
  4. 联合概率分布与边缘分布
  5. Java小项目(三)---汽车租赁系统(面向对象)
  6. win7/win10 密码忘了?没关系,利用5次shift,9步轻松破解密码
  7. ArcGIS Pro + PS 矢量化用地规划图
  8. 食品行业仓储条码管理系统解决方案
  9. Tomcat9安装及配置步骤
  10. S5p4418平台AP6212 WIFI稳定性的解决过程