GMT绘图笔记：绘制堆叠三维图

news2024/12/29 11:13:20

绘制三维图可以更好地表现物理场的空间分布和对比，也有很多的软件工具可以实现这一功能，本专栏主要强调使用GMT绘制。语法也非常简单，实例代码如下：

gmt begin tibet_moho_3D png,pdf E600
# 绘制底图
gmt set FORMAT_GEO_MAP = ddd:mm:ssF
gmt set MAP_FRAME_TYPE = plain
gmt set MAP_FRAME_PEN = 1p,black
gmt set FONT_ANNOT_PRIMARY = 14p,Helvetica,black
gmt set FONT_ANNOT_SECONDARY = 14p,Times-Roman,black
gmt set MAP_TICK_PEN = 1p,black
gmt set MAP_TICK_LENGTH = 1.5p
gmt set LABEL_FONT = 9p
gmt set FONT_LABEL = 9p
gmt set LABEL_FONT_SIZE = 9p
gmt set ANNOT_FONT_SIZE_PRIMARY = 9p
gmt set MAP_FRAME_WIDTH = 0.5p 
gmt set FRAM_WIDTH = 0.5c
gmt set MAP_ANNOT_OFFSET_PRIMARY = 0.2c
gmt set MAP_TICK_LENGTH_PRIMARY = 0.2c
gmt set FORMAT_ANNOT_PRIMARY = 2p,Times-Roman,black
gmt xyz2grd depthtomoho.txt -Rg -I1.1 -Gy1.grd
gmt grdsample y1.grd -Gb11.grd -I0.5
gmt xyz2grd csr.txt -Rg -I1.1 -Gm1.grd
gmt grdsample m1.grd -Ga11.grd -I0.5
gmt grdview b11.grd -JX8c/6c -R70/112/20/45/-70000/75000 -JZ5c -Qi -Cjet -I+d -Bxa10 -Bya10 -Bzaf+lm -p130/30 -Xc -Yc 
gmt grdview earth_relief_01m.grd -R70/112/20/45/-70000/75000 -JZ5c -Qi -Ctest.cpt -I+d -BSEz123+b -p130/30  -X0c -Y0c 
gmt grdview a11.grd -JX8c/6c -R70/112/20/45/-70000/75000 -JZ5c -Qi -Cjet -I+d -Bxa10 -Bya10 -Bzaf+lm -p130/30 -X0c -Y0c 
gmt end

结果图：这一图堆叠了（从底到顶）莫霍面深度、地形起伏和GRACE质量变化趋势信号