Synopsis
		sfstphase < in.rsf > out.rsf lambda= window_size=1 verb=n complex=n order=100 ref=1. hertz=n band=n

Parameters
		bool band=n [y/n] if y, compute instantaneous bandwidth   bool complex=n [y/n] if y, use complex-valued computations   bool hertz=n [y/n] if y, convert output to Hertz   float lambda=   int order=100 Hilbert transformer order   float ref=1. Hilbert transformer reference (0.5 < ref <= 1)   bool verb=n [y/n] verbosity   int window_size=1 window_size