Ok, now I know about the existence of EPSG codes, but apparently my rgdal installation can't find the list of them? Do I need to find a certain file in my system's PROJ4 installation, then tell that to rgdal?
> library(maps) > library(rgdal) > pnw <- map('state', region = c('washington', 'oregon', 'idaho'), plot=F) > epsg <- make_EPSG() Error in make_EPSG() : Error opening epsg file > pnw.spc <- spTransform(pnw, CRS=CRS("+init=epsg:102749")) Error in spTransform(pnw, CRS = CRS("+init=epsg:102749")) : error in evaluating the argument 'CRSobj' in selecting a method for function 'spTransform': Error in CRS("+init=epsg:102749") : no system list, errno: 13 --Scott > try following the instructions > here http://www.gpsmeter.com/index.php?page=faq14 to find an epsg code > then something along the lines of > spTransform(data, CRS=CRS("+init=epsg:xxx")) > should work? > > On 17 September 2013 00:09, Waichler, Scott R <scott.waich...@pnnl.gov> > wrote: > Erin, thanks, I'm getting closer. Rgdal's projInfo() lists UTM, which I > believe is closely related to the State Plane Coordinate system, but I > still don't see a direct transformation to SPC in R. --Scott Waichler > > > > Could you mean spTransform from rgdal, please? > > > > Hi, Does an R function exist to convert to and from the State Plane > > Coordinates system used in the U.S.? My search with Rseek was > fruitless. _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo