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@r-project.org > https://stat.ethz.ch/mailman/listinfo/r-sig-geo Josh London Wildlife Biologist National Marine Mammal Laboratory NOAA Seattle, Washington _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
