Synopsis
		sfpick < scn.rsf > pik.rsf vel0=o2 niter=100 an=1. gate=3 smooth=y rect#=(1,1,...) rect1=1 rect2=1 ...

rectN defines the size of the smoothing stencil in N-th dimension.

Theory in Appendix B of:
S. Fomel, 2009,
Velocity analysis using AB semblance: Geophysical Prospecting, v. 57, 311-321.
Reproducible version in RSFSRC/book/jsg/avo
http://ahay.org/RSF/book/jsg/avo/paper_html/

August 2012 program of the month:
http://ahay.org/rsflog/index.php?/archives/302-Program-of-the-month-sfpick.html

Parameters
		float an=1. axes anisotropy   int gate=3 picking gate   int niter=100 number of iterations   int rect#=(1,1,...) smoothing radius on #-th axis   bool smooth=y [y/n] if apply smoothing   float vel0=o2 surface velocity

Used In
		BEI          vela/vscan wvs/vscan   DATA          nankai/shots   GEE          prc/bob   GEO384W          hw3/synth hw6/blake   JSG          attr/vecta avo/avo avo/avo2 diffr/fault diffr/gom ltft/vecta oc3/bei oc3/fault pmig/mig reversible/nscnmo seislet/diplet seislet/gath shape/beivc simistack/bei simistack/flat4 timelapse/duri timelapse/long1 timelapse/timelapse   MILANO          taupvel/cmp   RSF          reprobook/exam3 scons/rsf usp/data   SEP          oc2/elfshot vc2/beivc vc2/elfvc vc2/sigvc