我有这段代码:
prf
我希望将此部分(如下所示)保存到csv文件中。我试图将png更改为csv,并使用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)
,但这并不容易......我是所有这些编码内容的新手。请帮帮忙?
.savetxt
我也尝试过使用熊猫,但它不起作用:
ppsd.plot_temporal([0.1, 1, 10],filename=(data.rsplit(".", 1)[0]) + "_05_PPSD temporal.png",show=False)
没有错误,但csv文件为空。
谢谢!