[R-sig-Geo] help with args to CRS(), spTransform

Michael Denslow michael.denslow at gmail.com
Thu Oct 29 14:38:03 CET 2009


Dear List,

I am attending to transform a set of points that are originally in
latitude longitude (unprojected) to a Lambert Conformal Conic
projection.
I am a bit confused about the arguments to CRS and would appreciate any help.

I have an old file that contains point distances that were measured
with the following projection information. This is from ArcMap which I
no longer use. I need to use this old file to compare with my new
distances.

Projected Coordinate System: NAD_1983_StatePlane_North_Carolina_FIPS_3200
Projection: Lambert_Conformal_Conic
False_Easting: 609601.22000000
False_Northing: 0.00000000
Central_Meridian: -79.00000000
Standard_Parallel_1: 34.33333333
Standard_Parallel_2: 36.16666667
Latitude_Of_Origin: 33.75000000
Linear Unit: Meter

Geographic Coordinate System: GCS_North_American_1983
Datum: D_North_American_1983
Prime Meridian: Greenwich
Angular Unit: Degree

I would like the new point distance measurements that I am doing in
SPATSTAT to be comparable to the old ones that used the projection
above (in meters?).

Here is what I have attempted...

library(sp)
library(rgdal)
# dput(herb.sp) I hope this works I have never used dput() before!

new("SpatialPoints"
    , coords = structure(c(-81.6809, -78.9455, -77.879, -83.1031, -78.6634,
-79.057, 36.2147, 35.9982, 34.227, 35.2571, 35.7853, 35.9097), .Dim = c(6L,
2L), .Dimnames = list(NULL, c("coords.x1", "coords.x2")))
    , bbox = structure(c(-83.1031, 34.227, -77.879, 36.2147), .Dim = c(2L,
2L), .Dimnames = list(NULL, c("min", "max")))
    , proj4string = new("CRS"
    , projargs = " +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"
)
)

herb.lcc <- spTransform(herb.sp,CRS("+proj=lcc"))


I am not sure that this contains all of the arguments needed. I tried
looking at the PROG.4 manual here:
ftp://ftp.remotesensing.org/proj/OF90-284.pdf , but I am still unclear
how to proceed.

Any help would be most appreciated!
Thanks again,
Michael



-- 
Michael Denslow

Graduate Student & Adjunct Instructor
I.W. Carpenter Jr. Herbarium [BOON]
Department of Biology
Appalachian State University
Boone, North Carolina U.S.A.
-- AND --
Communications Manager
Southeast Regional Network of Expertise and Collections
sernec.org

36.214177, -81.681480 +/- 3103 meters



More information about the R-sig-Geo mailing list