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
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…

RedWest
- 1
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 =…

Elena Iakimova
- 57
- 6
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