University of California, Los Angeles Department of Statistics. Universal kriging
|
|
- Estella Lambert
- 6 years ago
- Views:
Transcription
1 University of California, Los Angeles Department of Statistics Statistics C173/C273 Instructor: Nicolas Christou Universal kriging The Ordinary Kriging (OK) that was discussed earlier is based on the constant mean model given by Z(s) = µ + δ(s) where δ( ) has mean zero and variogram 2γ( ) Many times this is a too simple model to use The mean can be a function of the coordinates X, Y, in some linear, quadratic, or higher form For example the value of Z at location s can be expressed now as Z(s i ) = β 0 + β 1 X i + β 2 Y i + δ(s i ), linear Or Z(s i ) = β 0 + β 1 X i + β 2 Y i + β 3 X 2 i + β 4 X i Y i + β 5 Y 2 i + δ(s i ), quadratic, etc If this is the case then we say that there is a trend of the polynomial type We need to take this into account when we find the kriging weights The predicted value Z(s 0 ) at location s 0 will be again a linear combination of the observed Z(s i ), i = 1,, n values: where Ẑ(s 0 ) = w 1 Z(s 1 ) + w 2 Z(s 2 ) + + w n Z(s n ) = w i Z i w 1 + w w n = 1 Suppose now that a trend of the linear form is present Then the value Ẑ(s 0) can be expressed as or Ẑ(s 0 ) = Ẑ(s 0 ) = w i Z i = w i β 0 + w i β 1 X i + w i β 2 Y i + w i δ(s i ) n n w i Z i = β 0 + β 1 w i X i + β 2 w i Y i + w i δ(s i ) (1) But also the value of Z(s 0 ) can be expressed (based on the linear trend) as Z(s 0 ) = β 0 + β 1 X 0 + β 2 Y 0 + δ(s 0 ) (2) Compare (1) and (2) In order to ensure that we have an unbiased estimator we will need the following conditions: 1
2 and w i X i = X 0 w i Y i = Y 0 w i = 1 As with ordinary kriging, to find the weights when a trend is present we need to minimize the mean square error (MSE) min ( ) 2 Z(s 0 ) w i Z(s i ) subject to the above constraints This minimization will be unconstrained if we incorporate the 3 constraints in the objective function The result is a system of n + 3 equations for the linear trend If the trend is quadratic we will need n + 6 equations, and n + 10 equations for a threedimentional quadratic trend, etc On the next page the system of equations for the linear trend example is presented in matrix form Note: We should try to understand why the trend exists based on the nature of our data, use a simple form of the trend if possible, and avoid extrapolation beyond the available data Once we decided about which trend to use, we subtract this trend from the observed data to obtain the residuals We then use the residuals to compute the sample variogram, fit a model variogram to it, predict the values at the unsampled locations ( kriged the residuals), and finally add the kriged residuals back to the trend 2
3 The Kriging System with linear trend as a function of the x, y coordinates γ(s1 s1) γ(s1 s2) γ(s1 s3) γ(s1 sn) x1 y1 1 γ(s2 s1) γ(s2 s2) γ(s2 s3) γ(s2 sn) x2 y2 1 w1 w2 = γ(sn s1) γ(sn s2) γ(sn s3) γ(sn sn) xn yn 1 wn x1 x2 xn λ1 y1 y2 yn λ2 λ0 γ(s0 s1) γ(s0 s2) γ(s0 sn) x0 y0 1 3
4 Example of universal kriging: Let s return to the simple example with the 7 observed points as shown below We use the exponential semivariogram model with c 0 = 0, c 1 = 10, α = 333, γ(h) = 10(1 e h 333 ) s i x y z(s i ) s s s s s s s s ??? First we calculate the distance matrix as shown below: Distance matrix = s 0 s 1 s 2 s 3 s 4 s 5 s 6 s 7 s s s s s s s s Using the universal kriging equations (assuming that there is a linear trend in our data) the weights, the estimate, and its variance are computed as follows: Universal kriging using the variogram: W = Γ 1 γ =
5 The answer is: W = w 1 w 2 w 3 w 4 w 5 w 6 w 7 λ 1 λ 2 λ 0 = The last 3 elements of the W vector are the Lagrange multipliers, λ 0 = 23, λ 1 = 0066, λ 2 = 0023 We can verify that the sum of the elements 1 through 7 is equal to 1, as it should be The estimate is equal to: 7 ẑ(s 0 ) = w i z(s i ) = 0220(477) (783) = The variance of the estimate is: Or σ 2 (s 0 ) = w i γ(s 0 s i ) + λ 0 + λ 1 x λ k x k0 σ 2 (s 0 ) = 0220(7384) (9823) (65) (137) =
6 Universal kriging using the covariance: The same example is solved using the covariance function: C(h) = 10e h 333 Reminder: The corresponding semivariogram function was γ(h) = 10(1 e h 333 ) W 1 = C 1 1 c 1 = The answer is: W 1 = w 1 w 2 w 3 w 4 w 5 w 6 w 7 λ 1 λ 2 λ = The last 3 elements of the W 1 vector are the Lagrange multipliers, λ 0 = 23, λ 1 = 0067, λ 2 = 0023 We can verify that the sum of the elements 1 through 7 is equal to 1, as it should be The estimate is equal to: 7 ẑ(s 0 ) = w i z(s i ) = 0220(477) (783) = The variance of the estimate is: Or σ 2 (s 0 ) = C(0) w i C(s 0 s i ) + λ 0 + λ 1 x λ k x k0 σ 2 (s 0 ) = 10 [0220(2612) (0176)] (65) (137) =
7 Universal kriging using geor > a <- readtable(" kriging_11txt", header=true) > b <- asgeodata(a) Using the ksline function: > univ_kr <- ksline(b, covmodel="exp", covpars=c(10,333), nugget=0, locations=c(65,137), m0="kt", trend=1) ksline: kriging location: 1 out of 1 Kriging performed using global neighbourhood Warning message: assuming that there is only 1 prediction point in: checklocations(locations) And here is the estimate and its variance: > univ_kr$predict [1] > univ_kr$krigevar [1] Note: The argument trend=1 is required if the argument m0="kt" is used Using the krigeconv function: The locations argument in the krigeconv function must be a data frame You can create it as dd <- dataframe(65,137) Or simply use dataframe(65,137) as shown below: > kc <- krigeconv(b, locations=dd,krige=krigecontrol(typekrige="ok", covmodel="exp", covpars=c(10, 333), nugget=0, trendl="1st", trendd="1st")) kc <- krigeconv(b, locations=dataframe(65,137),krige=krigecontrol(typekrige="ok", covmodel="exp", covpars=c(10, 333), nugget=0, trendl="1st", trendd="1st")) Note: The arguments trendl, trendd specify the trend of the locations to be estimated and the trend of the data They must be the same A complete example on how to use universal kriging in geor and gstat is presented on the next page using the Wolfcamp data set 7
8 #Read the data: site=" a <- readtable(file=site, header=true) #Load the geor and gstat packages: library(geor) library(gstat) #Convert the data into a geodata object (for geor): b <- asgeodata(a) #Create the grid for spatial predictions: xrange <- asinteger(range(a[,1])) yrange <- asinteger(range(a[,2])) grd <- expandgrid(x=seq(from=xrange[1], to=xrange[2], by=2), y=seq(from=yrange[1], to=yrange[2], by=2)) #Perform universal kriging (geor using the function ksline): krig_trend <- ksline(b, covmodel="sph", covpars=c(30000, 60), nugget=10000, m0="kt", trend=1, locations=grd) #Perform universal kriging (geor using the function krigeconv): kc <- krigeconv(b, locations=grd,krige=krigecontrol(typekrige="ok", covmodel="sph", covpars=c(30000, 60), nugget=10000, trendl="1st", trendd="1st")) #Perform universal kriging (gstat - krige): pr_univk <- krige(id="level", level~x+y, locations=~x+y, model=vgm(30000, "Sph", 60, 10000), data=a, newdata=grd) #Compare the three methods: pred <- cbind(pr_univk$levelpred, krig_trend$pred, kc$predict) predvar <- cbind(pr_univk$levelvar, krig_trend$krigevar, kc$krigevar) #Display the first 5 rows: > pred[1:5,] [,1] [,2] [,3] [1,] [2,] [3,] [4,] [5,]
9 > predvar[1:5,] [,1] [,2] [,3] [1,] [2,] [3,] [4,] [5,] #Load the lattice package: library(lattice) #Construct a raster map using the kriged values: levelplot(pr_univk$levelpred~x+y, pr_univk, aspect ="iso", main="universal kriging predictions") #Construct a raster map using the variances of the kriged values: levelplot(pr_univk$levelvar~x+y, pr_univk, aspect ="iso", main="universal kriging variances") Here are the plots: Universal kriging predictions y x Universal kriging variances y x 9
10 Using the imageorig function: #Collapse the predicted values into a matrix: qqq <- matrix(pr_univk$levelpred, length(seq(from=xrange[1], to=xrange[2], by=2)), length(seq(from=yrange[1], to=yrange[2], by=2))) #Construct the raster map of the predicted values: imageorig(seq(from=xrange[1], to=xrange[2], by=2), seq(from=yrange[1],to=yrange[2], by=2), qqq, xlab="west to East",ylab="South to North", main="raster map of the predicted values") #Collapse the variances into a matrix: qqq1 <- matrix(pr_univk$levelvar, length(seq(from=xrange[1], to=xrange[2], by=2)), length(seq(from=yrange[1], to=yrange[2], by=2))) #Construct the raster map of the variances: imageorig(seq(from=xrange[1], to=xrange[2], by=2), seq(from=yrange[1],to=yrange[2], by=2), qqq1, xlab="west to East",ylab="South to North", main="raster map of the variances") Here are the plots: Raster map of the predicted values South to North West to East Raster map of the variances South to North West to East 10
University of California, Los Angeles Department of Statistics. Geostatistical data: variogram and kriging
University of California, Los Angeles Department of Statistics Statistics 403 Instructor: Nicolas Christou Geostatistical data: variogram and kriging Let Z(s) and Z(s + h) two random variables at locations
More informationPoint patterns. The average number of events (the intensity, (s)) is homogeneous
Point patterns Consider the point process {Z(s) :s 2 D R 2 }.Arealization of this process consists of a pattern (arrangement) of points in D. (D is a random set.) These points are called the events of
More informationUniversity of California, Los Angeles Department of Statistics. Data with trend
University of California, Los Angeles Department of Statistics Statistics C173/C273 Instructor: Nicolas Christou Data with trend As we discussed earlier, if the intrinsic stationarity assumption holds,
More informationUniversity of California, Los Angeles Department of Statistics
University of California, Los Angeles Department of Statistics Statistics C173/C273 Instructor: Nicolas Christou Working with gstat - An example We will use the data from the Meuse (Dutch Maas) river.
More informationUniversity of California, Los Angeles Department of Statistics. Effect of variogram parameters on kriging weights
University of California, Los Angeles Department of Statistics Statistics C173/C273 Instructor: Nicolas Christou Effect of variogram parameters on kriging weights We will explore in this document how the
More informationREN R 690 Lab Geostatistics Lab
REN R 690 Lab Geostatistics Lab The objective of this lab is to try out some basic geostatistical tools with R. Geostatistics is used primarily in the resource and environmental areas for estimation, uncertainty
More informationLecture 9: Introduction to Kriging
Lecture 9: Introduction to Kriging Math 586 Beginning remarks Kriging is a commonly used method of interpolation (prediction) for spatial data. The data are a set of observations of some variable(s) of
More informationKriging Luc Anselin, All Rights Reserved
Kriging Luc Anselin Spatial Analysis Laboratory Dept. Agricultural and Consumer Economics University of Illinois, Urbana-Champaign http://sal.agecon.uiuc.edu Outline Principles Kriging Models Spatial Interpolation
More informationUniversity of California, Los Angeles Department of Statistics. Introduction
University of California, Los Angeles Department of Statistics Statistics C173/C273 Instructor: Nicolas Christou Introduction What is geostatistics? Geostatistics is concerned with estimation and prediction
More informationData Break 8: Kriging the Meuse RiverBIOS 737 Spring 2004 p.1/27
Data Break 8: Kriging the Meuse River BIOS 737 Spring 2004 Data Break 8: Kriging the Meuse RiverBIOS 737 Spring 2004 p.1/27 Meuse River: Reminder library(gstat) Data included in gstat library. data(meuse)
More informationWhat s for today. All about Variogram Nugget effect. Mikyoung Jun (Texas A&M) stat647 lecture 4 September 6, / 17
What s for today All about Variogram Nugget effect Mikyoung Jun (Texas A&M) stat647 lecture 4 September 6, 2012 1 / 17 What is the variogram? Let us consider a stationary (or isotropic) random field Z
More informationPractical 12: Geostatistics
Practical 12: Geostatistics This practical will introduce basic tools for geostatistics in R. You may need first to install and load a few packages. The packages sp and lattice contain useful function
More informationPoint-Referenced Data Models
Point-Referenced Data Models Jamie Monogan University of Georgia Spring 2013 Jamie Monogan (UGA) Point-Referenced Data Models Spring 2013 1 / 19 Objectives By the end of these meetings, participants should
More informationIntroduction. Semivariogram Cloud
Introduction Data: set of n attribute measurements {z(s i ), i = 1,, n}, available at n sample locations {s i, i = 1,, n} Objectives: Slide 1 quantify spatial auto-correlation, or attribute dissimilarity
More informationTypes of Spatial Data
Spatial Data Types of Spatial Data Point pattern Point referenced geostatistical Block referenced Raster / lattice / grid Vector / polygon Point Pattern Data Interested in the location of points, not their
More informationPRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH
PRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH SURESH TRIPATHI Geostatistical Society of India Assumptions and Geostatistical Variogram
More informationExploring the World of Ordinary Kriging. Dennis J. J. Walvoort. Wageningen University & Research Center Wageningen, The Netherlands
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
More informationPractical 12: Geostatistics
Practical 12: Geostatistics This practical will introduce basic tools for geostatistics in R. You may need first to install and load a few packages. The packages sp and lattice contain useful function
More informationGRAD6/8104; INES 8090 Spatial Statistic Spring 2017
Lab #4 Semivariogram and Kriging (Due Date: 04/04/2017) PURPOSES 1. Learn to conduct semivariogram analysis and kriging for geostatistical data 2. Learn to use a spatial statistical library, gstat, in
More information11. Kriging. ACE 492 SA - Spatial Analysis Fall 2003
11. Kriging ACE 492 SA - Spatial Analysis Fall 2003 c 2003 by Luc Anselin, All Rights Reserved 1 Objectives The goal of this lab is to further familiarize yourself with ESRI s Geostatistical Analyst, extending
More informationInvestigation of Monthly Pan Evaporation in Turkey with Geostatistical Technique
Investigation of Monthly Pan Evaporation in Turkey with Geostatistical Technique Hatice Çitakoğlu 1, Murat Çobaner 1, Tefaruk Haktanir 1, 1 Department of Civil Engineering, Erciyes University, Kayseri,
More informationIntroduction to Geostatistics
Introduction to Geostatistics Abhi Datta 1, Sudipto Banerjee 2 and Andrew O. Finley 3 July 31, 2017 1 Department of Biostatistics, Bloomberg School of Public Health, Johns Hopkins University, Baltimore,
More informationSPATIAL ELECTRICAL LOADS MODELING USING THE GEOSTATISTICAL METHODS
19 th International CODATA Conference THE INFORMATION SOCIETY: NEW HORIZONS FOR SCIENCE Berlin, Germany 7-1 November 24 SPATIAL ELECTRICAL LOADS MODELING USING THE GEOSTATISTICAL METHODS Barbara Namysłowska-Wilczyńska
More informationChapter 4 - Fundamentals of spatial processes Lecture notes
TK4150 - Intro 1 Chapter 4 - Fundamentals of spatial processes Lecture notes Odd Kolbjørnsen and Geir Storvik January 30, 2017 STK4150 - Intro 2 Spatial processes Typically correlation between nearby sites
More informationBasics of Point-Referenced Data Models
Basics of Point-Referenced Data Models Basic tool is a spatial process, {Y (s), s D}, where D R r Chapter 2: Basics of Point-Referenced Data Models p. 1/45 Basics of Point-Referenced Data Models Basic
More informationBasics in Geostatistics 2 Geostatistical interpolation/estimation: Kriging methods. Hans Wackernagel. MINES ParisTech.
Basics in Geostatistics 2 Geostatistical interpolation/estimation: Kriging methods Hans Wackernagel MINES ParisTech NERSC April 2013 http://hans.wackernagel.free.fr Basic concepts Geostatistics Hans Wackernagel
More informationChapter 1. Summer School GEOSTAT 2014, Spatio-Temporal Geostatistics,
Chapter 1 Summer School GEOSTAT 2014, Geostatistics, 2014-06-19 sum- http://ifgi.de/graeler Institute for Geoinformatics University of Muenster 1.1 Spatial Data From a purely statistical perspective, spatial
More informationGeog 210C Spring 2011 Lab 6. Geostatistics in ArcMap
Geog 210C Spring 2011 Lab 6. Geostatistics in ArcMap Overview In this lab you will think critically about the functionality of spatial interpolation, improve your kriging skills, and learn how to use several
More information2.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics. References - geostatistics. References geostatistics (cntd.
.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics Spline interpolation was originally developed or image processing. In GIS, it is mainly used in visualization o spatial
More informationSpatial-Temporal Modeling of Active Layer Thickness
Spatial-Temporal Modeling of Active Layer Thickness Qian Chen Advisor : Dr. Tatiyana Apanasovich Department of Statistics The George Washington University Abstract The objective of this study is to provide
More informationGridding of precipitation and air temperature observations in Belgium. Michel Journée Royal Meteorological Institute of Belgium (RMI)
Gridding of precipitation and air temperature observations in Belgium Michel Journée Royal Meteorological Institute of Belgium (RMI) Gridding of meteorological data A variety of hydrologic, ecological,
More informationHandbook of Spatial Statistics Chapter 2: Continuous Parameter Stochastic Process Theory by Gneiting and Guttorp
Handbook of Spatial Statistics Chapter 2: Continuous Parameter Stochastic Process Theory by Gneiting and Guttorp Marcela Alfaro Córdoba August 25, 2016 NCSU Department of Statistics Continuous Parameter
More informationENGRG Introduction to GIS
ENGRG 59910 Introduction to GIS Michael Piasecki October 13, 2017 Lecture 06: Spatial Analysis Outline Today Concepts What is spatial interpolation Why is necessary Sample of interpolation (size and pattern)
More informationBeta-Binomial Kriging: An Improved Model for Spatial Rates
Available online at www.sciencedirect.com ScienceDirect Procedia Environmental Sciences 27 (2015 ) 30 37 Spatial Statistics 2015: Emerging Patterns - Part 2 Beta-Binomial Kriging: An Improved Model for
More informationSpatial Data Analysis in Archaeology Anthropology 589b. Kriging Artifact Density Surfaces in ArcGIS
Spatial Data Analysis in Archaeology Anthropology 589b Fraser D. Neiman University of Virginia 2.19.07 Spring 2007 Kriging Artifact Density Surfaces in ArcGIS 1. The ingredients. -A data file -- in.dbf
More informationNumerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error. 2- Fixed point
Numerical Analysis Solution of Algebraic Equation (non-linear equation) 1- Trial and Error In this method we assume initial value of x, and substitute in the equation. Then modify x and continue till we
More informationUmeå University Sara Sjöstedt-de Luna Time series analysis and spatial statistics
Umeå University 01-05-5 Sara Sjöstedt-de Luna Time series analysis and spatial statistics Laboration in ArcGIS Geostatistical Analyst These exercises are aiming at helping you understand ArcGIS Geostatistical
More information7 Geostatistics. Figure 7.1 Focus of geostatistics
7 Geostatistics 7.1 Introduction Geostatistics is the part of statistics that is concerned with geo-referenced data, i.e. data that are linked to spatial coordinates. To describe the spatial variation
More informationChapter 4 - Fundamentals of spatial processes Lecture notes
Chapter 4 - Fundamentals of spatial processes Lecture notes Geir Storvik January 21, 2013 STK4150 - Intro 2 Spatial processes Typically correlation between nearby sites Mostly positive correlation Negative
More informationgeorglm : a package for generalised linear spatial models introductory session
georglm : a package for generalised linear spatial models introductory session Ole F. Christensen & Paulo J. Ribeiro Jr. Last update: October 18, 2017 The objective of this page is to introduce the reader
More informationPackage STMedianPolish
Type Package Title Spatio-Temporal Median Polish Version 0.2 Date 2017-03-07 Package STMedianPolish March 8, 2017 Author William Martínez [aut, cre], Melo Carlos [aut],melo Oscar [aut]. Maintainer William
More informationGeostatistics in Hydrology: Kriging interpolation
Chapter Geostatistics in Hydrology: Kriging interpolation Hydrologic properties, such as rainfall, aquifer characteristics (porosity, hydraulic conductivity, transmissivity, storage coefficient, etc.),
More informationSpatial Data Mining. Regression and Classification Techniques
Spatial Data Mining Regression and Classification Techniques 1 Spatial Regression and Classisfication Discrete class labels (left) vs. continues quantities (right) measured at locations (2D for geographic
More information5. Geostatistics JEAN-MICHEL FLOCH INSEE. Abstract
5. Geostatistics JEAN-MICHEL FLOCH INSEE 5.1 Random functions 114 5.1.1 Definitions.............................................. 114 5.1.2 Stationarity............................................. 115
More informationREML Estimation and Linear Mixed Models 4. Geostatistics and linear mixed models for spatial data
REML Estimation and Linear Mixed Models 4. Geostatistics and linear mixed models for spatial data Sue Welham Rothamsted Research Harpenden UK AL5 2JQ December 1, 2008 1 We will start by reviewing the principles
More informationSpatial Linear Geostatistical Models
slm.geo{slm} Spatial Linear Geostatistical Models Description This function fits spatial linear models using geostatistical models. It can estimate fixed effects in the linear model, make predictions for
More informationSpatial analysis is the quantitative study of phenomena that are located in space.
c HYON-JUNG KIM, 2016 1 Introduction Spatial analysis is the quantitative study of phenomena that are located in space. Spatial data analysis usually refers to an analysis of the observations in which
More informationThe ProbForecastGOP Package
The ProbForecastGOP Package April 24, 2006 Title Probabilistic Weather Field Forecast using the GOP method Version 1.3 Author Yulia Gel, Adrian E. Raftery, Tilmann Gneiting, Veronica J. Berrocal Description
More informationChapter 14 Stein-Rule Estimation
Chapter 14 Stein-Rule Estimation The ordinary least squares estimation of regression coefficients in linear regression model provides the estimators having minimum variance in the class of linear and unbiased
More informationIntroduction to applied geostatistics. Short version. Overheads
Introduction to applied geostatistics Short version Overheads Department of Earth Systems Analysis International Institute for Geo-information Science & Earth Observation (ITC)
More informationAnomaly Density Estimation from Strip Transect Data: Pueblo of Isleta Example
Anomaly Density Estimation from Strip Transect Data: Pueblo of Isleta Example Sean A. McKenna, Sandia National Laboratories Brent Pulsipher, Pacific Northwest National Laboratory May 5 Distribution Statement
More informationAn Introduction to Spatial Statistics. Chunfeng Huang Department of Statistics, Indiana University
An Introduction to Spatial Statistics Chunfeng Huang Department of Statistics, Indiana University Microwave Sounding Unit (MSU) Anomalies (Monthly): 1979-2006. Iron Ore (Cressie, 1986) Raw percent data
More information11/8/2018. Spatial Interpolation & Geostatistics. Kriging Step 1
(Z i Z j ) 2 / 2 (Z i Zj) 2 / 2 Semivariance y 11/8/2018 Spatial Interpolation & Geostatistics Kriging Step 1 Describe spatial variation with Semivariogram Lag Distance between pairs of points Lag Mean
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry
More informationBayesian Transgaussian Kriging
1 Bayesian Transgaussian Kriging Hannes Müller Institut für Statistik University of Klagenfurt 9020 Austria Keywords: Kriging, Bayesian statistics AMS: 62H11,60G90 Abstract In geostatistics a widely used
More informationCopyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
Chapter 15. SPATIAL INTERPOLATION 15.1 Elements of Spatial Interpolation 15.1.1 Control Points 15.1.2 Type of Spatial Interpolation 15.2 Global Methods 15.2.1 Trend Surface Models Box 15.1 A Worked Example
More informationROeS Seminar, November
IASC Introduction: Spatial Interpolation Estimation at a certain location Geostatistische Modelle für Fließgewässer e.g. Air pollutant concentrations were measured at different locations. What is the concentration
More informationModels for spatial data (cont d) Types of spatial data. Types of spatial data (cont d) Hierarchical models for spatial data
Hierarchical models for spatial data Based on the book by Banerjee, Carlin and Gelfand Hierarchical Modeling and Analysis for Spatial Data, 2004. We focus on Chapters 1, 2 and 5. Geo-referenced data arise
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Department of Forestry & Department of Geography, Michigan State University, Lansing Michigan, U.S.A. 2 Biostatistics,
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 1 y Kriging Step 1 Describe spatial variation with Semivariogram (Z i Z j ) 2 / 2 Point cloud Map 3
More informationMUCK - A ovel Approach to Co-Kriging
Clark, W.V. Harper and K.L. Basinger, MUCK -- A novel approach to co-kriging in Proc. DOE/AECL 87 Conference on Geostatistical, Sensitivity and Uncertainty Methods for Ground-Water Flow and Radionuclide
More informationPracticum : Spatial Regression
: Alexandra M. Schmidt Instituto de Matemática UFRJ - www.dme.ufrj.br/ alex 2014 Búzios, RJ, www.dme.ufrj.br Exploratory (Spatial) Data Analysis 1. Non-spatial summaries Numerical summaries: Mean, median,
More informationRecent Developments in Biostatistics: Space-Time Models. Tuesday: Geostatistics
Recent Developments in Biostatistics: Space-Time Models Tuesday: Geostatistics Sonja Greven Department of Statistics Ludwig-Maximilians-Universität München (with special thanks to Michael Höhle for some
More informationLecture 2. The Simple Linear Regression Model: Matrix Approach
Lecture 2 The Simple Linear Regression Model: Matrix Approach Matrix algebra Matrix representation of simple linear regression model 1 Vectors and Matrices Where it is necessary to consider a distribution
More informationGstat: multivariable geostatistics for S
DSC 2003 Working Papers (Draft Versions) http://www.ci.tuwien.ac.at/conferences/dsc-2003/ Gstat: multivariable geostatistics for S Edzer J. Pebesma Dept. of Physical Geography, Utrecht University, P.O.
More informationSample questions for Fundamentals of Machine Learning 2018
Sample questions for Fundamentals of Machine Learning 2018 Teacher: Mohammad Emtiyaz Khan A few important informations: In the final exam, no electronic devices are allowed except a calculator. Make sure
More informationTime-lapse filtering and improved repeatability with automatic factorial co-kriging. Thierry Coléou CGG Reservoir Services Massy
Time-lapse filtering and improved repeatability with automatic factorial co-kriging. Thierry Coléou CGG Reservoir Services Massy 1 Outline Introduction Variogram and Autocorrelation Factorial Kriging Factorial
More informationIndex. Geostatistics for Environmental Scientists, 2nd Edition R. Webster and M. A. Oliver 2007 John Wiley & Sons, Ltd. ISBN:
Index Akaike information criterion (AIC) 105, 290 analysis of variance 35, 44, 127 132 angular transformation 22 anisotropy 59, 99 affine or geometric 59, 100 101 anisotropy ratio 101 exploring and displaying
More informationRegression with correlation for the Sales Data
Regression with correlation for the Sales Data Scatter with Loess Curve Time Series Plot Sales 30 35 40 45 Sales 30 35 40 45 0 10 20 30 40 50 Week 0 10 20 30 40 50 Week Sales Data What is our goal with
More informationI 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
Spatial analysis Huge topic! Key references Diggle (point patterns); Cressie (everything); Diggle and Ribeiro (geostatistics); Dormann et al (GLMMs for species presence/abundance); Haining; (Pinheiro and
More informationGeostatistics: Kriging
Geostatistics: Kriging 8.10.2015 Konetekniikka 1, Otakaari 4, 150 10-12 Rangsima Sunila, D.Sc. Background What is Geostatitics Concepts Variogram: experimental, theoretical Anisotropy, Isotropy Lag, Sill,
More informationGstat: Multivariable Geostatistics for S
New URL: http://www.r-project.org/conferences/dsc-2003/ Proceedings of the 3rd International Workshop on Distributed Statistical Computing (DSC 2003) March 20 22, Vienna, Austria ISSN 1609-395X Kurt Hornik,
More informationInterpolation and 3D Visualization of Geodata
Marek KULCZYCKI and Marcin LIGAS, Poland Key words: interpolation, kriging, real estate market analysis, property price index ABSRAC Data derived from property markets have spatial character, no doubt
More informationApplication and evaluation of universal kriging for optimal contouring of groundwater levels
Application and evaluation of universal kriging for optimal contouring of groundwater levels B V N P Kambhammettu 1,, Praveena Allena 2, and JamesPKing 1, 1 Civil Engineering Department, New Mexico State
More informationMETHODOLOGY WHICH APPLIES GEOSTATISTICS TECHNIQUES TO THE TOPOGRAPHICAL SURVEY
International Journal of Computer Science and Applications, 2008, Vol. 5, No. 3a, pp 67-79 Technomathematics Research Foundation METHODOLOGY WHICH APPLIES GEOSTATISTICS TECHNIQUES TO THE TOPOGRAPHICAL
More informationPackage ProbForecastGOP
Type Package Package ProbForecastGOP February 19, 2015 Title Probabilistic weather forecast using the GOP method Version 1.3.2 Date 2010-05-31 Author Veronica J. Berrocal , Yulia
More informationSimple Linear Regression
Simple Linear Regression MATH 282A Introduction to Computational Statistics University of California, San Diego Instructor: Ery Arias-Castro http://math.ucsd.edu/ eariasca/math282a.html MATH 282A University
More informationSpatio-temporal analysis and interpolation of PM10 measurements in Europe for 2009
Spatio-temporal analysis and interpolation of PM10 measurements in Europe for 2009 ETC/ACM Technical Paper 2012/8 March 2013 revised version Benedikt Gräler, Mirjam Rehr, Lydia Gerharz, Edzer Pebesma The
More informationWhat s for today. Introduction to Space-time models. c Mikyoung Jun (Texas A&M) Stat647 Lecture 14 October 16, / 19
What s for today Introduction to Space-time models c Mikyoung Jun (Texas A&M) Stat647 Lecture 14 October 16, 2012 1 / 19 Space-time Data So far we looked at the data that vary over space Now we add another
More informationChapter 4 - Spatial processes R packages and software Lecture notes
TK4150 - Intro 1 Chapter 4 - Spatial processes R packages and software Lecture notes Odd Kolbjørnsen and Geir Storvik February 13, 2017 STK4150 - Intro 2 Last time General set up for Kriging type problems
More informationLOGNORMAL ORDINARY KRIGING METAMODEL IN SIMULATION OPTIMIZATION
LOGNORMAL ORDINARY KRIGING METAMODEL IN SIMULATION OPTIMIZATION Muzaffer Balaban 1 and Berna Dengiz 2 1 Turkish Statistical Institute, Ankara, Turkey 2 Department of Industrial Engineering, Başkent University,
More informationPropagation of Errors in Spatial Analysis
Stephen F. Austin State University SFA ScholarWorks Faculty Presentations Spatial Science 2001 Propagation of Errors in Spatial Analysis Peter P. Siska I-Kuai Hung Arthur Temple College of Forestry and
More informationAssessing the covariance function in geostatistics
Statistics & Probability Letters 52 (2001) 199 206 Assessing the covariance function in geostatistics Ana F. Militino, M. Dolores Ugarte Departamento de Estadstica e Investigacion Operativa, Universidad
More informationRoger S. Bivand Edzer J. Pebesma Virgilio Gömez-Rubio. Applied Spatial Data Analysis with R. 4:1 Springer
Roger S. Bivand Edzer J. Pebesma Virgilio Gömez-Rubio Applied Spatial Data Analysis with R 4:1 Springer Contents Preface VII 1 Hello World: Introducing Spatial Data 1 1.1 Applied Spatial Data Analysis
More information(x + 1)(x 2) = 4. x
dvanced Integration Techniques: Partial Fractions The method of partial fractions can occasionally make it possible to find the integral of a quotient of rational functions. Partial fractions gives us
More informationOptimizing Sampling Schemes for Mapping and Dredging Polluted Sediment Layers
This file was created by scanning the printed publication. Errors identified by the software have been corrected; however, some errors may remain. Optimizing Sampling Schemes for Mapping and Dredging Polluted
More informationReport on Kriging in Interpolation
Tabor Reedy ENVS421 3/12/15 Report on Kriging in Interpolation In this project I explored use of the geostatistical analyst extension and toolbar in the process of creating an interpolated surface through
More informationAnalysing Spatial Data in R: Worked example: geostatistics
Analysing Spatial Data in R: Worked example: geostatistics Roger Bivand Department of Economics Norwegian School of Economics and Business Administration Bergen, Norway 17 April 2007 Worked example: geostatistics
More informationA Geostatistical Approach to Predict the Average Annual Rainfall of Bangladesh
Journal of Data Science 14(2016), 149-166 A Geostatistical Approach to Predict the Average Annual Rainfall of Bangladesh Mohammad Samsul Alam 1 and Syed Shahadat Hossain 1 1 Institute of Statistical Research
More informationCS1800: Sequences & Sums. Professor Kevin Gold
CS1800: Sequences & Sums Professor Kevin Gold Moving Toward Analysis of Algorithms Today s tools help in the analysis of algorithms. We ll cover tools for deciding what equation best fits a sequence of
More informationA prospective and retrospective spatial sampling scheme to characterize geochemicals in a mine tailings area
A prospective and retrospective spatial sampling scheme to characterize geochemicals in a mine tailings area Pravesh Debba CSIR Built Environment Logistics and Quantitative Methods (LQM) South Africa Presented
More informationWarm-Up. Simplify the following terms:
Warm-Up Simplify the following terms: 81 40 20 i 3 i 16 i 82 TEST Our Ch. 9 Test will be on 5/29/14 Complex Number Operations Learning Targets Adding Complex Numbers Multiplying Complex Numbers Rules for
More informationExperimental variogram of the residuals in the universal kriging (UK) model
Experimental variogram of the residuals in the universal kriging (UK) model Nicolas DESASSIS Didier RENARD Technical report R141210NDES MINES ParisTech Centre de Géosciences Equipe Géostatistique 35, rue
More informationAlgebra II Chapter 5: Polynomials and Polynomial Functions Part 1
Algebra II Chapter 5: Polynomials and Polynomial Functions Part 1 Chapter 5 Lesson 1 Use Properties of Exponents Vocabulary Learn these! Love these! Know these! 1 Example 1: Evaluate Numerical Expressions
More informationTABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1. Chapter Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9
TABLE OF CONTENTS INTRODUCTION, APPROXIMATION & ERRORS 1 Chapter 01.01 Introduction to numerical methods 1 Multiple-choice test 7 Problem set 9 Chapter 01.02 Measuring errors 11 True error 11 Relative
More informationStatistics Univariate Linear Models Gary W. Oehlert School of Statistics 313B Ford Hall
Statistics 5401 14. Univariate Linear Models Gary W. Oehlert School of Statistics 313B ord Hall 612-625-1557 gary@stat.umn.edu Linear models relate a target or response or dependent variable to known predictor
More informationGeographically Weighted Regression and Kriging: Alternative Approaches to Interpolation A Stewart Fotheringham
Geographically Weighted Regression and Kriging: Alternative Approaches to Interpolation A Stewart Fotheringham National Centre for Geocomputation National University of Ireland, Maynooth http://ncg.nuim.ie
More informationValidation of an offshore wind atlas using satellite data available at coastal regions of Portugal
Validation of an offshore wind atlas using satellite data available at coastal regions of Portugal R. Marujo, P. Costa, M. Fernandes, T. Simões, A. Estanqueiro European Seminar OWEMES 2012 SESSION 1: Wind
More informationOn block bootstrapping areal data Introduction
On block bootstrapping areal data Nicholas Nagle Department of Geography University of Colorado UCB 260 Boulder, CO 80309-0260 Telephone: 303-492-4794 Email: nicholas.nagle@colorado.edu Introduction Inference
More informationUnit 1 Notes. Polynomials
Unit 1 Notes 1 Day Number Date Topic Problem Set 1 Wed. Sept. Operations with Signed Numbers and Order of Operations Pre-Unit Review PS (1 ) Thurs. Sept. Working with Exponents Start Exponent Laws P.S.
More information