sflinear performs 1-D linear interpolation of irregularly spaced data.

The following example from rsf/su/rsflab4 shows a linearly interpolated velocity profile:

The input to sflinear contains coordinate-value pairs arranged so that the second dimension of the data is n2=2. The output contains regularly sampled values on the specified grid.

If the input coordinates are not in order and need sorting, use sort=y.

The output grid can be specified either by supplying it in a pattern file pattern= or by specifying the usual parameters n1=, o1=, d1=.

If the number of iterations specified by niter= is greater than zero, sflinear switches from simple linear interpolation to iterative interpolation by shaping regularization, which can produce a smooth output. The additional parameters to control this process are nw= (size of the local Lagrange interpolation filter for forward interpolation) and rect= (smoothing radius for shaping).

10 previous programs of the month: