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

Reply via email to