Stack Overflow
Stack Overflow

Questions tagged [gstat]

gstat is an R package for spatial and spatio-temporal geostatistical modelling, prediction and simulation

130 questions
0
votes
0 answers

Why does gstat.predict() function often return NaN values (GSTAT Package)? (R version 3.3.2, Windows 10)

I am trying to simulate a combination of two different random fields (yy1 and yy2 with different mean and correlation length) with an irregular boundary using Gstat package in R. I have attached the picture of my expected outcome. The code is not…
Saubhagya
  • 51
  • 1
  • 4
0
votes
0 answers

Confidence Intervals Variograms

I am currently using the gstat package to create variograms of oil and gas production in counties in West Virginia. Here is the link to my dataset of SpatialPointsDataFrame and the code to produce a variogram is…
user5425676
  • 1
  • 2
0
votes
1 answer

Plot a list of lines with R lattice package

I'm trying to show some variogram model fits on the top of xy.plots. The panel.plots work for all except for a list of lines I wanted to add to corresponding subplots. require(gstat) require(sp) data(meuse) names(meuse) #make directional…
ToNoY
  • 1,358
  • 2
  • 22
  • 43
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
0
votes
0 answers

Why the RMSE for Kriging is not strictly decreasing for increasing training data?

I am using gstat package for ordinary kriging and using the walker lake data (data size = 470). I have randomly taken 20 from that data in each trial and calculate the rmse for randomly chosen training dataset from 50-450 dataset. Then I have…
Chandan
  • 764
  • 2
  • 8
  • 21
0
votes
1 answer

IDW parameters in R

I want to perform IDW interpolation using R using the idw command from the gstat package. I have this data: #settings library(gstat) library(dplyr) library(sp) library(tidyr) id_rep <- rep(c(1,2), 20) f <- rep(c(930,930.2), each=20) perc <-…
Lince202
  • 143
  • 10
0
votes
1 answer

R gstat spatio-temporal variogram kriging

I am trying to use the function variogramST from the R package gstat to calculate a spatio-temporal variogram. There are 12 years of data with 20'000 data points at irregular points in space and time (no full grid or partial grid). I have to use the…
user823
  • 265
  • 2
  • 14
0
votes
1 answer

Adding an argument to a mathematical function

I'm trying to automate a process to produce interpolated geochemical maps. I've created a loop that essentially starts at column #13 and loops through until #67. However, part of the IDW code requires the header of the current column for the…
M. Aberle
  • 25
  • 1
  • 8
0
votes
0 answers

Create a gridded dataset in sp and gstat - R

I have trouble getting my data to create a grid using sp and gstat package in R. I have a relatively small data set with GPS coordinates and I want to interpolate my data and make a gridded map. My data is in a csv file and is currently not in a…
Karen Lykkebo
  • 1
0
votes
1 answer

Automap package: Negative values and variogram fitting

I started using the automap package for automated variogram modelling two days ago and played around with it (doing it manually was tedious and doesn't come with too many advantages...although I find it strange that automap seems really to love the…
GeoEki
  • 437
  • 1
  • 7
  • 20
0
votes
0 answers

How to find DEM elevation value at sill in semi-variogram

As semi-variograms gives semi-variance and lag distance for DEM. I want to know at what DEM elevation the sill value is flattening using penta spherical model.
sreedhar
  • 1
0
votes
0 answers

Local kriging reports error in `krigeST` with full data (`STFDF`)

I have been playing around with the krigeST function in the gstat package, and I tried to do local kriging with a full set of SpatialPolygons based data (STFDF). When I tried the following spatio_time_krige_sumProd = krigeST(conso~1, …
vermouth
  • 301
  • 2
  • 12
0
votes
1 answer

R: Bad variogram fitting, bad kriging results

I try to do a kriging in the Jakarta Bay. I have a set of measurement points with appropriated coordinates and attributes (pH, salinity,...) In order to do a kriging I first need to find a model for my variogram. When I use the "variogram" function…
ZKB
  • 101
  • 2
0
votes
1 answer

Customize inverse distance predictions within raster classes in R

I want to modify the following code so that for all values beyond the defined maximum distance I can use the average values of points which fall in the same soil order category. Any recommendations? # load packages and…
Geo-sp
  • 1,704
  • 3
  • 19
  • 42
0
votes
1 answer

R Gstat compute semivariogram on a grid

I have generated a 3D unconditional normal scores simulation using gstat package and would now like to check the variogram reproduction. How do a specify the calculation of the variogram on the grid? #-- Create the…
Markm0705
  • 1,340
  • 1
  • 13
  • 31
1 2 3
…
8
9