[R-sig-Geo] How to do Zonal Statistics after Kriging in R given Shapefile of Polygon

[Jesse Berman]

James,

Have you considered block kriging with gstat?  This will allow a smoothed
estimate of an area average (e.g. a county or town boundary) as opposed to
over a fine grid.

It's relatively simple, as you merely set the town boundary as your
'newdata,' although you will have to make sure it's in the same projection
as your points.

kr <- krige(POVERTY_INCIDENCE~1, muni.sp, newdata=town, model=v.fin)
spplot(kr, "var1.pred", col.regions = rev(topo.colors(20)))

Jesse




-----
--------------------------------------------------------
Jesse D Berman, PhD
Johns Hopkins Bloomberg School of Public Health
Department of Environmental Health Sciences
Post-Doc 
--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/How-to-do-Zonal-Statistics-after-Kriging-in-R-given-Shapefile-of-Polygon-tp7584330p7584332.html
Sent from the R-sig-geo mailing list archive at Nabble.com.