1000字范文,内容丰富有趣,学习的好帮手!
1000字范文 > GMT绘制子图 指北针 图例 比例尺

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

时间:2023-09-01 12:19:35

相关推荐

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-/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

效果图:

本内容不代表本网观点和政治立场,如有侵犯你的权益请联系我们处理。
网友评论
网友评论仅供其表达个人看法,并不表明网站立场。