PyGMT结合Shapefile绘图时图例生成失败求助（附代码）

核心问题定位与修复步骤 #

1. 缺失tempfile模块导入 #

代码中使用了tempfile.NamedTemporaryFile但未导入该模块，会直接导致临时图例文件无法创建，图例加载失败。需在代码开头添加：

import tempfile

2. 图例条目格式修正 #

PyGMT图例规范严格遵循GMT语法，现有条目存在格式错误：

• 海沟符号条目：原样式参数将画笔属性（pen）嵌入style字段，不符合规范，需拆分到独立字段：

  "S 0.2i f0.5c/0.3c+r+t+o0.3c 0.25p,black 0.3i black 0.4i Trench"
  

  格式说明：S <符号尺寸> <符号样式> <画笔参数> <符号与文本间距> <填充色> <文本缩进> <标签>
• 颜色条条目：G条目需补充颜色条标注格式，确保与GMT语法匹配：

  f"G 0.2i - 1i ./combined.cpt 0.4i -B+l\"Elevation (m)\" -Baf+u\"m\" {min_elevation}/{max_elevation}"
  

3. 临时文件路径与权限优化 #

避免系统临时目录的权限问题，显式指定临时文件存放在当前工作目录：

with tempfile.NamedTemporaryFile(delete=False, mode="w", suffix=".txt", dir="./") as legend_file:

4. 图例位置与图层顺序检查 #

• 调整图例偏移量，避免超出绘图区域：将o0.5i改为o0.2i测试
• 确保fig1.legend()调用在所有绘图命令之后，防止被其他图层覆盖

修正后的图例相关代码片段 #

import tempfile  # 新增导入

# ... 原有代码省略 ...

min_elevation = -5500
max_elevation = 4000

# 修正后的图例条目
legend_entries = [
    "S 0.2i f0.5c/0.3c+r+t+o0.3c 0.25p,black 0.3i black 0.4i Trench",
    "S 0.2i - 0.25p,black 0.3i - 0.4i Contours",
    f"G 0.2i - 1i ./combined.cpt 0.4i -B+l\"Elevation (m)\" -Baf+u\"m\" {min_elevation}/{max_elevation}"
]

# 显式指定临时文件存放目录
with tempfile.NamedTemporaryFile(delete=False, mode="w", suffix=".txt", dir="./") as legend_file:
    for entry in legend_entries:
        legend_file.write(entry + "\n")
    legend_file_name = legend_file.name

# 调整图例位置偏移
fig1.legend(position="JTR+jTR+o0.2i", box="+gwhite+p2p", spec=legend_file_name)

fig1.show()

额外排查建议 #

• 验证combined.cpt文件完整性：打开文件确认包含海洋、陆地的颜色映射数据，无格式错误
• 简化图例测试：先保留单个条目（如Contours），验证是否显示，逐步定位问题条目
• 启用调试日志：查看图例加载的错误信息

  import pygmt
  pygmt.set_log_level("debug")
  

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