Ich muss die Zentroide einer Reihe von räumlichen Zonen basierend auf einem separaten Populationsgitter-Dataset berechnen. Ich bin dankbar für eine Anleitung, wie dies im folgenden Beispiel erreicht werden kann.

Danke im Voraus.

require(raster)
require(spdep)
require(maptools)

dat <- raster(volcano)   # simulated population data
polys <- readShapePoly(system.file("etc/shapes/columbus.shp",package="spdep")[1])

# set consistent coordinate ref. systems and bounding boxes
proj4string(dat) <- proj4string(polys) <- CRS("+proj=longlat +datum=NAD27")
extent(dat) <- extent(polys)

# illustration plot
plot(dat, asp = TRUE)
plot(polys, add = TRUE)