绘制地形起伏剖面图

示例贡献者:

赵志远(作者)、陈箫翰(修订)


使用 GMT 绘制地形起伏剖面图的过程,可以简要概括为先利用 project 命令生成测线,然后使用 grdtrack 命令读取该测线上的地形数据起伏,最后使用 plot 命令绘制出地形剖面图。

下面的例子展示了如何基于以上方法绘制一条地形起伏剖面。例子中使用了精度为4弧分的地形数据,生成测线的数据点间隔为0.1度。读者应当根据自己的实际情况进行修改。

#!/usr/bin/env bash
#
# 绘制地形起伏剖面图

gmt begin profile

    # 绘制地形图
    gmt basemap -JM12c -R122/149/30/48 -Baf
    gmt grdimage @earth_relief_04m -Cglobe -I+d
    # 选取测线AB
    echo 126 42 A > tmp1
    echo 146 40 B >> tmp1
    # 绘制测线AB
    gmt plot tmp1 -W2p,red,-.-
    # 标注AB
    gmt text tmp1 -F+f15p -D0c/0.2c 
    
    # 沿测线AB绘制地形高度
    gmt basemap -R0/15/-8000/6000 -JX12c/3c -Bya4000+l"Elevation (m)" -Bxa2f1+l"Distance"+u"\260" -BWSrt -Y12c
    # 标注AB位置
    echo 0 7000 A | gmt text -F+f10p+jBC -N
    echo 15 7000 B | gmt text -F+f10p+jBC -N
    # 沿测线提取地形高度
    gmt project -C126/42 -E146/40 -G0.1 | gmt grdtrack -G@earth_relief_04m > tmp2
    # 将海平面以下填充为淡蓝色
    echo 0 0 > tmp3
    echo 15 0 >> tmp3
    gmt plot tmp3 -Wblack -Glightblue -L+y-8000
    # 将地形填充为灰色
    gmt plot tmp2 -i2,3 -Wblack -Ggray -L+y-8000 

    # 删除临时文件
    rm tmp*
gmt end show
../../_images/78a03ab725b96bca7227199575dffe9a.png