[转]百度坐标（BD09）、国测局坐标（火星坐标，GCJ02）、和WGS84坐标系之间的转换
来源：博客园
注：项目中面对不同的坐标体系，在地图上显示会有不小的偏差，下面是JavaScript的转换方法。 //定义一些常量 var x_PI = 3.79324 * 3000.0 / 180.0; var PI = 3.; var a = ; var ee = 0.; /** * 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换 * 即 百度 转 谷歌、高德 * @param bd_lon * @param bd_lat * @returns {*[]} */ function bd09togcj02 (bd_lon, bd_lat) { var x_pi = 3.79324 * 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 * 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 + var mglng = lng + return [mglng, mglat] } } /** * GCJ02 转换为 WGS84 * @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 * 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 + mglng = lng + 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
137.8347) || ((lat
55.8271) || false); } //使用示例 //gps坐标转火星坐标 var lng_lat_1 = wgs84togcj02(113..497629); console.log('火星坐标...',lng_lat_1); //火星坐标转百度坐标 var lng_lat_2 = gcj02tobd09(lng_lat_1[0], lng_lat_1[1]); console.log('百度坐标...',lng_lat_2); 附上在线地图经纬度查询链接： /maps.html
免责声明：本站部分内容、图片、文字、视频等来自于互联网，仅供大家学习与交流。相关内容如涉嫌侵犯您的知识产权或其他合法权益，请向本站发送有效通知，我们会及时处理。反馈邮箱&&&&。
学生服务号
在线咨询，奖学金返现，名师点评，等你来互动温馨提示！由于新浪微博认证机制调整，您的新浪微博帐号绑定已过期，请重新绑定！&&|&&
LOFTER精选
网易考拉推荐
用微信&&“扫一扫”
将文章分享到朋友圈。
用易信&&“扫一扫”
将文章分享到朋友圈。
参考:http://blog.csdn.net/coolypf/article/details/8686588http://blog.csdn.net/giswens/article/details/8775283/SourceControl/changeset/view/关于 GCJ-02 和 BD-09 ，请参考&&。算法代码如下，其中&bd_encrypt&将 GCJ-02 坐标转换成 BD-09 坐标，&bd_decrypt&反之。[cpp]&
http://blog.csdn.net/giswens/article/details/8775275其他参考/code/file/3872.htm?dir=/correct_lat_lon/src/com/hgq/correct/MapFix.java
阅读(1539)|
用微信&&“扫一扫”
将文章分享到朋友圈。
用易信&&“扫一扫”
将文章分享到朋友圈。
历史上的今天
loftPermalink:'',
id:'fks_',
blogTitle:'地球坐标系 (WGS-84) 到火星坐标系 (GCJ-02)百度坐标系 (BD-09) 的转换算法',
blogAbstract:'转自http://blog.csdn.net/winnyrain/article/details/地球坐标系 (WGS-84) 到火星坐标系 (GCJ-02) 的转换算法
{if x.moveFrom=='wap'}
{elseif x.moveFrom=='iphone'}
{elseif x.moveFrom=='android'}
{elseif x.moveFrom=='mobile'}
${a.selfIntro|escape}{if great260}${suplement}{/if}
{list a as x}
推荐过这篇日志的人：
{list a as x}
{if !!b&&b.length>0}
他们还推荐了：
{list b as y}
转载记录：
{list d as x}
{list a as x}
{list a as x}
{list a as x}
{list a as x}
{if x_index>4}{break}{/if}
${fn2(x.publishTime,'yyyy-MM-dd HH:mm:ss')}
{list a as x}
{if !!(blogDetail.preBlogPermalink)}
{if !!(blogDetail.nextBlogPermalink)}
{list a as x}
{if defined('newslist')&&newslist.length>0}
{list newslist as x}
{if x_index>7}{break}{/if}
{list a as x}
{var first_option =}
{list x.voteDetailList as voteToOption}
{if voteToOption==1}
{if first_option==false},{/if}&&“${b[voteToOption_index]}”&&
{if (x.role!="-1") },“我是${c[x.role]}”&&{/if}
&&&&&&&&${fn1(x.voteTime)}
{if x.userName==''}{/if}
网易公司版权所有&&
{list x.l as y}
{if defined('wl')}
{list wl as x}{/list}Javascript 是一种由Netscape的LiveScript发展而来的原型化继承的基于对象的动态类型的区分大小写的客户端脚本语言，主要目的是为了解决服务器端语言，比如Perl，遗留的速度问题，为客户提供更流畅的浏览效果。
这篇文章主要介绍了百度坐标（BD09）、国测局坐标（火星坐标，GCJ02）、和WGS84坐标系之间的转换的相关资料,需要的朋友可以参考下
在项目中面对不同的坐标体系，在地图上显示多多少少都会有点偏差，下面是使用javascript写的转换方法，具体代码如下所示：
　//定义一些常量
　　var x_PI = 3.79324 * 3000.0 / 180.0;
　　var PI = 3.;
　　var a = ;
　　var ee = 0.;
　　* 百度坐标系 (BD-09) 与 火星坐标系 (GCJ-02)的转换
　　* 即 百度 转 谷歌、高德
　　* @param bd_lon
　　* @param bd_lat
　　* @returns {*[]}
　　function bd09togcj02(bd_lon, bd_lat) {
　　var x_pi = 3.79324 * 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]
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 *
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 +
var mglng = lng +
return [mglng, mglat]
* GCJ02 转换为 WGS84
* @param lng
* @param lat
* @returns {*[]}
function gcj02towgs84(lng, lat) {
if (out_of_china(lng, lat)) {
return [lng, lat]
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 *
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 +
mglng = lng +
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);
　　//使用示例
　　//gps坐标转火星坐标
var lng_lat_1 = wgs84togcj02(113..497629);
console.log('火星坐标...',lng_lat_1);
//火星坐标转百度坐标
var lng_lat_2 = gcj02tobd09(lng_lat_1[0], lng_lat_1[1]);
console.log('百度坐标...',lng_lat_2);
Copyright &
All Rights Reserved &&&&&&SupportBy摘要：各种坐标体系之间如何转换？到底有哪些坐标体系？什么是火星坐标？为什么我的坐标，在地图上显示会有偏移？本文详细解答以上问题。最后给出坐标拾取工具。
---------------------------------------------------------------------------------------------------------------------------
一、坐标体系
首先我们要明白，开发者能接触到哪些坐标体系呢？
第一种分类：
1、& GPS，WGS-84，原始坐标体系。一般用国际标准的GPS记录仪记录下来的坐标，都是GPS的坐标。很可惜，在中国，任何一个地图产品都不允许使用GPS坐标，据说是为了保密。GPS坐标形式如图，度分秒形式的经纬度：
2、& GCJ-02，国测局02年发布的坐标体系。又称&火星坐标&。在中国，必须至少使用GCJ-02的坐标体系。比如谷歌，腾讯，高德都在用这个坐标体系。GCJ-02也是国内最广泛使用的坐标体系。
3、& 其他坐标体系。一般都是由GCJ-02进过偏移算法得到的。这种体系就根据每个公司的不同，坐标体系都不一样了。比如，百度和搜狗就使用自己的坐标体系，与其他坐标体系不兼容。
第二种分类：
首先明白，所有坐标体系的原点，都是非洲。
1、& 经纬度。这个是球面坐标，对于北京来说，就是(116..935961)这样的坐标。比如腾讯、高德、百度都是这样的经纬度坐标。谷歌是经纬度顺序写反的经纬度坐标。
如果是度分秒坐标，需要进行转换，才能得到这样的经纬度坐标。详见坐标转换。
2、& 墨卡托坐标。平面坐标，相当于是直线距离，数字一般都比较大，像这样的。(526.)
墨卡托坐标，主要用于程序的后台计算。直线距离嘛，加加减减几乎计算方便。
搜狗地图API就是直接使用的墨卡托坐标。
二、坐标转换
在各种web端平台，或者高德、腾讯、百度上取到的坐标，都不是GPS坐标，都是GCJ-02坐标，或者自己的偏移坐标系。
比如，你在谷歌地图API，高德地图API，腾讯地图API上取到的，都是GCJ-02坐标，他们三家都是通用的，也适用于大部分地图API产品，以及他们的地图产品。
例外，百度API上取到的，是BD-09坐标，只适用于百度地图相关产品。
例外，搜狗API上取到的，是搜狗坐标，只适用于搜狗地图相关产品。
例外，谷歌地球，google earth上取到的，是GPS坐标，而且是度分秒形式的经纬度坐标。在国内不允许使用。必须转换为GCJ-02坐标。
1、度分秒坐标转换为经纬度
比如，在GPS记录仪，或者google earth上采集到的是39&31'20.51，那么应该这样换算，31分就是31/60度，20.51秒就是20.51/3600度，结果就是39 + 31/60 + 20.51/3600 度。
2、& GPS转换为GCJ-02坐标
谷歌，高德，腾讯的地图API官网上，都不直接提供这样的坐标转换。如果要得到GCJ-02坐标，最好在他们的地图上直接取点，或者通过地址解析得到。（这个工具我后续会贴出来的。我就爱干这样的事情，哈哈。）
不过，在网上搜到了这样的接口，该接口的type=1就是GPS转到GCJ-02的墨卡托坐标。请大家对接口保密，哈哈。详见：
3、GCJ-02与BD-09之间互转
国测局GCJ-02坐标体系（谷歌、高德、腾讯），与百度坐标BD-09体系的转换，在CSDN上有很详细的讲解：
不过也有更简单的算法，线性算法（lat和lng是经纬度，球面坐标）：
To_B是转到百度，To_G是转到GCJ-02。
var TO_BLNG = function(lng){return
lng+0.0065;};
var TO_BLAT = function(lat){return
lat+0.0060;};
var TO_GLNG = function(lng){return
lng-0.0065;};
var TO_GLAT = function(lat){return
lat-0.0060;};
4、经纬纬度转成墨卡托
网上也有详细讲解：
5、各家API公司坐标转换接口的申请
一般需要将您的公司名称、项目名称、项目简介、联系人和联系方式，发邮件至地图API公司的商务部，经过申请，才能使用。
下面是他们的联系方式：
这个，国内没有公交接口，大家应该不会再使用了吧。　　
（大家发现没，高德是api，腾讯和百度是mapapi，说明什么？）
三、坐标偏移
如果您的坐标在转换之后，还有偏移，那么考虑以下几个方面。
A、原始坐标系弄错，比如以为自己是GPS坐标，但其实已经是GCJ-02坐标。
解决方案：请确保采集到的数据是哪个坐标体系，需要转换到哪个坐标系，再进行坐标转换。
B、原始坐标准确度不够
解决方案：如果您是GPS坐标，请确保采集GPS数据时，搜到至少4颗以上的卫星。并且GPS数据准不准，还取决于周围建筑物的高度，越高越不准，因为有遮挡。
如果本来就是GCJ-02坐标，在不同地图放大级别的时候，看到的地方可能不一样。比如你在地图级别4（国家）取到的坐标，放大到地图12级（街道）时，坐标就偏了。请确保在地图最大放大级别时，拾取坐标。
C、度分秒的概念混淆
比如，在google earth上采集到的是39&31'20.51，那么应该这样换算，31分就是31/60度，20.51秒就是20.51/3600度，结果就是39 + 31/60 + 20.51/3600 度。
D、经纬度顺序写反了
有些公司（比如高德，百度，腾讯）是先经度，再纬度，即Point(lng, lat)。但谷歌坐标的顺序恰好相反，是(lat, lng)。
四、坐标拾取工具
既然在国内必须至少使用GCJ-02的坐标系，而GCJ-02，&火星坐标&是在国内最广泛使用的坐标体系。那么，我们就来看看，如何直接获取到GCJ-02坐标呗。
请大家把这段代码保存到记事本里，然后后缀名改为.html，记得用UTF-8编码来保存。然后双击这个文件，就能打开网页了。
注册你自己的key，填在源代码里哦：&
&!DOCTYPE HTML&
&meta http-equiv="Content-Type" content="text/ charset=utf-8"&
&title&火星坐标拾取&/title&
&script language="javascript" src="/maps?v=1.2&key=您的key"&&/script&
#iMap{height:<span style="background-color: #f5f5f5; color: #px;width:<span style="background-color: #f5f5f5; color: #px;float:left;}
.info{float:left;margin:<span style="background-color: #f5f5f5; color: # 0 10px;}
label{width:80px;float:left;}
&body onLoad="mapInit()"&
&div id="iMap"&&/div&
&div class="info"&
&p&&label&火星坐标：&/label&&span id="lnglat"&&&/span&&/p&
&p&&label&地址：&/label&&span id="iAddress"&&&/span&&/p&
&p&说明：&/p&
&p&1、鼠标滚轮可以缩放地图，拖动地图。&/p&
&p&2、点击地图，即可获得GCJ-02的经纬度坐标，地址。&/p&
&script language="javascript"&
var lnglatXY;
//初始化地图
function mapInit(){
var opt = {
level: 13, //设置地图缩放级别
center: new AMap.LngLat(116.397428, 39.90923) //设置地图中心点
mapObj = new AMap.Map("iMap", opt);
AMap.event.addListener(mapObj,'click',getLnglat); //点击事件
function geocoder() {
//加载地理编码插件
mapObj.plugin(["AMap.Geocoder"], function() {
MGeocoder = new AMap.Geocoder({
radius: 1000,
extensions: "all"
//返回地理编码结果
AMap.event.addListener(MGeocoder, "complete", geocoder_CallBack);
//逆地理编码
MGeocoder.getAddress(lnglatXY);
var marker = new AMap.Marker({
map:mapObj,
icon: new AMap.Icon({
image: "/Public/images/js/mark.png",
size:new AMap.Size(58,30),
imageOffset: new AMap.Pixel(-32, -0)
position: lnglatXY,
offset: new AMap.Pixel(-5,-30)
mapObj.setFitView();
//回调函数
function geocoder_CallBack(data) {
//返回地址描述
address = data.regeocode.formattedA
//返回结果拼接输出
document.getElementById("iAddress").innerHTML =
//鼠标点击，获取经纬度坐标
function getLnglat(e){
mapObj.clearMap();
var x = e.lnglat.getLng();
var y = e.lnglat.getLat();
document.getElementById("lnglat").innerHTML = x + "," +
lnglatXY = new AMap.LngLat(x,y);
geocoder();
打开之后，是这样的：
&点击一下地图，即可获得GCJ-02坐标，地址：
demo地址：
高德官方坐标拾取工具：
阅读(...) 评论()