0%

grdview

发表于 更新于 阅读次数： Valine：

今天分享的GMT绘图比较常见：地形图。以印度洋的洋中脊的天成热液口为例子，简单介绍怎么绘制一副高质量的论文插图。

数据

绘图数据为多波束调查数据，来源于中国大洋科学考察航次。

代码

这个绘图比较简单，就不多写了。如果有问题，很可能是GMT版本造成，一般微调一下即可达到效果。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
rem 
gmt gmtset FORMAT_GEO_MAP = ddd:mmF
gmt gmtset MAP_FRAME_WIDTH=2p
gmt gmtset FONT_ANNOT_PRIMARY 7p,Helvetica,black FONT_LABEL 8p,Helvetica,black 

set ps=3d3.ps
set R=-R63.6/64.1/-28.1/-27.70
set R3=%R%/-5500/-2000
gmt grd2xyz 26VI_DBS02.grd %R% | gawk "!/NaN/ {print $1,$2,$3}" | gmt makecpt -Cseafloor -E200 -i2> t.cpt
gmt grdgradient 26VI_DBS02.grd -A0 -Gt_intens.nc -Nt0.75 %R%

gmt grdview 26VI_DBS02.grd -It_intens.nc %R3% -JM5i -p40/40 -JZ3i -Ct.cpt  -K -Y0.5i -Qs -B5m -N-5500+glightgray -BNESW > %ps%
echo 63.9233156 -27.8508410 JL-87| gmt grdtrack -G26VI_DBS02.grd  | gmt psxyz -R -J -JZ -Sc0.1c -Gblack -O -K -p>> %ps%
echo 63.9227	-27.8508 TianTang| gmt grdtrack -G26VI_DBS02.grd  | gmt psxyz -R -J -JZ -Sc0.1c -Gblack -O -K -p>> %ps%
echo 63.9218	-27.8504 TC-2| gmt grdtrack -G26VI_DBS02.grd  | gmt psxyz -R -J -JZ -Sc0.1c -Gblack -O -K -p>> %ps%
echo 4.25 5.5 T@#ian @# C@#heng@# | gmt pstext -R0/10/0/10 -Jx1i -F+f20p,4+jCB -O -K >> %ps%
gmt psscale %R% -JM5i -p -DjTC+o0/-0.7i+w3i/0.15i+h+e -Ct.cpt -I -O -Bx500 -By+lm  --FONT_LABEL=10p,Helvetica,black >> %ps%
gmt psconvert %ps% -A -P -Tg

结果

发送0301（微信公众号：海洋遥感数据共享），获得多波束网格数据。（数据分辨率较高，仅供学习交流，勿作他用。）

参考文献：[1]J. Sun et al., “Nearest vent, dearest friend: biodiversity of Tiancheng vent field reveals cross-ridge similarities in the Indian Ocean,” Royal Society Open Science, vol. 7, no. 3, p. 200110, 2020, doi: 10.1098/rsos.200110.