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.

Filter by
Sorted by
Tagged with
0
votes
0answers
21 views

Converting distance from km to degrees (WGS84)?

Is there an equation in which I can convert distance in kilometers into degrees? I am performing a Kriging analysis and when I am measuring the distance between the input points (for LAG and Major ...
0
votes
1answer
21 views

What are the formulas used in SAGA Ordinary Kriging?

can some one name the formulas avaible in the Ordinary Kriging Modul of SAGA: What ist the exponential, spherical, Gaussian-Model, etc.?
1
vote
1answer
79 views

R kriging raster returns NA for some prediction points?

I am trying to interpolate a raster layer of chlorophyll-a data using kriging (gstat::krige) in R, however the resulting interpolation is showing NA values in a peculiar straight line through the ...
0
votes
0answers
16 views

Cokriging cross validation variances too small?

I'm looking into kriging the Meuse dataset in R. I first performed ordinary kriging using copper, using krige.cv to calculate residuals. The distribution of prediction variances given by applying ...
0
votes
1answer
45 views

After I give kring command, this error comes up: data item in gstat object and newdata have different coordinate reference systems

> library(gstat) > library(readxl) > library(sp) > library(rgdal) #2013(FATALITIES) data <- read.csv("C:/Users/sbs/Desktop/Study/RA/Accident analysis/New folder/Iowa total crashes ...
1
vote
1answer
21 views

Estimation value of simple kriging are constant in geostatistical analyst (ArcGIS)

The objective of my study is to compare the accuracy of the model SK (simple kriging) and OK (ordinary kriging). When I use OK to apply on my data, there is a little bit variation on the prediction ...
0
votes
0answers
9 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
17 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
38 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
338 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)...
-2
votes
1answer
44 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
167 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
361 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
73 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
184 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
203 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
96 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
40 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
128 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
106 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
143 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
138 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
71 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
157 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
150 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
184 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
1k 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
136 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
253 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
692 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
18 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
54 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
407 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
80 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
68 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
501 views

How to calculate kriging weights?

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
87 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
183 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
33 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
86 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
574 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
2k 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
491 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
249 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
158 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
175 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
343 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
203 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 ...