sfhilbert (2.0-git)
Compute hilbert transform using different methods.

        sfhilbert < in.rsf > out.rsf order=100 ref=1. phase=90. type=
type=t -> time domain convolution
type=f -> frequency domain multiplication and use FFT
type=m -> closed-form design of maximally flat FIR hilbert transformer

int order=100
Hilbert transformer order if type=m
float phase=90.
phase shift (in degrees)
float ref=1.
Hilbert transformer reference (0.5 < ref <= 1) if type=m
string type=
Choosing hilbert transform method, type=t means time domain,
type=f means freqency domain, type=m means using more robust algorithm
the default is type=m