function to convert to/from U.S. State Plane Coordinates?
Whoops...I meant pnw -----Original Message----- From: r-sig-geo-bounces at r-project.org [mailto:r-sig-geo-bounces at r-project.org] On Behalf Of Hodgess, Erin Sent: Monday, September 16, 2013 7:12 PM To: Waichler, Scott R; 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 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
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo