Synopsis
		sfsphase < in.rsf > out.rsf niter=100 order=10 ref=1. rect#=(1,1,...)

Parameters
		int niter=100 number of iterations   int order=10 Hilbert transformer order   int rect#=(1,1,...) smoothing radius on #-th axis   float ref=1. Hilbert transformer reference (0.5 < ref <= 1)