Integrated trace in attributes

[EvanBianco]

A function that computes integrated trace of post-stack seismic amplitude data. Vectorized of course.

[kwinkunks]

Needs to have a window argument.

https://agilescientific.com/blog/2011/12/28/i-is-for-integrated-trace.html

Description of Landmark's old implementation:

Integration integrates the trace using Fourier transform. This is done by transforming the trace, dividing each frequency by (2Pifreq*sqrt(- 1)), and zeroing the zero frequency component (removing dc). This gives a result similar to the simple (trace(i) = trace(i) + trace (i-1)), but has a more accurate 90 degree phase response at all frequencies. The frequency domain method is also not dependent on dc biases caused by small amplitude variations of the first few samples.
Integration produces an output trace in which each output sample value is the sum of the original samples, including this original sample.

Another description here, from CREWES, look for BLIMP: https://crewes.org/ForOurSponsors/ResearchReports/1996/1996-21.pdf