利用GMT6.4及以上版本绘制DEM和地形起伏剖面图

首先，给大家介绍一个下载DEM的中科大镜像：Index of /gmtdata/ (ustc.edu.cn)，可以下载各种分辨率的DEM。本文是在手动下载好DEM的基础上进行的（主要是用GMT远程下载，总是不能下载或者没法将jp格式转为grd，哭了。。）

1、利用下载好的等高线数据，调整合适的CPT文件绘制DEM、测点A与B 和测线AB

gmt begin region png
set J=M10c
REM 以加利福尼亚地区为例
set R=235/244/33/43  
gmt set FONT_LABEL 20p,4,black
  gmt set FORMAT_GEO_MAP=ddd:mm:ssF
  gmt basemap -J%J% -R%R% -Baf
  gmt makecpt -Cglobe -T-10000/10000/200 -Z -D
  gmt grdimage earth_relief_30s.grd -Baf -BWSen -I+d
    REM # 选取测线AB
    echo 237 41 A > tmp1
    echo 241.5 34.2 B >> tmp1
    REM # 绘制测线AB
    gmt plot tmp1 -W2p,purple,-.-
    REM # 标注AB
    gmt text tmp1 -F+f15p -D0c/0.2c 
gmt end show

 2、绘制AB线所切得剖面的高程变化情况

gmt begin DEM_AB png
set J=M10cl
set R=235/244/33/43
gmt set FONT_LABEL 9p,4,black
    REM # 沿测线AB绘制地形高度
    gmt basemap -R0/7.65/-1000/2500 -JX12c/3c -Bya1000+l"Elevation (m)" -Bxaf1+l"Distance"+u"\260" -BWSrt -Y12c 
    REM # 标注AB位置
    echo 0 2800 A | gmt text -F+f10p+jBC -N
    echo 7.65 2800 B | gmt text -F+f10p+jBC -N
    REM # 沿测线提取地形高度
    gmt project -C237/41 -E241.5/34.2 -G0.1 | gmt grdtrack -Gearth_relief_04m.grd > tmp2
    REM # 将海平面以下填充为淡蓝色
    echo 0 0 > tmp3
    echo 7.65 0 >> tmp3
    gmt plot tmp3 -Wblack -Glightblue -L+y-1000
    REM # 将地形填充为灰色
    gmt plot tmp2 -i2,3 -Wblack -Ggray -L+y-1000 
gmt end show

参考信息：

1、GMT中文社区/中文手册：GMT 中文社区 (gmt-china.org)

2、社区绘图实例 -  绘制地形起伏剖面图 — GMT 中文手册 (gmt-china.org) 

致谢：特别感谢GMT中文社区志愿者的辛勤付出，你们辛苦了！有任何疑问，欢迎相互学习交流。