Hi All,
I am new to obspy and am trying to plot a PPSD of my miniseed data.
I keep getting an error saying no data accumulated.
I was able to plot the trace no problem.
The code I am currently running can be seen below ( I have added some unnecessary lines for trouble shooting).
from obspy import read
from obspy.signal import PPSD
st = read(“C:\Users\21352643\Desktop\S- Seismic Data\S- Garulp Sensor at South ITM\2017_04_07\01h00m_200_e.msd”)
tr = st.select(id=".CMG6…HHE")[0]
tr.plot()
paz = {‘gain’: 1983000, ‘poles’: [-0.02365+0.02365j, -0.02365-0.02365j, -393.011, -7.4904, -53.5979-21.7494j, -53.5979+21.7494j], ‘sensitivity’: 1259.73, ‘zeros’: [-5.03207, 0, 0]}
ppsd= PPSD(tr.stats, paz)
print(paz)
print(len(tr))
print(tr.stats)
ppsd.add(st)
#st = read(“C:\Users\21352643\Desktop\01h00m_200_e.msd”)
#print(ppsd.add(st))
print(ppsd.times_processed[:10])
print(ppsd.id)
#print(ppsd.id)
ppsd.plot()
The output of this is:

{‘zeros’: [-5.03207, 0, 0], ‘sensitivity’: 1259.73, ‘poles’: [(-0.02365+0.02365j), (-0.02365-0.02365j), -393.011, -7.4904, (-53.5979-21.7494j), (-53.5979+21.7494j)], ‘gain’: 1983000}
720000
network:
station: CMG6
location:
channel: HHE
starttime: 2017-04-07T01:00:00.000000Z
endtime: 2017-04-07T01:59:59.995000Z
sampling_rate: 200.0
delta: 0.005
npts: 720000
calib: 1.0
_format: MSEED
mseed: AttribDict({‘record_length’: 4096, ‘encoding’: u’STEIM1’, ‘filesize’: 1568768, u’dataquality’: u’D’, ‘number_of_records’: 383, ‘byteorder’: u’>’})
[]
.CMG6…HHE
Traceback (most recent call last):
File “”, line 1, in
runfile(‘C:/Users/21352643/untitled0.py’, wdir=‘C:/Users/21352643’)
File “C:\Users\21352643\AppData\Local\Continuum\Anaconda2\lib\site-packages\spyder\utils\site\sitecustomize.py”, line 866, in runfile
execfile(filename, namespace)
File “C:\Users\21352643\AppData\Local\Continuum\Anaconda2\lib\site-packages\spyder\utils\site\sitecustomize.py”, line 87, in execfile
exec(compile(scripttext, filename, ‘exec’), glob, loc)
File “C:/Users/21352643/untitled0.py”, line 18, in
ppsd.plot()
File “C:\Users\21352643\AppData\Local\Continuum\Anaconda2\envs\obspy\lib\site-packages\obspy\signal\spectral_estimation.py”, line 1541, in plot
self.__check_histogram()
File “C:\Users\21352643\AppData\Local\Continuum\Anaconda2\envs\obspy\lib\site-packages\obspy\signal\spectral_estimation.py”, line 758, in __check_histogram
raise Exception(msg)
Exception: No data accumulated
I have the calibration data for the sensor I am using (see attached).
Am I entering this correctly into the paz function?
Thanks all
T6294.doc (114 KB)