Converting Spatial Polygons to data frame for krige.conv
Thanks Tom,
On 30/05/2011, at 2:24 PM, Tom Gottfried wrote:
do you want the coordinates of all the vertices of your polygons?
Have a look at ?geometry. But this doesn't take you down to a mere list of coordinates. There might be a more elegant way, but you can try that little function to achieve this:
asText.SpatialPolygonsDataFrame <- function(obj){
coords <- matrix(double(), ncol=3)
for (Polygons in obj at polygons){
for (Polygon in Polygons at Polygons){
coords <- rbind(coords,
cbind(as.integer(Polygons at ID), Polygon at coords),
c(as.integer(Polygons at ID), NA, NA))
}
}
return(coords)
}
Perfect. It's where I thought I was going, but it probably would have taken me about two days! cheers Ben