在日常工作学习中常会涉及到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. Python七夕记

  2. 三、统一威胁管理（UTM）

    简介 统一威胁管理(Unified Threat Management),简称UTM. 2004年9月,IDC首度提出“统一威胁管理”的概念,即将防病毒.入侵检测和防火墙安全设备划归统一威胁管理(Un ...

  3. java基础(三)之面向对象编程

    对象的创建方法 语法: class 类名{ 属性; 方法; } 生成对象的方法 类名 对象名 = new 类名(); Dog dog = new Dog(); 对象的使用方法1.对象.变量;2.对象. ...

  4. android 获取颜色 getColor 方法 deprecated 过期

    可以使用下面代码代替: ContextCompat.getColor(getContext(), R.color.post_list_content_color) 需要引入: compile 'com ...

  5. css div布局示例2（head-main-footer

    <!doctype html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  6. Mysql主键外键操作

    外键: ·一对多 ·多对多 ·一对一 ·修改表 ·复制表   主键:   rimary key auto_increment(primary key是主键通常和auto_increment自动增加混合 ...

  7. 第三十二篇 玩转数据结构——AVL树（AVL Tree）

          1.. 平衡二叉树 平衡二叉树要求,对于任意一个节点,左子树和右子树的高度差不能超过1. 平衡二叉树的高度和节点数量之间的关系也是O(logn) 为二叉树标注节点高度并计算平衡因子 AVL ...

  8. js解释器

    又名js引擎 JavaScript是解释型语言,这就是说它无需编译,直接由JavaScript引擎直接执行. 既然说到了解释型语言,那么我们就来分别以下解释型语言和编译型语言的差别: 编译型语言:程序 ...

  9. 问题 I: 数字分组２

    问题 I: 数字分组2 时间限制: 1 Sec  内存限制: 128 MB[命题人:admin] 题目描述 已知一堆魔法石的重量,问如何分成两堆,使得它们质量和之差最大,但不能大于(可以等于)这些数中 ...

  10. Bridge(Ad Hoc)