Hi all,

I've had to extract coordinates from a SpatialPolygonsDataFrame for my
work as well.

# example from Sp Vignettes (Pebesma & Bivand)
 Sr1 = Polygon(cbind(c(2, 4, 4, 1, 2), c(2, 3, 5, 4, 2)))
 Sr2 = Polygon(cbind(c(5, 4, 2, 5), c(2, 3, 2, 2)))
 Sr3 = Polygon(cbind(c(4, 4, 5, 10, 4), c(5, 3, 2, 5, 5)))
 Sr4 = Polygon(cbind(c(5, 6, 6, 5, 5), c(4, 4, 3, 3, 4)), hole = TRUE)
 Srs1 = Polygons(list(Sr1), "s1")
 Srs2 = Polygons(list(Sr2), "s2")
 Srs3 = Polygons(list(Sr3, Sr4), "s3/4")
 SpP = SpatialPolygons(list(Srs1, Srs2, Srs3), 1:3)
 plot(SpP, col = 1:3, pbg = "white")
class(SpP)

# extracting coordinates
test<-...@polygons[[1]]
test2<-t...@polygons[[1]] # here '1' is the first polygon in your
shapefile (you can loop through them)
te...@coords

# Works for a SpatialPolygonsDataFrame as well. I've tried it.

Cheers,
Stephanie Melles

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Reply via email to