GPS坐标互轉：WGS-84(GPS)、GCJ-02(Google地圖)、BD-09(百度地圖)

WGS-84：是國際标準，GPS坐标（Google Earth使用、或者GPS子產品）

GCJ-02：中國坐标偏移标準，Google Map、高德、騰訊使用

BD-09：百度坐标偏移标準，Baidu Map使用

//WGS-84 to GCJ-02

GPS.gcj_encrypt();

//GCJ-02 to WGS-84 粗略

GPS.gcj_decrypt();

//GCJ-02 to WGS-84 精确(二分極限法)

// var threshold = 0.000000001; 目前設定的是精确到小數點後9位，這個值越小，越精确，但是javascript中，浮點運算本身就不太精确，九位在GPS裡也偏差不大了

GSP.gcj_decrypt_exact();

//GCJ-02 to BD-09

GPS.bd_encrypt();

//BD-09 to GCJ-02

GPS.bd_decrypt();

//求距離

GPS.distance();

示例：

document.write("GPS: 39.933676862706776,116.35608315379092<br />");

var arr2 = GPS.gcj_encrypt(39.933676862706776, 116.35608315379092);

document.write("中國:" + arr2['lat']+","+arr2['lon']+'<br />');

var arr3 = GPS.gcj_decrypt_exact(arr2['lat'], arr2['lon']);

document.write('逆算:' + arr3['lat']+","+arr3['lon']+' 需要和第一行相似（目前是小數點後9位相等）');

JavaScript代碼

<code> </code><code>delta :</code><code>function</code> <code>(lat, lon) {</code>

<code> </code><code>// Krasovsky 1940</code>

<code> </code><code>var</code> <code>a = 6378245.0;</code><code>// a: 衛星橢球坐标投影到平面地圖坐标系的投影因子。</code>

<code> </code><code>var</code> <code>ee = 0.00669342162296594323;</code><code>// ee: 橢球的偏心率。</code>

<code> </code><code>var</code> <code>dLat =</code><code>this</code><code>.transformLat(lon - 105.0, lat - 35.0);</code>

<code> </code><code>var</code> <code>dLon =</code><code>this</code><code>.transformLon(lon - 105.0, lat - 35.0);</code>

<code> </code><code>var</code> <code>radLat = lat / 180.0 *</code><code>this</code><code>.PI;</code>

<code> </code><code>var</code> <code>magic = Math.sin(radLat);</code>

<code> </code><code>magic = 1 - ee * magic * magic;</code>

<code> </code><code>var</code> <code>sqrtMagic = Math.sqrt(magic);</code>

<code> </code><code>dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) *</code><code>this</code><code>.PI);</code>

<code> </code><code>dLon = (dLon * 180.0) / (a / sqrtMagic * Math.cos(radLat) *</code><code>this</code><code>.PI);</code>

<code> </code><code>return</code> <code>{</code><code>'lat'</code><code>: dLat,</code><code>'lon'</code><code>: dLon};</code>

<code> </code><code>gcj_encrypt :</code><code>function</code> <code>(wgsLat, wgsLon) {</code>

<code> </code><code>if</code> <code>(</code><code>this</code><code>.outOfChina(wgsLat, wgsLon))</code>

<code> </code><code>return</code> <code>{</code><code>'lat'</code><code>: wgsLat,</code><code>'lon'</code><code>: wgsLon};</code>

<code> </code><code>var</code> <code>d =</code><code>this</code><code>.delta(wgsLat, wgsLon);</code>

<code> </code><code>return</code> <code>{</code><code>'lat'</code> <code>: wgsLat + d.lat,</code><code>'lon'</code> <code>: wgsLon + d.lon};</code>

<code> </code><code>gcj_decrypt :</code><code>function</code> <code>(gcjLat, gcjLon) {</code>

<code> </code><code>if</code> <code>(</code><code>this</code><code>.outOfChina(gcjLat, gcjLon))</code>

<code> </code><code>return</code> <code>{</code><code>'lat'</code><code>: gcjLat,</code><code>'lon'</code><code>: gcjLon};</code>

<code> </code><code>var</code> <code>d =</code><code>this</code><code>.delta(gcjLat, gcjLon);</code>

<code> </code><code>return</code> <code>{</code><code>'lat'</code><code>: gcjLat - d.lat,</code><code>'lon'</code><code>: gcjLon - d.lon};</code>

<code> </code><code>//GCJ-02 to WGS-84 exactly</code>

<code> </code><code>gcj_decrypt_exact :</code><code>function</code> <code>(gcjLat, gcjLon) {</code>

<code> </code><code>var</code> <code>initDelta = 0.01;</code>

<code> </code><code>var</code> <code>threshold = 0.000000001;</code>

<code> </code><code>var</code> <code>dLat = initDelta, dLon = initDelta;</code>

<code> </code><code>var</code> <code>mLat = gcjLat - dLat, mLon = gcjLon - dLon;</code>

<code> </code><code>var</code> <code>pLat = gcjLat + dLat, pLon = gcjLon + dLon;</code>

<code> </code><code>var</code> <code>wgsLat, wgsLon, i = 0;</code>

<code> </code><code>while</code> <code>(1) {</code>

<code> </code><code>wgsLat = (mLat + pLat) / 2;</code>

<code> </code><code>wgsLon = (mLon + pLon) / 2;</code>

<code> </code><code>var</code> <code>tmp =</code><code>this</code><code>.gcj_encrypt(wgsLat, wgsLon)</code>

<code> </code><code>dLat = tmp.lat - gcjLat;</code>

<code> </code><code>dLon = tmp.lon - gcjLon;</code>

<code> </code><code>break</code><code>;</code>

<code> </code><code>if</code> <code>(dLat > 0) pLat = wgsLat;</code><code>else</code> <code>mLat = wgsLat;</code>

<code> </code><code>if</code> <code>(dLon > 0) pLon = wgsLon;</code><code>else</code> <code>mLon = wgsLon;</code>

<code> </code><code>if</code> <code>(++i > 10000)</code><code>break</code><code>;</code>

<code> </code><code>//console.log(i);</code>

<code> </code><code>return</code> <code>{</code><code>'lat'</code><code>: wgsLat,</code><code>'lon'</code><code>: wgsLon};</code>

<code> </code><code>bd_encrypt :</code><code>function</code> <code>(gcjLat, gcjLon) {</code>

<code> </code><code>var</code> <code>x = gcjLon, y = gcjLat; </code>

<code> </code><code>var</code> <code>theta = Math.atan2(y, x) + 0.000003 * Math.cos(x *</code><code>this</code><code>.x_pi); </code>

<code> </code><code>bdLon = z * Math.cos(theta) + 0.0065; </code>

<code> </code><code>bdLat = z * Math.sin(theta) + 0.006;</code>

<code> </code><code>return</code> <code>{</code><code>'lat'</code> <code>: bdLat,</code><code>'lon'</code> <code>: bdLon};</code>

<code> </code><code>bd_decrypt :</code><code>function</code> <code>(bdLat, bdLon) {</code>

<code> </code><code>var</code> <code>x = bdLon - 0.0065, y = bdLat - 0.006; </code>

<code> </code><code>var</code> <code>theta = Math.atan2(y, x) - 0.000003 * Math.cos(x *</code><code>this</code><code>.x_pi); </code>

<code> </code><code>var</code> <code>gcjLon = z * Math.cos(theta); </code>

<code> </code><code>var</code> <code>gcjLat = z * Math.sin(theta);</code>

<code> </code><code>return</code> <code>{</code><code>'lat'</code> <code>: gcjLat,</code><code>'lon'</code> <code>: gcjLon};</code>

<code> </code><code>//WGS-84 to Web mercator</code>

<code> </code><code>//mercatorLat -> y mercatorLon -> x</code>

<code> </code><code>mercator_encrypt :</code><code>function</code><code>(wgsLat, wgsLon) {</code>

<code> </code><code>var</code> <code>x = wgsLon * 20037508.34 / 180.;</code>

<code> </code><code>var</code> <code>y = Math.log(Math.tan((90. + wgsLat) *</code><code>this</code><code>.PI / 360.)) / (</code><code>this</code><code>.PI / 180.);</code>

<code> </code><code>return</code> <code>{</code><code>'lat'</code> <code>: y,</code><code>'lon'</code> <code>: x};</code>

<code> </code><code>if ((Math.abs(wgsLon) > 180 || Math.abs(wgsLat) > 90))</code>

<code> </code><code>return null;</code>

<code> </code><code>var x = 6378137.0 * wgsLon * 0.017453292519943295;</code>

<code> </code><code>var a = wgsLat * 0.017453292519943295;</code>

<code> </code><code>return {'lat' : y, 'lon' : x};</code>

<code> </code><code>// Web mercator to WGS-84</code>

<code> </code><code>// mercatorLat -> y mercatorLon -> x</code>

<code> </code><code>mercator_decrypt :</code><code>function</code><code>(mercatorLat, mercatorLon) {</code>

<code> </code><code>var</code> <code>x = mercatorLon / 20037508.34 * 180.;</code>

<code> </code><code>var</code> <code>y = mercatorLat / 20037508.34 * 180.;</code>

<code> </code><code>if (Math.abs(mercatorLon) < 180 && Math.abs(mercatorLat) < 90)</code>

<code> </code><code>if ((Math.abs(mercatorLon) > 20037508.3427892) || (Math.abs(mercatorLat) > 20037508.3427892))</code>

<code> </code><code>var a = mercatorLon / 6378137.0 * 57.295779513082323;</code>

<code> </code><code>var x = a - (Math.floor(((a + 180.0) / 360.0)) * 360.0);</code>

<code> </code><code>var y = (1.5707963267948966 - (2.0 * Math.atan(Math.exp((-1.0 * mercatorLat) / 6378137.0)))) * 57.295779513082323;</code>

<code> </code><code>// two point's distance</code>

<code> </code><code>distance :</code><code>function</code> <code>(latA, lonA, latB, lonB) {</code>

<code> </code><code>var</code> <code>earthR = 6371000.;</code>

<code> </code><code>var</code> <code>alpha = Math.acos(s);</code>

<code> </code><code>var</code> <code>distance = alpha * earthR;</code>

<code> </code><code>return</code> <code>distance;</code>

<code> </code><code>outOfChina :</code><code>function</code> <code>(lat, lon) {</code>

<code> </code><code>return</code> <code>true</code><code>;</code>

<code> </code><code>return</code> <code>false</code><code>;</code>

<code> </code><code>transformLat :</code><code>function</code> <code>(x, y) {</code>

<code> </code><code>return</code> <code>ret;</code>

<code> </code><code>transformLon :</code><code>function</code> <code>(x, y) {</code>

若轉載請注明出處！若有疑問，請回複交流！

GPS坐标互轉：WGS-84(GPS)、GCJ-02(Google地圖)、BD-09(百度地圖)

繼續閱讀

賽維闆報丨動态相對GPS定位技術其他的GPS作業模式觀測資料需在測後處理，不僅無法實時地給出觀測站的定位結果，而且也無法

#前端#程式員#前端開發#前端面試自學過程中，除了學技術本身，更重要的是理清技術脈絡，用技術脈絡将技術串聯起來，形成系統

#哈羅單車我真的是服了，從頭到尾沒有一句說是自己定位問題，我拍的照片留的證據都不行，一直用哈羅，因為是支付寶裡面的小程式

不同場景UWB-TAG标簽方案的靈活應用，飛睿智能室内定位技術方案

賽維闆報丨GPS定位技術在地籍圖測繪中的應用同全站儀相比，RTK測量單點的時間需要幾秒到幾十秒，但是，它不要求通視，不需

「真的很少有老師，能就一個定位問題講上2個小時，老師實在是講得太全面周到了」「經過老師的剖析，我才發現自己的優勢所在」有

店鋪異常預警功能：電商物流/發貨/效期/單量預警(一)。大家好，今天來講解一下異常預警功能。現在淘寶、拼多多、抖音等各大

HarmonyOS在SDK9版本下FA模型geolocation無法定位問題解決

金航标kinghelm銷售總監鄧海峰介紹說，藍牙技術作為一種低成本的近距離無線連接配接技術，為固定和移動裝置建立了良好的通信

如何使用P-One的全鍊路壓測工具來定位問題節點和分析性能問題

室内定位技術的發展演變

IP地理位置定位技術原理

IP位址定位技術在各業中的應用

常用的北鬥高精度定位技術有哪些？

智能解決方案！無線傳感器網絡中的混合彙聚節點重新定位技術

#緬北#在某合網上，碰到了緬北做勒索詐騙的，用備用手機戲耍的騙子，手機裡放了30多個空号，分别備注上“總經理趙XX”、“