Re: [R-sig-Geo] plotting latlon axes of sp object of different projection

2011-07-08 Thread Roger Bivand
On Tue, 5 Jul 2011, Horacio Samaniego wrote: Thanks a lot Roger, This is a good recipe, I'ld never have figured by my own! To fine tune this, I need to place the labels on the other sides (north, east) of the graph. Modifying the pos column of the grdat object (output of gridat()) does not do

Re: [R-sig-Geo] plotting latlon axes of sp object of different projection

2011-07-06 Thread Josh London
Horacio See this thread from December: http://www.mail-archive.com/r-sig-geo@r-project.org/msg00154.html I haven't made much progress or further tested this approach since then. It seems to work for me and my needs. The example with the North Carolina map relies on pretty() to determine where

Re: [R-sig-Geo] plotting latlon axes of sp object of different projection

2011-07-05 Thread Horacio Samaniego
Thanks a lot Roger, This is a good recipe, I'ld never have figured by my own! To fine tune this, I need to place the labels on the other sides (north, east) of the graph. Modifying the pos column of the grdat object (output of gridat()) does not do the trick. I imagine that any modification shou

Re: [R-sig-Geo] plotting latlon axes of sp object of different projection

2011-07-05 Thread Roger Bivand
On Mon, 4 Jul 2011, Horacio Samaniego wrote: I need to add a lat/lon axes to an equal area plot of a SpatialLinesDataFrame object currently in units of meters. Is that achievable using sp/maptools? The .prj file where the object comes from is: PROJCS["Lambert_Azimuthal_Equal_Area",GEOGCS["GCS_u

[R-sig-Geo] plotting latlon axes of sp object of different projection

2011-07-04 Thread Horacio Samaniego
I need to add a lat/lon axes to an equal area plot of a SpatialLinesDataFrame object currently in units of meters. Is that achievable using sp/maptools? The .prj file where the object comes from is: PROJCS["Lambert_Azimuthal_Equal_Area",GEOGCS["GCS_unnamed ellipse",DATUM["D_unknown",SPHEROID["Unkn