[R-sig-Geo] temperature interpolation -- was: 3D kriging with gstat

[Edzer Pebesma]

On 01/26/2011 02:45 PM, piero campa wrote:
> 
> This is more about kriging theory actually,
> but I can't understand in which different ways the elevation information is
> taken into account in the prediction process when e.g.:
> 
>  - I apply ordinary cokriging on input data which include 3D (x,y,h) spatial
> coordinates;

If your temperature is (near) surface temperature, then no, because you
want to predict surface temperature, not temperature for arbitrary h values.

>  - I apply universal cokriging on input data with 2D (x,y) spatial
> coordinates and using the elevation h as external drift on the output grid.

This sounds more reasonable, yes.

-- 
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.de
http://www.52north.org/geostatistics      e.pebesma at wwu.de