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

[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 the lines are. Depending on your projection, you may find it better to specify the graticule lines.

cheers
Josh

On Jul 4, 2011, at 10:12 PM, 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_unnamed
> ellipse",DATUM["D_unknown",SPHEROID["Unknown",6370997,0]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["latitude_of_origin",45],PARAMETER["central_meridian",-100],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]
> 
> I'ld hate to install arcinfo to generate 1 plot! (plus the plot is
> already nice as it is!!! just need the axes!)
> 
> thanks,
> I'll summarize!
> 
> H
> 
> -- 
> Horacio Samaniego
> 
> Center for Nonlinear Studies
> Los Alamos National Laboratory
> MS B258
> Los Alamos, New Mexico 87545
> 
> Instituto de Silvicultura
> Universidad Austral de Chile
> Valdivia, Chile
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo

Josh London
Wildlife Biologist
National Marine Mammal Laboratory
NOAA
Seattle, Washington