坐标转换



import UIKit

import Foundation
import CoreLocation

@objcMembers
public class HTMCoorTransform: NSObject{
    
    static let shared = HTMCoorTransform.init()
    
    //WGS-84:是国际标准,GPS坐标(Google Earth使用、或者GPS模块)
    //GCJ-02:中国坐标偏移标准,Google Map、高德、腾讯使用
    //BD-09: 百度坐标偏移标准,Baidu Map使用
    
    let  a = 6378245.0;
    let  ee = 0.00669342162296594323;
    let  pi = 3.1415926535897932384626;
    let  xPi = Double.pi * 3000.0 / 180.0;
    
    //WGS-84 --> GCJ-02
    public static func transformFromWGSToGCJ(_ wgsLoc:CLLocationCoordinate2D)->CLLocationCoordinate2D{
        return shared.transformFromWGSToGCJ(wgsLoc)
    }
    func transformFromWGSToGCJ(_ wgsLoc:CLLocationCoordinate2D)->CLLocationCoordinate2D
    {
        var adjustLoc=CLLocationCoordinate2D();
        if( isLocationOutOfChina(location: wgsLoc))
        {
            adjustLoc = wgsLoc;
        }
        else
        {
            var adjustLat = transformLatWithX(wgsLoc.longitude - 105.0 ,wgsLoc.latitude - 35.0);
            var adjustLon = transformLonWithX(wgsLoc.longitude - 105.0 ,wgsLoc.latitude - 35.0);
            let radLat = wgsLoc.latitude / 180.0 * pi;
            var magic = sin(radLat);
            magic = 1 - ee * magic * magic;
            let sqrtMagic = sqrt(magic);
            adjustLat = (adjustLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * pi);
            adjustLon = (adjustLon * 180.0) / (a / sqrtMagic * cos(radLat) * pi);
            adjustLoc.latitude = wgsLoc.latitude + adjustLat;
            adjustLoc.longitude = wgsLoc.longitude + adjustLon;
        }
        return adjustLoc;
    }
    
    func transformLatWithX(_ x:Double,_ y:Double)->Double
    {
        var lat = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y ;
        lat += 0.2 * sqrt(fabs(x));
        
        lat += (20.0 * sin(6.0 * x * pi)) * 2.0 / 3.0;
        lat += (20.0 * sin(2.0 * x * pi)) * 2.0 / 3.0;
        lat += (20.0 * sin(y * pi)) * 2.0 / 3.0;
        lat += (40.0 * sin(y / 3.0 * pi)) * 2.0 / 3.0;
        lat += (160.0 * sin(y / 12.0 * pi)) * 2.0 / 3.0;
        lat += (320 * sin(y * pi / 30.0)) * 2.0 / 3.0;
        return lat;
    }
    
    func transformLonWithX(_ x:Double,_ y:Double)->Double
    {
        var lon = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y ;
        lon +=  0.1 * sqrt(fabs(x));
        lon += (20.0 * sin(6.0 * x * pi)) * 2.0 / 3.0;
        lon += (20.0 * sin(2.0 * x * pi)) * 2.0 / 3.0;
        lon += (20.0 * sin(x * pi)) * 2.0 / 3.0;
        lon += (40.0 * sin(x / 3.0 * pi)) * 2.0 / 3.0;
        lon += (150.0 * sin(x / 12.0 * pi)) * 2.0 / 3.0;
        lon += (300.0 * sin(x / 30.0 * pi)) * 2.0 / 3.0;
        return lon;
    }
    
    //GCJ-02 --> BD-09
    public static func transformFromGCJToBaidu(_ wgsLoc:CLLocationCoordinate2D)->CLLocationCoordinate2D{
        return shared.transformFromGCJToBaidu(wgsLoc)
    }
    func transformFromGCJToBaidu(_ p:CLLocationCoordinate2D) -> CLLocationCoordinate2D
    {
        let z = sqrt(p.longitude * p.longitude + p.latitude * p.latitude) + 0.00002 * sqrt(p.latitude * pi);
        let theta = atan2(p.latitude, p.longitude) + 0.000003 * cos(p.longitude * pi);
        var geoPoint=CLLocationCoordinate2D();
        geoPoint.latitude  = (z * sin(theta) + 0.006);
        geoPoint.longitude = (z * cos(theta) + 0.0065);
        return geoPoint;
    }
    
    
    //BD-09 --> GCJ-02
    public static func transformFromBaiduToGCJ(_ wgsLoc:CLLocationCoordinate2D)->CLLocationCoordinate2D{
        return shared.transformFromBaiduToGCJ(wgsLoc)
    }
    func transformFromBaiduToGCJ(_ p:CLLocationCoordinate2D)-> CLLocationCoordinate2D
    {
        let x = p.longitude - 0.0065, y = p.latitude - 0.006;
        let z = sqrt(x * x + y * y) - 0.00002 * sin(y * xPi);
        let theta = atan2(y, x) - 0.000003 * cos(x * xPi);
        var geoPoint = CLLocationCoordinate2D();
        geoPoint.latitude  = z * sin(theta);
        geoPoint.longitude = z * cos(theta);
        return geoPoint;
    }
    
    //GCJ-02 --> WGS-84
    public static func transformFromGCJToWGS(_ wgsLoc:CLLocationCoordinate2D)->CLLocationCoordinate2D{
        return shared.transformFromGCJToWGS(wgsLoc)
    }
    func transformFromGCJToWGS(_ p:CLLocationCoordinate2D) -> CLLocationCoordinate2D
    {
        let threshold = 0.00001;
        
        // The boundary
        var minLat = p.latitude - 0.5;
        var maxLat = p.latitude + 0.5;
        var minLng = p.longitude - 0.5;
        var maxLng = p.longitude + 0.5;
        
        var delta = 1.0;
        let maxIteration = 30;
        // Binary search
        while(true)
        {
            let leftBottom  = transformFromWGSToGCJ(CLLocationCoordinate2D(latitude: minLat,longitude: minLng));
            let rightBottom = transformFromWGSToGCJ(CLLocationCoordinate2D(latitude : minLat,longitude : maxLng));
            let leftUp      = transformFromWGSToGCJ(CLLocationCoordinate2D(latitude : maxLat,longitude : minLng));
            let midPoint    = transformFromWGSToGCJ(CLLocationCoordinate2D(latitude : ((minLat + maxLat) / 2),longitude : ((minLng + maxLng) / 2)));
            delta = fabs(midPoint.latitude - p.latitude) + fabs(midPoint.longitude - p.longitude);
            
            if(maxIteration <= 1 || delta <= threshold)
            {
                return CLLocationCoordinate2D(latitude: (minLat + maxLat) / 2, longitude: (minLng + maxLng) / 2);
                
            }
            
            if(isContains(p, p1: leftBottom, p2: midPoint))
            {
                maxLat = (minLat + maxLat) / 2;
                maxLng = (minLng + maxLng) / 2;
            }
            else if(isContains(p, p1: rightBottom, p2: midPoint))
            {
                maxLat = (minLat + maxLat) / 2;
                minLng = (minLng + maxLng) / 2;
            }
            else if(isContains(p, p1: leftUp, p2: midPoint))
            {
                minLat = (minLat + maxLat) / 2;
                maxLng = (minLng + maxLng) / 2;
            }
            else
            {
                minLat = (minLat + maxLat) / 2;
                minLng = (minLng + maxLng) / 2;
            }
        }
        
    }
    
    //WGS-84 --> BD-09
    public static func transformFromWGSToBaidu(_ wgsLoc:CLLocationCoordinate2D)->CLLocationCoordinate2D{
        return shared.transformFromWGSToBaidu(wgsLoc)
    }
    func transformFromWGSToBaidu(_ p:CLLocationCoordinate2D) -> CLLocationCoordinate2D
    {
        let gcj = transformFromWGSToGCJ(p);
        let bd = transformFromGCJToBaidu(gcj)
        return bd;
    }
    
    //BD-09 --> WGS-84
    public static func transformFromBaiduToWGS(_ wgsLoc:CLLocationCoordinate2D)->CLLocationCoordinate2D{
        return shared.transformFromBaiduToWGS(wgsLoc)
    }
    func transformFromBaiduToWGS(_ p:CLLocationCoordinate2D) -> CLLocationCoordinate2D
    {
        let gcj = transformFromBaiduToGCJ(p)
        let wgs = transformFromGCJToWGS(gcj)
        return wgs;
    }
    
    //判断点是否在p1和p2之间
    //point: 点
    //p1:    左上角
    //p2:    右下角
    func isContains(_ point:CLLocationCoordinate2D , p1:CLLocationCoordinate2D, p2:CLLocationCoordinate2D)->Bool
    {
        
        return (point.latitude >= min(p1.latitude, p2.latitude) && point.latitude <= max(p1.latitude, p2.latitude)) && (point.longitude >= min(p1.longitude,p2.longitude) && point.longitude <= max(p1.longitude, p2.longitude));
    }
    
    
    //是否在中国以外
    func isLocationOutOfChina(location:CLLocationCoordinate2D) -> Bool
    {
        if (location.longitude < 72.004 || location.longitude > 137.8347 || location.latitude < 0.8293 || location.latitude > 55.8271){
            return true;
        }else{
            return false;
        }
        
    }
    
    ///获取两点之间的距离
    public static func distanceByPoint(lat1:Double,lat2 :Double,lng1 :Double,lng2:Double)->Double{
        let dd = Double.pi/180;
        let x1=lat1*dd;
        let x2=lat2*dd;
        let y1=lng1*dd;
        let y2=lng2*dd;
        let R = 6371004;
        
        let temp = 2 - 2 * cos(x1) * cos(x2) * cos(y1-y2) - 2 * sin(x1) * sin(x2);
        
        let distance = Double(2) * Double(R) * asin(sqrt(temp)/2);
        
        //返回 m
        return   distance;
        
    }
    
    ///获取两点之间的距离
    public static  func distanceByPoint(point1:CLLocationCoordinate2D,point2:CLLocationCoordinate2D)->Double{
        return distanceByPoint(lat1: point1.latitude, lat2: point2.latitude, lng1: point1.longitude, lng2: point2.longitude);
        
    }
}


©著作权归作者所有,转载或内容合作请联系作者
禁止转载,如需转载请通过简信或评论联系作者。
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 217,406评论 6 503
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,732评论 3 393
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 163,711评论 0 353
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,380评论 1 293
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,432评论 6 392
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,301评论 1 301
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,145评论 3 418
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,008评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,443评论 1 314
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,649评论 3 334
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,795评论 1 347
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,501评论 5 345
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,119评论 3 328
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,731评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,865评论 1 269
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,899评论 2 370
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,724评论 2 354

推荐阅读更多精彩内容