地球座標，火星座標，百度座標轉換

         標準地球座標即GPS設備獲得的座標，該座標需要經過國家測繪局進行加密後形成火星座標（WGS-84 ），我們用的google座標和高德地圖座標也就是火星座標

         百度地圖，在火星座標的基礎上再進行一次加密，形成了百度地圖上的座標，因此，直接將標準地球座標顯示在百度地圖上是會有幾百米的偏差的。按照此原理，標準GPS座標經過兩步的轉換可得到百度座標。因爲在處理百度地圖時，例如查詢其POI都需要百度地圖上的座標，而不是標準座標，那麼這樣的轉換就是必要的了。下面是兩步轉換的C++ 程序，供開發的人做一個參考。

//============================================================================
// Name        : Convert_cpp.cpp
// Author      : roger
// Version     :
// Copyright   : Your copyright notice
// Description : Hello World in C++, Ansi-style
//============================================================================

#include <iostream>
#include <math.h>
#include <stdlib.h>
#include <iomanip>
using namespace std;

const double pi = 3.14159265358979324;
    const double a = 6378245.0;
    const double ee = 0.00669342162296594323;
    const  double x_pi = 3.14159265358979324 * 3000.0 / 180.0;

    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;
    }

     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;
    }

     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;
    }

    /**
     * 地球座標轉換爲火星座標
     * World Geodetic System ==> Mars Geodetic System
     *
     * @param wgLat  地球座標
     * @param wgLon
     *
     * mglat,mglon 火星座標
     */
     void transform2Mars(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;

    }

    /**
     * 火星座標轉換爲百度座標
     * @param gg_lat
     * @param gg_lon
     */
     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;

    }

    /**
     * 百度轉火星
     * @param bd_lat
     * @param bd_lon
     */
     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);

    }



int main() {
    double lat = 30.227607;
    double lon = 120.036565;

    //真實的經緯度轉化爲百度地圖上的經緯度，便於計算百度POI
    double marsLat = 0;
    double marsLon = 0;
    double resultLat = 0;
    double resultLon = 0;
    transform2Mars(lat,lon,marsLat,marsLon);
    bd_encrypt(marsLat,marsLon,resultLat,resultLon);

    //30.2193456 120.0348264
    cout<<setprecision(10)<<resultLat<<" "<<setprecision(10)<<resultLon<<endl;

}

         在自己的項目中，將標準GPS轉換爲百度座標後，再獲取百度地圖的POI，獲得的POI位置和真實的位置基本一致，可以使用。下面還有一些網上轉載的資料，比較豐富，值得一看。

-----------------------------------------------------------------------------------------------------

1. GCJ-02座標系統（火星座標）簡介：http://blog.csdn.net/giswens/article/details/8775121（存檔：http://mapbd.com/cms/2012/07/25/）
   
2. GCJ-02到真實座標反向變換的理論基礎：http://blog.csdn.net/giswens/article/details/8775213
   
3. WGS84座標與Web墨卡託座標互轉：http://blog.csdn.net/giswens/article/details/9634261
   
4. 地球座標系 (WGS-84) 到火星座標系 (GCJ-02) 的轉換算法：
   • C#代碼：https://on4wp7.codeplex.com/SourceControl/changeset/view/21483#353936（註解：http://blog.csdn.net/giswens/article/details/8775283）
   • Java代碼：http://emq.googlecode.com/svn/emq/src/Algorithm/Coords/Converter.java
   • iOS代碼：http://blog.csdn.net/giswens/article/details/8775183（存檔：http://www.keakon.net/2011/07/02/WGS84座標轉火星座標（iOS篇））
5. 火星座標系 (GCJ-02) 與百度座標系 (BD-09) 的轉換算法：http://blog.csdn.net/coolypf/article/details/8569813（解釋：國際經緯度座標標準爲WGS-84,國內必須至少使用國測局制定的GCJ-02,對地理位置進行首次加密。百度座標在此基礎上，進行了BD-09二次加密措施,更加保護了個人隱私。百度對外接口的座標系並不是GPS採集的真實經緯度，需要通過座標轉換接口進行轉換。）
   
6. 一種根據糾偏數據對火星座標進行完美擬合的方法：http://blog.sina.com.cn/s/blog_538036cf0100pxbl.html
   
7. 國內各地圖API座標系統比較：http://rovertang.com/archives/547
   
8. 關於百度地圖座標轉換接口的研究：http://rovertang.com/archives/24699

查詢過資料可得，通過程序進行 標準座標轉火星座標（google、高德），火星座標轉百度座標，百度座標轉火星座標，百度座標轉標準座標都是可行的，火星直接轉標準座標轉不了。