tile2lon：地图瓦片编号与经纬度的换算关系

在生成矢量切片之前，看到一个换算关系。。。就是切片的经纬度12/3000/4000之类的换算成12/120.111/30.222之类的。

https://www.jianshu.com/p/8902431e8f46

地图瓦片分割：

 

 目前所接触的瓦片都是从左上角开始进行编号的。。

左上角的经纬度是（-180°，85°），右下角的经纬度是（180°，-85°）

纬度的范围之所以是（-85°，85°）的原因是：为了保证地图是正方形。

 公式的具体推导过程见：https://en.wikipedia.org/wiki/Mercator_projection

 

整体是一个正方形，但是具体到每一小块，并不是正方形啊。。

所以是怎么切割成正方形的。。

在ArcGIS JS API 加载 TMS 地图瓦片 这篇笔记中提到过缩放等级 z 和每行（或每列）瓦片数量 n 的关系如下：

由上述可知投影后地图经度范围是[-180, 180]，在第 z 级别每行的瓦片数为 n。那么等级 z 下某一经度对应的 x 轴编号为：

而等级 z 下某一纬度对应的 y 轴编号则比较复杂（因为纬度分布不均匀）：

已知瓦片编号反算该瓦片左上角经纬度坐标公式如下：