[R-sig-Geo] Average kriging surface over region

[Paul Hiemstra]

Arnout Standaert wrote:
> Dear list,
>
> I'm into the process of learning R's geospatial framework. Currently I'm
> studying scripts that perform a kriging operation on heavy metal
> measurements on a EU scale, see this link:
>
> http://eusoils.jrc.ec.europa.eu/foregshmc/
>
> Now, I would like to derive averages of the fitted model surface for the EU
> member countries. So: take the intersection of the country boundaries and
> the interpolated surface, and derive an average value for that.
> Could someone provide me with some pointers to get started? What R lib would
> be suited for this?
>
> Thanks a bunch,
> Frans
>
> 	[[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,

You can use a polygon as a newdata object in the call to krige. This 
will give you predictions per polygon. Or alternatively, look at the 
overlay() function.

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