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

Plotting Probability Density Heatmap Over Time in R

Let's say I have the output of a monte-carlo simulation of one variable over several different iterations (think millions). For each iteration, I have the values of the variable at each point in time (ranging from t=1 to t=365). I would like to…
Alfalfa
  • 131
  • 7
3
votes
1 answer

how to change the lag distance while calculating a Variogram in R

I am trying to calculate an Experimental-Variogram value at different lag distances, so I am using the variogramm command variog1 <- variogram((Copper)~1,ds) but I can't know how to specify the needed lag distance. For example I want to get a…
Ahmed Emam
  • 33
  • 6
3
votes
1 answer

How to apply kriging for 3D arrays in Python?

I have a 3D numpy array with some elevation values. I would like to apply kriging interpolation method to them and get a full valued array with same given shape. My purpose is to create a surface actually. The values, I have, are hydrogelogical…
Mustafa Uçar
  • 442
  • 1
  • 6
  • 18
3
votes
0 answers

Variogram Iteration Warning

I would like to perform kriging using a variable call "Secchi.Disk" from a data set I have. However, when fitting the variogram I get the following warming message Warning message: In fit.variogram(v, vgm(0.5, "Sph", 40, 0.5)) : No convergence…
Perro
  • 348
  • 5
  • 14
3
votes
1 answer

Universal kriging using lat long gstat R

I'm new at R and I'm having some trouble to perform a universal kriging with gstat R. As Hengl et al. (2004) say "Universal kriging should be reserved for the case where the drift (or trend) is modelled as a function of the coordinates only". So, I…
Fjord
  • 144
  • 1
  • 9
3
votes
2 answers

Underlay a vector image to a grid for kriging in R

After searching around a lot, asking, and doing some code, I kinda got the bare minimum for doing kriging in R's gstat. Using 4 points (I know, totally bad), I kriged the unsampled points located between them. But in actuality, I don't need all of…
ace_01S
  • 387
  • 2
  • 5
  • 16
3
votes
1 answer

extrapolation using Gaussian processes regression or Kriging

Is there any way to estimate the extrapolation using kriging or Gaussian processes regression ? Gaussian processes work very well for interpolation of scattered data; however, I need to extrapolate a time series of variable in time. hoe can I…
shahrouzalimo
  • 41
  • 1
  • 2
3
votes
0 answers

Spatial interpolation for scattered data

I have a list of metars points (latitude, longitude, sea level pressure) and need to interpolate the pressure. Suppose I want to iterpolate the red point (on the kriging side of the image bellow) using the known measurements in the list of metars…
João Azevedo
  • 343
  • 1
  • 10
3
votes
1 answer

Time taken to krige in gstat package in R

The following R program creates an interpolated surface using 470 data points using walker Lake data in gstat package. source("D:/kriging/allfunctions.r") # Reads in all functions. source("D:/kriging/panel.gamma0.r") # Reads in…
Chandan
  • 764
  • 2
  • 8
  • 21
3
votes
1 answer

In R, kriging gives "system is exactly singular" error

I'm having an issue where my specific inputs generate an error. library(kriging) x <-…
Andrew
  • 7,201
  • 5
  • 25
  • 34
3
votes
2 answers

Could someone describe a 2d interpolation method that is better than bilinear interpolation?

I have a grid of data points that I currently use Bilinear interpolation on to find the missing points in the grid. I was pointed in the directions of Kriging aka thee best linear unbiased estimator, but I was unable to find good source code or an…
Bob
  • 1,065
  • 2
  • 16
  • 36
3
votes
0 answers

how to set the values of sill,nugget,range in ordinary kriging?

im work on temperature data ,and i want to interpolate the data to any location that i specify . my programming language is asp.net4 c#. i dont know how can i determine the sill,range,nugget for experimental variogram or theoratical variogram? i…
The_programmer
  • 37
  • 5
2
votes
1 answer

Get Sobol indices error from kriging with Openturns

I'm looking for a Python package which computes the Sobol' indices from a kriging model and which provides a confidence interval on these indices that takes into account both the meta-model and the Monte-Carlo errors. Is it possible to do this with…
Alexis
  • 21
  • 3
2
votes
1 answer

"CRS object has comment, which is lost in output" in krige

I am trying to krige for water quality dataset with Latitude, longitude which using CRS("+init=epsg:4326"). GGT <- read.csv("C:/Users/user/Data/newdata2019.csv") coordinates(GGT) = ~Lon+Lat GGT <- st_as_sf(GGT) st_crs(GGT) <- 4326 GGTgrid <-…
Heeone Lee
  • 43
  • 1
  • 6
2
votes
1 answer

spatial interpolation with gaussion process regression

I have a csv-file with 140.000 points(rows). It consists of: longitude value latitude value subsidence value at specific points. I assume that these points are spatially correlated. I want to perform a spatial interpolation analysis of the area of…
simonericmoon
  • 43
  • 5
1
2
3
…
12 13