从切片编号派生切片坐标

Question

我觉得这个问题以前已经回答过了，但不知何故我偏离了轨道。我有一个磁贴的URL (见下文)：我看到缩放级别是6，33个磁贴朝向东方，20个磁贴朝向南方。我的问题是，我能不能用上述数字(33 / 20)来推断这个瓦片的经纬度(web墨卡托)坐标？

https://heatmap-external-c.strava.com/tiles/all/hot/6/33/20@2x.png?v=19

​

​

Answer 1

正如问题的标签所述，链接的地图是通过mapbox GL呈现的。在mapbox docs中，我找到了对OSM Slippy Map Tilenames specification的引用，它提供了将tilename转换为后来的长坐标的方法，反之亦然。The following pseudocode回答了这个问题。

对于给定缩放级别(zoom)，将编号(xtile、ytile)平铺到经度/纬度

n = 2 ^ zoom
lon_deg = xtile / n * 360.0 - 180.0
lat_rad = arctan(sinh(π * (1 - 2 * ytile / n)))
lat_deg = lat_rad * 180.0 / π

顺便说一句:这种类型的问题在gis.stackexchange.com上更常见