Message-ID: <AANLkTimmp24RzE7dJwSvTf4gUw2lKt4L13oWH2DpZK37@mail.gmail.com>
Date: 2010-06-30T13:25:42Z
From: Barry Rowlingson
Subject: mapproject
In-Reply-To: <AANLkTinKL_lozMQ2mxaxOnpq2Marh_pSOaNJnMh0riYc@mail.gmail.com>
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