如何用JavaScript实现罗宾逊投影地图像素坐标转经纬度

罗宾逊投影属于伪圆柱投影，不存在通用的解析反解公式，你可以用迭代逼近法实现，3-4次迭代就能达到像素级精度，代码量极小完全符合你小于2KB的体积要求。

实现步骤 #

1. 统一参数基准
   先确定你正向投影的三个固定参数：中央经线（默认通用值为0°，如果你的正向实现自定义了该值，反解时保持一致即可）、投影缩放系数、画布偏移量，反解时先把输入的像素坐标减去偏移量、除以缩放系数，得到和正向输出对齐的归一化投影坐标(x, y)。
2. 迭代求解纬度
   罗宾逊投影的Y值仅和纬度相关，你可以直接复用正向投影里的纬度插值表，用二分法或者牛顿法迭代逼近对应的纬度，迭代3次精度就能到0.0001°级别，完全满足日常使用需求。
3. 计算经度
   得到纬度后，代入正向投影的经度计算公式反向推导，再加上中央经线的偏移值，就能得到最终经度。
4. 适配你的Y轴翻转逻辑
   你之前修改的latSign = 0-getSign(lat)只需要在反解时对Y值先取反再计算即可。

参考实现代码 #

// 罗宾逊投影基础系数表，和正向实现保持一致即可
const ROBINSON_COEFFS = [
  [0, 0.9986, 0],
  [5, 0.9986, 0.062],
  [10, 0.9986, 0.124],
  [15, 0.9986, 0.186],
  [20, 0.9986, 0.248],
  [25, 0.9985, 0.31],
  [30, 0.9983, 0.372],
  [35, 0.998, 0.434],
  [40, 0.9975, 0.4958],
  [45, 0.9969, 0.5571],
  [50, 0.996, 0.6176],
  [55, 0.9949, 0.6769],
  [60, 0.9932, 0.7346],
  [65, 0.9908, 0.7903],
  [70, 0.9874, 0.8435],
  [75, 0.9828, 0.8936],
  [80, 0.976, 0.9394],
  [85, 0.965, 0.9761],
  [90, 0.948, 1]
];
const RAD_TO_DEG = 180 / Math.PI;
const DEG_TO_RAD = Math.PI / 180;
// 你可以根据你的地图尺寸调整以下两个参数，和正向投影保持一致
const CENTRAL_MERIDIAN = 0; // 中央经线，默认0度
const SCALE = 1; // 投影缩放系数

// 正向投影（你已有的实现，这里仅做参数对齐示例）
function robinsonForward(lng, lat) {
  // 替换为你的原有正向逻辑，输出归一化(x,y)
}

// 反向投影实现
function robinsonInverse(x, y) {
  // 适配你之前的Y轴翻转逻辑，根据你的实际情况调整符号
  y = -y;

  // 二分法迭代求纬度
  let low = -90, high = 90, lat = 0;
  for (let i = 0; i < 4; i++) {
    lat = (low + high) / 2;
    const [, , yTest] = robinsonForward(0, lat);
    yTest < y ? low = lat : high = lat;
  }

  // 计算经度
  const [, xCoef] = ROBINSON_COEFFS[Math.floor(Math.abs(lat) / 5)];
  const lng = (x * 180) / (0.8487 * xCoef * SCALE) + CENTRAL_MERIDIAN;

  // 经纬度范围裁剪
  return [
    Math.max(-180, Math.min(180, lng)),
    Math.max(-90, Math.min(90, lat))
  ];
}

上述实现没有任何外部依赖，gzip压缩后体积不到1KB，完全符合你的体积要求。你只需要把robinsonForward替换成你已有的正向实现，调整SCALE和CENTRAL_MERIDIAN和你正向的参数对齐即可使用。

内容的提问来源于stack exchange，提问作者sc0ttj