PPSD workflow

Hi,
Does the PPSD class follow the same workflow as that described in the PQLX manual? i.e., is the day long file truncated to the next lowest power of two?

“First, to significantly improve the Fast Fourier Transform (FFT) speed ratio, by reducing the number of operations, the number of samples in the time series, N, is truncated to the next lowest power of two, 217, leaving N=131,072, thereby reducing the series length such that Th=3276.8s. Second, in order to further reduce the variance of the final PSD estimates, each roughly 1-hour time series record is divided into 13 segments, overlapping by 75 percent, where the length of each new time series segment is now, Tr=Th/4=819.2s with N=32,768=215.”

Thanks!

Januka

Hi Januka,

the approach very closely follows the processing workflows described in
the references papers / technical reports by McNamara et al.
http://docs.obspy.org/citations.html?highlight=mcnamara

By default PPSD works on 1-hour slices (like PQLX), but that can be
changed. For details please see the source code:

https://github.com/obspy/obspy/blob/4c83dfc6e4a6ae4c93c8681349661d67d48505bc/obspy/signal/spectral_estimation.py#L903-L1001

As you can see in there the psd is reused from matplotlib:
https://matplotlib.org/3.1.1/api/mlab_api.html#matplotlib.mlab.psd

cheers,
T

Thanks as always Tobias!

Januka Attanayake

Research Fellow | Earthquake Seismology
Homepage: https://sites.google.com/site/janukaattanayake

School of Earth Sciences | McCoy Bldg. 200
University of Melbourne | Parkville 3010 VIC
Australia

Hi Januka,

the approach very closely follows the processing workflows described in
the references papers / technical reports by McNamara et al.
http://docs.obspy.org/citations.html?highlight=mcnamara

By default PPSD works on 1-hour slices (like PQLX), but that can be
changed. For details please see the source code:

https://github.com/obspy/obspy/blob/4c83dfc6e4a6ae4c93c8681349661d67d48505bc/obspy/signal/spectral_estimation.py#L903-L1001

As you can see in there the psd is reused from matplotlib:
https://matplotlib.org/3.1.1/api/mlab_api.html#matplotlib.mlab.psd

cheers,
T