Stochastic Local Interaction Model: Geostatistics without Kriging

Classical geostatistical methods face serious computational challenges if they are confronted with large sets of spatially distributed data. We present a simplified stochastic local interaction (SLI) model for computationally efficient spatial prediction that can handle large data. The SLI method constructs a spatial interaction matrix (precision matrix) that accounts for the data values, their locations, and the sampling density variations without user input. We show that this precision matrix is strictly positive definite. The SLI approach does not require matrix inversion for parameter estimation, spatial prediction, and uncertainty estimation, leading to computational procedures that are significantly less intensive computationally than kriging. The precision matrix involves compact kernel functions (spherical, quadratic, etc.) which enable the application of sparse matrix methods, thus improving computational time and memory requirements. We investigate the proposed SLI method with a data set that includes approximately 11500 drill-hole data of coal thickness from Campbell County (Wyoming, USA). We also compare SLI with ordinary kriging (OK) in terms of estimation performance, using cross validation analysis, and computational time. According to the validation measures used, SLI performs slightly better in estimating seam thickness than OK in terms of cross-validation measures. In terms of computation time, SLI prediction is 3 to 25 times (depending on the size of the kriging neighborhood) faster than OK for the same grid size.

READ FULL TEXT

page 15

page 24

page 25

page 27

page 31

research
01/16/2015

Stochastic Local Interaction (SLI) Model: Interfacing Machine Learning and Geostatistics

Machine learning and geostatistics are powerful mathematical frameworks ...
research
02/20/2019

Stochastic Local Interaction Model with Sparse Precision Matrix for Space-Time Interpolation

The application of geostatistical and machine learning methods based on ...
research
02/16/2018

Bayesian cross-validation of geostatistical models

The problem of validating or criticising models for georeferenced data i...
research
03/13/2013

Estimation Stability with Cross Validation (ESCV)

Cross-validation (CV) is often used to select the regularization paramet...
research
06/09/2019

A Low Rank Gaussian Process Prediction Model for Very Large Datasets

Spatial prediction requires expensive computation to invert the spatial ...
research
01/06/2021

Cross-Validation and Uncertainty Determination for Randomized Neural Networks with Applications to Mobile Sensors

Randomized artificial neural networks such as extreme learning machines ...
research
03/23/2022

Application of Co-Kriging and Ordinary Kriging for Selecting Additional Well Locations

Well performance index (WPI), which is an indicator on how much producin...

Please sign up or login with your details

Forgot password? Click here to reset