现在地震数据的分发通常使用Miniseed格式数据,而如果使用SAC及其他软件处理Miniseed格式地震数据时需要进行格式转化,Python的Obspy可将任意格式转换成其他格式的地震数据。
利用glob可实现对文件夹里所有的mseed地震数据进行格式转换。
# convert miniseed to sac
# miniseed -> sac
import glob
from obspy import read
for mfile in glob.glob('*.mseed', recursive=True):
for (i, tr) in enumerate(read(mfile)):
tr.write(mfile.replace('.mseed', "{}.sac".format(i)))