I have two equidistant vectors (longitude and time) that make up a “grid” in which I want to plot velocity measurements. The velocity data points have corresponding longitudes and times, however not all of those will neccessarily be represented in the equidistant vectors, since I have more datapoints than the length of the equidistant vectors. Is there a function or package which would allow me to do this?
Thanks in advance!
If I understand your use case correctly, ScatteredInterpolation.jl should do the trick. I’ve personally used
scipy.interpolate.LinearNDInterpolator via PythonCall.jl, which uses linear rather than radial basis function interpolation.
The GMT gridit is a wrapper function to the several interpolation methods in GMT+GDAL
It could have a richer manual (that’s the intention) but given that it is a wrapper function, original methods have better documentation.