坐标系转换

const ee = 0.00669342162296594323; // ee: 椭球的偏心率。
const a = 6378245.0; // a: 卫星椭球坐标投影到平面地图坐标系的投影因子。
const PI = 3.14159265358979323846264338327950288;
const pi = (PI * 3000.0) / 180.0;
/**
 *  判断经纬度是否超出中国境内
 */
const isLocationOutOfChina = (latitude, longitude) =>
  !!(
    longitude < 72.004 ||
    longitude > 137.8347 ||
    latitude < 0.8293 ||
    latitude > 55.8271
  );

const transformLatWithXY = (x, y) => {
  var lat =
    -100.0 +
    2.0 * x +
    3.0 * y +
    0.2 * y * y +
    0.1 * x * y +
    0.2 * Math.sqrt(Math.abs(x));
  lat +=
    ((20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0) /
    3.0;
  lat +=
    ((20.0 * Math.sin(y * PI) + 40.0 * Math.sin((y / 3.0) * PI)) * 2.0) / 3.0;
  lat +=
    ((160.0 * Math.sin((y / 12.0) * PI) + 320 * Math.sin((y * PI) / 30.0)) *
      2.0) /
    3.0;
  return lat;
};

const transformLonWithXY = (x, y) => {
  var lon =
    300.0 +
    x +
    2.0 * y +
    0.1 * x * x +
    0.1 * x * y +
    0.1 * Math.sqrt(Math.abs(x));
  lon +=
    ((20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0) /
    3.0;
  lon +=
    ((20.0 * Math.sin(x * PI) + 40.0 * Math.sin((x / 3.0) * PI)) * 2.0) / 3.0;
  lon +=
    ((150.0 * Math.sin((x / 12.0) * PI) + 300.0 * Math.sin((x / 30.0) * PI)) *
      2.0) /
    3.0;
  return lon;
};

/**
 *  将WGS-84(国际标准)转为GCJ-02(火星坐标):
 */
export const wgs84_to_gcj02 = (lat, lon) => {
  if (isLocationOutOfChina(latitude, longitude))
    return { latitude: lat, longitude: lon };
  let latitude = transformLatWithXY(lon - 105.0, lat - 35.0);
  let longitude = transformLonWithXY(lon - 105.0, lat - 35.0);
  let radLat = (lat / 180.0) * PI;
  let magic = Math.sin(radLat);
  magic = 1 - ee * magic * magic;
  const sqrtMagic = Math.sqrt(magic);
  latitude = (latitude * 180.0) / (((a * (1 - ee)) / (magic * sqrtMagic)) * PI);
  longitude = (longitude * 180.0) / ((a / sqrtMagic) * Math.cos(radLat) * PI);
  return { latitude, longitude };
};

/**
 *  将GCJ-02(火星坐标)转为WGS-84:
 */
export const gcj02_to_wgs84 = (lat, lon) => {
  const d = wgs84_to_gcj02(lat, lon);
  return { latitude: lat - d.latitude, longitude: lon - d.longitude };
};

/**
 *  将GCJ-02(火星坐标)转为百度坐标:
 */
export const gcj02_to_baidu = (lat, lon) => {
  if (isLocationOutOfChina(latitude, longitude))
    return { latitude: lat, longitude: lon };
  const z = Math.sqrt(lon * lon + lat * lat) + 0.00002 * Math.sin(lat * pi);
  const theta = Math.atan2(lat, lon) + 0.000003 * Math.cos(lon * pi);
  return {
    latitude: z * Math.sin(theta) + 0.006,
    longitude: z * Math.cos(theta) + 0.0065,
  };
};

/**
 *  将百度坐标转为GCJ-02(火星坐标):
 */
export const baidu_to_gcj02 = (lat, lon) => {
  if (isLocationOutOfChina(latitude, longitude))
    return { latitude: lat, longitude: lon };
  const x = lon - 0.0065;
  const y = lat - 0.006;
  const z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * pi);
  const theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * pi);
  return { latitude: z * Math.sin(theta), longitude: z * Math.cos(theta) };
};

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值