地圖座標詳解及轉換方法

你們都知道，美國GPS使用的是WGS84的座標系統，以經緯度的形式來表示地球平面上的某一個位置，這應該是國際共識。但在我國，出於國家安全考慮，國內全部導航電子地圖必須使用國家測繪局制定的加密座標系統，即將一個真實的經緯度座標加密成一個不正確的經緯度座標，咱們在業內將前者稱之爲地球座標，後者稱之爲火星座標

下面是百度官方對百度座標爲什麼有偏移的解釋

　　國際經緯度座標標準爲WGS-84,國內必須至少使用國測局制定的GCJ-02,對地理位置進行首次加密。百度座標在此基礎上，進行了BD-09二次加密措施,更加保護了我的隱私。百度對外接口的座標系並非GPS採集的真實經緯度，須要經過座標轉換接口進行轉換。

1、座標體系

首先咱們要明白，開發者能接觸到哪些座標體系呢？

第一種分類：

一、  GPS，WGS-84，原始座標體系。

通常用國際標準的 GPS記錄儀 記錄下來的座標，都是GPS的座標。很惋惜，在中國，任何一個地圖產品都不容許使用GPS座標，聽說是爲了保密。GPS座標形式如圖，度分秒形式的經緯度：

二、  GCJ-02，國測局02年發佈的座標體系。又稱「火星座標」。在中國，必須至少使用GCJ-02的座標體系。好比谷歌，騰訊，高德都在用這個座標體系。GCJ-02也是國內最普遍使用的座標體系。

三、  其餘座標體系。通常都是由GCJ-02進過偏移算法獲得的。這種體系就根據每一個公司的不一樣，座標體系都不同了。好比，百度和搜狗就使用本身的座標體系，與其餘座標體系不兼容。

第二種分類：

首先明白，全部座標體系的原點，都是非洲。

一、  經緯度。這個是球面座標，對於北京來講，就是(116.38817139.935961)這樣的座標。

        好比騰訊、高德、百度都是這樣的經緯度座標。谷歌是經緯度順序寫反的經緯度座標。

        若是是度分秒座標，須要進行轉換，才能獲得這樣的經緯度座標。詳見座標轉換。

二、  墨卡託座標。平面座標，至關因而直線距離，數字通常都比較大，像這樣的。(215362.00021333335 99526.00034912192)

墨卡託座標，主要用於程序的後臺計算。直線距離嘛，加加減減幾乎計算方便。

搜狗地圖API就是直接使用的墨卡託座標。

2、座標轉換

在各類web端平臺，或者高德、騰訊、百度上取到的座標，都不是GPS座標，都是GCJ-02座標，或者本身的偏移座標系。

好比，你在谷歌地圖API，高德地圖API，騰訊地圖API上取到的，都是GCJ-02座標，他們三家都是通用的，也適用於大部分地圖API產品，以及他們的地圖產品。

例外，百度API上取到的，是BD-09座標，只適用於百度地圖相關產品。

例外，搜狗API上取到的，是搜狗座標，只適用於搜狗地圖相關產品。

例外，谷歌地球，google earth上取到的，是GPS座標，並且是度分秒形式的經緯度座標。在國內不容許使用。必須轉換爲GCJ-02座標。

一、度分秒座標轉換爲經緯度

好比，在GPS記錄儀，或者google earth上採集到的是39°31'20.51，那麼應該這樣換算，31分就是31/60度，20.51秒就是20.51/3600度，結果就是39 + 31/60 + 20.51/3600 度。

二、  GPS轉換爲GCJ-02座標

谷歌，高德，騰訊的地圖API官網上，都不直接提供這樣的座標轉換。若是要獲得GCJ-02座標，最好在他們的地圖上直接取點，或者經過地址解析獲得。（這個工具我後續會貼出來的。我就愛幹這樣的事情，哈哈。）

不過，在網上搜到了這樣的接口，該接口的type=1就是GPS轉到GCJ-02的墨卡託座標。請你們對接口保密，哈哈。詳見：

http://map.sogou.com/api/documentation/javascript/api2.5/interface_translate.html#late_intro

三、GCJ-02與BD-09之間互轉

國測局GCJ-02座標體系（谷歌、高德、騰訊），與百度座標BD-09體系的轉換，在CSDN上有很詳細的講解：

關於 GCJ-02 和 BD-09 ，請參考 http://developer.baidu.com/map/question.htm#qa0043 。

算法代碼以下，其中 bd_encrypt 將 GCJ-02 座標轉換成 BD-09 座標， bd_decrypt 反之

const double x_pi = 3.14159265358979324 * 3000.0 / 180.0;  
  
void bd_encrypt(double gg_lat, double gg_lon, double &bd_lat, double &bd_lon)  
{  
    double x = gg_lon, y = gg_lat;  
    double z = sqrt(x * x + y * y) + 0.00002 * sin(y * x_pi);  
    double theta = atan2(y, x) + 0.000003 * cos(x * x_pi);  
    bd_lon = z * cos(theta) + 0.0065;  
    bd_lat = z * sin(theta) + 0.006;  
}  
  
void bd_decrypt(double bd_lat, double bd_lon, double &gg_lat, double &gg_lon)  
{  
    double x = bd_lon - 0.0065, y = bd_lat - 0.006;  
    double z = sqrt(x * x + y * y) - 0.00002 * sin(y * x_pi);  
    double theta = atan2(y, x) - 0.000003 * cos(x * x_pi);  
    gg_lon = z * cos(theta);  
    gg_lat = z * sin(theta);  
}

不過也有更簡單的算法，線性算法（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.地球座標系 (WGS-84) 到火星座標系 (GCJ-02) 的轉換算法

WGS-84 到 GCJ-02 的轉換（即 GPS 加偏）算法是一個普通青年輕易沒法接觸到的「公開」的祕密。這個算法的代碼在互聯網上是公開的，詳情請使用 Google 搜索 "wgtochina_lb" 。

整理後的算法代碼請參考 https://on4wp7.codeplex.com/SourceControl/changeset/view/21483#353936 。代碼以下

using System;namespace Navi
{    class EvilTransform
    {        const double pi = 3.14159265358979324;        //
        // Krasovsky 1940
        //
        // a = 6378245.0, 1/f = 298.3
        // b = a * (1 - f)
        // ee = (a^2 - b^2) / a^2;
        const double a = 6378245.0;        const double ee = 0.00669342162296594323;        //
        // World Geodetic System ==> Mars Geodetic System
        public static void transform(double wgLat, double wgLon, out double mgLat, out double mgLon)
        {            if (outOfChina(wgLat, wgLon))
            {
                mgLat = wgLat;
                mgLon = wgLon;                return;
            }            double dLat = transformLat(wgLon - 105.0, wgLat - 35.0);            double dLon = transformLon(wgLon - 105.0, wgLat - 35.0);            double radLat = wgLat / 180.0 * pi;            double magic = Math.Sin(radLat);
            magic = 1 - ee * magic * magic;            double sqrtMagic = Math.Sqrt(magic);
            dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);
            dLon = (dLon * 180.0) / (a / sqrtMagic * Math.Cos(radLat) * pi);
            mgLat = wgLat + dLat;
            mgLon = wgLon + dLon;
        }        static bool outOfChina(double lat, double lon)
        {            if (lon < 72.004 || lon > 137.8347)                return true;            if (lat < 0.8293 || lat > 55.8271)                return true;            return false;
        }        static double transformLat(double x, double y)
        {            double 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 * pi) + 20.0 * Math.Sin(2.0 * x * pi)) * 2.0 / 3.0;
            ret += (20.0 * Math.Sin(y * pi) + 40.0 * Math.Sin(y / 3.0 * pi)) * 2.0 / 3.0;
            ret += (160.0 * Math.Sin(y / 12.0 * pi) + 320 * Math.Sin(y * pi / 30.0)) * 2.0 / 3.0;            return ret;
        }        static double transformLon(double x, double y)
        {            double 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 * pi) + 20.0 * Math.Sin(2.0 * x * pi)) * 2.0 / 3.0;
            ret += (20.0 * Math.Sin(x * pi) + 40.0 * Math.Sin(x / 3.0 * pi)) * 2.0 / 3.0;
            ret += (150.0 * Math.Sin(x / 12.0 * pi) + 300.0 * Math.Sin(x / 30.0 * pi)) * 2.0 / 3.0;            return ret;
        }
    }
}

知道了這個算法以後，就能夠離線進行 Google 地圖偏移校訂，沒必要像以前那麼麻煩。至於 GCJ-02 到 WGS-84 的轉換（即 GPS 糾偏），可使用二分法。

如何將非高德座標轉換爲高德座標系

Dear 開發者：

若是你苦悶於如何將非高德座標轉換爲高德座標系時，請看這裏。

當你在高德定位、規劃路徑等衆多優質服務面前徘徊，請必須看這裏！！

喜大普奔的消息來了！座標轉換服務在千呼萬喚中，終於bulingbuling 登場。

座標轉換後，使用JS API、SDK再無煩惱，大BOSS不再用擔憂個人工做，soeasy!

解決開發者的問題就是咱們的使命，一切開發者的障礙必須肅清。

座標轉換service:

一．適用對象：

苦苦糾結於如何將GPS、mapbar、baidu座標轉換爲高德座標系的開發者們。

二．轉換方法：

1.      首先固然仍是須要人手一個的百寶key。注意：權限君升級權限系統後，調用座標轉換接口的童鞋們須要從新申請。注意：綁定服務需選擇REST服務接口

2.      調用座標轉換接口，接口請求攻略以下：

1）  服務協議爲GET請求

2）  請求地址：http://restapi.amap.com/v3/assistant/coordinate/convert?parameters

3）  請求參數：

z轉換成功

將非百度座標轉換爲百度座標接口

http://developer.baidu.com/map/changeposition.htm

http://developer.baidu.com/map/question.htm#qa0043