IOS/JS坐标系转换工具

要淡定哦

于 2023-02-21 16:20:47 发布

阅读量310

点赞数

文章标签： swift javascript Powered by 金山文档

本文链接：https://blog.csdn.net/m0_73942546/article/details/129144824

版权

该文章提供了iOS和JavaScript的代码实现，用于处理不同地图坐标系之间的转换，包括从WGS84坐标系转换到GCJ02（火星坐标系）和BD09（百度坐标系）。这些转换在处理中国地区的地理位置数据时尤其重要，因为国内地图服务通常使用加密的坐标系。

摘要由CSDN通过智能技术生成

开发处理过一些地图类项目，对坐标系转换做一些工具类

iOS系统上通过定位服务CLLocation相关接口获取定位信息时，获取的经纬度坐标系是WGS84地球坐

WGS84坐标系：即地球坐标系，国际上通用的坐标系。设备一般包含GPS芯片或者北斗芯片获取的经纬度为WGS84地理坐标系,谷歌地图采用的是WGS84地理坐标系（中国范围除外）

GCJ02坐标系：即火星坐标系，是由中国国家测绘局制订的地理信息系统的坐标系统。由WGS84坐标系经加密后的坐标系。谷歌中国地图和搜搜中国地图采用的是GCJ02地理坐标系;

BD09坐标系：即百度坐标系，GCJ02坐标系经加密后的坐标系

搜狗坐标系、图吧坐标系等，估计也是在GCJ02基础上加密而成的

北斗芯片获取的经纬度为WGS84地理坐标

swift 版

class GPSTool {

    static let BAIDU_LBS_TYPE = "bd09ll"
    static let pi = 3.1415926535897932384626;
    static let a = 6378245.0
    static let ee = 0.00669342162296594323
    static let x_pi = pi * 3000.0 / 180.0
    
    static func gps84_To_Gcj02(lon:Double, lat:Double) ->CLLocationCoordinate2D?{
        if GPSTool.outOfChina(lon: lon, lat: lat) {
            return nil
        }
        var dLat = GPSTool.transformLat(x: lon - 105.0, y: lat - 35.0)
        var dLon = GPSTool.transformLon(x: lon - 105.0, y: lat - 35.0)
        let radLat = (lat / 180.0) * pi
        var magic = sin(radLat)
        magic = 1 - ee * magic * magic
        let sqrtMagic = sqrt(magic)
        dLat = (dLat * 180.0) / (((a * (1 - ee)) / (magic * sqrtMagic)) * pi)
        dLon = (dLon * 180.0) / ((a / sqrtMagic) * cos(radLat) * pi)
        let mgLat = lat + dLat
        let mgLon = lon + dLon
        return  CLLocationCoordinate2D(latitude: mgLat, longitude: mgLon)
    }
    
    static func gcj02_To_Bd09(gg_lon:Double, gg_lat:Double)->CLLocationCoordinate2D {
        let x = gg_lon
        let y = gg_lat
        let z = sqrt(x * x + y * y) + 0.00002 * sin(y * x_pi);
        let theta = atan2(y, x) + 0.000003 * cos(x * x_pi);
        let bd_lon = z * cos(theta) + 0.0065;
        let bd_lat = z * sin(theta) + 0.006;
        return CLLocationCoordinate2D(latitude: bd_lat, longitude: bd_lon)
    }
    
    static func gps84_To_Bd09(lon:Double, lat:Double) ->CLLocationCoordinate2D?{
        if GPSTool.outOfChina(lon: lon, lat: lat) {
            return nil
        }
        if let gcl = GPSTool.gps84_To_Gcj02(lon: lon, lat: lat) {
            let bd = GPSTool.gcj02_To_Bd09(gg_lon: gcl.longitude, gg_lat: gcl.latitude)
            return bd
        }else{
            return nil
        }
    }

    static func gcj_To_Gps84(lon:Double, lat:Double) ->CLLocationCoordinate2D{
        let gps = GPSTool.transform(lon:lon, lat:lat)
        let longitude = lon * 2 - gps.longitude
        let latitude = lat * 2 - gps.latitude
        return CLLocationCoordinate2D(latitude: latitude, longitude: longitude)
    }
    
    static func bd09_To_Gcj02(bd_lon:Double, bd_lat:Double) ->CLLocationCoordinate2D{
        let x = bd_lon - 0.0065
        let y = bd_lat - 0.006
        let z = sqrt(x * x + y * y) - 0.00002 * sin(y * x_pi)
        let theta = atan2(y, x) - 0.000003 * cos(x * x_pi)
        let gg_lon = z * cos(theta)
        let gg_lat = z * sin(theta)
        return CLLocationCoordinate2D(latitude: gg_lat, longitude: gg_lon)
    }
    
    static func bd09_To_Gps84(bd_lon:Double, bd_lat:Double) ->CLLocationCoordinate2D{
        let gcj02 = GPSTool.bd09_To_Gcj02(bd_lon: bd_lon, bd_lat: bd_lat)
      let map84 = GPSTool.gcj_To_Gps84(lon: gcj02.longitude, lat: gcj02.latitude)
      return map84
    }

    static func transform(lon:Double, lat:Double) ->CLLocationCoordinate2D{
        if (GPSTool.outOfChina(lon: lon, lat: lat)) {
          return CLLocationCoordinate2D(latitude: lat, longitude: lon)
        }
        var dLat = GPSTool.transformLat(x: lon - 105.0, y: lat - 35.0)
        var dLon = GPSTool.transformLon(x: lon - 105.0, y: lat - 35.0);
        let radLat = (lat / 180.0) * pi
        var magic = sin(radLat)
        magic = 1 - ee * magic * magic
        let sqrtMagic = sqrt(magic)
        dLat = (dLat * 180.0) / (((a * (1 - ee)) / (magic * sqrtMagic)) * pi)
        dLon = (dLon * 180.0) / ((a / sqrtMagic) * cos(radLat) * pi)
        let mgLat = lat + dLat
        let mgLon = lon + dLon
        return CLLocationCoordinate2D(latitude: mgLat, longitude: mgLon)
    }
    
    static func outOfChina(lon:Double, lat:Double) ->Bool{
        if (lon < 72.004 || lon > 137.8347) {
          return true
        }
        if (lat < 0.8293 || lat > 55.8271) {
          return true
        }
        return false
    }
    
    static 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 * 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
    }
    
    static 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 * 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
    }
}

js 版

gps 对象

class Gps {
  constructor(longitude, latitude) {
    this.longitude = longitude;
    this.latitude = latitude;
  }

  setWgLongitude(longitude) {
    this.longitude = longitude;
  }

  setWgLatitude(latitude) {
    this.latitude = latitude;
  }

  getWgLongitude(longitude) {
    return this.longitude;
  }

  getWgLatitude(latitude) {
    return this.latitude;
  }

  toString() {
    return this.longitude + ',' + this.latitude;
  }

  toJson(count = 12) {
    let result = {
      longitude: this.longitude.toFixed(count) * 1,
      latitude: this.latitude.toFixed(count) * 1,
    };
    return result;
  }
}
module.exports = Gps;

工具

const pi = 3.1415926535897932384626;
const a = 6378245.0;
const ee = 0.00669342162296594323;
const x_pi = pi * 3000.0 / 180.0;

class GPSTool {

  static gps84_To_Gcj02(lon, lat) {
    if (GPSTool.outOfChina(lon, lat)) {
      return null;
    }
    let dLat = GPSTool.transformLat(lon - 105.0, lat - 35.0);
    let dLon = GPSTool.transformLon(lon - 105.0, lat - 35.0);
    let radLat = (lat / 180.0) * pi;
    let magic = Math.sin(radLat);
    magic = 1 - ee * magic * magic;
    let sqrtMagic = Math.sqrt(magic);
    dLat = (dLat * 180.0) / (((a * (1 - ee)) / (magic * sqrtMagic)) * pi);
    dLon = (dLon * 180.0) / ((a / sqrtMagic) * Math.cos(radLat) * pi);
    let mgLat = lat + dLat;
    let mgLon = lon + dLon;
    return new Gps(mgLon, mgLat);
  }

  static gps84_To_Bd09(lon, lat) {
    let gcl = GPSTool.gps84_To_Gcj02(lon, lat);
    let bd = GPSTool.gcj02_To_Bd09(gcl.getWgLongitude(), gcl.getWgLatitude());
    return bd;
  }

  static gcj_To_Gps84(lon, lat) {
    let gps = GPSTool.transform(lon, lat);
    let longitude = lon * 2 - gps.getWgLongitude();
    let latitude = lat * 2 - gps.getWgLatitude();
    return new Gps(longitude, latitude);
  }

  static gcj02_To_Bd09(gg_lon, gg_lat) {
    let x = gg_lon,
      y = gg_lat;
    let z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * x_pi);
    let theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * x_pi);
    let bd_lon = z * Math.cos(theta) + 0.0065;
    let bd_lat = z * Math.sin(theta) + 0.006;
    let gps = new Gps(bd_lon, bd_lat);
    return gps;
  }

  static bd09_To_Gcj02(bd_lon, bd_lat) {
    let x = bd_lon - 0.0065,
      y = bd_lat - 0.006;
    let z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * x_pi);
    let theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_pi);
    let gg_lon = z * Math.cos(theta);
    let gg_lat = z * Math.sin(theta);
    return new Gps(gg_lon, gg_lat);
  }

  static bd09_To_Gps84(bd_lon, bd_lat) {
    let gcj02 = GPSTool.bd09_To_Gcj02(bd_lon, bd_lat);
    let map84 = GPSTool.gcj_To_Gps84(
      gcj02.getWgLongitude(),
      gcj02.getWgLatitude(),
    );
    return map84;
  }

  static outOfChina(lon, lat) {
    if (lon < 72.004 || lon > 137.8347) {
      return true;
    }
    if (lat < 0.8293 || lat > 55.8271) {
      return true;
    }
    return false;
  }

  static transform(lon, lat) {
    if (GPSTool.outOfChina(lon, lat)) {
      return new Gps(lon, lat);
    }
    let dLat = GPSTool.transformLat(lon - 105.0, lat - 35.0);
    let dLon = GPSTool.transformLon(lon - 105.0, lat - 35.0);
    let radLat = (lat / 180.0) * pi;
    let magic = Math.sin(radLat);
    magic = 1 - ee * magic * magic;
    let sqrtMagic = Math.sqrt(magic);
    dLat = (dLat * 180.0) / (((a * (1 - ee)) / (magic * sqrtMagic)) * pi);
    dLon = (dLon * 180.0) / ((a / sqrtMagic) * Math.cos(radLat) * pi);
    let mgLat = lat + dLat;
    let mgLon = lon + dLon;
    return new Gps(mgLon, mgLat);
  }

  static transformLat(x, y) {
    var ret =
      -100.0 +
      2.0 * x +
      3.0 * y +
      0.2 * y * y +
      0.1 * x * y +
      0.2 * Math.sqrt(Math.abs(x));
    ret +=
      ((20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0) /
      3.0;
    ret +=
      ((20.0 * Math.sin(y * pi) + 40.0 * Math.sin((y / 3.0) * pi)) * 2.0) / 3.0;
    ret +=
      ((160.0 * Math.sin((y / 12.0) * pi) + 320 * Math.sin((y * pi) / 30.0)) *
        2.0) /
      3.0;
    return ret;
  }

  static transformLon(x, y) {
    var ret =
      300.0 +
      x +
      2.0 * y +
      0.1 * x * x +
      0.1 * x * y +
      0.1 * Math.sqrt(Math.abs(x));
    ret +=
      ((20.0 * Math.sin(6.0 * x * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0) /
      3.0;
    ret +=
      ((20.0 * Math.sin(x * pi) + 40.0 * Math.sin((x / 3.0) * pi)) * 2.0) / 3.0;
    ret +=
      ((150.0 * Math.sin((x / 12.0) * pi) + 300.0 * Math.sin((x / 30.0) * pi)) *
        2.0) /
      3.0;
    return ret;
  }
}

export default GPSTool;