GMT绘制GPS速度场

#!/bin/bash
gmt begin
gmt figure gps pdf,png A0.2c,E600
data=gps_data.txt
# 绘制底图
gmt set FORMAT_GEO_MAP=ddd:mm:ssF
gmt set FONT_ANNOT_PRIMARY 8p,4
gmt set MAP_FRAME_PEN 0.5p
# utm.grd的作用仅为指定绘图范围，且让边框使用平面坐标格式
# 经纬度数据
gmt basemap -Rutm.grd -Ju+48/1:8000000 -Blrtb  #--MAP_FRAME_PEN=1p,red
#gmt grdcut @earth_relief_30s_p -R98/111/31/42 -Glon-lat.grd
#gmt grdgradient lon-lat.grd -Ne0.7 -A50 -Glon-lat_i.grd
gmt grdimage lon-lat.grd -Ccolor.cpt -Ilon-lat_i.grd -Bag2f1 --MAP_GRID_PEN_PRIMARY=0.4p,gray,-
# 绘制GPS速度场
# -Se后面三参数表示：速度值为1的矢量的长度/置信度0.95/文本的大小
# -W控制矢量以及误差椭圆的轮廓的宽度，颜色，线型
# -G矢量填充色
# -A控制矢量的属性，0.15c是矢量头的大小，+e表示在矢量尾端绘制箭头，+p0.75p矢量线段部分的宽度
gmt velo $data -Se0.05c/0.95/0 -A0.15c+e+p0.75p -Gblue -W0.2p,blue

# 绘制插值得到的边界条件
gmt velo boundary_conditions.txt -Se0.05c/0.95/0 -A0.15c+e+p0.75p -Gred -W0.2p,red


# utm
gmt basemap -Jx1:8000000 -Rutm.grd -Brtlb #--MAP_FRAME_PEN=0.4p,blue
# 研究区范围
gmt plot -W0.5p << EOF
126200  3655000
126200  4435000
780300  4435000
780300  3655000
126200  3655000
EOF
gmt end