I feel like I have done this before, so I am sorry if I have made a stupid mistake, however, I am out of ideas.
I am trying to make a shapefile out of a list of polygons. I want to include the a data frame as well, but I can't seem to create a SpatialPolygons-Class longer than one object. The one object has all the individual features when I map it, but it is only one object. Here is my code: xx = Polygons(polyList, "FUBAR") spPolys = SpatialPolygons(list(xx)) polyData = (The data frame I want to use) dfReturn = SpatialPolygonsDataFrame(spPolys, polyData, match.ID=F) Any ideas? If I replace the polyData with data.frame(ID = 1) stupid thing works! [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo