Stack Overflow
Stack Overflow

Questions tagged [gstat]

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

130 questions
2
votes
0 answers

What is the best way to interpolate (IDW) Spatial Data into raster ? [R]

I am working in flood forecasting and analysis with an intern 2D model which gives particular output of water level, not very good looking and not usable with meshers. Today I use "SMS" a mesher from AQUAVEO which can interpolate data on the mesh.…
Maxime auffret
  • 91
  • 5
2
votes
1 answer

How is variogramST dealing with distance and time?

I'm new in spatio-temporal analysis in r so I don't know if what I'm doing is right. I have a lot of spatio-temporal data and I would like to plot a spatio-temporal sample variogram. My problem is that I don't really understand how variogramST…
AliceW
  • 21
  • 1
2
votes
0 answers

How do I ensure a 3D synthetic random field reproduces its variogram using R gstat

I have attempted to generate a 3D synthetic spatial random field using unconditional simulation in R gstat; with the lines below. However, my attempts to check if the exhaustive data can reproduce the variogram used in generating the field have been…
TTOWG
  • 31
  • 2
2
votes
2 answers

How slow is too slow when kriging with gstat in R

I am trying to use the krige function in the gstat package of R to interpolate some spatial ocean depth data in R. I am finding for more than about ~1000 points, the function starts taking unreasonable amounts of time to finish (i.e., hours to days…
user3004015
  • 617
  • 1
  • 7
  • 14
2
votes
1 answer

Whats the difference between z~1 and z~LON+LAT in plotting variogram in R?

library(sp) data(meuse) coordinates(meuse)<-~x+y v<-variogram(log(zinc)~1,meuse) v np dist gamma dir.hor dir.ver id 1 57 79.29244 0.1234479 0 0 var1 2 299 163.97367 0.2162185 0 0 var1 3 …
Orpheus
  • 329
  • 6
  • 21
2
votes
1 answer

Finding nearest grid value in spatial dataframe in R

I am learning about spatial dataframes. I have created a simple example that uses inverse distance weighting (IDW) to interpolate values across a grid, based on an initial set of points. I would now like to obtain the interpolated value for each of…
user1885116
  • 1,757
  • 4
  • 26
  • 39
2
votes
1 answer

How to convert a dataframe into a gstat object in R? Specifically a gstatVariogram

I'd like to create a variogram fit from some semivariance data in R. I realize that a variogram fit requires an object generated by vgm, and raw.vgm gneerated by variogram, for example: x = seq(0,300,by=300/100) y = seq(0,0.9,by=(0.9/100)) Z =…
makansij
  • 9,303
  • 37
  • 105
  • 183
2
votes
1 answer

Logistic function for regression kriging

I would like to perform regression kriging (RK) for binary presence-absence and host grid data as a constant predictor. I have used logistic function to estimate relationship between binary outcome and predictor, however I think it is not passing…
mihoo
  • 157
  • 3
  • 12
2
votes
1 answer

Create spatial objects in R useful for coordinates() and spsample()

I'm trying to use this code, adapted from dataset meuse data<-list(var1,var2,x,y) coordinates(data)=~x+y grid = spsample(data, type = "regular", cellsize = c(0.05,0.05)) vt <- variogram(var1 ~ var2,data=data) vt.fit <- fit.variogram(vt, vgm(0.2,…
Darko
  • 1,448
  • 4
  • 27
  • 44
2
votes
2 answers

Spatial interpolation in R gstat without the message

When I'm doing interpolation in R using gstat package, the message like this '[inverse distance weighted interpolation]' or this '[ordinary or weighted least squares prediction] ' occurs.…
Jot eN
  • 6,120
  • 4
  • 40
  • 59
2
votes
2 answers

Is it possible to add lines/text/minor tick marks/personalized axes to a variogram plot?

I'm trying to add text, lines, minor tick marks, and/or personalized axes to a variogram using the plot.variogram function. When I try to add any of these using the traditional function calls (i.e. text("Text Here")) it returns the error of plot.new…
Jared Smith
  • 280
  • 3
  • 12
2
votes
1 answer

Spatial temporal kriging

I'm working on rainfall data. I have grid of rainfall points of 1 square kilometer with time step of 10 minutes. I want interpolate those point at every time step. also I want to check the temporal correlation of my rainfall points. The program I'm…
Herve
  • 21
  • 1
2
votes
2 answers

Assign and use data values of Spatial*DataFrame easily

Is it possible to somehow easily work with data in Spatial*DataFrame (* = Points, Lines, Polygons, Pixels, Grid,...)? In particular I have difficulties with assigning values to it and to operate with…
Tomas
  • 57,621
  • 49
  • 238
  • 373
1
vote
0 answers

STplot package: How to slice y-coordinate when using mode 'xt' for spatio-temporal kriging?

I'm trying to interpolate my dataset using spatio-temporal kriging approach from gstat and spacetime packages. For the plotting of the results, I use STplot function. I want to use the mode 'xt' which produce the space/time plots. But when I read…
Abdul Azim
  • 11
  • 2
1
vote
0 answers

Different results between kriging on residuals and universal kriging (R)

I have a dataset with 132 observations, 1 response variable (Lopend_gemiddelde), and 9 predictors on which I perform three methods: 1) multiple linear regression (MLR), 2) kriging on the residuals (based on a multiple linear regression), and 3)…
Foeke Boersma
  • 95
  • 6
1
2
3
…
8 9