GPS座標轉百度地圖座標的方法

首先需要認識一下GPS的座標系。GPS座標系遵循WGS-84標準，在這個標準下，GPS芯片可以發出不同的數據包格式。根據其數據幀幀頭的不同，GPS數據可以分類爲GPGGA、GPGSA、GPGSV、GPRMC等。這些幀頭標識了後續幀內數據的組成結構。通常情況下，我們所關心的定位數據如經緯度、速度、時間等均可以從GPRMC幀中獲取得到。

在次我不講解具體的幀的格式，在網上可以很容易搜到，也可以總芯片配套的接口文檔找到數據幀的格式。通過解讀GPRMC包，我們可以得到類似於“3040.8639,N,10405.7573,E”的位置數據，其字段含義可以在文檔中找到。這裏我解析的結果是：北緯30°40.8639'，東經104°5.7573'。

在得到這兩個數據後，我們可以試着去在地圖上定位一下，我們發現多數地圖定位用的是小數形式的座標。這裏我們就需要去講如何把一個GPS的原始數據，轉換爲小數形式的數據。

舉例：106°14'15"的轉換
因爲度分秒都是六十進制的
所以可以這樣轉換：
15/60=0.25分
(14+0.25)/60=0.2375度
106+0.2375=106.2375度

所以最後的結果是106.2375°

通過這個轉換之後，我的位置座標【北緯30°40.8639'，東經104°5.7573'】就可以變成【30.681065N,104.095955E】（這個就是WGS-84的座標）這個模樣。這個座標是GPS的物理定位，根據國際標準，需要對這個座標進行GCJ-02偏移轉換，轉換後的座標纔可以在google地圖、高德地圖、以及騰訊地圖上定位（以上三家遵循GCJ-02加密）。至於爲什麼要這樣做，是因爲加密座標的原因。因爲GCJ-02是不可逆的轉換。

這裏貼出由準尋WGS-84標準轉換爲GCJ-02的C++源碼：

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;
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 * sqrt(abs(x));
ret += (20.0 * sin(6.0 * x * pi) + 20.0 * sin(2.0 * x * pi)) * 2.0 / 3.0; 
ret += (20.0 * sin(y * pi) + 40.0 * sin(y / 3.0 * pi)) * 2.0 / 3.0; 
ret += (160.0 * sin(y / 12.0 * pi) + 320 * 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 * sqrt(abs(x)); 
ret += (20.0 * sin(6.0 * x * pi) + 20.0 * sin(2.0 * x * pi)) * 2.0 / 3.0; 
ret += (20.0 * sin(x * pi) + 40.0 * sin(x / 3.0 * pi)) * 2.0 / 3.0;
ret += (150.0 * sin(x / 12.0 * pi) + 300.0 * sin(x / 30.0 * pi)) * 2.0 / 3.0; 
return ret;
} 
/*
參數
wgLat:WGS-84緯度wgLon:WGS-84經度
返回值：
mgLat：GCJ-02緯度mgLon：GCJ-02經度
*/
void gps_transform( double wgLat, double wgLon, double& mgLat, 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 = sin(radLat);
magic = 1 - ee * magic * magic; double sqrtMagic = sqrt(magic);
dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);
dLon = (dLon * 180.0) / (a / sqrtMagic * cos(radLat) * pi);
mgLat = wgLat + dLat; mgLon = wgLon + dLon;
}

經過以上的轉換，我們的座標就可以在遵循GCJ-02標準的地圖刪定位了，這裏給出大家一個比較好用的測試定位的網址：http://www.gpsspg.com/maps.htm

百度地圖有些特殊，它在GCJ-02基礎上又進行了一次加密，百度把這個加密標準叫做BD-09，但是這個加密並沒有對外公開。但是萬能的網絡以及無所不能的網友總會有辦法解決這類問題。這裏我就貼出一個可行的轉換辦法，經過測試，定位很準確。

#include <math.h>    
const double x_pi = 3.14159265358979324 * 3000.0 / 180.0;  
//將 GCJ-02 座標轉換成 BD-09 座標
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);  
} 

這樣，我們的座標也可以在百度地圖上準確定位了，滿心歡喜~~
 另外再貼上一個比較容易解釋座標爲何需要轉換的圖例，以方便大家理解：
    國際通用的轉換做法：

    百度的做法：

    這也解釋了，爲什麼GPS座標需要經過轉換的原因。感謝貢獻這些資料到網上的人們！
參考：
http://m.blog.csdn.net/blog/wildboy2001/39497681
http://www.haodaima.net/art/2441684
http://www.gpsspg.com/maps.htm
http://wenda.haosou.com/q/1365472754066079
http://www.geekcome.com/content-10-1464-1.html