sfmyradon2 (2.0-git)
Linear/parabolic radon transform frequency domain implementation

        sfmyradon2 < in.rsf > out.rsf adj=y inv=n np= dp= p0= niter=100 eps=0.01 nx= ox= dx= parab=n x0=1. invmode=
Also referred to as high-resolution radon transform
Note: I borrowed a lot from /system/seismic/radon+Mradon.c. The distinction:
I am using FFTW because I am inexperienced in invoking kiss_fft.

bool adj=y [y/n]
if y, perform adjoint operation
float dp=
p sampling (if adj=y)
float dx=
sampling interval in x
float eps=0.01
regularization parameter
bool inv=n [y/n]
if y, perform inverse operation
string invmode=
inverse method: 'ls' if least-squares; 'toeplitz' if use FFT
int niter=100
number of CGLS iterations
int np=
number of p values (if adj=y)
int nx=
number of offsets (if adj=n)
float ox=
x origin
float p0=
p origin (if adj=y)
bool parab=n [y/n]
if y, parabolic Radon transform
float x0=1.
reference offset

Used In