> The maptools package has a function for computing the centroids of > polygons - install that and try help.search("centroid"). x <- read.shape(system.file("shapes/sids.shp", package="maptools") mycentroid <- get.Pcent(x)