瓦片坐标系基于Web Mercator投影，非线性映射且y轴反向；不能直接用像素除以瓦片尺寸计算，须先将经纬度转为Web Mercator米制坐标，再归一化并缩放至z级网格。
什么是瓦片坐标系，为什么不能直接用像素除以瓦片大小
瓦片坐标（如 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 翻转。
实操建议：
立即学习“C++免费学习笔记（深入）”；
-
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下危险）
简短示例：
组合式C++代码评审方案，融合静态分析、AI推理、多轮迭代评审和C++专项检查，适用于PR审查、增量代码审查、全项目评审和代码质量评分，触发词包括review cpp、cpp代码评审、C++review、代码审查。
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 << z); // 米/像素
int x = static_cast<int>(std::floor((mx + 20037508.34) / res));
int y = static_cast<int>(std::floor((20037508.34 - my) / res));
return {x, y};
}
已知瓦片坐标（z/x/y），怎么算出它左上角的经纬度
这是上一步的逆过程：先还原瓦片左上角的 Web Mercator 米坐标，再反解经纬度。关键点在于 y 要用 1 - y_norm，否则算出来是右下角或完全错位。
常见错误现象：
- 返回的纬度全是 NaN：没检查
y是否越界（y < 0 || y >= (1<<z)） - 结果偏移一整个瓦片：用了
y+1当作左上，实际应直接用y计算左上角 - 经度超 [-180,180]：反解后没做
fmod(lon + 180, 360) - 180归一化
参数差异：输入 y 是整数瓦片号，但计算时要转成归一化比例 y / (1.0 ，再映射回米值。
性能和兼容性要注意什么
Web Mercator 的 atanh(tan(...)) 计算开销不小，高频调用（如渲染循环）建议预计算缩放级对应的分辨率 res 并缓存；z 不变时，1 可提至外层避免重复位运算。
跨平台细节：
-
M_PI在某些标准（如 C++11 默认）下不定义，需加#define _USE_MATH_DEFINES或手动定义 -
std::clamp是 C++17 起支持，老项目用三元表达式替代：(lat < -85.0511 ? -85.0511 : (lat > 85.0511 ? 85.0511 : lat)) - Windows 上
std::log对负零可能返回 -inf，传入前确保my不是 NaN/inf（可用std::isfinite检查）
最易被忽略的是：瓦片坐标系中同一地理点在不同 z 下的 x/y 值不满足简单倍数关系——因为 Web Mercator 的非线性拉伸随纬度变化，不能靠插值或缩放推导，每次都要完整重算。
