Discrete natural neighbour interpolation with uncertainty using cross-validation error-distance fields

Interpolation techniques provide a method to convert point data of a geographic phenomenon into a continuous field estimate of that phenomenon, and have become a fundamental geocomputational technique of spatial and geographical analysts. Natural neighbour interpolation is one method of interpolation that has several useful properties: it is an exact interpolator, it creates a smooth surface free of any discontinuities, it is a local method, is spatially adaptive, requires no statistical assumptions, can be applied to small datasets, and is parameter free. However, as with any interpolation method, there will be uncertainty in how well the interpolated field values reflect actual phenomenon values. Using a method based on natural neighbour distance based rates of error calculated for data points via cross-validation, a cross-validation error-distance field can be produced to associate uncertainty with the interpolation. Virtual geography experiments demonstrate that given an appropriate number of data points and spatial-autocorrelation of the phenomenon being interpolated, the natural neighbour interpolation and cross-validation error-distance fields provide reliable estimates of value and error within the convex hull of the data points.While thismethod does not replace the need for analysts to use sound judgement in their interpolations, for those researchers for whomnatural neighbour interpolation is the best interpolation option the method presented provides a way to assess the uncertainty associated with natural neighbour interpolations. Subjects Computational Science, Spatial and Geographic Information Science

[1]  Donald S. Fussell,et al.  Computer rendering of stochastic models , 1982, Commun. ACM.

[2]  M. Goodchild,et al.  Uncertainty in geographical information , 2002 .

[3]  P. Burrough,et al.  Principles of geographical information systems , 1998 .

[4]  Alan E. Gelfand,et al.  Attaching uncertainty to deterministic spatial interpolations , 2012 .

[5]  Joel Nothman,et al.  SciPy 1.0-Fundamental Algorithms for Scientific Computing in Python , 2019, ArXiv.

[6]  Gaël Varoquaux,et al.  The NumPy Array: A Structure for Efficient Numerical Computation , 2011, Computing in Science & Engineering.

[7]  M. Zeri,et al.  Evaluation of methods of spatial interpolation for monthly rainfall data over the state of Rio de Janeiro, Brazil , 2018, Theoretical and Applied Climatology.

[8]  M. Sambridge,et al.  Geophysical parametrization and interpolation of irregular data using natural neighbours , 1995 .

[9]  C. Tomlin Geographic information systems and cartographic modeling , 1990 .

[10]  C. Willmott,et al.  Advantages of the mean absolute error (MAE) over the root mean square error (RMSE) in assessing average model performance , 2005 .

[11]  David O'Sullivan,et al.  Geographic Information Analysis , 2002 .

[12]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

[13]  Franziska Hoffmann,et al.  Spatial Tessellations Concepts And Applications Of Voronoi Diagrams , 2016 .

[14]  H. M. Yılmaz The effect of interpolation methods in surface definition: an experimental study , 2007 .

[15]  A. McEwen,et al.  Technical note: An evaluation of interpolation methods for Mars Orbiter Laser Altimeter (MOLA) data , 2004 .

[16]  W. Tobler A Computer Movie Simulating Urban Growth in the Detroit Region , 1970 .

[17]  Thomas R. Etherington,et al.  NLMpy: a python software package for the creation of neutral landscape models within a general numerical framework , 2015 .

[18]  R. Sibson,et al.  A brief description of natural neighbor interpolation , 1981 .

[19]  Bernd Hamann,et al.  Discrete Sibson interpolation , 2006, IEEE Transactions on Visualization and Computer Graphics.

[20]  S. Cole,et al.  CEH-GEAR: 1 km resolution daily and monthly areal rainfall estimates for the UK for hydrological and other applications , 2015 .

[21]  Kenji Matsuura,et al.  On the use of dimensioned measures of error to evaluate the performance of spatial interpolators , 2006, Int. J. Geogr. Inf. Sci..

[22]  Terry A. Slocum Thematic Cartography and Visualization , 1998 .

[23]  Dave Watson,et al.  The natural neighbor series manuals and source codes , 1999 .

[24]  C. Frei,et al.  Comparison of six methods for the interpolation of daily, European climate data , 2008 .

[25]  Fernando Pérez,et al.  Python: An Ecosystem for Scientific Computing , 2011, Computing in Science & Engineering.

[26]  Nicholas C. Coops,et al.  Evaluating error associated with lidar-derived DEM interpolation , 2009, Comput. Geosci..