iOS地图开发2-坐标系的转换(swift) 发表于 2015-07-25 | 分类于 ios iOS地图开发2(swift) 代码123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211//// ZJ_MapKits//// Created by PSVMC on 15/6/8.// Copyright (c) 2015年 PSVMC. All rights reserved.//import Foundationimport CoreLocationclass ZJ_MapKits{ //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.14159265358979324; let xPi = M_PI * 3000.0 / 180.0; //WGS-84 --> GCJ-02 func transformFromWGSToGCJ(wgsLoc:CLLocationCoordinate2D)->CLLocationCoordinate2D { var adjustLoc=CLLocationCoordinate2D(); if( isLocationOutOfChina(wgsLoc)) { adjustLoc = wgsLoc; } else { var adjustLat = transformLatWithX(wgsLoc.longitude - 105.0 ,y:wgsLoc.latitude - 35.0); var adjustLon = transformLonWithX(wgsLoc.longitude - 105.0 ,y:wgsLoc.latitude - 35.0); var radLat = wgsLoc.latitude / 180.0 * pi; var magic = sin(radLat); magic = 1 - ee * magic * magic; var 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 func transformFromGCJToBaidu(p:CLLocationCoordinate2D) -> CLLocationCoordinate2D { var z = sqrt(p.longitude * p.longitude + p.latitude * p.latitude) + 0.00002 * sqrt(p.latitude * pi); var 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 func transformFromBaiduToGCJ(p:CLLocationCoordinate2D)-> CLLocationCoordinate2D { var x = p.longitude - 0.0065, y = p.latitude - 0.006; var z = sqrt(x * x + y * y) - 0.00002 * sin(y * xPi); var 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 func transformFromGCJToWGS(p:CLLocationCoordinate2D) -> CLLocationCoordinate2D { var 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; var maxIteration = 30; // Binary search while(true) { var leftBottom = transformFromWGSToGCJ(CLLocationCoordinate2D(latitude: minLat,longitude: minLng)); var rightBottom = transformFromWGSToGCJ(CLLocationCoordinate2D(latitude : minLat,longitude : maxLng)); var leftUp = transformFromWGSToGCJ(CLLocationCoordinate2D(latitude : maxLat,longitude : minLng)); var 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 func transformFromWGSToBaidu(p:CLLocationCoordinate2D) -> CLLocationCoordinate2D { var gcj = transformFromWGSToGCJ(p); var bd = transformFromGCJToBaidu(gcj) return bd; } //BD-09 --> WGS-84 func transformFromBaiduToWGS(p:CLLocationCoordinate2D) -> CLLocationCoordinate2D { var gcj = transformFromBaiduToGCJ(p) var 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; } } ///获取两点之间的距离 static func distanceByPoint(lat1:Double,lat2 :Double,lng1 :Double,lng2:Double)->Double{ var dd = M_PI/180; var x1=lat1*dd; var x2=lat2*dd; var y1=lng1*dd; var y2=lng2*dd; var R = 6371004; var temp = 2 - 2 * cos(x1) * cos(x2) * cos(y1-y2) - 2 * sin(x1) * sin(x2); var distance = Double(2) * Double(R) * asin(sqrt(temp)/2); //返回 m return distance; } ///获取两点之间的距离 static func distanceByPoint(point1:CLLocationCoordinate2D,point2:CLLocationCoordinate2D)->Double{ return distanceByPoint(point1.latitude, lat2: point2.latitude, lng1: point1.longitude, lng2: point2.longitude); } }