GMT绘制子图、指北针、图例、比例尺

这里记录一下如何在GMT中添加子图、指北针、比例尺和图例,以防每次全网搜索找不齐全。下面是具体的内容:

添加子图代码

gmt psbasemap -R%Rg% -J%Jg% -D%R% -Bxa5 -Bya5 -Bg30 -BWeSn -Fd+p2p,red+c0.1 -X0c -Y0c -O>>%PS%


2.添加指北针:

https://docs.gmt-china.org/5.4/module/psbasemap/
gmt psbasemap -R%R% -J%J% -Bxa5 -Bya5 -BWeSn -X0c -Y0c -Tdg105/30+w0.5i+jCM -O -K>>%PS%

此处我使用的是GMT5,对于GMT6,有另外一种方法。参考:https://docs.gmt-china.org/latest/examples/ex007/#gmtplot-7adb3518e066c53838e13ab7c9582743
3.添加比例尺

pscoast -R%R% -J%J% -Dh -N1 -A1000  -I1 -Wthinnest,black -Lg70/12+c12+w500k+f+u -X0c -Y0c -O -K >>%ps%


4.添加图例

echo 65.2 10.2 > temp
echo 75 10.2 >> temp
echo 75 20 >> temp
echo 65.2 20 >> temp
gmt psxy temp -J%J% -R%R% -Gwhite -W0.1p,white -A -O -K>>%PS%

echo 67 18 > temp2
echo 70 18 >> temp2
echo 70 19 >> temp2
echo 67 19 >> temp2
gmt psxy temp2 -J%J% -R%R% -Gdarkyellow -W0.2p,darkyellow  -A -O -K>>%PS%
echo 72 18.4 Sediment | gmt pstext -J%J% -R%R% -F+f6p,4,black -O -K >>%PS%
echo 70.5 17 Chistrong_UCAS | gmt pstext -J%J% -R%R% -F+f6p,4,black -O -K >>%PS%


完整代码:

gmtset PLOT_DEGREE_FORMAT           ddd:mm:ssF
gmtset MAP_FRAME_TYPE plain
gmtset FRAM_WIDTH 0.04c
gmtset LABEL_FONT 2
gmtset LABEL_FONT_SIZE 8p
gmtset ANNOT_FONT_SIZE_PRIMARY 8p
gmtset MAP_FRAME_WIDTH 1p
gmtset MAP_ANNOT_OFFSET_PRIMARY 0.1c
set J=Cyl_stere/0/0/12c
set R=65/110/10/35
REM set J = N16c
set PS="hydro.ps"
rem %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
REM gmt grdcut ETOPO1_Bed_g_gdal.grd -R%R% -G1.grd
REM gmt grdsample 1.grd -G1.grd -I0.01/0.01
REM gmt grdgradient 1.grd -A270 -Ne1.2 -G1.int
gmt makecpt -Cterra -T-8000/8000/100 -Z -D>1.cpt
grdimage 1.grd  -R%R% -J%J% -B10f10/10f10wens  -C1.cpt -I1.int -Xc -Yc -K>%ps%
gmt psbasemap -R%R% -J%J% -Bxa5 -Bya5 -BWeSn  -X0c -Y0c -Tdg105/30+w0.5i+jCM -O -K>>%PS%
pscoast -R%R% -J%J% -Dh -N1 -A1000  -I1 -Wthinnest,black -Lg70/12+c12+w500k+f+u -X0c -Y0c -O -K >>%ps%
echo 65.2 10.2 > temp
echo 75 10.2 >> temp
echo 75 20 >> temp
echo 65.2 20 >> temp
gmt psxy temp -J%J% -R%R% -Gwhite -W0.1p,white -A -O -K>>%PS%
echo 67 18 > temp2
echo 70 18 >> temp2
echo 70 19 >> temp2
echo 67 19 >> temp2
gmt psxy temp2 -J%J% -R%R% -Gdarkyellow -W0.2p,darkyellow  -A -O -K>>%PS%
echo 72 18.4 Sediment | gmt pstext -J%J% -R%R% -F+f6p,4,black -O -K >>%PS%
echo 70.5 17 Chistrong_UCAS | gmt pstext -J%J% -R%R% -F+f6p,4,black -O -K >>%PS%
echo 105 28 N |gmt pstext -J%J% -R%R% -F+f6p,4,black -Gwhite -O -K >>%PS%
gmt psxy ganges.txt -J%J% -R%R% -Gdarkyellow -W0.1p,white -A -O -K>>%PS%
pscoast -R%R% -J%J% -Df -Wthinnest,black -Lg70/15+c40+w500k+f+u -X0c -Y0c -O -K >>%ps%
gmt psscale -Dx0.5c/-1.2c+jTC+w1.8c/0.15c+h+e --FONT_ANNOT_PRIMARY=8p,Helvetica,black -X0.8c -Y2c -O -K -C1.cpt -Ba5000 -By+lm>> %ps%

set Rg=-180/180/-90/90
set Jg=G90/21/3c
gmt pscoast -R%Rg% -J%Jg% -Glightyellow -Slightblue -N1 -A1000 -Wthinnest -O -K -X-1c -Y3c --MAP_FRAME_TYPE=plain>> %PS%
gmt psbasemap -R%Rg% -J%Jg% -D%R% -Bxa5 -Bya5 -Bg30 -BWeSn -Fd+p2p,red+c0.1 -X0c -Y0c -O>>%PS%
psconvert %PS% -A  -Tg  -E512


效果图:

  • 4
    点赞
  • 18
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

我是水怪的哥

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值