Stack Overflow
Stack Overflow

Questions tagged [geostatistics]

32 questions
6
votes
1 answer

Multi-output spatial statistics with gaussian processes

I've been investigating Gaussian processes lately. The perspective of probabilistic multi-output is promising in my field. In particular, spatial statistics. But I encountered three problems: multi-ouput overfitting and anisotropy. Let me run a…
essicolo
  • 743
  • 5
  • 12
3
votes
1 answer

how to change the lag distance while calculating a Variogram in R

I am trying to calculate an Experimental-Variogram value at different lag distances, so I am using the variogramm command variog1 <- variogram((Copper)~1,ds) but I can't know how to specify the needed lag distance. For example I want to get a…
Ahmed Emam
  • 33
  • 6
3
votes
2 answers

Underlay a vector image to a grid for kriging in R

After searching around a lot, asking, and doing some code, I kinda got the bare minimum for doing kriging in R's gstat. Using 4 points (I know, totally bad), I kriged the unsampled points located between them. But in actuality, I don't need all of…
ace_01S
  • 327
  • 2
  • 5
  • 12
2
votes
1 answer

GAM with "gp" smoother: how to retrieve the variogram parameters?

I am using the following geoadditive model library(gamair) library(mgcv) data(mack) mack$log.net.area <- log(mack$net.area) gm2 <- gam(egg.count ~ s(lon,lat,bs="gp",k=100,m=c(2,10,1)) + s(I(b.depth^.5)) + …
user3036416
  • 1,025
  • 2
  • 12
  • 25
1
vote
0 answers

Line segments coming out of the interior of a map of brazil using the ggrepel package

I am trying to reproduce the format of the map below using the geobr package, however, I have been facing some problems to insert these" lines "on the map. After a search, I found the ggrepel package that brings with it a series of functions capable…
Breno S.
  • 271
  • 1
  • 8
1
vote
1 answer

Is there an issue with how I'm reading in my shapefile and plotting it?

I've tried to read in my shape file and then plot it but it seems as though RStudio is stuck on something and will not finish running the plot function. Right now I have: library(rgdal) new_county_path <- paste(county_path,…
Buchlord
  • 11
  • 1
1
vote
1 answer

PyKrige witing a file with 3D kriged data

I am trying to krig soil volumes in 3D using pykrige. I am using ordinary kriging, but once I am done kriging the data, I am not sure how I can write it to an ascii file, Here is my code: from pykrige.ok3d import OrdinaryKriging3D from pykrige.uk3d…
Astro
  • 11
  • 1
1
vote
2 answers

Alternative to for-loop for large dataset to improve computational speed

For kriging, I need to compute large mesh arrays of length 20000. The code below works fine, especially for small mesh length (< 100), however, the computational time for such large mesh is very long (approx 45min). The length of data ranges between…
user2554925
  • 457
  • 2
  • 7
1
vote
0 answers

Nugget value is always zero while using scikit-geostat

import numpy as np import skgstat as skg coordinates = long_data values =N_data V = skg.Variogram(coordinates=coordinates, values=values ,maxlag=np.max(dists)/2.0 , use_nugget =True ) print(V) V.plot() enter image description here from the graph…
Nikhil Munna
  • 113
  • 1
  • 6
1
vote
0 answers

How do I fit directional variogram on gridded data, and for a given azimuth return the estimate range?

I have a digital elevation model of topographical features. My overall objective is to make a predictive model by regressing on the range of a directional variogram for a given azimuth. Are there any good python libraries which can accommodate this…
F.Dunbar
  • 37
  • 3
1
vote
0 answers

Get the lag vector from variogram in gstat

I want to compute the variogram from a set of data in R. I am using the function "variogram" from the gstat package. Now, I want to get the lag vector from the variogram. The problem is that myvariogram$dist returns the averages of the distances…
lgndrzzz
  • 274
  • 2
  • 12
1
vote
0 answers

Regression Kriging of binomial data in geoRglm R package

I am using binom.krige() function of the R package geoRglm for determining the spatial predictions of a binary (0, 1) response variable with several continuous as well as discrete covariates. Using glm() with binomial logit link function I found…
Asad Ali
  • 25
  • 6
1
vote
1 answer

How to merge plots (layers) in R Studio?

How can I merge two plots (layers)? On a first plot there are points with data, and on the second one there is a boundary of those points. Here is what i've got: data <- read.csv('data.csv') coordinates(data) <- ~x+y proj4string(data) <-…
Hubert
  • 27
  • 8
1
vote
1 answer

R: How to or should I drop an insignificant orthogonal polynomial basis in a linear model?

I have soil moisture data with x-, y- and z-coordinates like this: gue <- structure(list(x = c(311939.1507, 311935.4607, 311924.7316, 311959.553, 311973.5368, 311953.3743, 311957.9409, 311948.3151, 311946.7169, …
mattu
  • 156
  • 1
  • 16
1
vote
0 answers

R: Help needed to understand RMtrend (RandomFields)

When working with the R package RandomFields you specify your model describing the spatial distribution of your data by adding different components. E.g: Nearer points in your data are more familiar to each other than farer points are, so that they…
mattu
  • 156
  • 1
  • 16
1
2 3