百度官方本身就提供有谷歌坐标转百度坐标（BD-09）的API，而谷歌坐标中国部分采用的正是火星坐标（GCJ-02）。关于这个API的详细介绍可以参见这篇文章 http://gooderlee.iteye.com/blog/1186512。

本文我们主要来验证一下网上流传的火星坐标转百度坐标算法的准确性。

算法我最早是在这个地址看到的http://blog.csdn.net/coolypf/article/details/8569813，是C++的实现。

下面是对应的Python实现：

from __future__ import division
from math import pi, sqrt, sin, cos, atan2

# 火星坐标(GCJ-02)转百度坐标(BD-09)算法
X_PI = 3.14159265358979324 * 3000.0 / 180.0
def gcj02_to_bd09(lat, lng):
    z = sqrt(lng * lng + lat * lat) + 0.00002 * sin(lat * X_PI)
    theta = atan2(lat, lng) + 0.000003 * cos(lng * X_PI)
    bd_lng = z * cos(theta) + 0.0065
    bd_lat = z * sin(theta) + 0.006
    return bd_lat, bd_lng

if __name__ == '__main__':
    # 测试
    google_lat, google_lng = 31.2883064, 121.4904417
    print gc202_to_bd09(lat=google_lat, lng=google_lng)

验证一下准确性：

1）我们从大众点评网上选取一家商户http://www.dianping.com/shop/500000，其采用的是高德地图（也是火星坐标），其坐标值是：31.2883064, 121.4904417。

2）我们在谷歌地图上验证一下这个点，结果是吻合的。

3）使用上述算法转换后得到的百度坐标是：31.293981660287436, 121.49700476064442，我们在百度坐标拾取系统中进行验证，肉眼看这个点也是吻合的。详见下图所示（点击查看大图）。