Stack Overflow
Stack Overflow

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
2
votes
2 answers

Is it possible to add lines/text/minor tick marks/personalized axes to a variogram plot?

I'm trying to add text, lines, minor tick marks, and/or personalized axes to a variogram using the plot.variogram function. When I try to add any of these using the traditional function calls (i.e. text("Text Here")) it returns the error of plot.new…
Jared Smith
  • 280
  • 3
  • 12
2
votes
1 answer

Spatial temporal kriging

I'm working on rainfall data. I have grid of rainfall points of 1 square kilometer with time step of 10 minutes. I want interpolate those point at every time step. also I want to check the temporal correlation of my rainfall points. The program I'm…
Herve
  • 21
  • 1
1
vote
1 answer

Transformation and export from R of a ggplot map to raster file (tiff)

In R, and RStudio, I have run his code: #ordinary kriging mydata.krig <- krige(Clay~1, train_data, newdata=mask, vgm) names(mydata.krig) names(mydata.krig)[1] <- "Clay.pred" names(mydata.krig) min(mydata.krig$Clay.pred, na.rm=T);…
Dimitris K
  • 33
  • 5
1
vote
0 answers

Non-conformable arrays error in Regression Kriging in R

Grretings to all! In the Regression Kriging, i want to create a map for Standard Deviation of the Ordinary and the Regression Kriging. The pH.sd (from the OK, Ordinary Kriging) are in a matrix: double [250292] The pH.regkrig.sd (from the RK,…
Dimitris K
  • 33
  • 5
1
vote
0 answers

Cluster validation and kriging interpolation in R

I am stuck with the kriging interpolation analysis and colour-coding the silhouette output for correct representation of clusters. Any assistance in this regard is greatly appreciated For the kriging, I don't know how to proceed on this. For the…
Taiwo Fakorede
  • 11
  • 1
1
vote
0 answers

How to create a variogram

I'm a new user of R, which I use in my studies in digital soil mapping. I have my data, as you can see in my .csv file. There are columns SN (sample number), SAND, SILT, CLAY, OM (Organic Matter). I am trying to create a variogram for each property.…
Dimitris K
  • 33
  • 5
1
vote
0 answers

STplot package: How to slice y-coordinate when using mode 'xt' for spatio-temporal kriging?

I'm trying to interpolate my dataset using spatio-temporal kriging approach from gstat and spacetime packages. For the plotting of the results, I use STplot function. I want to use the mode 'xt' which produce the space/time plots. But when I read…
Abdul Azim
  • 11
  • 2
1
vote
0 answers

Different results between kriging on residuals and universal kriging (R)

I have a dataset with 132 observations, 1 response variable (Lopend_gemiddelde), and 9 predictors on which I perform three methods: 1) multiple linear regression (MLR), 2) kriging on the residuals (based on a multiple linear regression), and 3)…
Foeke Boersma
  • 95
  • 6
1
vote
0 answers

Universal kriging has a long processing time compared to other kriging methods - R

I have a dataset data_sp with 132 observations containing 9 predictors and 1 response variable. I perform a twenty-fold cross validation, using ordinary kriging and universal kriging. The raster grid100 on which I do my kriging, contains 61124…
Foeke Boersma
  • 95
  • 6
1
vote
1 answer

Missing output while kriging large point dataset

Long time reader/code parasite, first time question asker! I'm running into a relatively mysterious error while kriging a large point dataset that I'd love any insight people might be able to bring. In essence I'm trying to interpolate bathymetry…
Steven Brownlee
  • 21
  • 2
1
vote
0 answers

Kriging with gstat: including a random effect

I am getting acquainted with gstat and kriging more generally in R. The end goal is to interpolate temperature across my spatial grid, and I want to be able to do so for different time points (average daily temperature over several years). So, I've…
Johanna Fornberg
  • 11
  • 3
1
vote
0 answers

Universal kriging in r. The function predict.gstat return all values as NA

I'm trying to do a universal kriging to get a plot with the predicted values ​​inside a map of a data state. However, after fitting the model and calling the predict.gstat function, the var1.pred values ​​are all NA. It is likely that the error is…
Ivan Bezerra Allaman
  • 79
  • 1
  • 9
1
vote
0 answers

Scikit-learn: How to specify multiple input kernels in gaussian process regression?

This is what I am trying: kernel = C(np.ones(11), (1e-3, 1e3)) * RBF(length_scale=np.ones(11), length_scale_bounds=(1e-5, 1e3)) + WhiteKernel( noise_level=np.array(11*[1e-4]), noise_level_bounds=(1e-5, 1)) gp =…
win_rose
  • 11
  • 3
1
vote
2 answers

how to get better Kriging result graphs in openturns?

I performed spherical Kriging, but I can't seem to get good output graphs. The coordinates(x, and y) range from around around 51 latitude and around 6.5 as longitude my observations range from -70 to +10 here is my code : import openturns as…
simonericmoon
  • 43
  • 5
1
vote
0 answers

kriging approximation with ILNumerics F#

I need to estimate a multivariate function, known from discrete scattered data. I am using ILNumerics interpolation toolbox for that. I have the following code to test the library from F#: let inline (!) (x :RetArray<'T>) =…
pacta_sunt_servanda
  • 724
  • 4
  • 14
1 2 3
…
12 13