Сохранить ppsd-график в виде csv-файла в Python - PullRequest
0 голосов
/ 27 мая 2018

У меня есть этот код:

from obspy.core import read
from obspy.io.xseed import Parser
from obspy.signal import PPSD
from obspy.imaging.cm import pqlx
import glob

#read mseed data
for data in glob.glob("F:/SKRIPSI/MSEED/*.mseed"):
    st=read(data)
    tr=st.select()[0]

    #read dataless
    for dataless in glob.glob("F:/SKRIPSI/MSEED/*.dseed"):
        if (data.rsplit(".", 1)[0])==(dataless.rsplit(".", 1)[0]):
            parser=Parser(dataless)

            #plot ppsd
            ppsd=PPSD(tr.stats, metadata=parser)
            ppsd.add(st)
            ppsd.plot(filename=(data.rsplit(".", 1)[0]) + "_01_PPSD.png")
            ppsd.plot(filename=(data.rsplit(".", 1)[0]) + "_02_PPSD cumulative.png", cumulative=True)
            ppsd.plot(filename=(data.rsplit(".", 1)[0]) + "_03_PPSD pqlx.png", cmap=pqlx)
            ppsd.plot_spectrogram(filename=(data.rsplit(".", 1)[0]) + "_04_PPSD spectrogram.png",show=False)
            ppsd.plot_temporal([0.1, 1, 10],filename=(data.rsplit(".", 1)[0]) + "_05_PPSD temporal.png",show=False)

, и я хочу сохранить эту (как показано ниже) часть графика в CSV-файлы.Я пытался изменить png на csv и использовал .savetxt, но это не так просто, правда ... Я новичок во всем этом коде.Пожалуйста, помогите?

ppsd.plot_temporal([0.1, 1, 10],filename=(data.rsplit(".", 1)[0]) + "_05_PPSD temporal.png",show=False)

Я тоже пытался использовать панд, но это не работает:

import pandas as pd    

temp=ppsd.plot_temporal([0.1, 1, 10],filename=(data.rsplit(".", 1)[0]) + "_05_PPSD temporal.png",show=False)
csvtemp=pd.DataFrame(temp)
csvtemp.to_csv((data.rsplit(".", 1)[0]) + "_05_PPSD temporal.csv", index=False, header=False)

Нет ошибки, но файл csv пуст.

Спасибо!

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...