如何正确、统一的处理iOS地图漂移？

Question

已知：iOS5、6使用不同的数据源；iOS6会根据运营商和IP地址使用不同的数据源；国内地图存在一定的「漂移」。

问：怎么在iOS5和iOS6下，无论设备在墙内墙外：
1　获得并存储真实的经纬度（wgs84）
2　在任意类型地图（Google/AutoNavi/TomTom）上标记到正确的位置。

Answer 1

基于你的问题，我觉得描述的还不是很详细，我补充几点：

1. 常入眼的坐标系包括:wgs84,这是全球通用的标准（除天朝），你的调用手机的gps模块返回的坐标就是这个坐标系；gcj02,这个是我天朝测绘局基于wgs84进行的非线性偏移后的坐标系，国内的地图厂商如果想通过资质审核使用gcj02是最基本，还有些基于gcj02进一步进行了加密(百度的bd09坐标系)；

2. 由于在国内，所有的电子地图都是经过偏移的，所以你必须使用gcj02的坐标才能在地图上准确显示位置，如果使用wgs84那必定会有问题，实际上就是你用了真实的经纬度放在了另一个坐标系的电子地图上，所以会出现偏移，但是wgs84的经纬度绝对是真实的。

3、国内你能拿到的数据基本都是gcj02的，因为使用wgs84的数据是非法的，应该是保密法之类啥的吧

4、那怎么办？没辙了么？跪了么？当然不会，现在的问题就是如何把手机上得到wgs84坐标转换到gcj02上，答案是肯定可以的，wgs84->gcj02的转换一般是由国家测绘局提供的，不过我们一般拿不到，不过你可以使用第三方的地图资质厂商的平台提供的接口（这个肯定有，因为我之前用过，如果找不到我帮你找找）。

5、如果你不想对经纬度转换来转换去的，因为这样毕竟增加了不必要的代码，如果这些不是你的核心业务的话，你可以使用国内第三方的地图平台的整体解决方案，一般是在客户端集成他们的sdk，然后这个sdk帮你从定位到上图搞定一切，无需担心中间的坐标转换的问题，因为sdk已经内部做好了。

6、对于在“在任意类型地图（Google/AutoNavi/TomTom）上标记到正确的位置。”这个问题，怎么说呢，除非他们的坐标系是一致的，否则你会很恶心，AutoNavi是gcj02,Google国内是gcj02国外是wgs84,TomTom我还真不知道，总之，祝你好运吧..

Answer 2

#import "LocationTransform.h"
//
// Krasovsky 1940
//
// a = 6378245.0, 1/f = 298.3
// b = a * (1 - f)
// ee = (a^2 - b^2) / a^2;
const double pi = 3.14159265358979324;
const double a = 6378245.0;
const double ee = 0.00669342162296594323;

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

///
///  WGS-84 到 GCJ-02 的转换
///
Location transformFromWGSToGCJ(Location wgLoc)
{
    Location mgLoc;
    if (outOfChina(wgLoc.lat, wgLoc.lng))
    {
        mgLoc = wgLoc;
        return mgLoc;
    }
    double dLat = transformLat(wgLoc.lng - 105.0, wgLoc.lat - 35.0);
    double dLon = transformLon(wgLoc.lng - 105.0, wgLoc.lat - 35.0);
    double radLat = wgLoc.lat / 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);
    mgLoc.lat = wgLoc.lat + dLat;
    mgLoc.lng = wgLoc.lng + dLon;

    return mgLoc;
}

///
///  GCJ-02 坐标转换成 BD-09 坐标
///

const double x_pi = 3.14159265358979324 * 3000.0 / 180.0;
Location bd_encrypt(Location gcLoc)
{
    double x = gcLoc.lng, y = gcLoc.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);
    return LocationMake(z * cos(theta) + 0.0065, z * sin(theta) + 0.006);
}

///
///   BD-09 坐标转换成 GCJ-02坐标
///
///
Location bd_decrypt(Location bdLoc)
{
    double x = bdLoc.lng - 0.0065, y = bdLoc.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);
    return LocationMake(z * cos(theta), z * sin(theta));
}

看代码吧...不用谢