transform shapefiles in UTM to lat/long
On Tue, 21 Mar 2006, Fredy O. Mejia wrote:
Dear all: I have a shapefile in UTM coordinate system and I would like to transform it to Lat/Log coordinates (WSG84). The package PBSmapping contains function convUL to transform between the two coordinate systems when data is in the form of a data frame with attributes specifying the coordinate system. However, when shapefiles are imported using function read.shape (package maptools), a list instead of a data frame is generated. I could extract all the polygons in the list and transform them, individually, but I hope there is an easier way of doing that. Besides, the shapefile i'm trying to convert has two different UTM zones (is from Guatemala, and it covers zones 15 and 16). Any sugestions would be greatly appreciated.
My suggestion would be to read the shapefiles using the readOGR() function in the rgdal package, and then use the transform() method in that package to change from the appropriate UTM to geographical coordinates. Because I think your question involves more than that (I don't understand how a single shapefile can use two UTM zones), I suggest moving this thread to the R-sig-geo mailing list, and have cc-ed this reply there. Please also say which platform and version of R you are using.
Thanks, Fredy
Roger Bivand Economic Geography Section, Department of Economics, Norwegian School of Economics and Business Administration, Helleveien 30, N-5045 Bergen, Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43 e-mail: Roger.Bivand at nhh.no