(Linear) interpolation in three (or more) dimensions
package gstat offers simple 3D interpolation (idw, kriging), e.g. library(gstat) example(gstat3D) and also spatio-temporal kriging; some memory limitations there may be dealt with by using local kriging (see ?krigeST).
On 06/03/2014 08:27 AM, Mark Payne wrote:
Hi, I have a dataset that consists of satellite observations on a regular lon-lat grid at taken at regularly time pointts. Unfortunately due to cloud cover, there are gaps dispersed semi-randomly throughout the 3D data array. I would like to fill these gaps. The "best" solution I guess is space-time kriging, if memory limitations permit. But before I go that far, I would like to do a simple 3D linear interpolation to fill the gaps i.e interpolating the temperature at the missing (lon, lat, time) points. The akima package does an excellent job of this in 2D, but doesn't support higher dimensionality. Does anyone know of a function that supports three or more dimensions? Mark [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster Heisenbergstra?e 2, 48149 M?nster, Germany. Phone: +49 251 83 33081 http://ifgi.uni-muenster.de GPG key ID 0xAC227795 -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 490 bytes Desc: OpenPGP digital signature URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20140603/a1bcba40/attachment.bin>