Areal-weighted interpolation of polygon data
Areal-weighted interpolation of polygon data
st_interpolate_aw(x, to, extensive, ...)
x |
object of class |
to |
object of class |
extensive |
logical; if TRUE, the attribute variables are assumed to be spatially extensive (like population) and the sum is preserved, otherwise, spatially intensive (like population density) and the mean is preserved. |
... |
ignored |
nc = st_read(system.file("shape/nc.shp", package="sf")) g = st_make_grid(nc, n = c(20,10)) a1 = st_interpolate_aw(nc["BIR74"], g, extensive = FALSE) sum(a1$BIR74) / sum(nc$BIR74) # not close to one: property is assumed spatially intensive a2 = st_interpolate_aw(nc["BIR74"], g, extensive = TRUE) # verify mass preservation (pycnophylactic) property: sum(a2$BIR74) / sum(nc$BIR74) a1$intensive = a1$BIR74 a1$extensive = a2$BIR74 plot(a1[c("intensive", "extensive")], key.pos = 4)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.