Hi,
I'm trying to plot a set of geographical coordinates over the polygons in a
shapefile, so I need to transform the coordinates to the same projection.
I'm currently using the project function in the rgdal library, but I can't
seem to figure out how to format the proj string argument correctly. (I made
a couple attempts using http://www.remotesensing.org/geotiff/proj_list/ as a
reference.) I've never done worked with geographical data in R before so I
thought I'd see if someone more experienced could help. The text of the .prj
file is included below.
Thanks!
Patrick
PROJCS["NAD_1983_StatePlane_New_York_Long_Island_FIPS_3104_Feet",
GEOGCS["GCS_North_American_1983",
DATUM["D_North_American_1983",
SPHEROID["GRS_1980",6378137.0,298.257222101]],
PRIMEM["Greenwich",0.0],
UNIT["Degree",0.0174532925199433]],
PROJECTION["Lambert_Conformal_Conic"],
PARAMETER["False_Easting",984250.0],
PARAMETER["False_Northing",0.0],
PARAMETER["Central_Meridian",-74.0],
PARAMETER["Standard_Parallel_1",40.66666666666666],
PARAMETER["Standard_Parallel_2",41.03333333333333],
PARAMETER["Latitude_Of_Origin",40.16666666666666],
UNIT["Foot_US",0.3048006096012192]]
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
