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
1 answer
20 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 ...
0 votes
0 answers
33 views

PyKrige interpolation results very different to ArcGIS Pro

I had been using ArcGIS Pro's Kriging tool to apply Ordinary Kriging to my dataset. This produced the below results, which were exactly as I wanted (I have also included a screenshot of the parameters ...
0 votes
1 answer
11 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-...
1 vote
1 answer
37 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 ...
-1 votes
0 answers
30 views

Seeking Kriging Interpolation tool for Python

I have been using ESRI ArcGIS Pro's 3D Analyst Geoprocessing Tool Kriging to interpolate my data. More specifically, I have been using the Ordinary Kriging with a cell size of 0.5 on 36935 xyz points ...
0 votes
1 answer
43 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 ...
0 votes
1 answer
39 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 '...
  • 37
0 votes
0 answers
8 views

Change Kriging Output Scales using ArcGIS Pro

[![enter image description here][1]][1] I have made a series of Kriging maps based on different time periods, but I want them all to have the same scale. For example, all should be 33-69.5. Is this ...
  • 1
0 votes
0 answers
23 views

IDW vs. ordinary Kriging, spatial clusters

Using ArcMap, I have a dataset containing points with coordinates and value (cost of infrastructure). There is no even distribution, the points are mostly clustered around larger cities but with some ...
0 votes
0 answers
8 views

Number of Pairs in lag class of variogram using ArcGIS Geoststistical Analyst [duplicate]

My question has already been asked (but remains unanswered) as ArcGIS Geostatistical Analyst variogram lag pairs I'm looking to find out how many Point pairs are in a lag class in the semivariogram ...
0 votes
1 answer
55 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, ...
  • 561
0 votes
0 answers
11 views

Kriging output is blank

I have used kriging with default settings where I have five points with XY coordinates, ID and concentration. The concentration is ranging from 12 to 161900. The output I get from this is just a ...
0 votes
1 answer
67 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....
0 votes
1 answer
45 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
13 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
143 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 ...
  • 197
0 votes
0 answers
414 views

Cokriging Predict Function "Error in h(simpleError(msg, call))"

I tried to do CoKriging using gstat library in R. it works perfectly until the prediction function CK <- predict(g, Grid) It gave me this error: Error in h(simpleError(msg, call)) : error in ...
1 vote
1 answer
133 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 ...
1 vote
1 answer
78 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 ...
  • 463
0 votes
0 answers
72 views

QGIS Ordinary Kringing Python error

I am new to using the Smart-Map plugin. I added my data (which has one valid outlier) and did the ordinary Kriging interpolation. Just before it reaches 100% I get the following error message: 2022-03-...
  • 31
0 votes
0 answers
474 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, ...
  • 463
0 votes
0 answers
53 views

ArcGIS kriging unable to estimate semivariogram when processing extent is set

I'm trying to use ArcGIS kriging tool to interpolate some points to form a raster. The point covers approximately a polygon study are but not exactly the extent of the area. I set the processing ...
0 votes
1 answer
76 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
0 votes
0 answers
20 views

Appropriate kriging variogram model and parameter values in order not to loose information?

I have tried to apply kriging in order to interpolate the measurements to a regular grid. The function that generates the grid is the following: from pykrige.ok import OrdinaryKriging def ...
1 vote
1 answer
191 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
45 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 ...
0 votes
1 answer
133 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 ...
1 vote
0 answers
123 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
63 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 ...
0 votes
0 answers
71 views

Stein's parameterization matern model kriging SAGA in QGIS

I performed autovariogram in R. It gave me the output seen below. It suggest using a matern model-Steins parameterization (indicated by the "Ste"). I would like to try the kriging tools by ...
0 votes
0 answers
92 views

What do the kriging-parameters in SAGA QGIS mean?

I'm used to standard kriging-methods in an R environment. I'd like to know what each parameters in the kriging-tool from SAGA in QGIS mean. Parameters like: type of qualityindex Blockkrigen/Block ...
2 votes
1 answer
58 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
1 vote
1 answer
121 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(&...
1 vote
1 answer
430 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, ...
1 vote
0 answers
49 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 ...
0 votes
0 answers
196 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, ...
  • 91
1 vote
0 answers
382 views

Simple kriging in raster

I want to use a Python script for interpolation in 1st image to have the result of 2nd one, this result is from simple kriging saga. I tried the code below but I have two main problems. Need a lot of ...
0 votes
1 answer
41 views

Geo data normal assumption - effects

I have a basic question regarding Geo data: is the normal assumption mandatory? Given a dataset I identified outliers and inliers (using Moran local and plot test) so I excluded them from the original ...
  • 121
1 vote
0 answers
61 views

NAN's in the kriging predictions in R

I am trying to perform the kriging prediction of my data. I follow a pretty standard process. Firstly, I prep a grid to which the kriging predictions are supposed to be inserted. I prep data that I ...
  • 111
0 votes
0 answers
24 views

Warning message while kriging in R [duplicate]

krig_Gaussian<-krige(formula=(B)~1, locations=Amb_m, newdata=gd_m, model=model_Gaussian.variog_m) [using ordinary kriging] Warning messages: 1: In proj4string(d$data) : CRS object has comment, ...
  • 1
1 vote
0 answers
41 views

Kriging observations using additional knowledge from geophysical model

Here is a simplified explanation of the problem I am trying to solve : I have (1) a 100x100 grid with air pollutant concentrations in a city obtained from an air quality geophysical model, and (2) a ...
1 vote
1 answer
497 views

QGIS: Spatial extrapolation

I have few sample points with values on (agricultural) field and I want to interpolate and extrapolate values inside the field boundary. Ideally I would like to use Kriging (as it smooths local extrem ...
  • 355
1 vote
0 answers
77 views

Upper and lower bounds from gstat::krige() and gstat::idw()

Is it possible to obtain lower and upper bound values from a kriging estimation using gstat::krige() in the same way it can be obtained from a simple lm() model? The following lm() example comes from ...
3 votes
1 answer
598 views

How to create a grid to do kriging interpolation that correspond to areas in a map having actual observations?

I need to do kriging interpolation but wanted to have the interpolation in areas where I have actual data; not across the whole map. This is my example. Here is the base map (shp file in a .zip file): ...
2 votes
1 answer
305 views

Why krige function return all NA's values for the prediction when trying to interpolate a variable?

I am trying to build a map with kriging interpolation for a variable. When running the krige function the return is all NA values and the kriged map plot only one value. I couldn't find out why and I ...
0 votes
1 answer
438 views

Rule of thumb for complexity of 3d Kriging (with pykrige)

I've been using pykrige for some data interpolation. It's great and works well in 2d. I recently tried the 3D version though, and it takes forever to run (> a day before I killed it), just for ...
1 vote
0 answers
66 views

How to do a Correlation of various Raster Data and show the result in a Correlation Matrix

I am doing spatial correlation of soil Data sets. I did various Kiriging Interpolations with the krige function. The output was cut to a target raster to have the same raster size for all Kriging ...
  • 11
-2 votes
1 answer
106 views

Extracting time series trend analysis parameters from kriging using ArcGIS (versus R)

A grad student has asked me a question about extracting trend parameters for a time series of data using ordinary kriging ArcGIS tools. I'm not an expert in the spatial statistics options of ArcGIS, ...
  • 19
1 vote
0 answers
117 views

Forcing boundary conditions in interpolated/kriged data?

I am using pykrige.ok import OrdinaryKriging to generate a spatially continuous field from sparse point data (see image below). The data (glacier ice thickness) are such that there is a hard boundary (...
0 votes
0 answers
52 views

Kriging (or other interpolation) with integral constraints?

I would like to interpolate a continuous variable (from observed point locations) on a surface (gridded) but using an external constraint that the cumulative distribution function of the interpolated ...

1
2 3 4 5 6