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。

相关文章

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

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

《关于集合与数组转换实现方法》：本文主要介绍关于集合与数组转换实现方法,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1、Arrays.asList()1.1、方法作用1.2、内部实现1.3、修改元素的影响1.4、注意事项2、list.toArray()2.1、方
阅读更多...
Golang如何对cron进行二次封装实现指定时间执行定时任务

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

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

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

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

使用Python进行GRPC和Dubbo协议的高级测试

《使用Python进行GRPC和Dubbo协议的高级测试》GRPC（GoogleRemoteProcedureCall）是一种高性能、开源的远程过程调用（RPC）框架,Dubbo是一种高性能的分布式服... 目录01 GRPC测试安装gRPC编写.proto文件实现服务02 Dubbo测试1. 安装Dubb
阅读更多...
java Long 与long之间的转换流程

java Long 与long之间的转换流程

《javaLong与long之间的转换流程》Long类提供了一些方法,用于在long和其他数据类型（如String）之间进行转换,本文将详细介绍如何在Java中实现Long和long之间的转换,感... 目录概述流程步骤1：将long转换为Long对象步骤2：将Longhttp://www.cppcns.c
阅读更多...
在Java中将XLS转换为XLSX的实现方案

在Java中将XLS转换为XLSX的实现方案

《在Java中将XLS转换为XLSX的实现方案》在本文中,我们将探讨传统ExcelXLS格式与现代XLSX格式的结构差异,并为Java开发者提供转换方案,通过了解底层原理、性能优势及实用工具,您将掌握... 目录为什么升级XLS到XLSX值得投入？实际转换过程解析推荐技术方案对比Apache POI实现编程
阅读更多...
Linux使用scp进行远程目录文件复制的详细步骤和示例

Linux使用scp进行远程目录文件复制的详细步骤和示例

《Linux使用scp进行远程目录文件复制的详细步骤和示例》在Linux系统中,scp（安全复制协议）是一个使用SSH（安全外壳协议）进行文件和目录安全传输的命令,它允许在远程主机之间复制文件和目录,... 目录1. 什么是scp？2. 语法3. 示例示例 1: 复制本地目录到远程主机示例 2: 复制远程主
阅读更多...
C#代码实现解析WTGPS和BD数据

C#代码实现解析WTGPS和BD数据

《C#代码实现解析WTGPS和BD数据》在现代的导航与定位应用中,准确解析GPS和北斗（BD）等卫星定位数据至关重要,本文将使用C#语言实现解析WTGPS和BD数据,需要的可以了解下... 目录一、代码结构概览1. 核心解析方法2. 位置信息解析3. 经纬度转换方法4. 日期和时间戳解析5. 辅助方法二、L
阅读更多...
windows系统上如何进行maven安装和配置方式

windows系统上如何进行maven安装和配置方式

《windows系统上如何进行maven安装和配置方式》：本文主要介绍windows系统上如何进行maven安装和配置方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不... 目录1. Maven 简介2. maven的下载与安装2.1 下载 Maven2.2 Maven安装2.
阅读更多...
C/C++的OpenCV 进行图像梯度提取的几种实现

C/C++的OpenCV 进行图像梯度提取的几种实现

《C/C++的OpenCV进行图像梯度提取的几种实现》本文主要介绍了C/C++的OpenCV进行图像梯度提取的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的... 目录预www.chinasem.cn备知识1. 图像加载与预处理2. Sobel 算子计算 X 和 Y
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2