mapproject
On Wed, Jun 30, 2010 at 1:51 PM, Agustin Lobo <alobolistas at gmail.com> wrote:
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()?
Use spTransform from package:rgdal? That's what I use.
examples in ?spTransform include:
state.ll83 <- spTransform(states, CRS("+proj=longlat +ellps=GRS80"))
so things should just plug in.
Barry