[R-sig-Geo] mapproject

Agustin Lobo alobolistas at gmail.com
Wed Jun 30 14:51:14 CEST 2010


I have a SpPolDF in geographic coordinates (eugrd025DF) and must reproject it
to the projection of a raster object (Br):
> projection(eugrd025DF)
[1] "+proj=longlat +ellps=WGS84"

> projection(Br)
[1] "+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000
+ellps=GRS80 +units=m +no_defs"

Can I do this with mapproject (package mapproj) os is it better to
export to shape file, reproject
with a GIS and import back? If it's possible/convenient, how to write
the parameters argument
in mapproject()?

Thanks



More information about the R-sig-Geo mailing list