[R-sig-Geo] Help with project function in rgdal library

[Patrick Leahy]

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.) The
text of the .prj file is included below. Any help would be greatly
appreciated!

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]]