Hello all. I'm trying to write a new shapefile and am running into a weird issue. After I import the new shapefile into QGIS, I can see all the polygons in the attribute table, but none of them show up on the map no matter how I play with the projections. 

Here is some code that isolates the problem:
library("maptools")
spdf_1 <- readShapePoly("original.shp")
spdf_2 <- readShapePoly("original.shp", proj4string=CRS(as.character("+proj=utm +zone=19 +south +datum=WGS84 +units=m +no_defs")))
writeSpatialShape(spdf_1, "output_1")
writeSpatialShape(spdf_2, "output_2")

As you can see, I tried readShapePoly with and without CRS information, but both output files exhibit the same problem. I've used writeSpatialShape successfully in the past, and I don't know what I'm doing differently here.

Thanks in advance for your help.

Sam Rabin
     Graduate student
         Princeton University
         Ecology & Evolutionary Biology
	[[alternative HTML version deleted]]

