The Spatial task view on your nearest CRAN does say: "The rgdal package provides bindings to GDAL-supported raster formats and OGR-supported vector formats. It contains functions to write raster files in supported formats. The package also provides PROJ.4 projection support for vector objects ...". A Google search on "R spatial projection" reaches the task view as first hit.

Roger

---
Roger Bivand, NHH, Helleveien 30, N-5045 Bergen,
Roger.Bivand@nhh.no



-----Original Message-----
From: r-sig-geo-bounces@stat.math.ethz.ch on behalf of Agustin Lobo
Sent: Wed 2010-06-30 20:25
To: Michael Sumner; b.rowlingson@lancaster.ac.uk
Cc: r-sig-geo
Subject: Re: [R-sig-Geo] mapproject
 
True, thanks. I actually think that I have used spTransform in the past.

The problem was that I searched the help system for "reprojection"
and "projection", and spTransform did not come up.
Maybe both keywords can be included in the documentation of spTransform ?

Agus

2010/6/30 Michael Sumner <mdsumner@gmail.com>:
> 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@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@stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo


	[[alternative HTML version deleted]]

