This paper is review pending but the review hasn't started. Editor and reviewer assignments are happening over in GitHub over here »
xmsinterp is a spatial interpolation library. It is used in commercial applications (www.aquaveo.com - GMS, SMS, WMS) as well as research.
The library supports both 2D (x,y) and 3D (x,y,z) interpolation. Supported 2d interpolation methods include: linear, natural neighbor, clough tocher, and inverse distance weighted (IDW). IDW also supports 3D interpolation. IDW and natural neighbor support modeling data trends using gradient plane and quadratic nodal functions.
Example uses of the library for c++ can be found in the auto-generated documentation here: https://aquaveo.github.io/xmsinterp/md__interpolation__tutorial.html
The library also has a python distribution available through anaconda at https://anaconda.org/Aquaveo/xmsinterp.
An jupyter notebook demonstrating the python api can be found in the github repository in the examples folder and here https://aquaveo.github.io/examples/xmsinterp/2D%20Interpolation.html