深入了解XYZ瓦片:在线地图数据格式及实现原理

XYZ瓦片是一种在线地图数据格式,常见的地图底图如Google、OpenStreetMap等互联网的瓦片地图服务,都是XYZ瓦片,严格来说是ZXY规范的地图瓦片。

ZXY规范的地图瓦片规则如下:将地图全幅显示时的图片从左上角开始,往下和往右进行切割,切割的大小默认为 256*256 像素,左上角的格网行号为 0,列号为 0,往下和往右依次递增。

图片来源: ZXY标准瓦片 (supermap.com.cn)

从整体来说,XYZ瓦片数据结构是一种影像金字塔。

图片来源: 瓦片底图:在线地图的下载和使用 | Mars3D开发教程 ,此图仅供参考,图中的行列号和ZXY规范的地图瓦片的行列号编码存在差异。

XYZ瓦片与经纬度的计算以及原理

首先给出经纬度与XYZ行列号之间的计算公式。

图片来源: Slippy map tilenames - OpenStreetMap Wiki

现在解释一下原理

下面是一张OpenStreetMap在zoom等级为2时的瓦片示意图。

图片来源: Slippy map tilenames - OpenStreetMap Wiki

z 是当前的瓦片等级,就是缩放等级,由上面的图可以看出:z 等级时,共有 \(2^z\) 个瓦片,x范围为0-\(2^z-1\),y范围也是0-\(2^z-1\)。

首先 x 的计算很简单。

目的:将经度从-180度到180度,映射到0到\(2^z\)之间的整数列号上。

过程:先将经度加180度,使其从0到360度,然后除以360(归一化)再乘以\(2^z\)得到行号,最后向下取整数部分,得到最终的行号。

y 的计算就复杂多了。

目的:将纬度从-90度到90度,映射到0到\(2^z\)之间的整数行号上。

存在的问题:纬度分布不均匀,XYZ瓦片试图将地图展开为一个正方形(参考上图,本质上就是Web墨卡托投影),然而纬度是中间(赤道)长两极短,如果只是像 x 一样简单的映射,会导致两极的紧凑,赤道附近稀疏。

解决方案:将纬度通过一种映射,使其能均匀一点,然后就采用了下面的函数。

图片来源: Slippy map tilenames - OpenStreetMap Wiki

在浏览器端实现XYZ瓦片的加载示例

计算公式实现

根据上面的公式,很容易就把根据经纬度算行列号的函数写出来

function lon2tile(lon, zoom) { 
    return (Math.floor((lon + 180) / 360 * Math.pow(2, zoom)))
}
function lat2tile(lat, zoom) { 
    return (Math.floor((1 - Math.log(Math.tan(lat * Math.PI / 180) + 1 / Math.cos(lat * Math.PI / 180)) / Math.PI) / 2 * Math.pow(2, zoom)))
}

事实上,这个网站已经给出了这个公式的各种编程语言的实现: Slippy map tilenames - OpenStreetMap Wiki

核心代码

function lon2tile(lon, zoom) {
    return (Math.floor((lon + 180) / 360 * Math.pow(2, zoom)));
}
function lat2tile(lat, zoom) {
    return (Math.floor((1 - Math.log(Math.tan(lat * Math.PI / 180) + 1 / Math.cos(lat * Math.PI / 180)) / Math.PI) / 2 * Math.pow(2, zoom)));
}
const loadMapByBounds = (minLon, minLat, maxLon, maxLat, zoom) => {
    const minTileX = lon2tile(minLon, zoom);
    const minTileY = lat2tile(maxLat, zoom); // Y轴是反的,自上而下
    const maxTileX = lon2tile(maxLon, zoom);
    const maxTileY = lat2tile(minLat, zoom);
    for (let x = minTileX; x <= maxTileX; x++) {
        for (let y = minTileY; y <= maxTileY; y++) {
            loadTile(x, y, zoom); // 加载瓦片
        }
    }
}

完整实现

为了简单,这里使用img标签来加载瓦片图,并根据瓦片编号排列,设置对应的偏移值。

为了能拖动以浏览全图实现简单的交互,这里还设置了根据鼠标按压后拖动的偏移值来添加对应的偏移值。

实现效果如下:

图片来源: 瓦片底图:在线地图的下载和使用 | Mars3D开发教程

完整代码如下:

<!DOCTYPE html>
<html lang="en">
<head>
    <meta charset="UTF-8">
    <meta name="viewport" content="width=device-width, initial-scale=1.0">
    <title>Document</title>
    <style>
        img {
            width: 256px;
            height: 256px;
        }
        html,
        body {
            margin: 0;
            padding: 0;
            overflow: hidden;
            height: 100%;
            width: 100%;
        }
        #map {
            position: absolute;
            height: 100%;
            width: 100%;
            overflow: hidden;
            border: 1px solid #000;
        }
    </style>
</head>
<body>
    <div id="map"></div>
    <script>
        // 此处为代码内容,详见原文
    </script>
</body>
</html>

参考资料

[1] Slippy map tilenames - OpenStreetMap Wiki

[2] ZXY标准瓦片 (supermap.com.cn)

[3] 瓦片底图:在线地图的下载和使用 | Mars3D开发教程

标签:游戏攻略