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
2
votes
0 answers

Spatial interpolation with Copula

I was wondering if anyone knows about a library for spatial interpolation using Copulas in Python. I have already tried the copulalib and consequently ambhas, but none of them are suited for spatial interpolation. A reference to the problem can be…
2
votes
0 answers

How to mask data in scipy.interpolate.RectBivariateSpline

I have some curvilinear gridded data, and I'm trying to interpolate this data to a new grid. The data has unrealistically high variables for the masked data, so I tried masking using: temp = np.ma.masked_where(temp > 500, temp) if I do this…
user14241
  • 727
  • 1
  • 8
  • 27
2
votes
1 answer

How can I interpolate a point along a line string using DbGeography?

I have a linestring DbGeography object and I want to find some point along the linestring given a fractional value of the total distance. For example, if I give 0.5, I should get the midpoint of the line as a DbGeography object with a point WKT…
Chadley08
  • 559
  • 1
  • 8
  • 26
2
votes
1 answer

Interpolating Between Two Planes in 3d space

I'm developing a tool that lets you circle/enclose things on a 3d "volume." I want to save time by labelling "slices" 1 and 3, and "filling in" slice 2 from that information. Two easy solutions are to: 1. slice2 = slice1 AND slice3 (gets the overlap…
Frederick
  • 1,271
  • 1
  • 10
  • 29
2
votes
0 answers

interpolation across 6-D data, using MATLAB

I am attempting to interpolate across a 6-dimensional, non-monotonic, real domain. I'm using MATLAB. I am pretty sure that the function "griddatan" will do, but when I try to run it, I get the following error from one of its subroutines: "Error…
2
votes
1 answer

plotting issue with autofitVariogram in automap

I am a newbie in R. I am using autofitVariogram to daily rainfall data of 50 stations.The sample data is provided below.Some of stations have missing values represented by "NaN" values. My question is regarding the variogramfit. The variogram…
user1142937
  • 314
  • 5
  • 19
2
votes
1 answer

Method for interpolating value at the midpoint of a square, given the values, first and second derivatives at the corners?

All the practical examples of spatial interpolation I've been able to find work by sampling additional surrounding points to estimate the derivatives. But is there a simpler method if the derivatives are already known—and if you only need the value…
1
vote
2 answers

Vectorized method to insert arrays of different length at different spots in an ndarray

I'm trying to "expand" a 2D numpy array nadir of points (x,y,z), and fill gaps in space with interpolated points. Where there exist spatial gaps bigger than some tolerance dist, I want to use np.insert to insert the required number of nan rows to…
bt3
  • 156
  • 1
  • 10
1
vote
0 answers

Area Weighted Average Grid Interpolation

I would like to resize an n-dimensional grid similar to scipy.ndimage.affine_transform. For the stake of simplicity, I describe my problem with the 2-Dimensional case. I need to downsample an image with an arbitrary non integer factor by more than 2…
checkThisOut
  • 621
  • 1
  • 5
  • 18
1
vote
2 answers

Drawing something resembling contour lines for an annually advancing ray of fungi?

I would be happy for any tips for drawing something resembling contour lines for the growth advancement lines around a fungal fairy ring that advances annually uphill. This is a plot of the records of fruit bodies I have found from the years…
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)…
1
vote
0 answers

Universal kriging has a long processing time compared to other kriging methods - R

I have a dataset data_sp with 132 observations containing 9 predictors and 1 response variable. I perform a twenty-fold cross validation, using ordinary kriging and universal kriging. The raster grid100 on which I do my kriging, contains 61124…
1
vote
0 answers

Extracting spatially IDW interpolated z-point data using specified x, y inputs?

I was able to do an IDW ppp spatial interpolation in respect to a Datum level based on corresponding longitude, latitude coordinates for a grid I've specified. Although I achieve the desired plot, I am unsure how to input specific long-lat (x-y)…
1
vote
0 answers

Interpolate temperature values (with associated lat and lon values) at different depth intervals using Krig function in the 'fields' package

I have temperature data points in different depth intervals with associated lat long values across the study area. I want to make a raster and then interpolate between cells of the raster where there is no data. I can do it using Krig in the fields…
Leandri
  • 11
  • 2
1
vote
1 answer

xarray interp variable with two dimensional gird into one point

I have a dataset with two-dimension lon/lat and I want to calculate the value of one certain point. The ncfile can be download from ftp://ftp.awi.de/sea_ice/product/cryosat2_smos/v204/nh/LATEST/ and the variable can be retrieved by the following…