function to convert to/from U.S. State Plane Coordinates?

Hi Scott:

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