Stack Overflow
Stack Overflow

Questions tagged [geostatistics]

43 questions
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
  • 51
  • 9
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
  • 318
  • 2
  • 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
  • 318
  • 2
  • 16
1
vote
1 answer

R_number of pairs for each lag in a Variogram

I am using geoR package for spatial interpolation of rainfall. I have to tell that I am quite new to geostatistics. Thanks to some video tutorials in youtube, I understood (well, I think so) the theory behind variogram. As per my understanding, the…
rm167
  • 1,185
  • 2
  • 10
  • 26
0
votes
0 answers

Conditioning with External Drift Kriging using GSTools (and Pykrige...?)

I want to impute missing data on abortion rates of US-Counties. I´m comparing different Methods, one of them is Kriging. I used Pykrige for External Dirft Kriging: def Kriging_with_SGS(df, variable, kriged_var, variogram, lags, drifts=None): …
Lorenz
  • 1
  • 3
0
votes
0 answers

clustering using k-means/ k-means++, for data with geolocation

I need to define spatial domains over various types of data collected in my field of study. Each collection is performed at a georeferenced point. So I need to define the spatial domains through clustering. And generate a map with the domains…
Ferretto
  • 1
  • 1
0
votes
0 answers

software R: Error in (1 - h) * qs[i] : non-numeric argument to binary operator

I have three parameters with 899 points each, but one of them shows this message when I try to apply any command. Even the simplest command appears like this, for example: >basicStats(Values$Values.of.Th) Error in (1 - h) * qs[i] : non-numeric…
Victória
  • 1
0
votes
0 answers

Get URL of latest ONS Postcode/UPRN Lookup files

We'd like to use the UPRN and Postcode lookup files that the UK's Office for National Statistics (ONS) provide through their web portal: Postcode:…
Om95
  • 23
  • 4
0
votes
0 answers

How can i access to semivariogram in python with my information through scikit-gstat library?

I want to generate a semivariogram (exponential, gaussian, spherical) in python with regard to excel's information which has X, Y (need to calculate h), and EC in those mentioned places, please localize this in the Scikit-gstat…
Hosein
  • 1
  • 1
0
votes
1 answer

Extract values from netCDF to location by date

I am new to R and am now looking for a solution to my problem. I have a directory with files in .nc (netCDF) format that contain daily data on sea surface temperature. Each day during the period from December 1, 2019 to August 1, 2021 corresponds to…
Bowe22
  • 1
  • 1
0
votes
1 answer

How to restrict plot area in imshow plot?

Given: - 1- Data data (the entire data is to plotted). 2- A geostationary projection map_proj=crss.Geostationary(central_longitude=82.4) over which the data is to be plotted. 3- An area bounded within the geostationary projection where data is to be…
Pixel_Bear
  • 97
  • 1
  • 1
  • 11
0
votes
1 answer

Include the Matheron estimator in pykrige

Hi I have another question: how can I include the Matheron estimator in my program and still use Pykrige? My program looks like this: from pykrige.ok import OrdinaryKriging from pykrige.kriging_tools import write_asc_grid import…
Weiss
  • 176
  • 2
  • 16
0
votes
1 answer

autofitVariogram in automap: cutoff option doesn't work

I'd like to create a spatial separation distance up to which point pairs are included in semivariance estimates (cutoff function in variogram {gstat}), but using autofitVariogram in automap package. Despite the use of miscFitOptions function nothing…
Leprechault
  • 1,531
  • 12
  • 28
0
votes
1 answer

Cokrigging in R

I have tried to Cokrig based on the following code and everything seems to work upto fitting of coV and there was an error when I try to generate the interpolation. Here is the code that I accustomed to my…
Yon Balcha
  • 29
  • 6
0
votes
0 answers

How to quantify magnitude of change in dataset containing base values of 0?

I have a dataset with the low-water and high-water surface area of lakes/ponds within a delta for each year. These lakes can undergo substantial change from year to year, and sometimes can dry out completely. As such, surface area can have values of…
Nic
  • 3
  • 1
1
2
3