C++如何实现一个简单的2D地图瓦片坐标转换算法-C++-PHP中文网

瓦片坐标系基于web mercator投影，非线性映射且y轴反向；不能直接用像素除以瓦片尺寸计算，须先将经纬度转为web mercator米制坐标，再归一化并缩放至z级网格。

c++如何实现一个简单的2d地图瓦片坐标转换算法

什么是瓦片坐标系，为什么不能直接用像素除以瓦片大小

瓦片坐标（如 z/x/y）不是线性映射到经纬度的——它基于 Web Mercator 投影，且 y 轴方向与常规图像坐标相反（地图顶部对应 y=0，向下递增）。直接用像素位置除以瓦片尺寸（比如 256）会得到错误的 x、y，尤其在高纬度或跨赤道区域偏差极大。

真正可用的转换必须分两步：先将地理坐标（经纬度）转为 Web Mercator 米制坐标，再按当前缩放级 z 划分网格并取整。关键函数是 lonLatToMeters() 和 metersToTile()，不是简单四则运算。

C++里怎么把经纬度转成瓦片坐标（z/x/y）

核心逻辑是：经度 → Web Mercator X，纬度 → Web Mercator Y，再归一化到 [0,1] 区间，最后乘以 1 得整数瓦片索引。注意 <code>y 需用 1 - normalized_y 翻转。

实操建议：

std::atanh 和 std::sin 必须用弧度输入，别忘了 deg2rad 转换
纬度范围必须 clamp 在 (-85.0511, 85.0511)，否则 atanh(tan(...)) 会返回 inf 或 nan
瓦片 x 和 y 要用 floor 后强转 int，不是 round——瓦片编号永远向下取整
缩放级 z 通常为 0–19，超出时需提前截断，避免 1 溢出（<code>z > 30 在 int 下危险）

简短示例：

double deg2rad(double deg) { return deg * M_PI / 180.0; }
double lonLatToX(double lon) { return deg2rad(lon) * 20037508.34; }
double lonLatToY(double lat) {
    double latrad = deg2rad(lat);
    return std::log(std::tan(M_PI/4 + latrad/2)) * 20037508.34;
}
std::tuple<int int> latLonToTile(double lon, double lat, int z) {
    double mx = lonLatToX(lon);
    double my = lonLatToY(std::clamp(lat, -85.0511, 85.0511));
    double res = 40075016.6856 / (1 (std::floor((mx + 20037508.34) / res));
    int y = static_cast<int>(std::floor((20037508.34 - my) / res));
    return {x, y};
}
</int></int>