使用Python与Obspy对MSEED文件降采样遇编码匹配警告

问题描述 #

我正在尝试对MSEED格式文件进行降采样，使用的Python代码如下：

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)

运行代码时出现了以下警告：

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)

希望能解决这个编码与数据类型不匹配的警告问题。

问题原因 #

这个警告的核心原因很明确：你复制原始trace后，tr_new继承了原始MSEED文件的编码配置（比如常见的16位整数编码），但经过decimate降采样处理后，数据的数值类型（dtype）变成了浮点数，这就导致原本的编码设置和实际数据类型不匹配。ObsPy在写入文件时检测到这个矛盾，就会抛出警告，然后自动选择一个适配当前数据类型的编码。

解决方法 #

有两种简单的方式可以消除这个警告：

方式1：写入时显式指定编码 #

在调用write方法时，直接指定适合处理后数据类型的编码。因为降采样后的 data 是浮点数，你可以指定FLOAT32（32位浮点数）或者FLOAT64（64位浮点数）：

tr_new.write(sys.argv[1] + ".20sps", format="mseed", encoding='FLOAT32')

方式2：重置trace的编码配置 #

处理完数据后，手动清除或重置原始的编码设置，让ObsPy自动匹配合适的编码：

# 处理完降采样后，删除原有的编码设置
del tr_new.stats.mseed.encoding
# 或者将编码设置为None
# tr_new.stats.mseed.encoding = None

# 再执行写入操作
tr_new.write(sys.argv[1] + ".20sps", format="mseed")

额外优化建议 #

你的代码里手动fh.seek(512)的操作其实是多余的——ObsPy的obspy.read()可以直接读取MSEED文件，它会自动处理文件头部分。你可以把文件读取的代码简化为：

st = obspy.read(sys.argv[1])

这样不仅代码更简洁，还能避免手动跳过文件头可能带来的读取错误。

内容的提问来源于stack exchange，提问作者Miguel Sosa