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.
274 questions
0
votes
0
answers
10
views
Symbology and statistic issues with interpolation maps in Geostatistical Wizard of ArcGIS Pro
I have created a series of depth interpolation maps using ordinary kriging with the Geostatistical Wizard in ArcGIS Pro. To make the maps comparable to each other, I manually adjusted the class breaks ...
- 1
0
votes
0
answers
57
views
Groundwater interpolation using kriging with external drift
I have a bunch of groundwater measurements that I would like to interpolate taking into account the ground elevation. For this, I was advised to use kriging with external drift. I have tried different ...
- 11
0
votes
0
answers
91
views
QGIS SAGA Kriging tools not working with my data
I've been trying to create surfaces using Kriging from the geochemical sample data points I have, using SAGA ordinary kriging. But the tool doesn’t work.
Most of my data points have 50m*100m spacing. ...
- 76
-1
votes
1
answer
246
views
Performing empirical bayesian kriging in R
I am trying to apply empirical bayesian kriging in R, but I cannot find a package that provides a function for EBK.
Does anyone know of any or should I use ArcGIS?
- 1
0
votes
0
answers
80
views
Understanding Universal Kriging results
The results obtained using the Universal Kriging method with 41 data points appear as shown in the graph. I would like to inquire about the reasons behind the occurrence of the sawtooth pattern in the ...
- 1
0
votes
1
answer
78
views
Minimum datapoints for Kriging
Are the datapoints (34) for this Kriging operation sufficient?
Research Gate recommends > 30. [both -> here and here], while another source [-> here] mentions 30 to 50. All mention Webster ...
- 337
0
votes
0
answers
93
views
Using data for interpolation in Geostatistical Wizard of ArcGIS Pro
i do have soil temperature data (5 cm depth) from 1991 and 2021 for Germany, and wanted to interpolate the data for each year as a class of 3 seasonal months (winter = 12, 1, 2; summer = 6, 7, 8). ...
- 1
0
votes
0
answers
52
views
Clip a SpatialGrid using a SpatialPolygon
In order to perform kriging methods using automap::autoKrige I need to provide a SpatialGrid for which values can be predicted. I would like to reduce the rectangular grid to the actual area (in this ...
- 115
1
vote
1
answer
75
views
Improving kriging results from transect data
I have estimates of fish density along transects from a lake and would like to extrapolate these estimates for the entire lake surface, using kriging methods. Based on information from here, here and ...
- 115
1
vote
1
answer
365
views
What does ArcGIS Kriging use lag size for, if semivariogram model (range, sill, nugget) are prespecified
Using the Kriging tool in ArcMap I'm a bit confused on the meaning of parameters.
The interface asks for lag size, and optionally range, partial sill and nugget.
My understanding is that lag size is ...
- 2,343
0
votes
0
answers
60
views
Unexpected artifacts when creating a digital terrain model (DTM) using lidR and gstat libraries in R
I'm using the lidR, raster, sp, and gstat libraries in R to create a Digital Terrain Model (DTM) from LiDAR data. I'm encountering unexpected artifacts in my DTM and I'm unsure why.
Here is my code:
...
- 67
0
votes
1
answer
107
views
Kriging - spatial autocorrelation analysis
I'm reading a paper, where semivariogram analysis was done using kriging to calculate spatial autocorrelation of building hights.
The paper says:
As can be seen from Figure 3, the building in the ...
- 11
1
vote
0
answers
52
views
Generate simulations using variogram with given mean and standard deviation
I am trying to generate simulations using variogram using the R code below
library(sp)
library(gstat)
library(ggplot2)
library(tidyverse)
library(dplyr)
library(reshape2)
# Call data
points <- ...
- 77
0
votes
1
answer
262
views
Interpolate underwater quality data with QGIS as with Surfer with Kriging method [closed]
I'm trying to interpolate conductivity data form a coastal area about 4-6 square km from a series of 10-20 probe wells more or less equally distributed.
I have access to a licensed Surfer and to QGIS ...
- 31
1
vote
1
answer
77
views
Kriging: extract data for other unknown points
I have one month of data that I converted to STDFDF.
#set the spatial location of the regions
Data.loc <- data.frame(Stat = c("Firmo","Citta dei Ragazzi", "Rende",
...
0
votes
0
answers
583
views
SAGA GIS - regression kriging
I have a problem with regression kriging in SAGA under QGIS. It gives me an error - could not execute the tool. I want to create temperature maps from data set of meteorologic station and DEM.
I have ...
- 1
0
votes
1
answer
78
views
Exporting voxel sections to raster in ArcGIS Pro [closed]
I need to perform some algebra in voxel sections resulting from EBK 3D (NetCDF file).
I need to calculate a mean section and standard deviation, and to do that I need to export the sections to raster.
...
- 139
0
votes
0
answers
38
views
Generating vertical sections in 3D models in ArcGIS Pro
I have a model in NetCDF from EBK 3D in ArcGIS Pro and I need to take a vertical section, but ArcGIS only provides horizontal sections, has anyone faced this issue and know how to get vertical ...
- 139
0
votes
2
answers
992
views
QGIS Smart Map Plugin - Error Running Interpolation
I am new to using the Smart-Map plugin. I added my data and did the ordinary Kriging interpolation. Just before it reaches 100% I get the following error message:
An error has occurred while executing ...
1
vote
1
answer
85
views
Optimize summary stats PostGIS query for multiple polyline layers over a worldwide grid
I have two global PostGIS polyline layers, one of which is a subset of the other. The subset defines where there are issues with the system being mapped.
I am trying to map the likelihood of a ...
- 660
2
votes
0
answers
166
views
Standard deviation of interpolated data
I have a dataset with air-quality measurements in several points across one state, through several months. Thus, for each point I can calculate average, variance and so. My goal is to interpolate data ...
0
votes
0
answers
45
views
How to stack/merge geom_stars items/plots in one - in R
If I have plotted the Relative Error Maps for the Silt soil attribute for the following 4 different soil types: Cambisols, calcisols, leptosols and fluvisols that form my study case area:
ggplot() +
...
- 11
0
votes
0
answers
240
views
Spatial Outliers Detection and Remove in R?
I have run my data in R, and created a distribution mapof my data (Clay Soil Property).
I want to detect in a way the spatial outliers that i can see from the map, and remove them in order to execute ...
- 11
0
votes
1
answer
232
views
Spatio-temporal Kriging and how extract specific values
I have several years of data that needs to be converted to STDFDF.
Shapiro.test and Adf.test talk about non-stationarity and non-normal distribution.
particolato_May2021 <- read.delim(file = "...
0
votes
0
answers
241
views
How to create a Grid in R for Kriging
I have done my work for creating variograms etc, and i would like to run kriging. I've loaded my dataset, and i know that i should create a grid file, but i don't know how.
I have the shapefile from ...
- 11
0
votes
1
answer
202
views
Does Kriging change known values
Does the process of Kriging routinely change known values to accommodate other nearby known values?
You can see I know little to nothing of this Kriging stuff. My concern is that a known and accurate ...
- 11
1
vote
1
answer
227
views
Kriging with gstat and Gaussian model produces bad outputs
I am trying to produce kriging models from elevations points using the gstat package. I can fit models to the empirical variogram using exponential, spherical, and Gaussian curves. However, despite ...
- 11
0
votes
1
answer
167
views
Remove X and Y axis from PyKrige result
I am trying to use the PyKrige module in Python to interpolate via Ordinary Kriging my XYZ dataset.
I have the below script, code taken from https://colab.research.google.com/drive/1mZDVJ-...
- 25
2
votes
1
answer
727
views
PyKrige constantly running into memory errors even with small dataset
I have an xyz file in which I'd like to apply the interpolation method of kriging to via python.
In order to do this, I have come across the module PyKrige and duplicated this example perfectly on my ...
- 25
0
votes
1
answer
2k
views
Using PyKrige for ordinary kriging on an XYZ file
I've previously used ArcGIS for Ordinary Kriging, but I now need a Python alternative.
I found the PyKrige module which i decided to give a try, but I am bit confused by it.
So far, I have read my XYZ ...
- 71
0
votes
1
answer
394
views
Kriging Interpolation in ArcMap
I want to interpolate one field of the shapefile named 'Interpolation' by Kriging method. But want to change the processing environment to the extent of another file namely 'Polygon'.
I tried '...
- 141
0
votes
1
answer
94
views
Using a dynamic formula with parrallel kringing in R
I am trying to achieve parallel Kriging in R on several variables using a loop. Here is a reproducible example using data meuse and a code for parallel kriging that I found here. In the code below, ...
- 631
0
votes
1
answer
571
views
Ordinary Kriging: "Warning message: In fit.variogram(data.var, vgm("Mat")) : No convergence after 200 iterations: try different initial values?"
I am trying to fit a variogram for kriging using the sp and gstat packages in r, but when I fit the variogram using a mattern fit, I get this warning message:
Warning message:
In fit.variogram(data....
- 101
0
votes
1
answer
185
views
Error message with Krige function: Error in as.character.default(X[[i]], ...) : no method for coercing this S4 class to a vector
I am trying to run the Krige function but I keep receiving NA values for var1.pred and var1.var. After using the warning() command, it says: Error in as.character.default(X[[i]], ...) :
no method for ...
- 1
1
vote
0
answers
23
views
Calculating volume of IDW raster below Height of Plane, but it is failing
I am using ArcMap.
I am trying to calculate the volume of water using IDW Tool. When I reach a point to calculate, the statistics below height of plane cannot run, giving me 0.00,
What could be the ...
0
votes
1
answer
603
views
Using ArcGIS geostatistical wizard for Kriging in ArcPy
In ArcMap , I have many interpolations to create and need a way of doing batch processing. I would like to use Universal Cokriging for interpolation. I have point features in a regular 0.25 arc deg ...
- 329
1
vote
1
answer
491
views
R - Why are my kriged interpolation values so low when using ordinary kriging (gstat)?
I am using a synthetic dataset that I generated for the purposes of fulfilling an assignment. I am attempting to use ordinary kriging to interpolate point data across an entire field, but my ...
- 101
1
vote
1
answer
475
views
Kriging with gstat package in R produces empty interpolation points
I have kriged a dataset with the gstat package in R, but it has produced an empty variable as a result. I have 1,167 measurement points within a field, and I am trying to interpolate them across 3,464 ...
- 503
0
votes
0
answers
2k
views
Type error when kriging with python: TypeError - Input z must be at least a (2, 2) shaped array, but has shape (1, 1)
I am new to kriging with python, and I have a csv file that I am trying to interpolate and plot. However, I am getting the following error:
TypeError - Input z must be at least a (2, 2) shaped array, ...
- 503
0
votes
1
answer
193
views
Adding kriging results back to a GWR model in R?
How can I add the kriging results of a downscaling process back to a GWR model? Basically, I have 2 raster, 1 at 400m (let's call it x1) resolution and 1 at 100m resolution. I upscaled the latter to ...
- 91
1
vote
1
answer
385
views
Geographically weighted area-to-area Regression CoKriging for raster downscaling in R?
this will be probable a series of questions, so in this one I will try to set the big picture. I want to perform geographically weighted area-to-area regression Cokriging (GWATARCoK) to downscale 1 ...
- 91
0
votes
1
answer
115
views
Residual bubble plot of krige.cv function dataframe
I want to plot a bubble plot of residuals from the krige.cv function with the bubble function.
However I get an error:
Error in as(obj, "SpatialPointsDataFrame") :
no method or default for ...
- 161
1
vote
1
answer
462
views
Export kriging plot as GeoTIFF
I have performed kriging and plotted it with ggplot, now I want to load the plot in QGIS. For that to work I need a GeoTIFF file of my plot. How do I convert my plot into a GeoTIFF file?
EDIT: added ...
- 161
1
vote
0
answers
277
views
Ordinary kriging in QGIS with geodetic coordinates
I am using QGIS 3.16.9 Hannover. I want to create a surface (with kriging) from a list of points whose SRC is WGS84, so I am working with geodetic coordinates. They are in decimal degrees covering a ...
1
vote
0
answers
123
views
Changing value while converting Geostattistical layer to raster using Geostatistical analyst
I using geostatistical analyst for interpolation. when geostatistical layer create and I want to export it to the raster the range of min-max values changes. how could I export GA layer to raster ...
2
votes
1
answer
155
views
Using ArcPy to iterate feature classes for kriging
I have 12 feature classes and they are named from net_0 to net_11 but the order of feature classes In Arc catalog is,List_wells = [net_0,net_1,net_10,net_11,net_2,net_3,net_4,net_5,..net_9].
I used ...
- 23
2
votes
1
answer
773
views
Fitting a variogram for Kriging in R
I want to fit a variogram of my data for kriging.
Here is what I have so far:
library(sp)
library(rgdal)
require(gstat)
setwd("C:/Users/49151/Desktop")
#import data
data <- read.csv(&...
- 161
1
vote
1
answer
1k
views
Plotting a heatmap of Kriging with coordinate in R
I performed Kriging in R. But I have trouble getting a plot of a colored heatmap showing the resulted interpolated area along with a lon/lat values at the side and on the button.
Here is an example, ...
- 161
1
vote
0
answers
117
views
Why can't gstat::krige accept duplicate locations when universal kriging?
Let us create a simple simulation worthy of gstat::krige with universal kriging, in the language of linear mixed models, i.e.:
Where X is a n x n_vars matrix of n_vars additional variables, betas is ...
- 111
0
votes
0
answers
505
views
Ordinary kriging is not working in QGIS 3.20 & SAGA-7.2.0
I am using below data for making kriging interpolation.
structure(list(Serial = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L,
10L, 11L, 12L, 13L, 14L, 17L, 18L, 24L, 28L, 30L, 31L, 36L, 37L,
43L, 44L, 46L, ...
- 151