Lets say I have a several polygons stored in a SpatialPolygonsDataFrame
object. Using columbus¹ in maptools as an example:
setwd(system.file("shapes", package="maptools"))
columbus <- readShapeSpatial("columbus.shp²)
I want to extract a polygon 30 from columbus¹, and assign it to a new
SpatialPolygonsDataFrame object. I¹ve tried
newSPDF <- Polygons(list(columbus@polygons[[30]]), "29²)
but this didn¹t work. How can I do this?
Thanks
Ross
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo