[R-sig-Geo] function to convert to/from U.S. State Plane Coordinates?

[Waichler, Scott R]

I am working in Linux.  

I clicked on the ESRI:102749 option on the spatialreference.org website that Rob referred to, and got a spec that looks like it should work inside CRS(), but something seems to be amiss.

> library(rgdal)
> pnw.spc <- spTransform(pnw, CRS = CRS("+proj=lcc +lat_1=45.83333333333334 +lat_2=47.33333333333334 +lat_0=45.33333333333334 +lon_0=-120.5 +x_0=500000.0000000002 +y_0=0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs"))
Error in spTransform(pnw, CRS = CRS("+proj=lcc +lat_1=45.83333333333334 +lat_2=47.33333333333334 +lat_0=45.33333333333334 +lon_0=-120.5 +x_0=500000.0000000002 +y_0=0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs")) :
  load package rgdal for spTransform methods

> 
> Are you in Windows or a Linux type computer, please?
> 
> 
> 
> -----Original Message-----
> From: r-sig-geo-bounces at r-project.org [mailto:r-sig-geo-bounces at r-
> project.org] On Behalf Of Waichler, Scott R
> Sent: Monday, September 16, 2013 6:49 PM
> To: Rob Robinson
> Cc: r-sig-geo at r-project.org
> Subject: Re: [R-sig-Geo] function to convert to/from U.S. State Plane
> Coordinates?
> 
> 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.Waichler at 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 at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo