绘制三维图可以更好地表现物理场的空间分布和对比,也有很多的软件工具可以实现这一功能,本专栏主要强调使用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质量变化趋势信号
针对不同的坐标轴的展示,是由-BSEz123+b控制,在二维平面图中,通常仅使用-Bnews,而在三维图中,加上-BnewsZ,其中的1234代表Z的4条轴,不加上则表明不显示。而+b表明加上全部的轴。
❤️欢迎点赞收藏❤️