SUSPECK1K2 - 2D (K1,K2) Fourier SPECtrum of (x1,x2) data set

suspeck1k2 <infile >outfile [optional parameters]

Optional parameters:

d1=from header(d1) or 1.0 spatial sampling interval in first (fast)
dimension
d2=from header(d2) or 1.0 spatial sampling interval in second
(slow) dimension

verbose=0 verbose = 1 echoes information

tmpdir= if non-empty, use the value as a directory path
prefix for storing temporary files; else if the
the CWP_TMPDIR environment variable is set use
its value for the path; else use tmpfile()

Notes:
Because the data are assumed to be purely spatial (i.e. non-seismic),
the data are assumed to have trace id (30), corresponding to (z,x) data

To facilitate further processing, the sampling intervals in wavenumber
as well as the first frequency (0) and the first wavenumber are set in
the output header (as respectively d1, d2, f1, f2).

The relation: w = 2 pi F is well known for frequency, but there
doesn't seem to be a commonly used letter corresponding to F for the
spatial conjugate transform variables. We use K1 and K2 for this.
More specifically we assume a phase:
-i(k1 x1 + k2 x2) = -2 pi i(K1 x1 + K2 x2).
and K1, K2 define our respective wavenumbers.


Credits:
CWP: John Stockwell, 26 April 1995, based on original code by
Dave Hale and Jack Cohen

Trace header fields accessed: ns, d1, d2, trid
Trace header fields modified: tracl, ns, dt, trid, d1, f1, d2, f2