如何使Cartopy墨卡托投影地图保持1km×1km的等纵横比？

嘿，这个问题我之前帮人解决过好几次——Mercator投影默认是按经纬度（度数）来处理轴纵横比的，所以才会出现看起来“扁”或者“拉伸”的情况，尤其是在高纬度地区。要实现1km×1km的等比例显示，核心就是让轴的像素比例对应实际地面的米制距离，下面给你两个靠谱的方案：

方案1：手动计算纬度对应的米/度转换率 #

Mercator投影中，经度方向的每度距离会随纬度变化（纬度越高，每度经度对应的地面距离越短），而纬度方向的每度距离几乎固定在~111.32km。我们可以基于你关注区域的中心纬度，计算出准确的纵横比：

代码示例 #

import cartopy.crs as ccrs
import matplotlib.pyplot as plt
import numpy as np

# 初始化Mercator投影和绘图轴
proj = ccrs.Mercator()
fig, ax = plt.subplots(subplot_kw={'projection': proj})

# 设置你的目标经纬度范围（WGS84坐标系）
lon_min, lon_max = 110.0, 111.0
lat_min, lat_max = 30.0, 31.0
ax.set_extent([lon_min, lon_max, lat_min, lat_max], crs=ccrs.PlateCarree())

# 计算区域中心纬度
lat_center = (lat_min + lat_max) / 2
# 计算该纬度下，每度经度对应的米数（近似公式，足够日常使用）
meters_per_degree_lon = 111320 * np.cos(np.deg2rad(lat_center))
# 纬度方向每度固定为~111320米
meters_per_degree_lat = 111320

# 计算纵横比：确保图上1米的经度距离 = 1米的纬度距离
# matplotlib的aspect是height/width，所以要换算成(y方向总米数)/(x方向总米数)
aspect_ratio = ((lat_max - lat_min) * meters_per_degree_lat) / ((lon_max - lon_min) * meters_per_degree_lon)
ax.set_aspect(aspect_ratio, adjustable='box')

# 示例：绘制海岸线验证效果
ax.coastlines(resolution='10m')
plt.show()

方案2：直接用投影转换后的米制范围计算（更精确） #

如果你想避免手动计算的近似值，可以让Cartopy直接把经纬度范围转换成Mercator投影的米制坐标，然后基于实际的米数范围计算纵横比：

代码示例 #

import cartopy.crs as ccrs
import matplotlib.pyplot as plt

# 初始化投影和轴
proj = ccrs.Mercator()
fig, ax = plt.subplots(subplot_kw={'projection': proj})

# 设置经纬度范围
lon_min, lon_max = 110.0, 111.0
lat_min, lat_max = 30.0, 31.0
ax.set_extent([lon_min, lon_max, lat_min, lat_max], crs=ccrs.PlateCarree())

# 获取Mercator投影下的米制坐标范围
x_min, x_max, y_min, y_max = ax.get_extent(proj)

# 计算纵横比：保证图上x方向1米 = y方向1米
aspect_ratio = (y_max - y_min) / (x_max - x_min)
ax.set_aspect(aspect_ratio, adjustable='box')

# 绘制验证
ax.coastlines(resolution='10m')
plt.show()

关键注意点 #

• adjustable='box'：保持轴的位置不变，调整数据显示范围来匹配纵横比；如果希望轴的大小自动适配，可以改成adjustable='datalim'。
• 高纬度地区：Mercator投影在高纬度会严重拉伸，所以一定要基于你实际关注的区域计算纵横比，不要用全局平均值。
• 交互式调整：如果在Jupyter等交互式环境中动态调整地图范围，建议把纵横比计算逻辑封装成函数，每次更新范围后调用。

内容的提问来源于stack exchange，提问作者ChrisB