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

[Hodgess, Erin]

I think that you need to download some proj4 stuff.  If you're on Ubuntu, go to the Ubuntu package search page and you'll find it.

Also, what is the str(pnw.spc) please?

Thanks,
Erin


-----Original Message-----
From: Waichler, Scott R [mailto:Scott.Waichler at pnnl.gov] 
Sent: Monday, September 16, 2013 7:10 PM
To: Hodgess, Erin; 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?

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