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