velo
- 官方文档:
- 简介:
在地图上绘制三种符号：矢量（向量）、十字线和旋转饼图
-S选项
- -S
-S选项决定要绘制哪一种符号，及其输入数据的格式。-Se[velscale/]confidence[+ffont] 绘制矢量，并可选择用椭圆表示不确定度。输入数据中每列的含义为：
起点的经度
起点的纬度
矢量的东向分量
矢量的北向分量
东向分量的不确定度（1-sigma）
北向分量的不确定度（1-sigma）
东向分量和北向分量的correlation
名称（可选）
说明：
velscale控制了矢量的缩放比例，其对应了速度值为1的箭头的长度，默认单位由参数 PROJ_LENGTH_UNIT 决定，也可以自己添加长度单位i|c|pconfidence设置了椭圆的2D置信区间的极限，比如0.95代表95%置信区间font是名称的字体属性。默认值[9p,Helvetica,black]，使用 +f0 则不显示名称。椭圆的填充色由
-G选项控制椭圆的轮廓以及矢量的画笔属性由
-W选项控
-Sr[velscale/]confidence[+ffont] 与
-Se类似，区别在于输入文件不同：起点的经度
起点的纬度
矢量的东向分量
矢量的北向分量
椭圆的半长轴
椭圆的半短轴
椭圆的长轴相对于水平方向逆时针旋转的角度
名称（可选）
-Sn[barscale] 绘制各向异性矢量，
barscale是各向异性矢量的缩放比例。输入数据中每列的含义为：经度
纬度
各向异性矢量的东向分量
各向异性矢量的北向分量
-Sw[wedgescale/]wedgemag 绘制旋转饼图。输入数据中每列的含义为：
经度
纬度
旋转饼图的旋转角度（弧度制）
旋转饼图的不确定度（弧度制）
旋转饼图由两个同心扇形组成。一个扇形是旋转饼图本身（下称饼图扇形），另一个扇形则表示旋转饼图的不确定度（下称不确定度扇形）。 两个扇形的半径由
wedgescale给出。两个扇形的角度分别为旋转饼图的旋转角度和不确定度乘以wedgemag。 饼图的填充色由-G选项控制， 不确定度的颜色由-E控制。 饼图扇形从正北开始，顺时针旋转。不确定度扇形从旋转饼图的结束开始，向顺时针和逆时针两个方向旋转，两个方向旋转的角度均为不确定度的 2 倍。-Sx[cross_scale] 绘制应变十字线。
cross_scale控制十字线的大小。输入文件中每列的含义为：经度
纬度
eps1：应变张量的扩张本征值，扩张为正
eps2： 应变张量的压缩本征值，扩张为正
eps2的方位角（北方向顺时针旋转的角度）
选项
- -Aparameters
修改矢量的属性。注意默认情况下，矢量头的大小不随着矢量线的长度变化而变化。 当矢量线过短时，会出现矢量箭头不绘制的情况。此时可以添加 +nnorm 属性让箭头大小随矢量成比例变化。 详细说明见 矢量/箭头 一节。
- -Ccpt
提供一个CPT文件，并让通常由
-G设置的符号颜色根据 magnitude 来确定。 具体参见-Z。- -Dscale
控制
-Se和-Sr中速度不确定度以及-Sw中旋转不确定度的缩放。该参数可以与<confidence>一起控制不确定度的绘制- -Efill
设置用于填充不确定性楔形图（
-Sw）或速度误差椭圆（-Se或-Sr）的颜色或阴影。如果未指定-E，则不确定性区域将会是透明的。 注意：使用-C和-Z+e将会根据-Z中选择的度量标准 [如 magnitude error ] 更新不确定性填充颜色。- -Gfill
设置符号的填充色 注意：使用
-C和-Z会覆盖本选项。- -L[pen[+c[f|l]]]
绘制线段。椭圆和断层面会根据
-W绘制出轮廓。 使用 +cl 则会使用-C的设置使用CPT的颜色。 使用 +cf 则只会将CPT设置应用在 error 误差填充上。 使用 +c 则会将CPT设置应用在线条和填充上。- -N
不跳过区域范围外的符号
- -W[pen][+c[f|l]]
设置速度箭头、椭圆轮廓以及断层边界的颜色 [默认: width = 0.25p, color = black, style = solid]。 使用 +cl 则会使用
-C的设置使用CPT的颜色。 使用 +cf 则只会将CPT设置应用在符号填充上。 使用 +c 则会将CPT设置应用在线条和填充上。- -Z[m|e|n|u][+e]
- 选择与
-C一起使用的量，以设置填充颜色。可选择： m：magnitude（矢量magnitude或旋转magnitude）
e：东西方向速度
n：南北方向速度
u：用户提供的数据列（位于必需列之后提供）
- 选择与
如果希望使用相应的误差估计值（即矢量或旋转的不确定性）来查找颜色并用于绘制误差椭圆或楔形图，请添加后缀 +e。
- -U[label][+c][+jjust][+odx/dy] (more …)
在图上绘制GMT时间戳logo
- -t[transp] (more …)
设置图层透明度（百分比）。取值范围为0（不透明）到100（全透明）
示例
#!/usr/bin/env bash
gmt begin velo_ex1
# 左边三个符号
gmt velo -JM15c -R241/242/34/35 -B0.5 -BWeSn+tvelo \
-Sw3c/1.e7 -W1p -G60 -E210 -D2 << EOF
# lon lat spin(rad/yr) spin_sigma (rad/yr)
241.2 34.2 5.65E-08 1.17E-08
241.2 34.5 -4.85E-08 1.85E-08
241.2 34.7 4.46E-09 3.07E-08
EOF
# 右边三个符号
gmt velo -Sw3c/1.e7 -W1p -D2 << EOF
# lon lat spin(rad/yr) spin_sigma (rad/yr)
241.7 34.2 5.65E-08 1.17E-08
241.7 34.5 -4.85E-08 1.85E-08
241.7 34.7 4.46E-09 3.07E-08
EOF
gmt coast -W1p -Di
gmt end show
velo示例图1
#!/usr/bin/env bash
gmt begin velo_ex2
gmt velo -R-10/10/-10/10 -Wthin,red -Se0.2c/0.39/12 -B1g1 \
-BWeSn -Jx0.2i -Ggreen -Eblue -L -N -A1c+p3p+e << EOF
# Long. Lat. Evel Nvel Esig Nsig CorEN SITE
# (deg) (deg) (mm/yr) (mm/yr)
-8. 0. 5.0 0.0 4.0 6.0 0.500 4x6
-5. 0. 0.0 5.0 4.0 6.0 0.500 4x6
-5. 5. -5.0 0.0 4.0 6.0 0.500 4x6
-8. 5. 0.0 -5.0 0.0 0.0 0.500 4x6
-1. 5. 3.0 3.0 1.0 1.0 0.100 3x3
EOF
gmt velo -Se0.2c/0.39/14 -A0.25c+p0.25p+e << EOF
# Long. Lat. Evel Nvel Esig Nsig CorEN SITE
# (deg) (deg) (mm/yr) (mm/yr)
0. -8. 0.0 0.0 4.0 6.0 0.100 4x6
-8. 5. 3.0 3.0 0.0 0.0 0.200 3x3
0. 0. 4.0 6.0 4.0 6.0 0.300
-5. -5. 6.0 4.0 6.0 4.0 0.400 6x4
5. 0. -6.0 4.0 6.0 4.0 -0.300 -6x4
0. -5. 6.0 -4.0 6.0 4.0 -0.500 6x-4
EOF
gmt end show
velo示例图2
下面的示例展示了如何绘制应变十字线:
gmt begin test pdf
gmt velo -JX10c/10c -R0/10/0/10 -Sx10 -W1p << EOF
5 5 0.5 0.3 45
EOF
gmt end show