这里记录一下如何在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-/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-/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 > tempecho 75 10.2 >> tempecho 75 20 >> tempecho 65.2 20 >> tempgmt psxy temp -J%J% -R%R% -Gwhite -W0.1p,white -A -O -K>>%PS%echo 67 18 > temp2echo 70 18 >> temp2echo 70 19 >> temp2echo 67 19 >> temp2gmt 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_FORMATddd:mm:ssFgmtset MAP_FRAME_TYPE plaingmtset FRAM_WIDTH 0.04cgmtset LABEL_FONT 2gmtset LABEL_FONT_SIZE 8pgmtset ANNOT_FONT_SIZE_PRIMARY 8pgmtset MAP_FRAME_WIDTH 1pgmtset MAP_ANNOT_OFFSET_PRIMARY 0.1cset J=Cyl_stere/0/0/12cset R=65/110/10/35REM set J = N16cset PS="hydro.ps"rem %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%REM gmt grdcut ETOPO1_Bed_g_gdal.grd -R%R% -G1.grdREM gmt grdsample 1.grd -G1.grd -I0.01/0.01REM gmt grdgradient 1.grd -A270 -Ne1.2 -G1.intgmt makecpt -Cterra -T-8000/8000/100 -Z -D>1.cptgrdimage 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 > tempecho 75 10.2 >> tempecho 75 20 >> tempecho 65.2 20 >> tempgmt psxy temp -J%J% -R%R% -Gwhite -W0.1p,white -A -O -K>>%PS%echo 67 18 > temp2echo 70 18 >> temp2echo 70 19 >> temp2echo 67 19 >> temp2gmt 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/90set Jg=G90/21/3cgmt 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
效果图: