Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [kriging]

A group of geostatistical techniques to interpolate the value of a random field at an unobserved location from observations of its value at nearby locations.

0
votes
0answers
8 views

Cross validating part of spatial interpolation?

I am very new to ArcGIS Desktop . I am exploring some of the spatial interpolation techniques like IDW and Kriging. My query is on the cross validation part of the techniques. In the cross validation ...
1
vote
0answers
16 views

Using convex hull for interpolation of individual plots or lump plots together in ArcMap?

This question is inspired by Using kriging to extrapolate values outside of sampling polygon in ArcMap? I'm working on a project trying to understand the lasting effects of historic scientific ...
1
vote
0answers
29 views

Strange output of Kriging standard deviation

For my thesis I am interpolating some 4000 boreholepoints via Kriging to try and visualize a geological layer under the ground level. Since I know there is a strong trend in my data (SW-NE slope) I ...
2
votes
2answers
313 views

Making data table using linear kriging interpolation with variogram in R

I have a table of data with such structure: longitude_for_data = c(32,68,89,145,176, -14, -42) latitude_for_data = c(22,98,21,13 , 16,-134,-102) Z_for_data = c(30,40,60,70 , 20, 40, 100)...
-1
votes
1answer
43 views

Does the sampling design affect the variogram and mapping? [closed]

I am not sure if this is the right place to ask the question, but I don't know where to ask. So please excuse me if not so. I am considering two sampling strategies, the first is locating samples on ...
1
vote
0answers
98 views

QGIS 3.6 ordinary kriging produce one constant output

When I run kriging a+bx or any other model, I always get just one constant output. I have selected collect attribute. Prediction and quality layer has just one value. see output here:
0
votes
1answer
317 views

How kriging variance is calculated in R gstat?

How does exactly the krige function (a wrapper to gstat and predict functions) from package gstat calculate kriging variance (minimized estimation variance) in Ordinary Kriging? I wanted to reproduce ...
1
vote
1answer
41 views

Kriging with variable spatial precision?

I have a data-set comprised of hydraulic conductivity estimates measured at groundwater monitoring wells, as seen here: Exact GPS coordinates are available for some wells, while for others I only ...
1
vote
0answers
65 views

Error in Universal Kriging using gstat package in R [closed]

I am trying to perform universal kriging in R using the gstat package and running into the following error: krige(formula = CDD_LT_stC ~ Imperv1km + Elev_new, locations = sp.cal_data, newdata = sp....
0
votes
1answer
147 views

Choosing interpolation method for grid data?

I have a grid data of evenly spaced points that gives the number of people with a facebook account within a certain radius (approximately equal to the mid-distance between to points). From that I'd ...
0
votes
1answer
177 views

Kriging method vs. IDW in soil moisture interpolation

I have XY coordinates and soil moisture values for a grid of 10m and I used add XY data to import my excel data to ArcMap 10.5. I applied both IDW and kriging methods for interpolation. The results ...
0
votes
1answer
78 views

Which is the model semivariogram for binary data?

I have a binary data (1,0) to calculate spatial distribution of soil bacteria. I have to use indicator kriging (presence/absence) and fitting model of semivariogram. Then, which is the model ...
1
vote
1answer
38 views

Kriging for unevenly spaced points using R?

I started to read about Kriging and the implementation into R. In all the tutorials I found, a (regular) grid is used to do ordinary or universal kriging. My situation is as follows, I try to ...
1
vote
1answer
105 views

IDW/Kriging running slow in R

I have a dataset with 2157335 entries, and I am performing IDW/Kriging to interpolate. But the code is running very slow. I am not sure this is because of a large dataset or there is something wrong ...
0
votes
1answer
90 views

SAGA GIS Kriging in QGIS calculating wrong values

I have a point shapefile concerning rainfall for the whole country. The points represent the annual rainfall in mm and range from 675-3562. I have approximately 800 points. I am attempting to use ...
1
vote
0answers
119 views

Interpolation in QGIS shows wrong output result [closed]

I have done interpolation (IDW and Nearest Neighbour) in QGIS. But interestingly the output shows the wrong result. The area where I have a huge number of points results in a lower Z value (it is ...
5
votes
1answer
131 views

Cokriging: zero distance cross-semivariance in undersampled cases [gstat error]

Trying to cokrige two variables that are not perfectly colocated (one has sparser measurements than the other), I faced an issue that I'll illustrate with the following MRE. With the meuse dataset, ...
0
votes
0answers
66 views

Cokriging, Zero distance semivariance in cross-variogram error

Trying cokriging with simulated data, I faced a problem that did not seem one in the demo(cokriging) with the meuse dataset: I can't use fit.lmc() default but need to specify fit.method=1 as the ...
0
votes
1answer
145 views

Calculating coefficient of variation from gstat output to create uncertainty map in R

I was asked by my supervisor to create uncertainty maps of kriging interpolation based on the coefficient of variation ((sd/mean)*100%). I used krige function from gstat package to perform the ...
0
votes
0answers
124 views

Interpolate Multiple Data With Kriging In ModelBuilder?

I am having problems building a model for running interpolation on multiple data using the Kriging method. I have followed these steps: Convert multiple data to geodatabase file Used iteration on ...
5
votes
1answer
1k views

Kriging example with sf object

I am new to geo spacial analysis. I would like to learn kriging from a simple example. Say, we have four points with some z value. library(sf) library(gstat) library(tidyverse) dt<-tibble::...
2
votes
0answers
167 views

Empirical variogram plot in QGIS

I am trying to do kriging on mortality rates using QGIS (I have version 3.0 but I downloaded the 2.18.20 to be able to follow the web documentation and most of the answers here). From what I have ...
1
vote
0answers
930 views

Interpolation with QGIS 3

Resulting map is not generated after performing following steps: Add vector layer(shapefile of a country) Add demlimited text layer(.csv file of values/variables) Click on processing toolbox Type "...
0
votes
1answer
126 views

spatio temporal anisotropy in R: stAni argument in vgmST and krigeST

I have some problems in understanding the usage of argument 'stAni' of vgmST and krigeST functions, for the setting of spatio temporal anisotropy. How have I to interpret this argument, that is the ...
4
votes
1answer
224 views

Preparing dataset to perform co-kriging in R gstat?

This is the first time I'm using co-kriging in gstat. My problem is that I'm not sure how to prepare the data frame to supply to co-kriging when the variable of interest and auxiliary variables are ...
3
votes
0answers
623 views

Difficulty using ordinary kriging in QGIS 3

I have a shapefile with 65 data points. Each data point has a value for NO2 pollution. I would like to use kriging to create a grid of NO2 pollution values for the whole area. I have tried using ...
1
vote
0answers
17 views

Kriging with a non-uniform (not space filling with concentrated data in certain regions) sampling design?

I am using Kriging for optimization where local searches are often carried out. However, because of these local searches my sample design becomes poor and my Kriging prediction as well as variance ...
0
votes
1answer
51 views

Kriging: Large standard error very close to known value points

I am attempting to create surfaces of strata via Kriging in Arc GIS from borehole data. In addition to the predicted surface I also want surfaces showing 95% confidence above and below the predicted ...
-1
votes
1answer
373 views

Performing Cokriging in R? [closed]

I would like to perform cokriging of 60 observed data points. In the shapefile, there are 60 fields containing the observations as Time 1 to Time 60. I also have 60 rasters in tiff format. I want to ...
0
votes
1answer
69 views

Kriging interpolation with blocks as barriers for street level visualization in ArcMap?

I'm trying to do Kriging interpolation with blocks as barriers for street level visualization, but I can't realize how I can do this. I already have the barriers polygongs, data and all. But can I do ...
0
votes
1answer
61 views

Kriging with uncertain observational data

I'm currently trying to use ArcMap 10.5.1 for this analysis: Thus far, I have taken particulate matter trend estimates at monitoring sites across the U.S. and used Simple (Normal Score - due to ...
5
votes
2answers
439 views

How weights in Kriging are calculated, and predict beyond range of sample values?

In Esri's documentation on how to make a prediction in Kriging, I couldn't figure out the linkage between the semivariograms and the prediction process from reading the description therein: After ...
1
vote
0answers
85 views

Applying mask to Kriging shape in ArcGIS Desktop?

I'm trying to apply an mask to a Kriging image for after apply a Density Sampling Network, but unfortunately I'm not getting good results. I tried to apply on mask on Density Sampling Network and is ...
1
vote
1answer
167 views

how to predict on hold on test sample using krige.cv from gstat R toolbox?

I am using the 'krige.cv' function from the gstat toolbox in R. It does not take any 'newdata' as input for prediction ( as in universal kriging). How can I predict using the fitted model from this ...
1
vote
0answers
31 views

Is it possible to to have higher variances for predictions in Universal Kriging?

I am using gstat to do Universal kriging. I am finding that the variances for the predictions of a hold out sample (i.e., new spatial locations for which I want to determine the response variable) ...
1
vote
0answers
84 views

How to obtain nugget parameter from fit.variogram in gstat? [closed]

I have the following code where I fit a variogram using gstat. vA = variogram(dep_var~elevation=x+y, data, cressie=TRUE) vmodel = vgm("Sph") vB = fit.variogram(vA,vmodel) How can I get the nugget ...
1
vote
2answers
533 views

Unable to Export Kriging Prediction Map

I have performed Kriging using the autoKrige function, that is part of the automap package in R. I obtained this as my output I have realized that maybe the reason why I can not export the ...
1
vote
1answer
1k views

Fit 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,...
4
votes
1answer
447 views

Extracting the standard error map from autoKrige in R

autoKrige in R works fine and plotting the output gives a good overview of the interpolation with standard error. Following these steps: library(sp) library(automap) loadMeuse() demo(meuse) # ...
3
votes
0answers
224 views

Dealing with duplicate locations kriging - a method for replacing with a single location taking the mean of the duplicates? (R)

I understand that the existence of duplicate points (locations) when Kriging will result in a singular covariance matrix. The only way I've encountered people dealing with this is to keep the first of ...
0
votes
0answers
147 views

Appropriate method for transforming log-transformed kriging predictions back to original scale (R)

I am attempting to run some interpolations using kriging on transformed data using the transformation log (x + 1) - because my non-negative data has many zero values. The problem I have is with back-...
2
votes
0answers
167 views

What should be different variogram model in kriging [closed]

In ordinary kriging algorithm in QGIS (Saga algorithms provided in QGIS) there is field for Variogram model which by default is given as 'a + b * x' what does it signify and what are other options to ...
3
votes
1answer
315 views

kriging vs cokriging

I am attempting to interpolate sediments across an area with data from 3 sources using ArcGIS 10.3. All three sources divide sediments into percentages of mud, sand or gravel. Two sources come ...
1
vote
1answer
176 views

DEM, points values and elevation to create a map of temeperature

I am using SAGA - Universal Kriging (Global) for this, but I am stucked. I have grid (.tif) for Serbia with elevation values, and .shp file with points that have mean yearly temperature. I want to ...
2
votes
1answer
572 views

R kriging cross validation returns NA for all prediction points?

I'm trying to cross-validate the results of an [auto]Krige output but for some unknown reason the results of both autoKrige.cv (automap package) and krige.cv (gstat package) return NA for all values (...
1
vote
1answer
128 views

Creating prediction map by combining several rasters and vectors in ArcGIS Desktop?

I have three 'variables': 6000 GPS points that represent where my focal species (a carnivore) is present. Attached to these points is one value representing the activity level of the animal at the ...
1
vote
1answer
603 views

Error in Kriging function

I am using gstat::krige to interpolate a data set of arund 2700 points. However I get an error message "Error in seq.default(zrng[1], zrng[2], length.out = cuts + 2):'from' must be a finite number" ...
2
votes
1answer
364 views

How to make ordinary krigging with a spherical variogram in SAGA GIS?

I have a very special map to make: it's a sea bird density map. To do this, I interpolate the observation points with ordinary kriging with a variogram. As I've read, the variogram must be of "...
1
vote
1answer
43 views

Missing argument in SGS simulation of HPGL

Here is a simple example to test the Sequential Gaussian Simulation in open geostatistics library HPGL. I use the following file as hard data for the kriging: NEW_TEST_PROP.INC here is my simple ...
1
vote
0answers
171 views

How to do kriging on a country map from a linear model

I would like to use ordinary kriging to interpolate model predictions on a grid. In my code underneath, GPS is a 'SpatialPointsDataFrame' and EU is a 'SpatialPolygonDataFrame'. All covariates for ...