瓦片行列坐标转经纬度坐标

瓦片坐标说明：
列号x:对应经度
行号y:对应纬度
层级z:瓦片的等级（LV18/19/20）

两种转换函数：

#这种转成正常的wgs84的经纬度坐标
def XYZToLonLat(x,y,z):
    n = 2**z*1.0
    lon = x / n * 360.0 - 180.0 #经度
    lat = 180 - y / n * 360.0
    return lon,lat

#这种转成谷歌地图上的经纬度坐标
def XYZToLonLat(x,y,z):
    n = 2**z*1.0
    lon = x / n * 360.0 - 180.0 #经度
    lat = math.atan(math.sinh(math.pi * (1 - 2.0 * y / n)))
    lat = math.degrees(lat) #纬度
    return lon,lat

对于整个地球而言，基于 Web 墨卡托投影的地图左上角经纬度坐标为（180°，85.0511 °），右下角经纬度为（-180°，-85.0511°）。纬度范围是[-85.0511, 85.0511 ] 的原因是：保证整个地图为正方形。而谷歌地图采用的就是web墨卡托投影

作者：芒果香蕉_
链接：https://www.jianshu.com/p/8902431e8f46
来源：简书