obspy.signal.tf_misfit.fpg

fpg(st1, st2, dt=0.01, fmin=1.0, fmax=10.0, nf=100, w0=6, norm='global', st2_isref=True, a=10.0, k=1.0)[source]

Frequency-dependent Phase Goodness-of-Fit

See also

[Kristekova2009], Eq.(16)

Parameters:
  • st1 signal 1 of two signals to compare, type numpy.ndarray with shape (number of components, number of time samples) or (number of timesamples, ) for single component data
  • st2 signal 2 of two signals to compare, type and shape as st1
  • dt time step between two samples in st1 and st2
  • fmin minimal frequency to be analyzed
  • fmax maximal frequency to be analyzed
  • nf number of frequencies (will be chosen with logarithmic spacing)
  • w0 parameter for the wavelet, tradeoff between time and frequency resolution
  • norm ‘global’ or ‘local’ normalization of the misfit
  • st2_isref (bool) True if st2 is a reference signal, False if none is a reference
  • a Maximum value of Goodness-of-Fit for perfect agreement
  • k sensitivity of Goodness-of-Fit to the misfit
Returns:

frequency dependent Phase Goodness-of-Fit, type numpy.ndarray with shape (nf,) for single component data and (number of components, nf) for multicomponent data