plot画符号

plot画符号

要绘制符号,必须要有符号的经纬度坐标。

以下面坐标点为例 :

经度 纬度
104.39   29.90
 13.14   52.50
 19.99  -34.52
-77.15   38.89 
-52.47  -31.62
150.36  -33.79
 76.22  -69.22

下面的例子中,我们首先使用 UNIX 下的 cat 命令将数据写入到文件 stations.txt 中,然后使用 coast 模块绘制了一张底图,并使用 plot 模块绘制输入文件 stations.txt 中的数据。

下面的几行代码的作用是将两个 EOF 中间的三行数据保存到文件 stations.txt 中:

cat > stations.txt << EOF
104.39   29.90
 13.14   52.50
 19.99  -34.52
-77.15   38.89 
-52.47  -31.62
150.36  -33.79
 76.22  -69.22
EOF

下面的脚本中,采用 plot 模块绘制 stations.txt 给出的坐标点,并给 plot 添加了 -St0.5c 选项,绘制三角形符号,外接圆直径为0.5c,G` 设置符号的填充色为red

gmt begin figure1 png,pdf
    gmt coast -JH15c -Rg -Baf -W0.1p -Cblue -Ggreen -Slightblue
    gmt plot -St0.5c stations.txt -Gred
gmt end show

可以采用plot继续绘制坐标点:

cat > epicenter.txt << EOF
130.72 32.78
EOF

gmt begin figure1 png,pdf
    gmt coast -JH15c -Rg -Baf -W0.1p -Cblue -Ggreen -Slightblue
    gmt plot -St0.5c stations.txt -Gred
    gmt plot -Sa0.5c epicenter.txt -Gyellow -W1p,black
gmt end show

完整的绘图脚本如下:

#!/usr/bin/env bash
# GMT modern mode bash template
# Date:    2019-09-10T00:44:39
# User:    seisman
# Purpose: Purpose of this script
export GMT_SESSION_NAME=$$  # Set a unique session name

cat > stations.txt << EOF
104.39   29.90
 13.14   52.50
 19.99  -34.52
-77.15   38.89 
-52.47  -31.62
150.36  -33.79
 76.22  -69.22
EOF

cat > epicenter.txt << EOF
130.72 32.78
EOF

gmt begin figure1 png,pdf
    gmt coast -JH15c -Rg -Baf -W0.1p -Cblue -Ggreen -Slightblue
    gmt plot -St0.5c stations.txt -Gred
    gmt plot -Sa0.5c epicenter.txt -Gyellow -W1p,black
gmt end show

  • -S 选项相对复杂,与不同的选项连用,或者后面接不同的参数,所需要的输入数据的格式也不同。不管是什么符号,至少都需要给定符号的位置,即X和Y是必须的:
X   Y