在日常工作学习中常会涉及到WGS84、高德（火星/谷歌）、百度三种空间坐标系的坐标转换，本文将通过python、js两种语言实现坐标系的转换。

　坐标系说明：

  1. 　　wgs84:为一种大地坐标系，也是目前广泛使用的GPS全球卫星定位系统使用的坐标系。
  2. 　　高德/谷歌：采用GCJ－02坐标系，又名“火星坐标系”，是我国国测局独创的坐标体系，由WGS－84加密而成，在国内，必须至少使用GCJ－02坐标系，或者使用在GCJ－02加密后再进行加密的坐标系，如百度坐标系。高德和Google在国内都是使用GCJ－02坐标系，可以说，GCJ－02是国内最广泛使用的坐标系；
  3. 　　百度:为百度坐标系，在GCJ02坐标系基础上再次加密。其中bd09ll表示百度经纬度坐标，bd09mc表示百度墨卡托米制坐标。详情
  • python代码
# * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换

# * 即 百度 转 谷歌、高德

# * @param bd_lon

# * @param bd_lat

# * @returns {*[]}

# */

def bd09togcj02(bd_lon, bd_lat):

    x_pi = 3.14159265358979324 * 3000.0 / 180.0

    x = bd_lon - 0.0065

    y = bd_lat - 0.006

    z = math.sqrt(x * x + y * y) - 0.00002 * math.sin(y * x_pi)

    theta = math.atan2(y, x) - 0.000003 * math.cos(x * x_pi)

    gg_lng = z * math.cos(theta)

    gg_lat = z * math.sin(theta)

    return [gg_lng, gg_lat]

# * 火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换

# * 即谷歌、高德 转 百度

# */

def gcj02tobd09(lng, lat):

    x_PI = 3.14159265358979324 * 3000.0 / 180.0

    z = math.sqrt(lng * lng + lat * lat) + 0.00002 * math.sin(lat * x_PI)

    theta = math.atan2(lat, lng) + 0.000003 * math.cos(lng * x_PI)

    bd_lng = z * math.cos(theta) + 0.0065

    bd_lat = z * math.sin(theta) + 0.006

    return [bd_lng, bd_lat]

# wgs84转高德

def wgs84togcj02(lng, lat):

    PI = 3.1415926535897932384626

    ee = 0.00669342162296594323

    a = 6378245.0

    dlat = transformlat(lng - 105.0, lat - 35.0)

    dlng = transformlng(lng - 105.0, lat - 35.0)

    radlat = lat / 180.0 * PI

    magic = math.sin(radlat)

    magic = 1 - ee * magic * magic

    sqrtmagic = math.sqrt(magic)

    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI)

    dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * PI)

    mglat = lat + dlat

    mglng = lng + dlng

    return [mglng, mglat]

# GCJ02/谷歌、高德 转换为 WGS84 gcj02towgs84

def gcj02towgs84(localStr):

    lng = float(localStr.split(',')[0])

    lat = float(localStr.split(',')[1])

    PI = 3.1415926535897932384626

    ee = 0.00669342162296594323

    a = 6378245.0

    dlat = transformlat(lng - 105.0, lat - 35.0)

    dlng = transformlng(lng - 105.0, lat - 35.0)

    radlat = lat / 180.0 * PI

    magic = math.sin(radlat)

    magic = 1 - ee * magic * magic

    sqrtmagic = math.sqrt(magic)

    dlat = (dlat * 180.0) / ((a * (1 - ee)) / (magic * sqrtmagic) * PI)

    dlng = (dlng * 180.0) / (a / sqrtmagic * math.cos(radlat) * PI)

    mglat = lat + dlat

    mglng = lng + dlng

    return str(lng * 2 - mglng) + ',' + str(lat * 2 - mglat)

def transformlat(lng, lat):

    PI = 3.1415926535897932384626

    ret = -100.0 + 2.0 * lng + 3.0 * lat + 0.2 * lat * \

          lat + 0.1 * lng * lat + 0.2 * math.sqrt(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

def transformlng(lng, lat):

    PI = 3.1415926535897932384626

    ret = 300.0 + lng + 2.0 * lat + 0.1 * lng * lng + \

          0.1 * lng * lat + 0.1 * math.sqrt(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

　　

  • js代码
//定义一些常量

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 {*[]}

     */

     function bd09togcj02(bd_lon, bd_lat) {

     	var x_pi = 3.14159265358979324 * 3000.0 / 180.0;

     	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 {*[]}

     */

     function gcj02tobd09(lng, lat) {

     	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 {*[]}

     */

     function wgs84togcj02(lng, lat) {

     	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 gcj02towgs84

     * @param lng

     * @param lat

     * @returns {*[]}

     */

     function gcj02towgs84(lng, lat) {

     	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);

     		mglat = lat + dlat;

     		mglng = lng + dlng;

     		return [lng * 2 - mglng, lat * 2 - mglat]

     	}

     }

     function transformlat(lng, lat) {

     	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

     }

     function transformlng(lng, lat) {

     	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}

     */

     function out_of_china(lng, lat) {

     	return (lng < 72.004 || lng > 137.8347) || ((lat < 0.8293 || lat > 55.8271) || false)

     }

　　

python、js实现WGS84、高德（火星）、百度坐标转换的更多相关文章

  1. 使用Excel VBA编程将网点的百度坐标转换后标注到高德地图上

    公司网点表存储的坐标是百度坐标,现需要将网点位置标注到高德地图上,研究了一下高德地图的云图数据模版 http://lbs.amap.com/yuntu/reference/cloudstorage和坐 ...

  2. GPS坐标与百度坐标转换

    百度对外接口的坐标系,都是经过国家测绘局加密处理,符合国家测绘局对地理信息保密要求. 国际经纬度坐标标准为WGS-84,国内必须至少使用国测局制定的GCJ- 02,对地理位置进行首次加密.百度坐标在此 ...

  3. PHP+百度地图API+JAVASCRIPT实现GPS坐标与百度坐标转换的实例

    原文:PHP+百度地图API+JAVASCRIPT实现GPS坐标与百度坐标转换的实例 <!--小幅的坐标转换点位程序--> <!DOCTYPE html> <html&g ...

  4. Atitit.异步编程 java .net php python js 对照

    Atitit.异步编程 java .net php python js 的比較 1. 1.异步任务,异步模式,  APM模式,,  EAP模式, TAP 1 1.1.       APM模式: Beg ...

  5. js实现谷歌坐标转百度坐标

    js实现谷歌坐标转百度坐标 谷歌坐标转百度坐标 实现算法如下(以js为例,其他语言调整就行): //$lat 维度:$lng 经度  function GCJTobaidu($lat, $lng){ ...

  6. Atitit.异步编程 java .net php python js 的比较

    Atitit.异步编程 java .net php python js 的比较 1. 1.异步任务,异步模式,  APM模式,,  EAP模式, TAP 1 1.1.       APM模式: Beg ...

  7. PHP腾讯与百度坐标转换

    function coordinate_switch($a,$b){//百度转腾讯坐标转换 $a = Latitude , $b = Longitude $x = (double)$b - 0.006 ...

  8. Python爬虫实战之爬取百度贴吧帖子

    大家好,上次我们实验了爬取了糗事百科的段子,那么这次我们来尝试一下爬取百度贴吧的帖子.与上一篇不同的是,这次我们需要用到文件的相关操作. 本篇目标 对百度贴吧的任意帖子进行抓取 指定是否只抓取楼主发帖 ...

  9. 《python开发技术详解》|百度网盘免费下载|Python开发入门篇

    <python开发技术详解>|百度网盘免费下载|Python开发入门篇 提取码:2sby  内容简介 Python是目前最流行的动态脚本语言之一.本书共27章,由浅入深.全面系统地介绍了利 ...

随机推荐

  1. Js script type="text/template"的使用简单说明

    <script type="text/template" id="treeTableTpl"> <tr id="{{row.id}} ...

  2. csrf跨站点请求伪造

    什么是csrf(跨站请求伪造) 伪造请求的定义有很多种,我将不是用户本意发出的请求统称为伪造请求(在用户不知情的情况下执行某些操作)xss的通过用户对浏览器的信任造成的,csrf是通过服务器对浏览器的 ...

  3. 实现简单的 JS 模块加载器

    实现简单的 JS 模块加载器 1. 背景介绍 按需加载是前端性能优化的一个重要手段,按需加载的本质是从远程服务器加载一段JS代码(这里主要讨论JS,CSS或者其他资源大同小异),该JS代码就是一个模块 ...

  4. jsonp实现js跨域请求

    sonp是跨域通信的一个协议 具体来说jsonp实现跨域请求其实是使用js文件引用(js文件不一定是.js结尾)可跨域的性质,将请求的结果包裹在客户端需要调用的js方法内部.需要前后端配合使用. 前段 ...

  5. MySQL数据库、数据表和字段字符集查询、修改和配置

    一.设置编码 LINUX  修改vi/etc/my.cnf WINDOWS my.ini 在[client]下添加    default-character-set=utf8 在[mysqld]下添加 ...

  6. 2020 CCPC比赛

    https://ac.nowcoder.com/acm/contest/4010/A 这道题枚举区间长度的贡献值: 当区间长度为1时,就是所有元音数的个数: 当区间长度为2时,区间[2,n-1]的数贡 ...

  7. PHP SDK+Oss 上传文件流

    // Endpoint以杭州为例,其它Region请按实际情况填写. $endpoint = "http://oss-cn-hangzhou.aliyuncs.com"; // 云 ...

  8. python-excel读取-pyodbc

    https://github.com/mkleehammer/pyodbc/wiki/Cursor 利用pyodbc读取数据库,流程基本一样,就是配置connect对象时有所不同,下面是excel的: ...

  9. c#常用控件及简写

    c#常用控件及简写

  10. Java实现Word/Pdf/TXT转html

    引言: 最近公司在做一个教育培训学习及在线考试的项目,本人主要从事网络课程模块,主要做课程分类,课程,课件的创建及在线学习和统计的功能,因为课件涉及到多种类型,像视频,音频,图文,外部链接及文档类型. ...