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)