[R-sig-Geo] 3D kriging with gstat

[Doug Nychka]

Prodromos,

You might also consider the altitude entering just  as a linear term 
in elevation instead of a full spatial variable. Also don't forget that 
aspect (whether the slope faces north or South etc.) may be an important
covariate also. 

The easiest way to fit these kinds of models is to use back fitting ...



Doug


On Tue, 4 Jan 2005, Prodromos Zanis wrote:

> 
> Dear all,
> 
>  I want to do kriging with gstat for x,y,z coordinates (latitude, longtitude 
> and altitude) to interpolate irregularly distributed temperature station data 
> in order to produce a gridded dataset altidude corrected.
> 
> where can I find hints or examples on variogram modeling and kriging in 3D with 
> gstat?
> 
> I look forward for you help. 
> 
> thanks and regards,
> Prodromos Zanis
> 
> 

-- 
----------------------------------------------------------------------------- 
Doug Nychka
    
Institute for Mathematics Applied to Geosciences (IMAGe)             
Computing and Information Systems Laboratory           Email: nychka at ucar.edu 
National Center for Atmospheric Research               Voice: 303-497-1711 
PO Box 3000                                            FAX: 303-497-1333
Boulder, CO  80307-3000                                Web: www.cgd.ucar.edu/~nychka
                                              
Address for overnight mail: 1850 Table Mesa Drive, Boulder, CO 80305