DOCUMENTATION
SYNOPSIS
PERL PROGRAM NAME: SUICEPSTRUM - fft of complex log frequency traces to real time traces
AUTHOR: Juan Lorenzo (Perl module only)
DATE:
DESCRIPTION:
Version:
USE
NOTES
Examples
SEISMIC UNIX NOTES
SUICEPSTRUM - fft of complex log frequency traces to real time traces
suicepstrum <stdin >sdout sign2=-1
Required parameters:
none
Optional parameter:
sign1=1 sign in exponent of first fft
sign2=-1 sign in exponent of inverse fft
sym=0 =1 center output
dt=tr.dt time sampling interval (s) from header
if not set assumed to be .004s
Output traces are normalized by 1/N where N is the fft size.
Note:
The forward cepstral transform is the
F(t_c) = InvFFT[ln[FFT(F(t))]]
The inverse cepstral transform is the
F(t) = InvFFT[exp[FFT(F(t_c))]]
Here t_c is the cepstral time (quefrency) domain
Credits:
CWP: John Stockwell, Dec 2010 based on
suifft.c by:
CWP: Shuki Ronen, Chris Liner, Jack K. Cohen, c. 1989
Trace header fields accessed: ns, trid
Trace header fields modified: ns, trid
CHANGES and their DATES
sub Step
collects switches and assembles bash instructions by adding the program name
sub note
collects switches and assembles bash instructions by adding the program name
sub clear
sub dt
sub sign1
sub sign2
sub sym
sub get_max_index
max index = number of input variables -1