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

Ordinary kriging returns a prediction frame with identical predictions at every location

So I have two datasets with latitude, longitude and depth in them and I am trying to krige them so that I create a very high resolution bathymetric map for a study I am conducting. I am using the depth information that I collected to construct the…
G_man932
  • 1
  • 1
0
votes
0 answers

Conditioning with External Drift Kriging using GSTools (and Pykrige...?)

I want to impute missing data on abortion rates of US-Counties. I´m comparing different Methods, one of them is Kriging. I used Pykrige for External Dirft Kriging: def Kriging_with_SGS(df, variable, kriged_var, variogram, lags, drifts=None): …
Lorenz
  • 1
  • 3
0
votes
0 answers

Pykrige: Kriging with Specified Drifts

I can´t figure out how to include specified drifts in the interpolation. My code is the following: x = np.array([0, 1, 2, 3, 4]) y = np.array([0, 1, 2, 3, 4]) values = np.array([1, 2, 3, 5, 4]) drift = np.array([-14, 2, 6, 1, 3]) gridx =…
Lorenz
  • 1
  • 3
0
votes
0 answers

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

[![enter image description here][1]][1] .[![enter image description here][2]][2] Can you please help me with this situation? [1]: https://i.stack.imgur.com/Xxp7D.png [2]: https://i.stack.imgur.com/XUm8B.png
Dimitris K
  • 33
  • 5
0
votes
1 answer

Error in ST kriging of uneven spacetime data " the leading minor of order 2 is not positive definite"

I am trying to use universal kriging on some irregular data. I keep getting an error: Error in chol.default(A) : the leading minor of order 2 is not positive definite When I run the krigST command. Here is a sample of my…
Candy
  • 1
  • 4
0
votes
1 answer

Spatial Interpolation Map Temperature using ggplot in R

I would like to create 3 maps of spatial interpolation of a country (like belgium) using the GADM geopkg data that can be found there : https://gadm.org/download_country.html. The idea would be to create 3 interpolated maps : Voronoi…
0
votes
0 answers

gstat: simple vs. ordinary kriging

sorry for such a trivial question but I'm trying to find my way into the "kriging wild bush", a topic that is sometimes really confusing (at first); do you think my interpretation of the difference between simple (sk) vs. ordinary kriging (ok) and…
maxbre
  • 161
  • 9
0
votes
0 answers

MemoryError: Unable to allocate 7.59 GiB for an array with shape (1018651816,) and data type float64

Hi i tried to run the kriging code from pykrige. However, during running the code there is an error. Here is the code. p/s I am new to python **os.chdir("E:\\18_Jupyter") os.getcwd() pip install pykrige import matplotlib.pyplot as plt import numpy…
zam
  • 3
  • 3
0
votes
0 answers

ValueError: A 2-dimensional array must be passed

OK = OrdinaryKriging( x, y, sss, variogram_model='exponential', verbose=True, enable_plotting=True, ) when running it, i have problem/error. Here is the error ValueError: A 2-dimensional array must be passed. How can i clear the error
zam
  • 3
  • 3
0
votes
0 answers

Kriging and extract data

I have one month of data that I converted to STDFDF. Shapiro.test and Adf.test talk about stationarity and normal distribution of natural log transformed dataset. #set the spatial location of the regions Data.loc <- data.frame(Stat =…
0
votes
1 answer

Apply pykrige to interpolate a high-dimensional data

I am going to realize a black optimization problem, where one step is applying Kriging interpolation method to construct a approximate model. I have referred to the PyKrige Documentation and do not find any class supporting interpolation of…
Jeff Ju
  • 3
  • 1
0
votes
0 answers

scikit GPR returning 0 at points not used to train it

I am trying to use a gpr to model the behaviour of some data, but I am finding that after I have fitted the data, when I try to predict at certain points I get a result of 0, unless I am at/very close to a point used to train it. My GPR is given…
0
votes
0 answers

Wondering if my code is correct for Simple Kriging

I am a student and we learned about simple kriging with an example of 3 points with known elevations and 1 point with unknown elevation with the assumption that the empirical semivariogram is represented by the linear regression line with Y…
Frank D
  • 21
  • 1
0
votes
0 answers

Spatial interpolation using kriging in R

I have two datasets. The first one shows information about multiple weather phenomena in Brazil measured by weather stations in the country. I also have information regarding the latitude and longitude of these stations, and the weather data is…
Mateus Maciel
  • 151
  • 1
  • 1
  • 10
0
votes
0 answers

How to confine grid for kriging to shapefile or polygon boundaries? Help needed

I'm working on creating a grid for kriging in R. I am trying to limit my prediction grid to the water, and not include anything outside of it. I have a shape file for the water shoreline. I've only been able to create a convex shape from my data…
hydrotime
  • 1
  • 1