注意:本专栏的GMT脚本是在mac平台运行
(1)2D地形起伏图
#!/usr/bin/env bash
#Chistrong Wen
gmt begin hunza_dem png,pdf E600
gmt set FORMAT_GEO_MAP = ddd:mm:ssF
gmt set MAP_FRAME_TYPE = plain
gmt set MAP_FRAME_PEN = 1.2p,black
gmt set FONT_ANNOT_PRIMARY = 16p,Helvetica,black
gmt set FONT_ANNOT_SECONDARY = 16p,Times-Roman,black
gmt set MAP_TICK_PEN = 0.8p,black
gmt set MAP_TICK_LENGTH = 1.5p
gmt set LABEL_FONT = 16p
gmt set FONT_LABEL = 12p
gmt set LABEL_FONT_SIZE = 12p
gmt set ANNOT_FONT_SIZE_PRIMARY = 12p
gmt set MAP_FRAME_WIDTH = 1.5p
gmt set FRAM_WIDTH = 0.5c
gmt set MAP_ANNOT_OFFSET_PRIMARY = 0.25c
gmt set MAP_TICK_LENGTH_PRIMARY = 0.25c
gmt set FORMAT_ANNOT_PRIMARY = 2p,Times-Roman,black
# 绘制底图
# gmt grdview m1.grd -R74.52/74.68/36.25/36.4/1000/6000 -JZ5c -N-300+ggray -Qi -Ctest2.cpt -I+d -Ba -Bz -BwsenZ -p140/40 -Xc -Yc
# gmt basemap -p140/40/6000gmt basemap -R74.5/74.7/36.2/36.4 -Bf0.1a0.05 -JM16c -BWeSn -X22c -Yc
gmt grdimage m1.grd -Ctest2.cpt -I+d
gmt coast -A500 -W0.5p,black
gmt plot river.txt -W2p,blue
gmt plot road1.txt -W2p,yellow
gmt plot road2.txt -W2p,yellow
gmt plot 14_rgi60_SouthAsiaWest.shp -W2p,cyan -Gcyan
gmt plot path.txt -W2p,red# 沿测线提取地形高度
gmt basemap -R0/6000/2000/2800 -BeWsn -Byf200a200 -JX16c/3c -X0c -Y21c
gmt plot path_.txt -W2p,red -Ggray -L+y1000gmt end
需要注意的是,GMT目前只能绘制单一直线AB段的剖面。
如果要绘制沿曲线剖面(上图的红色曲线段),可以在matlab或者arcgis中提前得到path_.txt文件,再进行绘制
(2)3D地形起伏图
# Chistrong Wen
gmt begin dem_1 png,pdf E600
gmt set FORMAT_GEO_MAP = ddd:mm:ssF
gmt set MAP_FRAME_TYPE = plain
gmt set MAP_FRAME_PEN = 1.2p,black
gmt set FONT_ANNOT_PRIMARY = 16p,Helvetica,black
gmt set FONT_ANNOT_SECONDARY = 16p,Times-Roman,black
gmt set MAP_TICK_PEN = 0.8p,black
gmt set MAP_TICK_LENGTH = 1.5p
gmt set LABEL_FONT = 16p
gmt set FONT_LABEL = 12p
gmt set LABEL_FONT_SIZE = 12p
gmt set ANNOT_FONT_SIZE_PRIMARY = 12p
gmt set MAP_FRAME_WIDTH = 1.5p
gmt set FRAM_WIDTH = 0.5c
gmt set MAP_ANNOT_OFFSET_PRIMARY = 0.25c
gmt set MAP_TICK_LENGTH_PRIMARY = 0.25c
gmt set FORMAT_ANNOT_PRIMARY = 2p,Times-Roman,black
# 绘制底图
gmt grdview m1.grd -R74.52/74.68/36.25/36.4/1000/6000 -JZ5c -N-300+ggray -Qi -Ctest2.cpt -I+d -Ba -Bz -BwsenZ -p140/40 -Xc -Yc
gmt basemap -p140/40/6000
gmt end
绘制效果
代码中提到的数据见本专栏绑定的资源。
❤️欢迎点赞收藏❤️