function to convert to/from U.S. State Plane Coordinates?
On Tue, 17 Sep 2013, Waichler, Scott R wrote:
My system is Redhat, RHEL5.
The package does document all of these requirements:
file.show(system.file("README", package="rgdal"))
or:
https://r-forge.r-project.org/scm/viewvc.php/pkg/inst/README?view=markup&root=rgdal
or read the inst/README file from the source package tarball. Note that
the Fedora/RHEL case is covered specifically.
Nothing about projections is ever simple, nobody should assume that it is.
The assumption is very anachronistic, as clear standards have only been
available recently (say < 20 years, even < 10 years for GNSS/GPS), and
most maps are older; mapping agencies value continuity, and use their
legacy definitions, which vary greatly from place to place, and are often
not well defined.
Note also the link on the Spatial Task View to:
http://spatialreference.org/
as a resource for investigating which incantation may match your data. As
has been mentioned, State Plane contains hundreds of candidate
definitions, depending among other things on the part of the state
concerned and the datum in use at the time the data were created.
Searching this website for "Washington" - your location state - gives 18
hits in EPSG, and only you can know which is appropriate.
Hope this clarifies,
Roger
str(pnw)
List of 4 $ x : num [1:893] -117 -117 -117 -117 -117 ... $ y : num [1:893] 46 46 45.9 45.9 45.9 ... $ range: num [1:4] -125 -111 42 49 $ names: chr [1:7] "idaho" "oregon" "washington:san juan island" "washington:lopez island" ... - attr(*, "class")= chr "map" So, it sounds like there is extra proj4 stuff that I need to install beyond the base installation I already did, which also enabled installation of rgdal? Scott
-----Original Message----- From: Hodgess, Erin [mailto:HodgessE at uhd.edu] Sent: Monday, September 16, 2013 5: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
Roger Bivand Department of Economics, NHH Norwegian School of Economics, Helleveien 30, N-5045 Bergen, Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43 e-mail: Roger.Bivand at nhh.no