Go to the list of seismic processes.
Go to SIOSEIS introduction.
PROCESS T2F
------- ---
Parameters, alphabetically:
addwb coords etime fftlen fno ftr
lno lprint ltr stime window winlen
Document Date: 1 October 2009
Updates:
Oct. 2008 - Add coords AMPL
PROCESS T2F converts traces from the time domain to the frequency domain.
Each trace is windowed, padded with zeroes to a power of 2, transformed,
and optionally converted to polar coordinates. Traces left in rect-
angular coordinates are in complex form such that the imaginary part of
each frequency immediately follows the real (real-imaginary pairs).
Traces in polar form have the amplitude spectrum in the first half of the
output trace while the wrapped phase spectrum is in the second half of
the trace.
Each output trace will contain a power of 2 + 2 number of samples and is
the number of samples used in the FFT (fftlen+2). There are fftlen/2+1
frequencies in the output trace representing frequencies 0 (dc) to the
Nyquist frequency (1./(2*Sample interval)).
Any seismic process may be used on frequency data, but some post T2F
processes may not make geophysical sense. Frequency domain plots may be
made using polar coordinates and plotting the first half of the trace.
Bandpass filter and deconvolution should be done in the time domain using
processes filter and decon before the frequency domain processes.
A good review of the frequency domain is "A Guided Tour of the Fast
Fourier Transform" by G.D.Bergland in "IEEE Spectrum", July 1969.
THE PARAMETER DICTIONARY
--- --------- ----------
STIME - Start time, in seconds, of the data to transform. The use of
stime causes the the deep water delay to be set to 0. and the
first time sample on an inverse transform (PROCESS F2T) to be 0.
Preset = the delay of the trace. e.g. stime 3.3
ETIME - End time, in seconds, of the data to transform. The data after
etime will be padded with zeroes in order for the data length be
a power of 2.
Preset = the entire data trace. e.g. etime 4.0
ADDWB
WINDOW - The type of window to apply before computing the fft.
= hamm, Hamming
= hann, Hanning
= bart, Bartlett (triangular)
= rect, rectangular (box car - no window)
= blac, Blackman
= ebla, exact Blackman
= blha, Blackman-Harris
Preset=hamm e.g. window rect
WINLEN - The window length, in seconds. A window length of zero causes
the entire time domain gate to be windowed. A non zero length
indicates that winlen data will be modified at both ends of each
data gate.
Preset=0 e.g. winlen .2
COORDS - The coordinates of the output trace.
= rect, rectangular coordinates. The output trace will be complex.
The trace values will consist of real and imaginary pairs. The
use of rectangular coordinates is slightly faster than polar
coordinates, and is understood by frequency domain processes
such as f2t.
= polr, polar coordinates. The first half of the output trace
will be the full (-pi to +pi) amplitude spectrum and the second
half of the trace will be the phase spectrum.
= ampl, amplitude spectrum. Only the positive portion ( 0 to pi )
of the amplitude spectrum is output.
Preset = rect e.g. coords polr
FFTLEN - The number of points to use in the FFT. T2F will make this a
power of 2 if it is not.
Preset = the number of points in the trace. e.g. fftlen 1024
END - Terminates each parameter list.
Written and copyrighted (c) by:
Paul Henkart, Scripps Institution of Oceanography, January 1984
ALL RIGHTS RESERVED.