Итерация по списку трасс и создание сюжетов на одной фигуре (obspy) - PullRequest
Новая
       15

Итерация по списку трасс и создание сюжетов на одной фигуре (obspy)

0 голосов
/ 08 ноября 2019

Я переформатировал кодовую форму этой цепочки , чтобы иметь возможность перебирать список объектов с нежелательной трассировкой и отображать каждую трассу в своем собственном субплоте на одной фигуре.

Нечто странное происходит, когда создаются 4 фигуры, а на «Рис. 4» есть все четыре субплота, на рис. 3 - 3 субплота, нижний - не нанесен, 2 - 2 субплота и т. Д.

Я считаю, что проблема может заключаться в цикле for с zip, но я не знаю достаточно синтаксиса python на данный момент, чтобы решить эту проблему. Любая помощь будет принята с благодарностью, даже толчок в правильном направлении.

4 subplots

3 subplots

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()
...