Questions tagged [kriging]

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.

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…
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 +…
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 …
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…
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…
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…
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…
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…
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…