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

Josh London josh.london at noaa.gov
Tue Jul 5 23:26:53 CEST 2011


See this thread from December:


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.


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
Seattle, Washington

More information about the R-sig-Geo mailing list