一、坐標(biāo)簡(jiǎn)介
目前國(guó)內(nèi)主流坐標(biāo)系類型主要有三種：WGS84、GCJ02、BD09；
1.WGS84：為一種大地坐標(biāo)系，也是目前廣泛使用的GPS全球衛(wèi)星定位系統(tǒng)使用的坐標(biāo)系，如：iOS系統(tǒng)自帶地圖坐標(biāo)。

2. GCJ02：是由中國(guó)國(guó)家測(cè)繪局制訂的地理信息系統(tǒng)的坐標(biāo)系統(tǒng)，是由WGS84坐標(biāo)系經(jīng)加密后的坐標(biāo)系，如：高德地圖、google地圖、soso地圖、aliyun地圖、mapabc地圖和amap地圖所用坐標(biāo)
3. BD09：百度坐標(biāo)系，在GCJ02坐標(biāo)系基礎(chǔ)上再次加密。其中BD09ll表示百度經(jīng)緯度坐標(biāo)，BD09mc表示百度墨卡托米制坐標(biāo)，如：百度地圖所用坐標(biāo)

二、一般項(xiàng)目集成的sdk都支持別的坐標(biāo)轉(zhuǎn)成自己的坐標(biāo)，如百度地圖支持WGS84和GCJ02轉(zhuǎn)BD09ll，也支持BD09ll轉(zhuǎn)GCJ02，但是，下面就主要介紹一下GCJ02和BD09ll轉(zhuǎn)WGS84

二、BD09ll轉(zhuǎn)GCJ02

/**
 百度坐標(biāo)轉(zhuǎn)高德坐標(biāo)

 @param glat 緯度
 @param glon 經(jīng)度
 @return 數(shù)組【轉(zhuǎn)換后的緯度，轉(zhuǎn)換后的經(jīng)度】
 */
+ (NSArray *)bd09llToGCJLat:(double)glat lon:(double)glon {
    double X_PI = M_PI * 3000.0 / 180.0;
    double x = glon - 0.0065;
    double y = glat - 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);
    double lat = z * sin(theta);
    double lon = z * cos(theta);
    NSArray *latlon = @[[NSNumber numberWithDouble:lat],[NSNumber numberWithDouble:lon]];
    
    return latlon;
}

三、GCJ02轉(zhuǎn)WGS84（只支持國(guó)內(nèi)）

/**
 高德坐標(biāo)轉(zhuǎn)系統(tǒng)自帶坐標(biāo)

 @param array 數(shù)組 需要轉(zhuǎn)換的坐標(biāo)【緯度，經(jīng)度】
 @return 數(shù)組【轉(zhuǎn)換后的緯度，轉(zhuǎn)換后的經(jīng)度】
 */
+ (NSArray *)gcjToWGS:(NSArray *)array {
    double lat = [array.firstObject doubleValue];
    double lon = [array.lastObject doubleValue];
    if ([self outOfChinalat:lat lon:lon]){
        return @[[NSNumber numberWithDouble:lat],[NSNumber numberWithDouble:lon]];
    }
    NSArray *latlon = [self deltalat:lat lon:lon];
    return latlon;
}

+ (BOOL)outOfChinalat:(double)lat lon:(double)lon {
    if (lon < 72.004 || lon > 137.8347)
        return true;
    if (lat < 0.8293 || lat > 55.8271)
        return true;
    return false;
}

+ (NSArray *)deltalat:(double)wgLat lon:(double)wgLon {
    double OFFSET = 0.00669342162296594323;
    double AXIS = 6378245.0;
    
    double dLat = [self transformLatx:(wgLon - 105.0) lon:(wgLat - 35.0)];
    double dLon = [self transformLonx:(wgLon - 105.0) lon:(wgLat - 35.0)];
    double radLat = wgLat / 180.0 * M_PI;
    double magic = sin(radLat);
    magic = 1 - OFFSET * magic * magic;
    double sqrtMagic = sqrt(magic);
    dLat = (dLat * 180.0) / ((AXIS * (1 - OFFSET)) / (magic * sqrtMagic) * M_PI);
    dLon = (dLon * 180.0) / (AXIS / sqrtMagic * cos(radLat) * M_PI);
    return @[[NSNumber numberWithDouble:(wgLat - dLat)],[NSNumber numberWithDouble:(wgLon - dLon)]];
}

+ (double)transformLatx:(double)x lon:(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 * M_PI) + 20.0 * sin(2.0 * x * M_PI)) * 2.0 / 3.0;
    ret += (20.0 * sin(y * M_PI) + 40.0 * sin(y / 3.0 * M_PI)) * 2.0 / 3.0;
    ret += (160.0 * sin(y / 12.0 * M_PI) + 320 * sin(y * M_PI / 30.0)) * 2.0 / 3.0;
    return ret;
}

+ (double)transformLonx:(double)x lon:(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 * M_PI) + 20.0 * sin(2.0 * x * M_PI)) * 2.0 / 3.0;
    ret += (20.0 * sin(x * M_PI) + 40.0 * sin(x / 3.0 * M_PI)) * 2.0 / 3.0;
    ret += (150.0 * sin(x / 12.0 * M_PI) + 300.0 * sin(x / 30.0 * M_PI)) * 2.0 / 3.0;
    return ret;
}

四、BD09ll轉(zhuǎn)WGS84，目前的做法是先將BD09ll坐標(biāo)轉(zhuǎn)成GCJ02坐標(biāo)，然后再轉(zhuǎn)成WGS84坐標(biāo)。

五、經(jīng)實(shí)測(cè)，轉(zhuǎn)換后的坐標(biāo)，另GCJ02坐標(biāo)轉(zhuǎn)WGS84坐標(biāo)有一個(gè)更精確的算法，不過(guò)里面設(shè)計(jì)到了遞歸。