消失的极点：深度解析 Web 墨卡托投影与高纬度畸变

在之前的文章中，我们讨论了经纬度小数位数的实际影响-CSDN博客
以及深度解析WEBGIS中的“地图级别（Zoom Level）-CSDN博客
。当我们试图将这些数据在浏览器中渲染出来时，不可避免地会遇到制图学的终极难题：如何将一个球形的地球，摊平在一张正方形的屏幕上？

今天，我们就来聊聊现代互联网地图（Google Maps, 高德, Leaflet 等）的通用标准——Web 墨卡托投影（Web Mercator, 代码 EPSG:3857）。它为什么能统治互联网地图？又为何让格陵兰岛看起来和非洲一样大？

墨卡托投影（Mercator Projection）是由地理学家墨卡托在 1569 年提出的。它是一种等角圆柱投影。

想象一个透明的地球，在其内部中心放一盏灯。我们将一张纸卷成圆柱体套在地球外面，纸面与赤道相切。当地球上的地理信息被光影投射到这张纸上并展开，就得到了墨卡托投影图。

墨卡托投影最伟大的贡献是保持了等角性。这意味着在地图上任何一点，经线和纬线都是垂直相交的。对于航海者来说，这极其重要：如果你在地图上画一条直线，这条线与经线的夹角（方位角）在现实中是恒定的。

传统的墨卡托投影基于椭球体计算，公式复杂。为了提高互联网时代的渲染效率，Google Maps 在 2005 年对其进行了简化，形成了Web 墨卡托：

球体代替椭球体：它在投影计算时将地球视为完美的球体（R = 6378137 R=6378137R=6378137米），极大地降低了 CPU 的三角函数运算量。
投影成正方形：为了完美适配“瓦片金字塔”体系（即 Zoom 0 是一张正方形图片），Web 墨卡托对纬度范围进行了强制裁剪。

为了让地图的长宽相等，纬度的最大范围被限制在：

± 85.051129 ∘ \pm 85.051129^\circ±85.051129∘

这就是为什么你在大多数互联网地图上看不到南极点和北极点的原因——它们在数学上被投影到了无限远处。

墨卡托投影最被诟病的缺点就是面积畸变。这种变形源于它为了保持“角度不变”而付出的代价。

在球体上，经线在两极汇聚。但在墨卡托投影中，经线被强行平行化了。这意味着随着纬度升高，水平方向被强行拉伸了。在纬度ϕ \phiϕ处，水平拉伸倍数是1 / cos ⁡ ( ϕ ) 1/\cos(\phi)1/cos(ϕ)。

为了保持角度不变（即局部形状不失真），既然水平方向拉伸了，垂直方向也必须同步拉伸。

公式表达为：

y = ln ⁡ [ tan ⁡ ( π 4 + ϕ 2 ) ] y = \ln[\tan(\frac{\pi}{4} + \frac{\phi}{2})]y=ln[tan(4π+2ϕ)]

在赤道（ϕ = 0 ∘ \phi=0^\circϕ=0∘），cos ⁡ ( 0 ) = 1 \cos(0)=1cos(0)=1，没有任何变形。
在纬度60 ∘ 60^\circ60∘（如北欧），cos ⁡ ( 60 ∘ ) = 0.5 \cos(60^\circ)=0.5cos(60∘)=0.5，拉伸倍数是1 / 0.5 = 2 1/0.5 = 21/0.5=2倍。这意味着面积被放大了2 × 2 = 4 2 \times 2 = 42×2=4倍。
在纬度80 ∘ 80^\circ80∘附近，面积缩放比例达到了惊人的33 倍。

这就是著名的“格陵兰岛错觉”：在 Web 墨卡托地图上，格陵兰岛看起来和非洲大陆（实际面积是其 14 倍）差不多大。

尽管面积畸变严重，但 Web 墨卡托在互联网地图中地位不可动摇，原因有三：

作为 GIS 开发者或数据分析师，必须警惕投影带来的误导：

不要在 Web 墨卡托（EPSG:3857）下直接计算面积：如果你计算俄罗斯的面积，结果会比实际大得多。建议：先转换到等面积投影（如 EPSG:4490 或 WGS84 投影），或使用专门的测地线面积计算函数。
距离测量需修正：我们在上一篇提到的“地面分辨率”随纬度变化，正是因为投影拉伸。计算距离时应使用大圆距离（Haversine 公式）。
3D 地图的回归：近年来，Google Maps 和 Mapbox 在低级别（大尺度）下开始切换到球体投影（Globe View），只有放大到城市级别才切换回墨卡托，这完美解决了面积错觉问题。