Something like this should do the trick....

plot(shp) # Plot your shapefile in Lambert projection
plot(spTransform(pts, shp@proj4string), add = T)  # Add points converted to
Lambert

Cheers, Lyndon

On Sun, Feb 26, 2012 at 12:24 PM, <gregor.hochschild@gmx.de> wrote:

> Hi,
>
> I am not sure whether I am expressing myself correctly but I have the
> following problem:
> The shapefile I have is in Lambert Conformal Conic with coordinates like
> this (I not even sure whether the Lambert projection is the problem):
>       min       max
> x 913175.1 1067382.5
> y 120121.9  272844.3
>
> Plotting works fine but now I want to add points based on longitude and
> latitude (NYC area around -73 and 40).
>
> How can I convert my shapefile to longitude and latitude coordinates? Or
> alternatively, how can I convert my point locations?
>
>
> Thanks!
>
>
>
> The shapefile is from here:
> http://www.nyc.gov/html/dcp/html/bytes/dwndistricts.shtml
> (NYC Census Track 2010 file)
>
> Here are more information about the shapefile:
> http://www.nyc.gov/html/dcp/html/bytes/meta_dis_censustracts2010.shtml
>
> Map Projection Name: Lambert Conformal Conic
> Lambert Conformal Conic:
> Standard Parallel: 40.666667
> Standard Parallel: 41.033333
> Longitude of Central Meridian: -74.000000
> Latitude of Projection Origin: 40.166667
> False Easting: 984250.000000
> False Northing: 0.000000
> --
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

	[[alternative HTML version deleted]]

