在 Python 绘图（尤其是地理信息可视化）中添加指南针（指北针）通常涉及使用支持地图投影的库，如 `matplotlib` 的 `Basemap` 工具包或其继任者 `Cartopy`。由于 `Basemap` 已停止维护，现代实践更推荐使用 `Cartopy` [ref_5]。以下将分别介绍使用 `Cartopy` 和 `Basemap`（用于遗留代码参考）添加指南针的方法，并提供完整的实现代码。 ### 1. 使用 Cartopy 添加指南针 (推荐方法) `Cartopy` 是一个强大的地理空间数据处理库，与 `matplotlib` 深度集成。它提供了 `gridliner` 和 `add_compass`（通过 `geoaxes` 的扩展）等功能来添加指北针和刻度。 #### 方法一：使用 `gridliner` 显示指北针标签 `Cartopy` 的 `gridliner` 可以在地图边缘绘制经纬网格线，并通过设置 `xlocs` 和 `ylocs` 来控制标签的显示。虽然不直接绘制一个图形化的指南针，但通过正确标注经纬度，可以明确指示方向（通常地图上北下南）。 ```python import matplotlib.pyplot as plt import cartopy.crs as ccrs import cartopy.feature as cfeature # 创建图形和地理坐标轴，使用 PlateCarree 投影（最常见的经纬度投影） fig, ax = plt.subplots(figsize=(10, 8), subplot_kw={'projection': ccrs.PlateCarree()}) # 设置地图范围（西经、东经、南纬、北纬） ax.set_extent([-20, 60, -40, 40], crs=ccrs.PlateCarree()) # 添加地理特征：海岸线、陆地、海洋 ax.add_feature(cfeature.COASTLINE, linewidth=0.8) ax.add_feature(cfeature.LAND, facecolor='lightgray') ax.add_feature(cfeature.OCEAN, facecolor='lightblue') # 添加经纬网格线，并启用标签 gl = ax.gridlines(draw_labels=True, linewidth=1, color='gray', alpha=0.5, linestyle='--') # 自定义标签位置：顶部和右侧不显示标签，左侧和底部显示 gl.top_labels = False gl.right_labels = False gl.xlabel_style = {'size': 10, 'color': 'black'} gl.ylabel_style = {'size': 10, 'color': 'black'} ax.set_title('使用 Cartopy 绘制地图并显示经纬网格（指示方向）', fontsize=14, pad=20) plt.tight_layout() plt.show() ``` #### 方法二：自定义绘制图形化指南针 对于需要传统图形化指南针（如箭头和“N”标识）的场景，可以结合 `matplotlib` 的绘图功能在 `Cartopy` 地图的角落手动绘制一个。 ```python import matplotlib.pyplot as plt import cartopy.crs as ccrs import cartopy.feature as cfeature import numpy as np fig, ax = plt.subplots(figsize=(10, 8), subplot_kw={'projection': ccrs.PlateCarree()}) ax.set_extent([-20, 60, -40, 40], crs=ccrs.PlateCarree()) ax.add_feature(cfeature.COASTLINE) ax.add_feature(cfeature.LAND, facecolor='lightgray') ax.add_feature(cfeature.OCEAN, facecolor='lightblue') # 添加网格线（不带标签） gl = ax.gridlines(draw_labels=False, linewidth=0.5, color='gray', alpha=0.5, linestyle=':') # --- 自定义绘制指南针 --- # 定义指南针在图形坐标系中的位置和大小（左下角） compass_x, compass_y = 0.05, 0.05 # 图形坐标系下的位置 (0,0是左下角，1,1是右上角) compass_radius = 0.04 # 指南针半径（图形坐标系单位） # 在图形坐标系中添加一个子坐标轴用于绘制指南针 compass_ax = fig.add_axes([compass_x, compass_y, compass_radius*2, compass_radius*2]) compass_ax.set_aspect('equal') compass_ax.axis('off') # 隐藏坐标轴 # 绘制指南针圆盘 circle = plt.Circle((0.5, 0.5), 0.4, transform=compass_ax.transAxes, edgecolor='black', facecolor='white', linewidth=2) compass_ax.add_patch(circle) # 绘制指北针箭头（红色） # 使用箭头标注功能，调整起点和终点在坐标轴坐标系中的位置 compass_ax.annotate('', xy=(0.5, 0.9), xytext=(0.5, 0.5), xycoords='axes fraction', textcoords='axes fraction', arrowprops=dict(arrowstyle='->', color='red', lw=3, shrinkA=0, shrinkB=0)) # 添加“N”标签 compass_ax.text(0.5, 0.95, 'N', transform=compass_ax.transAxes, ha='center', va='center', fontsize=12, fontweight='bold', color='red') # 可选：添加其他方向标识（E, S, W） compass_ax.text(0.9, 0.5, 'E', transform=compass_ax.transAxes, ha='center', va='center', fontsize=10, color='black') compass_ax.text(0.5, 0.05, 'S', transform=compass_ax.transAxes, ha='center', va='center', fontsize=10, color='black') compass_ax.text(0.1, 0.5, 'W', transform=compass_ax.transAxes, ha='center', va='center', fontsize=10, color='black') ax.set_title('使用 Cartopy 绘制地图并添加自定义图形指南针', fontsize=14, pad=20) plt.tight_layout() plt.show() ``` ### 2. 使用 Basemap 添加指南针 (传统方法，已弃用) `Basemap` 提供了一个便捷的 `drawcompass()` 或 `drawmeridians()`/`drawparallels()` 结合自定义函数的方法来添加指南针 [ref_5]。但请注意，`Basemap` 自 2020 年起已不再维护，建议新项目使用 `Cartopy`。 #### 使用 `drawcompass()` 方法（如果可用） 在某些版本的 `Basemap` 中，可能存在 `drawcompass()` 方法，但官方文档并未广泛记载，且并非所有版本都支持。 #### 自定义绘制指南针函数（通用方法） 更可靠的方法是自定义一个绘制指南针的函数，并将其放置在地图的合适位置 [ref_5]。 ```python import matplotlib.pyplot as plt from mpl_toolkits.basemap import Basemap import numpy as np # 初始化地图 fig, ax = plt.subplots(figsize=(10, 8)) # 创建 Basemap 实例，使用等距圆柱投影 m = Basemap(projection='cyl', llcrnrlon=-20, llcrnrlat=-40, urcrnrlon=60, urcrnrlat=40, resolution='c', ax=ax) m.drawcoastlines(linewidth=0.8) m.fillcontinents(color='lightgray', lake_color='lightblue') m.drawmapboundary(fill_color='lightblue') # 绘制经纬网格线 m.drawmeridians(np.arange(-20, 61, 20), labels=[0,0,0,1], fontsize=10, linewidth=0.5, dashes=[1,0]) m.drawparallels(np.arange(-40, 41, 20), labels=[1,0,0,0], fontsize=10, linewidth=0.5, dashes=[1,0]) # --- 自定义函数：在地图角落添加指南针 --- def add_compass(m, x, y, size=0.05): """ 在地图上添加一个简单的指南针。 参数: m: Basemap 实例。 x, y: 指南针中心点在地图坐标系（经纬度）中的位置。 size: 指南针的大小（度）。 """ # 计算指南针四个方向的端点坐标 # 北方向 north_x, north_y = m(x, y + size) # 南方向 south_x, south_y = m(x, y - size) # 东方向 east_x, east_y = m(x + size, y) # 西方向 west_x, west_y = m(x - size, y) center_x, center_y = m(x, y) # 绘制十字线和圆圈 ax.plot([center_x, north_x], [center_y, north_y], 'r-', lw=2) # 北向红线 ax.plot([center_x, south_x], [center_y, south_y], 'k-', lw=1) # 南向黑线 ax.plot([center_x, east_x], [center_y, east_y], 'k-', lw=1) # 东向黑线 ax.plot([center_x, west_x], [center_y, west_y], 'k-', lw=1) # 西向黑线 # 绘制中心圆点 ax.scatter(center_x, center_y, color='black', s=30, zorder=5) # 添加方向标签 ax.text(north_x, north_y, 'N', fontsize=12, fontweight='bold', color='red', ha='center', va='bottom') ax.text(south_x, south_y, 'S', fontsize=10, ha='center', va='top') ax.text(east_x, east_y, 'E', fontsize=10, ha='left', va='center') ax.text(west_x, west_y, 'W', fontsize=10, ha='right', va='center') # 在地图的右下角添加指南针（例如：经度50，纬度-35） add_compass(m, x=50, y=-35, size=5) ax.set_title('使用 Basemap 绘制地图并添加自定义指南针', fontsize=14, pad=20) plt.tight_layout() plt.show() ``` ### 3. 方法对比与选择指南 下表对比了上述几种实现方式的主要特点和适用场景 [ref_5]： | 方法 | 所属库 | 核心原理 | 优点 | 缺点 | 适用场景 | | :--- | :--- | :--- | :--- | :--- | :--- | | **Cartopy 网格标签** | Cartopy | 通过绘制并标注经纬网格来指示方向。 | 1. 现代、活跃维护的库。<br>2. 与 matplotlib 集成好，代码简洁。<br>3. 方向指示符合专业地图规范。 | 1. 非传统图形化指南针样式。<br>2. 网格标签可能在地图边缘，不够突出。 | 需要标准、专业的地理可视化，方向由经纬度网格隐含指示。 | | **Cartopy 自定义图形** | Cartopy + Matplotlib | 在图形坐标系 (`fig.add_axes`) 中手动绘制箭头和文字。 | 1. 完全自定义指南针的外观和位置。<br>2. 不依赖地图投影，位置固定于图形角落。<br>3. 视觉上更直观、突出。 | 1. 需要额外的绘图代码。<br>2. 指南针方向是固定的（上北），不随地图旋转而旋转。 | 需要传统、醒目的图形化指南针，且地图为正北朝上（未旋转）的投影。 | | **Basemap 自定义函数** | Basemap | 在地图数据坐标系中计算位置并绘制几何图形。 | 1. 指南针位置基于地图坐标（经纬度），可随意放置。<br>2. 函数化，易于复用。 | 1. Basemap 已弃用，未来兼容性无保障。<br>2. 代码相对冗长。<br>3. 在复杂投影下，指南针形状可能失真。 | 维护遗留的 Basemap 代码库，或需要在地图特定经纬度位置放置指南针。 | **选择建议与实践技巧：** 1. **库的选择**：对于新项目，**强烈推荐使用 `Cartopy`** 替代 `Basemap`，因为它有更好的性能、更现代的 API 和持续的维护 [ref_5]。 2. **指南针样式选择**： * 如果您的读者熟悉地图阅读，使用 **Cartopy 的经纬网格标签** 是更专业和简洁的选择。 * 如果图表用于大众传播或需要非常直观的方向指示，**自定义图形化指南针** 更合适。 3. **位置与投影**： * 使用 `Cartopy` 自定义图形时，通过 `fig.add_axes([left, bottom, width, height])` 指定指南针在**图形**中的相对位置（0到1之间），这使其位置不受地图投影和范围的影响。 * 在 `Basemap` 或希望指南针随地图移动的场景中，需要计算**地图数据坐标**（经纬度）来定位。 4. **动态旋转（高级）**：对于非正北朝上的地图（如某些方位投影），真正的指北针方向可能与图形上方不一致。实现动态旋转的指南针需要根据地图投影在给定点的北方向进行计算，这涉及更复杂的地理变换，通常需要从投影对象获取该点的北向矢量，然后相应地旋转绘制的箭头。在 `Cartopy` 中，可以通过 `crs`（坐标参考系统）的相关方法进行计算，但这超出了基础指南针添加的范围。 **总结**：在 Python 绘图中添加指南针，最优方案是使用 `Cartopy` 库。根据需求，既可利用其网格线标注功能进行隐式指示，也可结合 `matplotlib` 的绘图功能在图形固定位置绘制一个美观的静态指南针图形。应避免在新项目中使用已弃用的 `Basemap` [ref_5]。