|
|
@@ -1,94 +0,0 @@
|
|
|
-//
|
|
|
-// Coordinate+E.swift
|
|
|
-// Pods
|
|
|
-//
|
|
|
-// Created by Groot on 2025/5/15.
|
|
|
-//
|
|
|
-
|
|
|
-import Foundation
|
|
|
-import CoreLocation
|
|
|
-
|
|
|
-extension CLLocationCoordinate2D {
|
|
|
- /// 地球坐标 WGC-84 转为为火星坐标 GCJ-02(限中国境内)
|
|
|
- /// - CLLocationManager 获取的是地球坐标,但地图上显示的需要火星坐标
|
|
|
- /// - 通过 CLLocationManager 获取到中国境内地理坐标,若需要在地图上显示,则需要转换为火星坐标
|
|
|
- var wgc84ToGCJ02: CLLocationCoordinate2D {
|
|
|
- let earthRadius = 6378245.0
|
|
|
- let gcj_ee = 0.00669342162296594323
|
|
|
-
|
|
|
- var adjustedLatitude = transformedLatitude
|
|
|
- var adjustedLongitude = transformedLongitude
|
|
|
-
|
|
|
- let radLat = latitude / 180.0 * .pi
|
|
|
- let magic: Double = {
|
|
|
- let temp = sin(radLat)
|
|
|
- let result = 1 - gcj_ee * temp * temp
|
|
|
- return result
|
|
|
- }()
|
|
|
- let sqrtMagic = sqrt(magic)
|
|
|
-
|
|
|
- adjustedLatitude = (adjustedLatitude * 180.0) / ((earthRadius * (1.0 - gcj_ee)) / (magic * sqrtMagic) * .pi)
|
|
|
- adjustedLongitude = (adjustedLongitude * 180.0) / (earthRadius / sqrtMagic * cos(radLat) * .pi)
|
|
|
-
|
|
|
- let result = CLLocationCoordinate2D(
|
|
|
- latitude: latitude + adjustedLatitude,
|
|
|
- longitude: longitude + adjustedLongitude
|
|
|
- )
|
|
|
-
|
|
|
- return result
|
|
|
- }
|
|
|
-
|
|
|
- var gcj02ToWGS84: CLLocationCoordinate2D {
|
|
|
- let a = 6378245.0
|
|
|
- let ee = 0.00669342162296594323
|
|
|
- let gcjLat = latitude
|
|
|
- let gcjLon = longitude
|
|
|
- var dLat = transformLat(x: gcjLon - 105.0, y: gcjLat - 35.0)
|
|
|
- var dLon = transformLon(x: gcjLon - 105.0, y: gcjLat - 35.0)
|
|
|
- let radLat = gcjLat / 180.0 * Double.pi
|
|
|
- var magic = sin(radLat)
|
|
|
- magic = 1 - ee * magic * magic
|
|
|
- let sqrtMagic = sqrt(magic)
|
|
|
- dLat = (dLat * 180.0) / ((a * (1 - ee)) / (magic * sqrtMagic) * Double.pi)
|
|
|
- dLon = (dLon * 180.0) / (a / sqrtMagic * cos(radLat) * Double.pi)
|
|
|
- let mgLat = gcjLat + dLat
|
|
|
- let mgLon = gcjLon + dLon
|
|
|
- return CLLocationCoordinate2D(latitude: gcjLat * 2 - mgLat, longitude: gcjLon * 2 - mgLon)
|
|
|
- }
|
|
|
-
|
|
|
- private var transformedLatitude: Double {
|
|
|
- let x = longitude - 105.0
|
|
|
- let y = latitude - 35.0
|
|
|
- var result = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * sqrt(abs(x))
|
|
|
- result += (20.0 * sin(6.0 * x * .pi) + 20.0 * sin(2.0 * x * .pi)) * 2.0 / 3.0
|
|
|
- result += (20.0 * sin(y * .pi) + 40.0 * sin(y / 3.0 * .pi)) * 2.0 / 3.0
|
|
|
- result += (160.0 * sin(y / 12.0 * .pi) + 320.0 * sin(y * .pi / 30.0)) * 2.0 / 3.0
|
|
|
- return result
|
|
|
- }
|
|
|
-
|
|
|
- private var transformedLongitude: Double {
|
|
|
- let x = longitude - 105.0
|
|
|
- let y = latitude - 35.0
|
|
|
- var result = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * sqrt(abs(x))
|
|
|
- result += (20.0 * sin(6.0 * x * .pi) + 20.0 * sin(2.0 * x * .pi)) * 2.0 / 3.0
|
|
|
- result += (20.0 * sin(x * .pi) + 40.0 * sin(x / 3.0 * .pi)) * 2.0 / 3.0
|
|
|
- result += (150.0 * sin(x / 12.0 * .pi) + 300.0 * sin(x / 30.0 * .pi)) * 2.0 / 3.0
|
|
|
- return result
|
|
|
- }
|
|
|
-
|
|
|
- func transformLat(x: Double, y: Double) -> Double {
|
|
|
- var 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 * Double.pi) + 20.0 * sin(2.0 * x * Double.pi)) * 2.0 / 3.0
|
|
|
- ret += (20.0 * sin(y * Double.pi) + 40.0 * sin(y / 3.0 * Double.pi)) * 2.0 / 3.0
|
|
|
- ret += (160.0 * sin(y / 12.0 * Double.pi) + 320 * sin(y * Double.pi / 30.0)) * 2.0 / 3.0
|
|
|
- return ret
|
|
|
- }
|
|
|
-
|
|
|
- func transformLon(x: Double, y: Double) -> Double {
|
|
|
- var 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 * Double.pi) + 20.0 * sin(2.0 * x * Double.pi)) * 2.0 / 3.0
|
|
|
- ret += (20.0 * sin(x * Double.pi) + 40.0 * sin(x / 3.0 * Double.pi)) * 2.0 / 3.0
|
|
|
- ret += (150.0 * sin(x / 12.0 * Double.pi) + 300.0 * sin(x / 30.0 * Double.pi)) * 2.0 / 3.0
|
|
|
- return ret
|
|
|
- }
|
|
|
-}
|