An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20140603/63b23b51/attachment.pl>
(Linear) interpolation in three (or more) dimensions
3 messages · Mark Payne, Pascal Oettli, Edzer Pebesma
Hello Mark, The DINEOF method is designed for problems as yours. There is an implementation in R here (http://menugget.blogspot.de/2012/10/dineof-data-interpolating-empirical.html) but it doesn't directly work on Raster* object. Hope this helps, Pascal
On Tue, Jun 3, 2014 at 3:27 PM, Mark Payne <markpayneatwork at gmail.com> 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
Pascal Oettli Project Scientist JAMSTEC Yokohama, Japan
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>