地震随震中距与反方位角的分布图

示例贡献者

lxrat(作者)、陈箫翰(修订)


在接收函数、S波分裂等地震学研究中,经常需要绘制以台站为中心的地震分布图。下面的例子展示了如何绘制地震随震中距与反方位角的分布图:

#!/usr/bin/env bash
#
# 地震随震中距与反方位角的分布图
#

# 假定台站位于 105E,30N 处
stlo=105
stla=30
# 使用 -JE 投影(四个参数: 中心经度/中心纬度/最大震中距/图片宽度)
J=E$stlo/$stla/90/10c

gmt begin map png,pdf
gmt set FORMAT_GEO_MAP=+D
gmt coast -J$J -Rg -A10000 -Ggrey

# 绘制台站位置(三角形)
echo $stlo $stla | gmt plot -St0.4c -Gblack -Bya180
# 绘制地震位置(五角星)
gmt plot -Sa0.25c -Gred << EOF
145 30
160 40
60 40
80 -10
EOF

# 绘制等震中距线:30度(直径 60 度)、60 度(直径 120 度)
echo $stlo $stla 60d | gmt plot -SE- -W1p,red
echo $stlo $stla 120d | gmt plot -SE- -W1p,red
# 添加文字:30、60、90 度处
gmt text -D0c/0.3c << EOF
$stlo 0 30\232
$stlo -30 60\232
$stlo -59.9 90\232
EOF

gmt end show
../../_images/index-gmtplot-06.png