wwa2psd (pyleoclim.utils.wavelet.wwa2psd)

pyleoclim.utils.wavelet.wwa2psd(wwa, ts, Neffs, freq=None, Neff=3, anti_alias=False, avgs=2)[source]

Return the power spectral density (PSD) using the weighted wavelet amplitude (WWA).

Parameters

  • wwa (array) – the weighted wavelet amplitude.

  • ts (array) – the time points, should be pre-truncated so that the span is exactly what is used for wwz

  • Neffs (array) – the matrix of effective number of points in the time-scale coordinates obtained from wwz

  • freq (array) – vector of frequency from wwz

  • Neff (int) – the threshold of the number of effective samples

  • anti_alias (bool) – whether to apply anti-alias filter

  • avgs (int) – flag for whether spectrum is derived from instantaneous point measurements (avgs<>1) OR from measurements averaged over each sampling interval (avgs==1)

Returns

psd – power spectral density

Return type

array

References

Kirchner’s C code for weighted psd calculation