【GMT5】用gmtselect命令筛选数据

gmtselect

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

六个空间准则包括：

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

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

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

结果：

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