У меня есть этот код:
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 пуст.
Спасибо!