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.
275
questions
44
votes
9
answers
30k
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 ...
- 16.3k
27
votes
2
answers
42k
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 ...
- 859
19
votes
2
answers
12k
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 ...
- 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 ...
- 153
11
votes
1
answer
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 .....
- 837
11
votes
1
answer
979
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). ...
- 501
10
votes
1
answer
870
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. ...
- 101
10
votes
3
answers
5k
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 ...
- 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 ...
- 205
9
votes
1
answer
3k
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::...
- 261
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 ...
$ ...
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 '...
- 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) ...
- 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 ...
- 861
8
votes
1
answer
461
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) ...
- 2,315
8
votes
1
answer
385
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....
- 2,546
8
votes
1
answer
219
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 ...
- 101
7
votes
2
answers
5k
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 ...
- 217
7
votes
2
answers
4k
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 ...
- 10.9k
7
votes
2
answers
19k
views
Performing Kriging in QGIS
Does anyone know how to perform kriging in QGIS? and not via SAGA but actually in the program itself.
- 171
7
votes
3
answers
2k
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 ...
- 71
7
votes
1
answer
403
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 ...
- 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 ...
- 20.9k
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?
- 69
6
votes
2
answers
4k
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 ...
- 61
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 ...
- 271
6
votes
1
answer
1k
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 ...
- 20.9k
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, ...
- 103
6
votes
1
answer
698
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, ...
5
votes
2
answers
1k
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 ...
- 2,546
5
votes
1
answer
255
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?
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 ...
- 30.2k
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 ...
- 2,546
5
votes
1
answer
158
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 ...
- 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)
# ...
- 136
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 ...
- 61
5
votes
1
answer
556
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 ...
- 538
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 ...
- 427
5
votes
1
answer
215
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 ...
5
votes
0
answers
535
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 ...
- 173
4
votes
2
answers
3k
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 ...
- 3,126
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 ...
- 71
4
votes
2
answers
1k
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 ...
- 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 = ...
- 189
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 ...
- 141
4
votes
1
answer
751
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 ...
- 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, ...
- 59
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"
- 43
4
votes
0
answers
195
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 ...
- 103
3
votes
1
answer
5k
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 ...
- 837