A statistical interpolation method, also known as Gaussian process regression, most used in geo-statistics. The goal is to map a surface given limited sample data. The process evaluates the variability of supplied data, then uses a weighted average of neighbouring points -- considering both distance and direction -- to interpolate the desired map points.
Questions tagged [kriging]
186 questions
0
votes
0 answers
GaussianProcessRegressor fitting perfectly but poor perfomance on test data?
I am trying to understand GPR, and I am testing it to predict some values. The response is the first component of a PCA, so it has relatively good data without outliers. The predictors also come from a PCA(n=2), and buth predictors columns has been…

ger.code
- 13
- 3
0
votes
0 answers
Prediction output for kriging in QGIS not working
I am fairly new to QGIS and am trying to use the kriging function. For context, I have a dataset for average traffic in Scotland, with traffic being counted at specific count points along all of the roads. I want to use the kriging function to…
0
votes
2 answers
Change krig plot with a continuous scale to discrete
I'm trying to solve a krigging/ggplot issue. I currently have the plot showing a continuous scale using:
Cu_DTPA_NL.kriged %>%
as.data.frame() %>%
ggplot(aes(x = x, y = y)) +
geom_raster(aes(fill = var1.pred), interpolate =…

RHT
- 21
- 3
0
votes
2 answers
Iterate over multiple dependent variables (columns) of an sp dataframe when using krige.cv
I have a SpatialPointsDataframe called rain and I would like to fit a variogram and perfom cross-validation for each one of its last 10 columns (dependent variables) like below:
fit.reg.vgm <- autofitVariogram(
column (dependent variable) ~ X +…

Georgios Boumis
- 25
- 4
0
votes
1 answer
autokriging and SpatialPixel grid file
I am trying to use autoKrige in automap package.
My data looks like
head(LL12)
Location Longitude Latitude percent
1 L1 16.48476 -28.60737 1.0000000
2 L3 16.90045 -29.24154 1.0000000
3 L4 16.99268 -28.44906 0.6363636
4 …

Maria Eugenia D'Amato
- 41
- 1
- 4
0
votes
0 answers
error wrong sign in 'by' argument while using grid.expand function in kriging (R)
I am building kriging model in R. I have plotted Lat Longs of city.With the help of locator function I have figured out boundries of city. Now I want to expand this thing which I can do using grid.expand fuction in which we do have "by" argument in…

Shruti Gote
- 3
- 6
0
votes
0 answers
Interpolation using matplotlib tricontourf function in python
I'm fairly new to python programming. The data I'm working with is generated by filtering process via a third party software, such that there are three arrays, x, y and z that don't need to be arranged in a 2d grid and its contour plot can be…

Chandan Dey
- 13
- 7
0
votes
0 answers
How to plot kriging analysis when packages presents conflicts and the option :: doesnt work?
I did a plot with set packages (geoR", "oceanmap","raster"); however, it was possible to build the figure below (the current figure is a print screen because it has a higher size than permitted, sorry by this); It is important to clarify that this…

Juan Carlos Rubio Polania
- 137
- 1
- 2
- 10
0
votes
1 answer
spatial interpolation using Ordinary Kriging
Two different institutions are using the same dataset for spatial interpolation using Ordinary Kriging. However, the resulting maps shows deviances.What are the potential causes of differences in maps?

SSB
- 1
- 2
0
votes
1 answer
Kriging in gstat
I have a problem with code that ran 1000 times without error.
Error message:
"Error in bb[, "max"] : subscript out of bounds"
I have no clue what bb means. I thought bbox but that is same for grid and…

Matevž U. Pavlič
- 1
- 1
0
votes
1 answer
Gstat prediciton with heterogenous measurement errors
I have a set of measurements for a variable Z over a non uniform grid in X and Y. For each Z value I have a mean plus a variance of the measurement. That is, I have a data set (Zmean_i,Zvariance_i,X_i,Y_i) i=1..N. In other words, some estimates of…
0
votes
0 answers
how to make a proper grid in case of a GAM + OK kriging method?
I am struggling with a Kriging + GAM problem since a while and hope someone will be able to help me.
I am trying to make interpolation of a pollutant into an area. To do so, I use ordinary kriging from gstat, in addition to GAM, based on this…
0
votes
1 answer
Exporting Kriged Values from R
Does anyone have any suggestions on how to export the values from your kriged data in R. I have them exported as a raster, but I need the actual cell values. The data that I'm working with are fish densities. I'm kriging the data and then converting…

Kelsey Spencer
- 57
- 1
- 1
- 9
0
votes
0 answers
Create a grid inside shapefile for kriging based on meuse.grid
I was trying to follow this kriging exercise and ended up going on a 6 hour adventure in an attempt to create a grid file like meuse.grid in the link. I'm very new to R and geostats in general so I'd really appreciate any help. I managed to create a…

savemeplease
- 1
- 1
0
votes
2 answers
What is internal drift and external drift?
For the Kriging model, I'm wondering what the difference is between internal and external drift? I understand the internal drift is part of the location (s) and the external drift is typically an additional variable that has a linear correlation…

user11591310
- 9
- 2