obspy.signal.array_analysis.array_processing

array_processing(stream, win_len, win_frac, sll_x, slm_x, sll_y, slm_y, sl_s, semb_thres, vel_thres, frqlow, frqhigh, stime, etime, prewhiten, verbose=False, coordsys='lonlat', timestamp='mlabday', method=0, store=<function nop at 0x518ce60>)[source]

Method for Seismic-Array-Beamforming/FK-Analysis/Capon

Parameters:
  • stream Stream object, the trace.stats dict like class must contain a obspy.core.util.AttribDict with ‘latitude’, ‘longitude’ (in degrees) and ‘elevation’ (in km), or ‘x’, ‘y’, ‘elevation’ (in km) items/attributes. See param coordsys
  • win_len (Float) Sliding window length in seconds
  • win_frac (Float) Fraction of sliding window to use for step
  • sll_x (Float) slowness x min (lower)
  • slm_x (Float) slowness x max
  • sll_y (Float) slowness y min (lower)
  • slm_y (Float) slowness y max
  • sl_s (Float) slowness step
  • semb_thres (Float) Threshold for semblance
  • vel_thres (Float) Threshold for velocity
  • frqlow (Float) lower frequency for fk/capon
  • frqhigh (Float) higher frequency for fk/capon
  • stime (UTCDateTime) Starttime of interest
  • etime (UTCDateTime) Endtime of interest
  • prewhiten (int) Do prewhitening, values: 1 or 0
  • coordsys valid values: ‘lonlat’ and ‘xy’, choose which stream attributes to use for coordinates
  • timestamp (string) valid values: ‘julsec’ and ‘mlabday’; ‘julsec’ returns the timestamp in secons 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 (int) the method to use 0 == bf, 1 == capon
  • store (int) a function which is called on each iteration with the relative power map and the time offset as argument. Usefull for storing or plotting the map for each iteration. For this purpose the dump and nop function of this module can be used.
Returns:

numpy.ndarray of timestamp, relative relpow, absolute relpow, backazimut, slowness

This Page