obspy.signal.polarization.polarization_analysis

polarization_analysis(stream, win_len, win_frac, frqlow, frqhigh, stime, etime, verbose=False, method='pm', var_noise=0.0, adaptive=True)[source]

Method carrying out polarization analysis with the [Flinn1965b], [Jurkevics1988], ParticleMotion, or [Vidale1986] algorithm.

Parameters:
  • stream (Stream) 3 component input data.
  • win_len (float) Sliding window length in seconds.
  • win_frac (float) Fraction of sliding window to use for step.
  • var_noise (float) resembles a sphere of noise in PM where the 3C is excluded
  • frqlow (float) lower frequency. Only used for method='vidale'.
  • frqhigh (float) higher frequency. Only used for method='vidale'.
  • stime (obspy.core.utcdatetime.UTCDateTime) Start time of interest
  • etime (obspy.core.utcdatetime.UTCDateTime) End time of interest
  • method (str) the method to use. one of "pm", "flinn" or "vidale".
  • adaptive (bool) switch for adaptive window estimation (defaults to True). If set to False, the window will be estimated as 3 * max(1/(fhigh-flow), 1/flow).
Return type:

dict

Returns:

Dictionary with keys "timestamp" (POSIX timestamp, can be used to initialize UTCDateTime objects), "azimuth", "incidence" (incidence angle) and additional keys depending on used method: "azimuth_error" and "incidence_error" (for method "pm"), "rectilinearity" and "planarity" (for methods "flinn" and "vidale") and "ellipticity" (for method "flinn"). Under each key a ndarray is stored, giving the respective values corresponding to the "timestamp" ndarray.