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结果:
黄色的是筛选前的点,绿色是筛选后的点,绿色的点把黄色的覆盖了一部分。
新闻热点
疑难解答