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
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…
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…
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
1 2
3
12 13