[R-sig-Geo] mapproject

[Michael Sumner]

I think rgdal package would be the preferable way. There are some old
and outdated map projection packages (I admittedly do not not the
details, and use rgdal exclusively).

library(rgdal)
eugrd025DF.laea <- spTransform(eugrd025DF, CRS(projection(Br)))

spTransform() handles all the details of back and forward transforming
(see inv argument in project) so that valid CRS() projection strings
can be used directly to specify reprojections.

I'm not intimate yet with the raster package so I might have missed
something that's possible in the way the CRS string is derived by
projection, but if proj4string(eugrd025DF) gives the same value as
projection(eugrd025DF) then the spTransform approach should be fine.

Cheers, Mike.

On Wed, Jun 30, 2010 at 10: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()?
>
> Thanks
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>