Python - 如何使用 obspy 模块将计数转换为 m/s

Question

我有一个带有单 channel 跟踪的小型种子文件，我假设数据是计数的(我如何检查跟踪的单位？)。我需要将其转换为 m/s。我已经检查了 obspy 教程，我的主要问题是我不知道如何从小型种子文件访问极点和零点以及放大系数。另外，我需要校准文件吗？

这是我的代码:

from obspy.core import *
st=read('/Users/guilhermew/Documents/Projecto/Dados sismicos 1 dia/2012_130_DOC01.mseed')
st.plot()

提前致谢，吉列尔梅

编辑:我终于明白如何转换数据了。 Obspy 有不同的方法来实现这一点，但这一切都归结为从波形数据中消除仪器响应。就像 @Robert Barsch 所说，我需要另一个文件来获取仪器响应元数据。所以我想出了以下代码:

parser=Parser("dir/parser/file")
for tr in stream_aux:
    stream_id=tr.stats.network+'.'+tr.stats.station+ '..' + tr.stats.channel
    paz=parser.getPAZ(stream_id, tr.stats.starttime)
    df = tr.stats.sampling_rate
    tr.data = seisSim(tr.data, df, paz_remove=paz)

我使用 seisSim 函数来转换数据。我现在的问题是输出看起来不正确(但我似乎无法发布图像)

Answer 1

这显然是一个应该问地震学社区而不是 StackOverflow 的问题!您写信给 ObsPy user mailinglist 怎么样？ ？

更新:我仍然认为答案是他/她应该直接在 ObsPy 邮件列表中询问。然而，为了对实际问题给出正确的答案:MiniSEED 是一种纯数据格式，不包含任何元信息，例如极点和零点或使用的单位。所以，是的，您将需要另一个文件，例如 RESP、SAC PAZ、Dataless SEED、Full SEED 等，以便获取电台特定的元数据。要应用地震仪校正，请阅读 http://docs.obspy.org/tutorial/code_snippets/seismometer_correction_simulation.html