Exploring the World of Ordinary Kriging. Dennis J. J. Walvoort. Wageningen University & Research Center Wageningen, The Netherlands

Similar documents
11/8/2018. Spatial Interpolation & Geostatistics. Kriging Step 1

Spatial Interpolation & Geostatistics

Introduction. Semivariogram Cloud

7 Geostatistics. Figure 7.1 Focus of geostatistics

Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.

Data Break 8: Kriging the Meuse RiverBIOS 737 Spring 2004 p.1/27

Geostatistics: Kriging

Spatial Data Analysis in Archaeology Anthropology 589b. Kriging Artifact Density Surfaces in ArcGIS

Spatiotemporal Analysis of Environmental Radiation in Korea

REN R 690 Lab Geostatistics Lab

Umeå University Sara Sjöstedt-de Luna Time series analysis and spatial statistics

Chapter 4 - Fundamentals of spatial processes Lecture notes

Soil Moisture Modeling using Geostatistical Techniques at the O Neal Ecological Reserve, Idaho

Gridding of precipitation and air temperature observations in Belgium. Michel Journée Royal Meteorological Institute of Belgium (RMI)

Spatial Data Mining. Regression and Classification Techniques

PRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH

An Introduction to Pattern Statistics

Basics of Point-Referenced Data Models

Lecture 5 Geostatistics

Relationships between Soil salinity and geopedological units in Saveh plain, Iran

A robust statistically based approach to estimating the probability of contamination occurring between sampling locations

Lab #3 Background Material Quantifying Point and Gradient Patterns

University of California, Los Angeles Department of Statistics. Effect of variogram parameters on kriging weights

SPATIAL-TEMPORAL TECHNIQUES FOR PREDICTION AND COMPRESSION OF SOIL FERTILITY DATA

I don t have much to say here: data are often sampled this way but we more typically model them in continuous space, or on a graph

OFTEN we need to be able to integrate point attribute information

Geostatistical Interpolation: Kriging and the Fukushima Data. Erik Hoel Colligium Ramazzini October 30, 2011

CREATION OF DEM BY KRIGING METHOD AND EVALUATION OF THE RESULTS

11. Kriging. ACE 492 SA - Spatial Analysis Fall 2003

Example: Identification

A MultiGaussian Approach to Assess Block Grade Uncertainty

University of California, Los Angeles Department of Statistics. Geostatistical data: variogram and kriging

ENVIRONMENTAL MONITORING Vol. II - Geostatistical Analysis of Monitoring Data - Mark Dowdall, John O Dea GEOSTATISTICAL ANALYSIS OF MONITORING DATA

Beers Law Instructor s Guide David T. Harvey

Understanding Anisotropy Computations 1

SPATIAL ELECTRICAL LOADS MODELING USING THE GEOSTATISTICAL METHODS

Beta-Binomial Kriging: An Improved Model for Spatial Rates

Investigation of Monthly Pan Evaporation in Turkey with Geostatistical Technique

Point-Referenced Data Models

What s for today. All about Variogram Nugget effect. Mikyoung Jun (Texas A&M) stat647 lecture 4 September 6, / 17

Improving Spatial Data Interoperability

Porosity prediction using cokriging with multiple secondary datasets

The Proportional Effect of Spatial Variables

Correcting Variogram Reproduction of P-Field Simulation

GEOSTATISTICAL ANALYSIS OF SPATIAL DATA. Goovaerts, P. Biomedware, Inc. and PGeostat, LLC, Ann Arbor, Michigan, USA

Chem 1 Kinetics. Objectives. Concepts

A kernel indicator variogram and its application to groundwater pollution

BD1. Reciprocal space and Fourier transforms. IB Mineral Sciences Module B: Reciprocal Space, Symmetry and Crystallography

Risk assessment: the global confidence uncertainty plume of SAKWeb

An Introduction to Spatial Autocorrelation and Kriging

The ProbForecastGOP Package

Introduction to Geostatistics

Applying MapCalc Map Analysis Software

Mapping Precipitation in Switzerland with Ordinary and Indicator Kriging

Influence of incorporating geometric anisotropy on the construction of thematic maps of simulated data and chemical attributes of soil

Conditional Distribution Fitting of High Dimensional Stationary Data

Practical 12: Geostatistics

Chapter 1. Summer School GEOSTAT 2014, Spatio-Temporal Geostatistics,

Activity 8b - Electric Field Exploration

ArcGIS for Geostatistical Analyst: An Introduction. Steve Lynch and Eric Krause Redlands, CA.

Optimal Interpolation

SIMULATION STUDY OF SPATIAL- POISSON DATA ASSESSING INCLUSION OF SPATIAL CORRELATION AND NON- NORMALITY IN THE ANALYSIS

through any three given points if and only if these points are not collinear.

Spatial Analysis II. Spatial data analysis Spatial analysis and inference

Introduction to applied geostatistics. Short version. Overheads

Toward an automatic real-time mapping system for radiation hazards

Index. Geostatistics for Environmental Scientists, 2nd Edition R. Webster and M. A. Oliver 2007 John Wiley & Sons, Ltd. ISBN:

Gstat: multivariable geostatistics for S

Linear inverse Gaussian theory and geostatistics a tomography example København Ø,

Space-time data. Simple space-time analyses. PM10 in space. PM10 in time

Student Exploration: Energy of a Pendulum

Influence of parameter estimation uncertainty in Kriging: Part 2 Test and case study applications

Lognormal Measurement Error in Air Pollution Health Effect Studies

Types of Spatial Data

Chapter 5: Ordinary Least Squares Estimation Procedure The Mechanics Chapter 5 Outline Best Fitting Line Clint s Assignment Simple Regression Model o

J. Bio. & Env. Sci. 2014

Evaluation of spatial heterogeneity for the design and layout of experimental sites

LOGNORMAL ORDINARY KRIGING METAMODEL

The purpose of this lab is to investigate phases of matter, temperature, and heat energy.

Application and evaluation of universal kriging for optimal contouring of groundwater levels

Geog 210C Spring 2011 Lab 6. Geostatistics in ArcMap

A Short Note on the Proportional Effect and Direct Sequential Simulation

Fluvial Variography: Characterizing Spatial Dependence on Stream Networks. Dale Zimmerman University of Iowa (joint work with Jay Ver Hoef, NOAA)

9.1 - Systems of Linear Equations: Two Variables

Scatterplots. 3.1: Scatterplots & Correlation. Scatterplots. Explanatory & Response Variables. Section 3.1 Scatterplots and Correlation

Spatial Variability of Penetration Resistance on Pseudogley

Gas Pressure and Temperature Relationships *

Rick Faber CE 513 Geostatistical Analyst Lab # 6 6/2/06

Simulating Spatial Distributions, Variability and Uncertainty of Soil Arsenic by Geostatistical Simulations in Geographic Information Systems

Practical application of drill hole spacing analysis in coal resource estimation

A Geostatistical Approach to Predict the Average Annual Rainfall of Bangladesh

Erik W. Gjesfjeld UNIVERSITY COLLEGE LONDON INSTITUTE OF ARCHAEOLOGY

Advanced analysis and modelling tools for spatial environmental data. Case study: indoor radon data in Switzerland

Solar Power Interpolation and Analysis using Spatial Autocorrelation

SPATIAL VARIABILITY MAPPING OF N-VALUE OF SOILS OF MUMBAI CITY USING ARCGIS

2.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics. References - geostatistics. References geostatistics (cntd.

Stepwise Conditional Transformation for Simulation of Multiple Variables 1

Harmonic Motion. Mass on a Spring. Physics 231: General Physics I Lab 6 Mar. 11, Goals:

adaptive prior information and Bayesian Partition Modelling.

BAYESIAN MODEL FOR SPATIAL DEPENDANCE AND PREDICTION OF TUBERCULOSIS

Transcription:

Exploring the World of Ordinary Kriging Wageningen University & Research Center Wageningen, The Netherlands July 2004 (version 0.2)

What is? What is it about? Potential Users a computer program for exploring ordinary kriging; interactive, self-explanatory, and easy to use. Hence, easy kriging and last but not least: it s freeware!

What is? What is it about? Potential Users a computer program for exploring ordinary kriging; interactive, self-explanatory, and easy to use. Hence, easy kriging and last but not least: it s freeware!

What is? What is it about? Potential Users a computer program for exploring ordinary kriging; interactive, self-explanatory, and easy to use. Hence, easy kriging and last but not least: it s freeware!

Who may benefit from it? What is it about? Potential Users 1 Students, who need to understand those mysterious kriging equations; 2 Lecturers, who want to explain kriging in an intuitive way; 3 Others, who just want to know what they are doing when using geostatistical software.

Who may benefit from it? What is it about? Potential Users 1 Students, who need to understand those mysterious kriging equations; 2 Lecturers, who want to explain kriging in an intuitive way; 3 Others, who just want to know what they are doing when using geostatistical software.

Who may benefit from it? What is it about? Potential Users 1 Students, who need to understand those mysterious kriging equations; 2 Lecturers, who want to explain kriging in an intuitive way; 3 Others, who just want to know what they are doing when using geostatistical software.

(GUI) Data Configuration Panel Semivariogram Panel Kriging Panel The GUI consists of three panels: 1 data configuration panel 2 semivariogram panel 3 kriging panel

(GUI) Data Configuration Panel Semivariogram Panel Kriging Panel The GUI consists of three panels: 1 data configuration panel 2 semivariogram panel 3 kriging panel

(GUI) Data Configuration Panel Semivariogram Panel Kriging Panel The GUI consists of three panels: 1 data configuration panel 2 semivariogram panel 3 kriging panel

(GUI) Data Configuration Panel Semivariogram Panel Kriging Panel The GUI consists of three panels: 1 data configuration panel 2 semivariogram panel 3 kriging panel

Data Configuration Panel Data Configuration Panel Semivariogram Panel Kriging Panel The prediction point is red and labelled 0.

Data Configuration Panel Data Configuration Panel Semivariogram Panel Kriging Panel The n sampling points are blue and labelled from 1 to n.

Data Configuration Panel Data Configuration Panel Semivariogram Panel Kriging Panel The configuration of points can be changed by means of the mouse (drag n drop).

Data Configuration Panel Data Configuration Panel Semivariogram Panel Kriging Panel The range is represented by a solid circle. A dashed circle represents the practical range. definition

Data Configuration Panel Data Configuration Panel Semivariogram Panel Kriging Panel The number of sampling points has to be specified here (see ellipse).

Data Configuration Panel Data Configuration Panel Semivariogram Panel Kriging Panel The values at the sampling locations have to be entered here.

Semivariogram Panel Data Configuration Panel Semivariogram Panel Kriging Panel The following semivariogram models can be selected: Spherical model details Exponential model details Gaussian model details

Semivariogram Panel Data Configuration Panel Semivariogram Panel Kriging Panel Use the sliders to set the semivariogram parameters: nugget (c 0 ) partial sill (c 1 ) range parameter (a)

Semivariogram Panel Data Configuration Panel Semivariogram Panel Kriging Panel This section gives the nugget variance, the (practical) range and the sill variance (c 0 + c 1 ).

Semivariogram Panel Data Configuration Panel Semivariogram Panel Kriging Panel The resulting semivariogram will be shown here.

Kriging Panel Introduction Data Configuration Panel Semivariogram Panel Kriging Panel The ordinary kriging weights are shown as bars. Click on the bars for numerical values.

Kriging Panel Introduction Data Configuration Panel Semivariogram Panel Kriging Panel Use the radiobuttons to switch between ordinary point kriging and ordinary block kriging.

Kriging Panel Introduction Data Configuration Panel Semivariogram Panel Kriging Panel The block size can be set by means of the sliders.

Kriging Panel Introduction Data Configuration Panel Semivariogram Panel Kriging Panel Check this checkbox to enforce square prediction blocks.

Kriging Panel Introduction Data Configuration Panel Semivariogram Panel Kriging Panel The ordinary kriging prediction and the associated variance of the prediction error are given here.

Kriging Panel Introduction Data Configuration Panel Semivariogram Panel Kriging Panel Press this button to get a glimpse of the underlying maths.

Effects of sampling point configuration Effect of semivariogram Distance effect Discover that points outside the range affect predictions differently than points within the range (cf. inverse squared distance interpolation).

Effects of sampling point configuration Effect of semivariogram Declustering effect Discover how ordinary kriging reduces the influence of clustered sampling points (cf. inverse squared distance interpolation).

Effects of sampling point configuration Effect of semivariogram Screening effect Recall that ordinary kriging is a non-convex interpolator, i.e., its predictions can be outside the data range. Explore data configurations and semivariogram settings that enhance this effect.

Effects of sampling point configuration Effect of semivariogram Effect of data values Discover how data values affect the weights, the prediction and the variance of the prediction error.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram model Study how model shape, nugget, sill and range affect the kriging results.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram model Study how model shape, nugget, sill and range affect the kriging results.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram model Study how model shape, nugget, sill and range affect the kriging results.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram model Study how model shape, nugget, sill and range affect the kriging results.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram model Study how model shape, nugget, sill and range affect the kriging results.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram scale Explore how multiplying the semivariogram by a positive value affects the kriging weights, the kriging prediction, and the variance of the prediction error.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram scale Explore how multiplying the semivariogram by a positive value affects the kriging weights, the kriging prediction, and the variance of the prediction error.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram scale Explore how multiplying the semivariogram by a positive value affects the kriging weights, the kriging prediction, and the variance of the prediction error.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram scale Explore how multiplying the semivariogram by a positive value affects the kriging weights, the kriging prediction, and the variance of the prediction error.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram scale Explore how multiplying the semivariogram by a positive value affects the kriging weights, the kriging prediction, and the variance of the prediction error.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram scale Explore how multiplying the semivariogram by a positive value affects the kriging weights, the kriging prediction, and the variance of the prediction error.

Effects of sampling point configuration Effect of semivariogram Effect of semivariogram scale Explore how multiplying the semivariogram by a positive value affects the kriging weights, the kriging prediction, and the variance of the prediction error.

Effects of sampling point configuration Effect of semivariogram Experiment with block size and see how it affects the kriging results.

Effects of sampling point configuration Effect of semivariogram Experiment with block size and see how it affects the kriging results.

Effects of sampling point configuration Effect of semivariogram Experiment with block size and see how it affects the kriging results.

Effects of sampling point configuration Effect of semivariogram Experiment with block size and see how it affects the kriging results.

Effects of sampling point configuration Effect of semivariogram Experiment with block size and see how it affects the kriging results.

Effects of sampling point configuration Effect of semivariogram Experiment with block size and see how it affects the kriging results.

Effects of sampling point configuration Effect of semivariogram Experiment with block size and see how it affects the kriging results.

Suggestions for further reading Further reading Geostatistical software License E.H. Isaaks and R.M. Srivastava. An Introduction to Applied Geostatistics. Oxford University Press, New York, 1989. P. Goovaerts. Geostatistics for Natural Resources Evaluation. Oxford University Press, New York, 1997.

Geostatistical software Further reading Geostatistical software License Vesper (www.usyd.edu.au/su/agric/acpa/vesper/vesper.html) GStat (www.gstat.org), GStat is also available as R-package (www.r-project.org)

Further reading Geostatistical software License License is freeware and provided as is without warranty of any kind, either express or implied.

Further reading Geostatistical software License Enjoy!

Definitions Spherical model γ s (h) = { 3h 2a 1 ( h ) 3 2 a h < a 1 h a γ(h) = c 0 + c 1 γ s (h) where: a : c 0 : c 1 : h : γ s γ range parameter nugget variance partial sill variance lag distance standardised semivariance semivariance Return

Definitions Exponential model γ s (h) = 1 exp ( h ) a γ(h) = c 0 + c 1 γ s (h) where: a : c 0 : c 1 : h : γ s γ range parameter nugget variance partial sill variance lag distance standardised semivariance semivariance Return

Definitions Gaussian model γ s (h) = 1 exp ( h2 a 2 γ(h) = c 0 + c 1 γ s (h) ) where: a : c 0 : c 1 : h : γ s γ range parameter nugget variance partial sill variance lag distance standardised semivariance semivariance Return

Definitions Practical range Lag h for which γ(h) = 0.95γ( ), i.e., that distance at which the semivariance is 95% of the sill. Return