火星坐标系（GCJ-02）与百度坐标系（BD-09）转换

算法代码如下，其中 bd_encrypt 将 GCJ-02 坐标转换成 BD-09 坐标， bd_decrypt 反之。

关于 GCJ-02 和 BD-09 ，请参考 http://developer.baidu.com/map/question.htm#qa0043 。

1. #include <math.h>
2. const double x_pi = 3.14159265358979324 * 3000.0 / 180.0;
3. void bd_encrypt(double gg_lat, double gg_lon, double &bd_lat, double &bd_lon)
4. {
5. double x = gg_lon, y = gg_lat;
6. double z = sqrt(x * x + y * y) + 0.00002 * sin(y * x_pi);
7. double theta = atan2(y, x) + 0.000003 * cos(x * x_pi);
8. bd_lon = z * cos(theta) + 0.0065;
9. bd_lat = z * sin(theta) + 0.006;
10. }
11. void bd_decrypt(double bd_lat, double bd_lon, double &gg_lat, double &gg_lon)
12. {
13. double x = bd_lon - 0.0065, y = bd_lat - 0.006;
14. double z = sqrt(x * x + y * y) - 0.00002 * sin(y * x_pi);
15. double theta = atan2(y, x) - 0.000003 * cos(x * x_pi);
16. gg_lon = z * cos(theta);
17. gg_lat = z * sin(theta);
18. }