主要参考网址见:https://seiscode.iris.washington.edu/projects/ws-fetch-scripts/wiki
首先利用FetchData获取metadata文件:
FetchData -N YP -S NE56 -C BH* -s 2011-01-01T00:00:00 -e 2011-01-01T01:00:00 -o YP.mseed -m YP.metadata
-N指台网名;-S台站名; -C是通道名(一般我们用三分量数据,所以用通配符即可);-s和-e分别是开始和结束时间;-o是需要获取的地震事件,一般直接从Jweed上直接下载的,所以大概率不用;-m 就是我们将metadata写入的文件。
然后利用mseed2sac给.mseed解压、转成sac格式、添加台站信息:
mseed2sac YP.mseed -m YP.metadata
接着再为原始数据添加事件的头段变量即可——利用自己写的 Addeqheaders_v1.1.pl 脚本