ObsPy Python stream图坐标轴限制不生效及指定数据读取问题

核心原因 #

你之前对坐标轴设置的范围、标签没有作用到实际绘制波形的轴上，根源是Obspy Stream.plot()方法的参数逻辑特性：
当向st.plot()传入fig=fig参数时，方法不会复用你在该fig下提前创建的ax对象，而是会在fig内自动生成一个全新的独立子图轴来承载波形。你所有对初始ax的配置，全作用在那个没有被用来绘图的空轴上，因此出现“坐标轴属性只对外框生效、数据在另一区域绘制”的现象。

正确实现方式 #

两种常用的可达到预期效果的写法：

• 直接将自定义绘图参数传入st.plot()方法，无需提前手动创建fig和ax，利用方法原生支持的参数完成配置：

st = obspy.read('Z7*')
st.plot(
    figsize=(14, 4),
    color='blue',
    ylim=(0, 0.01),
    xlabel="My x label",
    title="可自定义的图形标题"
)

• 保留手动创建ax的流程，将绘图目标参数从fig=fig替换为ax=ax。当st.plot()检测到传入了ax对象时，会直接把波形绘制到指定轴上，不会新建独立子图，你提前对ax设置的所有属性都会正常生效：

fig, ax = plt.subplots(figsize=(14, 4))
ax.set_ylim(0, 0.01)
ax.set_xlabel("My x label")
# 核心修改：传入ax=ax，而非fig=fig
st.plot(ax=ax, color='blue')

当前使用的obspy.read('Z7*')是通过通配符匹配当前目录下所有以Z7开头的文件，如果需要仅读取指定范围的文件，可根据场景选择以下方法：

• 精准读取少量指定文件时，直接将需要读取的文件路径组成列表，传入obspy.read()即可，例如仅读取Z7001、Z7002两个文件时，代码为st = obspy.read(["Z7001", "Z7002"])
• 需要批量筛选Z7开头的文件时，可先通过文件匹配逻辑生成符合条件的文件列表，再传入读取方法。例如仅筛选文件名包含2022年标识的Z7类文件：

import glob
matched_files = [f for f in glob.glob("Z7*") if "2022" in f]
st = obspy.read(matched_files)

• 已经完成全量文件读取后，也可直接在内存中对流对象做筛选，无需重复读取磁盘文件。例如仅保留时间范围在2023-01-01至2023-01-02的波形数据：

from obspy import UTCDateTime
st = st.select(
    starttime=UTCDateTime("2023-01-01"),
    endtime=UTCDateTime("2023-01-02")
)

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