GMT绘制GPS箭头代码
具体代码如下：其中changjaing.txt为流域边界数据，为经度、纬度两列数据，negative.txt为GPS站点数据，包含经度、纬度、东西向分量（为0）,南北向分量
REM 2021.3.29/ Z.Q.WEN ECJTU
set R=70/140/15/55
REM set J=j0.1i
set J = M8c
gmt set FORMAT_GEO_MAP=ddd:mm:ssF
gmt set MAP_FRAME_TYPE=plain
gmt set MAP_TICK_LENGTH=10p
gmt set FONT_LABEL 7p,Times-Roman
gmt set FONT_TITLE 7p,Times-Roman
gmt set MAP_TITLE_OFFSET=0.4c
gmt set PS_MEDIA=a3
gmt set FONT_ANNOT_PRIMARY 8p,Times-Roman
gmtset MAP_FRAME_TYPE plain
set PS="indus_ganges_.ps"
REM gmt grdcut ETOPO1_Bed_g_gdal.grd -R%R% -GcutTopo.grd
REM gmt grdgradient cutTopo.grd -Ne0.7 -A50 -GcutTopo_i.grd
REM gmt makecpt -Cglobe -T-8000/5000 -Di >t.cpt
gmt makecpt -Cgray -T1000/20000/100 -I >t.cpt
gmt grdimage cutTopo.grd -IcutTopo_i.grd -R%R% -J%J% -Xc -Yc -Ct.cpt -K >%PS%
gmt pscoast -R%R% -J%J% -Slightblue -A100 -W0.1p -O -K >>%PS%
gmt xyz2grd speed.txt -R%R% -I1 -Gtmp.grd
gmt grdsample tmp.grd -Gtmp.grd -I0.01
gmt makecpt -Crainbow -T-106/7.7/1 -D >tem.cpt
REM gmt psclip load.txt -J%J% -R%R% -Am -O -K >>%PS%
REM gmt grdimage tmp.grd -R%R% -J%J% -Ctem.cpt -B5f1/5f1WenS -t80 -O -K>>%PS%
gmt psxy gps.txt -J%J% -R%R% -Sc0.01c -Gblue -O -K>> %ps%
REM gmt psclip -C -O -K >>%PS%
gmt psxy changjiang.txt -R -J -O -K -W2p,gray44 >> %ps%
echo 108 30 1|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
gmt psxy hexi.txt -R -J -O -K -W2p,gray44>> %ps%
echo 100 40 2|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
gmt psxy huanghe.txt -R -J -O -K -W2p,gray44>> %ps%
echo 109 37 3|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
gmt psxy tp.txt -R -J -O -K -W2p,gray44>> %ps%
echo 86 33 4|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
gmt psxy neimenggu.txt -R -J -O -K -W2p,gray44>> %ps%
echo 113 43 5|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
gmt psxy southwest.txt -R -J -O -K -W2p,gray44 >> %ps%
echo 95 30 6|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
gmt psxy southeast.txt -R -J -O -K -W2p,gray44 >> %ps%
echo 119 27 7|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
gmt psxy tarim.txt -R -J -O -K -W2p,gray44 >> %ps%
echo 85 40 8|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
gmt psxy northeast.txt -R -J -O -K -W2p,gray44>> %ps%
echo 125 47 9|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
gmt psxy zj.txt -R -J -O -K -W2p,gray44>> %ps%
echo 110 24 10|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
gmt psxy zhungaer.txt -R -J -O -K -W2p,gray44 >> %ps%
echo 88 45 11|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman>> %ps%
gmt psxy haihe.txt -R -J -O -K-W2p,gray44>> %ps%
echo 116 38 12|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
gmt psxy huaihe1.txt -R -J -O -K -W2p,gray44>> %ps%
echo 118 34 13|gmt pstext -R%R% -J%J% -O -K -F+f10p,Times-Roman >> %ps%
REM gawk "{print $1,$2,$3}" kk.txt |
echo 70 15 > temp
echo 80 15 >> temp
echo 80 25 >> temp
echo 70 25 >> temp
gmt psxy temp -J -R -Gwhite -W0.1p,white -A -O -K >>%PS%
gmt psvelo positive.txt -R%R% -J%J% -Se0.2c/0.95/0 -A0.4c+e+p0.1p -Gred -W0.2p,red -O -K -N >>%PS%
gmt psvelo negative.txt -R%R% -J%J% -Se0.2c/0.95/0 -A0.4c+e+p0.1p -Gblue -W0.2p,blue -O -K -N >>%PS%
gmt psvelo nc.txt -R%R% -J%J% -Se0.05c/0.95/0 -A0.4c+e+p0.1p -Ggreen -W0.2p,green -O -K -N >>%PS%
echo 75 23 4 \261 1 mm/yr | gmt pstext -J -R -F+f7p,4,red -O -K >>%PS%
echo 75 18 -4 \261 1 mm/yr | gmt pstext -J -R -F+f7p,4,blue -O -K >>%PS%
echo 75 16 -16 \261 1 mm/yr | gmt pstext -J -R -F+f7p,4,green -O -K >>%PS%
REM echo 73 19 0 5 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0c+e+p0.1p -W0.1p,Gray45 -O -K >>%PS%
echo 73 19 0 4 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0.2c+e+p0.1p -Gred -W0.2p,red -O -K >>%PS%
REM echo 75 22 0 -5 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0c+e+p0.1p -W0.1p,Gray45 -O -K >>%PS%
echo 75 22 0 -4 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0.2c+e+p0.1p -Gblue -W0.2p,blue -O -K >>%PS%
REM echo 77 22 0 -5 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0c+e+p0.1p -W0.1p,Gray45 -O -K >>%PS%
echo 77 22 0 -4 0 0 0 | gmt psvelo -J -R -Se0.2c/0.95/0 -A0.2c+e+p0.1p -Ggreen -W0.2p,blue -O -K >>%PS%
gmt psbasemap -R%R% -J%J% -Bxa10 -Bya5 -BWesN -Xc -Yc -O >>%PS%
REM gmt psscale -Dx0.5c/-01.2c+jTC+w6c/0.4c+v+e --FONT_ANNOT_PRIMARY=14p,Helvetica,black -X1c -Y12c -O -Ctem.cpt -Ba25 -By+lmm/a>> %ps%
REM 绘制图中图
REM set Rg=-180/180/-90/90
REM set Jg=G78/32/5c
REM gmt grdimage ETOPO1_Bed_g_gdal.grd -R%Rg% -J%Jg% -Ct.cpt -X11.5c -Y-5.5c -O -K >>%PS%
REM gmt pscoast -R%Rg% -J%Jg% -N1 -A1000 -B0 -Wthinnest -O -K -X0c -Y0c --MAP_FRAME_TYPE=plain>> %PS%
REM gmt psbasemap -R%Rg% -J%Jg% -D%R% -Bf30a30 -Bg60 -Fd+p2p,red+c0.1 -O>>%PS%
gmt psconvert %PS% -A -Tg -P
