Questions tagged [gstat]

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

130 questions
1
vote
1 answer

Create variogram in R's gstat package

Suppose I have rainfall data taken at four weather stations over the span of 2004-2016. I fed the data into a database for retrieval in R. My goal is to take the data for every single day from that period, and krige using those values,…
ace_01S
  • 387
  • 2
  • 5
  • 16
1
vote
1 answer

length erros in gstat using cross validation krigecv in R

Today I am having a problem which I could not find a solution to. I want to use krigcv cross validation using gstat in R and i have an error message: 0%Error in model.frame.default(terms(formula), as(data, "data.frame"), na.action = na.fail) :…
Franck
  • 11
  • 2
1
vote
1 answer

Local Block Kriging with Local Variogram with gstat

I have been unable to find any information specific to local block kriging with a local variogram using the gstat package in R. There is freeware called VESPER from the Australian Center for Precision Agriculture that is able to do this, and from…
1
vote
1 answer

R: Autokrige.cv function in automap package generates NaNs

I’m fairly new to R and I am trying to make interpolations of temperature measurements that where gathered from different station across the Netherlands. I have data for about 35 stations that make measurements every 10 minutes covering a timespan…
T.Merkus
  • 11
  • 2
1
vote
3 answers

Duplicate data when using gstat or automap package in R

I am trying to using ordinary kriging to spatially predict data where an animal will occur based on predictor variables using the gstat or automap package in R. I have many (over 100) duplicate coordinate points, which I cannot throw out since those…
ms.elasmo
  • 13
  • 4
1
vote
0 answers

How do I fit a gstat model in regression kriging?

I am trying to fit a model using regression kriging using r gstat package. I have loaded some disease data called prevalence and some covariate rainfall. The rainfall raster data is loaded as a .grd file.The code that I am using to fit this model is…
Paulo
  • 41
  • 1
  • 5
1
vote
1 answer

How to plot multiple semi-variogram from a single dataset efficiently in R?

I have a dataframe named seoul1to7 contains the hourly PM10 concentration data from 1 march to 7 march,2012.please,download.In this dataset, time is in yyyymmddhr format.for example, 2012030101 means 1 march 2012, 1.00 a.m. Data is Look like: …
Orpheus
  • 329
  • 6
  • 21
1
vote
1 answer

Problems is projecting LON/LAT data in R

I have a dataframe contains PM10 concentration in air over seoul(capital city) in Korea. Please, take a look. I want to plot semivariogram from this data set. As LAT/LON data here, is in degree so I have project this data. I have projected data in…
Orpheus
  • 329
  • 6
  • 21
1
vote
2 answers

How to fit model with semivariogram using gstat in R?

I have a csv file file contains atmospheric PM10 concentration data of 1 march,12.00 pm.please, download. I want to draw a semivariogram using gstat package in R. I tried to write these code in R. but with these data, I cant fit the model. …
Orpheus
  • 329
  • 6
  • 21
1
vote
1 answer

R: Gstat universal cokriging resolution

I am trying to do universal cokriging in R with the Gstat package. I have a script that i was helped with, but now i'm stuck and can't ask assistance from the original source. The problem is that i can't change the output resolution of the cokriged…
Castle
  • 53
  • 6
1
vote
1 answer

can't install package gstat in R version 3.1.2

I am trying to install the gstat package on the latest version of R (3.1.2) on windows. It is unsuccessful, even when trying to change the 64-bit version to the 32-bit version. Here is the message I get: when installing the package, it says "package…
aphili01
  • 11
  • 2
1
vote
1 answer

R Error in chol.default(A) krigeST from gstat package

I am working with an hourly dataset of air temperature, recorded at ~200 stations over a relatively small area. I chose a space-time variogram (e.g. sum-metric) to fit my data and am now trying to make predictions over my same stations in order to…
1
vote
1 answer

How to ensure the value of spatial interpolation more than zero

I am using the automap package to calculate spatial interpolation because I do not want to specify the nugget, sill, or range in the krige function in gstat package. The value I want to do spatial interpolation is PM 2.5. Its value can not be…
Ben
  • 665
  • 1
  • 10
  • 27
1
vote
1 answer

r heatmap on ggmap from kriging data

I want to plot a heatmap on a ggmap. library(ggmap) turku<-get_map('turku', zoom=13) turkumap<-ggmap(turku, extent="device", legend="topleft") turkumap turkumap+geom_density2d(mapping=aes(x = lon, y = lat),data = test, ) We have made a…
Irene
  • 744
  • 1
  • 12
  • 36
1
vote
1 answer

idw() or krige() Error: dimensions do not match when missing values

Functions idw() and krige() from gstat package keep reporting errors when either response or predictor variable contains missing values (NA), even when na.action is set to na.omit: require(gstat) data(meuse) coordinates(meuse) =…
Tomas
  • 57,621
  • 49
  • 238
  • 373
1 2 3
8 9