obspy.signal.polarization.vidale_adapt

vidale_adapt(stream, noise_thres, fs, flow, fhigh, spoint, stime, etime, adaptive=True)[source]

Adaptive window polarization analysis after [Vidale1986] with the modification of adapted analysis window estimated by estimating the instantaneous frequency. It returns the azimuth, incidence, rectilinearity planarity and ellipticity.

Parameters:
  • stream (Stream) ZNE containing trace data
  • noise_thres (float) Variance of the noise sphere; data points are excluded when falling within the sphere of radius sqrt(noise_thres), Default = 0
  • fs (float) sampling rate
  • flow (float) lower frequency limit for analysis
  • fhigh (float) upper frequency limit for analysis
  • spoint (numpy.ndarray) array with traces’ individual start times in samples
  • stime (UTCDateTime) start time of the analysis
  • etime (UTCDateTime) end time for the analysis
  • 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).
Returns:

list of tuples containing azimuth, incidence, rectilinearity, planarity, and ellipticity