首页 > 学院 > 开发设计 > 正文

【GMT】用gmtselect命令筛选数据

2019-11-08 00:57:30
字体:
来源:转载
供稿:网友

gmtselect

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

六个空间准则包括:

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

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

R=0/360/-70/70ps=res.psJ=M5iraw=raw2.txtsub=subset.txtgmt 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

结果:

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


发表评论 共有条评论
用户名: 密码:
验证码: 匿名发表