Natural neighbor interpolation finds the closest subset of input samples to a query point and applies weights to them based on proportionate areas in order to interpolate a value (Sibson, 1981). It is also known as Sibson or "area-stealing" interpolation. Its basic properties are that it's local, using only a subset of samples that surround a query point, and that interpolated heights are guaranteed to be within the range of the samples used. It does not infer trends and will not produce peaks, pits, ridges or valleys that are not already represented by the input samples. The surface passes through the input samples and is smooth everywhere except at locations of the input samples. Breaklines may be used, in the case of TIN to Raster interpolation, to augment the surface, creating linear discontinuities where appropriate such as along roadsides and water bodies. It adapts locally to the structure of the input data, requiring no input from the user pertaining to search radius, sample count, or shape. It works equally well with regularly and irregularly distributed data (Watson, 1992).
The natural neighbors of any point are those associated with neighboring Voronoi (i.e. Thiessen) polygons. Initially, a Voronoi diagram is constructed of all the given points, represented by the olive colored polygons. A new Voronoi polygon, beige color, is then created around the interpolation point (red star). The proportion of overlap between this new polygon and the initial polygons are then used as the weights.