[R-sig-Geo] mapproject

[Barry Rowlingson]

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