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

效果图：