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:hodge...@uhd.edu]
Sent: Monday, September 16, 2013 5:12 PM
To: Waichler, Scott R; Rob Robinson
Cc: r-sig-geo@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.waich...@pnnl.gov]
Sent: Monday, September 16, 2013 7:10 PM
To: Hodgess, Erin; Rob Robinson
Cc: r-sig-geo@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-boun...@r-project.org [mailto:r-sig-geo-bounces@r-
project.org] On Behalf Of Waichler, Scott R
Sent: Monday, September 16, 2013 6:49 PM
To: Rob Robinson
Cc: r-sig-geo@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.waich...@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@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@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.biv...@nhh.no
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
