[R-sig-Geo] Average kriging surface over region
Paul Hiemstra
p.hiemstra at geo.uu.nl
Mon Oct 12 16:21:05 CEST 2009
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
More information about the R-sig-Geo
mailing list