# SpatialPolygonsDataFrame object from textfiles of coordinates
# ppt example of 5 regions 
R1_R5ppt.coord = read.table("R1_R5ppt_coord.txt")
# Extract the coordinates and create Polygons objects
R1ppt.poly = Polygon(R1_R5ppt.coord[1:6,])
R1ppt.polys = Polygons(list(Polygon(R1ppt.poly)), ID="R1")
R2ppt.polys = Polygons(list(Polygon(R1_R5ppt.coord[7:12,])), ID="R2")
R3ppt.polys = Polygons(list(Polygon(R1_R5ppt.coord[13:19,])), ID="R3")
R4ppt.polys = Polygons(list(Polygon(R1_R5ppt.coord[20:25,])), ID="R4")
R5ppt.polys = Polygons(list(Polygon(R1_R5ppt.coord[26:30,])), ID="R5")
# Converting the Polygons objects into a SpatialPolygons object
R1_R5ppt.SP = SpatialPolygons(list(R1ppt.polys,R2ppt.polys,R3ppt.polys,R4ppt.polys,R5ppt.polys))
# Drawing a map of the regional system
plot(R1_R5ppt.SP, axes=TRUE)