On Wed, Jun 30, 2010 at 1:51 PM, Agustin Lobo <alobolis...@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
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
