Skip to main content

Questions tagged [kriging]

A group of geostatistical techniques to interpolate the value of a random field at an unobserved location from observations of its value at nearby locations.

Filter by
Sorted by
Tagged with
44 votes
9 answers
32k views

Open source methods for kriging?

I have a point dataset which I'd like to Krige, ideally using an open-source software package. If possible, I'd also like to choose the semi-variogram model during the process to improve the ...
scw's user avatar
  • 16.4k
27 votes
2 answers
44k views

Choosing IDW vs Kriging Interpolation for DEM creation?

I'm trying to create a DEM using point data spaced very regularly about 10m apart. The area I am interpolating is an educational facility with many flat parking lots and football fields but still has ...
Calavin's user avatar
  • 869
19 votes
2 answers
13k views

Minimum number of samples for kriging interpolation

I get some data with number of specimens with a request to interpolate it using kriging method. After some investigation it appeared that kriging results (performed in ArcGIS Geostatistical Analyst ...
Marcin's user avatar
  • 1,821
15 votes
3 answers
20k views

Converting point data into gridded dataframe for histogram analysis using R?

I am very new at using GIS data and only modestly experienced with R. I've been reading about how to analyze spatial data using the spatial-analyst.net PDF book, so I'm not completely lost, but I ...
user1080253's user avatar
11 votes
2 answers
4k views

R gstat krige() - Covariance matrix singular at location [5.88,47.4,0]: skipping

When I want to perform the kriging it only works sometimes, depending on what values I use in my datatable. As a result of the krige function I get for var1.pred: NA NA NA ... and var1.var: NA NA NA .....
obrob's user avatar
  • 887
11 votes
1 answer
984 views

Kriging on GRASS+R - cannot allocate vector of size 146.5 Mb

I'm trying to do a groundwater elevation map with an Ordinary Kriging on R, but always resulting on the lack of memory cannot allocate vector of size 146.5 Mb allocation (windows xp 32bit). ...
Ulmi's user avatar
  • 501
10 votes
1 answer
873 views

How to grid unevenly sampled categorical data?

I am looking for a method for gridding categorical data. I have extracted from nautical charts and field sheets a collection of points specifying the nature of the surface of the bottom of the ocean. ...
Doug Hrynyk's user avatar
10 votes
3 answers
6k views

Kriging with barriers in ArcGIS Desktop?

I am trying to do interpolation in a coastal area with data points in several rivers emptying into a bay and on opposite sides of peninsulas. I have created a coastline layer to use it as a barrier ...
elayne927's user avatar
  • 171
9 votes
3 answers
2k views

Spatio-temporal block kriging with R package gstat?

In the documentation of sp I saw that the coords argument of the SpatialPolygons function only accepts 2 columns. Is it not possible to create 3-dimensional polygons? I am doing spatio-temporal ...
Nima's user avatar
  • 205
9 votes
1 answer
4k views

Kriging example with sf object

I am new to geo spacial analysis. I would like to learn kriging from a simple example. Say, we have four points with some z value. library(sf) library(gstat) library(tidyverse) dt<-tibble::...
Jacek Kotowski's user avatar
8 votes
1 answer
11k views

Plotting map resulted from kriging in R

I want to plot a map weighted by prices. I got a data frame of latitudes, longitudes and prices: > str(data) 'data.frame': 1253 obs. of 3 variables: $ lon : num 11.2 11.1 11 11 11.3 ... $ ...
Benedict Ernst's user avatar
8 votes
3 answers
2k views

Backtransformation of kriging predictions and variances

I'm not sure how to back-transform log-normal kriged results. This example using the "meuse" data shows how to make variogram and use it to get kriging predictions (and variances) using the popular '...
ToNoY's user avatar
  • 573
8 votes
2 answers
2k views

Interpolating Maps - Statistical Learning Techniques vs Spatial Statistics Techniques

I have recently begun experimenting with the interpolation of variables over a whole country. I have good sized training set where I have several dependent variables (continuous and categorical) ...
JEquihua's user avatar
  • 1,087
8 votes
1 answer
1k views

Hex grid value interpolation in POSTGIS

I have a table of tessellated hexagons that I created as a hex binning grid to cover an area. I've then overlaid a point set to give some of the hexagons a value, see diagram:- The shapes with value ...
user1331131's user avatar
8 votes
1 answer
495 views

Why is SAGA kriging slower in QGIS than in SAGA itself?

I'm using QGIS for a while. In my opinion the SAGA kriging tool (processing toolbox) is very slow. Is that an QGIS issue? I tried to interpolate a point layer with 300.000 points (I know it is big) ...
Pimpel's user avatar
  • 2,353
8 votes
1 answer
390 views

Interpolating points observed with different time duration

I am doing spatial interpolation with a group of points. Some points have complete observation of events: for a period of 24 hours we have observed N events in each point. These points have N events....
Gago-Silva's user avatar
  • 2,546
8 votes
1 answer
370 views

How does Empirical Bayesian Kriging of ArcGIS Pro handle sampling bias?

In ArcGIS Pro, simple kriging approaches allow you to de-cluster your data before applying a normal-score transformation. According to the Esri guideline, normal-score transformations are only ...
jslth's user avatar
  • 101
7 votes
2 answers
6k views

Strange spatial interpolation results from ordinary kriging

I'm using Ordinary Kriging to interpolate a set of spatial data. I'm working on different sets which represent different days of measurement. I get good results, apart from one specific set. Will post ...
natario's user avatar
  • 217
7 votes
2 answers
5k views

How to calculate kriging weights?

In Esri's documentation on how to make a prediction in Kriging, I couldn't figure out the linkage between the semivariograms and the prediction process from reading the description therein: After ...
tinlyx's user avatar
  • 11.1k
7 votes
2 answers
22k views

Performing Kriging in QGIS

Does anyone know how to perform kriging in QGIS? and not via SAGA but actually in the program itself.
C. Golani's user avatar
  • 181
7 votes
3 answers
3k views

Kriging with barriers in R?

I wonder if there is a way to interpolate xyz points in a raster within certain bounds using R? In ArcGIS this is called interpolation with barriers. I wonder if there is something similar in R? I ...
Camilo Mora's user avatar
7 votes
1 answer
417 views

Kriging to calculate forest biomass?

I am trying to create biomass map using Kriging from ground sample point. I have about 150 and 53 point plot for dense and sparse strata. The data is not normal so I have transformed it. After ...
Julia's user avatar
  • 71
7 votes
1 answer
2k views

Are special interpolation techniques required for weather data?

I've downloaded a series of weather observations from the local weather agency - it comprises around 1,000 points spread non-uniformly across Australia, with a greater density of points located near ...
Stephen Lead's user avatar
  • 21.1k
6 votes
2 answers
2k views

Kriging method to use for rainfall?

I am working on an rainfall analysis and I have to interpolate some datas with kriging method. Is co-kriging suitable for the rainfall analysis?
rblerk's user avatar
  • 69
6 votes
2 answers
5k views

Three-dimensional interpolation of temperature profiles using gstat in R?

I'm working on a problem in which I'd like to interpolate a set of oceanographic temperature profiles over a three-dimensional matrix of unsampled x, y, & z locations. The data to interpolate is ...
user30977's user avatar
6 votes
1 answer
5k views

Performing Ordinary Kriging in SAGA via QGIS SEXTANTE Plug-in?

I'm trying to learn how to apply the kriging analysis in QGIS using the SEXTANTE plug-in and SAGA tools. It appears that I have my SEXTANTE and SAGA folders configured properly according to the user ...
Patrick's user avatar
  • 271
6 votes
1 answer
2k views

Dealing with data which are not normally distributed when kriging?

I'm a novice with advanced spatial interpolation, and I'm trying to use the Geostatistical Analyst to interpolate a surface based on borehole data, using Ordinary Kriging. The help file explains how ...
Stephen Lead's user avatar
  • 21.1k
6 votes
1 answer
2k views

How do I implement 'weighted' kriging in R? [closed]

I'm trying to interpolate the values of an attribute (called "disease_rate") that i've sampled at a set of spatial points and i'd like to add 'weights' to the resulting kriging estimates. Basically, ...
Alberto's user avatar
  • 101
6 votes
1 answer
808 views

Cokriging: zero distance cross-semivariance in undersampled cases [gstat error]

Trying to cokrige two variables that are not perfectly colocated (one has sparser measurements than the other), I faced an issue that I'll illustrate with the following MRE. With the meuse dataset, ...
Yollanda Beetroot's user avatar
5 votes
2 answers
2k views

Method for estimating accuracy in spatial interpolation?

When using a spatial interpolation method (any method) without having a way to validate it, besides sampling some points out, is there a way to count as well with the mean distance between ...
Gago-Silva's user avatar
  • 2,546
5 votes
1 answer
283 views

Understanding regression function output in ArcGIS Kriging?

Below is "regression function" output in ArcGIS Kriging. I'm assuming the first value is my r2, but does the 2nd value represent error rate, residual, or something else?
Danielle Mitchell's user avatar
5 votes
1 answer
6k views

Variogram model - SAGA on QGIS Processing

QGIS 2.14 (Win) comes with SAGA 2.12, which Ordinary Kriging supports 6 variogram model options. However, if I run SAGA Ordinary Kriging from QGIS Processing Toolbox, there is no drop down menu. Only ...
Kazuhito's user avatar
  • 30.9k
5 votes
1 answer
1k views

Cross-validation of TPS and kriging outputs from the fields package in R

I am looking for existing package or functions that will allow me to cross-validate multiple outputs of Thin plate spline (TPS) and Kriging from the fields package. Does anyone know of any available ...
Gago-Silva's user avatar
  • 2,546
5 votes
1 answer
162 views

Can I "stack" semivariogram data?

I have several synoptic measurements (every few months) of a gravity field I would like to interpolate by kriging. Each group of measurements comprises about 45 measurement locations, and spatial ...
Jeff_K's user avatar
  • 53
5 votes
1 answer
1k views

Extracting the standard error map from autoKrige in R

autoKrige in R works fine and plotting the output gives a good overview of the interpolation with standard error. Following these steps: library(sp) library(automap) loadMeuse() demo(meuse) # ...
Martinyt's user avatar
  • 146
5 votes
1 answer
1k views

Spatio-temporal kriging - useful if not interpolating to new time points?

I have data measuring a variable Z say at several locations on a fixed network (each location at a fixed longitude and latitude). values of Z are measured daily at the same time. At the moment I am ...
piyushnz's user avatar
5 votes
1 answer
590 views

Validity of ordinary kriging of polygon centroids

I'm in a discussion with a colleague who believes that kriging should not be used to interpolate a surface from polygon based data (i.e. with kriging applied to centroids), because this effectively ...
geotheory's user avatar
  • 526
5 votes
1 answer
1k views

Preparing dataset to perform co-kriging in R gstat?

This is the first time I'm using co-kriging in gstat. My problem is that I'm not sure how to prepare the data frame to supply to co-kriging when the variable of interest and auxiliary variables are ...
rm167's user avatar
  • 469
5 votes
1 answer
226 views

Understanding ArcGIS Geostatistical kriging?

We use soil samples data for interpolating soil fertility parameters using ArcGIS geostatistical Kriging. While doing kriging, the resultant map values ranges are not matching with our input samples ...
Pradeep Balan's user avatar
5 votes
0 answers
536 views

Dealing with duplicate locations when kriging by replacing them with single location taking mean of duplicates using R [duplicate]

I understand that the existence of duplicate points (locations) when Kriging will result in a singular covariance matrix. The only way I've encountered people dealing with this is to keep the first of ...
Tom Newton's user avatar
4 votes
1 answer
6k views

R 'gstat - Warning: singular model in variogram fit

I am trying to perform a kriging interpolation. But when it comes to the point when i want to fit the variogram to the model I get the following warning: In fit.variogram(vario, vgm("Sph")) : singular ...
obrob's user avatar
  • 887
4 votes
2 answers
4k views

How to achieve parallel Kriging in R to speed up the process?

I have tried the gstat package in R to perform bathymetry predictions based on Kriging interpolation methods. However, the process is very time consuming due the large amount of sampling data I'm ...
Guz's user avatar
  • 3,176
4 votes
2 answers
5k views

Using kriging to extrapolate values outside of sampling polygon in ArcMap?

I have data points with values along 4 perpendicular transects. The kriging tool that I am currently using interpolates within the polygon created by the 4 transects, essentially a rectangle. I am ...
jcronin's user avatar
  • 71
4 votes
2 answers
2k views

Kriging values greater than the most extreme value of the input layer?

I'll admit I don't understand kriging to the fullest. Yet, I've read that the weights should add up to 1 as in IDW and Natural Neighbors. However, surrounding the biggest elevation points of my ...
Blerg's user avatar
  • 646
4 votes
1 answer
4k views

Error in In predict.gstat in R

I am trying to do kriging by gstat package through R. Here is my Code: krg<-krige(var ~ 1,data,grid,model=vgm,nmax=30) But each time I recive this warning: 1: In predict.gstat(g, newdata = ...
mohammad gerami's user avatar
4 votes
1 answer
1k views

Gstat krige() error: Covariance matrix singular at location [917300,3.6109e+06,0]: skipping

Recently, I am facing a strange problem when using krige() in gstat library. This error was not there a few months ago because I am using the same code with the same data that I successfully used ...
Asad Ali's user avatar
  • 141
4 votes
1 answer
840 views

Is co kriging appropriate for interpolation of three spatial co variates?

I want to estimate an interpolation model(s) for daily minimum temperature, daily maximum temperature and daily rainfall. I have samples of approximately the same size for all three variables at ...
emDiaz's user avatar
  • 159
4 votes
1 answer
2k views

Keeping unchanged observed values after kriging interpolation

Is there a way to force kriging to keep unchanged measured values of observation points? For example: if the observed value of a location is 1, after the interpolation it could be 0.9 or even 0.5, ...
e-falcon's user avatar
4 votes
1 answer
2k views

book/tutorial on using rpy2 in python for kriging?

I am proficient in Python, but have no knowledge of R. Can you suggest a good tutorial in using rpy2 or any other R library to perform Block/Cokriging "Kriging"
user1142937's user avatar
4 votes
0 answers
198 views

Interpolating species data from sample points in a small area [closed]

I have sample data from a small site (10 acres) which shows the number of occurrences of a particular species of snail over two sampling periods (fall and spring sampling). I want to try and ...
user37710's user avatar
  • 103

1
2 3 4 5 6