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
2
votes
0 answers
Kriging interpolation using GeoStats package in Julia
I am trying to build a model for kriging interpolation using GeoStats package in julia.
I have tried an example of 2D interpolations but the results are not accurate, as mentioned below.
Code for 2D interpolation:
using KrigingEstimators,…

Mohammad Saad
- 1,935
- 10
- 28
2
votes
1 answer
Kriging in Python
I have use pykrige for interpolation-2d.
A few data (x, y, z) which is location and elevation.
But the outcome is not good on ordinaryKriging (‘spherical’).
How can I adjust the parameters to get better outcome.
Or any suggestion for kriging…

Xiang
- 230
- 2
- 10
2
votes
1 answer
Geographic Interpolate in Python
I have 500 points with longitude x, latitude y, altitude z, and the value at these points.
On the other hand, I have other 200 points than I would like to interpolate, where the latitude, longitude, and altitude are known.
I would like to…

Francisco Gonzalez
- 437
- 1
- 3
- 15
2
votes
1 answer
How to interpolate 2D spatial data with kriging in Python?
I have a spatial 2D domain, say [0,1]×[0,1]. In this domain, there are 6 points where some scalar quantity of interest has been observed (e.g., temperature, mechanical stress, fluid density, etc.). How can I predict the quantity of interest at…

Michael Baudin
- 1,022
- 10
- 25
2
votes
0 answers
Issue in kriging some rainfall records using hydrokrige function in hydroTSM package
I have been trying to use hydrokrige function (hydroTSM package) to interpolate rainfall data. For some data it seems to be working fine while I can't make sense of result I am getting for some records. Below are two reproducible examples with and…

Dnstuladhar
- 31
- 5
2
votes
0 answers
object not a >= 2-column array
I am learning to krige in R and am stuck with this error "object not a >=2-column array".
I converted my initial data frame to a spatial data frame by assigning coordinates as follows:
coordinates(PB) <- cbind(PB$x, PB$y)
I created a grid…

Cary Lindsey
- 21
- 2
2
votes
0 answers
Cross-validation for kriging in R: how to include the trend while reestimating the variogram using xvalid?
I have a question very specific for the function xvalid (package geoR) in R which is used in spatial statistics only, so I hope it's not too specific for someone to be able to answer. In any case, suggestions for alternative functions/packages are…

Vera
- 29
- 4
2
votes
1 answer
How to apply universal kriging with custom prediction spatial grid using autoKrige in R
I want to apply universal kriging on a dataset using the autokrige function in R. I would like to create my own custom, spatial grid for the predicted points (for the new_data argument of autokrige). I am using R version 3.2.2 (64-bit) and RStudio…

Henry
- 23
- 4
2
votes
1 answer
Are there any softwares that implemented the multiple output gauss process?
I am trying to implement bayesian optimization using gauss process regression, and I want to try the multiple output GP firstly.
There are many softwares that implemented GP, like the fitrgp function in MATLAB and the ooDACE toolbox.
But I didn't…

Alaya
- 3,287
- 4
- 27
- 39
2
votes
1 answer
How to fit model with individual measurement error in DiceKriging, or can it?
I have a set of 5 data points (x=10,20,30,40,50 and its corresponding response values y and noise as s.d. of y). These data are obtained from stochastic computer experiments.
How can I use DiceKriging in R to fit a kriging model for these data?
x <-…

user2513881
- 69
- 1
- 7
2
votes
1 answer
set automatic colors levels in .filled.contour ? (kriging)
I have a hard time to automatically create colors using filled.contour displaying kriging plot, avoiding specifying levels().
I can plot my results and see the legend, but why my colors are repetitive? Thus, why interval 4.5 - 5.0 has the same color…

maycca
- 3,848
- 5
- 36
- 67
2
votes
1 answer
Interpolate TO a grid FROM random points
I have a series of coordinates that represent data in a 2D scalar field: x[i], y[i], F[i]. However, x[i] and y[i] are not on a rectilinear grid, they are "randomly" placed. I want to interpolate f[i] so I can estimate its value on a rectilinear…

user3533030
- 359
- 3
- 17
2
votes
1 answer
Why are my predicted values constant?
mydata
mydata_x <- rep(seq(200,600,by=200),times=1,each=4)
mydata_y <- rep(seq(600,1200,by=200),times=4,each=1)
mydata_z <- c(0,0,0,0,0,0,0,529,0,0,0,0)
mydata <- data.frame(x=mydata_x,y=mydata_y,z=mydata_z)
I've taken samples from a virtual field…

braxtonlewis
- 41
- 1
- 8
2
votes
1 answer
Kriging / Gaussian Process Conditional Simulations in Matlab
I would like to perform conditional simulations for Gaussian process (GP) models in Matlab. I have found a tutorial by Martin Kolář (http://mrmartin.net/?p=223).
sigma_f = 1.1251; %parameter of the squared exponential kernel
l = 0.90441;…

Nat
- 565
- 4
- 11
2
votes
1 answer
Logistic function for regression kriging
I would like to perform regression kriging (RK) for binary presence-absence and host grid data as a constant predictor. I have used logistic function to estimate relationship between binary outcome and predictor, however I think it is not passing…

mihoo
- 157
- 3
- 12