【GMT】用gmtselect命令筛选数据

gmtselect

该命令会从输入文件中读取前两列作为经度纬度,并判断经纬度是否符合一定的空间准则,以筛选出符合条件的记录。输入文件中仅前两列会被使用。

六个空间准则包括:

  1. 在矩形区域内(-R和-J)
  2. 与ptfile中的每个点的距离在dist之内
  3. 与linefile中的每条线的距离在distkm 之内
  4. 在polygonfile中的多边形内
  5. 在某个地理区域内(需要海岸线数据)
  6. z值在某个范围内

参考:《GMT5中文手册-seisman》

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
R=0/360/-70/70
ps=res.ps
J=M5i
raw=raw2.txt
sub=subset.txt
gmt gmtselect $raw -R$R -Dh -Nk/s >$sub
# ps文件头
gmt psxy -J$J -R$R -T -K > $ps
# 海岸线
gmt pscoast -R$R -J$J -Dl -Cblue -Baf -W1/0.25p,white -K -Slightblue -O >> $ps
#筛选前的数据,黄色
gmt psxy $raw -R$R -J$J -Sc0.03c -Gyellow -K -O >> $ps
#筛选后的数据,绿色
gmt psxy $sub -R$R -J$J -Sc0.03c -Ggreen -K -O >> $ps
# ps文件尾
gmt psxy -R$R -J$J -T -O >> $ps
# 保存为图片
gmt psconvert $ps -Tg -P -A

结果:

黄色的是筛选前的点,绿色是筛选后的点,绿色的点把黄色的覆盖了一部分。