首先,给大家介绍一个下载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中文社区志愿者的辛勤付出,你们辛苦了!有任何疑问,欢迎相互学习交流。