An interpolation method for multivariate data in a Delaunay triangulation. The value for an interpolation point is estimated using weighted values of the closest surrounding points in the triangulation. These points, the natural neighbors, are the ones the interpolation point would connect to if inserted into the triangulation.