[R-sig-Geo] Plotting a projected map in UTM coordinates; modifying CRS

[Edzer Pebesma]

Attempts to change a projection of spatial data by modifying the
proj4string, such as in:

On 03/25/2012 03:55 PM, Manuel Spínola wrote:
...
>> Am 25.03.2012 05:10, schrieb Manuel Spínola:
>>
>>> Dear list members,
>>>
>>> I am projecting a map to UTM but when I plot the map the axes display
>>> geographic coordinates, do you know why?
>>>
>>>  cr =
>>>>
>>> readOGR(dsn="/Users/**manuelspinola/Documents/**
>>> ProyectosRespacial/EvalHab_**reporte",
>>> layer="Cr_wgs84_meso")
...

>>> proj4string(cr)<-CRS("+proj=**tmerc +lon_0=-84 +lat_0=0 +x_0=500000
>>> +k=0.9999
>>> +datum=WGS84")

and up to now did not generate a warning. This will change; Roger
committed a change to sp with which such attempts now result in a
warning, and a pointer to rgdal::spTransform; see

http://lists.r-forge.r-project.org/pipermail/rspatial-devel/2012-March/000064.html

This warning can be turned off, but is on by default. I think this is an
improvement, but am open to hear other opinions.
-- 
Edzer Pebesma
Institute for Geoinformatics (ifgi), University of Münster
Weseler Straße 253, 48151 Münster, Germany. Phone: +49 251
8333081, Fax: +49 251 8339763  http://ifgi.uni-muenster.de
http://www.52north.org/geostatistics      e.pebesma at wwu.de