Synopsis
		sfepfbil < in.rsf > out.rsf twod=n sigmar=1.0 sigmad=1.0

Parameters
		float sigmad=1.0 sigma for domain filter   float sigmar=1.0 sigma for range filter   bool twod=n [y/n] y, 2D smoothing