obspy.signal.polarization.polarization_analysis¶
- polarization_analysis(stream, win_len, win_frac, frqlow, frqhigh, stime, etime, verbose=False, timestamp=u'mlabday', method=u'pm', var_noise=0.0)[source]¶
Method carrying out polarization analysis with the Flinn, Jurkevics, ParticleMotion, or Vidale 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 for PM
- frqhigh (float) higher frequency for PM
- stime (obspy.core.utcdatetime.UTCDateTime) Start time of interest
- etime (obspy.core.utcdatetime.UTCDateTime) End time of interest
- timestamp (str) valid values: "julsec" and "mlabday"; "julsec" returns the timestamp in seconds since 1970-01-01T00:00:00, "mlabday" returns the timestamp in days (decimals represent hours, minutes and seconds) since 0001-01-01T00:00:00 as needed for matplotlib date plotting (see e.g. matplotlibs num2date)
- method (str) the method to use. one of "pm", "flinn" or "vidale".
Returns: Dictionary with azimuth, incidence angle, errors, rectilinearity, planarity, and/or ellipticity (the returned values depend on the used method).