1、  GPS，WGS-84，原始坐标体系。一般用国际标准的GPS记录仪记录下来的坐标，

都是GPS的坐标。很可惜，在中国，任何一个地图产品都不允许使用GPS坐标，

据说是为了保密。GPS坐标形式如图，度分秒形式的经纬度。

2、  GCJ-02，国测局02年发布的坐标体系。又称“火星坐标”。

在中国，必须至少使用GCJ-02的坐标体系。比如谷歌，腾讯，高德都在

用这个坐标体系。GCJ-02也是国内最广泛使用的坐标体系。

3、  其他坐标体系。一般都是由GCJ-02进过偏移算法得到的。

这种体系就根据每个公司的不同，坐标体系都不一样了。比如，

百度和搜狗就使用自己的坐标体系，与其他坐标体系不兼容。

首先明白，所有坐标体系的原点，都是非洲。

1、  经纬度。这个是球面坐标，对于北京来说，就是(116.38817139.935961)这

样的坐标。比如腾讯、高德、百度都是这样的经纬度坐标。

谷歌是经纬度顺序写反的经纬度坐标。如果是度分秒坐标，需要进行转换

，才能得到这样的经纬度坐标。详见坐标转换。

2、  墨卡托坐标。平面坐标，相当于是直线距离，数字一般都比较大，像这样的。

(215362.00021333335 99526.00034912192)墨卡托坐标，

主要用于程序的后台计算。直线距离嘛，加加减减几乎计算方便。

搜狗地图API就是直接使用的墨卡托坐标。

在各种web端平台，或者高德、腾讯、百度上取到的坐标，都不是GPS坐标，

都是GCJ-02坐标，或者自己的偏移坐标系。

比如，你在谷歌地图API，高德地图API，腾讯地图API上取到的，

都是GCJ-02坐标，他们三家都是通用的，也适用于大部分地图API产品，

以及他们的地图产品。

例外，百度API上取到的，是BD-09坐标，只适用于百度地图相关产品。

例外，搜狗API上取到的，是搜狗坐标，只适用于搜狗地图相关产品。

例外，谷歌地球，google earth上取到的，是GPS坐标，

而且是度分秒形式的经纬度坐标。在国内不允许使用。必须转换为GCJ-02坐标。

百度地图坐标与苹果自带地图经纬度之间的相互转换方法：

/// 百度坐标转高德坐标

+ (CLLocationCoordinate2D)GCJ02FromBD09:(CLLocationCoordinate2D)coor

{

    CLLocationDegrees x_pi = 3.14159265358979324 * 3000.0 / 180.0;

    CLLocationDegrees x = coor.longitude - 0.0065, y = coor.latitude - 0.006;

    CLLocationDegrees z = sqrt(x * x + y * y) - 0.00002 * sin(y * x_pi);

    CLLocationDegrees theta = atan2(y, x) - 0.000003 * cos(x * x_pi);

    CLLocationDegrees gg_lon = z * cos(theta);

    CLLocationDegrees gg_lat = z * sin(theta);

    return CLLocationCoordinate2DMake(gg_lat, gg_lon);

}

// 高德坐标转百度坐标

+ (CLLocationCoordinate2D)BD09FromGCJ02:(CLLocationCoordinate2D)coor

{

    CLLocationDegrees x_pi = 3.14159265358979324 * 3000.0 / 180.0;

    CLLocationDegrees x = coor.longitude, y = coor.latitude;

    CLLocationDegrees z = sqrt(x * x + y * y) + 0.00002 * sin(y * x_pi);

    CLLocationDegrees theta = atan2(y, x) + 0.000003 * cos(x * x_pi);

    CLLocationDegrees bd_lon = z * cos(theta) + 0.0065;

    CLLocationDegrees bd_lat = z * sin(theta) + 0.006;

    return CLLocationCoordinate2DMake(bd_lat, bd_lon);

}

# -*- coding: utf-8 -*-

import math

x_pi = 3.14159265358979324 * 3000.0 / 180.0

def bd_encrypt(gg):

    x = gg["gg_lon"]

    y = gg["gg_lat"]

    z = math.sqrt(x * x + y * y) + 0.00002 * math.sin(y * x_pi)

    theta = math.atan2(y, x) + 0.000003 * math.cos(x * x_pi)

    bd_lon = z * math.cos(theta) + 0.0065

    bd_lat = z * math.sin(theta) + 0.006

    return {"bd_lon":bd_lon, "bd_lat":bd_lat}

def bd_decrypt(bd):

    x = bd["bd_lon"] - 0.0065

    y = bd["bd_lat"] - 0.006;

    z = sqrt(x * x + y * y) - 0.00002 * sin(y * x_pi);

    theta = atan2(y, x) - 0.000003 * cos(x * x_pi);

    gg_lon = z * cos(theta);

    gg_lat = z * sin(theta);

    return {"gg_lon":gg_lon, "gg_lat":gg_lat}