C#中纬度与归一化墨卡托Y坐标反向转换精度问题（SuperCluster Mapbox）

问题重现 #

实现了LatY（纬度转归一化墨卡托y坐标）和YLat（y坐标转回纬度）函数，但逆转换结果与原纬度偏差极大：

转换代码 #

LatY函数：

double LatY(double lat)
{
    var sin = Math.Sin(lat * Math.PI / 180);
    
    var y = 0.5 - 0.25 * Math.Log((1 + sin) / (1 - sin)) / Math.PI;

    return y;
}

YLat函数：

double YLat(double y)
{
    var exp = -(Math.PI / 0.25) * (y - 0.5);

    var antiLog = Math.Pow(10, exp);

    var sin = (antiLog - 1) / (antiLog + 1);

    var lat2 = Math.Asin(sin) * (180 / Math.PI);
    return lat2;
}

测试案例 #

var y = LatY(-112.484379);
var lat = YLat(y);

预期结果：lat = -112.484379，实际结果：-87.22371494918505

错误原因 #

核心问题是对数与指数函数不匹配：

• LatY中使用的Math.Log是自然对数（以e为底，ln）
• YLat中错误地使用了以10为底的指数函数Math.Pow(10, exp)，而非自然指数函数Math.Exp(exp)

数学上，自然对数ln(x)的反函数是e^x，而非10^x。两者的不匹配导致antiLog的计算完全错误，进而使得sin值偏离真实值，最终arcsin得到的纬度与原数值差距极大。

修复方案 #

将YLat中的Math.Pow(10, exp)替换为Math.Exp(exp)，修正后的函数如下：

double YLat(double y)
{
    var exp = -(Math.PI / 0.25) * (y - 0.5);
    var antiLog = Math.Exp(exp); // 使用自然指数替代10的幂
    var sin = (antiLog - 1) / (antiLog + 1);
    var lat2 = Math.Asin(sin) * (180 / Math.PI);
    return lat2;
}

验证结果 #

修复后执行测试案例，返回的lat值为-112.48437899999999，与原纬度几乎完全一致（仅存在浮点数运算的微小精度误差）。

额外说明 #

墨卡托投影的有效纬度范围是±85.0511287798°（超出此范围时，归一化y坐标会超出0-1区间，且数学上趋向无穷），但修复后的函数在超出该范围时仍能保证双向转换的一致性。

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