使用.xml文件去仪器响应

使用.xml文件去仪器响应

June 7, 2022
地震学
python, 数据处理, 去仪器响应

使用BQmail申请地震数据后，要对地震数据进行去仪器响应。现在的数据格式为miniseed。

去仪器响应代码 #

代码如下：

import obspy
from obspy import read, read_inventory
import glob
import os


# 搜索mseed列表

mseedlst = glob.glob(os.path.join("*.mseed"))
# 循环去响应
for file in mseedlst:
    st = read('TEPE.YL.mseed')
    # 分割字符串为列表
    strlst = file.split(".")
    # inventory  详尽目录
    # 一般为.xml格式 自动匹配
    inv = read_inventory('IRISDMC-' + strlst[0] + '.YL.xml')
    # 滤波预处理 防止反褶积产生的噪声
    pre_filt = (0.005, 0.006, 15, 19)

    # output
    #        :param output: Output units. One of:

    #            ``"DISP"``
    #                displacement, output unit is meters
    #            ``"VEL"``
    #               velocity, output unit is meters/second
    #            ``"ACC"``
    #               acceleration, output unit is meters/second**2
    #
    st.remove_response(inventory=inv, output="VEL", pre_filt=pre_filt)
    # 写入新文件中
    st.write('norespmseed' + '\\' + strlst[0] + '.YL.mseed')
返回