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.