13.55 psvelo

官方文档:psvelo
简介:在地图上绘制三种符号：矢量（向量）、十字线和旋转饼图

13.55.1 -S选项

-S 选项决定要绘制哪一种符号，及其输入数据的格式。

-Se<velscale>/<confidence>/<fontsize>

绘制矢量，并可选择用椭圆表示不确定度。输入数据中每列的含义为：

  1. 起点的经度
  2. 起点的纬度
  3. 矢量的东向分量
  4. 矢量的北向分量
  5. 东向分量的不确定度（1-sigma）
  6. 北向分量的不确定度（1-sigma）
  7. 东向分量和北向分量的correlation
  8. 名称（可选）

说明：

  1. <velscale> 控制了矢量的缩放比例，其对应了速度值为1的箭头的长度，默认单位由参数 PROJ_LENGTH_UNIT 决定，也可以自己添加长度单位 i|c|p
  2. <confidence> 设置了椭圆的2D置信区间的极限，比如0.95代表95%置信区间
  3. <fontsize> 是名称的文本字号
  4. 椭圆的填充色由 -G 选项控制
  5. 椭圆的轮廓以及矢量的画笔属性由 -W 选项控
-Sr<velscale>/<confidence>/<fontsize>

与 -Se 类似，区别在于输入文件不同：

  1. 起点的经度
  2. 起点的纬度
  3. 矢量的东向分量
  4. 矢量的北向分量
  5. 椭圆的半长轴
  6. 椭圆的半短轴
  7. 椭圆的长轴相对于水平方向逆时针旋转的角度
  8. 名称（可选）
-Sn<barscale>

绘制各向异性矢量， <barscale> 是各向异性矢量的缩放比例。输入数据中每列的含义为：

  1. 经度
  2. 纬度
  3. 各向异性矢量的东向分量
  4. 各向异性矢量的北向分量
-Sw<wedge_scale>/<wedge_mag>

绘制旋转饼图。输入数据中每列的含义为：

  1. 经度
  2. 纬度
  3. 旋转饼图的旋转角度（弧度制）
  4. 旋转饼图的不确定度（弧度制）

旋转饼图由两个同心扇形组成。一个扇形是旋转饼图本身（下称饼图扇形），另一个扇形则表示旋转饼图的不确定度（下称不确定度扇形）。两个扇形的半径由 <wedge_scale> 给出。两个扇形的角度分别为旋转饼图的旋转角度和不确定度乘以 <wedge_mag> 。饼图的填充色由 -G 选项控制， 不确定度的颜色由 -E 控制。饼图扇形从正北开始，顺时针旋转。不确定度扇形从旋转饼图的结束开始，向顺时针和逆时针两个方向旋转，两个方向旋转的角度均为不确定度的 2 倍。

-Sx<cross_scale>

绘制应变十字线。 <cross_scale> 控制十字线的大小。输入文件中每列的含义为：

  1. 经度
  2. 纬度
  3. eps1：应变张量的扩张本征值，扩张为正
  4. eps2： 应变张量的压缩本征值，扩张为正
  5. eps2的方位角（北方向顺时针旋转的角度）

13.55.2 选项

-A<params>
修改矢量的属性，见 绘制矢量/箭头 一节
-D<sigma_scale>
控制 -Se 和 -Sr 中速度不确定度以及 -Sw 中旋转不确定度的缩放。该参数可以与 <confidence> 一起控制不确定度的绘制
-F<fill>
控制边框和标注的颜色
-E<fill>
控制 -Sw 中饼图不确定度的填充
-G<fill>
控制颜色？
-L
绘制线段。使用该选择，则椭圆和断层面会根据 -W 选项绘制轮廓
-N
不跳过区域范围外的符号
-W<pen>
设置速度箭头、椭圆轮廓以及断层边界的颜色

13.55.3 示例

#!/bin/bash
ps=psvelo_ex1.ps

# 左边三个符号
gmt psvelo << EOF -Xc -JM15c -R241/242/34/35 -B0.5 -BWeSn+tpsvelo \
    -Sw3c/1.e7 -W1p -G60 -E210 -D2 -P -K > $ps
# 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 psvelo -J -R -Sw3c/1.e7 -W1p -D2 -O -K << EOF >> $ps
# 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 pscoast -O -R -J -W1p -Di >> $ps
rm gmt.*
../../_images/psvelo_ex1.png

psvelo示例图1

#!/bin/bash
ps=psvelo_ex2.ps

gmt psvelo -R-10/10/-10/10 -Wthin,red \
	-Se0.2c/0.39/12 -B1g1 -BWeSn -Jx0.2i -Ggreen -Eblue -L -N \
	-A1c+p3p+e -K << EOF > $ps
# 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 psvelo -Se0.2c/0.39/14 -R -J -A0.25c+p0.25p+e -O << EOF >> $ps
# 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
rm gmt.*
../../_images/psvelo_ex2.png

psvelo示例图2

下面的示例展示了如何绘制应变十字线:

gmt psvelo -JX10c/10c -R0/10/0/10 -Sx10 -W1p > test.ps << EOF
5 5 0.5 0.3 45
EOF