WGS转墨卡托坐标转换函数输出不符预期的修正方法

问题本质 #

你当前使用的是标准正轴等角墨卡托投影的归一化Y计算公式，这个投影的固有特性是高纬度（两极附近）区域被拉伸，和你需要的「两极纬度分段更窄、赤道纬度分段更宽」的需求正好相反，才会出现45°纬度位置和预期对齐位置偏差的问题。

现有公式偏差验证 #

你当前的代码实现：

WgsToMercator(coord)
{
           yRadian = coord.y * Math.PI / 180.0;
           sinLat = Math.Sin(yRadian);
           y = 0.5 - Math.Log((1 + sinLat) / (1 - sinLat)) / (Math.PI * 4); 
           return y;
}

代入45°纬度计算时，sin(45°)≈0.707，对数项计算结果约为0.14，最终输出y≈0.36，45°位置距离赤道（y=0.5）仅0.14单位，距离北极（y=0）有0.36单位，属于典型的高纬度拉伸效果，和你的预期正好反向。

调整后实现方案 #

不需要保留原公式的对数变换，直接基于等距圆柱投影做幂次曲率调整，就能完全匹配你给出的预期输出，同时满足赤道宽、两极窄的分段要求：

WgsToMercator(coord)
{
    // 纬度转弧度，北纬为正，范围[-π/2, π/2]
    double yRadian = coord.y * Math.PI / 180.0;
    // 纬度线性归一化到[-0.5, 0.5]区间
    double normLinearLat = yRadian / Math.PI;
    // 曲率参数，取0.47时可精准匹配45°对应0.14/0.86的校准要求
    const double curvature = 0.47;
    // 幂次调整曲线，实现赤道拉伸、两极压缩
    double adjustedOffset = Math.Sign(normLinearLat) * 0.5 * Math.Pow(Math.Abs(normLinearLat) * 2, curvature);
    // 映射到[0,1]归一化区间
    double y = 0.5 - adjustedOffset;
    return y;
}

输出结果验证 #

代入测试值计算结果完全匹配预期：

• WgsToMercator(90) = 0
• WgsToMercator(45) ≈ 0.14
• WgsToMercator(0) = 0.5
• WgsToMercator(-45) ≈ 0.86
• WgsToMercator(-90) = 1

参数微调说明 #

如果后续渲染对齐存在微小偏差，可以直接调整curvature参数适配：

• 参数值越小，两极区域分段越窄、赤道区域分段越宽
• 参数值=1时为标准等距圆柱投影，所有纬度分段宽度一致
• 参数值>1时会呈现和原墨卡托公式类似的两极拉伸效果

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