Easier, better, faster, stronger
The reflectivity equations in
module now work on arrays for the Vp, Vs, and rho values, as well as the theta values. This is about 10 times faster than running a loop over elements; the Zoeppritz solution is 100× faster.
- The various Zoeppritz solutions and the Aki–Richards approximations now return the complex reflectivity and therefore show post-critical amplitudes correctly.
A new reflection coefficient series function,
, makes it easier to compute offset reflectivities from logs.
Several new linear and non-linear filters are in
(good for seismic horizons),
(good for waveform classification), symmetric nearest-neighbours or
(aka Klauder) and
wavelet now all work for a sequence of frequencies, returning a wavelet bank. Also added a
wavelet, with a taper option to attenuate the sidelobes.
, and other density and velocity transforms, to
(average velocity) and
(naïve Backus average). These all operate in a 'cumulative' average-down-to sense.
- Added a coordinate transformation to translate between arbitrarily oriented (x,y) and (inline, line) coordinates.