WGS-84 ,GCJ-02与BD-09的坐标进行转换

2024-01-28 15:58
文章标签 进行 02 坐标 转换 09 84 bd gcj wgs

本文主要是介绍WGS-84 ,GCJ-02与BD-09的坐标进行转换,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

Java实现

package com.zehin.map.util;public class GPSUtil {  public static double pi = 3.1415926535897932384626;  public static double x_pi = 3.14159265358979324 * 3000.0 / 180.0;  public static double a = 6378245.0;  public static double ee = 0.00669342162296594323;  public static double transformLat(double x, double y) {  double ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y  + 0.2 * Math.sqrt(Math.abs(x));  ret += (20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0 / 3.0;  ret += (20.0 * Math.sin(y * pi) + 40.0 * Math.sin(y / 3.0 * pi)) * 2.0 / 3.0;  ret += (160.0 * Math.sin(y / 12.0 * pi) + 320 * Math.sin(y * pi / 30.0)) * 2.0 / 3.0;  return ret;  }  public static double transformLon(double x, double y) {  double ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1  * Math.sqrt(Math.abs(x));  ret += (20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0 / 3.0;  ret += (20.0 * Math.sin(x * pi) + 40.0 * Math.sin(x / 3.0 * pi)) * 2.0 / 3.0;  ret += (150.0 * Math.sin(x / 12.0 * pi) + 300.0 * Math.sin(x / 30.0  * pi)) * 2.0 / 3.0;  return ret;  }  public static double[] transform(double lat, double lon) {  if (outOfChina(lat, lon)) {  return new double[]{lat,lon};  }  double dLat = transformLat(lon - 105.0, lat - 35.0);  double dLon = transformLon(lon - 105.0, lat - 35.0);  double radLat = lat / 180.0 * pi;  double magic = Math.sin(radLat);  magic = 1 - ee * magic * magic;  double sqrtMagic = Math.sqrt(magic);  dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);  dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);  double mgLat = lat + dLat;  double mgLon = lon + dLon;  return new double[]{mgLat,mgLon};  }  public static boolean outOfChina(double lat, double lon) {  if (lon < 72.004 || lon > 137.8347)  return true;  if (lat < 0.8293 || lat > 55.8271)  return true;  return false;  }  /** * 84 to 火星坐标系 (GCJ-02) World Geodetic System ==> Mars Geodetic System * * @param lat * @param lon * @return */  public static double[] gps84_To_Gcj02(double lat, double lon) {  if (outOfChina(lat, lon)) {  return new double[]{lat,lon};  }  double dLat = transformLat(lon - 105.0, lat - 35.0);  double dLon = transformLon(lon - 105.0, lat - 35.0);  double radLat = lat / 180.0 * pi;  double magic = Math.sin(radLat);  magic = 1 - ee * magic * magic;  double sqrtMagic = Math.sqrt(magic);  dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);  dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * pi);  double mgLat = lat + dLat;  double mgLon = lon + dLon;  return new double[]{mgLat, mgLon};  }  /** * * 火星坐标系 (GCJ-02) to 84 * * @param lon * @param lat * @return * */  public static double[] gcj02_To_Gps84(double lat, double lon) {  double[] gps = transform(lat, lon);  double lontitude = lon * 2 - gps[1];  double latitude = lat * 2 - gps[0];  return new double[]{latitude, lontitude};  }  /** * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 将 GCJ-02 坐标转换成 BD-09 坐标 * * @param lat * @param lon */  public static double[] gcj02_To_Bd09(double lat, double lon) {  double x = lon, 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 tempLon = z * Math.cos(theta) + 0.0065;  double tempLat = z * Math.sin(theta) + 0.006;  double[] gps = {tempLat,tempLon};  return gps;  }  /** * * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 * * 将 BD-09 坐标转换成GCJ-02 坐标 * * @param * bd_lat * @param bd_lon * @return */  public static double[] bd09_To_Gcj02(double lat, double lon) {  double x = lon - 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 tempLon = z * Math.cos(theta);  double tempLat = z * Math.sin(theta);  double[] gps = {tempLat,tempLon};  return gps;  }  /**将gps84转为bd09 * @param lat * @param lon * @return */  public static double[] gps84_To_bd09(double lat,double lon){  double[] gcj02 = gps84_To_Gcj02(lat,lon);  double[] bd09 = gcj02_To_Bd09(gcj02[0],gcj02[1]);  return bd09;  }  public static double[] bd09_To_gps84(double lat,double lon){  double[] gcj02 = bd09_To_Gcj02(lat, lon);  double[] gps84 = gcj02_To_Gps84(gcj02[0], gcj02[1]);  //保留小数点后六位  gps84[0] = retain6(gps84[0]);  gps84[1] = retain6(gps84[1]);  return gps84;  }  /**保留小数点后六位 * @param num * @return */  private static double retain6(double num){  String result = String .format("%.6f", num);  return Double.valueOf(result);  }  } 

js实现

var GPS = {PI : 3.14159265358979324,x_pi : 3.14159265358979324 * 3000.0 / 180.0,delta : function (lat, lon) {// Krasovsky 1940//// a = 6378245.0, 1/f = 298.3// b = a * (1 - f)// ee = (a^2 - b^2) / a^2;var a = 6378245.0; //  a: 卫星椭球坐标投影到平面地图坐标系的投影因子。var ee = 0.00669342162296594323; //  ee: 椭球的偏心率。var dLat = this.transformLat(lon - 105.0, lat - 35.0);var dLon = this.transformLon(lon - 105.0, lat - 35.0);var radLat = lat / 180.0 * this.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) * this.PI);dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) * this.PI);return {'lat': dLat, 'lon': dLon};},//WGS-84 to GCJ-02gcj_encrypt : function (wgsLat, wgsLon) {if (this.outOfChina(wgsLat, wgsLon))return {'lat': wgsLat, 'lon': wgsLon};var d = this.delta(wgsLat, wgsLon);return {'lat' : wgsLat + d.lat,'lon' : wgsLon + d.lon};},//GCJ-02 to WGS-84gcj_decrypt : function (gcjLat, gcjLon) {if (this.outOfChina(gcjLat, gcjLon))return {'lat': gcjLat, 'lon': gcjLon};var d = this.delta(gcjLat, gcjLon);return {'lat': gcjLat - d.lat, 'lon': gcjLon - d.lon};},//GCJ-02 to WGS-84 exactlygcj_decrypt_exact : function (gcjLat, gcjLon) {var initDelta = 0.01;var threshold = 0.000000001;var dLat = initDelta, dLon = initDelta;var mLat = gcjLat - dLat, mLon = gcjLon - dLon;var pLat = gcjLat + dLat, pLon = gcjLon + dLon;var wgsLat, wgsLon, i = 0;while (1) {wgsLat = (mLat + pLat) / 2;wgsLon = (mLon + pLon) / 2;var tmp = this.gcj_encrypt(wgsLat, wgsLon)dLat = tmp.lat - gcjLat;dLon = tmp.lon - gcjLon;if ((Math.abs(dLat) < threshold) && (Math.abs(dLon) < threshold))break;if (dLat > 0) pLat = wgsLat; else mLat = wgsLat;if (dLon > 0) pLon = wgsLon; else mLon = wgsLon;if (++i > 10000) break;}//console.log(i);return {'lat': wgsLat, 'lon': wgsLon};},//GCJ-02 to BD-09bd_encrypt : function (gcjLat, gcjLon) {var x = gcjLon, y = gcjLat;var z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * this.x_pi);var theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * this.x_pi);bdLon = z * Math.cos(theta) + 0.0065;bdLat = z * Math.sin(theta) + 0.006;return {'lat' : bdLat,'lon' : bdLon};},//BD-09 to GCJ-02bd_decrypt : function (bdLat, bdLon) {var x = bdLon - 0.0065, y = bdLat - 0.006;var z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * this.x_pi);var theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * this.x_pi);var gcjLon = z * Math.cos(theta);var gcjLat = z * Math.sin(theta);return {'lat' : gcjLat, 'lon' : gcjLon};},//WGS-84 to Web mercator//mercatorLat -> y mercatorLon -> xmercator_encrypt : function(wgsLat, wgsLon) {var x = wgsLon * 20037508.34 / 180.;var y = Math.log(Math.tan((90. + wgsLat) * this.PI / 360.)) / (this.PI / 180.);y = y * 20037508.34 / 180.;return {'lat' : y, 'lon' : x};/*if ((Math.abs(wgsLon) > 180 || Math.abs(wgsLat) > 90))return null;var x = 6378137.0 * wgsLon * 0.017453292519943295;var a = wgsLat * 0.017453292519943295;var y = 3189068.5 * Math.log((1.0 + Math.sin(a)) / (1.0 - Math.sin(a)));return {'lat' : y, 'lon' : x};//*/},// Web mercator to WGS-84// mercatorLat -> y mercatorLon -> xmercator_decrypt : function(mercatorLat, mercatorLon) {var x = mercatorLon / 20037508.34 * 180.;var y = mercatorLat / 20037508.34 * 180.;y = 180 / this.PI * (2 * Math.atan(Math.exp(y * this.PI / 180.)) - this.PI / 2);return {'lat' : y, 'lon' : x};/*if (Math.abs(mercatorLon) < 180 && Math.abs(mercatorLat) < 90)return null;if ((Math.abs(mercatorLon) > 20037508.3427892) || (Math.abs(mercatorLat) > 20037508.3427892))return null;var a = mercatorLon / 6378137.0 * 57.295779513082323;var x = a - (Math.floor(((a + 180.0) / 360.0)) * 360.0);var y = (1.5707963267948966 - (2.0 * Math.atan(Math.exp((-1.0 * mercatorLat) / 6378137.0)))) * 57.295779513082323;return {'lat' : y, 'lon' : x};//*/},// two point's distancedistance : function (latA, lonA, latB, lonB) {var earthR = 6371000.;var x = Math.cos(latA * this.PI / 180.) * Math.cos(latB * this.PI / 180.) * Math.cos((lonA - lonB) * this.PI / 180);var y = Math.sin(latA * this.PI / 180.) * Math.sin(latB * this.PI / 180.);var s = x + y;if (s > 1) s = 1;if (s < -1) s = -1;var alpha = Math.acos(s);var distance = alpha * earthR;return distance;},outOfChina : function (lat, lon) {if (lon < 72.004 || lon > 137.8347)return true;if (lat < 0.8293 || lat > 55.8271)return true;return false;},transformLat : function (x, y) {var ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x));ret += (20.0 * Math.sin(6.0 * x * this.PI) + 20.0 * Math.sin(2.0 * x * this.PI)) * 2.0 / 3.0;ret += (20.0 * Math.sin(y * this.PI) + 40.0 * Math.sin(y / 3.0 * this.PI)) * 2.0 / 3.0;ret += (160.0 * Math.sin(y / 12.0 * this.PI) + 320 * Math.sin(y * this.PI / 30.0)) * 2.0 / 3.0;return ret;},transformLon : function (x, y) {var ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.sqrt(Math.abs(x));ret += (20.0 * Math.sin(6.0 * x * this.PI) + 20.0 * Math.sin(2.0 * x * this.PI)) * 2.0 / 3.0;ret += (20.0 * Math.sin(x * this.PI) + 40.0 * Math.sin(x / 3.0 * this.PI)) * 2.0 / 3.0;ret += (150.0 * Math.sin(x / 12.0 * this.PI) + 300.0 * Math.sin(x / 30.0 * this.PI)) * 2.0 / 3.0;return ret;}
};

这篇关于WGS-84 ,GCJ-02与BD-09的坐标进行转换的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/654096

相关文章

如何使用Lombok进行spring 注入

《如何使用Lombok进行spring注入》本文介绍如何用Lombok简化Spring注入,推荐优先使用setter注入,通过注解自动生成getter/setter及构造器,减少冗余代码,提升开发效... Lombok为了开发环境简化代码,好处不用多说。spring 注入方式为2种,构造器注入和setter

MySQL进行数据库审计的详细步骤和示例代码

《MySQL进行数据库审计的详细步骤和示例代码》数据库审计通过触发器、内置功能及第三方工具记录和监控数据库活动,确保安全、完整与合规,Java代码实现自动化日志记录,整合分析系统提升监控效率,本文给大... 目录一、数据库审计的基本概念二、使用触发器进行数据库审计1. 创建审计表2. 创建触发器三、Java

Kotlin Map映射转换问题小结

《KotlinMap映射转换问题小结》文章介绍了Kotlin集合转换的多种方法,包括map(一对一转换)、mapIndexed(带索引)、mapNotNull(过滤null)、mapKeys/map... 目录Kotlin 集合转换:map、mapIndexed、mapNotNull、mapKeys、map

MySQL深分页进行性能优化的常见方法

《MySQL深分页进行性能优化的常见方法》在Web应用中,分页查询是数据库操作中的常见需求,然而,在面对大型数据集时,深分页(deeppagination)却成为了性能优化的一个挑战,在本文中,我们将... 目录引言:深分页,真的只是“翻页慢”那么简单吗?一、背景介绍二、深分页的性能问题三、业务场景分析四、

SpringBoot结合Docker进行容器化处理指南

《SpringBoot结合Docker进行容器化处理指南》在当今快速发展的软件工程领域,SpringBoot和Docker已经成为现代Java开发者的必备工具,本文将深入讲解如何将一个SpringBo... 目录前言一、为什么选择 Spring Bootjavascript + docker1. 快速部署与

linux解压缩 xxx.jar文件进行内部操作过程

《linux解压缩xxx.jar文件进行内部操作过程》:本文主要介绍linux解压缩xxx.jar文件进行内部操作,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、解压文件二、压缩文件总结一、解压文件1、把 xxx.jar 文件放在服务器上,并进入当前目录#

SpringBoot中如何使用Assert进行断言校验

《SpringBoot中如何使用Assert进行断言校验》Java提供了内置的assert机制,而Spring框架也提供了更强大的Assert工具类来帮助开发者进行参数校验和状态检查,下... 目录前言一、Java 原生assert简介1.1 使用方式1.2 示例代码1.3 优缺点分析二、Spring Fr

关于集合与数组转换实现方法

《关于集合与数组转换实现方法》:本文主要介绍关于集合与数组转换实现方法,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1、Arrays.asList()1.1、方法作用1.2、内部实现1.3、修改元素的影响1.4、注意事项2、list.toArray()2.1、方

Golang如何对cron进行二次封装实现指定时间执行定时任务

《Golang如何对cron进行二次封装实现指定时间执行定时任务》:本文主要介绍Golang如何对cron进行二次封装实现指定时间执行定时任务问题,具有很好的参考价值,希望对大家有所帮助,如有错误... 目录背景cron库下载代码示例【1】结构体定义【2】定时任务开启【3】使用示例【4】控制台输出总结背景

利用Python脚本实现批量将图片转换为WebP格式

《利用Python脚本实现批量将图片转换为WebP格式》Python语言的简洁语法和库支持使其成为图像处理的理想选择,本文将介绍如何利用Python实现批量将图片转换为WebP格式的脚本,WebP作为... 目录简介1. python在图像处理中的应用2. WebP格式的原理和优势2.1 WebP格式与传统