[R-sig-Geo] (no subject)

harrykim at stat.Berkeley.EDU harrykim at stat.Berkeley.EDU
Fri Dec 1 19:08:32 CET 2006


Hello R-sig-geo list,

   I've been trying to figure this out by myself but i am getting kind of
lost so i was hoping if you could point me to the right direction.
   I have a shape file for map of berkeley--the one in
california--formatted in long/lat coordinates. I also have set of
locations in long/lat coordinates which i wish to plot over the map of
berkeley.
   So far i have

library("maptools")
plot(readShapePoly("V1 Land.shp"))
points(tmp$long,tmp$lat)

which does the job. But I want to convert the shapefile and the locations
into UTM coordinates and ultimately convert them into meters (or
kilometers); I want to work with the coordinates in meters.
   I have been using readShapePoly() in maptools package to read in the
shapefile which returns spatial polygons data frame. I've read that the
spatial polygons data frame can be converted into UTM using
sptransform() in Rgdal package but i am not so sure what i need to
exactly. I tried to read up on using CRS class (a class used by proj4)
but i was getting lost in the documentation.

Thank you very much
Harry
p.s:The UTM zone for berkeley is 10S




More information about the R-sig-Geo mailing list