如何在墨卡托地图（JPEG）上从x，y获取纬度和经度？

Question

我有一个墨卡托投影映射为JPEG，我想知道如何将给定的x，y坐标与其经度和纬度关联。我看过Gudermannian函数，但我真的不明白如何使用该函数并应用它。也就是说，它期望什么输入？我发现的实现（JavaScript）似乎需要-PI和PI之间的范围，但是我的y值（以像素为单位）和该范围之间的相关性如何？

此外，我发现这个函数需要一个纬度并返回谷歌地图，也使用墨卡托瓷砖。似乎如果我知道如何反转这个功能，我会非常接近我的答案。

/*<summary>Get the vertical tile number from a latitude 
using Mercator projection formula</summary>*/ 

    private int getMercatorLatitude(double lati) 
    { 
     double maxlat = Math.PI; 

     double lat = lati; 

     if (lat > 90) lat = lat - 180; 
     if (lat < -90) lat = lat + 180; 

     // conversion degre=>radians 
     double phi = Math.PI * lat/180; 

     double res; 
     //double temp = Math.Tan(Math.PI/4 - phi/2); 
     //res = Math.Log(temp); 
     res = 0.5 * Math.Log((1 + Math.Sin(phi))/(1 - Math.Sin(phi))); 
     double maxTileY = Math.Pow(2, zoom); 
     int result = (int)(((1 - res/maxlat)/2) * (maxTileY)); 

     return (result); 
    } 

Answer 1

这里有一些代码给你...让我知道如果你需要更多的解释。

/// <summary> 
    /// Calculates the Y-value (inverse Gudermannian function) for a latitude. 
    /// <para><see cref="http://en.wikipedia.org/wiki/Gudermannian_function"/></para> 
    /// </summary> 
    /// <param name="latitude">The latitude in degrees to use for calculating the Y-value.</param> 
    /// <returns>The Y-value for the given latitude.</returns> 
    public static double GudermannianInv(double latitude) 
    { 
     double sign = Math.Sign(latitude); 
     double sin = Math.Sin(latitude * RADIANS_PER_DEGREE * sign); 
     return sign * (Math.Log((1.0 + sin)/(1.0 - sin))/2.0); 
    } 

    /// <summary> 
    /// Returns the Latitude in degrees for a given Y. 
    /// </summary> 
    /// <param name="y">Y is in the range of +PI to -PI.</param> 
    /// <returns>Latitude in degrees.</returns> 
    public static double Gudermannian(double y) 
    { 
     return Math.Atan(Math.Sinh(y)) * DEGREES_PER_RADIAN; 
    } 

Answer 2

谷歌等等，使用“球墨卡托”，即墨卡托投影使用球形地球模型而不是更慢和更复杂的椭圆方程。

中的变换是可作为的的OpenLayers代码的一部分：

http://docs.openlayers.org/library/spherical_mercator.html

Answer 3

我已经做了类似的东西。特别是如果你有来自世界的一部分的图像。裁切过的地图，或只是不是一个完整的世界地图：https://stackoverflow.com/a/10401734/730823

Answer 4

进行逆时一个重要的注意的是，有作为“的墨卡托地图”没有这样的事情是与其他大多数地图投影的情况。存在的每个墨卡托地图根据输入的phi值而不同。根据维基百科，谷歌使用85.051129，其他地图提供商使用85.05113。因此，Gudermannian的输入值必须根据例如GudermannianInv（85.05113）。

Answer 5

Erich Mirabal's answer是完全正确的（如果不是完全的话）。

我刚刚使用'理论上的256x256墨卡托瓷砖'（谷歌单一瓷砖版本的世界地图）对它进行了测试。

这里有一个小更多的代码（JavaScript的，但易于遵循）澄清。

我住在澳大利亚，纬度约为-33°。

convertRange(
    GudermannianInv(-33), 
    [Math.PI, - Math.PI], 
    [0, 256] 
); 

152.88327883810192

如果算上152个像素从瓦的顶部向下，你会发现澳大利亚。 我也通过将结果与已知的良好函数进行比较来验证此答案是正确的。

可以肯定的，我们可以扭转计算：

Gudermannian(
    convertRange(
     152.88, 
     [0, 256], 
     [Math.PI, - Math.PI] 
)); 

现在我们又回到-32.99613291758226。

棘手的部分不处于古德曼函数，但在两个尺度之间的转换。

幸运的是，是相当懒了，恨这些类型的结垢问题，我已经有了一个小功能来做到这一点凌乱的转换我。

/** 
    * convert number from _n_ of r1[0] .. r1[1] to _n_ of r2[0] .. r2[1] 
    * @example `convertRange(5, [0, 10], [0, 100]) === 50` 
    * 
    * @param {number} value 
    * @param {array<number>} r1 old range 
    * @param {array<number>} r2 new range 
    * @returns {number} value adjusted for new range 
    */ 
    function convertRange(value, r1, r2) { 
     return (value - r1[0]) 
      * (r2[1] - r2[0]) 
      /(r1[1] - r1[0]) 
      + r2[0]; 
    } 

和JavaScript版本的原有功能自然是：

function Gudermannian(y) { 
    return Math.atan(Math.sinh(y)) * (180/Math.PI) 
} 

function GudermannianInv(latitude) 
{ 
    var sign = Math.sign(latitude); 
    var sin = Math.sin(
          latitude 
         * (Math.PI/180) 
         * sign 
    ); 
    return sign * (
     Math.log(
      (1 + sin)/(1 - sin) 
     )/2 
    ); 
}