array_processing.algorithms.fk_freq module

array_processing.algorithms.fk_freq.fk_freq(data, fs, rij, vmin, vmax, fmin, fmax, nvel, ntheta)[source]

\(f\)–\(k\) beamforming with loop over frequency bands.

Parameters:
  • data – (m, n) array; time series with m samples from n traces as columns
  • rij – (d, n) array; n sensor coordinates as [northing, easting, {elevation}] column vectors in d dimensions
  • fs (int or float) – Sample rate [Hz]
  • vmin (int or float) – Min velocity in km/s, suggest 0.25
  • vmax (int or float) – Max velocity in km/s, suggest 0.45
  • fmin (int or float) – Minimum frequency in Hz
  • fmax (int or float) – Maximum frequency in Hz
  • nvel (int or float) – Number of velocity iterations, suggest 100–200
  • ntheta (int or float) – Number of azimuth iterations, suggest 100–200
Returns:

(ntheta, nvel) array; beamformed slowness map, not normalized. Can find max using

ix, iy = np.unravel_index(bmpwr.argmax(), bmpwr.shape)