GPS、谷歌、百度、高德坐标相互转换

news/2024/10/24 4:48:41/

GPS、谷歌、百度、高德坐标相互转换

一、在进行地图开发过程中，我们一般能接触到以下三种类型的地图坐标系：

1.WGS－84原始坐标系，一般用国际GPS纪录仪记录下来的经纬度，通过GPS定位拿到的原始经纬度，Google和高德地图定位的的经纬度（国外）都是基于WGS－84坐标系的；但是在国内是不允许直接用WGS84坐标系标注的，必须经过加密后才能使用；

2.GCJ－02坐标系，又名“火星坐标系”，是我国国测局独创的坐标体系，由WGS－84加密而成，在国内，必须至少使用GCJ－02坐标系，或者使用在GCJ－02加密后再进行加密的坐标系，如百度坐标系。高德和Google在国内都是使用GCJ－02坐标系，可以说，GCJ－02是国内最广泛使用的坐标系；

3.百度坐标系:bd-09，百度坐标系是在GCJ－02坐标系的基础上再次加密偏移后形成的坐标系，只适用于百度地图。(目前百度API提供了从其它坐标系转换为百度坐标系的API，但却没有从百度坐标系转为其他坐标系的API)

二、为什么会发生偏移？

1.由于坐标系之间不兼容，如在百度地图上定位的经纬度拿到高德地图上直接描点就肯定会发生偏移；只考虑国内的情况，高德地图和Google地图是可以不经过转换也能够准确显示的（在国内用的都是GCJ－02坐标系）；下面是收录了网上的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);  }  
}

c#：

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 bool 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[] Gps84ToGcj02(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[] Gcj02ToGps84(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 };}/// <summary>/// 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 将 GCJ-02 坐标转换成 BD-09 坐标 /// 高德谷歌转为百度/// </summary>/// <param name="lat"></param>/// <param name="lon"></param>/// <returns></returns>public static double[] Gcj02ToBd09(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;}/// <summary>/// 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法 * * 将 BD-09 坐标转换成GCJ-02 坐标/// 百度坐标转为高德谷歌坐标/// </summary>/// <param name="lat"></param>/// <param name="lon"></param>/// <returns></returns>public static double[] Bd09ToGcj02(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;}/// <summary>/// gps84转为bd09/// GPS坐标转为百度坐标/// </summary>/// <param name="lat"></param>/// <param name="lon"></param>/// <returns></returns>public static double[] Gps84ToBd09(double lat, double lon){double[] gcj02 = Gps84ToGcj02(lat, lon);double[] bd09 = Gcj02ToBd09(gcj02[0], gcj02[1]);return bd09;}/// <summary>/// 百度坐标转成GPS坐标/// </summary>/// <param name="lat"></param>/// <param name="lon"></param>/// <returns></returns>public static double[] Bd09ToGps84(double lat, double lon){double[] gcj02 = Bd09ToGcj02(lat, lon);double[] gps84 = Gcj02ToGps84(gcj02[0], gcj02[1]);//保留小数点后六位  gps84[0] = Retain6(gps84[0]);gps84[1] = Retain6(gps84[1]);return gps84;}/// <summary>/// 保留小数点后六位 /// </summary>/// <param name="num"></param>/// <returns></returns>private static double Retain6(double num){String result = String.Format("%.6f", num);return Double.Parse(result);}
}

作者：青檬可乐_
链接：https://www.jianshu.com/p/d3dd4149bb0b
來源：简书
简书著作权归作者所有，任何形式的转载都请联系作者获得授权并注明出处。


http://www.ppmy.cn/news/659580.html

相关文章

高德地图通过经纬度获取地区城市

高德地图通过经纬度获取地区城市

php curl 请求接口 /*** 模拟post进行url请求 * param string $url * param array $postData */ function https_request($url, $data null) {$curl curl_init();//初始化curl_setopt($curl, CURLOPT_URL, $url);curl_setopt($curl, CURLOPT_TIMEOUT, 30);//允许 cURL 函数执…
阅读更多...
高德地图的逆地理编码 | 将经纬度坐标转化为对应的地理位置

高德地图的逆地理编码 | 将经纬度坐标转化为对应的地理位置

官网 地理/逆地理编码-API文档-开发指南-Web服务 API | 高德地图API 官方解释&#xff1a;地理编码/逆地理编码 API 是通过 HTTP/HTTPS 协议访问远程服务的接口&#xff0c;提供结构化地址与经纬度之间的相互转化的能力。 通读易懂也就是将经纬度坐标转化为对应的地理位置信息…
阅读更多...
高德地图获取手机定位以及经纬度组件

高德地图获取手机定位以及经纬度组件

父组件代码&#xff1a; html代码 <div><get-position :position.sync"currentPosition" :place.sync"currentPlace"/> </div>import GetPosition from "../../components/GetPosition" //引入GetPosition组件 export defau…
阅读更多...
uniapp 打包app用不了腾讯地图|高德地图

uniapp 打包app用不了腾讯地图|高德地图

原因1:好像app只支持高德|百度地图 好像app只支持高德|百度地图的,申请一个高德地图的key就可 SHA1&#xff1a; 1.dcolud开发平台地址:开发者中心 2.查看SHA1&#xff08;点击应用名称进入&#xff09; 3.查看证书 包名: key: 原因2:配置权限及设置key 1.配置key 2.配置权…
阅读更多...
高德地图获取可视区域内四角坐标（东北，东南，西南，西北）

高德地图获取可视区域内四角坐标（东北，东南，西南，西北）

高德地图获取可视区域内四角坐标&#xff08;东北&#xff0c;东南&#xff0c;西南&#xff0c;西北&#xff09; 场景&#xff1a;当需要实现查询某坐标区域内的所有共享单车&#xff0c;那么常规操作就是把所选区域内的左上和右下两点经纬度给到后端。获取方法如下&#xff…
阅读更多...
百度经纬度与高德经纬度互转

百度经纬度与高德经纬度互转

使用不同的地图&#xff0c;在不同的地图上根据经纬度展示会有出入&#xff0c; 这边记录下 百度和高德经纬度的互转 // 将百度地图经纬度转换为腾讯/高德地图经纬度 用于大屏bMapTransQQMap(lng, lat) {let x_pi 3.14159265358979324 * 3000.0 / 180.0;let x lng - …
阅读更多...
知道经纬度来调高德地图的官网API来获取所在的位置(逆地理编码)

知道经纬度来调高德地图的官网API来获取所在的位置(逆地理编码)

第一步&#xff1a;首先要有一个key值 高德的参考 https://lbs.amap.com/api/webservice/guide/api/georegeo/?sug_index2 第二步&#xff1a;调接口 https://restapi.amap.com/v3/geocode/regeo?key"你的key值"&location经度,纬度 注&#xff1a;key值要…
阅读更多...
高德地图一些常用的地理位置获取方法

高德地图一些常用的地理位置获取方法

1、首先引入plugin //webapi.amap.com/maps?v1.4.15&key000&pluginAMap.CustomLayer,AMap.Geocoder 2、mixin import { uuid } from /utils/utils; export default {data() {return {map: null,mapStyle: amap://styles/bcb8c05d5306c3e9004a513721f93262, //自定义…
阅读更多...
最新文章

Copyright @ 2022~2023