Я переформатировал кодовую форму этой цепочки , чтобы иметь возможность перебирать список объектов с нежелательной трассировкой и отображать каждую трассу в своем собственном субплоте на одной фигуре.
Нечто странное происходит, когда создаются 4 фигуры, а на «Рис. 4» есть все четыре субплота, на рис. 3 - 3 субплота, нижний - не нанесен, 2 - 2 субплота и т. Д.
Я считаю, что проблема может заключаться в цикле for с zip
, но я не знаю достаточно синтаксиса python на данный момент, чтобы решить эту проблему. Любая помощь будет принята с благодарностью, даже толчок в правильном направлении.
import obspy
st = obspy.read('sac_file.sac') #vertical component
tr = st[0]
index = 0
windowed_tr = []
for windowed in tr.slide(window_length=500, step=100):
windowed_tr.append(windowed)
fig, axes = plt.subplots(nrows=4, ncols=1)
for tr_, ax in zip(windowed_tr, axes.reshape(-1)):
ax.plot(windowed.times(), tr.data, color='k')
print(tr_)
plt.show()