Stack Overflow
Stack Overflow

Questions tagged [kriging]

A statistical interpolation method, also known as Gaussian process regression, most used in geo-statistics. The goal is to map a surface given limited sample data. The process evaluates the variability of supplied data, then uses a weighted average of neighbouring points -- considering both distance and direction -- to interpolate the desired map points.

186 questions
0
votes
0 answers

Interpolating (kriging?) Grid Data in MATLAB

I have Lat, Lon, and NO2 data that are gridded. I'm trying to select a few points from this data and interpolated out some data/a plot, comparing how this interpolated data/plot improves as I interpolate it. My data is as follows: Lat is size…
Taylor
  • 89
  • 1
  • 11
0
votes
1 answer

Interpolating values of a vector with Gaussian Process / Kriging in Python

I'm attempting to interpolate wind speed and direction values on a map given some lat/long coordinates and then compare those values to my observed values. A couple papers suggest that Gaussian Process / Kriging are effective methods for this, but I…
markcoder1980
  • 13
  • 1
  • 3
0
votes
1 answer

Problem in exporting a kriging map as Raster

I created a kriging map which was created by using the kriging() and image() functions of the kriging package (table is the data with the coords and values): krig <- kriging(table@coords[ ,1], table@coords[ ,2], response = table@data$Joined.l_8,…
Stathis G.
  • 145
  • 2
  • 11
0
votes
1 answer

Problems with kriging function in R

I have a dataframe named as kef, consisted of 512 rows, and the fields x, y (referring to coordinates) and v (refering to a certain numeric value for each cell). I also have a map layer named as grecia.map, loaded in R through the readOGR command…
Stathis G.
  • 145
  • 2
  • 11
0
votes
1 answer

How to extract specific values with point coordinates from Kriging interpolations made in R?

By using R version 3.4.2 and the library "geoR", I made kriging interpolations for different variables (bellow I give an example of my process). I also made a matrix with the coordinates for 305 trees with distinct marks (species, DBH, Height) that…
Carlos SH
  • 1
  • 2
0
votes
1 answer

Kriging at a single spatial point?

This question may have to do with my poor knowledge of kriging - is it possible to compute kriged value at a single spatial location? As I understand it, a typical kriging method uses the spatial correlation embedded in a sparse data to interpolate…
Manojit
  • 611
  • 1
  • 8
  • 18
0
votes
1 answer

R_Co-kriging when the variable of interest and auxiliary variable(s) are not measured at the same locations

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 not measured at the same locations. Data frame with…
rm167
  • 1,185
  • 2
  • 10
  • 26
0
votes
1 answer

geoR kriging error: solve.default

I am trying to do kriging in geoR for a fairly large area (~1 million km^2). It is for my thesis, so unfortunately I cannot share the data. I have already checked for duplicated in the coordinates and the actual data to see if this is why there is…
mcgill
  • 1
  • 1
0
votes
0 answers

Is Monte Carlo sampling in gstat: sequential Gaussian Simulation duplicitous?

How can I control over sampling in sequential Gaussian simulation? For example in the following code, how can I guarantee that the Monte Carlo samples are not duplicitous? library(sp) library(gstat) data(meuse) coordinates(meuse) =…
Mohammad
  • 67
  • 1
  • 9
0
votes
1 answer

why does the kriging give the same values as the observed?

I did kriging using spPredict from spBayes package for Bayesian kriging and krige from gstat package for non-Bayesian kriging. I didn't use any covariates (only constant mean term) and used 1283 points for kriging. Out of 1283 points, 1012 are new…
Donna
  • 11
  • 1
  • 5
0
votes
1 answer

R 'automap' how to create prediction grid to use with AutoKrige (e.g. meuse.grid)?

I'm having a lot of difficulty creating a prediction grid (for the new_data argument) to use with the autoKrige function in the automap package. I've already tried following the steps in this post (How to subset SpatialGrid using SpatialPolygon)…
Tom Newton
  • 199
  • 2
  • 14
0
votes
0 answers

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

I would like to use ordinary kriging get model predictions on a grid. In my code underneath, GPS is a 'SpatialPointsDataFrame' and EU is a 'SpatialPolygonDataFrame'. All covariates for estimating the model are contained in 'All_data'. These…
Esmee
  • 21
  • 2
0
votes
1 answer

How to write a formula in a loop for autoKrige {automap} in R

I have a large matrix and the column names are as follows: colid=vector(length = 60) for(i in 1"60) { colid[i]=paste0("V",i) } When I use the autoKrige function in automap, a formula must be written in such format: V1~1. When I try to do this…
Yang Yang
  • 858
  • 3
  • 26
  • 49
0
votes
0 answers

Perform Kriging Interpolation using Arcpy

I have list of point feature class. I am trying to write a python script to perform Krigging interpolation. I am getting error massage in this code "Point_Num" is not defined, Below script i am working import arcpy from arcpy import env from…
Pruthvi
  • 1
  • 1
0
votes
0 answers

What is the difference between the R code for ordinary kriging and block kriging used gstat package and the krige() function?

i have some problems with block kriging. First I create a grid (5000*5000 m). Than i exclude stations, which are very close to each other. The next step is creating the variogram and the fit.variogram. I assign the projection and the coordinates.…
Kira
  • 1
  • 1
1 2 3
…
12 13