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: [EMAIL PROTECTED] ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
