GeoServer生成EPSG:3857切片转经纬度后纬度为负的原因

我用GeoServer生成了EPSG:3857投影的切片，希望通过切片文件名中的x、y、z参数转换得到经纬度，使用了OpenStreetMap相关的JavaScript转换函数：

function tile2long(x,z) {
  return (x/Math.pow(2,z)*360-180);
}
function tile2lat(y,z) {
  var n=Math.PI-2*Math.PI*y/Math.pow(2,z);
  return (180/Math.PI*Math.atan(0.5*(Math.exp(n)-Math.exp(-n))));
}

转换后的经度符合预期，但纬度始终为负值，请问这是什么原因？

原因分析 #

核心问题是切片坐标系的y轴原点定义不一致：

• OpenStreetMap的切片规范中，y轴原点位于地图左上角（对应北纬约85°），y值从0开始向下递增；
• GeoServer生成EPSG:3857切片时，默认y轴原点位于地图左下角（对应南纬约85°），y值从0开始向上递增。

由于两者y轴方向完全相反，直接用OSM的公式代入GeoServer的y值，计算出的纬度自然会呈现负值。

解决方法 #

只需要将GeoServer的y值转换为符合OSM规范的y值，再代入原公式即可：

// 转换GeoServer的y值为OSM规范的y值
function tile2lat_geoserver(y, z) {
  const osmY = Math.pow(2, z) - 1 - y;
  const n = Math.PI - 2 * Math.PI * osmY / Math.pow(2, z);
  return (180 / Math.PI) * Math.atan(0.5 * (Math.exp(n) - Math.exp(-n)));
}

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