GMT6.4绘图进阶：给你的地形剖面图加上高程填充与海平面标识
GMT6.4地形剖面图高阶美化：从基础线条到专业出版级图表
当一张地形剖面图从单调的黑白线条演变为带有高程填充色、清晰海平面标识的彩色图表时，其传达信息的效率会呈指数级提升。作为地球科学领域的标准工具，GMT6.4提供了强大的图形定制能力，能让你的剖面图在学术报告或期刊出版物中脱颖而出。
1. 数据准备与环境配置
在开始绘制前，确保已获取适当分辨率的DEM数据。对于区域尺度研究，30弧秒（约1公里）分辨率的earth_relief_30s.grd通常足够；而局部精细分析可能需要15弧秒或更高分辨率数据。
# 检查GMT版本及数据可用性 gmt --version gmt which @earth_relief_30s --show建议：将常用参数设为变量便于后续调用：
REGION="235/244/33/43" # 研究区域经纬度范围 PROFILE_START="237/41" # 剖面起点A坐标 PROFILE_END="241.5/34.2" # 剖面终点B坐标注意：使用远程DEM数据时可能遇到下载中断问题，可预先通过
gmt grdconvert将数据转换为本地GRD格式保存。
2. 基础剖面绘制与高程提取
剖面线生成是后续美化的基础。通过gmt grdtrack提取的高程数据需要合理处理海平面参考：
# 生成AB剖面线并提取高程 gmt project -C${PROFILE_START} -E${PROFILE_END} -G0.1 \ | gmt grdtrack -G@earth_relief_30s > profile_raw.txt原始数据格式示例：
| 经度 | 纬度 | 距离(km) | 高程(m) |
|---|---|---|---|
| 237.0 | 41.0 | 0.0 | 1280 |
| 237.1 | 40.9 | 12.3 | 1150 |
关键参数说明：
-G0.1：设置采样间隔为0.1公里-i2,3：指定使用第3、4列数据（距离和高程）
3. 智能填充技术实现
3.1 海平面上下分色填充
利用-G和-L参数组合实现分色填充：
# 创建海平面参考线 echo "0 0" > sea_level.txt echo "$(tail -1 profile_raw.txt | awk '{print $3}') 0" >> sea_level.txt # 绘制填充效果 gmt plot sea_level.txt -W0.5p,blue -Glightblue -L+y-10000 gmt plot profile_raw.txt -i2,3 -W1p,darkgray -Ggray -L+y-10000颜色搭配建议：
- 海平面以下：浅蓝色系（
lightblue/azure） - 海平面以上：土色系（
gray/tan/bisque） - 特殊地貌：可通过
-G+z实现高程渐变
3.2 多图层叠加技巧
通过控制-X/-Y参数实现多图组合：
gmt begin profile_final png # 主剖面图 gmt plot [...] -JX15c/5c -R0/800/-2000/3000 # 插入小地图定位 gmt inset begin -DjTR+w4c -F+gwhite gmt coast -R${REGION} -Glightbrown -Slightblue -Baf gmt plot -W1p,red << EOF ${PROFILE_START} ${PROFILE_END} EOF gmt inset end gmt end show4. 专业标注与出版级优化
4.1 关键点标注增强
# 标注剖面端点 echo "0 3200 A" | gmt text -F+f12p,Helvetica-Bold,red+jBC -N echo "800 3200 B" | gmt text -F+f12p,Helvetica-Bold,red+jBC -N # 添加比例尺 gmt plot -Sv0.5c+e -W1p << EOF 400 2500 0 500 EOF echo "400 2300 500 km" | gmt text -F+f10p,black+jTC4.2 字体与样式规范
学术出版推荐设置：
gmt set FONT_ANNOT_PRIMARY 12p,Helvetica,black gmt set FONT_LABEL 14p,Helvetica-Bold,black gmt set MAP_FRAME_TYPE plain gmt set MAP_TICK_LENGTH_PRIMARY 5p4.3 输出格式优化
不同用途的输出建议：
- 期刊投稿：EPS或PDF矢量格式
- 网页展示：PNG（600dpi）
- 海报打印：TIFF（1200dpi）
# 高质量PNG输出示例 gmt begin my_profile png gmt set PS_MEDIA A4 gmt set FORMAT_GEO_MAP ddd.xx [...绘图命令...] gmt end show在实际项目制图中，我发现最影响最终效果的是颜色对比度和标注清晰度。建议在正式输出前，先用gmt preview快速检查各元素的可读性，特别是当图表需要黑白印刷时，要确保灰度对比足够明显。