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
1 answer

Understanding gstat's variogram

I am new to using gstat and variograms and hoping someone can help explain some details to me, that don't appear to be in the documentation. I have used the default data and an example from one of the help pages to fit a spherical…
Sarah
  • 3,022
  • 1
  • 19
  • 40
0
votes
0 answers

How to create moving window filter of semivariogram outputs in R using focal area function?

I'm trying to create a raster filled with semivariogram outputs such as the sill and range of an area which describes the spatial autocorrelation. I wanted to try using the focal function in R as a way to scan an area, and the variogram function…
Melinda Martinez
  • 39
  • 1
0
votes
0 answers

How to extract values from Kriging interpolation

I can't extract ‘var1.pred’, ‘var1.var’ from the predicted values in kri, they are all shown as NA. data=read.csv('D:/R/R程序/2016.12.31/对比气温数据/2016年均温.csv') bound=readOGR('F:/ArcGis/RACGIS数据库/江苏省行政区/泰州.shp') dsp <- SpatialPoints(rev(data[,1:2]),…
wzp
  • 1
  • 2
0
votes
0 answers

Dynamic variable names in for loops combined with gdal

I have the following problem trying to somehow solve: I created 10 arrays of different IDW interpolations using the gstat/gdal package. Now I have the following variables, which I´m trying to export to my harddrive: DI.IDW.SAND.P.0.1,…
Microwave_Woodchuck
  • 13
  • 5
0
votes
1 answer

Pass dataframe function to gstat formula argument

I'm trying to create inverse distance weighted rasters using gstat() and raster::interpolate(). I'm running into issues passing a column name to the formula argument of the gstat function. Hardcoding the column name works just fine: gs <-…
philiporlando
  • 941
  • 4
  • 19
  • 31
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
0 answers

Variogram plot for Sill, Nugget, Range

I'm working on Temporal-spatio data. Data is in STFDF structure. I'm at the stage where applying variogram and try to determine Nugget, Sill, and Range from variogram plot that should look like this . something like this: However when I plot my…
Reta
  • 363
  • 3
  • 4
  • 15
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

Kriging simulation using ff package

I'm trying to understand the way I can use the ff package to overcome the error "Error: cannot allocate vector of size 1.1 Mb" while using kriging/ gaussian simulation. I don't know how to change the input data. Is there any idea to help me do…
Mohammad
  • 67
  • 1
  • 9
0
votes
1 answer

Extracting fitted values of a variogram model

How can I extract the fitted values (not parameters!) of a variogram model? I intend to draw the theoretical variogram in the Excel software. I used the fit.variogram in gstat package as follows: library(sp) data(meuse) coordinates(meuse) =…
Mohammad
  • 67
  • 1
  • 9
0
votes
1 answer

Defining new correlation model in gstat package in R?

Is it possible to define and use the new correlation model in vgm function of gstat package in R, in order to simulate random fields? I need to use correlation function (as suggested here ) for the log-transformed exponentially correlated spatial…
Saubhagya
  • 51
  • 1
  • 4
0
votes
1 answer

Temporal Variogram

I would like to build a temporal variogram, without specifying spatial coordinates as I have only one gauge that I'm studing temporally. Here is a sample of my data: > head(test) POSTE DATE RR6 4485 63113001 1358609040 0.0 4486…
I.Seck
  • 11
  • 5
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 identify the common coordinates in two different SpatialPixelDataFrames

I have two SpatialPixelDataFrames (a and b) which have few common coordinates/locations. I need to identify these common coordinates and remove them from both the dataframes. Though I have been able to identify and separate common coordinate points…
Saubhagya
  • 51
  • 1
  • 4
1 2 3
…
8 9