Stack Overflow
Stack Overflow

Questions tagged [spatial-interpolation]

Estimating new data points based on existing spatial observations

Spatial interpolation is a method of estimating new data points based on existing spatial interpolation. This tag focusses on spatial interpolation (data with a geographic component), although interpolation in a general also applies to e.g. timeseries.

Commonly, spatial interpolation involves using of the surrounding observations to estimate a new data point. This can be done using some kind of weighted mean (inverse distance weighted interpolation), by fitting a mathematical function through the existing points (splines), or a combination of both (kriging with external drift).

More information regarding which packages can be used to perform interpolation can be found on the Spatial Task View on CRAN, but in general the gstat, automap, Fields, and geoR packages are a good start for a wide range of spatial interpolation methods.

119 questions
1
vote
0 answers

R Spatial Interpolation of irregular grid data in R

Probably it is a newbie question, but I search for a while and unfortunately I did not find anything, which fits to my problem. I want to do a spatial interpolation between one finer irregular grid to one coarser irregular gird. The Problem: I have…
Boje
  • 11
  • 3
1
vote
1 answer

Interpolation technique for weirdly spaced point data

I have a spatial dataset that consists of a large number of point measurements (n=10^4) that were taken along regular grid lines (500m x 500m) and some arbitrary lines and blocks in between. Single measurements taken with a spacing of about 0.3-1.0m…
Leo
  • 11
  • 2
1
vote
1 answer

Spatial Overlay for Kriging Predictions in R

I have data that are structured as follows: winter.dat<-structure(list(ELON = c(-98.02325, -96.66909, -99.33808, -98.70974, -98.29216, -97.08568, -99.90308, -100.53012, -99.05847, -95.86621, -97.25452, -102.49713, -96.63121, -97.69394, -96.35404,…
Joe Ripberger
  • 569
  • 1
  • 7
  • 19
1
vote
0 answers

Vectorize the generation of a 3D matrix made within a loop up to now

Up to now, I used 3 loops to build a matrix, but it takes such a long time that I would like to vectorize it. But either it's not possible, or I took the wrong way to do it. Anyway, I am in front of a wall, and I need your help. Here are the loops…
Nenuneix
  • 11
  • 1
1
vote
0 answers

scipy.interpolate.griddata bombs out when trying to interpolate over a rectangular array

I mark the bad pixels on an image with a boolean array. A bad pixel has a value of True and a good one a value of False. I use this routine for various images, none with the same dimensions. The code I use is: image =…
user2739303
  • 31
  • 4
1
vote
0 answers

Algorithm for grid interpolation/regression given points with holes

I have an image containg a set of identical objects which are disposed into a regular grid fashon. The number of these objects is unknown. I need to find in a really accurate way the bounding rect of that set of objects and then generate the grid…
Alfatau
  • 233
  • 2
  • 8
1
vote
1 answer

Can regression and IDW spatial interpolation be done in one model in R?

I do spatial modelling of variable T (temperature). I use what is commonly used in literature - perform regression (using variables like altitude etc.) and then spatially interpolate the residuals using IDW. R package gstat seems to have this…
Tomas
  • 57,621
  • 49
  • 238
  • 373
1
vote
1 answer

Issue with cross-validation using the automap package

I want to do a cross-validation for the ca20-Dataset from the geoR package. With for example the meuse-dataset, this works fine, but for this dataset, I encounter a strange problem with the dimensions of the …
P_M
  • 328
  • 1
  • 7
1
vote
2 answers

grid point interpolation of wrf data to lat lon in python with numpy and scipy

Im new to python and I'm trying to write a function that will take a numpy array from a netcdf file with dimensions [time,height,longitude,latitude] and interpolate the function to a specified lat and lon. I have looked into scipy.interpolate but…
Barry Baker
  • 19
  • 1
  • 3
1
vote
1 answer

What are my options when I want to find a good interpolation for a coordinate in an irregular grid of values?

I have a list points and their values that cover an area like so: #x, y, value list = [ (1.123, 1.232, 15), (1.121, 2.323, 12), (1.124, 3.451, 17), (1.120, 4.593, 06), (1.121, 5.638, 09), (2.123, 1.232, 19), (2.121, 2.323, 72), (2.124, 3.451,…
AME
  • 2,499
  • 5
  • 29
  • 45
1
vote
2 answers

How I can do Spline Interpolation in R

I am a newbie in R and trying to do interpolation in R(sorry if it seems a lame question) . I have some Points(Lat,Lan and one integer value for each point) that I need to interpolate them. I am a GIS user I did it in ARCGIS but ArcGIS doesn't give…
Hank
  • 70
  • 2
  • 8
0
votes
1 answer

Inverse Distance Weighting with gstat

I am attempting to execute an idw interpolation with gstat but am unable to. text…
arkriger
  • 207
  • 2
  • 7
0
votes
1 answer

Spatial Interpolation Map Temperature using ggplot in R

I would like to create 3 maps of spatial interpolation of a country (like belgium) using the GADM geopkg data that can be found there : https://gadm.org/download_country.html. The idea would be to create 3 interpolated maps : Voronoi…
RedWest
  • 1
0
votes
1 answer

2D slice from 3D volume (interpolation)

I have a 3D volume (black) with each voxel having a known 3D coordinates (x,y,z), and a corresponding intensity value f(x,y,z). Now I like to get a slice (blue) from the 3D volume. All (xq, yq, zq) coordinates of each voxel of the slice is know, and…
Nick X Tsui
  • 2,737
  • 6
  • 39
  • 73
0
votes
0 answers

Unable to create raster files using writeGDAL function

I am trying to perform IDW interpolation on precipitation data in R using the 'gstat' and 'raster' packages. However, I am having trouble creating raster files from my interpolated data using the 'writeGDAL' function from the 'rgdal' package. Here…
Rana
  • 1
  • 1
1 2 3 4
5
6 7 8