使用js对WGS-84 ,GCJ-02与BD-09的坐标进行转换

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

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)之间的相互转换的方法,经测试,是转换后相对准确可用的:

作者:青檬可乐_

链接:http://08643.cn/p/d3dd4149bb0b

來源:简书

著作权归作者所有。商业转载请联系作者获得授权,非商业转载请注明出处。

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-02??

gcj_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-84??

gcj_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?exactly??

gcj_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-09??

bd_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-02??

bd_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?->?x??

mercator_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?->?x??

mercator_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?distance??

distance: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;??

????}??

};??

最后编辑于
?著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 213,992评论 6 493
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,212评论 3 388
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 159,535评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,197评论 1 287
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,310评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,383评论 1 292
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,409评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,191评论 0 269
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,621评论 1 306
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,910评论 2 328
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,084评论 1 342
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,763评论 4 337
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,403评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,083评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,318评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 46,946评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 43,967评论 2 351