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
1
vote
1answer
35 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
1answer
36 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
0answers
34 views

How would a Kriging model react to different z values at same point

Say I have a very large dataset of values collected at points over many years. For the sake of an example, we'll pretend it's the altitude at certain lat/longs. If some of those points appeared twice ...
0
votes
0answers
17 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
0answers
16 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 ...
-2
votes
1answer
30 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, ...
0
votes
0answers
19 views

Issues fitting Variogram Model in R

I have taken a shapefile from Open NYC Data and performed the following method. My end goal is to predict Taxi trip_duration at various points across the city of NYC in March at given time segments. I ...
0
votes
0answers
58 views

How to perform regression Kriging with R studio using gstat?

I'm struggling with regression kriging in R. I want to interpolate the temperature of a watershed under consideration of elevation, but get always the same error when I want to perform the last step, ...
1
vote
0answers
27 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
0answers
4 views

Is it possible to access fields from a relationship class to use as a kriging value

I have a point dataset of well locations related to a table of calculated production figures by year through a simple relationship class. I am trying to use the production figures as the z value for ...
0
votes
0answers
21 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 ...
0
votes
1answer
62 views

0…0.kriging error messages, grid file problem ? (automap, gstat)

I am trying to conduct kriging of variable called percent in South Africa. head(LL12) Location Longitude Latitude percent 1 L1 16.48476 -28.60737 1.0000000 2 L3 16....
2
votes
0answers
103 views

Spatial interpolation of climate data using latitude, longitude and elevation in r

I have downloaded E-OBS gridded daily climate data on 0.1-degree grid. For each grid point, I have longitude, latitude, elevation and climate data, i.e. precipitation sum / mean temperature. Now, I ...
4
votes
0answers
65 views

How does Empirical Bayesian Kriging of ArcGIS Pro handle sampling bias?

In ArcGIS Pro, simple kriging approaches allow you to de-cluster your data before applying a normal-score transformation. According to the Esri guideline, normal-score transformations are only ...
0
votes
0answers
24 views

Kriging Interpolation using ArcMap

I am using ArcMap. I am using Kriging interpolation in order to show the heat-map on my study area. I have total of 6 variables which ranges between 0.119~0.300 when I select kriging from arc toolbox ...
0
votes
0answers
101 views

Interpolating SpatialPointsDataFrame to RasterLayer using autoKrige in R

I have a SpatialPointsDataFrame with elevation data and 51'257 elements. I also have a RasterLayer object of dimensions c(114,210) covering roughly the same area (all the spatial points are within the ...
0
votes
1answer
78 views

R function krige.cv() returns NAs

I recently completed a geostatistics course, and now trying to apply my knewly acquired knowledge to a real case. I have a dataset with 2 122 groundwater level measurement points. I would like to ...
3
votes
1answer
195 views

Interpolating DEM using Kriging in R

I have several DEMs as geotiff files which I load into R using raster::raster. Now some of these DEMs contain patches of missing values, represented by NA. I would like to interpolate these missing ...
0
votes
0answers
35 views

Spatio-temporal kriging

I have a geospatial dataset of air quality and I want to interpolate the air quality value across time and space for unmeasured location and time instace. I found gstat package in R could do the ...
1
vote
0answers
50 views

Interpolation with barriers in Google Earth Engine [closed]

Does anyone have a code example or helpful information to conduct Empirical Bayesian kriging on Google Earth Engine or of an interpolation with barriers technique in Google Earth Engine?
0
votes
1answer
33 views

Kriging for Soil Nutrient Analysis

I have collected soil analysis results on a gridded basis within fields. These grids were generated on a 35m based upon digitised polygons of the field extents. I have imported my soil analysis ...
0
votes
1answer
111 views

Removing raster values using a threshold in QGIS and interpolating

I have an imagery (plant cell density) comprising vine rows and grass inter-rows.I'm trying to permanently remove all the values of inter-row grass through setting up a threshold value. The vines ...
1
vote
1answer
286 views

Interpolation for unevenly spaced data points

I am currently having some trouble finding a good interpolation technique for unevenly distributed spatial data points. Let's give a quick example : I have a 1000 points dataset. As several surveys ...
0
votes
0answers
29 views

How to add additional GRASS modules to QGIS 3.2?

I like to use the v.kriging tool from GRASS (https://grass.osgeo.org/grass7/manuals/addons/v.kriging.html) in QGIS 3.2 but its not in the processing toolbox. Is this possible to import this tool and ...
1
vote
1answer
390 views

QGIS SAGA “Kriging” tools doesn't work

I try out different methods of point interpolation for a digital terrain model. For this purpose I tried to run the "ordinary kriging tool" from the SAGA Tools with QGIS 2.18 and QGIS 3.2. When I run ...
-1
votes
1answer
147 views

Seeking open source tool for 3D Interpolation

Empirical Bayesian Kriging 3D (EBK3D) looks to be the tool I need within ArcGIS. It requires a GeoStatistical Analyst license. Is there is a similar tool or library that is free and open source? I'...
0
votes
1answer
147 views

Latin Hypercube Sampling (LHS) in R

I am new to R and would like to generate Latin Hypercube samples (LHS) for three variables within the following limts: 1.25 ≤ x1 ≤ 5 , 1.25 ≤ x2 ≤ 5 7 ≤ x1 ≤ 25 with total 3*10 = 30 samples. What ...
0
votes
0answers
26 views

Joining z-values (from point file) to map generated by Kriging in ArcGIS Pro

I am using ArcGIS Pro 2.4.0. I am making a fertilizer application map that changes the rates, based on the location of the field. The overall goal of the operation is to create a grid of squares, ...
1
vote
1answer
25 views

Defining one numeric variable by defining another string variable (via parameter selection)

I'm creating an ArcGIS tool that will allow users to figure out how much nitrogen that a field needs based on interpolated points (Kriging). I will have the following parameters: #numeric variables ...
1
vote
1answer
76 views

Kriging error when writing ArcPy script (ERROR 000622 & ERROR 000628)

I am trying to write an ArcPy script that uses Kriging to interpolate various nutrients in the soil (i.e. nitrogen) using ArcGIS Pro 2.4.0. These nutrients were collected by sampling the soil, so it ...
3
votes
1answer
166 views

Variogram fit won't converge

I am having some difficulty fitting a variogram to the data provided at the end of this post. For context, this is baseball data where the x and z coordinates are for the location of the pitch and the ...
2
votes
0answers
167 views

Error while performing kriging interpolation in QGIS 3.8.3

I am using QGIS 3.8.3 -Zanzibar in order to perform kriging interpolation. I get this error massage: "algorithm 'Ordinary kriging' starting… Input parameters: { 'BLOCK' : True, 'DBLOCK' : 100, 'FIELD'...
1
vote
1answer
107 views

Kriging in QGIS 3.4 Problem with Variogram Model

In QGIS 3.4.13 - Madeira there seems to be a bug in the Ordinary Kriging Modul (Processing Toolbox > SAGA). When I enter the modell "1600000*x+19202*x^2" QGIS will add a second "^" leading to wrong ...
1
vote
0answers
555 views

Building Semivariogram/Variogram in QGIS

I am trying to Interpolate Point Data in QGIS using the Ordinary Kriging Modul (SAGA) from the processing menu. Is it possible to build a Semivariogram in QGIS or call the SAGA Variogramm function ...
2
votes
1answer
97 views

Define Range when calling SAGA Kriging from QGIS

Is there a possibility to define the range-value or display the Variogram when you call the SAGA-Kriging-Module from QGIS? Maybe by using the python console?
3
votes
0answers
49 views

Performing meaningful kriging in R (averaging nearby points) [closed]

I have a set of data on the quantitative composition of the particle size distribution of the sedimentary layer. I have to build a map of the distribution of certain fractions across the entire study ...
1
vote
0answers
58 views

Raster distortion from parallel kriging in R

I am relatively new to parallel computing in R. I am using a high performance cluster to perform kriging for soil data points. I found this previous post which I used as a template: How to achieve ...
0
votes
0answers
23 views

Kriging along perpendicular transects [duplicate]

I collected data along two perpendicular transects (in a "+" shape) and would like to create a prediction surface map with kriging. What form of kriging (method, parameters, etc.) is appropriate to ...
0
votes
1answer
128 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.?
2
votes
1answer
506 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
54 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
379 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
48 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 ...
1
vote
0answers
72 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
77 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
562 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
60 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
450 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:
1
vote
1answer
780 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
49 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
2 3 4 5