[R-sig-Geo] ordinary kriging and elevation with GSTAT

[Paul Hiemstra]

Tobin Cara wrote:
> Hello,
>
> If I generate a variogram with an outside program on the x-y plane, would this variogram be valid for ordinary kriging in gstat? I am giving gstat the x,y,z coordinates and I am trying to predict at the elevations of a digital elevation model.
>
> Thank you in advance for taking the time to answer such a basic question.
>
> Cheers,
>
> Cara
>
> 	[[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>   
Hi,

I would say there is no reason why not. But what do you mean by 
variogram, I think the variogram model. To use the model in R you only 
need to put it into the correct format.

library(automap)
data(meuse)
coordinates(meuse) = ~x+y
kr = autoKrige(meuse)
class(kr$var_model)
kr$var_model
str(kr$var_model)

cheers,
Paul

-- 
Drs. Paul Hiemstra
Department of Physical Geography
Faculty of Geosciences
University of Utrecht
Heidelberglaan 2
P.O. Box 80.115
3508 TC Utrecht
Phone:  +3130 274 3113 Mon-Tue
Phone:  +3130 253 5773 Wed-Fri
http://intamap.geo.uu.nl/~paul