Question

我正在尝试对mseed格式文件进行抽取。现在我正在使用这段代码。

＆＃13;

import numpy as np
import matplotlib.pyplot as plt
import obspy
import sys

with open(sys.argv[1],"rb") as fh:
        fh.seek(512)
        st = obspy.read(fh)


tr = st[0]
tr_new = tr.copy()
tr_new.decimate(factor=5, strict_length=False)
tr_new.write(sys.argv[1] + ".20sps",format="mseed")

tr_filt = tr.copy()
tr_filt.filter('lowpass', freq=0.4 * tr.stats.sampling_rate / 4.0)

t = np.arange(0, tr.stats.npts / tr.stats.sampling_rate, tr.stats.delta)
t_new = np.arange(0, tr_new.stats.npts / tr_new.stats.sampling_rate,
tr_new.stats.delta)

＆＃13;

但是当我运行它时，我收到以下错误。

＆＃13;

Analizando: /home/miguel/Documentos/mseed_decimacion/HHZ.D/caig.ig.hhz.d.2018.107.0000
/usr/lib/python2.7/dist-packages/obspy/io/mseed/core.py:772: UserWarning: The encoding specified in trace.stats.mseed.encoding does not match the dtype of the data.
A suitable encoding will be chosen.
  warnings.warn(msg, UserWarning)

＆＃13;

Answer 1

原始MiniSEED数据很可能存储为整数，并使用仅整数压缩（STEIM）进行压缩。读取数据后，此元信息将与Trace对象一起保留。在抽取过程中，将应用一个过滤器，该过滤器将数据转换为浮点数。然后，当您最后将数据写入MiniSEED时，ObsPy会查看仍具有仅整数压缩方案的元信息，并会意识到它无法使用该压缩方案写入浮点数据。然后回到编写未压缩的浮点MiniSEED，所以这里真的没有什么可担心的。毕竟只是警告，没有错误。

使用python和obspy对文件进行抽取

1 个答案: