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]]
_______________________________________________
R-sig-Geo mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
