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).
code
more code
~~~~